ZMuMuAnalyzer_cynematics Class Reference

Inheritance diagram for ZMuMuAnalyzer_cynematics:

Public Member Functions
	ZMuMuAnalyzer_cynematics (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
void	endJob () override
bool	isContained (const Candidate &, const Candidate &)

Private Attributes
double	binEta
double	binPhi
double	binPt
double	etacut_
double	etamax
double	etamin
TH2D *	h_OneSta_mass
TH1D *	h_zMuMu_mass_
TH1D *	h_zMuMu_numberOfCand
TH1D *	h_zMuMu_numberOfCand_etapassed
TH1D *	h_zMuMu_numberOfCand_isopassed
TH1D *	h_zMuMu_numberOfCand_masspassed
TH1D *	h_zMuMu_numberOfCand_passed
TH1D *	h_zMuMu_numberOfCand_ptetamasspassed
TH1D *	h_zMuMu_numberOfCand_ptetapassed
TH1D *	h_zMuMu_numberOfCand_ptpassed
TH1D *	h_zMuMuMatched_mass_
TH1D *	h_zMuSingleStandAlone_mass_
TH1D *	h_zMuSingleStandAloneMatched_mass_
TH1D *	h_zMuSingleStandAloneOverlap_mass_
TH1D *	h_zMuSingleStandAloneOverlapMatched_mass_
TH1D *	h_zMuSingleTrack_mass_
TH1D *	h_zMuSingleTrackMatched_mass_
TH1D *	h_zMuSta_MCmatched_numberOfCand_passed
TH1D *	h_zMuSta_numberOfCand
TH1D *	h_zMuSta_numberOfCand_notcontained
TH1D *	h_zMuSta_numberOfCand_passed
TH1D *	h_zMuTrack_MCmatched_numberOfCand_passed
TH1D *	h_zMuTrack_numberOfCand
TH1D *	h_zMuTrack_numberOfCand_notcontained
TH1D *	h_zMuTrack_numberOfCand_passed
vector< TH1D * >	hmumu_eta
vector< TH1D * >	hmumu_phi
vector< TH1D * >	hmumu_pt
vector< TH1D * >	hmusta_eta
vector< TH1D * >	hmusta_phi
vector< TH1D * >	hmusta_pt
vector< TH1D * >	hmutrack_eta
vector< TH1D * >	hmutrack_phi
vector< TH1D * >	hmutrack_pt
double	isocut_
double	maxZmass_
double	minZmass_
EDGetTokenT< IsolationCollection >	muIsoToken_
int	numberOfIntervals
OverlapChecker	overlap_
double	phimax
double	phimin
double	ptcut_
double	ptmax
double	ptmin
double	ptSTAcut_
EDGetTokenT< IsolationCollection >	standAloneIsoToken_
EDGetTokenT< IsolationCollection >	trackIsoToken_
EDGetTokenT< CandMatchMap >	zMuMuMapToken_
EDGetTokenT< CandidateCollection >	zMuMuToken_
EDGetTokenT< CandMatchMap >	zMuStandAloneMapToken_
EDGetTokenT< CandidateCollection >	zMuStandAloneToken_
EDGetTokenT< CandMatchMap >	zMuTrackMapToken_
EDGetTokenT< CandidateCollection >	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 ZMuMuAnalyzer_cynematics.cc.

Constructor & Destructor Documentation

ZMuMuAnalyzer_cynematics::ZMuMuAnalyzer_cynematics

(

const edm::ParameterSet &

pset

)

Definition at line 88 of file ZMuMuAnalyzer_cynematics.cc.

References a, b, binEta, binPhi, binPt, etacut_, etamax, etamin, h_OneSta_mass, h_zMuMu_mass_, h_zMuMu_numberOfCand, h_zMuMu_numberOfCand_etapassed, h_zMuMu_numberOfCand_isopassed, h_zMuMu_numberOfCand_masspassed, h_zMuMu_numberOfCand_passed, h_zMuMu_numberOfCand_ptetamasspassed, h_zMuMu_numberOfCand_ptetapassed, h_zMuMu_numberOfCand_ptpassed, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, h_zMuSta_MCmatched_numberOfCand_passed, h_zMuSta_numberOfCand, h_zMuSta_numberOfCand_notcontained, h_zMuSta_numberOfCand_passed, h_zMuTrack_MCmatched_numberOfCand_passed, h_zMuTrack_numberOfCand, h_zMuTrack_numberOfCand_notcontained, h_zMuTrack_numberOfCand_passed, hmumu_eta, hmumu_phi, hmumu_pt, hmusta_eta, hmusta_phi, hmusta_pt, hmutrack_eta, hmutrack_phi, hmutrack_pt, mps_fire::i, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), numberOfIntervals, phimax, phimin, ptcut_, ptmax, and ptmin.

    : zMuMuToken_(consumes<CandidateCollection>(pset.getParameter<InputTag>("zMuMu"))),
      zMuTrackToken_(consumes<CandidateCollection>(pset.getParameter<InputTag>("zMuTrack"))),
      zMuStandAloneToken_(consumes<CandidateCollection>(pset.getParameter<InputTag>("zMuStandAlone"))),
      muIsoToken_(consumes<IsolationCollection>(pset.getParameter<InputTag>("muIso"))),
      trackIsoToken_(consumes<IsolationCollection>(pset.getParameter<InputTag>("trackIso"))),
      standAloneIsoToken_(consumes<IsolationCollection>(pset.getParameter<InputTag>("standAloneIso"))),
      zMuMuMapToken_(mayConsume<CandMatchMap>(pset.getParameter<InputTag>("zMuMuMap"))),
      zMuTrackMapToken_(mayConsume<CandMatchMap>(pset.getParameter<InputTag>("zMuTrackMap"))),
      zMuStandAloneMapToken_(mayConsume<CandMatchMap>(pset.getParameter<InputTag>("zMuStandAloneMap"))),
      isocut_(pset.getParameter<double>("isocut")),
      etacut_(pset.getParameter<double>("etacut")),
      ptcut_(pset.getParameter<double>("ptcut")),
      ptSTAcut_(pset.getParameter<double>("ptSTAcut")),
      minZmass_(pset.getParameter<double>("minZmass")),
      maxZmass_(pset.getParameter<double>("maxZmass")) {
  Service<TFileService> fs;
  h_zMuMu_numberOfCand = fs->make<TH1D>("ZMuMunumberOfCand", "number of ZMuMu cand", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_passed =
      fs->make<TH1D>("ZMuMunumberOfCandpassed", "number of ZMuMu cand selected", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_ptpassed =
      fs->make<TH1D>("ZMuMunumberOfCandptpassed", "number of ZMuMu cand after pt cut selected", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_etapassed =
      fs->make<TH1D>("ZMuMunumberOfCandetapassed", "number of ZMuMu cand after eta cut selected", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_masspassed =
      fs->make<TH1D>("ZMuMunumberOfCandmasspassed", "number of ZMuMu cand after mass cut selected", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_isopassed =
      fs->make<TH1D>("ZMuMunumberOfCandisopassed", "number of ZMuMu cand after iso cut selected", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_ptetapassed =
      fs->make<TH1D>("ZMuMunumberOfCandptetapassed", "number of ZMuMu cand after pt & eta cut selected", 10, -.5, 9.5);
  h_zMuMu_numberOfCand_ptetamasspassed = fs->make<TH1D>(
      "ZMuMunumberOfCandptetamaspassed", "number of ZMuMu cand after pt & eta & mass cut selected", 10, -.5, 9.5);

  h_zMuMu_mass_ = fs->make<TH1D>("ZMuMumass", "ZMuMu mass(GeV)", 200, 0., 200.);
  h_zMuSingleTrack_mass_ = fs->make<TH1D>("ZMuSingleTrackmass", "ZMuSingleTrack mass(GeV)", 100, 0., 200.);
  h_zMuSingleStandAlone_mass_ =
      fs->make<TH1D>("ZMuSingleStandAlonemass", "ZMuSingleStandAlone mass(GeV)", 50, 0., 200.);
  h_zMuSingleStandAloneOverlap_mass_ =
      fs->make<TH1D>("ZMuSingleStandAloneOverlapmass", "ZMuSingleStandAloneOverlap  mass(GeV)", 50, 0., 200.);

  h_zMuMuMatched_mass_ = fs->make<TH1D>("ZMuMuMatchedmass", "ZMuMu Matched  mass(GeV)", 200, 0., 200.);
  h_zMuSingleTrackMatched_mass_ =
      fs->make<TH1D>("ZMuSingleTrackmassMatched", "ZMuSingleTrackMatched mass(GeV)", 100, 0., 200.);
  h_zMuSingleStandAloneMatched_mass_ =
      fs->make<TH1D>("ZMuSingleStandAlonemassMatched", "ZMuSingleStandAloneMatched mass(GeV)", 50, 0., 200.);
  h_zMuSingleStandAloneOverlapMatched_mass_ = fs->make<TH1D>(
      "ZMuSingleStandAloneOverlapmassMatched", "ZMuSingleStandAloneMatched Overlap  mass(GeV)", 50, 0., 200.);

  h_zMuSta_numberOfCand =
      fs->make<TH1D>("ZMuStanumberOfCand", "number of ZMuSta cand (if ZMuMu not selected)", 10, -.5, 9.5);
  h_OneSta_mass =
      fs->make<TH2D>("ZOneMuStaMass",
                     "inv. mass of ZMuSta1 vs ZMuSta2 when one ZMuSta has been found (if ZMuMu not selected)",
                     100,
                     0.,
                     400,
                     100,
                     0.,
                     400.);
  h_zMuSta_numberOfCand_notcontained = fs->make<TH1D>(
      "ZMuStanumberOfCandnotcontained", "number of independent ZMuSta cand (if ZMuMu not selected)", 10, -.5, 9.5);
  h_zMuSta_numberOfCand_passed = fs->make<TH1D>(
      "ZMuStanumberOfCandpassed", "number of ZMuSta cand selected (if ZMuMu not selected)", 10, -.5, 9.5);
  h_zMuSta_MCmatched_numberOfCand_passed =
      fs->make<TH1D>("ZMuStaMCmatchedNumberOfCandpassed",
                     "number of ZMuSta MC matched cand selected (if ZMuMu not selected)",
                     10,
                     -.5,
                     9.5);
  h_zMuTrack_numberOfCand =
      fs->make<TH1D>("ZMuTranumberOfCand", "number of ZMuTrack cand (if ZMuMu and ZMuSTa not selected)", 10, -.5, 9.5);
  h_zMuTrack_numberOfCand_notcontained =
      fs->make<TH1D>("ZMuTranumberOfCandnotcontaind",
                     "number of indeendent ZMuTrack cand (if ZMuMu and ZMuSTa not selected)",
                     10,
                     -.5,
                     9.5);
  h_zMuTrack_numberOfCand_passed = fs->make<TH1D>(
      "ZMuTranumberOfCandpassed", "number of ZMuTrack cand selected (if ZMuMu and ZMuSta not selected)", 10, -.5, 9.5);
  h_zMuTrack_MCmatched_numberOfCand_passed =
      fs->make<TH1D>("ZMuTraMCmacthedNumberOfCandpassed",
                     "number of ZMuTrack MC matched cand selected (if ZMuMu and ZMuSta not selected)",
                     10,
                     -.5,
                     9.5);

  // creating histograms for each Pt, eta, phi interval

  etamin = -etacut_;
  etamax = etacut_;
  phimin = -3.1415;
  phimax = 3.1415;
  ptmin = ptcut_;
  ptmax = 100;
  numberOfIntervals = 8;  // number of intervals in which to divide cynematic variables
  binEta = (etamax - etamin) / numberOfIntervals;
  binPhi = (phimax - phimin) / numberOfIntervals;
  binPt = (ptmax - ptmin) / numberOfIntervals;
  TFileDirectory etaDirectory =
      fs->mkdir("etaIntervals");  // in this directory will be saved all the histos of different eta intervals
  TFileDirectory phiDirectory =
      fs->mkdir("phiIntervals");  // in this directory will be saved all the histos of different phi intervals
  TFileDirectory ptDirectory =
      fs->mkdir("ptIntervals");  // in this directory will be saved all the histos of different pt intervals

  // eta histograms creation

  for (int i = 0; i < numberOfIntervals; i++) {
    double range0 = etamin + i * binEta;
    double range1 = range0 + binEta;
    char a[30], b[50];
    sprintf(a, "zmumu_etaRange%d", i);
    sprintf(b, "zmumu mass eta Range %f to %f", range0, range1);
    hmumu_eta.push_back(etaDirectory.make<TH1D>(a, b, 200, 0., 200.));
    char asta[30], bsta[50];
    sprintf(asta, "zmusta_etaRange%d", i);
    sprintf(bsta, "zmusta mass eta Range %f to %f", range0, range1);
    hmusta_eta.push_back(etaDirectory.make<TH1D>(asta, bsta, 50, 0., 200.));
    char atk[30], btk[50];
    sprintf(atk, "zmutrack_etaRange%d", i);
    sprintf(btk, "zmutrack mass eta Range %f to %f", range0, range1);
    hmutrack_eta.push_back(etaDirectory.make<TH1D>(atk, btk, 100, 0., 200.));
  }

  // phi histograms creation

  for (int i = 0; i < numberOfIntervals; i++) {
    double range0 = phimin + i * binPhi;
    double range1 = range0 + binPhi;
    char a[30], b[50];
    sprintf(a, "zmumu_phiRange%d", i);
    sprintf(b, "zmumu mass phi Range %f to %f", range0, range1);
    hmumu_phi.push_back(phiDirectory.make<TH1D>(a, b, 200, 0., 200.));
    char asta[30], bsta[50];
    sprintf(asta, "zmusta_phiRange%d", i);
    sprintf(bsta, "zmusta mass phi Range %f to %f", range0, range1);
    hmusta_phi.push_back(phiDirectory.make<TH1D>(asta, bsta, 50, 0., 200.));
    char atk[30], btk[50];
    sprintf(atk, "zmutrack_phiRange%d", i);
    sprintf(btk, "zmutrack mass phi Range %f to %f", range0, range1);
    hmutrack_phi.push_back(phiDirectory.make<TH1D>(atk, btk, 100, 0., 200.));
  }

  // pt histograms creation

  for (int i = 0; i < numberOfIntervals; i++) {
    double range0 = ptmin + i * binPt;
    double range1 = range0 + binPt;
    char a[30], b[50];
    sprintf(a, "zmumu_ptRange%d", i);
    sprintf(b, "zmumu mass pt Range %f to %f", range0, range1);
    hmumu_pt.push_back(ptDirectory.make<TH1D>(a, b, 200, 0., 200.));
    char asta[30], bsta[50];
    sprintf(asta, "zmusta_ptRange%d", i);
    sprintf(bsta, "zmusta mass pt Range %f to %f", range0, range1);
    hmusta_pt.push_back(ptDirectory.make<TH1D>(asta, bsta, 50, 0., 200.));
    char atk[30], btk[50];
    sprintf(atk, "zmutrack_ptRange%d", i);
    sprintf(btk, "zmutrack mass pt Range %f to %f", range0, range1);
    hmutrack_pt.push_back(ptDirectory.make<TH1D>(atk, btk, 100, 0., 200.));
  }
}

Member Function Documentation

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

overrideprivate

Definition at line 251 of file ZMuMuAnalyzer_cynematics.cc.

References binEta, binPhi, binPt, edm::RefToBase< T >::castTo(), gather_cfg::cout, reco::Candidate::daughter(), edm::AssociationMap< Tag >::end(), etacut_, etamin, HcalObjRepresent::Fill(), edm::AssociationMap< Tag >::find(), edm::Ref< C, T, F >::get(), h_zMuMu_mass_, h_zMuMu_numberOfCand, h_zMuMu_numberOfCand_etapassed, h_zMuMu_numberOfCand_isopassed, h_zMuMu_numberOfCand_masspassed, h_zMuMu_numberOfCand_passed, h_zMuMu_numberOfCand_ptetamasspassed, h_zMuMu_numberOfCand_ptetapassed, h_zMuMu_numberOfCand_ptpassed, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, h_zMuSta_MCmatched_numberOfCand_passed, h_zMuSta_numberOfCand, h_zMuSta_numberOfCand_notcontained, h_zMuSta_numberOfCand_passed, h_zMuTrack_MCmatched_numberOfCand_passed, h_zMuTrack_numberOfCand, h_zMuTrack_numberOfCand_notcontained, h_zMuTrack_numberOfCand_passed, hmumu_eta, hmumu_phi, hmumu_pt, hmusta_eta, hmusta_phi, hmusta_pt, hmutrack_eta, hmutrack_phi, hmutrack_pt, mps_fire::i, edm::Ref< C, T, F >::id(), globals_cff::id1, globals_cff::id2, isContained(), isocut_, dqmiolumiharvest::j, edm::Ref< C, T, F >::key(), RazorAnalyzer::lep1, RazorAnalyzer::lep2, visualization-live-secondInstance_cfg::m, reco::Candidate::mass(), reco::Candidate::masterClone(), maxZmass_, minZmass_, AnalysisDataFormats_SUSYBSMObjects::ms, TtSemiLepEvtBuilder_cfi::mt, muIsoToken_, numberOfIntervals, phimin, ptcut_, ptmin, ptSTAcut_, edm::OwnVector< T, P >::size(), standAloneIsoToken_, trackIsoToken_, ZMuMuAnalysisNtupler_cff::zMuMu, zMuMuMapToken_, zMuMuToken_, zMuStandAloneMapToken_, zMuStandAloneToken_, zMuTrackMapToken_, and zMuTrackToken_.

                                                                                        {
  Handle<CandidateCollection> zMuMu;
  event.getByToken(zMuMuToken_, zMuMu);
  Handle<CandidateCollection> zMuTrack;
  event.getByToken(zMuTrackToken_, zMuTrack);
  Handle<CandidateCollection> zMuStandAlone;
  event.getByToken(zMuStandAloneToken_, zMuStandAlone);

  unsigned int nZMuMu = zMuMu->size();
  unsigned int nZTrackMu = zMuTrack->size();
  unsigned int nZStandAloneMu = zMuStandAlone->size();
  //  static const double zMass = 91.1876; // PDG Z mass

  cout << "++++++++++++++++++++++++++" << endl;
  cout << "nZMuMu = " << nZMuMu << endl;
  cout << "nZTrackMu = " << nZTrackMu << endl;
  cout << "nZStandAloneMu = " << nZStandAloneMu << endl;
  cout << "++++++++++++++++++++++++++" << endl;

  // ZMuMu counters

  int ZMuMu_passed = 0;
  int ZMuMu_ptcut_counter = 0;
  int ZMuMu_etacut_counter = 0;
  int ZMuMu_masscut_counter = 0;
  int ZMuMu_isocut_counter = 0;
  int ZMuMu_ptetacut_counter = 0;
  int ZMuMu_ptetamasscut_counter = 0;
  int ZMuMu_allcut_counter = 0;

  // ZMuTrack counters

  int ZMuTrack_passed = 0;
  int ZMuTrack_notcontained = 0;
  int ZMuTrack_MCmatched_passed = 0;

  // ZMuStandalone counters
  int ZMuStandalone_notcontained = 0;
  int ZMuStandalone_passed = 0;
  int ZMuStandalone_MCmatched_passed = 0;

  Handle<CandMatchMap> zMuMuMap;
  if (nZMuMu > 0) {
    event.getByToken(zMuMuMapToken_, zMuMuMap);
  }

  Handle<CandMatchMap> zMuTrackMap;
  if (nZTrackMu > 0) {
    event.getByToken(zMuTrackMapToken_, zMuTrackMap);
  }

  Handle<CandMatchMap> zMuStandAloneMap;
  if (nZStandAloneMu > 0) {
    event.getByToken(zMuStandAloneMapToken_, zMuStandAloneMap);
  }

  Handle<IsolationCollection> muIso;
  event.getByToken(muIsoToken_, muIso);
  ProductID muIsoId = muIso->keyProduct().id();
  Handle<IsolationCollection> trackIso;
  event.getByToken(trackIsoToken_, trackIso);
  ProductID trackIsoId = trackIso->keyProduct().id();

  Handle<IsolationCollection> standAloneIso;
  event.getByToken(standAloneIsoToken_, standAloneIso);
  ProductID standAloneIsoId = standAloneIso->keyProduct().id();

  if (nZMuMu > 0) {
    // double mass = 1000000.;
    for (unsigned int i = 0; i < nZMuMu; i++) {
      bool ptcutAccept = false;
      bool etacutAccept = false;
      bool masscutAccept = false;
      bool isocutAccept = false;
      const Candidate& zmmCand = (*zMuMu)[i];
      CandidateRef CandRef(zMuMu, i);
      CandidateRef lep1 = zmmCand.daughter(0)->masterClone().castTo<CandidateRef>();
      CandidateRef lep2 = zmmCand.daughter(1)->masterClone().castTo<CandidateRef>();

      const double iso1 = muIso->value(lep1.key());
      const double iso2 = muIso->value(lep2.key());

      double m = zmmCand.mass();
      // check single cuts

      if (lep1->pt() > ptcut_ && lep2->pt() > ptcut_)
        ptcutAccept = true;
      if (fabs(lep1->eta()) < etacut_ && fabs(lep2->eta()) < etacut_)
        etacutAccept = true;
      if (m > minZmass_ && m < maxZmass_)
        masscutAccept = true;
      if (iso1 < isocut_ && iso2 < isocut_)
        isocutAccept = true;

      if (ptcutAccept)
        ZMuMu_ptcut_counter++;
      if (etacutAccept)
        ZMuMu_etacut_counter++;
      if (masscutAccept)
        ZMuMu_masscut_counter++;
      if (isocutAccept)
        ZMuMu_isocut_counter++;

      // check sequencial cuts

      if (ptcutAccept && etacutAccept) {
        ZMuMu_ptetacut_counter++;
        if (masscutAccept) {
          ZMuMu_ptetamasscut_counter++;
          if (isocutAccept) {
            ZMuMu_passed++;
          }
        }
      }

      if (ptcutAccept && etacutAccept && masscutAccept && isocutAccept) {
        ZMuMu_allcut_counter++;
        h_zMuMu_mass_->Fill(m);

        // check the cynematics to fill correct histograms
        for (int j = 0; j < numberOfIntervals; j++) {
          bool statusBinEta = false;
          bool statusBinPhi = false;
          bool statusBinPt = false;
          double range0 = etamin + j * binEta;
          double range1 = range0 + binEta;
          double range0phi = phimin + j * binPhi;
          double range1phi = range0phi + binPhi;
          double range0pt = ptmin + j * binPt;
          double range1pt = range0pt + binPt;
          // eta histograms
          if (lep1->eta() >= range0 && lep1->eta() < range1) {
            hmumu_eta[j]->Fill(m);
            statusBinEta = true;
          }
          if (lep2->eta() >= range0 && lep2->eta() < range1 && !statusBinEta) {
            hmumu_eta[j]->Fill(m);  // If eta1 is in the same bin of eta2 fill just once
          }
          // phi histograms
          if (lep1->phi() >= range0phi && lep1->phi() < range1phi) {
            hmumu_phi[j]->Fill(m);
            statusBinPhi = true;
          }
          if (lep2->phi() >= range0phi && lep2->phi() < range1phi && !statusBinPhi) {
            hmumu_phi[j]->Fill(m);  // If phi1 is in the same bin of phi2 fill just once
          }
          // pt histograms
          if (lep1->pt() >= range0pt && lep1->pt() < range1pt) {
            hmumu_pt[j]->Fill(m);
            statusBinPt = true;
          }
          if (lep2->pt() >= range0pt && lep2->pt() < range1pt && !statusBinPt) {
            hmumu_pt[j]->Fill(m);  // If pt1 is in the same bin of pt2 fill just once
          }
        }

        CandMatchMap::const_iterator m0 = zMuMuMap->find(CandRef);
        if (m0 != zMuMuMap->end()) {  // the Z is matched to MC thruth
          h_zMuMuMatched_mass_->Fill(m);
        }
      }
    }
  }

  h_zMuMu_numberOfCand->Fill(nZMuMu);                              // number of Z cand found per event
  h_zMuMu_numberOfCand_passed->Fill(ZMuMu_allcut_counter);         // number of Z cand after all cuts found per event
  h_zMuMu_numberOfCand_ptpassed->Fill(ZMuMu_ptcut_counter);        // number of Z cand afer pt cut found per event
  h_zMuMu_numberOfCand_etapassed->Fill(ZMuMu_etacut_counter);      // number of Z cand afer eta cut found per event
  h_zMuMu_numberOfCand_masspassed->Fill(ZMuMu_masscut_counter);    // number of Z cand afer mass cut found per event
  h_zMuMu_numberOfCand_isopassed->Fill(ZMuMu_isocut_counter);      // number of Z cand afer iso cut found per event
  h_zMuMu_numberOfCand_ptetapassed->Fill(ZMuMu_ptetacut_counter);  // number of Z cand afer pt&eta cut found per event
  h_zMuMu_numberOfCand_ptetamasspassed->Fill(
      ZMuMu_ptetamasscut_counter);  // number of Z cand afer pt&eta&mass cut found per event

  //ZmuSingleStandAlone (check MuStandalone if MuMu has not been selected by cuts)
  //  cout << "ZMuMuanalyzer : n of zMuMu " << nZMuMu << " passed " << ZMuMu_passed << "     n. of zStaMu " << nZStandAloneMu << endl;

  if (ZMuMu_passed == 0 && nZStandAloneMu > 0) {
    //      unsigned int index = 1000;
    for (unsigned int j = 0; j < nZStandAloneMu; j++) {
      const Candidate& zsmCand = (*zMuStandAlone)[j];
      bool skipZ = false;
      for (unsigned int i = 0; i < nZMuMu; i++) {  // chek if the ZMuSTandalone is contained in a ZMuMu
        const Candidate& zmmCand = (*zMuMu)[i];    // if yes .. the event has to be skipped
        if (isContained(zmmCand, zsmCand))
          skipZ = true;
      }
      if (!skipZ) {  // ZSMuSTandalone not contained in a ZMuMu
        ZMuStandalone_notcontained++;
        CandidateRef CandRef(zMuStandAlone, j);
        CandidateRef lep1 = zsmCand.daughter(0)->masterClone().castTo<CandidateRef>();
        CandidateRef lep2 = zsmCand.daughter(1)->masterClone().castTo<CandidateRef>();

        ProductID id1 = lep1.id();
        ProductID id2 = lep2.id();
        double iso1 = -1;
        double iso2 = -1;

        if (id1 == muIsoId)
          iso1 = muIso->value(lep1.key());
        else if (id1 == standAloneIsoId)
          iso1 = standAloneIso->value(lep1.key());

        if (id2 == muIsoId)
          iso2 = muIso->value(lep2.key());
        else if (id2 == standAloneIsoId)
          iso2 = standAloneIso->value(lep2.key());

        double ms = zsmCand.mass();
        if (lep1->pt() > ptSTAcut_ && lep2->pt() > ptSTAcut_ && fabs(lep1->eta()) < etacut_ &&
            fabs(lep2->eta()) < etacut_ && ms > minZmass_ && ms < maxZmass_ && iso1 < isocut_ && iso2 < isocut_) {
          h_zMuSingleStandAlone_mass_->Fill(ms);
          ZMuStandalone_passed++;
          // check the cynematics to fill correct histograms
          for (int j = 0; j < numberOfIntervals; j++) {
            double range0 = etamin + j * binEta;
            double range1 = range0 + binEta;
            double range0phi = phimin + j * binPhi;
            double range1phi = range0phi + binPhi;
            double range0pt = ptmin + j * binPt;
            double range1pt = range0pt + binPt;

            // check which muon is a standalone (standalone means that there is a reference missing.)
            if ((lep1->get<TrackRef, reco::StandAloneMuonTag>()).isNull()) {
              if (lep1->eta() >= range0 && lep1->eta() < range1)
                hmusta_eta[j]->Fill(ms);
              if (lep1->phi() >= range0phi && lep1->phi() < range1phi)
                hmusta_phi[j]->Fill(ms);
              if (lep1->pt() >= range0pt && lep1->pt() < range1pt)
                hmusta_pt[j]->Fill(ms);
            }
            if ((lep2->get<TrackRef, reco::StandAloneMuonTag>()).isNull()) {
              if (lep2->eta() >= range0 && lep2->eta() < range1)
                hmusta_eta[j]->Fill(ms);
              if (lep2->phi() >= range0phi && lep2->phi() < range1phi)
                hmusta_phi[j]->Fill(ms);
              if (lep2->pt() >= range0pt && lep2->pt() < range1pt)
                hmusta_pt[j]->Fill(ms);
            }
          }
          CandMatchMap::const_iterator m0 = zMuStandAloneMap->find(CandRef);
          if (m0 != zMuStandAloneMap->end()) {
            ZMuStandalone_MCmatched_passed++;
            h_zMuSingleStandAloneMatched_mass_->Fill(ms);
          }
        }
      }
    }
    h_zMuSta_numberOfCand->Fill(
        nZStandAloneMu);  // number of ZMuStandalone cand found per event (no higher priority Z selected)
    h_zMuSta_numberOfCand_notcontained->Fill(ZMuStandalone_notcontained);
    h_zMuSta_numberOfCand_passed->Fill(
        ZMuStandalone_passed);  // number of ZMuSTa cand after all cuts found per event (no higher prioriy Z selected)
    h_zMuSta_MCmatched_numberOfCand_passed->Fill(
        ZMuStandalone_MCmatched_passed);  // number of ZMuSTa MC matched cand after all cuts found per event (no higher prioriy Z selected)
  }

  //ZmuSingleTRack  (check MuTrack if MuMu has not been selected)
  if (ZMuMu_passed == 0 && ZMuStandalone_passed == 0 && nZTrackMu > 0) {
    for (unsigned int j = 0; j < nZTrackMu; j++) {
      const Candidate& ztmCand = (*zMuTrack)[j];
      bool skipZ = false;
      for (unsigned int i = 0; i < nZMuMu; i++) {  // chek if the ZMuTrack is contained in a ZMuMu
        const Candidate& zmmCand = (*zMuMu)[i];    // if yes .. the event has to be skipped
        if (isContained(zmmCand, ztmCand))
          skipZ = true;
      }
      if (!skipZ) {
        ZMuTrack_notcontained++;
        CandidateRef CandRef(zMuTrack, j);
        CandidateRef lep1 = ztmCand.daughter(0)->masterClone().castTo<CandidateRef>();
        CandidateRef lep2 = ztmCand.daughter(1)->masterClone().castTo<CandidateRef>();

        ProductID id1 = lep1.id();
        ProductID id2 = lep2.id();
        double iso1 = -1;
        double iso2 = -1;

        if (id1 == muIsoId)
          iso1 = muIso->value(lep1.key());
        else if (id1 == trackIsoId)
          iso1 = trackIso->value(lep1.key());

        if (id2 == muIsoId)
          iso2 = muIso->value(lep2.key());
        else if (id2 == trackIsoId)
          iso2 = trackIso->value(lep2.key());

        double mt = ztmCand.mass();
        if (lep1->pt() > ptcut_ && lep2->pt() > ptcut_ && fabs(lep1->eta()) < etacut_ && fabs(lep2->eta()) < etacut_ &&
            mt > minZmass_ && mt < maxZmass_ && iso1 < isocut_ && iso2 < isocut_) {
          h_zMuSingleTrack_mass_->Fill(mt);
          ZMuTrack_passed++;

          // check the cynematics to fill correct histograms
          for (int j = 0; j < numberOfIntervals; j++) {
            double range0 = etamin + j * binEta;
            double range1 = range0 + binEta;
            double range0phi = phimin + j * binPhi;
            double range1phi = range0phi + binPhi;
            double range0pt = ptmin + j * binPt;
            double range1pt = range0pt + binPt;

            // check which muon is a track only (track only means that there is a reference missing.)
            if ((lep1->get<TrackRef, reco::StandAloneMuonTag>()).isNull()) {
              if (lep1->eta() >= range0 && lep1->eta() < range1)
                hmutrack_eta[j]->Fill(mt);
              if (lep1->phi() >= range0phi && lep1->phi() < range1phi)
                hmutrack_phi[j]->Fill(mt);
              if (lep1->pt() >= range0pt && lep1->pt() < range1pt)
                hmutrack_pt[j]->Fill(mt);
            }
            if ((lep2->get<TrackRef, reco::StandAloneMuonTag>()).isNull()) {
              if (lep2->eta() >= range0 && lep2->eta() < range1)
                hmutrack_eta[j]->Fill(mt);
              if (lep2->phi() >= range0phi && lep2->phi() < range1phi)
                hmutrack_phi[j]->Fill(mt);
              if (lep2->pt() >= range0pt && lep2->pt() < range1pt)
                hmutrack_pt[j]->Fill(mt);
            }
          }
          CandMatchMap::const_iterator m0 = zMuTrackMap->find(CandRef);
          if (m0 != zMuTrackMap->end()) {
            ZMuTrack_MCmatched_passed++;
            h_zMuSingleTrackMatched_mass_->Fill(mt);
          }
        }
      }
    }
    h_zMuTrack_numberOfCand->Fill(
        nZTrackMu);  // number of ZMuTrack cand found per event (no higher priority Z selected)
    h_zMuTrack_numberOfCand_notcontained->Fill(
        ZMuTrack_notcontained);  // number of ZMuTrack cand not cntained in ZMuMu (no higher priority Z selected)

    h_zMuTrack_numberOfCand_passed->Fill(
        ZMuTrack_passed);  // number of ZMuTrack cand after all cuts found per event (no higher priority Z selected)

    h_zMuTrack_MCmatched_numberOfCand_passed->Fill(ZMuTrack_MCmatched_passed);
  }
}

void ZMuMuAnalyzer_cynematics::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 690 of file ZMuMuAnalyzer_cynematics.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, h_zMuMu_mass_, h_zMuMu_numberOfCand, h_zMuMu_numberOfCand_etapassed, h_zMuMu_numberOfCand_isopassed, h_zMuMu_numberOfCand_masspassed, h_zMuMu_numberOfCand_passed, h_zMuMu_numberOfCand_ptetamasspassed, h_zMuMu_numberOfCand_ptetapassed, h_zMuMu_numberOfCand_ptpassed, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, h_zMuSta_MCmatched_numberOfCand_passed, h_zMuSta_numberOfCand, h_zMuSta_numberOfCand_notcontained, h_zMuSta_numberOfCand_passed, h_zMuTrack_MCmatched_numberOfCand_passed, h_zMuTrack_numberOfCand, h_zMuTrack_numberOfCand_notcontained, and h_zMuTrack_numberOfCand_passed.

                                      {
  // candidate analysis
  // ZMuMu
  double Nzmmc = h_zMuMu_numberOfCand->GetEntries();
  double Nzmmc_0Z = h_zMuMu_numberOfCand->GetBinContent(1);
  double Nzmmc_1Z = h_zMuMu_numberOfCand->GetBinContent(2);
  double Nzmmc_moreZ = Nzmmc - Nzmmc_0Z - Nzmmc_1Z;
  double Nzmmc_passed_0Z = h_zMuMu_numberOfCand_passed->GetBinContent(1);
  double Nzmmc_passed_1Z = h_zMuMu_numberOfCand_passed->GetBinContent(2);
  double Nzmmc_passed_moreZ = Nzmmc - Nzmmc_passed_0Z - Nzmmc_passed_1Z;
  double Nzmmc_ptpassed_0Z = h_zMuMu_numberOfCand_ptpassed->GetBinContent(1);
  double Nzmmc_ptpassed_1Z = h_zMuMu_numberOfCand_ptpassed->GetBinContent(2);
  double Nzmmc_etapassed_0Z = h_zMuMu_numberOfCand_etapassed->GetBinContent(1);
  double Nzmmc_etapassed_1Z = h_zMuMu_numberOfCand_etapassed->GetBinContent(2);
  double Nzmmc_masspassed_0Z = h_zMuMu_numberOfCand_masspassed->GetBinContent(1);
  double Nzmmc_masspassed_1Z = h_zMuMu_numberOfCand_masspassed->GetBinContent(2);
  double Nzmmc_isopassed_0Z = h_zMuMu_numberOfCand_isopassed->GetBinContent(1);
  double Nzmmc_isopassed_1Z = h_zMuMu_numberOfCand_isopassed->GetBinContent(2);
  double Nzmmc_ptetapassed_0Z = h_zMuMu_numberOfCand_ptetapassed->GetBinContent(1);
  double Nzmmc_ptetapassed_1Z = h_zMuMu_numberOfCand_ptetapassed->GetBinContent(2);
  double Nzmmc_ptetamasspassed_0Z = h_zMuMu_numberOfCand_ptetamasspassed->GetBinContent(1);
  double Nzmmc_ptetamasspassed_1Z = h_zMuMu_numberOfCand_ptetamasspassed->GetBinContent(2);
  double Nzmmc_ptpassed_moreZ = Nzmmc - Nzmmc_ptpassed_0Z - Nzmmc_ptpassed_1Z;
  double Nzmmc_etapassed_moreZ = Nzmmc - Nzmmc_etapassed_0Z - Nzmmc_etapassed_1Z;
  double Nzmmc_masspassed_moreZ = Nzmmc - Nzmmc_masspassed_0Z - Nzmmc_masspassed_1Z;
  double Nzmmc_isopassed_moreZ = Nzmmc - Nzmmc_isopassed_0Z - Nzmmc_isopassed_1Z;
  double Nzmmc_ptetapassed_moreZ = Nzmmc - Nzmmc_ptetapassed_0Z - Nzmmc_ptetapassed_1Z;
  double Nzmmc_ptetamasspassed_moreZ = Nzmmc - Nzmmc_ptetamasspassed_0Z - Nzmmc_ptetamasspassed_1Z;
  double Nzmsc = h_zMuSta_numberOfCand->GetEntries();
  double Nzmsc_0Z = h_zMuSta_numberOfCand->GetBinContent(1);
  double Nzmsc_1Z = h_zMuSta_numberOfCand->GetBinContent(2);
  double Nzmsc_moreZ = Nzmsc - Nzmsc_0Z - Nzmsc_1Z;
  double Nzmsc_notcontained_0Z = h_zMuSta_numberOfCand_notcontained->GetBinContent(1);
  double Nzmsc_notcontained_1Z = h_zMuSta_numberOfCand_notcontained->GetBinContent(2);
  double Nzmsc_notcontained_moreZ = Nzmsc - Nzmsc_notcontained_0Z - Nzmsc_notcontained_1Z;
  double Nzmsc_passed_0Z = h_zMuSta_numberOfCand_passed->GetBinContent(1);
  double Nzmsc_passed_1Z = h_zMuSta_numberOfCand_passed->GetBinContent(2);
  double Nzmsc_passed_moreZ = Nzmsc - Nzmsc_passed_0Z - Nzmsc_passed_1Z;
  double Nzmsc_MCmatched_passed_0Z = h_zMuSta_MCmatched_numberOfCand_passed->GetBinContent(1);
  double Nzmsc_MCmatched_passed_1Z = h_zMuSta_MCmatched_numberOfCand_passed->GetBinContent(2);
  double Nzmsc_MCmatched_passed_moreZ = Nzmsc - Nzmsc_MCmatched_passed_0Z - Nzmsc_MCmatched_passed_1Z;
  double Nzmtc = h_zMuTrack_numberOfCand->GetEntries();
  double Nzmtc_0Z = h_zMuTrack_numberOfCand->GetBinContent(1);
  double Nzmtc_1Z = h_zMuTrack_numberOfCand->GetBinContent(2);
  double Nzmtc_moreZ = Nzmtc - Nzmtc_0Z - Nzmtc_1Z;
  double Nzmtc_notcontained_0Z = h_zMuTrack_numberOfCand_notcontained->GetBinContent(1);
  double Nzmtc_notcontained_1Z = h_zMuTrack_numberOfCand_notcontained->GetBinContent(2);
  double Nzmtc_notcontained_moreZ = Nzmtc - Nzmtc_notcontained_0Z - Nzmtc_notcontained_1Z;
  double Nzmtc_passed_0Z = h_zMuTrack_numberOfCand_passed->GetBinContent(1);
  double Nzmtc_passed_1Z = h_zMuTrack_numberOfCand_passed->GetBinContent(2);
  double Nzmtc_passed_moreZ = Nzmtc - Nzmtc_passed_0Z - Nzmtc_passed_1Z;
  double Nzmtc_MCmatched_passed_0Z = h_zMuTrack_MCmatched_numberOfCand_passed->GetBinContent(1);
  double Nzmtc_MCmatched_passed_1Z = h_zMuTrack_MCmatched_numberOfCand_passed->GetBinContent(2);
  double Nzmtc_MCmatched_passed_moreZ = Nzmtc - Nzmtc_MCmatched_passed_0Z - Nzmtc_MCmatched_passed_1Z;

  cout << "--------------- Statistics ----------------------------------------------------------" << endl;
  cout << "n of ZMuMu entries   ...................................................... " << Nzmmc << endl;
  cout << "n of ZMuMu events with 0 cand ............................................. " << Nzmmc_0Z << endl;
  cout << "n of ZMuMu events with 1 cand ............................................. " << Nzmmc_1Z << endl;
  cout << "n of ZMuMu events with 2 or more cand ..................................... " << Nzmmc_moreZ << endl << endl;

  cout << "n of ZMuMu events not selected by cuts .................................... " << Nzmmc_passed_0Z << endl;
  cout << "n of ZMuMu events with 1 cand selected by cuts ............................ " << Nzmmc_passed_1Z << endl;
  cout << "n of ZMuMu events with 2 or more cand elected by cuts ..................... " << Nzmmc_passed_moreZ << endl
       << endl;

  cout << "n of ZMuMu events not selected by pt cut .................................. " << Nzmmc_ptpassed_0Z << endl;
  cout << "n of ZMuMu events with 1 cand selected by pt cut .......................... " << Nzmmc_ptpassed_1Z << endl;
  cout << "n of ZMuMu events with 2 or more cand elected by pt cut ................... " << Nzmmc_ptpassed_moreZ << endl
       << endl;

  cout << "n of ZMuMu events not selected by eta cut ................................. " << Nzmmc_etapassed_0Z << endl;
  cout << "n of ZMuMu events with 1 cand selected by eta cut ......................... " << Nzmmc_etapassed_1Z << endl;
  cout << "n of ZMuMu events with 2 or more cand elected by eta cut .................. " << Nzmmc_etapassed_moreZ
       << endl
       << endl;

  cout << "n of ZMuMu events not selected by mass cut ................................ " << Nzmmc_masspassed_0Z << endl;
  cout << "n of ZMuMu events with 1 cand selected by mass cut ........................ " << Nzmmc_masspassed_1Z << endl;
  cout << "n of ZMuMu events with 2 or more cand elected by mass cut ................. " << Nzmmc_masspassed_moreZ
       << endl
       << endl;

  cout << "n of ZMuMu events not selected by iso cut ................................. " << Nzmmc_isopassed_0Z << endl;
  cout << "n of ZMuMu events with 1 cand selected iso cut ............................ " << Nzmmc_isopassed_1Z << endl;
  cout << "n of ZMuMu events with 2 or more cand elected iso cut ..................... " << Nzmmc_isopassed_moreZ
       << endl
       << endl;

  cout << "n of ZMuMu events not selected by pt and eta cut .......................... " << Nzmmc_ptetapassed_0Z
       << endl;
  cout << "n of ZMuMu events with 1 cand selected by pt and eta cut .................. " << Nzmmc_ptetapassed_1Z
       << endl;
  cout << "n of ZMuMu events with 2 or more cand elected by pt and eta cut ........... " << Nzmmc_ptetapassed_moreZ
       << endl
       << endl;

  cout << "n of ZMuMu events not selected by pt and eta and mass cut ................. " << Nzmmc_ptetamasspassed_0Z
       << endl;
  cout << "n of ZMuMu events with 1 cand selected by pt and eta and mass cut ......... " << Nzmmc_ptetamasspassed_1Z
       << endl;
  cout << "n of ZMuMu events with 2 or more cand elected by pt and eta and mass cut .. " << Nzmmc_ptetamasspassed_moreZ
       << endl
       << endl;

  cout << "................When No ZMuMu are selected.................................." << endl;
  cout << "n of ZMuSta entries ....................................................... " << Nzmsc << endl;
  cout << "n of ZMuSta events with 0 cand ............................................ " << Nzmsc_0Z << endl;
  cout << "n of ZMuSta events with 1 cand ............................................ " << Nzmsc_1Z << endl;
  cout << "n of ZMuSta events with 2 or more cand .................................... " << Nzmsc_moreZ << endl << endl;

  cout << "n of ZMuSta not contained events with 0 cand .............................. " << Nzmsc_notcontained_0Z
       << endl;
  cout << "n of ZMuSta events not contained with 1 cand .............................. " << Nzmsc_notcontained_1Z
       << endl;
  cout << "n of ZMuSta events no contained with 2 or more cand ....................... " << Nzmsc_notcontained_moreZ
       << endl
       << endl;

  cout << "n of ZMuSta cand not selectd by cuts ...................................... " << Nzmsc_passed_0Z << endl;
  cout << "n of ZMuSta events with 1 cand selected by cuts ........................... " << Nzmsc_passed_1Z << endl;
  cout << "n of ZMuSta events with 2 or more cand selected by cuts ................... " << Nzmsc_passed_moreZ << endl
       << endl;

  cout << "n of ZMuSta MCmatched cand not selectd by cuts ............................ " << Nzmsc_MCmatched_passed_0Z
       << endl;
  cout << "n of ZMuSta MCmatched events with 1 cand selected by cuts ................. " << Nzmsc_MCmatched_passed_1Z
       << endl;
  cout << "n of ZMuSta MCmatched events with 2 or more cand selected by cuts ......... " << Nzmsc_MCmatched_passed_moreZ
       << endl
       << endl;

  cout << "...............When no ZMuMu and ZMuSta are selcted........................." << endl;
  cout << "n of ZMuTrack entries   ................................................... " << Nzmtc << endl;
  cout << "n of ZMuTrack events with 0 cand ..........................................." << Nzmtc_0Z << endl;
  cout << "n of ZMuTrack events with 1 cand ..........................................." << Nzmtc_1Z << endl;
  cout << "n of ZMuTrack events with 2 or more cand ..................................." << Nzmtc_moreZ << endl << endl;

  cout << "n of ZMuTrack not contained events with 0 cand ............................ " << Nzmtc_notcontained_0Z
       << endl;
  cout << "n of ZMuTrack events not contained with 1 cand ............................ " << Nzmtc_notcontained_1Z
       << endl;
  cout << "n of ZMuTrack events no contained with 2 or more cand ..................... " << Nzmtc_notcontained_moreZ
       << endl
       << endl;

  cout << "n of ZMuTrack cand not selectd by cuts ....................................." << Nzmtc_passed_0Z << endl;
  cout << "n of ZMuTrack events with 1 cand selected by cuts .........................." << Nzmtc_passed_1Z << endl;
  cout << "n of ZMuTrack events with 2 or more cand selected by cuts .................." << Nzmtc_passed_moreZ << endl
       << endl;

  cout << "n of ZMuTrack MCmatched cand not selectd by cuts .......................... " << Nzmtc_MCmatched_passed_0Z
       << endl;
  cout << "n of ZMuTrcak MCmatched events with 1 cand selected by cuts ............... " << Nzmtc_MCmatched_passed_1Z
       << endl;
  cout << "n of ZMuTrack MCmatched events with 2 or more cand selected by cuts ....... " << Nzmtc_MCmatched_passed_moreZ
       << endl;

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

  double Nzmm = h_zMuMu_mass_->GetEntries();
  double Nzsm = h_zMuSingleStandAlone_mass_->GetEntries();
  double Nzsnom = h_zMuSingleStandAloneOverlap_mass_->GetEntries();
  double Nztm = h_zMuSingleTrack_mass_->GetEntries();

  double NzmmMatch = h_zMuMuMatched_mass_->GetEntries();
  double NzsmMatch = h_zMuSingleStandAloneMatched_mass_->GetEntries();
  double NzsnomMatch = h_zMuSingleStandAloneOverlapMatched_mass_->GetEntries();
  double NztmMatch = h_zMuSingleTrackMatched_mass_->GetEntries();

  cout << "-- N SingleTrackMu = " << Nztm << endl;
  cout << "-----N SinglStandAloneMu = " << Nzsm << endl;
  cout << "-----N SingleStandAloneOverlapMu = " << Nzsnom << endl;
  cout << "------- N MuMu = " << Nzmm << endl;

  cout << "-- N SingleTrackMuMatched = " << NztmMatch << endl;
  cout << "-----N SinglStandAloneMuMatched = " << NzsmMatch << endl;
  cout << "-----N SingleStandAloneOverlapMuMatched  = " << NzsnomMatch << endl;
  cout << "------- N MuMu Matched  = " << NzmmMatch << endl;
}

bool ZMuMuAnalyzer_cynematics::isContained ( const Candidate &  obj1, const Candidate &  obj2  )

private

Definition at line 592 of file ZMuMuAnalyzer_cynematics.cc.

References gather_cfg::cout, reco::Candidate::daughter(), RemoveAddSevLevel::flag, reco::Candidate::get(), mps_fire::i, edm::Ref< C, T, F >::id(), dqmiolumiharvest::j, crabWrapper::key, edm::Ref< C, T, F >::key(), muonTagProbeFilters_cff::matched, and reco::Candidate::numberOfDaughters().

Referenced by analyze().

                                                                                       {
  // check if a candidate obj2 is different from obj1  (assume that obj1 is a ZMuMu and obj2 is any other type)
  // (for example a Z can be done with two global muons, or with a global muon plus a standalone muon.
  // if the standalone muon is part of the second global muon in fact this is the same Z)

  const int maxd = 10;
  const Candidate* daughters1[maxd];
  const Candidate* daughters2[maxd];
  TrackRef trackerTrack1[maxd];
  TrackRef stAloneTrack1[maxd];
  TrackRef globalTrack1[maxd];
  TrackRef trackerTrack2[maxd];
  TrackRef stAloneTrack2[maxd];
  TrackRef globalTrack2[maxd];
  bool flag;
  unsigned int nd1 = obj1.numberOfDaughters();
  unsigned int nd2 = obj2.numberOfDaughters();
  unsigned int matched = 0;

  for (unsigned int i = 0; i < nd1; ++i) {
    daughters1[i] = obj1.daughter(i);
    trackerTrack1[i] = daughters1[i]->get<TrackRef>();
    stAloneTrack1[i] = daughters1[i]->get<TrackRef, reco::StandAloneMuonTag>();
    globalTrack1[i] = daughters1[i]->get<TrackRef, reco::CombinedMuonTag>();

    /*********************************************** just used for debug ********************
    if (trackerTrack1[i].isNull())
      cout << "in ZMuMu daughter " << i << " tracker ref non found " << endl;
    else
      cout << "in ZMuMu daughter " << i << " tracker ref FOUND"
       << " id: " << trackerTrack1[i].id() << ", index: " << trackerTrack1[i].key()
       << endl;
    if (stAloneTrack1[i].isNull())
      cout << "in ZMuMu daughter " << i << " stalone ref non found " << endl;
    else
      cout << "in ZMuMu daughter " << i << " stalone ref FOUND"
       << " id: " << stAloneTrack1[i].id() << ", index: " << stAloneTrack1[i].key()
       << endl;

    if (globalTrack1[i].isNull())
      cout << "in ZMuMu daughter " << i << " global ref non found " << endl;
    else
      cout << "in ZMuMu daughter " << i << " global ref FOUND"
       << " id: " << globalTrack1[i].id() << ", index: " << globalTrack1[i].key()
       << endl;
    */
  }
  for (unsigned int i = 0; i < nd2; ++i) {
    daughters2[i] = obj2.daughter(i);
    trackerTrack2[i] = daughters2[i]->get<TrackRef>();
    stAloneTrack2[i] = daughters2[i]->get<TrackRef, reco::StandAloneMuonTag>();
    globalTrack2[i] = daughters2[i]->get<TrackRef, reco::CombinedMuonTag>();

    /******************************************** just used for debug ************
    if (trackerTrack2[i].isNull())
      cout << "in ZMuSta daughter " << i << " tracker ref non found " << endl;
    else
      cout << "in ZMuSta daughter " << i << " tracker ref FOUND"
       << " id: " << trackerTrack2[i].id() << ", index: " << trackerTrack2[i].key()
       << endl;
    if (stAloneTrack2[i].isNull())
      cout << "in ZMuSta daughter " << i << " standalone ref non found " << endl;
    else
      cout << "in ZMuSta daughter " << i << " standalone ref FOUND"
       << " id: " << stAloneTrack2[i].id() << ", index: " << stAloneTrack2[i].key()
       << endl;

    if (globalTrack2[i].isNull())
      cout << "in ZMuSta daughter " << i << " global ref non found " << endl;
    else
      cout << "in ZMuSta daughter " << i << " global ref FOUND"
       << " id: " << globalTrack2[i].id() << ", index: " << globalTrack2[i].key()
       << endl;

   */
  }
  if (nd1 != nd2) {
    cout << "ZMuMuAnalyzer::isContained WARNING n.of daughters different " << nd1 << "  " << nd2 << endl;
  } else {
    for (unsigned int i = 0; i < nd1; i++) {
      flag = false;
      for (unsigned int j = 0; j < nd2;
           j++) {  // if the obj2 is a standalone the trackref is alwais in the trackerTRack position
        if (((trackerTrack2[i].id() == trackerTrack1[j].id()) && (trackerTrack2[i].key() == trackerTrack1[j].key())) ||
            ((trackerTrack2[i].id() == stAloneTrack1[j].id()) && (trackerTrack2[i].key() == stAloneTrack1[j].key()))) {
          flag = true;
        }
      }
      if (flag)
        matched++;
    }
  }
  if (matched == nd1)  // return true if all the childrens of the ZMuMu have a children matched in ZMuXX
    return true;
  else
    return false;
}

Member Data Documentation

double ZMuMuAnalyzer_cynematics::binEta

private

Definition at line 82 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::binPhi

private

Definition at line 82 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::binPt

private

Definition at line 82 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::etacut_

private

Definition at line 68 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::etamax

private

Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::etamin

private

Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

TH2D* ZMuMuAnalyzer_cynematics::h_OneSta_mass

private

Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by ZMuMuAnalyzer_cynematics().

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D* ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_etapassed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_isopassed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_masspassed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_passed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_ptetamasspassed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_ptetapassed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_ptpassed

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuMuMatched_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAlone_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAloneMatched_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAloneOverlap_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by endJob(), and ZMuMuAnalyzer_cynematics().

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAloneOverlapMatched_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by endJob(), and ZMuMuAnalyzer_cynematics().

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleTrack_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleTrackMatched_mass_

private

Definition at line 69 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSta_MCmatched_numberOfCand_passed

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D* ZMuMuAnalyzer_cynematics::h_zMuSta_numberOfCand

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSta_numberOfCand_notcontained

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuSta_numberOfCand_passed

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_MCmatched_numberOfCand_passed

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_numberOfCand

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_numberOfCand_notcontained

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_numberOfCand_passed

private

Definition at line 75 of file ZMuMuAnalyzer_cynematics.cc.

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

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmumu_eta

private

Definition at line 83 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmumu_phi

private

Definition at line 84 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmumu_pt

private

Definition at line 85 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmusta_eta

private

Definition at line 83 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmusta_phi

private

Definition at line 84 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmusta_pt

private

Definition at line 85 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmutrack_eta

private

Definition at line 83 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmutrack_phi

private

Definition at line 84 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

vector<TH1D*> ZMuMuAnalyzer_cynematics::hmutrack_pt

private

Definition at line 85 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::isocut_

private

Definition at line 68 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

double ZMuMuAnalyzer_cynematics::maxZmass_

private

Definition at line 68 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

double ZMuMuAnalyzer_cynematics::minZmass_

private

Definition at line 68 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<IsolationCollection> ZMuMuAnalyzer_cynematics::muIsoToken_

private

Definition at line 62 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

int ZMuMuAnalyzer_cynematics::numberOfIntervals

private

Definition at line 81 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

OverlapChecker ZMuMuAnalyzer_cynematics::overlap_

private

Definition at line 58 of file ZMuMuAnalyzer_cynematics.cc.

double ZMuMuAnalyzer_cynematics::phimax

private

Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::phimin

private

Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::ptcut_

private

Definition at line 68 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::ptmax

private

Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::ptmin

private

Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze(), and ZMuMuAnalyzer_cynematics().

double ZMuMuAnalyzer_cynematics::ptSTAcut_

private

Definition at line 68 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<IsolationCollection> ZMuMuAnalyzer_cynematics::standAloneIsoToken_

private

Definition at line 64 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<IsolationCollection> ZMuMuAnalyzer_cynematics::trackIsoToken_

private

Definition at line 63 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<CandMatchMap> ZMuMuAnalyzer_cynematics::zMuMuMapToken_

private

Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<CandidateCollection> ZMuMuAnalyzer_cynematics::zMuMuToken_

private

Definition at line 59 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<CandMatchMap> ZMuMuAnalyzer_cynematics::zMuStandAloneMapToken_

private

Definition at line 67 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<CandidateCollection> ZMuMuAnalyzer_cynematics::zMuStandAloneToken_

private

Definition at line 61 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<CandMatchMap> ZMuMuAnalyzer_cynematics::zMuTrackMapToken_

private

Definition at line 66 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().

EDGetTokenT<CandidateCollection> ZMuMuAnalyzer_cynematics::zMuTrackToken_

private

Definition at line 60 of file ZMuMuAnalyzer_cynematics.cc.

Referenced by analyze().