Inheritance diagram for ZMuMu_MCanalyzer:

Public Member Functions
	ZMuMu_MCanalyzer (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 Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

bool	check_ifZmumu (const Candidate dauGen0, const Candidate dauGen1, const Candidate *dauGen2)

void	endJob () override

float	getParticleEta (const int ipart, const Candidate dauGen0, const Candidate dauGen1, const Candidate *dauGen2)

Particle::LorentzVector	getParticleP4 (const int ipart, const Candidate dauGen0, const Candidate dauGen1, const Candidate *dauGen2)

float	getParticlePhi (const int ipart, const Candidate dauGen0, const Candidate dauGen1, const Candidate *dauGen2)

float	getParticlePt (const int ipart, const Candidate dauGen0, const Candidate dauGen1, const Candidate *dauGen2)

Private Attributes
double	alpha_

double	beta_

bool	bothMuons_

double	dREcal_

double	dRHcal_

double	dRTrk_

double	dRVetoTrk_

double	etamax_

double	etamin_

double	etEcalThreshold_

double	etHcalThreshold_

EDGetTokenT< GenParticleCollection >	genParticlesToken_

reco::CandidateBaseRef	globalMuonCandRef_

TH1D *	h_ProbeOk_eta

TH1D *	h_ProbeOk_pt

TH1D *	h_staProbe_eta

TH1D *	h_staProbe_pt

TH1D *	h_trackProbe_eta

TH1D *	h_trackProbe_pt

string	hltPath_

double	isoMax_

double	massMax_

double	massMin_

EDGetTokenT< CandidateView >	muonsToken_

int	n2GlobalMuonsMatched_passedIso

int	n2GlobalMuonsMatched_passedIso2Trg

int	nGlobalMuonsMatched_passed

int	nGlobalMuonsMatched_passedIso

int	nMu0onlyTriggered

int	nMu1onlyTriggered

int	nStaMuonsMatched_passedIso

int	nTracksMuonsMatched_passedIso

int	nZMuMu_matched

int	nZMuSta_matched

int	nZMuTrk_matched

OverlapChecker	overlap_

double	ptmin_

double	ptThreshold_

bool	relativeIsolation_

reco::CandidateBaseRef	standAloneMuonCandRef_

reco::CandidateBaseRef	trackMuonCandRef_

EDGetTokenT< CandidateView >	tracksToken_

EDGetTokenT< GenParticleMatch >	zMuMuMatchMapToken_

EDGetTokenT< CandidateView >	zMuMuToken_

EDGetTokenT< GenParticleMatch >	zMuStandAloneMatchMapToken_

EDGetTokenT< CandidateView >	zMuStandAloneToken_

EDGetTokenT< GenParticleMatch >	zMuTrackMatchMapToken_

EDGetTokenT< CandidateView >	zMuTrackToken_

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 49 of file ZMuMu_MCanalyzer.cc.

Constructor & Destructor Documentation

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

Definition at line 166 of file ZMuMu_MCanalyzer.cc.

References MuonErrorMatrixValues_cff::etaRange, h_ProbeOk_eta, h_ProbeOk_pt, h_staProbe_eta, h_staProbe_pt, h_trackProbe_eta, h_trackProbe_pt, TFileService::make(), n2GlobalMuonsMatched_passedIso, n2GlobalMuonsMatched_passedIso2Trg, nGlobalMuonsMatched_passed, nGlobalMuonsMatched_passedIso, nMu0onlyTriggered, nMu1onlyTriggered, nStaMuonsMatched_passedIso, nTracksMuonsMatched_passedIso, nZMuMu_matched, nZMuSta_matched, nZMuTrk_matched, and MuonErrorMatrixValues_cff::ptRange.

                                                            :
   zMuMuToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuMu"))),
   zMuMuMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuMuMatchMap"))),
   zMuStandAloneToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuStandAlone"))),
   zMuStandAloneMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuStandAloneMatchMap"))),
   zMuTrackToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuTrack"))),
   zMuTrackMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuTrackMatchMap"))),
   muonsToken_(consumes<CandidateView>(pset.getParameter<InputTag>("muons"))),
   tracksToken_(consumes<CandidateView>(pset.getParameter<InputTag>("tracks"))),
   genParticlesToken_(consumes<GenParticleCollection>(pset.getParameter<InputTag>("genParticles"))),
 
   bothMuons_(pset.getParameter<bool>("bothMuons")),
   etamin_(pset.getUntrackedParameter<double>("etamin")),
   etamax_(pset.getUntrackedParameter<double>("etamax")),
   ptmin_(pset.getUntrackedParameter<double>("ptmin")),
   massMin_(pset.getUntrackedParameter<double>("zMassMin")),
   massMax_(pset.getUntrackedParameter<double>("zMassMax")),
   isoMax_(pset.getUntrackedParameter<double>("isomax")),
   ptThreshold_(pset.getUntrackedParameter<double>("ptThreshold")),
   etEcalThreshold_(pset.getUntrackedParameter<double>("etEcalThreshold")),
   etHcalThreshold_(pset.getUntrackedParameter<double>("etHcalThreshold")),
   dRVetoTrk_(pset.getUntrackedParameter<double>("deltaRVetoTrk")),
   dRTrk_(pset.getUntrackedParameter<double>("deltaRTrk")),
   dREcal_(pset.getUntrackedParameter<double>("deltaREcal")),
   dRHcal_(pset.getUntrackedParameter<double>("deltaRHcal")),
   alpha_(pset.getUntrackedParameter<double>("alpha")),
   beta_(pset.getUntrackedParameter<double>("beta")),
   relativeIsolation_(pset.getUntrackedParameter<bool>("relativeIsolation")),
   hltPath_(pset.getUntrackedParameter<std::string >("hltPath")) {
   Service<TFileService> fs;
 
   // binning of entries array (at moment defined by hand and not in cfg file)
   double etaRange[8] = {-2.5,-2.,-1.2,-0.8,0.8,1.2,2.,2.5};
   double ptRange[4] = {20.,40.,60.,100.};
 
   // general histograms
   h_trackProbe_eta = fs->make<TH1D>("trackProbeEta","Eta of tracks",7,etaRange);
   h_trackProbe_pt = fs->make<TH1D>("trackProbePt","Pt of tracks",3,ptRange);
   h_staProbe_eta = fs->make<TH1D>("standAloneProbeEta","Eta of standAlone",7,etaRange);
   h_staProbe_pt = fs->make<TH1D>("standAloneProbePt","Pt of standAlone",3,ptRange);
   h_ProbeOk_eta = fs->make<TH1D>("probeOkEta","Eta of probe Ok",7,etaRange);
   h_ProbeOk_pt = fs->make<TH1D>("probeOkPt","Pt of probe ok",3,ptRange);
 
   // clear global counters
   nGlobalMuonsMatched_passed = 0;
   nGlobalMuonsMatched_passedIso = 0;
   n2GlobalMuonsMatched_passedIso = 0;
   nStaMuonsMatched_passedIso = 0;
   nTracksMuonsMatched_passedIso = 0;
   n2GlobalMuonsMatched_passedIso2Trg = 0;
   nMu0onlyTriggered = 0;
   nMu1onlyTriggered = 0;
   nZMuMu_matched = 0;
   nZMuSta_matched = 0;
   nZMuTrk_matched = 0;
 }

Member Function Documentation

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

overrideprivate

Definition at line 223 of file ZMuMu_MCanalyzer.cc.

References funct::abs(), alpha_, beta_, candidateIsolation(), reco::Candidate::daughter(), dREcal_, dRHcal_, dRTrk_, dRVetoTrk_, edm::View< T >::empty(), reco::Candidate::eta(), etamax_, etamin_, etEcalThreshold_, etHcalThreshold_, GenHFHadronMatcher_cfi::genParticles, genParticlesToken_, h_ProbeOk_eta, h_ProbeOk_pt, h_staProbe_eta, h_staProbe_pt, h_trackProbe_eta, h_trackProbe_pt, hltPath_, mps_fire::i, edm::Ref< C, T, F >::isNonnull(), isoMax_, RazorAnalyzer::lep1, ResonanceBuilder::mass, reco::Candidate::mass(), massMax_, massMin_, reco::Candidate::masterClone(), extraflags_cff::muons, muonsToken_, n2GlobalMuonsMatched_passedIso, n2GlobalMuonsMatched_passedIso2Trg, nGlobalMuonsMatched_passed, nGlobalMuonsMatched_passedIso, nMu0onlyTriggered, nMu1onlyTriggered, nStaMuonsMatched_passedIso, nTracksMuonsMatched_passedIso, nZMuMu_matched, nZMuSta_matched, nZMuTrk_matched, reco::Candidate::pt(), hiDetachedQuadStep_cff::pt1, ptmin_, ptThreshold_, edm::View< T >::refAt(), relativeIsolation_, edm::View< T >::size(), l1t::tracks, tracksToken_, pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), zMuMuMatchMapToken_, zMuMuToken_, zMuStandAloneMatchMapToken_, zMuStandAloneToken_, zMuTrackMatchMapToken_, and zMuTrackToken_.

                                                                           {
   Handle<CandidateView> zMuMu;
   Handle<GenParticleMatch> zMuMuMatchMap; //Map of Z made by Mu global + Mu global
   Handle<CandidateView> zMuStandAlone;
   Handle<GenParticleMatch> zMuStandAloneMatchMap; //Map of Z made by Mu + StandAlone
   Handle<CandidateView> zMuTrack;
   Handle<GenParticleMatch> zMuTrackMatchMap; //Map of Z made by Mu + Track
   Handle<CandidateView> muons; //Collection of Muons
   Handle<CandidateView> tracks; //Collection of Tracks
 
   Handle<GenParticleCollection> genParticles;  // Collection of Generatd Particles
 
   event.getByToken(zMuMuToken_, zMuMu);
   event.getByToken(zMuStandAloneToken_, zMuStandAlone);
   event.getByToken(zMuTrackToken_, zMuTrack);
   event.getByToken(genParticlesToken_, genParticles);
   event.getByToken(muonsToken_, muons);
   event.getByToken(tracksToken_, tracks);
 
   /*
   cout << "*********  zMuMu         size : " << zMuMu->size() << endl;
   cout << "*********  zMuStandAlone size : " << zMuStandAlone->size() << endl;
   cout << "*********  zMuTrack      size : " << zMuTrack->size() << endl;
   cout << "*********  muons         size : " << muons->size() << endl;
   cout << "*********  tracks        size : " << tracks->size() << endl;
   */
   //      std::cout<<"Run-> "<<event.id().run()<<std::endl;
   //      std::cout<<"Event-> "<<event.id().event()<<std::endl;
 
 
   bool zMuMu_found = false;
 
   // loop on ZMuMu
   if (!zMuMu->empty() ) {
     event.getByToken(zMuMuMatchMapToken_, zMuMuMatchMap);
     for(unsigned int i = 0; i < zMuMu->size(); ++i) { //loop on candidates
       const Candidate & zMuMuCand = (*zMuMu)[i]; //the candidate
       CandidateBaseRef zMuMuCandRef = zMuMu->refAt(i);
 
       const Candidate * lep0 = zMuMuCand.daughter( 0 );
       const Candidate * lep1 = zMuMuCand.daughter( 1 );
       const pat::Muon & muonDau0 = dynamic_cast<const pat::Muon &>(*lep0->masterClone());
       //      double trkiso0 = muonDau0.trackIso();
       const pat::Muon & muonDau1 = dynamic_cast<const pat::Muon &>(*lep1->masterClone());
       //double trkiso1 = muonDau1.trackIso();
 
       double iso0 = candidateIsolation(lep0,ptThreshold_, etEcalThreshold_, etHcalThreshold_,dRVetoTrk_, dRTrk_, dREcal_ , dRHcal_,  alpha_, beta_, relativeIsolation_);
       double iso1 = candidateIsolation(lep1,ptThreshold_, etEcalThreshold_, etHcalThreshold_,dRVetoTrk_, dRTrk_, dREcal_ , dRHcal_,  alpha_, beta_, relativeIsolation_);
 
       double pt0 = zMuMuCand.daughter(0)->pt();
       double pt1 = zMuMuCand.daughter(1)->pt();
       double eta0 = zMuMuCand.daughter(0)->eta();
       double eta1 = zMuMuCand.daughter(1)->eta();
       double mass = zMuMuCand.mass();
 
       // HLT match
       const pat::TriggerObjectStandAloneCollection mu0HLTMatches =
     muonDau0.triggerObjectMatchesByPath( hltPath_ );
       const pat::TriggerObjectStandAloneCollection mu1HLTMatches =
     muonDau1.triggerObjectMatchesByPath( hltPath_ );
 
       bool trig0found = false;
       bool trig1found = false;
       if( !mu0HLTMatches.empty() )
     trig0found = true;
       if( !mu1HLTMatches.empty() )
     trig1found = true;
 
       GenParticleRef zMuMuMatch = (*zMuMuMatchMap)[zMuMuCandRef];
       if(zMuMuMatch.isNonnull()) {  // ZMuMu matched
     zMuMu_found = true;
     nZMuMu_matched++;
     if (pt0>ptmin_ && pt1>ptmin_ && abs(eta0)>etamin_ && abs(eta1) >etamin_ && abs(eta0)<etamax_ && abs(eta1) <etamax_ && mass >massMin_ && mass < massMax_ && (trig0found || trig1found)) { // kinematic and trigger cuts passed
       nGlobalMuonsMatched_passed++; // first global Muon passed kine cuts
       nGlobalMuonsMatched_passed++; // second global muon passsed kine cuts
       if (iso0<isoMax_) nGlobalMuonsMatched_passedIso++;       // first global muon passed the iso cut
       if (iso1<isoMax_) nGlobalMuonsMatched_passedIso++;       // second global muon passed the iso cut
       if (iso0<isoMax_ && iso1<isoMax_) {
         n2GlobalMuonsMatched_passedIso++;  // both muons passed iso cut
         if (trig0found && trig1found) n2GlobalMuonsMatched_passedIso2Trg++;  // both muons have HLT
         if (trig0found && !trig1found) nMu0onlyTriggered++;
         if (trig1found && !trig0found) nMu1onlyTriggered++;
         // histograms vs eta and pt
         if (trig1found) {         // check efficiency of muon0 not imposing the trigger on it
           h_trackProbe_eta->Fill(eta0);
           h_trackProbe_pt->Fill(pt0);
           h_staProbe_eta->Fill(eta0);
           h_staProbe_pt->Fill(pt0);
           h_ProbeOk_eta->Fill(eta0);
           h_ProbeOk_pt->Fill(pt0);
         }
         if (trig0found) {         // check efficiency of muon1 not imposing the trigger on it
           h_trackProbe_eta->Fill(eta1);
           h_staProbe_eta->Fill(eta1);
           h_trackProbe_pt->Fill(pt1);
           h_staProbe_pt->Fill(pt1);
           h_ProbeOk_eta->Fill(eta1);
           h_ProbeOk_pt->Fill(pt1);
         }
       }
     }
       } // end MC match
 
     }  // end loop on ZMuMu cand
   }    // end if ZMuMu size > 0
 
   // loop on ZMuSta
   bool zMuSta_found = false;
   if (!zMuMu_found && !zMuStandAlone->empty() ) {
     event.getByToken(zMuStandAloneMatchMapToken_, zMuStandAloneMatchMap);
     for(unsigned int i = 0; i < zMuStandAlone->size(); ++i) { //loop on candidates
       const Candidate & zMuStandAloneCand = (*zMuStandAlone)[i]; //the candidate
       CandidateBaseRef zMuStandAloneCandRef = zMuStandAlone->refAt(i);
       GenParticleRef zMuStandAloneMatch = (*zMuStandAloneMatchMap)[zMuStandAloneCandRef];
 
       const Candidate * lep0 = zMuStandAloneCand.daughter( 0 );
       const Candidate * lep1 = zMuStandAloneCand.daughter( 1 );
       const pat::Muon & muonDau0 = dynamic_cast<const pat::Muon &>(*lep0->masterClone());
       //double trkiso0 = muonDau0.trackIso();
       //      const pat::Muon & muonDau1 = dynamic_cast<const pat::Muon &>(*lep1->masterClone());
       //double trkiso1 = muonDau1.trackIso();
 
       double iso0 = candidateIsolation(lep0,ptThreshold_, etEcalThreshold_, etHcalThreshold_,dRVetoTrk_, dRTrk_, dREcal_ , dRHcal_,  alpha_, beta_, relativeIsolation_);
       double iso1 = candidateIsolation(lep1,ptThreshold_, etEcalThreshold_, etHcalThreshold_,dRVetoTrk_, dRTrk_, dREcal_ , dRHcal_,  alpha_, beta_, relativeIsolation_);
 
       double pt0 = zMuStandAloneCand.daughter(0)->pt();
       double pt1 = zMuStandAloneCand.daughter(1)->pt();
       double eta0 = zMuStandAloneCand.daughter(0)->eta();
       double eta1 = zMuStandAloneCand.daughter(1)->eta();
       double mass = zMuStandAloneCand.mass();
 
       // HLT match (check just dau0 the global)
       const pat::TriggerObjectStandAloneCollection mu0HLTMatches =
     muonDau0.triggerObjectMatchesByPath( hltPath_ );
 
       bool trig0found = false;
       if( !mu0HLTMatches.empty() )
     trig0found = true;
 
       if(zMuStandAloneMatch.isNonnull()) {  // ZMuStandAlone matched
     zMuSta_found = true;
     nZMuSta_matched++;
     if (pt0>ptmin_ && pt1>ptmin_ && abs(eta0)>etamin_ && abs(eta1)>etamin_ && abs(eta0)<etamax_ && abs(eta1) <etamax_ && mass >massMin_ &&
         mass < massMax_ && iso0<isoMax_ && iso1 < isoMax_ && trig0found) { // all cuts and trigger passed
       nStaMuonsMatched_passedIso++;
       // histograms vs eta and pt
       h_staProbe_eta->Fill(eta1);
       h_staProbe_pt->Fill(pt1);
     }
       } // end MC match
     }  // end loop on ZMuStandAlone cand
   }    // end if ZMuStandAlone size > 0
 
 
   // loop on ZMuTrack
   if (!zMuMu_found && !zMuSta_found && !zMuTrack->empty() ) {
     event.getByToken(zMuTrackMatchMapToken_, zMuTrackMatchMap);
     for(unsigned int i = 0; i < zMuTrack->size(); ++i) { //loop on candidates
       const Candidate & zMuTrackCand = (*zMuTrack)[i]; //the candidate
       CandidateBaseRef zMuTrackCandRef = zMuTrack->refAt(i);
       const Candidate * lep0 = zMuTrackCand.daughter( 0 );
       const Candidate * lep1 = zMuTrackCand.daughter( 1 );
       const pat::Muon & muonDau0 = dynamic_cast<const pat::Muon &>(*lep0->masterClone());
       //double trkiso0 = muonDau0.trackIso();
       //const pat::GenericParticle & trackDau1 = dynamic_cast<const pat::GenericParticle &>(*lep1->masterClone());
       //double trkiso1 = trackDau1.trackIso();
       double iso0 = candidateIsolation(lep0,ptThreshold_, etEcalThreshold_, etHcalThreshold_,dRVetoTrk_, dRTrk_, dREcal_ , dRHcal_,  alpha_, beta_, relativeIsolation_);
       double iso1 = candidateIsolation(lep1,ptThreshold_, etEcalThreshold_, etHcalThreshold_,dRVetoTrk_, dRTrk_, dREcal_ , dRHcal_,  alpha_, beta_, relativeIsolation_);
 
 
       double pt0 = zMuTrackCand.daughter(0)->pt();
       double pt1 = zMuTrackCand.daughter(1)->pt();
       double eta0 = zMuTrackCand.daughter(0)->eta();
       double eta1 = zMuTrackCand.daughter(1)->eta();
       double mass = zMuTrackCand.mass();
 
       // HLT match (check just dau0 the global)
       const pat::TriggerObjectStandAloneCollection mu0HLTMatches =
     muonDau0.triggerObjectMatchesByPath( hltPath_ );
 
       bool trig0found = false;
       if( !mu0HLTMatches.empty() )
     trig0found = true;
 
       GenParticleRef zMuTrackMatch = (*zMuTrackMatchMap)[zMuTrackCandRef];
       if(zMuTrackMatch.isNonnull()) {  // ZMuTrack matched
     nZMuTrk_matched++;
     if (pt0>ptmin_ && pt1>ptmin_ && abs(eta0)>etamin_ && abs(eta1)>etamin_ && abs(eta0)<etamax_ && abs(eta1) <etamax_ && mass >massMin_ &&
         mass < massMax_ && iso0<isoMax_ && iso1 < isoMax_ && trig0found) { // all cuts and trigger passed
       nTracksMuonsMatched_passedIso++;
       // histograms vs eta and pt
       h_trackProbe_eta->Fill(eta1);
       h_trackProbe_pt->Fill(pt1);
     }
       }  // end MC match
     }  // end loop on ZMuTrack cand
   }    // end if ZMuTrack size > 0
 
 }       // end analyze

