#include <Calibration/IsolatedParticles/plugins/StudyHLT.cc>

Inheritance diagram for StudyHLT:

Public Member Functions
	StudyHLT (const edm::ParameterSet &)

	~StudyHLT () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

void	beginJob () override

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	clear ()

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

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

void	fillEnergy (int, int, double, double, double, double, double)

void	fillIsolation (int, double, double, double)

void	fillTrack (int, double, double, double, double)

int	trackPID (const reco::Track *, const edm::Handle< reco::GenParticleCollection > &)

std::string	truncate_str (const std::string &)

Private Attributes
bool	changed_

const bool	doTree_

int	etaBin_ [nEtaBin_+1]

bool	firstEvent_

edm::Service< TFileService >	fs_

TH1D *	h_ediff [2]

TH1D *	h_ene1 [2]

TH1D *	h_ene2 [2]

TH1D *	h_energy [nPVBin_+8][nPBin_][nEtaBin_][6]

TH1D *	h_eta [nGen_+2]

TH1I *	h_goodPV

TH1I *	h_goodRun

TH1I *	h_HLTAccept

std::vector< TH1I * >	h_HLTAccepts

TH1I *	h_HLTCorr

TH1D *	h_maxNearP [2]

TH1I *	h_nHLT

TH2I *	h_nHLTvsRN

TH1I *	h_ntrk [2]

TH1I *	h_numberPV

TH1D *	h_p [nGen_+2]

TH1D *	h_phi [nGen_+2]

TH1D *	h_pt [nGen_+2]

HLTConfigProvider	hltConfig_

std::vector< std::string >	HLTNames_

const bool	isItAOD_

const double	maxTrackEta_

const double	minTrackP_

const std::vector< std::string >	newNames_

int	nRun_

double	pBin_ [nPBin_+1]

const std::vector< double >	puWeights_

int	pvBin_ [nPVBin_+1]

spr::trackSelectionParameters	selectionParameters_

const std::string	theTrackQuality_

const edm::InputTag	theTriggerResultsLabel_

const double	tMaxE_

const double	tMaxH_

const double	tMinE_

const double	tMinH_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< LumiDetails >	tok_lumi

edm::EDGetTokenT< reco::GenParticleCollection >	tok_parts_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes

double	tr_eventWeight

std::vector< double >	tr_FE11x11P

std::vector< double >	tr_FE15x15P

std::vector< double >	tr_FE7x7P

int	tr_goodPV

int	tr_goodRun

std::vector< double >	tr_H3x3

std::vector< double >	tr_H5x5

std::vector< double >	tr_H7x7

std::vector< int >	tr_iEta

std::vector< double >	tr_MaxNearHcalP7x7

std::vector< double >	tr_MaxNearP31X31

std::vector< bool >	tr_SE11x11P

std::vector< bool >	tr_SE15x15P

std::vector< bool >	tr_SE7x7P

std::vector< std::string >	tr_TrigName

std::vector< double >	tr_TrkEta

std::vector< int >	tr_TrkID

std::vector< double >	tr_TrkP

std::vector< double >	tr_TrkPhi

std::vector< double >	tr_TrkPt

TTree *	tree_

const edm::InputTag	triggerEvent_

const std::vector< std::string >	trigNames_

const int	verbosity_

Static Private Attributes
static const int	nEtaBin_ =4

static const int	nGen_ =(nPVBin_+12)

static const int	nPBin_ =10

static const int	nPVBin_ =4

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Description: Studies single particle response measurements in data/MC

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

Definition at line 82 of file StudyHLT.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 157 of file StudyHLT.cc.

                                                  : 
   verbosity_(iConfig.getUntrackedParameter<int>("verbosity",0)),
   trigNames_(iConfig.getUntrackedParameter<std::vector<std::string> >("triggers")),
   newNames_(iConfig.getUntrackedParameter<std::vector<std::string> >("newNames")),
   theTrackQuality_(iConfig.getUntrackedParameter<std::string>("trackQuality","highPurity")),
   minTrackP_(iConfig.getUntrackedParameter<double>("minTrackP",1.0)),
   maxTrackEta_(iConfig.getUntrackedParameter<double>("maxTrackEta",2.5)),
   tMinE_(iConfig.getUntrackedParameter<double>("timeMinCutECAL",-500.)),
   tMaxE_(iConfig.getUntrackedParameter<double>("timeMaxCutECAL",500.)),
   tMinH_(iConfig.getUntrackedParameter<double>("timeMinCutHCAL",-500.)),
   tMaxH_(iConfig.getUntrackedParameter<double>("timeMaxCutHCAL",500.)),
   isItAOD_(iConfig.getUntrackedParameter<bool>("isItAOD",false)),
   doTree_(iConfig.getUntrackedParameter<bool>("doTree",false)),
   puWeights_(iConfig.getUntrackedParameter<std::vector<double> >("puWeights")),
   triggerEvent_(edm::InputTag("hltTriggerSummaryAOD","","HLT")),
   theTriggerResultsLabel_(edm::InputTag("TriggerResults","","HLT")), nRun_(0) {
 
   usesResource(TFileService::kSharedResource);
 
   reco::TrackBase::TrackQuality trackQuality=reco::TrackBase::qualityByName(theTrackQuality_);
   selectionParameters_.minPt          = iConfig.getUntrackedParameter<double>("minTrackPt", 10.0);
   selectionParameters_.minQuality     = trackQuality;
   selectionParameters_.maxDxyPV       = iConfig.getUntrackedParameter<double>("maxDxyPV", 0.2);
   selectionParameters_.maxDzPV        = iConfig.getUntrackedParameter<double>("maxDzPV",  5.0);
   selectionParameters_.maxChi2        = iConfig.getUntrackedParameter<double>("maxChi2",  5.0);
   selectionParameters_.maxDpOverP     = iConfig.getUntrackedParameter<double>("maxDpOverP",  0.1);
   selectionParameters_.minOuterHit    = iConfig.getUntrackedParameter<int>("minOuterHit", 4);
   selectionParameters_.minLayerCrossed= iConfig.getUntrackedParameter<int>("minLayerCrossed", 8);
   selectionParameters_.maxInMiss      = iConfig.getUntrackedParameter<int>("maxInMiss", 0);
   selectionParameters_.maxOutMiss     = iConfig.getUntrackedParameter<int>("maxOutMiss", 0);
 
   // define tokens for access
   tok_lumi      = consumes<LumiDetails, edm::InLumi>(edm::InputTag("lumiProducer"));
   tok_trigEvt   = consumes<trigger::TriggerEvent>(triggerEvent_);
   tok_trigRes   = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
   tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
   tok_recVtx_   = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
   tok_parts_    = consumes<reco::GenParticleCollection>(iConfig.getParameter<edm::InputTag>("particleSource"));
 
   if (isItAOD_) {
     tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
     tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));
     tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("reducedHcalRecHits", "hbhereco"));
   } else {
     tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
     tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
     tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
   }
   tok_ew_       = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
   
   edm::LogInfo("IsoTrack") << "Verbosity " << verbosity_ << " with " 
                << trigNames_.size() << " triggers:";
   for (unsigned int k=0; k<trigNames_.size(); ++k)
     edm::LogInfo("IsoTrack") << " [" << k << "] " << trigNames_[k];
   edm::LogInfo("IsoTrack") << "TrackQuality " << theTrackQuality_ << " Minpt "
                << selectionParameters_.minPt << " maxDxy " 
                << selectionParameters_.maxDxyPV << " maxDz "
                << selectionParameters_.maxDzPV << " maxChi2 "
                << selectionParameters_.maxChi2 << " maxDp/p "
                << selectionParameters_.maxDpOverP << " minOuterHit "
                << selectionParameters_.minOuterHit << " minLayerCrossed "
                << selectionParameters_.minLayerCrossed << " maxInMiss "
                << selectionParameters_.maxInMiss << " maxOutMiss " 
                << selectionParameters_.maxOutMiss << " minTrackP "
                << minTrackP_ << " maxTrackEta " << maxTrackEta_
                << " tMinE_ " << tMinE_ << " tMaxE " << tMaxE_ 
                << " tMinH_ " << tMinH_ << " tMaxH_ " << tMaxH_ 
                << " isItAOD " << isItAOD_ << " doTree " << doTree_;
 
   double pBins[nPBin_+1] = {1.0,2.0,3.0,4.0,5.0,6.0,7.0,9.0,11.0,15.0,20.0};
   int    etaBins[nEtaBin_+1] = {1, 7, 13, 17, 23};
   int    pvBins[nPVBin_+1] = {1, 2, 3, 5, 100};
   for (int i=0; i<=nPBin_; ++i)   pBin_[i]   = pBins[i];
   for (int i=0; i<=nEtaBin_; ++i) etaBin_[i] = etaBins[i];
   for (int i=0; i<=nPVBin_; ++i)  pvBin_[i]  = pvBins[i];
 
   firstEvent_ = true;
   changed_    = false;
 }

