L1Analysis::L1AnalysisGenerator Class Reference

#include <L1AnalysisGenerator.h>

Public Member Functions
L1AnalysisGeneratorDataFormat *	getData ()
	L1AnalysisGenerator ()
void	Reset ()
void	Set (const edm::Event &e)
	~L1AnalysisGenerator ()

Private Attributes
L1AnalysisGeneratorDataFormat	generator_

Detailed Description

Definition at line 15 of file L1AnalysisGenerator.h.

Constructor & Destructor Documentation

L1Analysis::L1AnalysisGenerator::L1AnalysisGenerator

(

)

Definition at line 10 of file L1AnalysisGenerator.cc.

{
}

L1Analysis::L1AnalysisGenerator::~L1AnalysisGenerator

(

)

Definition at line 14 of file L1AnalysisGenerator.cc.

{
}

Member Function Documentation

L1AnalysisGeneratorDataFormat* L1Analysis::L1AnalysisGenerator::getData ( void  )

inline

Definition at line 22 of file L1AnalysisGenerator.h.

References generator_.

Referenced by L1NtupleProducer::L1NtupleProducer().

{return &generator_;}

void L1Analysis::L1AnalysisGenerator::Reset ( void  )

inline

Definition at line 20 of file L1AnalysisGenerator.h.

References generator_, and L1Analysis::L1AnalysisGeneratorDataFormat::Reset().

Referenced by L1NtupleProducer::analyzeGenerator().

{generator_.Reset();}

void L1Analysis::L1AnalysisGenerator::Set

(

const edm::Event &

e

)

Definition at line 18 of file L1AnalysisGenerator.cc.

References funct::abs(), reco::LeafCandidate::energy(), genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByLabel(), i, reco::CompositeRefCandidateT< D >::mother(), reco::CompositeRefCandidateT< D >::numberOfMothers(), AlCaHLTBitMon_ParallelJobs::p, SingleMuPt40Fwdv2_cfi_GEN_SIM::pdgId, reco::LeafCandidate::pdgId(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), and reco::LeafCandidate::status().

{
  edm::Handle<reco::GenParticleCollection> genParticles;
   e.getByLabel("genParticles", genParticles);
   for(size_t i = 0; i < genParticles->size(); ++ i) {
     const GenParticle & p = (*genParticles)[i];
     int id = p.pdgId();
     //int st = p.status();  
        if (abs(id) == 13) {
            unsigned int nMo=p.numberOfMothers();
//          std::cout << "id " << id << "; st " << st 
//                          << "; nMo " << nMo << std::endl;
            for(unsigned int i=0;i<nMo;++i){
//              int thisParentID = dynamic_cast <const reco::GenParticle*>(p.mother(i))->pdgId();
//              std::cout << "   mother ID " << thisParentID << std::endl;
            }
        }
            
//
// See if the parent was interesting
        int parentID = -10000;
        unsigned int nMo=p.numberOfMothers();
        for(unsigned int i=0;i<nMo;++i){
            int thisParentID = dynamic_cast
                    <const reco::GenParticle*>(p.mother(i))->pdgId();
//
// Is this a bottom hadron?
            int hundredsIndex = abs(thisParentID)/100;
            int thousandsIndex = abs(thisParentID)/1000;
            if ( ((abs(thisParentID) >= 23) && 
                        (abs(thisParentID) <= 25)) ||
                        (abs(thisParentID) == 6) ||
                        (hundredsIndex == 5) ||
                        (hundredsIndex == 4) ||
                        (thousandsIndex == 5) ||
                        (thousandsIndex == 4) 
                    )
                parentID = thisParentID;
        }
        if ((parentID == -10000) && (nMo > 0)) 
            parentID = dynamic_cast
                    <const reco::GenParticle*>(p.mother(0))->pdgId();
//
// If the parent of this particle is interesting, store all of the info
        if ((parentID != p.pdgId()) &&
            ((parentID > -9999) 
               || (abs(id) == 11)
               || (abs(id) == 13)
               || (abs(id) == 23)
               || (abs(id) == 24)
               || (abs(id) == 25)
               || (abs(id) == 4)
               || (abs(id) == 5)
               || (abs(id) == 6))
            )
        {
         generator_.id.push_back(p.pdgId());
            generator_.status.push_back(p.status());
            generator_.px.push_back(p.px());
            generator_.py.push_back(p.py());
            generator_.pz.push_back(p.pz());
            generator_.e.push_back(p.energy());
            generator_.parent_id.push_back(parentID);
        }
   }

}

Member Data Documentation

L1AnalysisGeneratorDataFormat L1Analysis::L1AnalysisGenerator::generator_

private

Definition at line 25 of file L1AnalysisGenerator.h.

Referenced by getData(), and Reset().