Inheritance diagram for MCEfficiencyAnalyzer:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

void	endJob () override

	MCEfficiencyAnalyzer (const edm::ParameterSet &pset)

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
double	deltaRStacut_

double	etacut_

EDGetTokenT< GenParticleCollection >	genParticlesToken_

EDGetTokenT< CandMatchMap >	MuonsMapToken_

EDGetTokenT< CandidateCollection >	MuonsToken_

int	nMuMC

int	nMureco

int	nNotMuMatching

int	nSta

int	nTrk

double	ptcut_

EDGetTokenT< CandMatchMap >	StandAloneMapToken_

EDGetTokenT< CandidateCollection >	StandAloneToken_

EDGetTokenT< CandMatchMap >	TracksMapToken_

EDGetTokenT< CandidateCollection >	TracksToken_

vector< double >	v_

EDGetTokenT< CandidateCollection >	zMuMuToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 29 of file MC_Efficiency_Analyzer.cc.

Constructor & Destructor Documentation

MCEfficiencyAnalyzer::MCEfficiencyAnalyzer ( const edm::ParameterSet & pset )

inline

Definition at line 32 of file MC_Efficiency_Analyzer.cc.

                                                     :
     zMuMuToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "zMuMu" ) ) ),
     MuonsToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "Muons" ) ) ),
     MuonsMapToken_( consumes<CandMatchMap>( pset.getParameter<InputTag>( "MuonsMap" ) ) ),
     TracksToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "Tracks" ) ) ),
     TracksMapToken_( consumes<CandMatchMap>( pset.getParameter<InputTag>( "TracksMap" ) ) ),
     genParticlesToken_( consumes<GenParticleCollection>( pset.getParameter<InputTag>( "genParticles" ) ) ),
     StandAloneToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "StandAlone" ) ) ),
     StandAloneMapToken_( consumes<CandMatchMap>( pset.getParameter<InputTag>( "StandAloneMap" ) ) ),
     etacut_( pset.getParameter<double>( "etacut" ) ),
     ptcut_( pset.getParameter<double>( "ptcut" ) ),
     deltaRStacut_( pset.getParameter<double>( "deltaRStacut" ) )
 
   {
     nMuMC =0; nMureco =0; nTrk=0; nSta=0; nNotMuMatching =0 ;
   }

Member Function Documentation

void MCEfficiencyAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

inlineoverride

Definition at line 49 of file MC_Efficiency_Analyzer.cc.

