#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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 >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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_FULL_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);
 }