Inheritance diagram for ZMuMu_efficiencyAnalyzer:

Public Member Functions
	ZMuMu_efficiencyAnalyzer (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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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 &	operator= (EDConsumerBase &&)=default

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

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
bool	bothMuons_

unsigned int	etaBins

double	etamax_

double	etaRange [7]

EDGetTokenT< GenParticleCollection >	genParticlesToken_

reco::CandidateBaseRef	globalMuonCandRef_

TH1D *	h_zmm1HLTminus_mass

TH1D *	h_zmm1HLTplus_mass

TH1D *	h_zmm2HLT_mass

TH1D *	h_zmm_mass

TH1D *	h_zmmNotIsominus_mass

TH1D *	h_zmmNotIsoplus_mass

TH1D *	h_zmsminus_mass

TH1D *	h_zmsplus_mass

TH1D *	h_zmtminus_mass

TH1D *	h_zmtplus_mass

vector< TH1D * >	hmumu1HLTminus_eta

vector< TH1D * >	hmumu1HLTminus_pt

vector< TH1D * >	hmumu1HLTplus_eta

vector< TH1D * >	hmumu1HLTplus_pt

vector< TH1D * >	hmumu2HLTminus_eta

vector< TH1D * >	hmumu2HLTminus_pt

vector< TH1D * >	hmumu2HLTplus_eta

vector< TH1D * >	hmumu2HLTplus_pt

vector< TH1D * >	hmumuNotIsominus_eta

vector< TH1D * >	hmumuNotIsominus_pt

vector< TH1D * >	hmumuNotIsoplus_eta

vector< TH1D * >	hmumuNotIsoplus_pt

vector< TH1D * >	hmustaminus_eta

vector< TH1D * >	hmustaminus_pt

vector< TH1D * >	hmustaplus_eta

vector< TH1D * >	hmustaplus_pt

vector< TH1D * >	hmutrackminus_eta

vector< TH1D * >	hmutrackminus_pt

vector< TH1D * >	hmutrackplus_eta

vector< TH1D * >	hmutrackplus_pt

double	isoMax_

double	massMax_

double	massMin_

EDGetTokenT< CandidateView >	muonsToken_

int	nGlobalMuonsMatched_passed

OverlapChecker	overlap_

EDGetTokenT< VertexCollection >	primaryVerticesToken_

unsigned int	ptBins

double	ptmin_

double	ptRange [5]

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	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Detailed Description

Definition at line 47 of file ZMuMu_efficiencyAnalyzer.cc.

Constructor & Destructor Documentation

◆ ZMuMu_efficiencyAnalyzer()

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

Definition at line 110 of file ZMuMu_efficiencyAnalyzer.cc.

     : 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"))),
       primaryVerticesToken_(consumes<VertexCollection>(pset.getParameter<InputTag>("primaryVertices"))),
  
       bothMuons_(pset.getParameter<bool>("bothMuons")),
  
       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")) {
   Service<TFileService> fs;
  
   // general histograms
   h_zmm_mass = fs->make<TH1D>("zmm_mass", "zmumu mass", 100, 0., 200.);
   h_zmm2HLT_mass = fs->make<TH1D>("zmm2HLT_mass", "zmumu 2HLT mass", 100, 0., 200.);
   h_zmm1HLTplus_mass = fs->make<TH1D>("zmm1HLTplus_mass", "zmumu 1HLT plus mass", 100, 0., 200.);
   h_zmmNotIsoplus_mass = fs->make<TH1D>("zmmNotIsoplus_mass", "zmumu a least One Not Iso plus mass", 100, 0., 200.);
   h_zmsplus_mass = fs->make<TH1D>("zmsplus_mass", "zmusta plus mass", 100, 0., 200.);
   h_zmtplus_mass = fs->make<TH1D>("zmtplus_mass", "zmutrack plus mass", 100, 0., 200.);
   h_zmm1HLTminus_mass = fs->make<TH1D>("zmm1HLTminus_mass", "zmumu 1HLT minus mass", 100, 0., 200.);
   h_zmmNotIsominus_mass = fs->make<TH1D>("zmmNotIsominus_mass", "zmumu a least One Not Iso minus mass", 100, 0., 200.);
   h_zmsminus_mass = fs->make<TH1D>("zmsminus_mass", "zmusta minus mass", 100, 0., 200.);
   h_zmtminus_mass = fs->make<TH1D>("zmtminus_mass", "zmutrack minus mass", 100, 0., 200.);
  
   cout << "primo" << endl;
   // creating histograms for each Pt, eta interval
  
   TFileDirectory etaDirectory =
       fs->mkdir("etaIntervals");  // in this directory will be saved all the histos of different eta intervals
   TFileDirectory ptDirectory =
       fs->mkdir("ptIntervals");  // in this directory will be saved all the histos of different pt intervals
  
   // binning of entries array (at moment defined by hand and not in cfg file)
   etaBins = 6;
   ptBins = 4;
   double etaRangeTmp[7] = {-2., -1.2, -0.8, 0., 0.8, 1.2, 2.};
   double ptRangeTmp[5] = {20., 40., 60., 80., 100.};
   for (unsigned int i = 0; i <= etaBins; i++)
     etaRange[i] = etaRangeTmp[i];
   for (unsigned int i = 0; i <= ptBins; i++)
     ptRange[i] = ptRangeTmp[i];
  
   // eta histograms creation
   cout << "eta istograms creation " << endl;
  
   for (unsigned int i = 0; i < etaBins; i++) {
     cout << " bin eta plus  " << i << endl;
     // muon plus
     double range0 = etaRange[i];
     double range1 = etaRange[i + 1];
     std::string ap, bp;
     ap = "zmumu2HLTplus_etaRange" + std::to_string(i);
     bp = "zmumu2HLT plus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap << "   " << bp << endl;
     hmumu2HLTplus_eta.push_back(etaDirectory.make<TH1D>(ap.c_str(), bp.c_str(), 200, 0., 200.));
     ap = "zmumu1HLTplus_etaRange" + std::to_string(i);
     bp = "zmumu1HLT plus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap << "   " << bp << endl;
     hmumu1HLTplus_eta.push_back(etaDirectory.make<TH1D>(ap.c_str(), bp.c_str(), 200, 0., 200.));
     ap = "zmustaplus_etaRange" + std::to_string(i);
     bp = "zmusta plus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap << "   " << bp << endl;
     hmustaplus_eta.push_back(etaDirectory.make<TH1D>(ap.c_str(), bp.c_str(), 50, 0., 200.));
     ap = "zmutrackplus_etaRange" + std::to_string(i);
     bp = "zmutrack plus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap << "   " << bp << endl;
     hmutrackplus_eta.push_back(etaDirectory.make<TH1D>(ap.c_str(), bp.c_str(), 100, 0., 200.));
     ap = "zmumuNotIsoplus_etaRange" + std::to_string(i);
     bp = "zmumuNotIso plus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap << "   " << bp << endl;
     hmumuNotIsoplus_eta.push_back(etaDirectory.make<TH1D>(ap.c_str(), bp.c_str(), 100, 0., 200.));
     // muon minus
     cout << " bin eta minus  " << i << endl;
     std::string am, bm;
     am = "zmumu2HLTminus_etaRange" + std::to_string(i);
     bm = "zmumu2HLT minus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am << "   " << bm << endl;
     hmumu2HLTminus_eta.push_back(etaDirectory.make<TH1D>(am.c_str(), bm.c_str(), 200, 0., 200.));
     am = "zmumu1HLTminus_etaRange" + std::to_string(i);
     bm = "zmumu1HLT minus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am << "   " << bm << endl;
     hmumu1HLTminus_eta.push_back(etaDirectory.make<TH1D>(am.c_str(), bm.c_str(), 200, 0., 200.));
     am = "zmustaminus_etaRange" + std::to_string(i);
     bm = "zmusta minus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am << "   " << bm << endl;
     hmustaminus_eta.push_back(etaDirectory.make<TH1D>(am.c_str(), bm.c_str(), 50, 0., 200.));
     am = "zmutrackminus_etaRange" + std::to_string(i);
     bm = "zmutrack minus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am << "   " << bm << endl;
     hmutrackminus_eta.push_back(etaDirectory.make<TH1D>(am.c_str(), bm.c_str(), 100, 0., 200.));
     am = "zmumuNotIsominus_etaRange" + std::to_string(i);
     bm = "zmumuNotIso minus mass eta Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am << "   " << bm << endl;
     hmumuNotIsominus_eta.push_back(etaDirectory.make<TH1D>(am.c_str(), bm.c_str(), 100, 0., 200.));
   }
  
   // pt histograms creation
   cout << "pt istograms creation " << endl;
  
   for (unsigned int i = 0; i < ptBins; i++) {
     double range0 = ptRange[i];
     double range1 = ptRange[i + 1];
     // muon plus
     cout << " bin pt plus  " << i << endl;
     std::string ap1, bp1;
     ap1 = "zmumu2HLTplus_ptRange" + std::to_string(i);
     bp1 = "zmumu2HLT plus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap1 << "   " << bp1 << endl;
     hmumu2HLTplus_pt.push_back(ptDirectory.make<TH1D>(ap1.c_str(), bp1.c_str(), 200, 0., 200.));
     ap1 = "zmumu1HLTplus_ptRange" + std::to_string(i);
     bp1 = "zmumu1HLT plus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap1 << "   " << bp1 << endl;
     hmumu1HLTplus_pt.push_back(ptDirectory.make<TH1D>(ap1.c_str(), bp1.c_str(), 200, 0., 200.));
     ap1 = "zmustaplus_ptRange" + std::to_string(i);
     bp1 = "zmusta plus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap1 << "   " << bp1 << endl;
     hmustaplus_pt.push_back(ptDirectory.make<TH1D>(ap1.c_str(), bp1.c_str(), 50, 0., 200.));
     ap1 = "zmutrackplus_ptRange" + std::to_string(i);
     bp1 = "zmutrack plus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap1 << "   " << bp1 << endl;
     hmutrackplus_pt.push_back(ptDirectory.make<TH1D>(ap1.c_str(), bp1.c_str(), 100, 0., 200.));
     ap1 = "zmumuNotIsoplus_ptRange" + std::to_string(i);
     bp1 = "zmumuNotIso plus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << ap1 << "   " << bp1 << endl;
     hmumuNotIsoplus_pt.push_back(ptDirectory.make<TH1D>(ap1.c_str(), bp1.c_str(), 100, 0., 200.));
     // muon minus
     cout << " bin pt minus  " << i << endl;
     std::string am1, bm1;
     am1 = "zmumu2HLTminus_ptRange" + std::to_string(i);
     bm1 = "zmumu2HLT minus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am1 << "   " << bm1 << endl;
     hmumu2HLTminus_pt.push_back(ptDirectory.make<TH1D>(am1.c_str(), bm1.c_str(), 200, 0., 200.));
     am1 = "zmumu1HLTminus_ptRange" + std::to_string(i);
     bm1 = "zmumu1HLT minus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am1 << "   " << bm1 << endl;
     hmumu1HLTminus_pt.push_back(ptDirectory.make<TH1D>(am1.c_str(), bm1.c_str(), 200, 0., 200.));
     am1 = "zmustaminus_ptRange" + std::to_string(i);
     bm1 = "zmusta minus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am1 << "   " << bm1 << endl;
     hmustaminus_pt.push_back(ptDirectory.make<TH1D>(am1.c_str(), bm1.c_str(), 50, 0., 200.));
     am1 = "zmutrackminus_ptRange" + std::to_string(i);
     bm1 = "zmutrack minus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am1 << "   " << bm1 << endl;
     hmutrackminus_pt.push_back(ptDirectory.make<TH1D>(am1.c_str(), bm1.c_str(), 100, 0., 200.));
     am1 = "zmumuNotIsominus_ptRange" + std::to_string(i);
     bm1 = "zmumuNotIso minus mass pt Range " + std::to_string(range0) + " to " + std::to_string(range1);
     cout << am1 << "   " << bm1 << endl;
     hmumuNotIsominus_pt.push_back(ptDirectory.make<TH1D>(am1.c_str(), bm1.c_str(), 100, 0., 200.));
   }
  
   // clear global counters
   nGlobalMuonsMatched_passed = 0;
 }