References funct::abs(), edm::AssociationMap< Tag >::begin(), gather_cfg::cout, edmIntegrityCheck::d, reco::Candidate::daughter(), boostedElectronIsolation_cff::deltaR, MillePedeFileConverter_cfg::e, edm::AssociationMap< Tag >::end(), reco::Candidate::eta(), GenHFHadronMatcher_cfi::genParticles, mps_fire::i, gen::k, crabWrapper::key, checklumidiff::l, CaloTowersParam_cfi::mc, Reconstruction_hiPF_cff::Muons, gen::n, reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::Candidate::pt(), mps_update::status, reco::Candidate::status(), gather_cfg::Tracks, and heppy_batch::val.

   {
 
     Handle<CandidateCollection> zMuMu;
     event.getByToken(zMuMuToken_, zMuMu);
 
 
     Handle<CandidateCollection> Muons;
     event.getByToken(MuonsToken_, Muons);
 
     Handle<CandidateCollection> Tracks;
     event.getByToken(TracksToken_, Tracks);
 
     Handle<CandidateCollection> StandAlone;
     event.getByToken(StandAloneToken_, StandAlone);
 
     Handle<CandMatchMap> MuonsMap;
     event.getByToken(MuonsMapToken_, MuonsMap);
 
     Handle<CandMatchMap> TracksMap;
     event.getByToken(TracksMapToken_, TracksMap);
 
     Handle<CandMatchMap> StandAloneMap;
     event.getByToken(StandAloneMapToken_, StandAloneMap);
 
     Handle<GenParticleCollection> genParticles;
     event.getByToken(genParticlesToken_, genParticles);
 
     //Getting muons from Z MC
     for( unsigned int k = 0; k < genParticles->size(); k++ )
       {
     const Candidate & ZCand = (*genParticles)[ k ];
     int status = ZCand.status();
 
     if (ZCand.pdgId()==23&& status==3 )
       {
         // positive muons
         const Candidate  *  muCand1 =   ZCand.daughter(0);
         if (muCand1->status()==3)
           {
         for (unsigned int d =0 ; d< muCand1->numberOfDaughters(); d++)
           {
             const Candidate * muCandidate = muCand1->daughter(d);
             if (muCandidate->pdgId() == muCand1->pdgId() )
               {
             muCand1 = muCand1->daughter(d);
               }
           }
           }
         // negative muons
         const Candidate  *  muCand2 =   ZCand.daughter(1);
         if (muCand2->status()==3)
           {
         for (unsigned int e =0 ; e< muCand2->numberOfDaughters(); e++)
           {
             const Candidate * muCandidate = muCand2->daughter(e);
             if (muCandidate->pdgId() == muCand2->pdgId() )
               {
             muCand2 = muCand2->daughter(e);
               }
           }
           }
 
         double deltaR_Mu_Sta =0;
         int nMurecoTemp = nMureco;
         // getting mu matched
         CandMatchMap::const_iterator i;
         for(i = MuonsMap->begin(); i != MuonsMap->end(); i++ )
           {
         const Candidate/* & reco = * i -> key,*/ &  mc =  * i -> val;
         if ((muCand1 == &mc) && (mc.pt()>ptcut_) && (std::abs(mc.eta())<etacut_))
           {
             nMuMC++;
             nMureco++;
             break;
           }
           }
         if (nMureco == nMurecoTemp ) // I.E. MU RECO NOT FOUND!!!!
           {
                 int nTrkTemp = nTrk;
         int nStaTemp = nSta;
 
         // getting tracks matched and doing the same, CONTROLLING IF MU IS RECONSTRUCTED AS A TRACK AND NOT AS A MU
         CandMatchMap::const_iterator l;
         for(l = TracksMap->begin(); l != TracksMap->end(); l++ )
           {
             const Candidate /* & Trkreco = * l -> key, */ &  Trkmc =  * l -> val;
             if (( muCand1 == & Trkmc) && (Trkmc.pt()>ptcut_) && (std::abs(Trkmc.eta())<etacut_))
               {
                 nMuMC++;
             nTrk++;
             break;
               }
           }
         // the same for standalone
         CandMatchMap::const_iterator n;
         for(n = StandAloneMap->begin(); n != StandAloneMap->end(); n++ )
           {
             const Candidate & Stareco = * n -> key, &  Stamc =  * n -> val;
             if ((muCand1 == &Stamc ) && (Stamc.pt()>ptcut_) && (std::abs(Stamc.eta())<etacut_))
               {
             nMuMC++;
             nSta++;
             deltaR_Mu_Sta = deltaR(Stareco, *muCand1);
 
             //  cout<<"Ho trovato un sta reco "<<endl;
             break;
               }
           }
         // controlling if sta and trk are reconstrucetd both, and if so the get the deltaR beetween muon MC and reco sta, to controll correlation to this happening
         if ((nSta == nStaTemp + 1) && (nTrk == nTrkTemp + 1 ) )
           {
             nNotMuMatching ++;
             if ((deltaR_Mu_Sta< deltaRStacut_))
               {
             v_.push_back(deltaR_Mu_Sta) ;
                    cout << "Not matching from trk and sta matched to MC mu, to reconstruct a recoMU" << endl;
               }
           }
           }
       }
       }
   }

void MCEfficiencyAnalyzer::endJob ( void )

inlineoverridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 174 of file MC_Efficiency_Analyzer.cc.

References gather_cfg::cout, and mps_fire::i.

   {
 
 
     cout <<"--- nMuMC == "<<nMuMC<<endl;
     cout <<"--- nMureco == "<<nMureco<<endl;
     cout <<"--- nSta == "<<nSta<<endl;
     cout <<"--- nTrk == "<<nTrk<<endl;
     cout <<"--- nNotMuMatching from a trk and sta matched to a Mu MC == "<<nNotMuMatching<<endl;
     if (nMuMC!=0 )
       {
     cout<<" effMu == "<<(double) nMureco/nMuMC<<endl;
     cout<<" effTrk == "<< (double)(nTrk + nMureco) /nMuMC<<endl;
     cout<<" effSta == "<< (double)(nSta + nMureco) / nMuMC<<endl;
       }
 
     for (unsigned int i =0 ; i < v_.size(); ++i )
       {
     cout<<" delta R Mu Sta == "<< v_[i]<<endl;
         }
   }

Member Data Documentation

double MCEfficiencyAnalyzer::deltaRStacut_

private

Definition at line 210 of file MC_Efficiency_Analyzer.cc.

double MCEfficiencyAnalyzer::etacut_

private

Definition at line 210 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<GenParticleCollection> MCEfficiencyAnalyzer::genParticlesToken_

private

Definition at line 207 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::MuonsMapToken_

private

Definition at line 204 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::MuonsToken_

private

Definition at line 203 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nMuMC

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nMureco

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nNotMuMatching

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nSta

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nTrk

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

double MCEfficiencyAnalyzer::ptcut_

private

Definition at line 210 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::StandAloneMapToken_

private

Definition at line 209 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::StandAloneToken_

private

Definition at line 208 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::TracksMapToken_

private

Definition at line 206 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::TracksToken_

private

Definition at line 205 of file MC_Efficiency_Analyzer.cc.

vector<double> MCEfficiencyAnalyzer::v_

private

Definition at line 212 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::zMuMuToken_

private

Definition at line 202 of file MC_Efficiency_Analyzer.cc.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation