ZMuMu_Radiative_analyzer Class Reference

Inheritance diagram for ZMuMu_Radiative_analyzer:

Public Member Functions
	ZMuMu_Radiative_analyzer (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
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
void	endJob () override

Private Attributes
double	dRTrk_
double	dRVeto_
int	evntcounter
bool	FSR_mu
bool	FSR_mu0
bool	FSR_mu1
bool	FSR_tk
TH1D *	h_Iso_
TH3D *	h_Iso_3D_
TH1D *	h_Iso_FSR_
TH3D *	h_Iso_FSR_3D_
TH1D *	h_ProbeOk_eta_FSR
TH1D *	h_ProbeOk_eta_no_FSR
TH2D *	h_ProbeOk_pt_eta_FSR_
TH2D *	h_ProbeOk_pt_eta_no_FSR_
TH1D *	h_ProbeOk_pt_FSR
TH1D *	h_ProbeOk_pt_no_FSR
TH1D *	h_staProbe_eta_FSR
TH1D *	h_staProbe_eta_no_FSR
TH2D *	h_staProbe_pt_eta_FSR_
TH2D *	h_staProbe_pt_eta_no_FSR_
TH1D *	h_staProbe_pt_FSR
TH1D *	h_staProbe_pt_no_FSR
TH1D *	h_trackProbe_eta_FSR
TH1D *	h_trackProbe_eta_no_FSR
TH1D *	h_trackProbe_pt_FSR
TH1D *	h_trackProbe_pt_no_FSR
TH1D *	h_zmass_FSR
TH1D *	h_zmass_no_FSR
TH1D *	h_zMuSamass_FSR
TH1D *	h_zMuSamass_no_FSR
TH1D *	h_zMuTkmass_FSR
TH1D *	h_zMuTkmass_no_FSR
double	ptThreshold_
bool	trig0found
bool	trig1found
int	zmmcounter
int	zmscounter
int	zmtcounter
EDGetTokenT< GenParticleMatch >	zMuMuMatchMapToken_
EDGetTokenT< CandidateView >	zMuMuToken_
EDGetTokenT< GenParticleMatch >	zMuSaMatchMapToken_
EDGetTokenT< CandidateView >	zMuSaToken_
EDGetTokenT< GenParticleMatch >	zMuTkMatchMapToken_
EDGetTokenT< CandidateView >	zMuTkToken_

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
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>
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)
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 51 of file ZMuMu_radiative_analysis.cc.

Constructor & Destructor Documentation

ZMuMu_Radiative_analyzer()

ZMuMu_Radiative_analyzer::ZMuMu_Radiative_analyzer

(

const edm::ParameterSet &

pset

)

Definition at line 86 of file ZMuMu_radiative_analysis.cc.

    : zMuMuToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuMu"))),
      zMuMuMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuMuMatchMap"))),
      zMuTkToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuTk"))),
      zMuTkMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuTkMatchMap"))),
      zMuSaToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuSa"))),
      zMuSaMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuSaMatchMap"))),
      dRVeto_(pset.getUntrackedParameter<double>("veto")),
      dRTrk_(pset.getUntrackedParameter<double>("deltaRTrk")),
      ptThreshold_(pset.getUntrackedParameter<double>("ptThreshold")) {
  zmmcounter = 0;
  zmscounter = 0;
  zmtcounter = 0;
  evntcounter = 0;
  Service<TFileService> fs;
 
  // general histograms
  h_zmass_FSR = fs->make<TH1D>("h_zmass_FRS", "Invariant Z mass distribution", 200, 0, 200);
  h_zmass_no_FSR = fs->make<TH1D>("h_zmass_no_FSR", "Invariant Z mass distribution", 200, 0, 200);
  h_zMuSamass_FSR = fs->make<TH1D>("h_zMuSamass_FRS", "Invariant Z mass distribution", 200, 0, 200);
  h_zMuSamass_no_FSR = fs->make<TH1D>("h_zMuSamass_no_FSR", "Invariant Z mass distribution", 200, 0, 200);
  h_zMuTkmass_FSR = fs->make<TH1D>("h_zMuTkmass_FRS", "Invariant Z mass distribution", 200, 0, 200);
  h_zMuTkmass_no_FSR = fs->make<TH1D>("h_zMuTkmass_no_FSR", "Invariant Z mass distribution", 200, 0, 200);
  h_Iso_ = fs->make<TH1D>("h_iso", "Isolation distribution of muons without FSR", 100, 0, 20);
  h_Iso_FSR_ = fs->make<TH1D>("h_iso_FSR", "Isolation distribution of muons with FSR ", 100, 0, 20);
  h_Iso_3D_ = fs->make<TH3D>(
      "h_iso_3D", "Isolation distribution of muons without FSR", 100, 20, 100, 100, -2.0, 2.0, 100, 0, 20);
  h_Iso_FSR_3D_ = fs->make<TH3D>(
      "h_iso_FSR_3D", "Isolation distribution of muons with FSR ", 100, 20, 100, 100, -2.0, 2.0, 100, 0, 20);
  h_staProbe_pt_eta_no_FSR_ =
      fs->make<TH2D>("h_staProbe_pt_eta_no_FSR", "Pt vs Eta StandAlone without FSR ", 100, 20, 100, 100, -2.0, 2.0);
  h_staProbe_pt_eta_FSR_ =
      fs->make<TH2D>("h_staProbe_pt_eta_FSR", "Pt vs Eta StandAlone with FSR ", 100, 20, 100, 100, -2.0, 2.0);
  h_ProbeOk_pt_eta_no_FSR_ =
      fs->make<TH2D>("h_ProbeOk_pt_eta_no_FSR", "Pt vs Eta probeOk without FSR ", 100, 20, 100, 100, -2.0, 2.0);
  h_ProbeOk_pt_eta_FSR_ =
      fs->make<TH2D>("h_ProbeOk_pt_eta_FSR", "Pt vs Eta probeOk with FSR ", 100, 20, 100, 100, -2.0, 2.0);
 
  h_trackProbe_eta_no_FSR = fs->make<TH1D>("trackProbeEta_no_FSR", "Eta of tracks", 100, -2.0, 2.0);
  h_trackProbe_pt_no_FSR = fs->make<TH1D>("trackProbePt_no_FSR", "Pt of tracks", 100, 0, 100);
  h_staProbe_eta_no_FSR = fs->make<TH1D>("standAloneProbeEta_no_FSR", "Eta of standAlone", 100, -2.0, 2.0);
  h_staProbe_pt_no_FSR = fs->make<TH1D>("standAloneProbePt_no_FSR", "Pt of standAlone", 100, 0, 100);
  h_ProbeOk_eta_no_FSR = fs->make<TH1D>("probeOkEta_no_FSR", "Eta of probe Ok", 100, -2.0, 2.0);
  h_ProbeOk_pt_no_FSR = fs->make<TH1D>("probeOkPt_no_FSR", "Pt of probe ok", 100, 0, 100);
 
  h_trackProbe_eta_FSR = fs->make<TH1D>("trackProbeEta_FSR", "Eta of tracks", 100, -2.0, 2.0);
  h_trackProbe_pt_FSR = fs->make<TH1D>("trackProbePt_FSR", "Pt of tracks", 100, 0, 100);
  h_staProbe_eta_FSR = fs->make<TH1D>("standAloneProbeEta_FSR", "Eta of standAlone", 100, -2.0, 2.0);
  h_staProbe_pt_FSR = fs->make<TH1D>("standAloneProbePt_FSR", "Pt of standAlone", 100, 0, 100);
  h_ProbeOk_eta_FSR = fs->make<TH1D>("probeOkEta_FSR", "Eta of probe Ok", 100, -2.0, 2.0);
  h_ProbeOk_pt_FSR = fs->make<TH1D>("probeOkPt_FSR", "Pt of probe ok", 100, 0, 100);
}

References evntcounter, h_Iso_, h_Iso_3D_, h_Iso_FSR_, h_Iso_FSR_3D_, h_ProbeOk_eta_FSR, h_ProbeOk_eta_no_FSR, h_ProbeOk_pt_eta_FSR_, h_ProbeOk_pt_eta_no_FSR_, h_ProbeOk_pt_FSR, h_ProbeOk_pt_no_FSR, h_staProbe_eta_FSR, h_staProbe_eta_no_FSR, h_staProbe_pt_eta_FSR_, h_staProbe_pt_eta_no_FSR_, h_staProbe_pt_FSR, h_staProbe_pt_no_FSR, h_trackProbe_eta_FSR, h_trackProbe_eta_no_FSR, h_trackProbe_pt_FSR, h_trackProbe_pt_no_FSR, h_zmass_FSR, h_zmass_no_FSR, h_zMuSamass_FSR, h_zMuSamass_no_FSR, h_zMuTkmass_FSR, h_zMuTkmass_no_FSR, TFileService::make(), zmmcounter, zmscounter, and zmtcounter.

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 139 of file ZMuMu_radiative_analysis.cc.

                                                                                  {
  evntcounter++;
  Handle<CandidateView> zMuMu;             //Collection of Z made by  Mu global + Mu global
  Handle<GenParticleMatch> zMuMuMatchMap;  //Map of Z made by Mu global + Mu global with MC
  event.getByToken(zMuMuToken_, zMuMu);
  Handle<CandidateView> zMuTk;  //Collection of Z made by  Mu global + Track
  Handle<GenParticleMatch> zMuTkMatchMap;
  event.getByToken(zMuTkToken_, zMuTk);
  Handle<CandidateView> zMuSa;  //Collection of Z made by  Mu global + Sa
  Handle<GenParticleMatch> zMuSaMatchMap;
  event.getByToken(zMuSaToken_, zMuSa);
  cout << "**********  New Event  ***********" << endl;
  // 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);
 
      CandidateBaseRef dau0 = zMuMuCand.daughter(0)->masterClone();
      CandidateBaseRef dau1 = zMuMuCand.daughter(1)->masterClone();
      const pat::Muon& mu0 = dynamic_cast<const pat::Muon&>(*dau0);  //cast in patMuon
      const pat::Muon& mu1 = dynamic_cast<const pat::Muon&>(*dau1);
 
      double zmass = zMuMuCand.mass();
      double pt0 = mu0.pt();
      double pt1 = mu1.pt();
      double eta0 = mu0.eta();
      double eta1 = mu1.eta();
      if (pt0 > 20 && pt1 > 20 && abs(eta0) < 2 && abs(eta1) < 2 && zmass > 20 && zmass < 200) {
        GenParticleRef zMuMuMatch = (*zMuMuMatchMap)[zMuMuCandRef];
        if (zMuMuMatch.isNonnull()) {  // ZMuMu matched
          zmmcounter++;
          cout << "         Zmumu cuts && matched" << endl;
          FSR_mu0 = false;
          FSR_mu1 = false;
 
          //Isodeposit
          const pat::IsoDeposit* mu0TrackIso = mu0.isoDeposit(pat::TrackIso);
          const pat::IsoDeposit* mu1TrackIso = mu1.isoDeposit(pat::TrackIso);
          Direction mu0Dir = Direction(mu0.eta(), mu0.phi());
          Direction mu1Dir = Direction(mu1.eta(), mu1.phi());
 
          reco::IsoDeposit::AbsVetos vetos_mu0;
          vetos_mu0.push_back(new ConeVeto(mu0Dir, dRVeto_));
          vetos_mu0.push_back(new ThresholdVeto(ptThreshold_));
 
          reco::IsoDeposit::AbsVetos vetos_mu1;
          vetos_mu1.push_back(new ConeVeto(mu1Dir, dRVeto_));
          vetos_mu1.push_back(new ThresholdVeto(ptThreshold_));
 
          double Tracker_isovalue_mu0 = mu0TrackIso->sumWithin(dRTrk_, vetos_mu0);
          double Tracker_isovalue_mu1 = mu1TrackIso->sumWithin(dRTrk_, vetos_mu1);
 
          //trigger study
          trig0found = false;
          trig1found = false;
          const pat::TriggerObjectStandAloneCollection mu0HLTMatches = mu0.triggerObjectMatchesByPath("HLT_Mu9");
          const pat::TriggerObjectStandAloneCollection mu1HLTMatches = mu1.triggerObjectMatchesByPath("HLT_Mu9");
          if (!mu0HLTMatches.empty())
            trig0found = true;
          if (!mu1HLTMatches.empty())
            trig1found = true;
 
          //MonteCarlo Study
          const reco::GenParticle* muMc0 = mu0.genLepton();
          const reco::GenParticle* muMc1 = mu1.genLepton();
          const Candidate* motherMu0 = muMc0->mother();
          int num_dau_muon0 = motherMu0->numberOfDaughters();
          const Candidate* motherMu1 = muMc1->mother();
          int num_dau_muon1 = motherMu1->numberOfDaughters();
          cout << "         muone0" << endl;
          cout << "         num di daughters = " << num_dau_muon0 << endl;
          if (num_dau_muon0 > 1) {
            for (int j = 0; j < num_dau_muon0; ++j) {
              int id = motherMu0->daughter(j)->pdgId();
              cout << "         dauther[" << j << "] pdgId = " << id << endl;
              if (id == 22)
                FSR_mu0 = true;
            }
          }  //end check of gamma
 
          cout << "         muone1" << endl;
          cout << "         num di daughters = " << num_dau_muon1 << endl;
          if (num_dau_muon1 > 1) {
            for (int j = 0; j < num_dau_muon1; ++j) {
              int id = motherMu1->daughter(j)->pdgId();
              cout << "         dauther[" << j << "] pdgId = " << id << endl;
              if (id == 22)
                FSR_mu1 = true;
            }
          }  //end check of gamma
 
          if (FSR_mu0 || FSR_mu1)
            h_zmass_FSR->Fill(zmass);
          else
            h_zmass_no_FSR->Fill(zmass);
 
          if (trig1found) {  // check efficiency of muon0 not imposing the trigger on it
            cout << "muon 1 is triggered " << endl;
            if (FSR_mu0) {
              cout << "and muon 0 does FSR" << endl;
              h_trackProbe_eta_FSR->Fill(eta0);
              h_trackProbe_pt_FSR->Fill(pt0);
              h_staProbe_eta_FSR->Fill(eta0);
              h_staProbe_pt_FSR->Fill(pt0);
              h_staProbe_pt_eta_FSR_->Fill(pt0, eta0);
              h_ProbeOk_eta_FSR->Fill(eta0);
              h_ProbeOk_pt_FSR->Fill(pt0);
              h_ProbeOk_pt_eta_FSR_->Fill(pt0, eta0);
            } else {
              cout << "and muon 0 doesn't FSR" << endl;
              h_trackProbe_eta_no_FSR->Fill(eta0);
              h_trackProbe_pt_no_FSR->Fill(pt0);
              h_staProbe_eta_no_FSR->Fill(eta0);
              h_staProbe_pt_no_FSR->Fill(pt0);
              h_staProbe_pt_eta_no_FSR_->Fill(pt0, eta0);
              h_ProbeOk_eta_no_FSR->Fill(eta0);
              h_ProbeOk_pt_no_FSR->Fill(pt0);
              h_ProbeOk_pt_eta_no_FSR_->Fill(pt0, eta0);
            }
            if (FSR_mu0) {
              h_Iso_FSR_->Fill(Tracker_isovalue_mu0);
              h_Iso_FSR_3D_->Fill(pt0, eta0, Tracker_isovalue_mu0);
            } else {
              h_Iso_->Fill(Tracker_isovalue_mu0);
              h_Iso_3D_->Fill(pt0, eta0, Tracker_isovalue_mu0);
            }
          }
          if (trig0found) {  // check efficiency of muon1 not imposing the trigger on it
            cout << "muon 0 is triggered" << endl;
            if (FSR_mu1) {
              cout << "and muon 1  does FSR" << endl;
              h_trackProbe_eta_FSR->Fill(eta1);
              h_staProbe_eta_FSR->Fill(eta1);
              h_trackProbe_pt_FSR->Fill(pt1);
              h_staProbe_pt_FSR->Fill(pt1);
              h_ProbeOk_eta_FSR->Fill(eta1);
              h_ProbeOk_pt_FSR->Fill(pt1);
              h_staProbe_pt_eta_FSR_->Fill(pt1, eta1);
              h_ProbeOk_pt_eta_FSR_->Fill(pt1, eta1);
 
            } else {
              cout << "and muon 1 doesn't FSR" << endl;
              h_trackProbe_eta_no_FSR->Fill(eta1);
              h_staProbe_eta_no_FSR->Fill(eta1);
              h_trackProbe_pt_no_FSR->Fill(pt1);
              h_staProbe_pt_no_FSR->Fill(pt1);
              h_ProbeOk_eta_no_FSR->Fill(eta1);
              h_ProbeOk_pt_no_FSR->Fill(pt1);
              h_staProbe_pt_eta_no_FSR_->Fill(pt1, eta1);
              h_ProbeOk_pt_eta_no_FSR_->Fill(pt1, eta1);
            }
            if (FSR_mu1) {
              h_Iso_FSR_->Fill(Tracker_isovalue_mu1);
              h_Iso_FSR_3D_->Fill(pt1, eta1, Tracker_isovalue_mu1);
            } else {
              h_Iso_->Fill(Tracker_isovalue_mu1);
              h_Iso_3D_->Fill(pt1, eta1, Tracker_isovalue_mu1);
            }
          }
        }  // end MC match
      }    //end of cuts
    }      // end loop on ZMuMu cand
  }        // end if ZMuMu size > 0
 
  // ZMuSa
  if (!zMuSa->empty()) {
    event.getByToken(zMuSaMatchMapToken_, zMuSaMatchMap);
    for (unsigned int i = 0; i < zMuSa->size(); ++i) {  //loop on candidates
 
      const Candidate& zMuSaCand = (*zMuSa)[i];  //the candidate
      CandidateBaseRef zMuSaCandRef = zMuSa->refAt(i);
      const Candidate* lep0 = zMuSaCand.daughter(0);
      const Candidate* lep1 = zMuSaCand.daughter(1);
      const CandidateBaseRef& dau0 = lep0->masterClone();
      const CandidateBaseRef& dau1 = lep1->masterClone();
      const pat::Muon& mu0 = dynamic_cast<const pat::Muon&>(*dau0);  //cast in patMuon
      const pat::Muon& mu1 = dynamic_cast<const pat::Muon&>(*dau1);
 
      double zmass = zMuSaCand.mass();
      double pt0 = mu0.pt();
      double pt1 = mu1.pt();
      double eta0 = mu0.eta();
      double eta1 = mu1.eta();
      if (pt0 > 20 && pt1 > 20 && abs(eta0) < 2 && abs(eta1) < 2 && zmass > 20 && zmass < 200) {
        GenParticleRef zMuSaMatch = (*zMuSaMatchMap)[zMuSaCandRef];
        if (zMuSaMatch.isNonnull()) {  // ZMuSa matched
          cout << "         Zmusa cuts && matched" << endl;
          FSR_mu0 = false;
          FSR_mu1 = false;
          zmscounter++;
          //Isodeposit
          const pat::IsoDeposit* mu0TrackIso = mu0.isoDeposit(pat::TrackIso);
          const pat::IsoDeposit* mu1TrackIso = mu1.isoDeposit(pat::TrackIso);
          Direction mu0Dir = Direction(mu0.eta(), mu0.phi());
          Direction mu1Dir = Direction(mu1.eta(), mu1.phi());
 
          reco::IsoDeposit::AbsVetos vetos_mu0;
          vetos_mu0.push_back(new ConeVeto(mu0Dir, dRVeto_));
          vetos_mu0.push_back(new ThresholdVeto(ptThreshold_));
 
          reco::IsoDeposit::AbsVetos vetos_mu1;
          vetos_mu1.push_back(new ConeVeto(mu1Dir, dRVeto_));
          vetos_mu1.push_back(new ThresholdVeto(ptThreshold_));
 
          double Tracker_isovalue_mu0 = mu0TrackIso->sumWithin(dRTrk_, vetos_mu0);
          double Tracker_isovalue_mu1 = mu1TrackIso->sumWithin(dRTrk_, vetos_mu1);
 
          // HLT match (check just dau0 the global)
          const pat::TriggerObjectStandAloneCollection mu0HLTMatches = mu0.triggerObjectMatchesByPath("HLT_Mu9");
          const pat::TriggerObjectStandAloneCollection mu1HLTMatches = mu1.triggerObjectMatchesByPath("HLT_Mu9");
          trig0found = false;
          trig1found = false;
          if (!mu0HLTMatches.empty())
            trig0found = true;
          if (!mu1HLTMatches.empty())
            trig1found = true;
 
          //MonteCarlo Study
          const reco::GenParticle* muMc0 = mu0.genLepton();
          const reco::GenParticle* muMc1 = mu1.genLepton();
          const Candidate* motherMu0 = muMc0->mother();
          const Candidate* motherMu1 = muMc1->mother();
          int num_dau_muon0 = motherMu0->numberOfDaughters();
          int num_dau_muon1 = motherMu1->numberOfDaughters();
          cout << "         muone0" << endl;
          cout << "         num di daughters = " << num_dau_muon0 << endl;
          if (num_dau_muon0 > 1) {
            for (int j = 0; j < num_dau_muon0; ++j) {
              int id = motherMu0->daughter(j)->pdgId();
              cout << "         dauther[" << j << "] pdgId = " << id << endl;
              if (id == 22)
                FSR_mu0 = true;
            }
          }  //end check of gamma
 
          cout << "         muone1" << endl;
          cout << "         num di daughters = " << num_dau_muon1 << endl;
          if (num_dau_muon1 > 1) {
            for (int j = 0; j < num_dau_muon1; ++j) {
              int id = motherMu1->daughter(j)->pdgId();
              cout << "         dauther[" << j << "] pdgId = " << id << endl;
              if (id == 22)
                FSR_mu1 = true;
            }
          }  //end check of gamma
          if (FSR_mu0 || FSR_mu1)
            h_zMuSamass_FSR->Fill(zmass);
          else
            h_zMuSamass_no_FSR->Fill(zmass);
          if (lep0->isGlobalMuon() && trig0found) {
            if (FSR_mu1) {
              h_staProbe_eta_FSR->Fill(eta1);
              h_staProbe_pt_FSR->Fill(pt1);
              h_staProbe_pt_eta_FSR_->Fill(pt1, eta1);
 
            } else {
              h_staProbe_eta_no_FSR->Fill(eta1);
              h_staProbe_pt_no_FSR->Fill(pt1);
              h_staProbe_pt_eta_no_FSR_->Fill(pt1, eta1);
            }
            if (FSR_mu1) {
              h_Iso_FSR_->Fill(Tracker_isovalue_mu1);
              h_Iso_FSR_3D_->Fill(pt1, eta1, Tracker_isovalue_mu1);
            } else {
              h_Iso_->Fill(Tracker_isovalue_mu1);
              h_Iso_3D_->Fill(pt1, eta1, Tracker_isovalue_mu1);
            }
          }
          if (lep1->isGlobalMuon() && trig1found) {
            if (FSR_mu0) {
              h_staProbe_eta_FSR->Fill(eta0);
              h_staProbe_pt_FSR->Fill(pt0);
              h_staProbe_pt_eta_FSR_->Fill(pt0, eta0);
 
            } else {
              h_staProbe_eta_no_FSR->Fill(eta0);
              h_staProbe_pt_no_FSR->Fill(pt0);
              h_staProbe_pt_eta_FSR_->Fill(pt0, eta0);
            }
            if (FSR_mu0) {
              h_Iso_FSR_->Fill(Tracker_isovalue_mu0);
              h_Iso_FSR_3D_->Fill(pt0, eta0, Tracker_isovalue_mu0);
            } else {
              h_Iso_->Fill(Tracker_isovalue_mu0);
              h_Iso_3D_->Fill(pt0, eta0, Tracker_isovalue_mu0);
            }
          }
        }  // end MC match
      }    //end of cuts
    }      // end loop on ZMuSa cand
  }        // end if ZMuSa size > 0
 
  //ZMuTk
  if (!zMuTk->empty()) {
    event.getByToken(zMuTkMatchMapToken_, zMuTkMatchMap);
    for (unsigned int i = 0; i < zMuTk->size(); ++i) {  //loop on candidates
      const Candidate& zMuTkCand = (*zMuTk)[i];         //the candidate
      CandidateBaseRef zMuTkCandRef = zMuTk->refAt(i);
      CandidateBaseRef dau0 = zMuTkCand.daughter(0)->masterClone();
      CandidateBaseRef dau1 = zMuTkCand.daughter(1)->masterClone();
      const pat::Muon& mu0 = dynamic_cast<const pat::Muon&>(*dau0);  //cast in patMuon
      const pat::GenericParticle& mu1 = dynamic_cast<const pat::GenericParticle&>(*dau1);
 
      double zmass = zMuTkCand.mass();
      double pt0 = mu0.pt();
      double pt1 = mu1.pt();
      double eta0 = mu0.eta();
      double eta1 = mu1.eta();
      if (pt0 > 20 && pt1 > 20 && abs(eta0) < 2 && abs(eta1) < 2 && zmass > 20 && zmass < 200) {  //kinematical cuts
        GenParticleRef zMuTkMatch = (*zMuTkMatchMap)[zMuTkCandRef];
        if (zMuTkMatch.isNonnull()) {  // ZMuTk matched
          FSR_mu = false;
          FSR_tk = false;
          cout << "          ZmuTk cuts && matched" << endl;
          zmtcounter++;
          //Isodeposit
          const pat::IsoDeposit* muTrackIso = mu0.isoDeposit(pat::TrackIso);
          const pat::IsoDeposit* tkTrackIso = mu1.isoDeposit(pat::TrackIso);
          Direction muDir = Direction(mu0.eta(), mu0.phi());
          Direction tkDir = Direction(mu1.eta(), mu1.phi());
 
          IsoDeposit::AbsVetos vetos_mu;
          vetos_mu.push_back(new ConeVeto(muDir, dRVeto_));
          vetos_mu.push_back(new ThresholdVeto(ptThreshold_));
 
          reco::IsoDeposit::AbsVetos vetos_tk;
          vetos_tk.push_back(new ConeVeto(tkDir, dRVeto_));
          vetos_tk.push_back(new ThresholdVeto(ptThreshold_));
 
          double Tracker_isovalue_mu = muTrackIso->sumWithin(dRTrk_, vetos_mu);
          double Tracker_isovalue_tk = tkTrackIso->sumWithin(dRTrk_, vetos_tk);
 
          //MonteCarlo Study
          const reco::GenParticle* muMc0 = mu0.genLepton();
          const reco::GenParticle* muMc1 = mu1.genParticle();
          const Candidate* motherMu0 = muMc0->mother();
          const Candidate* motherMu1 = muMc1->mother();
          int num_dau_muon0 = motherMu0->numberOfDaughters();
          int num_dau_muon1 = motherMu1->numberOfDaughters();
          cout << "numero di figli muone0 = " << num_dau_muon0 << endl;
          cout << "numero di figli muone1 = " << num_dau_muon1 << endl;
 
          cout << "         muon" << endl;
          cout << "         num di daughters = " << num_dau_muon0 << endl;
          if (num_dau_muon0 > 1) {
            for (int j = 0; j < num_dau_muon0; ++j) {
              int id = motherMu0->daughter(j)->pdgId();
              cout << "         dau[" << j << "] pdg ID = " << id << endl;
              if (id == 22) {
                FSR_mu = true;
              }
            }
          }  //end check of gamma
          else
            cout << "         dau[0] pdg ID = " << motherMu0->daughter(0)->pdgId() << endl;
          cout << "         traccia" << endl;
          cout << "         num di daughters = " << num_dau_muon1 << endl;
          if (num_dau_muon1 > 1) {
            for (int j = 0; j < num_dau_muon1; ++j) {
              int id = motherMu1->daughter(j)->pdgId();
              cout << "         dau[" << j << "] pdg ID = " << id << endl;
              if (id == 22) {
                FSR_tk = true;
              }
            }
          }  //end check of gamma
          else
            cout << "         dau[0] pdg ID = " << motherMu1->daughter(0)->pdgId() << endl;
          cout << "Mu Isolation = " << Tracker_isovalue_mu << endl;
          cout << "Track Isolation = " << Tracker_isovalue_tk << endl;
          if (FSR_mu) {
            h_Iso_FSR_->Fill(Tracker_isovalue_mu);
            h_Iso_FSR_3D_->Fill(pt0, eta0, Tracker_isovalue_mu);
          } else {
            h_Iso_->Fill(Tracker_isovalue_mu);
            h_Iso_3D_->Fill(pt0, eta0, Tracker_isovalue_mu);
          }
          if (FSR_tk) {
            h_Iso_FSR_->Fill(Tracker_isovalue_tk);
            h_Iso_FSR_3D_->Fill(pt1, eta1, Tracker_isovalue_tk);
            h_trackProbe_eta_FSR->Fill(eta1);
            h_trackProbe_pt_FSR->Fill(pt1);
          } else {
            h_Iso_->Fill(Tracker_isovalue_tk);
            h_Iso_3D_->Fill(pt1, eta1, Tracker_isovalue_tk);
            h_trackProbe_eta_no_FSR->Fill(eta1);
            h_trackProbe_pt_no_FSR->Fill(pt1);
          }
        }  // end MC match
      }    //end Kine-cuts
    }      // end loop on ZMuTk cand
  }        // end if ZMuTk size > 0
}  // end analyze

References funct::abs(), gather_cfg::cout, reco::Candidate::daughter(), dRTrk_, dRVeto_, edm::View< T >::empty(), reco::LeafCandidate::eta(), HLT_2018_cff::eta1, evntcounter, FSR_mu, FSR_mu0, FSR_mu1, FSR_tk, pat::Lepton< LeptonType >::genLepton(), pat::PATObject< ObjectType >::genParticle(), h_Iso_, h_Iso_3D_, h_Iso_FSR_, h_Iso_FSR_3D_, h_ProbeOk_eta_FSR, h_ProbeOk_eta_no_FSR, h_ProbeOk_pt_eta_FSR_, h_ProbeOk_pt_eta_no_FSR_, h_ProbeOk_pt_FSR, h_ProbeOk_pt_no_FSR, h_staProbe_eta_FSR, h_staProbe_eta_no_FSR, h_staProbe_pt_eta_FSR_, h_staProbe_pt_eta_no_FSR_, h_staProbe_pt_FSR, h_staProbe_pt_no_FSR, h_trackProbe_eta_FSR, h_trackProbe_eta_no_FSR, h_trackProbe_pt_FSR, h_trackProbe_pt_no_FSR, h_zmass_FSR, h_zmass_no_FSR, h_zMuSamass_FSR, h_zMuSamass_no_FSR, mps_fire::i, reco::Candidate::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), pat::Lepton< LeptonType >::isoDeposit(), pat::GenericParticle::isoDeposit(), dqmiolumiharvest::j, RazorAnalyzer::lep1, reco::Candidate::mass(), reco::Candidate::masterClone(), reco::CompositeRefCandidateT< D >::mother(), reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), HLT_2018_cff::pt1, ptThreshold_, edm::View< T >::refAt(), edm::View< T >::size(), reco::IsoDeposit::sumWithin(), pat::TrackIso, trig0found, trig1found, pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), zmmcounter, zmscounter, zmtcounter, ZMuMuAnalysisNtupler_cff::zMuMu, zMuMuMatchMapToken_, zMuMuToken_, ZMuMuAnalysisNtupler_cff::zMuSa, zMuSaMatchMapToken_, zMuSaToken_, zMuTkMatchMapToken_, and zMuTkToken_.

endJob()

void ZMuMu_Radiative_analyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 536 of file ZMuMu_radiative_analysis.cc.

                                      {
  cout << " ============= Summary ==========" << endl;
  cout << " Numero di eventi " << evntcounter << endl;
  cout << " 1)Numero di ZMuMu matched dopo i tagli cinematici = " << zmmcounter << endl;
  cout << " 2)Numero di ZMuSa matched dopo i tagli cinematici = " << zmscounter << endl;
  cout << " 3)Numero di ZMuTk matched dopo i tagli cinematici = " << zmtcounter << endl;
  double n1 = h_Iso_FSR_->Integral();
  double icut1 = h_Iso_FSR_->Integral(0, 15);
  double eff_iso_FSR = (double)icut1 / (double)n1;
  double err_iso_FSR = sqrt(eff_iso_FSR * (1 - eff_iso_FSR) / n1);
  double n2 = h_Iso_->Integral();
  double icut2 = h_Iso_->Integral(0, 15);
  double eff_iso = (double)icut2 / (double)n2;
  double err_iso = sqrt(eff_iso * (1 - eff_iso) / n2);
  cout << " ============= Isolation Efficiecy ==========" << endl;
  cout << "Isolation Efficiency  = " << eff_iso << " +/- " << err_iso << endl;
  cout << "Isolation Efficiency with FSR = " << eff_iso_FSR << " +/- " << err_iso_FSR << endl;
}

References gather_cfg::cout, evntcounter, h_Iso_, h_Iso_FSR_, mathSSE::sqrt(), zmmcounter, zmscounter, and zmtcounter.

Member Data Documentation

dRTrk_

double ZMuMu_Radiative_analyzer::dRTrk_

private

Definition at line 65 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

dRVeto_

double ZMuMu_Radiative_analyzer::dRVeto_

private

Definition at line 65 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

evntcounter

int ZMuMu_Radiative_analyzer::evntcounter

private

Definition at line 81 of file ZMuMu_radiative_analysis.cc.

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

FSR_mu

bool ZMuMu_Radiative_analyzer::FSR_mu

private

Definition at line 78 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

FSR_mu0

bool ZMuMu_Radiative_analyzer::FSR_mu0

private

Definition at line 78 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

FSR_mu1

bool ZMuMu_Radiative_analyzer::FSR_mu1

private

Definition at line 78 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

FSR_tk

bool ZMuMu_Radiative_analyzer::FSR_tk

private

Definition at line 78 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

h_Iso_

TH1D* ZMuMu_Radiative_analyzer::h_Iso_

private

Definition at line 70 of file ZMuMu_radiative_analysis.cc.

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

h_Iso_3D_

TH3D* ZMuMu_Radiative_analyzer::h_Iso_3D_

private

Definition at line 71 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_Iso_FSR_

TH1D * ZMuMu_Radiative_analyzer::h_Iso_FSR_

private

Definition at line 70 of file ZMuMu_radiative_analysis.cc.

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

h_Iso_FSR_3D_

TH3D * ZMuMu_Radiative_analyzer::h_Iso_FSR_3D_

private

Definition at line 71 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_ProbeOk_eta_FSR

TH1D * ZMuMu_Radiative_analyzer::h_ProbeOk_eta_FSR

private

Definition at line 75 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_ProbeOk_eta_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_ProbeOk_eta_no_FSR

private

Definition at line 73 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_ProbeOk_pt_eta_FSR_

TH2D * ZMuMu_Radiative_analyzer::h_ProbeOk_pt_eta_FSR_

private

Definition at line 72 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_ProbeOk_pt_eta_no_FSR_

TH2D * ZMuMu_Radiative_analyzer::h_ProbeOk_pt_eta_no_FSR_

private

Definition at line 72 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_ProbeOk_pt_FSR

TH1D * ZMuMu_Radiative_analyzer::h_ProbeOk_pt_FSR

private

Definition at line 75 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_ProbeOk_pt_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_ProbeOk_pt_no_FSR

private

Definition at line 73 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_staProbe_eta_FSR

TH1D * ZMuMu_Radiative_analyzer::h_staProbe_eta_FSR

private

Definition at line 75 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_staProbe_eta_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_staProbe_eta_no_FSR

private

Definition at line 73 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_staProbe_pt_eta_FSR_

TH2D * ZMuMu_Radiative_analyzer::h_staProbe_pt_eta_FSR_

private

Definition at line 72 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_staProbe_pt_eta_no_FSR_

TH2D* ZMuMu_Radiative_analyzer::h_staProbe_pt_eta_no_FSR_

private

Definition at line 72 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_staProbe_pt_FSR

TH1D * ZMuMu_Radiative_analyzer::h_staProbe_pt_FSR

private

Definition at line 75 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_staProbe_pt_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_staProbe_pt_no_FSR

private

Definition at line 73 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_trackProbe_eta_FSR

TH1D* ZMuMu_Radiative_analyzer::h_trackProbe_eta_FSR

private

Definition at line 75 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_trackProbe_eta_no_FSR

TH1D* ZMuMu_Radiative_analyzer::h_trackProbe_eta_no_FSR

private

Definition at line 73 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_trackProbe_pt_FSR

TH1D * ZMuMu_Radiative_analyzer::h_trackProbe_pt_FSR

private

Definition at line 75 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_trackProbe_pt_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_trackProbe_pt_no_FSR

private

Definition at line 73 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_zmass_FSR

TH1D* ZMuMu_Radiative_analyzer::h_zmass_FSR

private

Definition at line 67 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_zmass_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_zmass_no_FSR

private

Definition at line 67 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_zMuSamass_FSR

TH1D* ZMuMu_Radiative_analyzer::h_zMuSamass_FSR

private

Definition at line 68 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_zMuSamass_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_zMuSamass_no_FSR

private

Definition at line 68 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze(), and ZMuMu_Radiative_analyzer().

h_zMuTkmass_FSR

TH1D* ZMuMu_Radiative_analyzer::h_zMuTkmass_FSR

private

Definition at line 69 of file ZMuMu_radiative_analysis.cc.

Referenced by ZMuMu_Radiative_analyzer().

h_zMuTkmass_no_FSR

TH1D * ZMuMu_Radiative_analyzer::h_zMuTkmass_no_FSR

private

Definition at line 69 of file ZMuMu_radiative_analysis.cc.

Referenced by ZMuMu_Radiative_analyzer().

ptThreshold_

double ZMuMu_Radiative_analyzer::ptThreshold_

private

Definition at line 65 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

trig0found

bool ZMuMu_Radiative_analyzer::trig0found

private

Definition at line 79 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

trig1found

bool ZMuMu_Radiative_analyzer::trig1found

private

Definition at line 79 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

zmmcounter

int ZMuMu_Radiative_analyzer::zmmcounter

private

Definition at line 81 of file ZMuMu_radiative_analysis.cc.

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

zmscounter

int ZMuMu_Radiative_analyzer::zmscounter

private

Definition at line 81 of file ZMuMu_radiative_analysis.cc.

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

zmtcounter

int ZMuMu_Radiative_analyzer::zmtcounter

private

Definition at line 81 of file ZMuMu_radiative_analysis.cc.

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

zMuMuMatchMapToken_

EDGetTokenT<GenParticleMatch> ZMuMu_Radiative_analyzer::zMuMuMatchMapToken_

private

Definition at line 60 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

zMuMuToken_

EDGetTokenT<CandidateView> ZMuMu_Radiative_analyzer::zMuMuToken_

private

Definition at line 59 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

zMuSaMatchMapToken_

EDGetTokenT<GenParticleMatch> ZMuMu_Radiative_analyzer::zMuSaMatchMapToken_

private

Definition at line 64 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

zMuSaToken_

EDGetTokenT<CandidateView> ZMuMu_Radiative_analyzer::zMuSaToken_

private

Definition at line 63 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

zMuTkMatchMapToken_

EDGetTokenT<GenParticleMatch> ZMuMu_Radiative_analyzer::zMuTkMatchMapToken_

private

Definition at line 62 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().

zMuTkToken_

EDGetTokenT<CandidateView> ZMuMu_Radiative_analyzer::zMuTkToken_

private

Definition at line 61 of file ZMuMu_radiative_analysis.cc.

Referenced by analyze().