testAnalyzer Class Reference

Inheritance diagram for testAnalyzer:

Public Member Functions
	testAnalyzer (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

Private Attributes
int	DoubleTrigger_
TH1D *	h_EtaDist_Pt80_
TH1D *	h_numberMuon_
TH1D *	h_numberMuon_ptStudy_
TH1D *	h_numberTrigMuon_
TH1D *	h_numberTrigMuon_ptStudy_
TH1D *	h_pt_distribution_
double	maxEta_
double	maxPt_
double	minEta_
double	minPt_
int	nbinsEta_
int	nbinsEtaPt_
int	nbinsPt_
int	NoTrigger_
string	pathName_
EDGetTokenT< vector< pat::Muon > >	selectMuonToken_
int	SingleTrigger_
vector< double >	vectorEta
vector< double >	vectorPt
int	zmumuIncrement_
EDGetTokenT< CandidateView >	zMuMuToken_

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 31 of file triggerEfficiency.cc.

Constructor & Destructor Documentation

testAnalyzer::testAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 54 of file triggerEfficiency.cc.

References DoubleTrigger_, h_EtaDist_Pt80_, h_numberMuon_, h_numberMuon_ptStudy_, h_numberTrigMuon_, h_numberTrigMuon_ptStudy_, h_pt_distribution_, TFileService::make(), maxEta_, maxPt_, minEta_, minPt_, nbinsEta_, nbinsEtaPt_, nbinsPt_, NoTrigger_, SingleTrigger_, and zmumuIncrement_.

    : selectMuonToken_(consumes<vector<pat::Muon> >(pset.getParameter<InputTag>("selectMuon"))),
      zMuMuToken_(consumes<CandidateView>(pset.getParameter<InputTag>("ZMuMu"))),
      pathName_(pset.getParameter<string>("pathName")),
      nbinsEta_(pset.getParameter<int>("EtaBins")),
      minEta_(pset.getParameter<double>("minEta")),
      maxEta_(pset.getParameter<double>("maxEta")),
      nbinsPt_(pset.getParameter<int>("PtBins")),
      minPt_(pset.getParameter<double>("minPt")),
      maxPt_(pset.getParameter<double>("maxPt")),
      nbinsEtaPt_(pset.getParameter<int>("EtaPt80Bins")) {
  SingleTrigger_ = 0;
  DoubleTrigger_ = 0;
  NoTrigger_ = 0;
  zmumuIncrement_ = 0;
  Service<TFileService> fs;
  h_pt_distribution_ = fs->make<TH1D>("PtResolution ", "Pt Resolution", 200, -4., 4.);
  h_numberMuon_ = fs->make<TH1D>("Denominatore", "Number of Muons vs Eta", nbinsEta_, minEta_, maxEta_);
  h_numberTrigMuon_ =
      fs->make<TH1D>("NumeratoreTrigMuon", "Number of Triggered Muons vs Eta", nbinsEta_, minEta_, maxEta_);
  h_numberMuon_ptStudy_ = fs->make<TH1D>("DenominatorePtStudy", "Number of Muons vs Pt", nbinsPt_, minPt_, maxPt_);
  h_numberTrigMuon_ptStudy_ =
      fs->make<TH1D>("NumeratoreTrigMuonPtStudy", "Number of Triggered Muons vs Pt", nbinsPt_, minPt_, maxPt_);
  h_EtaDist_Pt80_ = fs->make<TH1D>("EtaDistr", "Eta distribution (Pt>80)", nbinsEtaPt_, minEta_, maxEta_);
}

Member Function Documentation

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

overrideprivate

Definition at line 80 of file triggerEfficiency.cc.

References reco::Candidate::daughter(), DoubleTrigger_, reco::LeafCandidate::eta(), h_EtaDist_Pt80_, h_numberMuon_, h_numberMuon_ptStudy_, h_numberTrigMuon_, h_numberTrigMuon_ptStudy_, h_pt_distribution_, mps_fire::i, dqmiolumiharvest::j, reco::Candidate::masterClone(), NoTrigger_, pathName_, reco::LeafCandidate::pt(), selectMuonToken_, SingleTrigger_, edm::View< T >::size(), pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), vectorEta, vectorPt, ZMuMuAnalysisNtupler_cff::zMuMu, zmumuIncrement_, and zMuMuToken_.

                                                                            {
  Handle<vector<pat::Muon> > selectMuon;
  event.getByToken(selectMuonToken_, selectMuon);
  Handle<CandidateView> zMuMu;
  event.getByToken(zMuMuToken_, zMuMu);
  int zmumuSize = zMuMu->size();
  if (zmumuSize > 0) {
    for (int i = 0; i < zmumuSize; ++i) {
      bool singleTrigFlag0 = false;
      bool singleTrigFlag1 = false;
      zmumuIncrement_++;
      const Candidate& zMuMuCand = (*zMuMu)[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);
      const pat::TriggerObjectStandAloneCollection mu0HLTMatches = mu0.triggerObjectMatchesByPath(pathName_);
      const pat::TriggerObjectStandAloneCollection mu1HLTMatches = mu1.triggerObjectMatchesByPath(pathName_);
      double EtaPatMu0 = mu0.eta();
      double EtaPatMu1 = mu1.eta();
      double PtPatMu0 = mu0.pt();
      double PtPatMu1 = mu1.pt();
      h_numberMuon_->Fill(EtaPatMu0);
      h_numberMuon_->Fill(EtaPatMu1);
      h_numberMuon_ptStudy_->Fill(PtPatMu0);
      h_numberMuon_ptStudy_->Fill(PtPatMu1);
      int dimTrig0 = mu0HLTMatches.size();
      int dimTrig1 = mu1HLTMatches.size();
      if (dimTrig0 != 0) {
        for (int j = 0; j < dimTrig0; ++j) {
          singleTrigFlag0 = true;
          h_numberTrigMuon_->Fill(EtaPatMu0);
          h_numberTrigMuon_ptStudy_->Fill(PtPatMu0);
          double PtTrig = mu0HLTMatches[j].pt();
          double PtDif = PtTrig - PtPatMu0;
          h_pt_distribution_->Fill(PtDif);
        }
      } else {
        if (PtPatMu0 > 80) {
          h_EtaDist_Pt80_->Fill(EtaPatMu0);
          vectorPt.push_back(PtPatMu0);
          vectorEta.push_back(EtaPatMu0);
        }
      }
      if (dimTrig1 != 0) {
        for (int j = 0; j < dimTrig1; ++j) {
          singleTrigFlag1 = true;
          h_numberTrigMuon_->Fill(EtaPatMu1);
          h_numberTrigMuon_ptStudy_->Fill(PtPatMu1);
          double PtTrig = mu0HLTMatches[j].pt();
          double PtDif = PtTrig - PtPatMu1;
          h_pt_distribution_->Fill(PtDif);
        }
      } else {
        if (PtPatMu0 > 80) {
          h_EtaDist_Pt80_->Fill(EtaPatMu1);
          vectorPt.push_back(PtPatMu0);
          vectorEta.push_back(EtaPatMu0);
        }
      }

      if (singleTrigFlag0 && singleTrigFlag1)
        DoubleTrigger_++;
      if (singleTrigFlag0 && !singleTrigFlag1)
        SingleTrigger_++;
      if (!singleTrigFlag0 && singleTrigFlag1)
        SingleTrigger_++;
      if (!singleTrigFlag0 && !singleTrigFlag1)
        NoTrigger_++;

    }  //end loop on ZMuMu candidates
  }    //end check on ZMuMu
}

void testAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 154 of file triggerEfficiency.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, DoubleTrigger_, mps_fire::i, NoTrigger_, SingleTrigger_, vectorEta, vectorPt, and zmumuIncrement_.

                          {
  cout << "DoubleTrigger = " << DoubleTrigger_ << " , SingleTrigger = " << SingleTrigger_
       << " , NoTrigger = " << NoTrigger_ << " ,zmumuIncrement = " << zmumuIncrement_ << endl;
  double OneTrig = (double)SingleTrigger_ / (double)zmumuIncrement_;
  double DoubleTrig = (double)DoubleTrigger_ / (double)zmumuIncrement_;
  cout << "eps^2 = " << DoubleTrig << endl;
  cout << "2eps(1 - eps) = " << OneTrig << endl;
  int dimVec = vectorPt.size();
  for (int i = 0; i < dimVec; ++i)
    cout << "Pt = " << vectorPt[i] << " ==> Eta = " << vectorEta[i] << endl;
}

Member Data Documentation

int testAnalyzer::DoubleTrigger_

private

Definition at line 38 of file triggerEfficiency.cc.

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

TH1D* testAnalyzer::h_EtaDist_Pt80_

private

Definition at line 50 of file triggerEfficiency.cc.

Referenced by analyze(), and testAnalyzer().

TH1D * testAnalyzer::h_numberMuon_

private

Definition at line 48 of file triggerEfficiency.cc.

Referenced by analyze(), and testAnalyzer().

TH1D * testAnalyzer::h_numberMuon_ptStudy_

private

Definition at line 49 of file triggerEfficiency.cc.

Referenced by analyze(), and testAnalyzer().

TH1D* testAnalyzer::h_numberTrigMuon_

private

Definition at line 48 of file triggerEfficiency.cc.

Referenced by analyze(), and testAnalyzer().

TH1D* testAnalyzer::h_numberTrigMuon_ptStudy_

private

Definition at line 49 of file triggerEfficiency.cc.

Referenced by analyze(), and testAnalyzer().

TH1D* testAnalyzer::h_pt_distribution_

private

Definition at line 47 of file triggerEfficiency.cc.

Referenced by analyze(), and testAnalyzer().

double testAnalyzer::maxEta_

private

Definition at line 43 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

double testAnalyzer::maxPt_

private

Definition at line 45 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

double testAnalyzer::minEta_

private

Definition at line 43 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

double testAnalyzer::minPt_

private

Definition at line 45 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

int testAnalyzer::nbinsEta_

private

Definition at line 42 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

int testAnalyzer::nbinsEtaPt_

private

Definition at line 46 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

int testAnalyzer::nbinsPt_

private

Definition at line 44 of file triggerEfficiency.cc.

Referenced by testAnalyzer().

int testAnalyzer::NoTrigger_

private

Definition at line 38 of file triggerEfficiency.cc.

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

string testAnalyzer::pathName_

private

Definition at line 41 of file triggerEfficiency.cc.

Referenced by analyze().

EDGetTokenT<vector<pat::Muon> > testAnalyzer::selectMuonToken_

private

Definition at line 39 of file triggerEfficiency.cc.

Referenced by analyze().

int testAnalyzer::SingleTrigger_

private

Definition at line 38 of file triggerEfficiency.cc.

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

vector<double> testAnalyzer::vectorEta

private

Definition at line 51 of file triggerEfficiency.cc.

Referenced by analyze(), and endJob().

vector<double> testAnalyzer::vectorPt

private

Definition at line 51 of file triggerEfficiency.cc.

Referenced by analyze(), and endJob().

int testAnalyzer::zmumuIncrement_

private

Definition at line 38 of file triggerEfficiency.cc.

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

EDGetTokenT<CandidateView> testAnalyzer::zMuMuToken_

private

Definition at line 40 of file triggerEfficiency.cc.

Referenced by analyze().