References gather_cfg::cout, etaBins, etaRange, h_zmm1HLTminus_mass, h_zmm1HLTplus_mass, h_zmm2HLT_mass, h_zmm_mass, h_zmmNotIsominus_mass, h_zmmNotIsoplus_mass, h_zmsminus_mass, h_zmsplus_mass, h_zmtminus_mass, h_zmtplus_mass, hmumu1HLTminus_eta, hmumu1HLTminus_pt, hmumu1HLTplus_eta, hmumu1HLTplus_pt, hmumu2HLTminus_eta, hmumu2HLTminus_pt, hmumu2HLTplus_eta, hmumu2HLTplus_pt, hmumuNotIsominus_eta, hmumuNotIsominus_pt, hmumuNotIsoplus_eta, hmumuNotIsoplus_pt, hmustaminus_eta, hmustaminus_pt, hmustaplus_eta, hmustaplus_pt, hmutrackminus_eta, hmutrackminus_pt, hmutrackplus_eta, hmutrackplus_pt, mps_fire::i, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), nGlobalMuonsMatched_passed, ptBins, ptRange, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 273 of file ZMuMu_efficiencyAnalyzer.cc.

                                                                                   {
   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
   Handle<VertexCollection> primaryVertices;    // Collection of primary Vertices
  
   event.getByToken(zMuMuToken_, zMuMu);
   event.getByToken(zMuStandAloneToken_, zMuStandAlone);
   event.getByToken(zMuTrackToken_, zMuTrack);
   event.getByToken(genParticlesToken_, genParticles);
   event.getByToken(primaryVerticesToken_, primaryVertices);
   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;
   cout << "*********  vertices      size : " << primaryVertices->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()) {
     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();
  
       // kinemtic variables
       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 charge0 = zMuMuCand.daughter(0)->charge();
       double charge1 = zMuMuCand.daughter(1)->charge();
       double mass = zMuMuCand.mass();
  
       // HLT match
       const pat::TriggerObjectStandAloneCollection mu0HLTMatches = muonDau0.triggerObjectMatchesByPath("HLT_Mu9");
       const pat::TriggerObjectStandAloneCollection mu1HLTMatches = muonDau1.triggerObjectMatchesByPath("HLT_Mu9");
  
       bool trig0found = false;
       bool trig1found = false;
       if (!mu0HLTMatches.empty())
         trig0found = true;
       if (!mu1HLTMatches.empty())
         trig1found = true;
  
       // kinematic selection
  
       bool checkOppositeCharge = false;
       if (charge0 != charge1)
         checkOppositeCharge = true;
       if (pt0 > ptmin_ && pt1 > ptmin_ && abs(eta0) < etamax_ && abs(eta1) < etamax_ && mass > massMin_ &&
           mass < massMax_ && checkOppositeCharge) {
         if (trig0found || trig1found) {                  // at least one muon match HLT
           zMuMu_found = true;                            // Z found as global-global (so don't check Zms and Zmt)
           if (trkiso0 < isoMax_ && trkiso1 < isoMax_) {  // both muons are isolated
             if (trig0found && trig1found) {
               // ******************** category zmm 2 HLT ****************
  
               h_zmm2HLT_mass->Fill(mass);
               h_zmm_mass->Fill(mass);
  
               // check the cynematics to fill correct histograms
  
               for (unsigned int j = 0; j < etaBins; j++) {  // eta Bins loop
                 double range0 = etaRange[j];
                 double range1 = etaRange[j + 1];
  
                 // eta histograms
  
                 if (eta0 >= range0 && eta0 < range1) {
                   if (charge0 < 0)
                     hmumu2HLTminus_eta[j]->Fill(mass);  // mu- in bin eta
                   if (charge0 > 0)
                     hmumu2HLTplus_eta[j]->Fill(mass);  // mu+ in bin eta
                 }
                 if (eta1 >= range0 && eta1 < range1) {
                   if (charge1 < 0)
                     hmumu2HLTminus_eta[j]->Fill(mass);  // mu- in bin eta
                   if (charge1 > 0)
                     hmumu2HLTplus_eta[j]->Fill(mass);  // mu+ in bin eta
                 }
               }  // end loop etaBins
  
               for (unsigned int j = 0; j < ptBins; j++) {  // pt Bins loop
                 double range0pt = ptRange[j];
                 double range1pt = ptRange[j + 1];
                 // pt histograms
                 if (pt0 >= range0pt && pt0 < range1pt) {
                   if (charge0 < 0)
                     hmumu2HLTminus_pt[j]->Fill(mass);  // mu- in bin eta
                   if (charge0 > 0)
                     hmumu2HLTplus_pt[j]->Fill(mass);  // mu+ in bin eta
                 }
                 if (pt1 >= range0pt && pt1 < range1pt) {
                   if (charge1 < 0)
                     hmumu2HLTminus_pt[j]->Fill(mass);  // mu- in bin eta
                   if (charge1 > 0)
                     hmumu2HLTplus_pt[j]->Fill(mass);  // mu+ in bin eta
                 }
               }  // end loop  ptBins
  
             }  // ******************* end category zmm 2 HLT ****************
  
             if (!trig0found || !trig1found) {
               // ****************** category zmm 1 HLT ******************
               h_zmm_mass->Fill(mass);
               double eta = 9999;
               double pt = 9999;
               double charge = 0;
               if (trig0found) {
                 eta = eta1;  // check  muon not HLT matched
                 pt = pt1;
                 charge = charge1;
               } else {
                 eta = eta0;
                 pt = pt0;
                 charge = charge0;
               }
               if (charge < 0)
                 h_zmm1HLTminus_mass->Fill(mass);
               if (charge > 0)
                 h_zmm1HLTplus_mass->Fill(mass);
  
               for (unsigned int j = 0; j < etaBins; j++) {  // eta Bins loop
                 double range0 = etaRange[j];
                 double range1 = etaRange[j + 1];
                 // eta histograms fill the bin of the muon not HLT matched
                 if (eta >= range0 && eta < range1) {
                   if (charge < 0)
                     hmumu1HLTminus_eta[j]->Fill(mass);
                   if (charge > 0)
                     hmumu1HLTplus_eta[j]->Fill(mass);
                 }
               }                                            // end loop etaBins
               for (unsigned int j = 0; j < ptBins; j++) {  // pt Bins loop
                 double range0 = ptRange[j];
                 double range1 = ptRange[j + 1];
                 // pt histograms
                 if (pt >= range0 && pt < range1) {
                   if (charge < 0)
                     hmumu1HLTminus_pt[j]->Fill(mass);
                   if (charge > 0)
                     hmumu1HLTplus_pt[j]->Fill(mass);
                 }
               }  // end loop ptBins
  
             }  // ****************** end category zmm 1 HLT ***************
  
           } else {  // one or both muons are not isolated
             // *****************  category zmumuNotIso **************** (per ora non studio iso vs eta e pt da capire meglio)
  
           }  // end if both muons isolated
  
         }  // end if at least 1 HLT trigger found
       }    // end if kinematic selection
  
     }  // 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 pt0 = zMuStandAloneCand.daughter(0)->pt();
       double pt1 = zMuStandAloneCand.daughter(1)->pt();
       double eta0 = zMuStandAloneCand.daughter(0)->eta();
       double eta1 = zMuStandAloneCand.daughter(1)->eta();
       double charge0 = zMuStandAloneCand.daughter(0)->charge();
       double charge1 = zMuStandAloneCand.daughter(1)->charge();
       double mass = zMuStandAloneCand.mass();
  
       // HLT match
       const pat::TriggerObjectStandAloneCollection mu0HLTMatches = muonDau0.triggerObjectMatchesByPath("HLT_Mu9");
       const pat::TriggerObjectStandAloneCollection mu1HLTMatches = muonDau1.triggerObjectMatchesByPath("HLT_Mu9");
  
       bool trig0found = false;
       bool trig1found = false;
       if (!mu0HLTMatches.empty())
         trig0found = true;
       if (!mu1HLTMatches.empty())
         trig1found = true;
  
       // check HLT match of Global muon and save eta, pt of second muon (standAlone)
       bool trigGlbfound = false;
       double pt = 999.;
       double eta = 999.;
       double charge = 0;
       if (muonDau0.isGlobalMuon()) {
         trigGlbfound = trig0found;
         pt = pt1;
         eta = eta1;
         charge = charge1;
       }
       if (muonDau1.isGlobalMuon()) {
         trigGlbfound = trig1found;
         pt = pt0;
         eta = eta0;
         charge = charge0;
       }
  
       bool checkOppositeCharge = false;
       if (charge0 != charge1)
         checkOppositeCharge = true;
  
       if (checkOppositeCharge && trigGlbfound && pt0 > ptmin_ && pt1 > ptmin_ && abs(eta0) < etamax_ &&
           abs(eta1) < etamax_ && mass > massMin_ && mass < massMax_ && trkiso0 < isoMax_ &&
           trkiso1 < isoMax_) {  // global mu match HLT + kinematic cuts + opposite charge
  
         if (charge < 0)
           h_zmsminus_mass->Fill(mass);
         if (charge > 0)
           h_zmsplus_mass->Fill(mass);
  
         for (unsigned int j = 0; j < etaBins; j++) {  // eta Bins loop
           double range0 = etaRange[j];
           double range1 = etaRange[j + 1];
           // eta histograms
           if (eta >= range0 && eta < range1) {
             if (charge < 0)
               hmustaminus_eta[j]->Fill(mass);
             if (charge > 0)
               hmustaplus_eta[j]->Fill(mass);
           }
         }                                            // end loop etaBins
         for (unsigned int j = 0; j < ptBins; j++) {  // pt Bins loop
           double range0 = ptRange[j];
           double range1 = ptRange[j + 1];
           // pt histograms
           if (pt >= range0 && pt < range1) {
             if (charge < 0)
               hmustaminus_pt[j]->Fill(mass);
             if (charge > 0)
               hmustaplus_pt[j]->Fill(mass);
           }
         }  // end loop ptBins
  
       }  // end if trigGlbfound + kinecuts + OppostieCharge
     }    // end loop on ZMuStandAlone cand
   }      // end if ZMuStandAlone size > 0
  
   // loop on ZMuTrack
   //  bool zMuTrack_found = false;
   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 pt0 = zMuTrackCand.daughter(0)->pt();
       double pt1 = zMuTrackCand.daughter(1)->pt();
       double eta0 = zMuTrackCand.daughter(0)->eta();
       double eta1 = zMuTrackCand.daughter(1)->eta();
       double charge0 = zMuTrackCand.daughter(0)->charge();
       double charge1 = zMuTrackCand.daughter(1)->charge();
       double mass = zMuTrackCand.mass();
  
       // HLT match (check just dau0 the global)
       const pat::TriggerObjectStandAloneCollection mu0HLTMatches = muonDau0.triggerObjectMatchesByPath("HLT_Mu9");
  
       bool trig0found = false;
       if (!mu0HLTMatches.empty())
         trig0found = true;
  
       bool checkOppositeCharge = false;
       if (charge0 != charge1)
         checkOppositeCharge = true;
  
       if (checkOppositeCharge && trig0found && pt0 > ptmin_ && pt1 > ptmin_ && abs(eta0) < etamax_ &&
           abs(eta1) < etamax_ && mass > massMin_ && mass < massMax_ && trkiso0 < isoMax_ &&
           trkiso1 < isoMax_) {  // global mu match HLT + kinematic cuts + opposite charge
  
         if (charge1 < 0)
           h_zmtminus_mass->Fill(mass);
         if (charge1 > 0)
           h_zmtplus_mass->Fill(mass);
  
         for (unsigned int j = 0; j < etaBins; j++) {  // eta Bins loop
           double range0 = etaRange[j];
           double range1 = etaRange[j + 1];
           // eta histograms
           if (eta1 >= range0 && eta1 < range1) {
             if (charge1 < 0)
               hmutrackminus_eta[j]->Fill(mass);  // just check muon1 (mu0 is global by definition)
             if (charge1 > 0)
               hmutrackplus_eta[j]->Fill(mass);  // just check muon1 (mu0 is global by definition)
           }
         }                                            // end loop etaBins
         for (unsigned int j = 0; j < ptBins; j++) {  // pt Bins loop
           double range0 = ptRange[j];
           double range1 = ptRange[j + 1];
           // pt histograms
           if (pt1 >= range0 && pt1 < range1) {
             if (charge1 < 0)
               hmutrackminus_pt[j]->Fill(mass);  // just check muon1 (mu0 is global by definition)
             if (charge1 > 0)
               hmutrackplus_pt[j]->Fill(mass);  // just check muon1 (mu0 is global by definition)
           }
         }  // end loop ptBins
  
       }  // end if trig0found
  
     }  // end loop on ZMuTrack cand
   }    // end if ZMuTrack size > 0
  
 }  // end analyze

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, reco::Candidate::charge(), reco::Candidate::daughter(), edm::View< T >::empty(), PVValHelper::eta, reco::Candidate::eta(), HLT_FULL_cff::eta1, etaBins, etamax_, etaRange, genParticles2HepMC_cfi::genParticles, genParticlesToken_, h_zmm1HLTminus_mass, h_zmm1HLTplus_mass, h_zmm2HLT_mass, h_zmm_mass, h_zmsminus_mass, h_zmsplus_mass, h_zmtminus_mass, h_zmtplus_mass, hmumu1HLTminus_eta, hmumu1HLTminus_pt, hmumu1HLTplus_eta, hmumu1HLTplus_pt, hmumu2HLTminus_eta, hmumu2HLTminus_pt, hmumu2HLTplus_eta, hmumu2HLTplus_pt, hmustaminus_eta, hmustaminus_pt, hmustaplus_eta, hmustaplus_pt, hmutrackminus_eta, hmutrackminus_pt, hmutrackplus_eta, hmutrackplus_pt, mps_fire::i, reco::Muon::isGlobalMuon(), isoMax_, dqmiolumiharvest::j, RazorAnalyzer::lep1, EgHLTOffHistBins_cfi::mass, reco::Candidate::mass(), massMax_, massMin_, reco::Candidate::masterClone(), PDWG_BPHSkim_cff::muons, muonsToken_, zMuMuMuonUserData::primaryVertices, primaryVerticesToken_, DiDispStaMuonMonitor_cfi::pt, reco::Candidate::pt(), HLT_FULL_cff::pt1, ptBins, ptmin_, ptRange, edm::View< T >::refAt(), edm::View< T >::size(), pat::Muon::trackIso(), pat::GenericParticle::trackIso(), PDWG_EXOHSCP_cff::tracks, tracksToken_, pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), ZMuMuAnalysisNtupler_cff::zMuMu, zMuMuToken_, zMuStandAloneMatchMapToken_, zMuStandAloneToken_, zMuTrackMatchMapToken_, and zMuTrackToken_.

◆ check_ifZmumu()

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

private

Definition at line 617 of file ZMuMu_efficiencyAnalyzer.cc.

                                                                        {
   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);
 }

References reco::Candidate::pdgId().

◆ endJob()

void ZMuMu_efficiencyAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 775 of file ZMuMu_efficiencyAnalyzer.cc.

775 {}

◆ getParticleEta()

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

private

Definition at line 670 of file ZMuMu_efficiencyAnalyzer.cc.

                                                                          {
   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;
 }

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

◆ getParticleP4()

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

private

Definition at line 740 of file ZMuMu_efficiencyAnalyzer.cc.

                                                                                           {
   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;
 }

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

◆ getParticlePhi()

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

private

Definition at line 705 of file ZMuMu_efficiencyAnalyzer.cc.

                                                                          {
   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;
 }

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

◆ getParticlePt()

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

private

Definition at line 635 of file ZMuMu_efficiencyAnalyzer.cc.

                                                                         {
   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;
 }

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

Member Data Documentation

◆ bothMuons_

bool ZMuMu_efficiencyAnalyzer::bothMuons_

private

Definition at line 74 of file ZMuMu_efficiencyAnalyzer.cc.

◆ etaBins

unsigned int ZMuMu_efficiencyAnalyzer::etaBins