StudyHLT::~StudyHLT ( )

inlineoverride

Definition at line 86 of file StudyHLT.cc.

References analyze(), beginJob(), beginRun(), endRun(), and fillDescriptions().

86 { }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 270 of file StudyHLT.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationEcal(), spr::chargeIsolationHcal(), clear(), doTree_, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, spr::eECALmatrix(), spr::eHCALmatrix(), edm::EventID::event(), fillEnergy(), fillIsolation(), fillTrack(), firstEvent_, particlelevel_cff::genEventInfo, GenHFHadronMatcher_cfi::genParticles, edm::EventSetup::get(), edm::Event::getByToken(), h_goodPV, h_goodRun, h_HLTAccept, h_HLTAccepts, h_HLTCorr, h_nHLT, h_nHLTvsRN, h_ntrk, h_numberPV, photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, HLTNames_, mps_fire::i, edm::EventBase::id(), createfilelist::int, edm::HandleBase::isValid(), edm::InputTag::label(), mix_E10TeV_L13E31_BX432_cfi::Lumi, edm::EventBase::luminosityBlock(), M_PI, match(), maxTrackEta_, spr::trackSelectionParameters::minQuality, minTrackP_, reco::TrackBase::momentum(), newNames_, nGen_, nPVBin_, nRun_, convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), p1, edm::Handle< T >::product(), edm::ESHandle< T >::product(), spr::propagateCALO(), reco::TrackBase::pt(), hiDetachedQuadStep_cff::pt1, puWeights_, pvBin_, jets_cff::quality, reco::TrackBase::quality(), edm::EventID::run(), selectionParameters_, edm::HLTGlobalStatus::size(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality_, tMaxE_, tMaxH_, tMinE_, tMinH_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_parts_, tok_recVtx_, tok_trigEvt, tok_trigRes, tr_eventWeight, tr_FE11x11P, tr_FE15x15P, tr_FE7x7P, tr_goodPV, tr_goodRun, tr_H3x3, tr_H5x5, tr_H7x7, tr_iEta, tr_MaxNearHcalP7x7, tr_MaxNearP31X31, tr_SE11x11P, tr_SE15x15P, tr_SE7x7P, tr_TrigName, tr_TrkEta, tr_TrkID, tr_TrkP, tr_TrkPhi, tr_TrkPt, trackPID(), tree_, triggerEvent_, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, trigNames_, truncate_str(), verbosity_, and GenEventInfoProduct::weight().

Referenced by ~StudyHLT().

                                                                           {
   clear();
   if (verbosity_ > 0) 
     edm::LogInfo("IsoTrack") << "Event starts===================================="; 
   int RunNo = iEvent.id().run();
   int EvtNo = iEvent.id().event();
   int Lumi  = iEvent.luminosityBlock();
   int Bunch = iEvent.bunchCrossing();
   
   std::vector<int>  newAccept(newNames_.size()+1,0);
   float mybxlumi=-1;
   /*
   edm::Handle<LumiDetails> Lumid;
   iEvent.getLuminosityBlock().getByToken(tok_lumi,Lumid);
   if (Lumid.isValid()) mybxlumi=Lumid->lumiValue(LumiDetails::kOCC1,iEvent.bunchCrossing())*6.37;
   */
   if (verbosity_ > 0)
     edm::LogInfo("IsoTrack") << "RunNo " << RunNo << " EvtNo " << EvtNo 
                  << " Lumi " << Lumi << " Bunch " << Bunch 
                  << " mybxlumi " << mybxlumi;
   
   trigger::TriggerEvent triggerEvent;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   iEvent.getByToken(tok_trigEvt,triggerEventHandle);
   
   bool ok(false);
   if (!triggerEventHandle.isValid()) {
     edm::LogWarning("IsoTrack") << "Error! Can't get the product "
                 << triggerEvent_.label();
   } else {
     triggerEvent = *(triggerEventHandle.product());
     
     edm::Handle<edm::TriggerResults> triggerResults;
     iEvent.getByToken(tok_trigRes, triggerResults);
 
     if (triggerResults.isValid()) {
       h_nHLT->Fill(triggerResults->size());
       h_nHLTvsRN->Fill(RunNo, triggerResults->size());
 
       const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);      
       const std::vector<std::string> & triggerNames_ = triggerNames.triggerNames();
       for (unsigned int iHLT=0; iHLT<triggerResults->size(); iHLT++) {
     //        unsigned int triggerindx = hltConfig_.triggerIndex(triggerNames_[iHLT]);
     //        const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
         int ipos=-1;
     std::string newtriggerName = truncate_str(triggerNames_[iHLT]);
     for (unsigned int i=0; i<HLTNames_.size(); ++i) {
       if (newtriggerName == HLTNames_[i]) {
         ipos = i+1;
         break;
       }
     }
     if (ipos < 0) {
       HLTNames_.push_back(newtriggerName);
       ipos = (int)(HLTNames_.size());
       if (ipos <= h_HLTAccept->GetNbinsX())
         h_HLTAccept->GetXaxis()->SetBinLabel(ipos,newtriggerName.c_str());
     }
     if ((int)(iHLT+1) > h_HLTAccepts[nRun_]->GetNbinsX()) {
       edm::LogInfo("IsoTrack") << "Wrong trigger " << RunNo << " Event " 
                    << EvtNo << " Hlt " << iHLT;
     } else {
       if (firstEvent_)  h_HLTAccepts[nRun_]->GetXaxis()->SetBinLabel(iHLT+1, newtriggerName.c_str());
     }
     int hlt    = triggerResults->accept(iHLT);
     if (hlt) {
       h_HLTAccepts[nRun_]->Fill(iHLT+1);
       h_HLTAccept->Fill(ipos);
     }
     if (trigNames_.empty()) {
       ok = true;
     } else {
       for (unsigned int i=0; i<trigNames_.size(); ++i) {
         if (newtriggerName.find(trigNames_[i]) != std::string::npos) {
           if (verbosity_%10 > 0)  
         edm::LogInfo("IsoTrack") << newtriggerName;
           if (hlt > 0) {
         ok = true;
         tr_TrigName.push_back(newtriggerName);
           }
         }
       }
       for (unsigned int i=0; i<newNames_.size(); ++i) {
         if (newtriggerName.find(newNames_[i]) != std::string::npos) {
           if (verbosity_%10 > 0)
         edm::LogInfo("IsoTrack") << "[" << i << "] " << newNames_[i] 
                      << " : " << newtriggerName;
           if (hlt > 0) newAccept[i] = 1;
         }
       }
     }
       }
       int iflg(0), indx(1);
       for (unsigned int i=0; i<newNames_.size(); ++i) {
     iflg += (indx*newAccept[i]); indx *= 2;
       }
       h_HLTCorr->Fill(iflg);
     }
   }
 
   //Look at the tracks  
   if (ok) {
     h_goodRun->Fill(RunNo);
     tr_goodRun = RunNo;
     // get handles to calogeometry and calotopology
     edm::ESHandle<CaloGeometry> pG;
     iSetup.get<CaloGeometryRecord>().get(pG);
     const CaloGeometry* geo = pG.product();
   
     edm::ESHandle<CaloTopology> theCaloTopology;
     iSetup.get<CaloTopologyRecord>().get(theCaloTopology); 
     const CaloTopology *caloTopology = theCaloTopology.product();
   
     edm::ESHandle<HcalTopology> htopo;
     iSetup.get<HcalRecNumberingRecord>().get(htopo);
     const HcalTopology* theHBHETopology = htopo.product();
  
     edm::ESHandle<MagneticField> bFieldH;
     iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
     const MagneticField *bField = bFieldH.product();
 
     edm::ESHandle<EcalChannelStatus> ecalChStatus;
     iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
     const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
 
     edm::Handle<reco::VertexCollection> recVtxs;
     iEvent.getByToken(tok_recVtx_,recVtxs);
     int       ntrk(0), ngoodPV(0), nPV(-1);
     int       nvtxs = (int)(recVtxs->size());
     for (int ind=0; ind<nvtxs; ind++) {
       if (!((*recVtxs)[ind].isFake()) && (*recVtxs)[ind].ndof() > 4) ngoodPV++;
     }
     for (int i=0; i<nPVBin_; ++i) {
       if (ngoodPV >= pvBin_[i] && ngoodPV < pvBin_[i+1]) {
     nPV = i; break;
       }
     }
 
     tr_eventWeight = 1.0;
     edm::Handle<GenEventInfoProduct> genEventInfo;
     iEvent.getByToken(tok_ew_, genEventInfo);
     if (genEventInfo.isValid()) tr_eventWeight = genEventInfo->weight();  
 
     if ((verbosity_/10)%10 > 0) 
       edm::LogInfo("IsoTrack") << "Number of vertices: " << nvtxs
                    << " Good " << ngoodPV << " Bin " << nPV
                    << " Event weight " << tr_eventWeight;
     h_numberPV->Fill(nvtxs,tr_eventWeight);
     h_goodPV->Fill(ngoodPV,tr_eventWeight);
     tr_goodPV   = ngoodPV;
     
     if (!puWeights_.empty()) {
       int npbin = h_goodPV->FindBin(ngoodPV);
       if (npbin > 0 && npbin <= (int)(puWeights_.size())) 
     tr_eventWeight *= puWeights_[npbin-1];
       else  
     tr_eventWeight  = 0;
     }
 
     edm::Handle<reco::TrackCollection> trkCollection;
     iEvent.getByToken(tok_genTrack_, trkCollection);
 
     //=== genParticle information
     edm::Handle<reco::GenParticleCollection> genParticles;
     iEvent.getByToken(tok_parts_, genParticles);
     if (genParticles.isValid()) {
       for (const auto & p : (reco::GenParticleCollection)(*genParticles)) {
     double pt1         = p.momentum().Rho();
     double p1          = p.momentum().R();
     double eta1        = p.momentum().Eta();
     double phi1        = p.momentum().Phi();
     fillTrack(nGen_, pt1,p1,eta1,phi1);
     bool match(false);
     double phi2(phi1);
     if (phi2 < 0) phi2 += 2.0*M_PI;
     for (const auto& trk : (reco::TrackCollection)(*trkCollection)) {
       bool quality = trk.quality(selectionParameters_.minQuality);
       if (quality) {
         double dEta = trk.eta() - eta1;
         double phi0 = trk.phi();
         if (phi0 < 0) phi0 += 2.0*M_PI;
         double dPhi = phi0-phi2;
         if (dPhi > M_PI)       dPhi -= 2.*M_PI;
         else if (dPhi < -M_PI) dPhi += 2.*M_PI;
         double dR   = sqrt(dEta*dEta+dPhi*dPhi);
         if (dR < 0.01) { match = true; break;}
       }
     }
     if (match) fillTrack(nGen_+1, pt1,p1,eta1,phi1);
       }
     }
 
     reco::TrackCollection::const_iterator trkItr;
     for (trkItr=trkCollection->begin(); trkItr != trkCollection->end(); ++trkItr,++ntrk) {
       const reco::Track* pTrack = &(*trkItr);
       double pt1         = pTrack->pt();
       double p1          = pTrack->p();
       double eta1        = pTrack->momentum().eta();
       double phi1        = pTrack->momentum().phi();
       bool quality       = pTrack->quality(selectionParameters_.minQuality);
       fillTrack(0, pt1,p1,eta1,phi1);
       if (quality) fillTrack(1, pt1,p1,eta1,phi1);
     }
     h_ntrk[0]->Fill(ntrk,tr_eventWeight);
 
     std::vector<spr::propagatedTrackID> trkCaloDets;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, ((verbosity_/100)%10 > 0));
     std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
     for (trkDetItr = trkCaloDets.begin(),ntrk=0; trkDetItr != trkCaloDets.end(); trkDetItr++,ntrk++) {
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
       double pt1         = pTrack->pt();
       double p1          = pTrack->p();
       double eta1        = pTrack->momentum().eta();
       double phi1        = pTrack->momentum().phi();
       if ((verbosity_/10)%10 > 0) 
     edm::LogInfo("IsoTrack") << "track: p " << p1 << " pt " << pt1 
                  << " eta " << eta1 << " phi " << phi1 
                  << " okEcal " << trkDetItr->okECAL;
       fillTrack(2, pt1,p1,eta1,phi1);
       if (pt1>minTrackP_ && std::abs(eta1)<maxTrackEta_ && trkDetItr->okECAL) { 
     fillTrack(3, pt1,p1,eta1,phi1);
     double maxNearP31x31 = spr::chargeIsolationEcal(ntrk, trkCaloDets, geo, caloTopology, 15, 15, ((verbosity_/1000)%10 > 0));
     
     edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
     iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
   
     edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
     edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
     iEvent.getByToken(tok_EB_,barrelRecHitsHandle);
     iEvent.getByToken(tok_EE_,endcapRecHitsHandle);
     // get ECal Tranverse Profile
     std::pair<double, bool>  e7x7P, e11x11P, e15x15P;
     const DetId isoCell = trkDetItr->detIdECAL;
     e7x7P   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3, 0.030, 0.150, tMinE_,tMaxE_, ((verbosity_/10000)%10 > 0));
     e11x11P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, 0.030, 0.150, tMinE_,tMaxE_, ((verbosity_/10000)%10 > 0));
     e15x15P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, 0.030, 0.150, tMinE_,tMaxE_, ((verbosity_/10000)%10 > 0));
 
     double  maxNearHcalP7x7 = spr::chargeIsolationHcal(ntrk, trkCaloDets, theHBHETopology, 3,3, ((verbosity_/1000)%10 > 0));
     int    ieta(0);
     double h3x3(0), h5x5(0), h7x7(0);
     fillIsolation(0, maxNearP31x31,e11x11P.first,e15x15P.first);
     if ((verbosity_/10)%10 > 0) 
       edm::LogInfo("IsoTrack") << "Accepted Tracks reaching Ecal maxNearP31x31 " 
                    << maxNearP31x31 << " e11x11P " 
                    << e11x11P.first << " e15x15P " 
                    << e15x15P.first << " okHCAL " 
                    << trkDetItr->okHCAL;
 
     int trackID = trackPID(pTrack,genParticles);
     if (trkDetItr->okHCAL) {
       edm::Handle<HBHERecHitCollection> hbhe;
       iEvent.getByToken(tok_hbhe_, hbhe);
       const DetId ClosestCell(trkDetItr->detIdHCAL);
       ieta = ((HcalDetId)(ClosestCell)).ietaAbs();
       h3x3 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,1,1, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_, ((verbosity_/10000)%10 > 0));  
       h5x5 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,2,2, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_, ((verbosity_/10000)%10 > 0) );  
       h7x7 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,3,3, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_, ((verbosity_/10000)%10 > 0) );  
       fillIsolation(1, maxNearHcalP7x7,h5x5,h7x7);
       if ((verbosity_/10)%10 > 0)
         edm::LogInfo("IsoTrack") << "Tracks Reaching Hcal maxNearHcalP7x7/h5x5/h7x7 " 
                      << maxNearHcalP7x7 << "/" << h5x5 << "/" << h7x7;
       tr_TrkPt.push_back(pt1);
       tr_TrkP.push_back(p1);
       tr_TrkEta.push_back(eta1);
       tr_TrkPhi.push_back(phi1);
       tr_TrkID.push_back(trackID);
       tr_MaxNearP31X31.push_back(maxNearP31x31);
       tr_MaxNearHcalP7x7.push_back(maxNearHcalP7x7);
       tr_FE7x7P.push_back(e7x7P.first);
       tr_FE11x11P.push_back(e11x11P.first);
       tr_FE15x15P.push_back(e15x15P.first);
       tr_SE7x7P.push_back(e7x7P.second);
       tr_SE11x11P.push_back(e11x11P.second);
       tr_SE15x15P.push_back(e15x15P.second);
       tr_iEta.push_back(ieta);
       tr_H3x3.push_back(h3x3);
       tr_H5x5.push_back(h5x5);
       tr_H7x7.push_back(h7x7);
     
       if (maxNearP31x31 < 0) {
         fillTrack(4, pt1,p1,eta1,phi1);
         fillEnergy(0,ieta,p1,e7x7P.first,h3x3,e11x11P.first,h5x5);
         if (maxNearHcalP7x7 < 0) {
           fillTrack(5, pt1,p1,eta1,phi1);
           fillEnergy(1,ieta,p1,e7x7P.first,h3x3,e11x11P.first,h5x5);
           if ((e11x11P.second) && (e15x15P.second) &&
           (e15x15P.first-e11x11P.first)<2.0) {
         fillTrack(6, pt1,p1,eta1,phi1);
         fillEnergy(2,ieta,p1,e7x7P.first,h3x3,e11x11P.first,h5x5);
         if (h7x7-h5x5 < 2.0) {
           fillTrack(7, pt1,p1,eta1,phi1);
           fillEnergy(3,ieta,p1,e7x7P.first,h3x3,e11x11P.first,h5x5);
           if (nPV >= 0) {
             fillTrack(nPV+8, pt1,p1,eta1,phi1);
             fillEnergy(nPV+4,ieta,p1,e7x7P.first,h3x3,e11x11P.first,h5x5);
           }
           if (trackID > 0) {
             fillTrack(nPVBin_+trackID+7, pt1,p1,eta1,phi1);
             fillEnergy(nPVBin_+trackID+3,ieta,p1,e7x7P.first,h3x3,e11x11P.first,h5x5);
           }
         }
           }
         }
       }
     }
       }
     }
     h_ntrk[1]->Fill(ntrk,tr_eventWeight);
     if ((!tr_TrkPt.empty()) && doTree_) tree_->Fill();
   }
   firstEvent_ = false;
 }

void StudyHLT::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 584 of file StudyHLT.cc.

References doTree_, etaBin_, fs_, h_ediff, h_ene1, h_ene2, h_energy, h_eta, h_goodPV, h_goodRun, h_HLTAccept, h_HLTCorr, h_maxNearP, h_nHLT, h_nHLTvsRN, h_ntrk, h_numberPV, h_p, h_phi, h_pt, mps_fire::i, TFileService::make(), nEtaBin_, nGen_, nPBin_, nPVBin_, pBin_, pvBin_, AlCaHLTBitMon_QueryRunRegistry::string, tr_eventWeight, tr_FE11x11P, tr_FE15x15P, tr_FE7x7P, tr_goodPV, tr_goodRun, tr_H3x3, tr_H5x5, tr_H7x7, tr_iEta, tr_MaxNearHcalP7x7, tr_MaxNearP31X31, tr_SE11x11P, tr_SE15x15P, tr_SE7x7P, tr_TrigName, tr_TrkEta, tr_TrkID, tr_TrkP, tr_TrkPhi, tr_TrkPt, and tree_.

Referenced by ~StudyHLT().

                         {
   // Book histograms
   h_nHLT        = fs_->make<TH1I>("h_nHLT" , "size of trigger Names", 1000, 0, 1000);
   h_HLTAccept   = fs_->make<TH1I>("h_HLTAccept", "HLT Accepts for all runs", 500, 0, 500);
   for (int i=1; i<=500; ++i) h_HLTAccept->GetXaxis()->SetBinLabel(i," ");
   h_nHLTvsRN    = fs_->make<TH2I>("h_nHLTvsRN" , "size of trigger Names vs RunNo", 2168, 190949, 193116, 100, 400, 500);
   h_HLTCorr     = fs_->make<TH1I>("h_HLTCorr", "Correlation among different paths", 100, 0, 100);
   h_numberPV    = fs_->make<TH1I>("h_numberPV", "Number of Primary Vertex", 100, 0, 100);
   h_goodPV      = fs_->make<TH1I>("h_goodPV", "Number of good Primary Vertex", 100, 0, 100);
   h_goodRun     = fs_->make<TH1I>("h_goodRun","Number of accepted events for Run", 4000, 190000, 1940000);
   char hname[60], htit[200];
   std::string CollectionNames[2] = {"Reco", "Propagated"};
   for (unsigned int i=0; i<2; i++) {
     sprintf(hname, "h_nTrk_%s", CollectionNames[i].c_str());
     sprintf(htit, "Number of %s tracks", CollectionNames[i].c_str());
     h_ntrk[i] = fs_->make<TH1I>(hname, htit, 500, 0, 500);
   }
   std::string TrkNames[8]       = {"All", "Quality", "NoIso", "okEcal", "EcalCharIso", "HcalCharIso", "EcalNeutIso", "HcalNeutIso"};
   std::string particle[4]       = {"Electron", "Pion", "Kaon", "Proton"};
   for (unsigned int i=0; i<=nGen_+1; i++) {
     if (i < 8) {
       sprintf(hname, "h_pt_%s", TrkNames[i].c_str());
       sprintf(htit, "p_{T} of %s tracks", TrkNames[i].c_str());
     } else if (i < 8+nPVBin_) {
       sprintf(hname, "h_pt_%s_%d", TrkNames[7].c_str(), i-8);
       sprintf(htit, "p_{T} of %s tracks (PV=%d:%d)", TrkNames[7].c_str(), pvBin_[i-8], pvBin_[i-7]-1);
     } else if (i >= nGen_) {
       sprintf(hname, "h_pt_%s_%d", TrkNames[0].c_str(), i-nGen_);
       sprintf(htit, "p_{T} of %s Generator tracks", TrkNames[0].c_str());
     } else {
       sprintf(hname, "h_pt_%s_%s", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
       sprintf(htit, "p_{T} of %s tracks (%s)", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
     }
     h_pt[i]   = fs_->make<TH1D>(hname, htit, 400, 0, 200.0);
     h_pt[i]->Sumw2();
 
     if (i < 8) {
       sprintf(hname, "h_p_%s", TrkNames[i].c_str());
       sprintf(htit, "Momentum of %s tracks", TrkNames[i].c_str());
     } else if (i < 8+nPVBin_) {
       sprintf(hname, "h_p_%s_%d", TrkNames[7].c_str(), i-8);
       sprintf(htit, "Momentum of %s tracks (PV=%d:%d)", TrkNames[7].c_str(), pvBin_[i-8], pvBin_[i-7]-1);
     } else if (i >= nGen_) {
       sprintf(hname, "h_p_%s_%d", TrkNames[0].c_str(), i-nGen_);
       sprintf(htit, "Momentum of %s Generator tracks", TrkNames[0].c_str());
     } else {
       sprintf(hname, "h_p_%s_%s", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
       sprintf(htit, "Momentum of %s tracks (%s)", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
     }
     h_p[i]    = fs_->make<TH1D>(hname, htit, 400, 0, 200.0);
     h_p[i]->Sumw2();
 
     if (i < 8) {
       sprintf(hname, "h_eta_%s", TrkNames[i].c_str());
       sprintf(htit, "Eta of %s tracks", TrkNames[i].c_str());
     } else if (i < 8+nPVBin_) {
       sprintf(hname, "h_eta_%s_%d", TrkNames[7].c_str(), i-8);
       sprintf(htit, "Eta of %s tracks (PV=%d:%d)", TrkNames[7].c_str(), pvBin_[i-8], pvBin_[i-7]-1);
     } else if (i >= nGen_) {
       sprintf(hname, "h_eta_%s_%d", TrkNames[0].c_str(), i-nGen_);
       sprintf(htit, "Eta of %s Generator tracks", TrkNames[0].c_str());
     } else {
       sprintf(hname, "h_eta_%s_%s", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
       sprintf(htit, "Eta of %s tracks (%s)", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
     }
     h_eta[i]  = fs_->make<TH1D>(hname, htit, 60, -3.0, 3.0);
     h_eta[i]->Sumw2();
 
     if (i < 8) {
       sprintf(hname, "h_phi_%s", TrkNames[i].c_str());
       sprintf(htit, "Phi of %s tracks", TrkNames[i].c_str());
     } else if (i < 8+nPVBin_) {
       sprintf(hname, "h_phi_%s_%d", TrkNames[7].c_str(), i-8);
       sprintf(htit, "Phi of %s tracks (PV=%d:%d)", TrkNames[7].c_str(), pvBin_[i-8], pvBin_[i-7]-1);
     } else if (i >= nGen_) {
       sprintf(hname, "h_phi_%s_%d", TrkNames[0].c_str(), i-nGen_);
       sprintf(htit, "Phi of %s Generator tracks", TrkNames[0].c_str());
     } else {
       sprintf(hname, "h_phi_%s_%s", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
       sprintf(htit, "Phi of %s tracks (%s)", TrkNames[7].c_str(), particle[i-8-nPVBin_].c_str());
     }
     h_phi[i]  = fs_->make<TH1D>(hname, htit, 100, -3.15, 3.15);
     h_phi[i]->Sumw2();
   }
   std::string IsolationNames[2] = {"Ecal", "Hcal"};
   for (unsigned int i=0; i<2; i++) {
     sprintf(hname, "h_maxNearP_%s", IsolationNames[i].c_str());
     sprintf(htit, "Energy in ChargeIso region for %s", IsolationNames[i].c_str());
     h_maxNearP[i] = fs_->make<TH1D>(hname, htit, 120, -1.5, 10.5);
     h_maxNearP[i]->Sumw2(); 
 
     sprintf(hname, "h_ene1_%s", IsolationNames[i].c_str());
     sprintf(htit, "Energy in smaller cone for %s", IsolationNames[i].c_str());
     h_ene1[i]     = fs_->make<TH1D>(hname, htit, 400, 0.0, 200.0);
     h_ene1[i]->Sumw2();
 
     sprintf(hname, "h_ene2_%s", IsolationNames[i].c_str());
     sprintf(htit, "Energy in bigger cone for %s", IsolationNames[i].c_str());
     h_ene2[i]     = fs_->make<TH1D>(hname, htit, 400, 0.0, 200.0);
     h_ene2[i]->Sumw2(); 
 
     sprintf(hname, "h_ediff_%s", IsolationNames[i].c_str());
     sprintf(htit, "Energy in NeutralIso region for %s", IsolationNames[i].c_str());
     h_ediff[i]      = fs_->make<TH1D>(hname, htit, 100, -0.5, 19.5);
     h_ediff[i]->Sumw2();
   }
   std::string energyNames[6]={"E_{7x7}", "H_{3x3}", "(E_{7x7}+H_{3x3})",
                   "E_{11x11}", "H_{5x5}", "{E_{11x11}+H_{5x5})"};
   for (int i=0; i<4+nPVBin_+4; ++i) {
     for (int ip=0; ip<nPBin_; ++ip) {
       for (int ie=0; ie<nEtaBin_; ++ie) {
     for (int j=0; j<6; ++j) {
       sprintf(hname, "h_energy_%d_%d_%d_%d", i, ip, ie, j);
       if (i < 4) {
         sprintf(htit,"%s/p (p=%4.1f:%4.1f; i#eta=%d:%d) for tracks with %s",
             energyNames[j].c_str(),pBin_[ip],pBin_[ip+1],etaBin_[ie],
             (etaBin_[ie+1]-1), TrkNames[i+4].c_str());
       } else if (i < 4+nPVBin_) {
         sprintf(htit,"%s/p (p=%4.1f:%4.1f, i#eta=%d:%d, PV=%d:%d) for tracks with %s",
             energyNames[j].c_str(),pBin_[ip],pBin_[ip+1],etaBin_[ie],
             (etaBin_[ie+1]-1), pvBin_[i-4], pvBin_[i-3],
             TrkNames[7].c_str());
       } else {
         sprintf(htit,"%s/p (p=%4.1f:%4.1f, i#eta=%d:%d %s) for tracks with %s",
             energyNames[j].c_str(),pBin_[ip],pBin_[ip+1],etaBin_[ie],
             (etaBin_[ie+1]-1), particle[i-4-nPVBin_].c_str(),
             TrkNames[7].c_str());
       }
       h_energy[i][ip][ie][j] = fs_->make<TH1D>(hname, htit, 500, -0.1, 4.9);
       h_energy[i][ip][ie][j]->Sumw2();
     }
       }
     }
   }
 
   // Now the tree
   if (doTree_) {
     tree_ = fs_->make<TTree>("testTree", "new HLT Tree");
     tree_->Branch("tr_goodRun",         &tr_goodRun,        "tr_goodRun/I");
     tree_->Branch("tr_goodPV",          &tr_goodPV,         "tr_goodPV/I");
     tree_->Branch("tr_eventWeight",     &tr_eventWeight,    "tr_eventWeight/D");
     tree_->Branch("tr_tr_TrigName",     &tr_TrigName);
     tree_->Branch("tr_TrkPt",           &tr_TrkPt);
     tree_->Branch("tr_TrkP",            &tr_TrkP);
     tree_->Branch("tr_TrkEta",          &tr_TrkEta);
     tree_->Branch("tr_TrkPhi",          &tr_TrkPhi);
     tree_->Branch("tr_TrkID",           &tr_TrkID);
     tree_->Branch("tr_MaxNearP31X31",   &tr_MaxNearP31X31);   
     tree_->Branch("tr_MaxNearHcalP7x7", &tr_MaxNearHcalP7x7);
     tree_->Branch("tr_FE7x7P",          &tr_FE7x7P);
     tree_->Branch("tr_FE11x11P",        &tr_FE11x11P);
     tree_->Branch("tr_FE15x15P",        &tr_FE15x15P);
     tree_->Branch("tr_SE7x7P",          &tr_SE7x7P);
     tree_->Branch("tr_SE11x11P",        &tr_SE11x11P);
     tree_->Branch("tr_SE15x15P",        &tr_SE15x15P);
     tree_->Branch("tr_H3x3",            &tr_H3x3);
     tree_->Branch("tr_H5x5",            &tr_H5x5);
     tree_->Branch("tr_H7x7",            &tr_H7x7);
     tree_->Branch("tr_iEta",            &tr_iEta);
   }
 }

virtual void StudyHLT::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

inlineprivatevirtual

Definition at line 95 of file StudyHLT.cc.

95 { }

void StudyHLT::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprivate

Definition at line 747 of file StudyHLT.cc.

References changed_, firstEvent_, fs_, h_HLTAccepts, hltConfig_, mps_fire::i, HLTConfigProvider::init(), TFileService::make(), nRun_, and edm::RunBase::run().

Referenced by ~StudyHLT().

                                                                        {
   char hname[100], htit[400];
   edm::LogInfo("IsoTrack")  << "Run[" << nRun_ << "] " << iRun.run() << " hltconfig.init " 
                 << hltConfig_.init(iRun,iSetup,"HLT",changed_);
   sprintf(hname, "h_HLTAccepts_%i", iRun.run());
   sprintf(htit, "HLT Accepts for Run No %i", iRun.run());
   TH1I *hnew = fs_->make<TH1I>(hname, htit, 500, 0, 500);
   for (int i=1; i<=500; ++i) hnew->GetXaxis()->SetBinLabel(i," ");
   h_HLTAccepts.push_back(hnew);
   edm::LogInfo("IsoTrack") << "beginrun " << iRun.run();
   firstEvent_ = true;
   changed_    = false;
 }

void StudyHLT::clear ( void )

private

Definition at line 767 of file StudyHLT.cc.

References tr_FE11x11P, tr_FE15x15P, tr_FE7x7P, tr_H3x3, tr_H5x5, tr_H7x7, tr_iEta, tr_MaxNearHcalP7x7, tr_MaxNearP31X31, tr_SE11x11P, tr_SE15x15P, tr_SE7x7P, tr_TrigName, tr_TrkEta, tr_TrkID, tr_TrkP, tr_TrkPhi, and tr_TrkPt.

Referenced by analyze(), Vispa.Views.WidgetView.WidgetView::closeEvent(), Vispa.Views.BoxDecayView.BoxDecayView::closeEvent(), endLuminosityBlock(), Vispa.Share.FindAlgorithm.FindAlgorithm::findUsingFindDialog(), Vispa.Views.LineDecayView.LineDecayView::setDataObjects(), Vispa.Views.WidgetView.WidgetView::setDataObjects(), BeautifulSoup.Tag::setString(), Vispa.Views.TreeView.TreeView::updateContent(), Vispa.Views.TableView.TableView::updateContent(), Vispa.Views.BoxDecayView.BoxDecayView::updateContent(), and Vispa.Views.PropertyView.PropertyView::updateContent().

                      {
   tr_TrigName.clear();
   tr_TrkPt.clear();  tr_TrkP.clear();  tr_TrkEta.clear();  tr_TrkPhi.clear();
   tr_TrkID.clear();  tr_MaxNearP31X31.clear();  tr_MaxNearHcalP7x7.clear(); 
   tr_FE7x7P.clear(); tr_FE11x11P.clear();  tr_FE15x15P.clear();
   tr_SE7x7P.clear(); tr_SE11x11P.clear();  tr_SE15x15P.clear();
   tr_H3x3.clear();   tr_H5x5.clear();      tr_H7x7.clear();
   tr_iEta.clear();
 }

virtual void StudyHLT::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

inlineprivatevirtual

Definition at line 96 of file StudyHLT.cc.

References clear(), fillEnergy(), fillIsolation(), fillTrack(), AlCaHLTBitMon_QueryRunRegistry::string, trackPID(), and truncate_str().

96 { }

void StudyHLT::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &
	)

overrideprivate

Definition at line 762 of file StudyHLT.cc.

References nRun_, and edm::RunBase::run().

Referenced by ~StudyHLT().

                                                               {
   ++nRun_;
   edm::LogInfo("IsoTrack") << "endrun[" << nRun_ << "] " << iRun.run();
 }

void StudyHLT::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 237 of file StudyHLT.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), AlCaHLTBitMon_QueryRunRegistry::string, and create_public_lumi_plots::weights.

Referenced by ~StudyHLT().

                                                                           {
   std::vector<std::string> trig;
   std::vector<double>      weights;
   std::vector<std::string> newNames = {"HLT", "PixelTracks_Multiplicity",
                        "HLT_Physics_","HLT_JetE",
                        "HLT_ZeroBias"};
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("particleSource",edm::InputTag("genParticles"));
   desc.addUntracked<int>("verbosity",         0);
   desc.addUntracked<std::vector<std::string> >("triggers",trig);
   desc.addUntracked<std::vector<std::string> >("newNames",newNames);
   desc.addUntracked<std::string>("trackQuality","highPurity");
   desc.addUntracked<double>("minTrackPt",     1.0);
   desc.addUntracked<double>("maxDxyPV",       0.02);
   desc.addUntracked<double>("maxDzPV",        0.02);
   desc.addUntracked<double>("maxChi2",        5.0);
   desc.addUntracked<double>("maxDpOverP",     0.1);
   desc.addUntracked<int>("minOuterHit",       4);
   desc.addUntracked<int>("minLayerCrossed",   8);
   desc.addUntracked<int>("maxInMiss",         0);
   desc.addUntracked<int>("maxOutMiss",        0);
   desc.addUntracked<double>("minTrackP",      1.0);
   desc.addUntracked<double>("maxTrackEta",    2.6);
   desc.addUntracked<double>("timeMinCutECAL",-500.0);
   desc.addUntracked<double>("timeMaxCutECAL", 500.0);
   desc.addUntracked<double>("timeMinCutHCAL",-500.0);
   desc.addUntracked<double>("timeMaxCutHCAL", 500.0);
   desc.addUntracked<bool>("isItAOD",          false);
   desc.addUntracked<bool>("doTree",           false);
   desc.addUntracked<std::vector<double> >("puWeights", weights);
   descriptions.add("studyHLT",desc);
 }

void StudyHLT::fillEnergy	(	int	flag,
		int	ieta,
		double	p,
		double	enEcal1,
		double	enHcal1,
		double	enEcal2,
		double	enHcal2
	)

private

Definition at line 791 of file StudyHLT.cc.

References etaBin_, RemoveAddSevLevel::flag, h_energy, mps_fire::i, nEtaBin_, nPBin_, pBin_, and tr_eventWeight.

Referenced by analyze(), and endLuminosityBlock().

                                                               {
   int ip(-1), ie(-1);
   for (int i=0; i<nPBin_; ++i) {
     if (p >= pBin_[i] && p < pBin_[i+1]) { ip = i; break; }
   }
   for (int i=0; i<nEtaBin_; ++i) {
     if (ieta >= etaBin_[i] && ieta < etaBin_[i+1]) { ie = i; break; }
   }
   if (ip >= 0 && ie >= 0 && enEcal1 > 0.02 && enHcal1 > 0.1) {
     h_energy[flag][ip][ie][0]->Fill(enEcal1/p,tr_eventWeight);
     h_energy[flag][ip][ie][1]->Fill(enHcal1/p,tr_eventWeight);
     h_energy[flag][ip][ie][2]->Fill((enEcal1+enHcal1)/p,tr_eventWeight);
     h_energy[flag][ip][ie][3]->Fill(enEcal2/p,tr_eventWeight);
     h_energy[flag][ip][ie][4]->Fill(enHcal2/p,tr_eventWeight);
     h_energy[flag][ip][ie][5]->Fill((enEcal2+enHcal2)/p,tr_eventWeight);
   }
 }

void StudyHLT::fillIsolation	(	int	i,
		double	emaxnearP,
		double	eneutIso1,
		double	eneutIso2
	)

private

Definition at line 784 of file StudyHLT.cc.

References h_ediff, h_ene1, h_ene2, h_maxNearP, mps_fire::i, and tr_eventWeight.

Referenced by analyze(), and endLuminosityBlock().

                                                                                        {
   h_maxNearP[i]->Fill(emaxnearP,tr_eventWeight);
   h_ene1[i]->Fill(eneutIso1,tr_eventWeight);
   h_ene2[i]->Fill(eneutIso2,tr_eventWeight);
   h_ediff[i]->Fill(eneutIso2-eneutIso1,tr_eventWeight);
 }

void StudyHLT::fillTrack	(	int	i,
		double	pt,
		double	p,
		double	eta,
		double	phi
	)

private

Definition at line 777 of file StudyHLT.cc.

References h_eta, h_p, h_phi, h_pt, mps_fire::i, and tr_eventWeight.

Referenced by analyze(), and endLuminosityBlock().

                                                                           {
   h_pt[i]->Fill(pt,tr_eventWeight);
   h_p[i]->Fill(p,tr_eventWeight);
   h_eta[i]->Fill(eta,tr_eventWeight);
   h_phi[i]->Fill(phi,tr_eventWeight);
 }

int StudyHLT::trackPID	(	const reco::Track *	pTrack,
		const edm::Handle< reco::GenParticleCollection > &	genParticles
	)

private

Definition at line 822 of file StudyHLT.cc.

References funct::abs(), DEFINE_FWK_MODULE, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::TrackBase::eta(), triggerObjects_cff::id, training_settings::idx, edm::HandleBase::isValid(), M_PI, AlCaHLTBitMon_ParallelJobs::p, common_cff::pdgId, reco::TrackBase::phi(), and mathSSE::sqrt().

Referenced by analyze(), and endLuminosityBlock().

                                                                                                          {
   int id(0);
   if (genParticles.isValid()) {
     unsigned int indx;
     reco::GenParticleCollection::const_iterator p;
     double mindR(999.9);
     for (p=genParticles->begin(),indx=0; p!=genParticles->end(); ++p,++indx) {
       int pdgId  = std::abs(p->pdgId());
       int idx    = (pdgId == 11) ? 1 :
     ((pdgId == 211) ? 2 : ((pdgId == 321) ? 3 : ((pdgId == 2212) ? 4 : 0)));
       if (idx > 0) {
     double dEta = pTrack->eta() - p->momentum().Eta();
     double phi1 = pTrack->phi();
     double phi2 = p->momentum().Phi();
     if (phi1 < 0) phi1 += 2.0*M_PI;
     if (phi2 < 0) phi2 += 2.0*M_PI;
     double dPhi = phi1-phi2;
     if (dPhi > M_PI)       dPhi -= 2.*M_PI;
     else if (dPhi < -M_PI) dPhi += 2.*M_PI;
     double dR   = sqrt(dEta*dEta+dPhi*dPhi);
     if (dR < mindR) {
       mindR = dR; id = idx;
     }
       }
     }
   }
   return id;
 }

std::string StudyHLT::truncate_str ( const std::string & str )

private

Definition at line 810 of file StudyHLT.cc.

References mps_fire::i, AlCaHLTBitMon_QueryRunRegistry::string, and z.

Referenced by analyze(), and endLuminosityBlock().

                                                      {
   std::string truncated_str(str);
   int length = str.length();
   for (int i=0; i<length-2; i++){
     if (str[i]=='_' && str[i+1]=='v' && isdigit(str.at(i+2))) {
       int z = i+1;
       truncated_str = str.substr(0,z);
     } 
   }
   return(truncated_str);
 }

Member Data Documentation

bool StudyHLT::changed_

private

Definition at line 121 of file StudyHLT.cc.

Referenced by beginRun(), and StudyHLT().

const bool StudyHLT::doTree_

private

Definition at line 116 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and StudyHLT().

int StudyHLT::etaBin_[nEtaBin_+1]

private

Definition at line 144 of file StudyHLT.cc.

Referenced by beginJob(), fillEnergy(), and StudyHLT().

bool StudyHLT::firstEvent_

private

Definition at line 121 of file StudyHLT.cc.

Referenced by analyze(), beginRun(), and StudyHLT().

edm::Service<TFileService> StudyHLT::fs_

private

Definition at line 110 of file StudyHLT.cc.

Referenced by beginJob(), and beginRun().

TH1D * StudyHLT::h_ediff[2]

private

Definition at line 141 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1D * StudyHLT::h_ene1[2]

private

Definition at line 141 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1D * StudyHLT::h_ene2[2]

private

Definition at line 141 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1D* StudyHLT::h_energy[nPVBin_+8][nPBin_][nEtaBin_][6]

private

Definition at line 142 of file StudyHLT.cc.

Referenced by beginJob(), and fillEnergy().

TH1D* StudyHLT::h_eta[nGen_+2]

private

Definition at line 139 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

TH1I* StudyHLT::h_goodPV

private

Definition at line 135 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I * StudyHLT::h_goodRun

private

Definition at line 135 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I * StudyHLT::h_HLTAccept

private

Definition at line 134 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

std::vector<TH1I*> StudyHLT::h_HLTAccepts

private

Definition at line 137 of file StudyHLT.cc.

Referenced by analyze(), and beginRun().

TH1I * StudyHLT::h_HLTCorr

private

Definition at line 134 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1D* StudyHLT::h_maxNearP[2]

private

Definition at line 141 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1I* StudyHLT::h_nHLT

private

Definition at line 134 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH2I* StudyHLT::h_nHLTvsRN

private

Definition at line 136 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I* StudyHLT::h_ntrk[2]

private

Definition at line 140 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I * StudyHLT::h_numberPV

private

Definition at line 134 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1D* StudyHLT::h_p[nGen_+2]

private

Definition at line 138 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

TH1D * StudyHLT::h_phi[nGen_+2]

private

Definition at line 139 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

TH1D * StudyHLT::h_pt[nGen_+2]

private

Definition at line 138 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

HLTConfigProvider StudyHLT::hltConfig_

private

Definition at line 109 of file StudyHLT.cc.

Referenced by beginRun().

std::vector<std::string> StudyHLT::HLTNames_

private

Definition at line 120 of file StudyHLT.cc.

Referenced by analyze().

const bool StudyHLT::isItAOD_

private

Definition at line 116 of file StudyHLT.cc.

Referenced by StudyHLT().

const double StudyHLT::maxTrackEta_

private

Definition at line 114 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const double StudyHLT::minTrackP_

private

Definition at line 114 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const int StudyHLT::nEtaBin_ =4

staticprivate

Definition at line 107 of file StudyHLT.cc.

Referenced by beginJob(), fillEnergy(), and StudyHLT().

const std::vector<std::string> StudyHLT::newNames_

private

Definition at line 112 of file StudyHLT.cc.

Referenced by analyze().

const int StudyHLT::nGen_ =(nPVBin_+12)

staticprivate

Definition at line 108 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

const int StudyHLT::nPBin_ =10

staticprivate

Definition at line 107 of file StudyHLT.cc.

Referenced by beginJob(), fillEnergy(), and StudyHLT().

const int StudyHLT::nPVBin_ =4

staticprivate

Definition at line 107 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and StudyHLT().

int StudyHLT::nRun_

private

Definition at line 144 of file StudyHLT.cc.

Referenced by analyze(), beginRun(), and endRun().

double StudyHLT::pBin_[nPBin_+1]

private

Definition at line 145 of file StudyHLT.cc.

Referenced by beginJob(), fillEnergy(), and StudyHLT().

const std::vector<double> StudyHLT::puWeights_

private

Definition at line 117 of file StudyHLT.cc.

Referenced by analyze().

int StudyHLT::pvBin_[nPVBin_+1]

private

Definition at line 144 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and StudyHLT().

spr::trackSelectionParameters StudyHLT::selectionParameters_

private

Definition at line 119 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const std::string StudyHLT::theTrackQuality_

private

Definition at line 113 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const edm::InputTag StudyHLT::theTriggerResultsLabel_

private

Definition at line 118 of file StudyHLT.cc.

Referenced by StudyHLT().

const double StudyHLT::tMaxE_

private

Definition at line 115 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const double StudyHLT::tMaxH_

private

Definition at line 115 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const double StudyHLT::tMinE_

private

Definition at line 115 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const double StudyHLT::tMinH_

private

Definition at line 115 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<EcalRecHitCollection> StudyHLT::tok_EB_

private

Definition at line 129 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<EcalRecHitCollection> StudyHLT::tok_EE_

private

Definition at line 130 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<GenEventInfoProduct> StudyHLT::tok_ew_

private

Definition at line 132 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<reco::TrackCollection> StudyHLT::tok_genTrack_

private

Definition at line 127 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<HBHERecHitCollection> StudyHLT::tok_hbhe_

private

Definition at line 131 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<LumiDetails> StudyHLT::tok_lumi

private

Definition at line 123 of file StudyHLT.cc.

Referenced by StudyHLT().

edm::EDGetTokenT<reco::GenParticleCollection> StudyHLT::tok_parts_

private

Definition at line 126 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<reco::VertexCollection> StudyHLT::tok_recVtx_

private

Definition at line 128 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<trigger::TriggerEvent> StudyHLT::tok_trigEvt

private

Definition at line 124 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<edm::TriggerResults> StudyHLT::tok_trigRes

private

Definition at line 125 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

double StudyHLT::tr_eventWeight

private

Definition at line 147 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), fillEnergy(), fillIsolation(), and fillTrack().

std::vector<double> StudyHLT::tr_FE11x11P

private

Definition at line 152 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_FE15x15P

private

Definition at line 152 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_FE7x7P

private

Definition at line 152 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

int StudyHLT::tr_goodPV

private

Definition at line 146 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

int StudyHLT::tr_goodRun

private

Definition at line 146 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

std::vector<double> StudyHLT::tr_H3x3

private

Definition at line 151 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_H5x5

private

Definition at line 151 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_H7x7

private

Definition at line 151 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<int> StudyHLT::tr_iEta

private

Definition at line 154 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_MaxNearHcalP7x7

private

Definition at line 150 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_MaxNearP31X31

private

Definition at line 150 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<bool> StudyHLT::tr_SE11x11P

private

Definition at line 153 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<bool> StudyHLT::tr_SE15x15P

private

Definition at line 153 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<bool> StudyHLT::tr_SE7x7P

private

Definition at line 153 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<std::string> StudyHLT::tr_TrigName

private

Definition at line 148 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_TrkEta

private

Definition at line 149 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<int> StudyHLT::tr_TrkID

private

Definition at line 154 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_TrkP

private

Definition at line 149 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_TrkPhi

private

Definition at line 149 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

std::vector<double> StudyHLT::tr_TrkPt

private

Definition at line 149 of file StudyHLT.cc.

Referenced by analyze(), beginJob(), and clear().

TTree* StudyHLT::tree_

private

Definition at line 143 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

const edm::InputTag StudyHLT::triggerEvent_

private

Definition at line 118 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const std::vector<std::string> StudyHLT::trigNames_

private

Definition at line 112 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const int StudyHLT::verbosity_

private

Definition at line 111 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation