ZMuMu_efficiencyAnalyzer Class Reference

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 ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup) override
bool	check_ifZmumu (const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
virtual 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
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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 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 46 of file ZMuMu_efficiencyAnalyzer.cc.

Constructor & Destructor Documentation

ZMuMu_efficiencyAnalyzer::ZMuMu_efficiencyAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 105 of file ZMuMu_efficiencyAnalyzer.cc.

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, i, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), nGlobalMuonsMatched_passed, ptBins, and 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"))),
  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];
    char ap[30], bp[50];
    sprintf(ap,"zmumu2HLTplus_etaRange%d",i);
    sprintf(bp,"zmumu2HLT plus mass eta Range %f to %f",range0,range1);
    cout << ap << "   " << bp << endl;
    hmumu2HLTplus_eta.push_back(etaDirectory.make<TH1D>(ap,bp,200,0.,200.));
    sprintf(ap,"zmumu1HLTplus_etaRange%d",i);
    sprintf(bp,"zmumu1HLT plus mass eta Range %f to %f",range0,range1);
    cout << ap << "   " << bp << endl;
    hmumu1HLTplus_eta.push_back(etaDirectory.make<TH1D>(ap,bp,200,0.,200.));
    sprintf(ap,"zmustaplus_etaRange%d",i);
    sprintf(bp,"zmusta plus mass eta Range %f to %f",range0,range1);
    cout << ap << "   " << bp << endl;
    hmustaplus_eta.push_back(etaDirectory.make<TH1D>(ap,bp,50,0.,200.));
    sprintf(ap,"zmutrackplus_etaRange%d",i);
    sprintf(bp,"zmutrack plus mass eta Range %f to %f",range0,range1);
    cout << ap << "   " << bp << endl;
    hmutrackplus_eta.push_back(etaDirectory.make<TH1D>(ap,bp,100,0.,200.));
    sprintf(ap,"zmumuNotIsoplus_etaRange%d",i);
    sprintf(bp,"zmumuNotIso plus mass eta Range %f to %f",range0,range1);
    cout << ap << "   " << bp << endl;
    hmumuNotIsoplus_eta.push_back(etaDirectory.make<TH1D>(ap,bp,100,0.,200.));
    // muon minus
    cout << " bin eta minus  " << i << endl;
    char am[30], bm[50];
    sprintf(am,"zmumu2HLTminus_etaRange%d",i);
    sprintf(bm,"zmumu2HLT minus mass eta Range %f to %f",range0,range1);
    cout << am << "   " << bm << endl;
    hmumu2HLTminus_eta.push_back(etaDirectory.make<TH1D>(am,bm,200,0.,200.));
    sprintf(am,"zmumu1HLTminus_etaRange%d",i);
    sprintf(bm,"zmumu1HLT minus mass eta Range %f to %f",range0,range1);
    cout << am << "   " << bm << endl;
    hmumu1HLTminus_eta.push_back(etaDirectory.make<TH1D>(am,bm,200,0.,200.));
    sprintf(am,"zmustaminus_etaRange%d",i);
    sprintf(bm,"zmusta minus mass eta Range %f to %f",range0,range1);
    cout << am << "   " << bm << endl;
    hmustaminus_eta.push_back(etaDirectory.make<TH1D>(am,bm,50,0.,200.));
    sprintf(am,"zmutrackminus_etaRange%d",i);
    sprintf(bm,"zmutrack minus mass eta Range %f to %f",range0,range1);
    cout << am << "   " << bm << endl;
    hmutrackminus_eta.push_back(etaDirectory.make<TH1D>(am,bm,100,0.,200.));
    sprintf(am,"zmumuNotIsominus_etaRange%d",i);
    sprintf(bm,"zmumuNotIso minus mass eta Range %f to %f",range0,range1);
    cout << am << "   " << bm << endl;
    hmumuNotIsominus_eta.push_back(etaDirectory.make<TH1D>(am,bm,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;
    char ap1[30], bp1[50];
    sprintf(ap1,"zmumu2HLTplus_ptRange%d",i);
    sprintf(bp1,"zmumu2HLT plus mass pt Range %f to %f",range0,range1);
    cout << ap1 << "   " << bp1 << endl;
    hmumu2HLTplus_pt.push_back(ptDirectory.make<TH1D>(ap1,bp1,200,0.,200.));
    sprintf(ap1,"zmumu1HLTplus_ptRange%d",i);
    sprintf(bp1,"zmumu1HLT plus mass pt Range %f to %f",range0,range1);
    cout << ap1 << "   " << bp1 << endl;
    hmumu1HLTplus_pt.push_back(ptDirectory.make<TH1D>(ap1,bp1,200,0.,200.));
    sprintf(ap1,"zmustaplus_ptRange%d",i);
    sprintf(bp1,"zmusta plus mass pt Range %f to %f",range0,range1);
    cout << ap1 << "   " << bp1 << endl;
    hmustaplus_pt.push_back(ptDirectory.make<TH1D>(ap1,bp1,50,0.,200.));
    sprintf(ap1,"zmutrackplus_ptRange%d",i);
    sprintf(bp1,"zmutrack plus mass pt Range %f to %f",range0,range1);
    cout << ap1 << "   " << bp1 << endl;
    hmutrackplus_pt.push_back(ptDirectory.make<TH1D>(ap1,bp1,100,0.,200.));
    sprintf(ap1,"zmumuNotIsoplus_ptRange%d",i);
    sprintf(bp1,"zmumuNotIso plus mass pt Range %f to %f",range0,range1);
    cout << ap1 << "   " << bp1 << endl;
    hmumuNotIsoplus_pt.push_back(ptDirectory.make<TH1D>(ap1,bp1,100,0.,200.));
    // muon minus
    cout << " bin pt minus  " << i << endl;
    char am1[30], bm1[50];
    sprintf(am1,"zmumu2HLTminus_ptRange%d",i);
    sprintf(bm1,"zmumu2HLT minus mass pt Range %f to %f",range0,range1);
    cout << am1 << "   " << bm1 << endl;
    hmumu2HLTminus_pt.push_back(ptDirectory.make<TH1D>(am1,bm1,200,0.,200.));
    sprintf(am1,"zmumu1HLTminus_ptRange%d",i);
    sprintf(bm1,"zmumu1HLT minus mass pt Range %f to %f",range0,range1);
    cout << am1 << "   " << bm1 << endl;
    hmumu1HLTminus_pt.push_back(ptDirectory.make<TH1D>(am1,bm1,200,0.,200.));
    sprintf(am1,"zmustaminus_ptRange%d",i);
    sprintf(bm1,"zmusta minus mass pt Range %f to %f",range0,range1);
    cout << am1 << "   " << bm1 << endl;
    hmustaminus_pt.push_back(ptDirectory.make<TH1D>(am1,bm1,50,0.,200.));
    sprintf(am1,"zmutrackminus_ptRange%d",i);
    sprintf(bm1,"zmutrack minus mass pt Range %f to %f",range0,range1);
    cout << am1 << "   " << bm1 << endl;
    hmutrackminus_pt.push_back(ptDirectory.make<TH1D>(am1,bm1,100,0.,200.));
    sprintf(am1,"zmumuNotIsominus_ptRange%d",i);
    sprintf(bm1,"zmumuNotIso minus mass pt Range %f to %f",range0,range1);
    cout << am1 << "   " << bm1 << endl;
    hmumuNotIsominus_pt.push_back(ptDirectory.make<TH1D>(am1,bm1,100,0.,200.));
  }

  // clear global counters
  nGlobalMuonsMatched_passed = 0;
}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 264 of file ZMuMu_efficiencyAnalyzer.cc.

