MuTriggerAnalyzer Class Reference

Inheritance diagram for MuTriggerAnalyzer:

Public Member Functions
	MuTriggerAnalyzer (const edm::ParameterSet &pset)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
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
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup) override
virtual void	endJob () override
bool	IsMuMatchedToHLTMu (const reco::Muon &, std::vector< reco::Particle > &, double, double)

Private Attributes
TH1D *	deltaPtOverPt_
TH1D *	deltaR_
double	etaMuCut_
edm::Handle< trigger::TriggerEvent >	handleTriggerEvent_
std::string	hltPath_
TH1D *	hTrigMuonPtDenS_
TH1D *	hTrigMuonPtNumS_
std::string	L3FilterName_
double	maxDeltaR_
double	maxDPtRel_
edm::EDGetTokenT< View< Muon > >	muonToken_
const int	nbins_
const double	ptMax_
double	ptMuCut_
edm::InputTag	trigEvTag_
edm::EDGetTokenT< trigger::TriggerEvent >	trigEvToken_
edm::Handle< edm::TriggerResults >	triggerResults_
edm::TriggerNames const *	trigNames_
edm::EDGetTokenT< TriggerResults >	trigToken_

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 &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 49 of file MuTriggerAnalyzer.cc.

Constructor & Destructor Documentation

MuTriggerAnalyzer::MuTriggerAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 113 of file MuTriggerAnalyzer.cc.

References deltaPtOverPt_, deltaR_, hTrigMuonPtDenS_, hTrigMuonPtNumS_, TFileService::make(), maxDeltaR_, maxDPtRel_, nbins_, and ptMax_.

                                                                :
  trigToken_(consumes<TriggerResults>(cfg.getParameter<edm::InputTag> ("TrigTag"))),
  trigEvTag_(cfg.getParameter<edm::InputTag> ("triggerEvent")),
  trigEvToken_(consumes< trigger::TriggerEvent >(trigEvTag_)),
  muonToken_(consumes<View<Muon> > (cfg.getUntrackedParameter<edm::InputTag>("muons"))),
  ptMuCut_(cfg.getUntrackedParameter<double>("ptMuCut")),
  etaMuCut_(cfg.getUntrackedParameter<double>("etaMuCut")),
  hltPath_(cfg.getParameter<std::string >("hltPath")),
  L3FilterName_(cfg.getParameter<std::string >("L3FilterName")),
  maxDPtRel_(cfg.getParameter<double>("maxDPtRel")),
  maxDeltaR_(cfg.getParameter<double>("maxDeltaR")),
  nbins_(cfg.getParameter<double>("ptMax_")),
  ptMax_(cfg.getParameter<double>("ptMax_")){
   Service<TFileService> fs;

 hTrigMuonPtNumS_ = fs->make<TH1D>("hTrigMuonPtNumS", "hTrigMuonPtNumS", nbins_ + 1, 0, ptMax_);
  hTrigMuonPtDenS_ = fs->make<TH1D>("hTrigMuonPtDenS", "hTrigMuonPtDenS", nbins_ +1 , 0, ptMax_);
  deltaR_ = fs->make<TH1D>("deltaR", "deltaR", nbins_+1, 0, maxDeltaR_);
  deltaPtOverPt_ = fs->make<TH1D>("deltaPtOverPt", "deltaPtOverPt", nbins_ + 1, 0, maxDPtRel_);

}

Member Function Documentation

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

overrideprivatevirtual