private

Definition at line 79 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), QGLikelihoodCalculator.QGLikelihoodCalculator::findEntry(), QGLikelihoodCalculator.QGLikelihoodCalculator::init(), QGLikelihoodCalculator.QGLikelihoodCalculator::isValidRange(), and ZMuMu_efficiencyAnalyzer().

◆ etamax_

double ZMuMu_efficiencyAnalyzer::etamax_

private

Definition at line 76 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ etaRange

double ZMuMu_efficiencyAnalyzer::etaRange[7]

private

Definition at line 81 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ genParticlesToken_

EDGetTokenT<GenParticleCollection> ZMuMu_efficiencyAnalyzer::genParticlesToken_

private

Definition at line 71 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ globalMuonCandRef_

reco::CandidateBaseRef ZMuMu_efficiencyAnalyzer::globalMuonCandRef_

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

◆ h_zmm1HLTminus_mass

TH1D* ZMuMu_efficiencyAnalyzer::h_zmm1HLTminus_mass

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmm1HLTplus_mass

TH1D* ZMuMu_efficiencyAnalyzer::h_zmm1HLTplus_mass

private

Definition at line 89 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmm2HLT_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmm2HLT_mass

private

Definition at line 88 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmm_mass

TH1D* ZMuMu_efficiencyAnalyzer::h_zmm_mass

private

Definition at line 88 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmmNotIsominus_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmmNotIsominus_mass

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ h_zmmNotIsoplus_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmmNotIsoplus_mass

private

Definition at line 89 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ h_zmsminus_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmsminus_mass

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmsplus_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmsplus_mass

private

Definition at line 89 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmtminus_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmtminus_mass

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ h_zmtplus_mass

TH1D * ZMuMu_efficiencyAnalyzer::h_zmtplus_mass

private

Definition at line 89 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu1HLTminus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTminus_eta

private

Definition at line 97 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu1HLTminus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTminus_pt

private

Definition at line 98 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu1HLTplus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTplus_eta

private

Definition at line 95 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu1HLTplus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTplus_pt

private

Definition at line 96 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu2HLTminus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTminus_eta

private

Definition at line 97 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu2HLTminus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTminus_pt

private

Definition at line 98 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu2HLTplus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTplus_eta

private

Definition at line 95 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumu2HLTplus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTplus_pt

private

Definition at line 96 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmumuNotIsominus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsominus_eta

private

Definition at line 97 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ hmumuNotIsominus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsominus_pt

private

Definition at line 98 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ hmumuNotIsoplus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsoplus_eta

private

Definition at line 95 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ hmumuNotIsoplus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsoplus_pt

private

Definition at line 96 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ hmustaminus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaminus_eta

private

Definition at line 97 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmustaminus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaminus_pt

private

Definition at line 98 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmustaplus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaplus_eta

private

Definition at line 95 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmustaplus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaplus_pt

private

Definition at line 96 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmutrackminus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackminus_eta

private

Definition at line 97 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmutrackminus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackminus_pt

private

Definition at line 98 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmutrackplus_eta

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackplus_eta

private

Definition at line 95 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ hmutrackplus_pt

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackplus_pt

private

Definition at line 96 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ isoMax_

double ZMuMu_efficiencyAnalyzer::isoMax_

private

Definition at line 76 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ massMax_

double ZMuMu_efficiencyAnalyzer::massMax_

private

Definition at line 76 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ massMin_

double ZMuMu_efficiencyAnalyzer::massMin_

private

Definition at line 76 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ muonsToken_

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::muonsToken_

private

Definition at line 69 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ nGlobalMuonsMatched_passed

int ZMuMu_efficiencyAnalyzer::nGlobalMuonsMatched_passed

private

Definition at line 93 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

◆ overlap_

OverlapChecker ZMuMu_efficiencyAnalyzer::overlap_

private

Definition at line 85 of file ZMuMu_efficiencyAnalyzer.cc.

