#include <DQMOffline/CalibCalo/src/DQMHcalIsoTrackAlCaReco.cc>

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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

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)

virtual	~EDConsumerBase () noexcept(false)

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 edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

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

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

Definition at line 46 of file DQMHcalIsoTrackAlCaReco.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 28 of file DQMHcalIsoTrackAlCaReco.cc.

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

                                                                                {
   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";
 }

DQMHcalIsoTrackAlCaReco::~DQMHcalIsoTrackAlCaReco ( )

override

Definition at line 52 of file DQMHcalIsoTrackAlCaReco.cc.

52 {}

Member Function Documentation

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

override

Definition at line 54 of file DQMHcalIsoTrackAlCaReco.cc.

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

                                                                                          {
 
   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";
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 187 of file DQMHcalIsoTrackAlCaReco.cc.

References DQMStore::IBooker::book1D(), etaRange_, folderName_, hEnEcal_, hHltEta_, hHltP_, hHltPhi_, hIeta_, hIphi_, hL1Eta_, hL1phi_, hL1Pt_, hL3Dr_, hL3Rat_, hltFilterTag_, hMaxP_, hOffP_, ifL3_, indexH_, gen::k, checklumidiff::l, l1FilterTag_, LogDebug, dataset::name, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, and type_.

                                                  {
 
   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);
 }