Definition at line 144 of file MuTriggerAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), gather_cfg::cout, DEFINE_FWK_MODULE, edm::InputTag::encode(), stringResolutionProvider_cfi::eta, reco::LeafCandidate::eta(), etaMuCut_, trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), hltPath_, mps_fire::i, reco::Muon::isGlobalMuon(), IsMuMatchedToHLTMu(), gen::k, L3FilterName_, maxDeltaR_, maxDPtRel_, RPCpg::mu, electronCleaner_cfi::muons, muonToken_, dataset::name, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::p4(), EnergyCorrector::pt, reco::LeafCandidate::pt(), ptMuCut_, edm::HLTGlobalStatus::size(), trigger::TriggerEvent::sizeFilters(), AlCaHLTBitMon_QueryRunRegistry::string, trigEvTag_, trigEvToken_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), triggerResults, trigNames, trigNames_, and trigToken_.

                                                                     {
          //      bool overlap = false;

      // Trigger
      Handle<TriggerResults> triggerResults;
      TriggerNames trigNames;
      if (!ev.getByToken(trigToken_, triggerResults)) {
    LogWarning("") << ">>> TRIGGER collection does not exist !!!";
    return;
      }
      // trigNames.init(*triggerResults);
      trigNames_ = &ev.triggerNames(*triggerResults);
      //bool trigger_fired = false;
      for (unsigned int i=0; i<triggerResults->size(); i++) {


        std::string trigName = trigNames_->triggerName(i);
    if ( trigName == hltPath_ && triggerResults->accept(i)) {
      //trigger_fired = true;
    }
      }
      edm::Handle< trigger::TriggerEvent > handleTriggerEvent;
      //  LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << hltPath_ << ")";
      if ( ! ev.getByToken( trigEvToken_, handleTriggerEvent ))  {
    LogWarning( "errorTriggerEventValid" ) << "trigger::TriggerEvent product with InputTag " << trigEvTag_.encode() << " not in event";
    return;
      }
      const trigger::TriggerObjectCollection & toc(handleTriggerEvent->getObjects());
      size_t nMuHLT =0;
      std::vector<reco::Particle>  HLTMuMatched;
      for ( size_t ia = 0; ia < handleTriggerEvent->sizeFilters(); ++ ia) {
    std::string fullname = handleTriggerEvent->filterTag(ia).encode();
    std::string name;
    size_t p = fullname.find_first_of(':');
    if ( p != std::string::npos) {
      name = fullname.substr(0, p);
    }
    else {
      name = fullname;
    }
    if ( toc.size() != 0 ) {
      const trigger::Keys & k = handleTriggerEvent->filterKeys(ia);
      for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
        if (name == L3FilterName_  ) {
          HLTMuMatched.push_back(toc[*ki].particle());
          nMuHLT++;
        }
      }
    }
      }

    //  looping on muon....
      Handle<View<Muon> >   muons;
      if (!ev.getByToken(muonToken_, muons)) {
    LogError("") << ">>> muon collection does not exist !!!";
    return;
      }
      //saving only muons with pt> ptMuCut and eta<etaMuCut
      std::vector<reco::Muon>  highPtGlbMuons;

      for (unsigned int i=0; i<muons->size(); i++ ){
        const reco::Muon & mu = muons->at(i);
    double pt = mu.pt();
    double eta = mu.eta();
    if (pt> ptMuCut_ && fabs(eta)< etaMuCut_) {
      if (mu.isGlobalMuon()) highPtGlbMuons.push_back(mu);
    }
      }
      unsigned int nHighPtGlbMu = highPtGlbMuons.size();
      std::cout << "I've got " << nHighPtGlbMu << " nHighPtGlbMu" << std::endl;
      // unsigned int nHighPtStaMu = highPtStaMuons.size();

    // stop the loop after 10 cicles....
        (nHighPtGlbMu> 10)?   nHighPtGlbMu=10 : 1;

    if (nHighPtGlbMu>0 ){

           for(unsigned int i =0 ; i < nHighPtGlbMu ; i++) {
        reco::Muon muon1 = highPtGlbMuons[i];
        math::XYZTLorentzVector mu1(muon1.p4());
        //      double pt1= muon1.pt();

        /* bool singleTrigFlag1 =*/ IsMuMatchedToHLTMu ( muon1,  HLTMuMatched ,maxDeltaR_, maxDPtRel_ );

      }

    }

}

void MuTriggerAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 135 of file MuTriggerAnalyzer.cc.

References gather_cfg::cout, hTrigMuonPtDenS_, hTrigMuonPtNumS_, mps_fire::i, and nbins_.

                               {
  for (int i = 0; i < nbins_+1; ++i){
    std::cout << "number of reco muon in bin " << i << " = " << hTrigMuonPtDenS_->GetBinContent(i) << std::endl;
    std::cout << "number of hlt muon in bin " << i << " = " << hTrigMuonPtNumS_->GetBinContent(i) << std::endl;


  }
}

