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 31 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 52 of file MC_Efficiency_Analyzer.cc.

References funct::abs(), edm::AssociationMap< Tag >::begin(), gather_cfg::cout, ztail::d, reco::Candidate::daughter(), PbPb_ZMuSkimMuonDPG_cff::deltaR, MillePedeFileConverter_cfg::e, edm::AssociationMap< Tag >::end(), reco::Candidate::eta(), genParticles2HepMC_cfi::genParticles, mps_fire::i, dqmdumpme::k, edm::helpers::KeyVal< K, V >::key, cmsLHEtoEOSManager::l, CaloTowersParam_cfi::mc, HLTTauDQMOffline_cfi::Muons, dqmiodumpmetadata::n, reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::Candidate::pt(), mps_update::status, reco::Candidate::status(), gather_cfg::Tracks, edm::helpers::KeyVal< K, V >::val, and ZMuMuAnalysisNtupler_cff::zMuMu.

                                                                            {
     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 157 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 183 of file MC_Efficiency_Analyzer.cc.

double MCEfficiencyAnalyzer::etacut_

private

Definition at line 183 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<GenParticleCollection> MCEfficiencyAnalyzer::genParticlesToken_

private

Definition at line 180 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::MuonsMapToken_

private

Definition at line 177 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::MuonsToken_

private

Definition at line 176 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nMuMC

private

Definition at line 184 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nMureco

private

Definition at line 184 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nNotMuMatching

private

Definition at line 184 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nSta

private

Definition at line 184 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nTrk

private

Definition at line 184 of file MC_Efficiency_Analyzer.cc.

double MCEfficiencyAnalyzer::ptcut_

private

Definition at line 183 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::StandAloneMapToken_

private

Definition at line 182 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::StandAloneToken_

private

Definition at line 181 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::TracksMapToken_

private

Definition at line 179 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::TracksToken_

private

Definition at line 178 of file MC_Efficiency_Analyzer.cc.

vector<double> MCEfficiencyAnalyzer::v_

private

Definition at line 185 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::zMuMuToken_

private

Definition at line 175 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