◆ primaryVerticesToken_

EDGetTokenT<VertexCollection> ZMuMu_efficiencyAnalyzer::primaryVerticesToken_

private

Definition at line 72 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ ptBins

unsigned int ZMuMu_efficiencyAnalyzer::ptBins

private

Definition at line 80 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), QGLikelihoodCalculator.QGLikelihoodCalculator::findEntry(), QGLikelihoodCalculator.QGLikelihoodCalculator::init(), QGLikelihoodCalculator.QGLikelihoodCalculator::isValidRange(), and ZMuMu_efficiencyAnalyzer().

◆ ptmin_

double ZMuMu_efficiencyAnalyzer::ptmin_

private

Definition at line 76 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ ptRange

double ZMuMu_efficiencyAnalyzer::ptRange[5]

private

Definition at line 82 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

◆ standAloneMuonCandRef_

reco::CandidateBaseRef ZMuMu_efficiencyAnalyzer::standAloneMuonCandRef_

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

◆ trackMuonCandRef_

reco::CandidateBaseRef ZMuMu_efficiencyAnalyzer::trackMuonCandRef_

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

◆ tracksToken_

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::tracksToken_

private

Definition at line 70 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ zMuMuMatchMapToken_

EDGetTokenT<GenParticleMatch> ZMuMu_efficiencyAnalyzer::zMuMuMatchMapToken_

private

Definition at line 64 of file ZMuMu_efficiencyAnalyzer.cc.

◆ zMuMuToken_

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::zMuMuToken_

private

Definition at line 63 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ zMuStandAloneMatchMapToken_

EDGetTokenT<GenParticleMatch> ZMuMu_efficiencyAnalyzer::zMuStandAloneMatchMapToken_

private

Definition at line 66 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ zMuStandAloneToken_

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::zMuStandAloneToken_

private

Definition at line 65 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ zMuTrackMatchMapToken_

EDGetTokenT<GenParticleMatch> ZMuMu_efficiencyAnalyzer::zMuTrackMatchMapToken_

private

Definition at line 68 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

◆ zMuTrackToken_

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::zMuTrackToken_

private

Definition at line 67 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ZMuMu_efficiencyAnalyzer()

Member Function Documentation

◆ analyze()

◆ check_ifZmumu()

◆ endJob()

◆ getParticleEta()

◆ getParticleP4()

◆ getParticlePhi()

◆ getParticlePt()

Member Data Documentation

◆ bothMuons_

◆ etaBins

◆ etamax_

◆ etaRange

◆ genParticlesToken_

◆ globalMuonCandRef_

◆ h_zmm1HLTminus_mass

◆ h_zmm1HLTplus_mass

◆ h_zmm2HLT_mass

◆ h_zmm_mass

◆ h_zmmNotIsominus_mass

◆ h_zmmNotIsoplus_mass

◆ h_zmsminus_mass

◆ h_zmsplus_mass

◆ h_zmtminus_mass

◆ h_zmtplus_mass

◆ hmumu1HLTminus_eta

◆ hmumu1HLTminus_pt

◆ hmumu1HLTplus_eta

◆ hmumu1HLTplus_pt

◆ hmumu2HLTminus_eta

◆ hmumu2HLTminus_pt

◆ hmumu2HLTplus_eta

◆ hmumu2HLTplus_pt

◆ hmumuNotIsominus_eta

◆ hmumuNotIsominus_pt

◆ hmumuNotIsoplus_eta

◆ hmumuNotIsoplus_pt

◆ hmustaminus_eta

◆ hmustaminus_pt

◆ hmustaplus_eta

◆ hmustaplus_pt

◆ hmutrackminus_eta

◆ hmutrackminus_pt

◆ hmutrackplus_eta

◆ hmutrackplus_pt

◆ isoMax_

◆ massMax_

◆ massMin_

◆ muonsToken_

◆ nGlobalMuonsMatched_passed

◆ overlap_

◆ primaryVerticesToken_

◆ ptBins

◆ ptmin_

◆ ptRange

◆ standAloneMuonCandRef_

◆ trackMuonCandRef_

◆ tracksToken_

◆ zMuMuMatchMapToken_

◆ zMuMuToken_

◆ zMuStandAloneMatchMapToken_

◆ zMuStandAloneToken_

◆ zMuTrackMatchMapToken_

◆ zMuTrackToken_