bool MuTriggerAnalyzer::IsMuMatchedToHLTMu ( const reco::Muon &  mu, std::vector< reco::Particle > &  HLTMu, double  DR, double  DPtRel  )

private

Definition at line 80 of file MuTriggerAnalyzer.cc.

References gather_cfg::cout, deltaR(), deltaR_, reco::LeafCandidate::eta(), HcalObjRepresent::Fill(), gen::k, nJettinessAdder_cfi::nPass, EnergyCorrector::pt, and reco::LeafCandidate::pt().

Referenced by analyze().

                                                                                                                            {
  size_t dim =  HLTMu.size();
  size_t nPass=0;

  // filling the denumerator;
  double muRecoPt= mu.pt();
  hTrigMuonPtDenS_-> Fill(muRecoPt);

  if (dim==0) return false;
  for (size_t k =0; k< dim; k++ ) {
    if (  (deltaR(HLTMu[k], mu) < DR)   && (fabs(HLTMu[k].pt() - mu.pt())/ HLTMu[k].pt()<DPtRel)){
      nPass++ ;
      std::cout<< " matching a muon " << std::endl;
      std::cout  << "muon reco  pt : " << muRecoPt<< std::endl;
      std::cout  << "muon reco  eta " <<  mu.eta() << std::endl;
      std::cout << "muon trigger  pt "<< HLTMu[k].pt() << std::endl;
      // filling the numerator, at the same bin as the denum.....
      hTrigMuonPtNumS_-> Fill(muRecoPt);
      deltaR_-> Fill(deltaR(HLTMu[k], mu));
      deltaPtOverPt_-> Fill(fabs(HLTMu[k].pt() - mu.pt())/ HLTMu[k].pt());


      std::cout << "muon trigger  eta : "<< HLTMu[k].eta() << std::endl;
      std::cout  <<"deltaR((HLTMu[k], mu)): "<< deltaR(HLTMu[k], mu) << std::endl;
      std::cout  <<"deltaPtOverPt: "<< fabs(HLTMu[k].pt() - mu.pt())/ HLTMu[k].pt() << std::endl;
    }
  }

  return (nPass>0);
}

Member Data Documentation

TH1D* MuTriggerAnalyzer::deltaPtOverPt_

private

Definition at line 75 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

TH1D* MuTriggerAnalyzer::deltaR_

private

Definition at line 74 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

double MuTriggerAnalyzer::etaMuCut_

private

Definition at line 62 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::Handle< trigger::TriggerEvent > MuTriggerAnalyzer::handleTriggerEvent_

private

Definition at line 67 of file MuTriggerAnalyzer.cc.

std::string MuTriggerAnalyzer::hltPath_

private

Definition at line 63 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

TH1D* MuTriggerAnalyzer::hTrigMuonPtDenS_

private

Definition at line 73 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

TH1D* MuTriggerAnalyzer::hTrigMuonPtNumS_

private

Definition at line 72 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

std::string MuTriggerAnalyzer::L3FilterName_

private

Definition at line 64 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

double MuTriggerAnalyzer::maxDeltaR_

private

Definition at line 69 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

double MuTriggerAnalyzer::maxDPtRel_

private

Definition at line 69 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

edm::EDGetTokenT<View<Muon> > MuTriggerAnalyzer::muonToken_

private

Definition at line 61 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

const int MuTriggerAnalyzer::nbins_

private

Definition at line 70 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

const double MuTriggerAnalyzer::ptMax_

private

Definition at line 71 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

double MuTriggerAnalyzer::ptMuCut_

private

Definition at line 62 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::InputTag MuTriggerAnalyzer::trigEvTag_

private

Definition at line 59 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT< trigger::TriggerEvent > MuTriggerAnalyzer::trigEvToken_

private

Definition at line 60 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::Handle<edm::TriggerResults> MuTriggerAnalyzer::triggerResults_

private

Definition at line 65 of file MuTriggerAnalyzer.cc.

edm::TriggerNames const* MuTriggerAnalyzer::trigNames_

private

Definition at line 66 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<TriggerResults> MuTriggerAnalyzer::trigToken_

private

Definition at line 58 of file MuTriggerAnalyzer.cc.

Referenced by analyze().