bool ZMuMu_MCanalyzer::check_ifZmumu	(	const Candidate *	dauGen0,
		const Candidate *	dauGen1,
		const Candidate *	dauGen2
	)

private

Definition at line 423 of file ZMuMu_MCanalyzer.cc.

References reco::Candidate::pdgId().

 {
   int partId0 = dauGen0->pdgId();
   int partId1 = dauGen1->pdgId();
   int partId2 = dauGen2->pdgId();
   bool muplusFound=false;
   bool muminusFound=false;
   bool ZFound=false;
   if (partId0==13 || partId1==13 || partId2==13) muminusFound=true;
   if (partId0==-13 || partId1==-13 || partId2==-13) muplusFound=true;
   if (partId0==23 || partId1==23 || partId2==23) ZFound=true;
   return (muplusFound && muminusFound && ZFound);
 }

void ZMuMu_MCanalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 571 of file ZMuMu_MCanalyzer.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, n2GlobalMuonsMatched_passedIso, n2GlobalMuonsMatched_passedIso2Trg, nGlobalMuonsMatched_passed, nGlobalMuonsMatched_passedIso, nMu0onlyTriggered, nMu1onlyTriggered, nStaMuonsMatched_passedIso, nTracksMuonsMatched_passedIso, nZMuMu_matched, nZMuSta_matched, nZMuTrk_matched, and mathSSE::sqrt().

                               {
 
 
   double eff_Iso = double(nGlobalMuonsMatched_passedIso)/nGlobalMuonsMatched_passed;
   double err_effIso = sqrt(eff_Iso*(1-eff_Iso)/nGlobalMuonsMatched_passed);
 
   double n1_afterIso = 2*n2GlobalMuonsMatched_passedIso2Trg+nMu0onlyTriggered+nMu1onlyTriggered+nTracksMuonsMatched_passedIso;
   double n2_afterIso = 2*n2GlobalMuonsMatched_passedIso2Trg+nMu0onlyTriggered+nMu1onlyTriggered+nStaMuonsMatched_passedIso;
   double nGLB_afterIso = 2*n2GlobalMuonsMatched_passedIso2Trg+nMu0onlyTriggered+nMu1onlyTriggered;
   double effSta_afterIso = (2*n2GlobalMuonsMatched_passedIso2Trg+nMu0onlyTriggered+nMu1onlyTriggered)/n1_afterIso;
   double effTrk_afterIso = (2*n2GlobalMuonsMatched_passedIso2Trg+nMu0onlyTriggered+nMu1onlyTriggered)/n2_afterIso;
   double effHLT_afterIso = (2.* n2GlobalMuonsMatched_passedIso2Trg)/(2.* n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered);
   double err_effHLT_afterIso= sqrt( effHLT_afterIso * (1 - effHLT_afterIso)/nGLB_afterIso);
   double err_effsta_afterIso = sqrt(effSta_afterIso*(1-effSta_afterIso)/n1_afterIso);
   double err_efftrk_afterIso = sqrt(effTrk_afterIso*(1-effTrk_afterIso)/n2_afterIso);
 
 
   cout << "------------------------------------  Counters  --------------------------------" << endl;
 
   cout << "number of events zMuMu matched " << nZMuMu_matched << endl;
   cout << "number of events zMuSta matched " << nZMuSta_matched << endl;
   cout << "number of events zMuTk matched " << nZMuTrk_matched << endl;
   cout << "number of events zMuMu with mu0 only triggered " << nMu0onlyTriggered << endl;
   cout << "number of events zMuMu with mu1 only triggered " << nMu1onlyTriggered << endl;
   cout << "=========================================" << endl;
   cout << "n. of global muons MC matched and passing cuts:           " << nGlobalMuonsMatched_passed << endl;
   cout << "n. of global muons MC matched and passing also Iso cut:       " << nGlobalMuonsMatched_passedIso << endl;
   cout << "n. of Z -> 2 global muons MC matched and passing ALL cuts:    " << n2GlobalMuonsMatched_passedIso << endl;
   cout << "n. of ZMuSta MC matched and passing ALL cuts:    " << nStaMuonsMatched_passedIso << endl;
   cout << "n. of ZmuTrck MC matched and passing ALL cuts:   " << nTracksMuonsMatched_passedIso << endl;
   cout << "n. of Z -> 2 global muons MC matched and passing ALL cuts and both triggered: " << n2GlobalMuonsMatched_passedIso2Trg << endl;
   cout << "=================================================================================" << endl;
   cout << "Iso efficiency: " << eff_Iso << " +/- " << err_effIso << endl;
   cout << "HLT efficiency: " << effHLT_afterIso << " +/- " << err_effHLT_afterIso << endl;
   cout << "eff StandAlone (after Isocut) : " << effSta_afterIso << "+/-" << err_effsta_afterIso << endl;
   cout << "eff Tracker (after Isocut)    : " << effTrk_afterIso << "+/-" << err_efftrk_afterIso << endl;
 
 }

float ZMuMu_MCanalyzer::getParticleEta	(	const int	ipart,
		const Candidate *	dauGen0,
		const Candidate *	dauGen1,
		const Candidate *	dauGen2
	)

private

Definition at line 470 of file ZMuMu_MCanalyzer.cc.

References funct::abs(), reco::Candidate::daughter(), reco::Candidate::eta(), gen::k, reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), and reco::Candidate::status().

 {
   int partId0 = dauGen0->pdgId();
   int partId1 = dauGen1->pdgId();
   int partId2 = dauGen2->pdgId();
   float etapart=0.;
   if (partId0 == ipart) {
     for(unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen0->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     etapart = dauMuGen->eta();
       }
     }
   }
   if (partId1 == ipart) {
     for(unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen1->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     etapart = dauMuGen->eta();
       }
     }
   }
   if (partId2 == ipart) {
     for(unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen2->daughter(k);
       if(abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() ==1) {
     etapart = dauMuGen->eta();
       }
     }
   }
   return etapart;
 }

Particle::LorentzVector ZMuMu_MCanalyzer::getParticleP4	(	const int	ipart,
		const Candidate *	dauGen0,
		const Candidate *	dauGen1,
		const Candidate *	dauGen2
	)

private

Definition at line 536 of file ZMuMu_MCanalyzer.cc.

References funct::abs(), reco::Candidate::daughter(), gen::k, reco::Candidate::numberOfDaughters(), reco::Candidate::p4(), reco::Candidate::pdgId(), and reco::Candidate::status().

 {
   int partId0 = dauGen0->pdgId();
   int partId1 = dauGen1->pdgId();
   int partId2 = dauGen2->pdgId();
   Particle::LorentzVector p4part(0.,0.,0.,0.);
   if (partId0 == ipart) {
     for(unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen0->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     p4part = dauMuGen->p4();
       }
     }
   }
   if (partId1 == ipart) {
     for(unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen1->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     p4part = dauMuGen->p4();
       }
     }
   }
   if (partId2 == ipart) {
     for(unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen2->daughter(k);
       if(abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() ==1) {
     p4part = dauMuGen->p4();
       }
     }
   }
   return p4part;
 }

float ZMuMu_MCanalyzer::getParticlePhi	(	const int	ipart,
		const Candidate *	dauGen0,
		const Candidate *	dauGen1,
		const Candidate *	dauGen2
	)

private

Definition at line 503 of file ZMuMu_MCanalyzer.cc.

References funct::abs(), reco::Candidate::daughter(), gen::k, reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::Candidate::phi(), and reco::Candidate::status().

 {
   int partId0 = dauGen0->pdgId();
   int partId1 = dauGen1->pdgId();
   int partId2 = dauGen2->pdgId();
   float phipart=0.;
   if (partId0 == ipart) {
     for(unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen0->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     phipart = dauMuGen->phi();
       }
     }
   }
   if (partId1 == ipart) {
     for(unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen1->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     phipart = dauMuGen->phi();
       }
     }
   }
   if (partId2 == ipart) {
     for(unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen2->daughter(k);
       if(abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() ==1) {
     phipart = dauMuGen->phi();
       }
     }
   }
   return phipart;
 }

float ZMuMu_MCanalyzer::getParticlePt	(	const int	ipart,
		const Candidate *	dauGen0,
		const Candidate *	dauGen1,
		const Candidate *	dauGen2
	)

private

Definition at line 437 of file ZMuMu_MCanalyzer.cc.

References funct::abs(), reco::Candidate::daughter(), gen::k, reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::Candidate::pt(), and reco::Candidate::status().

 {
   int partId0 = dauGen0->pdgId();
   int partId1 = dauGen1->pdgId();
   int partId2 = dauGen2->pdgId();
   float ptpart=0.;
   if (partId0 == ipart) {
     for(unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen0->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     ptpart = dauMuGen->pt();
       }
     }
   }
   if (partId1 == ipart) {
     for(unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen1->daughter(k);
       if(dauMuGen->pdgId() == ipart && dauMuGen->status() ==1) {
     ptpart = dauMuGen->pt();
       }
     }
   }
   if (partId2 == ipart) {
     for(unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
       const Candidate * dauMuGen = dauGen2->daughter(k);
       if(abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() ==1) {
     ptpart = dauMuGen->pt();
       }
     }
   }
   return ptpart;
 }

Member Data Documentation

double ZMuMu_MCanalyzer::alpha_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::beta_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

bool ZMuMu_MCanalyzer::bothMuons_

