ZMuMu_MCanalyzer Class Reference

Inheritance diagram for ZMuMu_MCanalyzer:

Public Member Functions
	ZMuMu_MCanalyzer (const edm::ParameterSet &pset)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override
bool	check_ifZmumu (const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
void	endJob () override
float	getParticleEta (const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
Particle::LorentzVector	getParticleP4 (const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
float	getParticlePhi (const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
float	getParticlePt (const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)

Private Attributes
double	alpha_
double	beta_
bool	bothMuons_
double	dREcal_
double	dRHcal_
double	dRTrk_
double	dRVetoTrk_
double	etamax_
double	etamin_
double	etEcalThreshold_
double	etHcalThreshold_
EDGetTokenT< GenParticleCollection >	genParticlesToken_
reco::CandidateBaseRef	globalMuonCandRef_
TH1D *	h_ProbeOk_eta
TH1D *	h_ProbeOk_pt
TH1D *	h_staProbe_eta
TH1D *	h_staProbe_pt
TH1D *	h_trackProbe_eta
TH1D *	h_trackProbe_pt
string	hltPath_
double	isoMax_
double	massMax_
double	massMin_
EDGetTokenT< CandidateView >	muonsToken_
int	n2GlobalMuonsMatched_passedIso
int	n2GlobalMuonsMatched_passedIso2Trg
int	nGlobalMuonsMatched_passed
int	nGlobalMuonsMatched_passedIso
int	nMu0onlyTriggered
int	nMu1onlyTriggered
int	nStaMuonsMatched_passedIso
int	nTracksMuonsMatched_passedIso
int	nZMuMu_matched
int	nZMuSta_matched
int	nZMuTrk_matched
OverlapChecker	overlap_
double	ptmin_
double	ptThreshold_
bool	relativeIsolation_
reco::CandidateBaseRef	standAloneMuonCandRef_
reco::CandidateBaseRef	trackMuonCandRef_
EDGetTokenT< CandidateView >	tracksToken_
EDGetTokenT< GenParticleMatch >	zMuMuMatchMapToken_
EDGetTokenT< CandidateView >	zMuMuToken_
EDGetTokenT< GenParticleMatch >	zMuStandAloneMatchMapToken_
EDGetTokenT< CandidateView >	zMuStandAloneToken_
EDGetTokenT< GenParticleMatch >	zMuTrackMatchMapToken_
EDGetTokenT< CandidateView >	zMuTrackToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 49 of file ZMuMu_MCanalyzer.cc.

Constructor & Destructor Documentation

ZMuMu_MCanalyzer::ZMuMu_MCanalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 207 of file ZMuMu_MCanalyzer.cc.

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

    : zMuMuToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuMu"))),
      zMuMuMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuMuMatchMap"))),
      zMuStandAloneToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuStandAlone"))),
      zMuStandAloneMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuStandAloneMatchMap"))),
      zMuTrackToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuTrack"))),
      zMuTrackMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuTrackMatchMap"))),
      muonsToken_(consumes<CandidateView>(pset.getParameter<InputTag>("muons"))),
      tracksToken_(consumes<CandidateView>(pset.getParameter<InputTag>("tracks"))),
      genParticlesToken_(consumes<GenParticleCollection>(pset.getParameter<InputTag>("genParticles"))),

      bothMuons_(pset.getParameter<bool>("bothMuons")),
      etamin_(pset.getUntrackedParameter<double>("etamin")),
      etamax_(pset.getUntrackedParameter<double>("etamax")),
      ptmin_(pset.getUntrackedParameter<double>("ptmin")),
      massMin_(pset.getUntrackedParameter<double>("zMassMin")),
      massMax_(pset.getUntrackedParameter<double>("zMassMax")),
      isoMax_(pset.getUntrackedParameter<double>("isomax")),
      ptThreshold_(pset.getUntrackedParameter<double>("ptThreshold")),
      etEcalThreshold_(pset.getUntrackedParameter<double>("etEcalThreshold")),
      etHcalThreshold_(pset.getUntrackedParameter<double>("etHcalThreshold")),
      dRVetoTrk_(pset.getUntrackedParameter<double>("deltaRVetoTrk")),
      dRTrk_(pset.getUntrackedParameter<double>("deltaRTrk")),
      dREcal_(pset.getUntrackedParameter<double>("deltaREcal")),
      dRHcal_(pset.getUntrackedParameter<double>("deltaRHcal")),
      alpha_(pset.getUntrackedParameter<double>("alpha")),
      beta_(pset.getUntrackedParameter<double>("beta")),
      relativeIsolation_(pset.getUntrackedParameter<bool>("relativeIsolation")),
      hltPath_(pset.getUntrackedParameter<std::string>("hltPath")) {
  Service<TFileService> fs;

  // binning of entries array (at moment defined by hand and not in cfg file)
  double etaRange[8] = {-2.5, -2., -1.2, -0.8, 0.8, 1.2, 2., 2.5};
  double ptRange[4] = {20., 40., 60., 100.};

  // general histograms
  h_trackProbe_eta = fs->make<TH1D>("trackProbeEta", "Eta of tracks", 7, etaRange);
  h_trackProbe_pt = fs->make<TH1D>("trackProbePt", "Pt of tracks", 3, ptRange);
  h_staProbe_eta = fs->make<TH1D>("standAloneProbeEta", "Eta of standAlone", 7, etaRange);
  h_staProbe_pt = fs->make<TH1D>("standAloneProbePt", "Pt of standAlone", 3, ptRange);
  h_ProbeOk_eta = fs->make<TH1D>("probeOkEta", "Eta of probe Ok", 7, etaRange);
  h_ProbeOk_pt = fs->make<TH1D>("probeOkPt", "Pt of probe ok", 3, ptRange);

  // clear global counters
  nGlobalMuonsMatched_passed = 0;
  nGlobalMuonsMatched_passedIso = 0;
  n2GlobalMuonsMatched_passedIso = 0;
  nStaMuonsMatched_passedIso = 0;
  nTracksMuonsMatched_passedIso = 0;
  n2GlobalMuonsMatched_passedIso2Trg = 0;
  nMu0onlyTriggered = 0;
  nMu1onlyTriggered = 0;
  nZMuMu_matched = 0;
  nZMuSta_matched = 0;
  nZMuTrk_matched = 0;
}

Member Function Documentation

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

overrideprivate

Definition at line 264 of file ZMuMu_MCanalyzer.cc.

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

                                                                          {
  Handle<CandidateView> zMuMu;
  Handle<GenParticleMatch> zMuMuMatchMap;  //Map of Z made by Mu global + Mu global
  Handle<CandidateView> zMuStandAlone;
  Handle<GenParticleMatch> zMuStandAloneMatchMap;  //Map of Z made by Mu + StandAlone
  Handle<CandidateView> zMuTrack;
  Handle<GenParticleMatch> zMuTrackMatchMap;  //Map of Z made by Mu + Track
  Handle<CandidateView> muons;                //Collection of Muons
  Handle<CandidateView> tracks;               //Collection of Tracks

  Handle<GenParticleCollection> genParticles;  // Collection of Generatd Particles

  event.getByToken(zMuMuToken_, zMuMu);
  event.getByToken(zMuStandAloneToken_, zMuStandAlone);
  event.getByToken(zMuTrackToken_, zMuTrack);
  event.getByToken(genParticlesToken_, genParticles);
  event.getByToken(muonsToken_, muons);
  event.getByToken(tracksToken_, tracks);

  /*
  cout << "*********  zMuMu         size : " << zMuMu->size() << endl;
  cout << "*********  zMuStandAlone size : " << zMuStandAlone->size() << endl;
  cout << "*********  zMuTrack      size : " << zMuTrack->size() << endl;
  cout << "*********  muons         size : " << muons->size() << endl;
  cout << "*********  tracks        size : " << tracks->size() << endl;
  */
  //      std::cout<<"Run-> "<<event.id().run()<<std::endl;
  //      std::cout<<"Event-> "<<event.id().event()<<std::endl;

  bool zMuMu_found = false;

  // loop on ZMuMu
  if (!zMuMu->empty()) {
    event.getByToken(zMuMuMatchMapToken_, zMuMuMatchMap);
    for (unsigned int i = 0; i < zMuMu->size(); ++i) {  //loop on candidates
      const Candidate& zMuMuCand = (*zMuMu)[i];         //the candidate
      CandidateBaseRef zMuMuCandRef = zMuMu->refAt(i);

      const Candidate* lep0 = zMuMuCand.daughter(0);
      const Candidate* lep1 = zMuMuCand.daughter(1);
      const pat::Muon& muonDau0 = dynamic_cast<const pat::Muon&>(*lep0->masterClone());
      //      double trkiso0 = muonDau0.trackIso();
      const pat::Muon& muonDau1 = dynamic_cast<const pat::Muon&>(*lep1->masterClone());
      //double trkiso1 = muonDau1.trackIso();

      double iso0 = candidateIsolation(lep0,
                                       ptThreshold_,
                                       etEcalThreshold_,
                                       etHcalThreshold_,
                                       dRVetoTrk_,
                                       dRTrk_,
                                       dREcal_,
                                       dRHcal_,
                                       alpha_,
                                       beta_,
                                       relativeIsolation_);
      double iso1 = candidateIsolation(lep1,
                                       ptThreshold_,
                                       etEcalThreshold_,
                                       etHcalThreshold_,
                                       dRVetoTrk_,
                                       dRTrk_,
                                       dREcal_,
                                       dRHcal_,
                                       alpha_,
                                       beta_,
                                       relativeIsolation_);

      double pt0 = zMuMuCand.daughter(0)->pt();
      double pt1 = zMuMuCand.daughter(1)->pt();
      double eta0 = zMuMuCand.daughter(0)->eta();
      double eta1 = zMuMuCand.daughter(1)->eta();
      double mass = zMuMuCand.mass();

      // HLT match
      const pat::TriggerObjectStandAloneCollection mu0HLTMatches = muonDau0.triggerObjectMatchesByPath(hltPath_);
      const pat::TriggerObjectStandAloneCollection mu1HLTMatches = muonDau1.triggerObjectMatchesByPath(hltPath_);

      bool trig0found = false;
      bool trig1found = false;
      if (!mu0HLTMatches.empty())
        trig0found = true;
      if (!mu1HLTMatches.empty())
        trig1found = true;

      GenParticleRef zMuMuMatch = (*zMuMuMatchMap)[zMuMuCandRef];
      if (zMuMuMatch.isNonnull()) {  // ZMuMu matched
        zMuMu_found = true;
        nZMuMu_matched++;
        if (pt0 > ptmin_ && pt1 > ptmin_ && abs(eta0) > etamin_ && abs(eta1) > etamin_ && abs(eta0) < etamax_ &&
            abs(eta1) < etamax_ && mass > massMin_ && mass < massMax_ &&
            (trig0found || trig1found)) {  // kinematic and trigger cuts passed
          nGlobalMuonsMatched_passed++;    // first global Muon passed kine cuts
          nGlobalMuonsMatched_passed++;    // second global muon passsed kine cuts
          if (iso0 < isoMax_)
            nGlobalMuonsMatched_passedIso++;  // first global muon passed the iso cut
          if (iso1 < isoMax_)
            nGlobalMuonsMatched_passedIso++;  // second global muon passed the iso cut
          if (iso0 < isoMax_ && iso1 < isoMax_) {
            n2GlobalMuonsMatched_passedIso++;  // both muons passed iso cut
            if (trig0found && trig1found)
              n2GlobalMuonsMatched_passedIso2Trg++;  // both muons have HLT
            if (trig0found && !trig1found)
              nMu0onlyTriggered++;
            if (trig1found && !trig0found)
              nMu1onlyTriggered++;
            // histograms vs eta and pt
            if (trig1found) {  // check efficiency of muon0 not imposing the trigger on it
              h_trackProbe_eta->Fill(eta0);
              h_trackProbe_pt->Fill(pt0);
              h_staProbe_eta->Fill(eta0);
              h_staProbe_pt->Fill(pt0);
              h_ProbeOk_eta->Fill(eta0);
              h_ProbeOk_pt->Fill(pt0);
            }
            if (trig0found) {  // check efficiency of muon1 not imposing the trigger on it
              h_trackProbe_eta->Fill(eta1);
              h_staProbe_eta->Fill(eta1);
              h_trackProbe_pt->Fill(pt1);
              h_staProbe_pt->Fill(pt1);
              h_ProbeOk_eta->Fill(eta1);
              h_ProbeOk_pt->Fill(pt1);
            }
          }
        }
      }  // end MC match

    }  // end loop on ZMuMu cand
  }    // end if ZMuMu size > 0

  // loop on ZMuSta
  bool zMuSta_found = false;
  if (!zMuMu_found && !zMuStandAlone->empty()) {
    event.getByToken(zMuStandAloneMatchMapToken_, zMuStandAloneMatchMap);
    for (unsigned int i = 0; i < zMuStandAlone->size(); ++i) {   //loop on candidates
      const Candidate& zMuStandAloneCand = (*zMuStandAlone)[i];  //the candidate
      CandidateBaseRef zMuStandAloneCandRef = zMuStandAlone->refAt(i);
      GenParticleRef zMuStandAloneMatch = (*zMuStandAloneMatchMap)[zMuStandAloneCandRef];

      const Candidate* lep0 = zMuStandAloneCand.daughter(0);
      const Candidate* lep1 = zMuStandAloneCand.daughter(1);
      const pat::Muon& muonDau0 = dynamic_cast<const pat::Muon&>(*lep0->masterClone());
      //double trkiso0 = muonDau0.trackIso();
      //      const pat::Muon & muonDau1 = dynamic_cast<const pat::Muon &>(*lep1->masterClone());
      //double trkiso1 = muonDau1.trackIso();

      double iso0 = candidateIsolation(lep0,
                                       ptThreshold_,
                                       etEcalThreshold_,
                                       etHcalThreshold_,
                                       dRVetoTrk_,
                                       dRTrk_,
                                       dREcal_,
                                       dRHcal_,
                                       alpha_,
                                       beta_,
                                       relativeIsolation_);
      double iso1 = candidateIsolation(lep1,
                                       ptThreshold_,
                                       etEcalThreshold_,
                                       etHcalThreshold_,
                                       dRVetoTrk_,
                                       dRTrk_,
                                       dREcal_,
                                       dRHcal_,
                                       alpha_,
                                       beta_,
                                       relativeIsolation_);

      double pt0 = zMuStandAloneCand.daughter(0)->pt();
      double pt1 = zMuStandAloneCand.daughter(1)->pt();
      double eta0 = zMuStandAloneCand.daughter(0)->eta();
      double eta1 = zMuStandAloneCand.daughter(1)->eta();
      double mass = zMuStandAloneCand.mass();

      // HLT match (check just dau0 the global)
      const pat::TriggerObjectStandAloneCollection mu0HLTMatches = muonDau0.triggerObjectMatchesByPath(hltPath_);

      bool trig0found = false;
      if (!mu0HLTMatches.empty())
        trig0found = true;

      if (zMuStandAloneMatch.isNonnull()) {  // ZMuStandAlone matched
        zMuSta_found = true;
        nZMuSta_matched++;
        if (pt0 > ptmin_ && pt1 > ptmin_ && abs(eta0) > etamin_ && abs(eta1) > etamin_ && abs(eta0) < etamax_ &&
            abs(eta1) < etamax_ && mass > massMin_ && mass < massMax_ && iso0 < isoMax_ && iso1 < isoMax_ &&
            trig0found) {  // all cuts and trigger passed
          nStaMuonsMatched_passedIso++;
          // histograms vs eta and pt
          h_staProbe_eta->Fill(eta1);
          h_staProbe_pt->Fill(pt1);
        }
      }  // end MC match
    }    // end loop on ZMuStandAlone cand
  }      // end if ZMuStandAlone size > 0

  // loop on ZMuTrack
  if (!zMuMu_found && !zMuSta_found && !zMuTrack->empty()) {
    event.getByToken(zMuTrackMatchMapToken_, zMuTrackMatchMap);
    for (unsigned int i = 0; i < zMuTrack->size(); ++i) {  //loop on candidates
      const Candidate& zMuTrackCand = (*zMuTrack)[i];      //the candidate
      CandidateBaseRef zMuTrackCandRef = zMuTrack->refAt(i);
      const Candidate* lep0 = zMuTrackCand.daughter(0);
      const Candidate* lep1 = zMuTrackCand.daughter(1);
      const pat::Muon& muonDau0 = dynamic_cast<const pat::Muon&>(*lep0->masterClone());
      //double trkiso0 = muonDau0.trackIso();
      //const pat::GenericParticle & trackDau1 = dynamic_cast<const pat::GenericParticle &>(*lep1->masterClone());
      //double trkiso1 = trackDau1.trackIso();
      double iso0 = candidateIsolation(lep0,
                                       ptThreshold_,
                                       etEcalThreshold_,
                                       etHcalThreshold_,
                                       dRVetoTrk_,
                                       dRTrk_,
                                       dREcal_,
                                       dRHcal_,
                                       alpha_,
                                       beta_,
                                       relativeIsolation_);
      double iso1 = candidateIsolation(lep1,
                                       ptThreshold_,
                                       etEcalThreshold_,
                                       etHcalThreshold_,
                                       dRVetoTrk_,
                                       dRTrk_,
                                       dREcal_,
                                       dRHcal_,
                                       alpha_,
                                       beta_,
                                       relativeIsolation_);

      double pt0 = zMuTrackCand.daughter(0)->pt();
      double pt1 = zMuTrackCand.daughter(1)->pt();
      double eta0 = zMuTrackCand.daughter(0)->eta();
      double eta1 = zMuTrackCand.daughter(1)->eta();
      double mass = zMuTrackCand.mass();

      // HLT match (check just dau0 the global)
      const pat::TriggerObjectStandAloneCollection mu0HLTMatches = muonDau0.triggerObjectMatchesByPath(hltPath_);

      bool trig0found = false;
      if (!mu0HLTMatches.empty())
        trig0found = true;

      GenParticleRef zMuTrackMatch = (*zMuTrackMatchMap)[zMuTrackCandRef];
      if (zMuTrackMatch.isNonnull()) {  // ZMuTrack matched
        nZMuTrk_matched++;
        if (pt0 > ptmin_ && pt1 > ptmin_ && abs(eta0) > etamin_ && abs(eta1) > etamin_ && abs(eta0) < etamax_ &&
            abs(eta1) < etamax_ && mass > massMin_ && mass < massMax_ && iso0 < isoMax_ && iso1 < isoMax_ &&
            trig0found) {  // all cuts and trigger passed
          nTracksMuonsMatched_passedIso++;
          // histograms vs eta and pt
          h_trackProbe_eta->Fill(eta1);
          h_trackProbe_pt->Fill(pt1);
        }
      }  // end MC match
    }    // end loop on ZMuTrack cand
  }      // end if ZMuTrack size > 0

}  // end analyze

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

private

Definition at line 526 of file ZMuMu_MCanalyzer.cc.

References reco::Candidate::pdgId().

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

void ZMuMu_MCanalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 682 of file ZMuMu_MCanalyzer.cc.

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

                              {
  double eff_Iso = double(nGlobalMuonsMatched_passedIso) / nGlobalMuonsMatched_passed;
  double err_effIso = sqrt(eff_Iso * (1 - eff_Iso) / nGlobalMuonsMatched_passed);

  double n1_afterIso =
      2 * n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered + nTracksMuonsMatched_passedIso;
  double n2_afterIso =
      2 * n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered + nStaMuonsMatched_passedIso;
  double nGLB_afterIso = 2 * n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered;
  double effSta_afterIso =
      (2 * n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered) / n1_afterIso;
  double effTrk_afterIso =
      (2 * n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered) / n2_afterIso;
  double effHLT_afterIso = (2. * n2GlobalMuonsMatched_passedIso2Trg) /
                           (2. * n2GlobalMuonsMatched_passedIso2Trg + nMu0onlyTriggered + nMu1onlyTriggered);
  double err_effHLT_afterIso = sqrt(effHLT_afterIso * (1 - effHLT_afterIso) / nGLB_afterIso);
  double err_effsta_afterIso = sqrt(effSta_afterIso * (1 - effSta_afterIso) / n1_afterIso);
  double err_efftrk_afterIso = sqrt(effTrk_afterIso * (1 - effTrk_afterIso) / n2_afterIso);

  cout << "------------------------------------  Counters  --------------------------------" << endl;

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

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

private

Definition at line 577 of file ZMuMu_MCanalyzer.cc.

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

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

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

private

Definition at line 647 of file ZMuMu_MCanalyzer.cc.

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

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

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

private

Definition at line 612 of file ZMuMu_MCanalyzer.cc.

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

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

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

private

Definition at line 542 of file ZMuMu_MCanalyzer.cc.

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

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

Member Data Documentation

double ZMuMu_MCanalyzer::alpha_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::beta_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

bool ZMuMu_MCanalyzer::bothMuons_

private

Definition at line 75 of file ZMuMu_MCanalyzer.cc.

double ZMuMu_MCanalyzer::dREcal_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::dRHcal_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::dRTrk_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::dRVetoTrk_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etamax_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etamin_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etEcalThreshold_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::etHcalThreshold_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleCollection> ZMuMu_MCanalyzer::genParticlesToken_

private

Definition at line 73 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

reco::CandidateBaseRef ZMuMu_MCanalyzer::globalMuonCandRef_

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

TH1D * ZMuMu_MCanalyzer::h_ProbeOk_eta

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_ProbeOk_pt

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_staProbe_eta

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_staProbe_pt

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D* ZMuMu_MCanalyzer::h_trackProbe_eta

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

TH1D * ZMuMu_MCanalyzer::h_trackProbe_pt

private

Definition at line 87 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze(), and ZMuMu_MCanalyzer().

string ZMuMu_MCanalyzer::hltPath_

private

Definition at line 82 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::isoMax_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::massMax_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::massMin_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::muonsToken_

private

Definition at line 71 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

int ZMuMu_MCanalyzer::n2GlobalMuonsMatched_passedIso

private

Definition at line 92 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::n2GlobalMuonsMatched_passedIso2Trg

private

Definition at line 95 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nGlobalMuonsMatched_passed

private

Definition at line 90 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nGlobalMuonsMatched_passedIso

private

Definition at line 91 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nMu0onlyTriggered

private

Definition at line 96 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nMu1onlyTriggered

private

Definition at line 97 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nStaMuonsMatched_passedIso

private

Definition at line 93 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nTracksMuonsMatched_passedIso

private

Definition at line 94 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nZMuMu_matched

private

Definition at line 99 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nZMuSta_matched

private

Definition at line 100 of file ZMuMu_MCanalyzer.cc.

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

int ZMuMu_MCanalyzer::nZMuTrk_matched

private

Definition at line 101 of file ZMuMu_MCanalyzer.cc.

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

OverlapChecker ZMuMu_MCanalyzer::overlap_

private

Definition at line 84 of file ZMuMu_MCanalyzer.cc.

double ZMuMu_MCanalyzer::ptmin_

private

Definition at line 77 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

double ZMuMu_MCanalyzer::ptThreshold_

private

Definition at line 79 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

bool ZMuMu_MCanalyzer::relativeIsolation_

private

Definition at line 81 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

reco::CandidateBaseRef ZMuMu_MCanalyzer::standAloneMuonCandRef_

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

reco::CandidateBaseRef ZMuMu_MCanalyzer::trackMuonCandRef_

private

Definition at line 83 of file ZMuMu_MCanalyzer.cc.

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::tracksToken_

private

Definition at line 72 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_MCanalyzer::zMuMuMatchMapToken_

private

Definition at line 66 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::zMuMuToken_

private

Definition at line 65 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_MCanalyzer::zMuStandAloneMatchMapToken_

private

Definition at line 68 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::zMuStandAloneToken_

private

Definition at line 67 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<GenParticleMatch> ZMuMu_MCanalyzer::zMuTrackMatchMapToken_

private

Definition at line 70 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> ZMuMu_MCanalyzer::zMuTrackToken_

private

Definition at line 69 of file ZMuMu_MCanalyzer.cc.

Referenced by analyze().