Inheritance diagram for StudyHLT:

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

	~StudyHLT ()

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

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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
virtual void	analyze (edm::Event const &, edm::EventSetup const &) override

virtual void	beginJob () override

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

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

void	clear ()

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

virtual 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_

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 [nPVBin+12]

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 [nPVBin+12]

TH1D *	h_phi [nPVBin+12]

TH1D *	h_pt [nPVBin+12]

HLTConfigProvider	hltConfig_

std::vector< std::string >	HLTNames_

bool	isItAOD_

double	maxTrackEta_

double	minTrackP_

int	nRun

double	pBin [nPBin+1]

std::vector< double >	puWeights_

int	pvBin [nPVBin+1]

spr::trackSelectionParameters	selectionParameters_

std::string	theTrackQuality_

edm::InputTag	theTriggerResultsLabel_

double	tMaxE_

double	tMaxH_

double	tMinE_

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_

edm::InputTag	triggerEvent_

std::vector< std::string >	trigNames_

int	verbosity_

Static Private Attributes
static const int	nEtaBin =4

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

Definition at line 64 of file StudyHLT.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 136 of file StudyHLT.cc.

                                                  : nRun(0) {
 
   usesResource("TFileService");
 
   verbosity_                          = iConfig.getUntrackedParameter<int>("Verbosity",0);
   trigNames_                          = iConfig.getUntrackedParameter<std::vector<std::string> >("Triggers");
   theTrackQuality_                    = iConfig.getUntrackedParameter<std::string>("TrackQuality","highPurity");
   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);
   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");
 
   // 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 ( )

Definition at line 215 of file StudyHLT.cc.

215 {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 226 of file StudyHLT.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationEcal(), spr::chargeIsolationHcal(), clear(), doTree_, spr::eECALmatrix(), spr::eHCALmatrix(), edm::EventID::event(), fillEnergy(), fillIsolation(), fillTrack(), firstEvent_, 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, HLTNames_, mps_fire::i, edm::EventBase::id(), createfilelist::int, edm::HandleBase::isValid(), edm::InputTag::label(), mix_E10TeV_L13E31_BX432_cfi::Lumi, edm::EventBase::luminosityBlock(), maxTrackEta_, spr::trackSelectionParameters::minQuality, minTrackP_, reco::TrackBase::momentum(), nPVBin, nRun, convertSQLiteXML::ok, reco::TrackBase::p(), p1, edm::Handle< T >::product(), edm::ESHandle< T >::product(), spr::propagateCALO(), reco::TrackBase::pt(), CAHitQuadrupletGenerator_cfi::pt1, puWeights_, pvBin, btvTracks_cfi::quality, reco::TrackBase::quality(), edm::EventID::run(), selectionParameters_, edm::HLTGlobalStatus::size(), 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().

                                                                           {
   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::string newNames[5]={"HLT","PixelTracks_Multiplicity","HLT_Physics_","HLT_JetE","HLT_ZeroBias"};
   int         newAccept[5];
   for (int i=0; i<5; ++i) newAccept[i] = 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_.size() < 1) {
       ok = true;
     } else {
       for (unsigned int i=0; i<trigNames_.size(); ++i) {
         if (newtriggerName.find(trigNames_[i].c_str())!=std::string::npos) {
           if (verbosity_%10 > 0)  
         edm::LogInfo("IsoTrack") << newtriggerName;
           if (hlt > 0) {
         ok = true;
         tr_TrigName.push_back(newtriggerName);
           }
         }
       }
       for (int i=0; i<5; ++i) {
         if (newtriggerName.find(newNames[i].c_str())!=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 (int i=0; i<5; ++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_.size() > 0) {
       int npbin = h_goodPV->FindBin(ngoodPV);
       if (npbin > 0 && npbin <= (int)(puWeights_.size())) 
     tr_eventWeight *= puWeights_[npbin-1];
       else  
     tr_eventWeight  = 0;
     }
 
     //=== genParticle information
     edm::Handle<reco::GenParticleCollection> genParticles;
     iEvent.getByToken(tok_parts_, genParticles);
 
     edm::Handle<reco::TrackCollection> trkCollection;
     iEvent.getByToken(tok_genTrack_, trkCollection);
     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.size() > 0 && doTree_) tree_->Fill();
   }
   firstEvent_ = false;
 }

void StudyHLT::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 515 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, 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_.

                         {
   // 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[50], 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<8+nPVBin+4; 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 {
       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 {
       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 {
       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 {
       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);
   }
 }

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

privatevirtual

Definition at line 687 of file StudyHLT.cc.

687 {}

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

overrideprivatevirtual

Definition at line 666 of file StudyHLT.cc.

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

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

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

privatevirtual

Definition at line 689 of file StudyHLT.cc.

689 {}

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

overrideprivatevirtual

Definition at line 681 of file StudyHLT.cc.

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

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

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

static

Definition at line 218 of file StudyHLT.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                           {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

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

private

Definition at line 715 of file StudyHLT.cc.

References etaBin, RemoveAddSevLevel::flag, h_energy, mps_fire::i, nEtaBin, nPBin, pBin, and tr_eventWeight.

Referenced by analyze().

                                                               {
   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 708 of file StudyHLT.cc.

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

Referenced by analyze().

                                                                                        {
   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 701 of file StudyHLT.cc.

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

Referenced by analyze().

                                                                           {
   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 746 of file StudyHLT.cc.

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

Referenced by analyze().

                                                                                                          {
   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 734 of file StudyHLT.cc.

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

Referenced by analyze().

                                                      {
   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 99 of file StudyHLT.cc.

Referenced by beginRun(), and StudyHLT().

bool StudyHLT::doTree_

private

Definition at line 99 of file StudyHLT.cc.

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

int StudyHLT::etaBin[nEtaBin+1]

private

Definition at line 123 of file StudyHLT.cc.

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

bool StudyHLT::firstEvent_

private

Definition at line 99 of file StudyHLT.cc.

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

edm::Service<TFileService> StudyHLT::fs_

private

Definition at line 91 of file StudyHLT.cc.

Referenced by beginJob(), and beginRun().

TH1D * StudyHLT::h_ediff[2]

private

Definition at line 120 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1D * StudyHLT::h_ene1[2]

private

Definition at line 120 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1D * StudyHLT::h_ene2[2]

private

Definition at line 120 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1D* StudyHLT::h_energy[nPVBin+8][nPBin][nEtaBin][6]

private

Definition at line 121 of file StudyHLT.cc.

Referenced by beginJob(), and fillEnergy().

TH1D* StudyHLT::h_eta[nPVBin+12]

private

Definition at line 118 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

TH1I* StudyHLT::h_goodPV

private

Definition at line 114 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I * StudyHLT::h_goodRun

private

Definition at line 114 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I * StudyHLT::h_HLTAccept

private

Definition at line 113 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 116 of file StudyHLT.cc.

Referenced by analyze(), and beginRun().

TH1I * StudyHLT::h_HLTCorr

private

Definition at line 113 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1D* StudyHLT::h_maxNearP[2]

private

Definition at line 120 of file StudyHLT.cc.

Referenced by beginJob(), and fillIsolation().

TH1I* StudyHLT::h_nHLT

private

Definition at line 113 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH2I* StudyHLT::h_nHLTvsRN

private

Definition at line 115 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I* StudyHLT::h_ntrk[2]

private

Definition at line 119 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1I * StudyHLT::h_numberPV

private

Definition at line 113 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

TH1D* StudyHLT::h_p[nPVBin+12]

private

Definition at line 117 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

TH1D * StudyHLT::h_phi[nPVBin+12]

private

Definition at line 118 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

TH1D * StudyHLT::h_pt[nPVBin+12]

private

Definition at line 117 of file StudyHLT.cc.

Referenced by beginJob(), and fillTrack().

HLTConfigProvider StudyHLT::hltConfig_

private

Definition at line 90 of file StudyHLT.cc.

Referenced by beginRun().

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

private

Definition at line 94 of file StudyHLT.cc.

Referenced by analyze().

bool StudyHLT::isItAOD_

private

Definition at line 99 of file StudyHLT.cc.

Referenced by StudyHLT().

double StudyHLT::maxTrackEta_

private

Definition at line 97 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

double StudyHLT::minTrackP_

private

Definition at line 97 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

const int StudyHLT::nEtaBin =4

staticprivate

Definition at line 89 of file StudyHLT.cc.

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

const int StudyHLT::nPBin =10

staticprivate

Definition at line 89 of file StudyHLT.cc.

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

const int StudyHLT::nPVBin =4

staticprivate

Definition at line 89 of file StudyHLT.cc.

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

int StudyHLT::nRun

private

Definition at line 123 of file StudyHLT.cc.

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

double StudyHLT::pBin[nPBin+1]

private

Definition at line 124 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::puWeights_

private

Definition at line 96 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

int StudyHLT::pvBin[nPVBin+1]

private

Definition at line 123 of file StudyHLT.cc.

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

spr::trackSelectionParameters StudyHLT::selectionParameters_

private

Definition at line 93 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

std::string StudyHLT::theTrackQuality_

private

Definition at line 95 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::InputTag StudyHLT::theTriggerResultsLabel_

private

Definition at line 101 of file StudyHLT.cc.

Referenced by StudyHLT().

double StudyHLT::tMaxE_

private

Definition at line 98 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

double StudyHLT::tMaxH_

private

Definition at line 98 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

double StudyHLT::tMinE_

private

Definition at line 98 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

double StudyHLT::tMinH_

private

Definition at line 98 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<EcalRecHitCollection> StudyHLT::tok_EB_

private

Definition at line 108 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<EcalRecHitCollection> StudyHLT::tok_EE_

private

Definition at line 109 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<GenEventInfoProduct> StudyHLT::tok_ew_

private

Definition at line 111 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

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

private

Definition at line 106 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<HBHERecHitCollection> StudyHLT::tok_hbhe_

private

Definition at line 110 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

edm::EDGetTokenT<LumiDetails> StudyHLT::tok_lumi

private

Definition at line 102 of file StudyHLT.cc.

Referenced by StudyHLT().

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

private

Definition at line 105 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

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

private

Definition at line 107 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

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

private

Definition at line 103 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

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

private

Definition at line 104 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

double StudyHLT::tr_eventWeight

private

Definition at line 126 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_FE11x11P

private

Definition at line 131 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_FE15x15P

private

Definition at line 131 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_FE7x7P

private

Definition at line 131 of file StudyHLT.cc.

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

int StudyHLT::tr_goodPV

private

Definition at line 125 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

int StudyHLT::tr_goodRun

private

Definition at line 125 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

std::vector<double> StudyHLT::tr_H3x3

private

Definition at line 130 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_H5x5

private

Definition at line 130 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_H7x7

private

Definition at line 130 of file StudyHLT.cc.

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

std::vector<int> StudyHLT::tr_iEta

private

Definition at line 133 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_MaxNearHcalP7x7

private

Definition at line 129 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_MaxNearP31X31

private

Definition at line 129 of file StudyHLT.cc.

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

std::vector<bool> StudyHLT::tr_SE11x11P

private

Definition at line 132 of file StudyHLT.cc.

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

std::vector<bool> StudyHLT::tr_SE15x15P

private

Definition at line 132 of file StudyHLT.cc.

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

std::vector<bool> StudyHLT::tr_SE7x7P

private

Definition at line 132 of file StudyHLT.cc.

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

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

private

Definition at line 127 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_TrkEta

private

Definition at line 128 of file StudyHLT.cc.

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

std::vector<int> StudyHLT::tr_TrkID

private

Definition at line 133 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_TrkP

private

Definition at line 128 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_TrkPhi

private

Definition at line 128 of file StudyHLT.cc.

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

std::vector<double> StudyHLT::tr_TrkPt

private

Definition at line 128 of file StudyHLT.cc.

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

TTree* StudyHLT::tree_

private

Definition at line 122 of file StudyHLT.cc.

Referenced by analyze(), and beginJob().

edm::InputTag StudyHLT::triggerEvent_

private

Definition at line 101 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

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

private

Definition at line 94 of file StudyHLT.cc.

Referenced by analyze(), and StudyHLT().

int StudyHLT::verbosity_

private

Definition at line 92 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