private

Definition at line 71 of file ZMuMu_MCanalyzer.cc.

double ZMuMu_MCanalyzer::dREcal_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::dRHcal_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::dRTrk_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::dRVetoTrk_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etamax_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etamin_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etEcalThreshold_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etHcalThreshold_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleCollection> ZMuMu_MCanalyzer::genParticlesToken_

private

Definition at line 69 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

reco::CandidateBaseRef ZMuMu_MCanalyzer::globalMuonCandRef_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

TH1D * ZMuMu_MCanalyzer::h_ProbeOk_eta

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_ProbeOk_pt

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_staProbe_eta

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_staProbe_pt

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D* ZMuMu_MCanalyzer::h_trackProbe_eta

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_trackProbe_pt

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

string ZMuMu_MCanalyzer::hltPath_

private

Definition at line 78 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::isoMax_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::massMax_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::massMin_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::muonsToken_

private

Definition at line 67 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

int ZMuMu_MCanalyzer::n2GlobalMuonsMatched_passedIso

private

Definition at line 88 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::n2GlobalMuonsMatched_passedIso2Trg

private

Definition at line 91 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nGlobalMuonsMatched_passed

private

Definition at line 86 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nGlobalMuonsMatched_passedIso

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nMu0onlyTriggered

private

Definition at line 92 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nMu1onlyTriggered