References funct::abs(), DeDxDiscriminatorTools::charge(), reco::Candidate::charge(), reco::Candidate::daughter(), reco::Candidate::eta(), eta(), etaBins, etamax_, etaRange, genParticleCandidates2GenParticles_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, i, reco::Muon::isGlobalMuon(), isoMax_, j, reco::Candidate::mass(), massMax_, massMin_, reco::Candidate::masterClone(), patZpeak::muons, muonsToken_, primaryVerticesToken_, RecoTauCleanerPlugins::pt, reco::Candidate::pt(), ptBins, ptmin_, ptRange, pat::Muon::trackIso(), pat::GenericParticle::trackIso(), testEve_cfg::tracks, tracksToken_, pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), ZMuMuAnalysisNtupler_cff::zMuMu, 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
  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->size() > 0 ) {
    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.size()>0 )
    trig0found = true;
      if( mu1HLTMatches.size()>0 )
    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->size() > 0 ) {
    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.size()>0 )
    trig0found = true;
      if( mu1HLTMatches.size()>0 )
    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->size() > 0 ) {
    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.size()>0 )
    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

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

private

Definition at line 591 of file ZMuMu_efficiencyAnalyzer.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_efficiencyAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 739 of file ZMuMu_efficiencyAnalyzer.cc.

                                      {



}

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

private

Definition at line 638 of file ZMuMu_efficiencyAnalyzer.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_efficiencyAnalyzer::getParticleP4 ( const int  ipart, const Candidate *  dauGen0, const Candidate *  dauGen1, const Candidate *  dauGen2  )

private

Definition at line 704 of file ZMuMu_efficiencyAnalyzer.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_efficiencyAnalyzer::getParticlePhi ( const int  ipart, const Candidate *  dauGen0, const Candidate *  dauGen1, const Candidate *  dauGen2  )

private

Definition at line 671 of file ZMuMu_efficiencyAnalyzer.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_efficiencyAnalyzer::getParticlePt ( const int  ipart, const Candidate *  dauGen0, const Candidate *  dauGen1, const Candidate *  dauGen2  )

private

Definition at line 605 of file ZMuMu_efficiencyAnalyzer.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

bool ZMuMu_efficiencyAnalyzer::bothMuons_

private

Definition at line 69 of file ZMuMu_efficiencyAnalyzer.cc.

unsigned int ZMuMu_efficiencyAnalyzer::etaBins

private

Definition at line 74 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

double ZMuMu_efficiencyAnalyzer::etamax_

private

Definition at line 71 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

double ZMuMu_efficiencyAnalyzer::etaRange[7]

private

Definition at line 76 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

EDGetTokenT<GenParticleCollection> ZMuMu_efficiencyAnalyzer::genParticlesToken_

private

Definition at line 66 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

reco::CandidateBaseRef ZMuMu_efficiencyAnalyzer::globalMuonCandRef_

private

Definition at line 79 of file ZMuMu_efficiencyAnalyzer.cc.

TH1D* ZMuMu_efficiencyAnalyzer::h_zmm1HLTminus_mass

private

Definition at line 85 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D* ZMuMu_efficiencyAnalyzer::h_zmm1HLTplus_mass

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmm2HLT_mass

private

Definition at line 83 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D* ZMuMu_efficiencyAnalyzer::h_zmm_mass

private

Definition at line 83 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmmNotIsominus_mass

private

Definition at line 85 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmmNotIsoplus_mass

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmsminus_mass

private

Definition at line 85 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmsplus_mass

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmtminus_mass

private

Definition at line 85 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

TH1D * ZMuMu_efficiencyAnalyzer::h_zmtplus_mass

private

Definition at line 84 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTminus_eta

private

Definition at line 92 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTminus_pt

private

Definition at line 93 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTplus_eta

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu1HLTplus_pt

private

Definition at line 91 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTminus_eta

private

Definition at line 92 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTminus_pt

private

Definition at line 93 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTplus_eta

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumu2HLTplus_pt

private

Definition at line 91 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsominus_eta

private

Definition at line 92 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsominus_pt

private

Definition at line 93 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsoplus_eta

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmumuNotIsoplus_pt

private

Definition at line 91 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaminus_eta

private

Definition at line 92 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaminus_pt

private

Definition at line 93 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaplus_eta

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmustaplus_pt

private

Definition at line 91 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackminus_eta

private

Definition at line 92 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackminus_pt

private

Definition at line 93 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackplus_eta

private

Definition at line 90 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

vector<TH1D *> ZMuMu_efficiencyAnalyzer::hmutrackplus_pt

private

Definition at line 91 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

double ZMuMu_efficiencyAnalyzer::isoMax_

private

Definition at line 71 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

double ZMuMu_efficiencyAnalyzer::massMax_

private

Definition at line 71 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

double ZMuMu_efficiencyAnalyzer::massMin_

private

Definition at line 71 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::muonsToken_

private

Definition at line 64 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

int ZMuMu_efficiencyAnalyzer::nGlobalMuonsMatched_passed

private

Definition at line 88 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by ZMuMu_efficiencyAnalyzer().

OverlapChecker ZMuMu_efficiencyAnalyzer::overlap_

private

Definition at line 80 of file ZMuMu_efficiencyAnalyzer.cc.

EDGetTokenT<VertexCollection> ZMuMu_efficiencyAnalyzer::primaryVerticesToken_

private

Definition at line 67 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

unsigned int ZMuMu_efficiencyAnalyzer::ptBins

private

Definition at line 75 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

double ZMuMu_efficiencyAnalyzer::ptmin_

private

Definition at line 71 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

double ZMuMu_efficiencyAnalyzer::ptRange[5]

private

Definition at line 77 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze(), and ZMuMu_efficiencyAnalyzer().

reco::CandidateBaseRef ZMuMu_efficiencyAnalyzer::standAloneMuonCandRef_

private

Definition at line 79 of file ZMuMu_efficiencyAnalyzer.cc.

reco::CandidateBaseRef ZMuMu_efficiencyAnalyzer::trackMuonCandRef_

private

Definition at line 79 of file ZMuMu_efficiencyAnalyzer.cc.

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::tracksToken_

private

Definition at line 65 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_efficiencyAnalyzer::zMuMuMatchMapToken_

private

Definition at line 59 of file ZMuMu_efficiencyAnalyzer.cc.

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::zMuMuToken_

private

Definition at line 58 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_efficiencyAnalyzer::zMuStandAloneMatchMapToken_

private

Definition at line 61 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::zMuStandAloneToken_

private

Definition at line 60 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_efficiencyAnalyzer::zMuTrackMatchMapToken_

private

Definition at line 63 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_efficiencyAnalyzer::zMuTrackToken_

private

Definition at line 62 of file ZMuMu_efficiencyAnalyzer.cc.

Referenced by analyze().