Using only final particle in jet clustering algorithms

kgautam · August 16, 2021, 1:48pm

Hi people!

I’m trying to use the jet clustering algorithms using MC particles, I can run the algorithms by referring to other examples, but I can’t find any command to filter the particles I want to use for jet clustering in the analysis.py file. Is there a way that I can use only final particles (status = 1) as inputs in the clustering algos while running analysis.py? Or is the only way to do it is to loop over the jet constituents and keep only the final particles?
Without any interference, my jets have all the particles and therefore energies and momenta far exceeding the logical limits.
Any help is appreciated. Thank you

Cheers,
Kunal

clement.helsens · August 17, 2021, 6:33am

hello,

I believe this is what you need:

github.com

HEP-FCC/FCCAnalyses/blob/master/analyzers/dataframe/MCParticle.h#L38#L43

    
      
            bool  operator() (ROOT::VecOps::RVec<edm4hep::MCParticleData> in);
          };
          
          
/// select MCParticles with transverse momentum greater than a minimum value [GeV]
          struct sel_pt {
            sel_pt(float arg_min_pt);
            float m_min_pt = 20; //> transverse momentum threshold [GeV]
            ROOT::VecOps::RVec<edm4hep::MCParticleData>  operator() (ROOT::VecOps::RVec<edm4hep::MCParticleData> in);
          };
          
          
/// select MCParticles with their status
          struct sel_genStatus {
            sel_genStatus(int arg_status);
            float m_status = 1; //> Generator status
            ROOT::VecOps::RVec<edm4hep::MCParticleData>  operator() (ROOT::VecOps::RVec<edm4hep::MCParticleData> in);
          };
          
          
/// select MCParticles with their PDG id
          struct sel_pdgID {
            sel_pdgID(int arg_pdg, bool arg_chargeconjugate);
            int m_pdg = 13;

Best,
Clement

kgautam · August 17, 2021, 7:59am

Hi Clement,

This worked perfectly. Thank you so much.

Cheers,
Kunal