#include <DQMHcalIsoTrackAlCaReco.h>

Inheritance diagram for DQMHcalIsoTrackAlCaReco:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	DQMHcalIsoTrackAlCaReco (const edm::ParameterSet &)

	~DQMHcalIsoTrackAlCaReco () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
std::vector< double >	etaRange_

std::string	folderName_

MonitorElement *	hEnEcal_

std::vector< MonitorElement * >	hHltEta_

std::vector< MonitorElement * >	hHltP_

std::vector< MonitorElement * >	hHltPhi_

MonitorElement *	hIeta_

MonitorElement *	hIphi_

std::vector< MonitorElement * >	hL1Eta_

std::vector< MonitorElement * >	hL1phi_

std::vector< MonitorElement * >	hL1Pt_

MonitorElement *	hL3Dr_

MonitorElement *	hL3Rat_

std::vector< std::string >	hltFilterTag_

MonitorElement *	hMaxP_

std::vector< MonitorElement * >	hOffP_

std::vector< bool >	ifL3_

std::vector< unsigned int >	indexH_

std::vector< std::string >	l1FilterTag_

edm::InputTag	labelTrack_

edm::InputTag	labelTrigger_

int	nHLTaccepts_

int	nTotal_

double	pThr_

edm::EDGetTokenT< reco::HcalIsolatedTrackCandidateCollection >	tokTrack_

edm::EDGetTokenT< trigger::TriggerEvent >	tokTrigger_

std::vector< int >	type_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQMOffline/CalibCalo/src/DQMHcalIsoTrackAlCaReco.cc

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 45 of file DQMHcalIsoTrackAlCaReco.h.

Constructor & Destructor Documentation

◆ DQMHcalIsoTrackAlCaReco()

DQMHcalIsoTrackAlCaReco::DQMHcalIsoTrackAlCaReco ( const edm::ParameterSet & iConfig )

explicit

Definition at line 29 of file DQMHcalIsoTrackAlCaReco.cc.

                                                                                {
   folderName_ = iConfig.getParameter<std::string>("FolderName");
   l1FilterTag_ = iConfig.getParameter<std::vector<std::string>>("L1FilterLabel");
   hltFilterTag_ = iConfig.getParameter<std::vector<std::string>>("HltFilterLabels");
   type_ = iConfig.getParameter<std::vector<int>>("TypeFilter");
   labelTrigger_ = iConfig.getParameter<edm::InputTag>("TriggerLabel");
   labelTrack_ = iConfig.getParameter<edm::InputTag>("TracksLabel");
   pThr_ = iConfig.getUntrackedParameter<double>("pThrL3", 0);
  
   nTotal_ = nHLTaccepts_ = 0;
   tokTrigger_ = consumes<trigger::TriggerEvent>(labelTrigger_);
   tokTrack_ = consumes<reco::HcalIsolatedTrackCandidateCollection>(labelTrack_);
   LogDebug("HcalIsoTrack") << "Folder " << folderName_ << " Input Tag for Trigger " << labelTrigger_ << " track "
                            << labelTrack_ << " threshold " << pThr_ << " with " << l1FilterTag_.size()
                            << " level 1 and " << hltFilterTag_.size() << " hlt filter tags"
                            << "\n";
   for (unsigned int k = 0; k < l1FilterTag_.size(); ++k)
     LogDebug("HcalIsoTrack") << "L1FilterTag[" << k << "] " << l1FilterTag_[k] << "\n";
   for (unsigned int k = 0; k < hltFilterTag_.size(); ++k)
     LogDebug("HcalIsoTrack") << "HLTFilterTag[" << k << "] " << hltFilterTag_[k] << "\n";
 }

References folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltFilterTag_, dqmdumpme::k, l1FilterTag_, labelTrack_, labelTrigger_, LogDebug, nHLTaccepts_, nTotal_, pThr_, AlCaHLTBitMon_QueryRunRegistry::string, tokTrack_, tokTrigger_, and type_.

◆ ~DQMHcalIsoTrackAlCaReco()

DQMHcalIsoTrackAlCaReco::~DQMHcalIsoTrackAlCaReco ( )

override

Definition at line 51 of file DQMHcalIsoTrackAlCaReco.cc.

51 {}

Member Function Documentation

◆ analyze()

void DQMHcalIsoTrackAlCaReco::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 53 of file DQMHcalIsoTrackAlCaReco.cc.

                                                                                          {
   nTotal_++;
   bool accept(false);
   edm::Handle<trigger::TriggerEvent> trEv;
   iEvent.getByToken(tokTrigger_, trEv);
  
   edm::Handle<reco::HcalIsolatedTrackCandidateCollection> recoIsoTracks;
   iEvent.getByToken(tokTrack_, recoIsoTracks);
   LogDebug("HcalIsoTrack") << "Gets Trigger information with " << trEv.isValid() << " and offline tracks with "
                            << recoIsoTracks.isValid() << "\n";
  
   if (trEv.isValid()) {
     const trigger::TriggerObjectCollection &TOCol(trEv->getObjects());
     const trigger::size_type nFilt(trEv->sizeFilters());
     // plots for L1 trigger
     for (unsigned int k = 0; k < l1FilterTag_.size(); ++k) {
       trigger::Keys KEYSl1;
       double etaTrigl1(-10000), phiTrigl1(-10000), ptMaxl1(0);
       for (trigger::size_type iFilt = 0; iFilt != nFilt; iFilt++) {
         LogDebug("HcalIsoTrack") << trEv->filterTag(iFilt).label() << " find for " << l1FilterTag_[k] << " gives "
                                  << (trEv->filterTag(iFilt).label()).find(l1FilterTag_[k]) << "\n";
         if ((trEv->filterTag(iFilt).label()).find(l1FilterTag_[k]) != std::string::npos) {
           KEYSl1 = trEv->filterKeys(iFilt);
           trigger::size_type nRegl1 = KEYSl1.size();
           LogDebug("HcalIsoTrack") << "# of objects " << nRegl1 << "\n";
           for (trigger::size_type iReg = 0; iReg < nRegl1; iReg++) {
             const trigger::TriggerObject &TObj(TOCol[KEYSl1[iReg]]);
             LogDebug("HcalIsoTrack") << "Object[" << iReg << "] with pt " << TObj.pt() << " " << TObj.eta() << " "
                                      << TObj.phi() << "\n";
             if (TObj.pt() > ptMaxl1) {
               etaTrigl1 = TObj.eta();
               phiTrigl1 = TObj.phi();
               ptMaxl1 = TObj.pt();
             }
           }
         }
       }
       LogDebug("HcalIsoTrack") << "For L1 trigger type " << k << " pt " << ptMaxl1 << " eta " << etaTrigl1 << " phi "
                                << phiTrigl1 << "\n";
       if (ptMaxl1 > 0) {
         hL1Pt_[k]->Fill(ptMaxl1);
         hL1Eta_[k]->Fill(etaTrigl1);
         hL1phi_[k]->Fill(phiTrigl1);
       }
     }
     // Now make plots for hlt objects
     trigger::Keys KEYS;
     for (unsigned l = 0; l < hltFilterTag_.size(); l++) {
       for (trigger::size_type iFilt = 0; iFilt != nFilt; iFilt++) {
         LogDebug("HcalIsoTrack") << trEv->filterTag(iFilt).label() << " find for " << hltFilterTag_[l] << " gives "
                                  << (trEv->filterTag(iFilt).label()).find(hltFilterTag_[l]) << "\n";
         if ((trEv->filterTag(iFilt).label()).find(hltFilterTag_[l]) != std::string::npos) {
           KEYS = trEv->filterKeys(iFilt);
           trigger::size_type nReg = KEYS.size();
           LogDebug("HcalIsoTrack") << "# of objects for HLT " << nReg << "\n";
           // checks with IsoTrack trigger results
           for (trigger::size_type iReg = 0; iReg < nReg; iReg++) {
             const trigger::TriggerObject &TObj(TOCol[KEYS[iReg]]);
             LogDebug("HcalIsoTrack") << "HLT Filter Tag " << l << " trigger " << iFilt << " object " << iReg << " p "
                                      << TObj.p() << " pointer " << indexH_[l] << ":" << hHltP_[indexH_[l]] << ":"
                                      << hHltEta_[indexH_[l]] << ":" << hHltPhi_[indexH_[l]] << "\n";
             if (TObj.p() > pThr_) {
               hHltP_[indexH_[l]]->Fill(TObj.p());
               hHltEta_[indexH_[l]]->Fill(TObj.eta());
               hHltPhi_[indexH_[l]]->Fill(TObj.phi());
               if (ifL3_[l])
                 accept = true;
               if (recoIsoTracks.isValid() && ifL3_[l]) {
                 double minRecoL3dist(1000), pt(1000);
                 reco::HcalIsolatedTrackCandidateCollection::const_iterator mrtr;
                 for (mrtr = recoIsoTracks->begin(); mrtr != recoIsoTracks->end(); mrtr++) {
                   double R = deltaR(mrtr->eta(), mrtr->phi(), TObj.eta(), TObj.phi());
                   if (R < minRecoL3dist) {
                     minRecoL3dist = R;
                     pt = mrtr->pt();
                   }
                 }
                 LogDebug("HcalIsoTrack") << "Minimum R " << minRecoL3dist << " pt " << pt << ":" << TObj.pt() << "\n";
                 hL3Dr_->Fill(minRecoL3dist);
                 if (minRecoL3dist < 0.02)
                   hL3Rat_->Fill(TObj.pt() / pt);
               }
             }
           }
         }
       }
     }
   }
  
   // general distributions
   if (recoIsoTracks.isValid()) {
     for (reco::HcalIsolatedTrackCandidateCollection::const_iterator itr = recoIsoTracks->begin();
          itr != recoIsoTracks->end();
          itr++) {
       hMaxP_->Fill(itr->maxP());
       hEnEcal_->Fill(itr->energyEcal());
       std::pair<int, int> etaphi = itr->towerIndex();
       hIeta_->Fill(etaphi.first);
       hIphi_->Fill(etaphi.second);
       LogDebug("HcalIsoTrack") << "Reco track p " << itr->p() << " eta|phi " << etaphi.first << "|" << etaphi.second
                                << " maxP " << itr->maxP() << " EcalE " << itr->energyEcal() << " pointers " << hHltP_[3]
                                << ":" << hHltEta_[3] << ":" << hHltPhi_[3] << "\n";
       if (itr->p() >= pThr_) {
         hHltP_[3]->Fill(itr->p());
         hHltEta_[3]->Fill(itr->eta());
         hHltPhi_[3]->Fill(itr->phi());
       }
       double etaAbs = std::abs(itr->eta());
       hOffP_[0]->Fill(itr->p());
       for (unsigned int l = 1; l < etaRange_.size(); l++) {
         if (etaAbs >= etaRange_[l - 1] && etaAbs < etaRange_[l]) {
           LogDebug("HcalIsoTrack") << "Range " << l << " p " << itr->p() << " pointer " << hOffP_[l];
           hOffP_[l]->Fill(itr->p());
           break;
         }
       }
     }
   }
  
   if (accept)
     nHLTaccepts_++;
   LogDebug("HcalIsoTrack") << "Accept " << accept << "\n";
 }

References funct::abs(), accept(), PbPb_ZMuSkimMuonDPG_cff::deltaR, trigger::TriggerObject::eta(), fastmath::etaphi(), etaRange_, dqm::impl::MonitorElement::Fill(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), spr::find(), trigger::TriggerEvent::getObjects(), hEnEcal_, hHltEta_, hHltP_, hHltPhi_, hIeta_, hIphi_, hL1Eta_, hL1phi_, hL1Pt_, hL3Dr_, hL3Rat_, hltFilterTag_, hMaxP_, hOffP_, iEvent, ifL3_, indexH_, edm::HandleBase::isValid(), dqmdumpme::k, cmsLHEtoEOSManager::l, l1FilterTag_, edm::InputTag::label(), LogDebug, HLT_2018_cff::nFilt, nHLTaccepts_, nTotal_, trigger::TriggerObject::p(), trigger::TriggerObject::phi(), DiDispStaMuonMonitor_cfi::pt, trigger::TriggerObject::pt(), pThr_, dttmaxenums::R, trigger::TriggerEvent::sizeFilters(), tokTrack_, and tokTrigger_.

◆ bookHistograms()

void DQMHcalIsoTrackAlCaReco::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 177 of file DQMHcalIsoTrackAlCaReco.cc.

                                                                                                             {
   iBooker.setCurrentFolder(folderName_);
   LogDebug("HcalIsoTrack") << "Set the folder to " << folderName_ << "\n";
   char name[100], title[200];
   for (unsigned int k = 0; k < l1FilterTag_.size(); ++k) {
     sprintf(name, "hp%s", l1FilterTag_[k].c_str());
     sprintf(title, "p_T of L1 object for %s", l1FilterTag_[k].c_str());
     hL1Pt_.push_back(iBooker.book1D(name, title, 1000, 0, 1000));
     hL1Pt_[k]->setAxisTitle("p_T (GeV)", 1);
     sprintf(name, "heta%s", l1FilterTag_[k].c_str());
     sprintf(title, "#eta of L1 object for %s", l1FilterTag_[k].c_str());
     hL1Eta_.push_back(iBooker.book1D(name, title, 100, -2.5, 2.5));
     hL1Eta_[k]->setAxisTitle("#eta", 1);
     sprintf(name, "hphi%s", l1FilterTag_[k].c_str());
     sprintf(title, "#phi of L1 object for %s", l1FilterTag_[k].c_str());
     hL1phi_.push_back(iBooker.book1D(name, title, 100, -3.2, 3.2));
     hL1phi_[k]->setAxisTitle("#phi", 1);
   }
  
   std::string types[4] = {"L2", "L2x", "L3", "Off"};
   for (unsigned int l = 0; l < 4; l++) {
     sprintf(name, "hp%s", types[l].c_str());
     sprintf(title, "Momentum of %s object", types[l].c_str());
     hHltP_.push_back(iBooker.book1D(name, title, 200, 0, 1000));
     hHltP_[l]->setAxisTitle("p (GeV)", 1);
     sprintf(name, "heta%s", types[l].c_str());
     sprintf(title, "#eta of %s object", types[l].c_str());
     hHltEta_.push_back(iBooker.book1D(name, title, 16, -2, 2));
     hHltEta_[l]->setAxisTitle("#eta", 1);
     sprintf(name, "hphi%s", types[l].c_str());
     sprintf(title, "#phi of %s object", types[l].c_str());
     hHltPhi_.push_back(iBooker.book1D(name, title, 16, -3.2, 3.2));
     hHltPhi_[l]->setAxisTitle("#phi", 1);
   }
   sprintf(title, "Distance of offline track from L3 object");
   hL3Dr_ = (iBooker.book1D("hDRL3", title, 40, 0, 0.2));
   hL3Dr_->setAxisTitle("R(#eta,#phi)", 1);
   sprintf(title, "Ratio of p L3/Offline");
   hL3Rat_ = (iBooker.book1D("hRatL3", title, 500, 0, 3));
   indexH_.clear();
   ifL3_.clear();
   for (unsigned int l = 0; l < hltFilterTag_.size(); l++) {
     unsigned int indx = (type_[l] >= 0 && type_[l] < 3) ? type_[l] : 0;
     indexH_.push_back(indx);
     ifL3_.push_back(indx == 2);
     LogDebug("HcalIsoTrack") << "Filter[" << l << "] " << hltFilterTag_[l] << " type " << type_[l] << " index "
                              << indexH_[l] << " L3? " << ifL3_[l] << "\n";
   }
  
   double etaV[6] = {0.0, 0.5, 1.0, 1.5, 2.0, 2.5};
   for (unsigned int k = 0; k < 6; ++k) {
     sprintf(name, "hOffP%d", k);
     if (k == 0) {
       sprintf(title, "p of AlCaReco object (All)");
     } else {
       sprintf(title, "p of AlCaReco object (%3.1f < |#eta| < %3.1f)", etaV[k - 1], etaV[k]);
     }
     etaRange_.push_back(etaV[k]);
     hOffP_.push_back(iBooker.book1D(name, title, 1000, 0, 1000));
     hOffP_[k]->setAxisTitle("E (GeV)", 1);
   }
   hMaxP_ = iBooker.book1D("hChgIsol", "Energy for charge isolation", 110, -10, 100);
   hMaxP_->setAxisTitle("p (GeV)", 1);
   hEnEcal_ = iBooker.book1D("hEnEcal", "Energy in ECAL", 100, 0, 20);
   hEnEcal_->setAxisTitle("E (GeV)", 1);
   hIeta_ = iBooker.book1D("hIEta", "i#eta for HCAL tower", 90, -45, 45);
   hIeta_->setAxisTitle("i#eta", 1);
   hIphi_ = iBooker.book1D("hIPhi", "i#phi for HCAL tower", 72, 0, 72);
   hIphi_->setAxisTitle("i#phi", 1);
 }