private

Definition at line 93 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nStaMuonsMatched_passedIso

private

Definition at line 89 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nTracksMuonsMatched_passedIso

private

Definition at line 90 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nZMuMu_matched

private

Definition at line 95 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nZMuSta_matched

private

Definition at line 96 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

int ZMuMu_MCanalyzer::nZMuTrk_matched

private

Definition at line 97 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), endJob(), and ZMuMu_MCanalyzer().

OverlapChecker ZMuMu_MCanalyzer::overlap_

private

Definition at line 80 of file ZMuMu_MCanalyzer.cc.

double ZMuMu_MCanalyzer::ptmin_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::ptThreshold_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

bool ZMuMu_MCanalyzer::relativeIsolation_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

reco::CandidateBaseRef ZMuMu_MCanalyzer::standAloneMuonCandRef_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

reco::CandidateBaseRef ZMuMu_MCanalyzer::trackMuonCandRef_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::tracksToken_

private

Definition at line 68 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_MCanalyzer::zMuMuMatchMapToken_

private

Definition at line 62 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::zMuMuToken_

private

Definition at line 61 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_MCanalyzer::zMuStandAloneMatchMapToken_

private

Definition at line 64 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::zMuStandAloneToken_

private

Definition at line 63 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_MCanalyzer::zMuTrackMatchMapToken_

private

Definition at line 66 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::zMuTrackToken_

private

Definition at line 65 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation