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

Inheritance diagram for StudyCaloResponse:

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

	~StudyCaloResponse () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

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 double	cutEcal_

const double	cutMuon_

const double	cutRatio_

const bool	doTree_

int	etaBin_ [nEtaBin_+1]

bool	firstEvent_

edm::Service< TFileService >	fs_

TH1D *	h_counter [8]

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

TH1D *	h_h_pNew [8]

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 edm::InputTag	labelGenTrack_

const edm::InputTag	labelMuon_

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::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< LumiDetails >	tok_lumi

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_magField_

edm::EDGetTokenT< reco::MuonCollection >	tok_Muon_

edm::EDGetTokenT< reco::GenParticleCollection >	tok_parts_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_topo_

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_

const bool	vetoEcal_

const bool	vetoMuon_

const bool	vetoTrigger_

Static Private Attributes
static const int	nEtaBin_ = 4

static const int	nGen_ = (nPVBin_ + 12)

static const int	nPBin_ = 15

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: Studies single particle response measurements in data/MC

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

Definition at line 85 of file StudyCaloResponse.cc.

Constructor & Destructor Documentation

◆ StudyCaloResponse()

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

explicit

Definition at line 168 of file StudyCaloResponse.cc.

     : verbosity_(iConfig.getUntrackedParameter<int>("verbosity", 0)),
       trigNames_(iConfig.getUntrackedParameter<std::vector<std::string> >("triggers")),
       newNames_(iConfig.getUntrackedParameter<std::vector<std::string> >("newNames")),
       labelMuon_(iConfig.getUntrackedParameter<edm::InputTag>("labelMuon")),
       labelGenTrack_(iConfig.getUntrackedParameter<edm::InputTag>("labelTrack")),
       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)),
       vetoTrigger_(iConfig.getUntrackedParameter<bool>("vetoTrigger", false)),
       doTree_(iConfig.getUntrackedParameter<bool>("doTree", false)),
       vetoMuon_(iConfig.getUntrackedParameter<bool>("vetoMuon", true)),
       vetoEcal_(iConfig.getUntrackedParameter<bool>("vetoEcal", false)),
       cutMuon_(iConfig.getUntrackedParameter<double>("cutMuon", 0.1)),
       cutEcal_(iConfig.getUntrackedParameter<double>("cutEcal", 2.0)),
       cutRatio_(iConfig.getUntrackedParameter<double>("cutRatio", 0.90)),
       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_Muon_ = consumes<reco::MuonCollection>(labelMuon_);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   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"));
  
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
   tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
   tok_topo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
   tok_magField_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
   tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
   tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
  
   edm::LogVerbatim("IsoTrack") << "Verbosity " << verbosity_ << " with " << trigNames_.size() << " triggers:";
   for (unsigned int k = 0; k < trigNames_.size(); ++k)
     edm::LogVerbatim("IsoTrack") << " [" << k << "] " << trigNames_[k];
   edm::LogVerbatim("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_ << " vetoTrigger " << vetoTrigger_ << " vetoMuon "
                                << vetoMuon_ << ":" << cutMuon_ << " vetoEcal " << vetoEcal_ << ":" << cutEcal_ << ":"
                                << cutRatio_;
  
   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, 25.0, 30.0, 40.0, 60.0, 100.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;
 }

◆ ~StudyCaloResponse()

StudyCaloResponse::~StudyCaloResponse ( )

inlineoverride

Definition at line 88 of file StudyCaloResponse.cc.

88 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 303 of file StudyCaloResponse.cc.

                                                                                    {
   clear();
   int counter0[1000] = {0};
   int counter1[1000] = {0};
   int counter2[1000] = {0};
   int counter3[1000] = {0};
   int counter4[1000] = {0};
   int counter5[1000] = {0};
   int counter6[1000] = {0};
   int counter7[1000] = {0};
   if (verbosity_ > 0)
     edm::LogVerbatim("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);
   if (verbosity_ > 0)
     edm::LogVerbatim("IsoTrack") << "RunNo " << RunNo << " EvtNo " << EvtNo << " Lumi " << Lumi << " Bunch " << Bunch;
  
   trigger::TriggerEvent triggerEvent;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   iEvent.getByToken(tok_trigEvt, triggerEventHandle);
  
   bool ok(false);
   std::string triggerUse("None");
   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++) {
         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::LogVerbatim("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::LogVerbatim("IsoTrack") << newtriggerName;
               if (hlt > 0) {
                 if (!ok)
                   triggerUse = newtriggerName;
                 ok = true;
                 tr_TrigName.push_back(newtriggerName);
               }
             }
           }
           if (vetoTrigger_)
             ok = !ok;
           for (unsigned int i = 0; i < newNames_.size(); ++i) {
             if (newtriggerName.find(newNames_[i]) != std::string::npos) {
               if (verbosity_ % 10 > 0)
                 edm::LogVerbatim("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);
     }
   }
   if ((verbosity_ / 10) % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "Trigger check gives " << ok << " with " << triggerUse;
  
   //Look at the tracks
   if (ok) {
     h_goodRun->Fill(RunNo);
     tr_goodRun = RunNo;
     // get handles to calogeometry and calotopology
     const CaloGeometry* geo = &iSetup.getData(tok_geom_);
     const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
     const HcalTopology* theHBHETopology = &iSetup.getData(tok_topo_);
     const MagneticField* bField = &iSetup.getData(tok_magField_);
     const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
  
     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::LogVerbatim("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);
  
     edm::Handle<reco::MuonCollection> muonEventHandle;
     iEvent.getByToken(tok_Muon_, muonEventHandle);
  
     //=== 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);
       if (p1 < 1000) {
         h_h_pNew[0]->Fill(p1);
         ++counter0[(int)(p1)];
       }
     }
     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::LogVerbatim("IsoTrack") << "track: p " << p1 << " pt " << pt1 << " eta " << eta1 << " phi " << phi1
                                      << " okEcal " << trkDetItr->okECAL;
       fillTrack(2, pt1, p1, eta1, phi1);
  
       bool vetoMuon(false);
       double chiGlobal(0), dr(0);
       bool goodGlob(false);
       if (vetoMuon_) {
         if (muonEventHandle.isValid()) {
           for (reco::MuonCollection::const_iterator recMuon = muonEventHandle->begin();
                recMuon != muonEventHandle->end();
                ++recMuon) {
             if (((recMuon->isPFMuon()) && (recMuon->isGlobalMuon() || recMuon->isTrackerMuon())) &&
                 (recMuon->innerTrack()->validFraction() > 0.49) && (recMuon->innerTrack().isNonnull())) {
               chiGlobal = ((recMuon->globalTrack().isNonnull()) ? recMuon->globalTrack()->normalizedChi2() : 999);
               goodGlob = (recMuon->isGlobalMuon() && chiGlobal < 3 &&
                           recMuon->combinedQuality().chi2LocalPosition < 12 && recMuon->combinedQuality().trkKink < 20);
               if (muon::segmentCompatibility(*recMuon) > (goodGlob ? 0.303 : 0.451)) {
                 const reco::Track* pTrack0 = (recMuon->innerTrack()).get();
                 dr = reco::deltaR(pTrack0->eta(), pTrack0->phi(), pTrack->eta(), pTrack->phi());
                 if (dr < cutMuon_) {
                   vetoMuon = true;
                   break;
                 }
               }
             }
           }
         }
       }
       if ((verbosity_ / 10) % 10 > 0)
         edm::LogVerbatim("IsoTrack") << "vetoMuon: " << vetoMuon_ << ":" << vetoMuon << " chi:good:dr " << chiGlobal
                                      << ":" << goodGlob << ":" << dr;
       if (pt1 > minTrackP_ && std::abs(eta1) < maxTrackEta_ && trkDetItr->okECAL && (!vetoMuon)) {
         fillTrack(3, pt1, p1, eta1, phi1);
         double maxNearP31x31 =
             spr::chargeIsolationEcal(ntrk, trkCaloDets, geo, caloTopology, 15, 15, ((verbosity_ / 1000) % 10 > 0));
  
         const EcalSeverityLevelAlgo* sevlv = &iSetup.getData(tok_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,
                                  3,
                                  3,
                                  0.030,
                                  0.150,
                                  tMinE_,
                                  tMaxE_,
                                  ((verbosity_ / 10000) % 10 > 0));
         e11x11P = spr::eECALmatrix(isoCell,
                                    barrelRecHitsHandle,
                                    endcapRecHitsHandle,
                                    *theEcalChStatus,
                                    geo,
                                    caloTopology,
                                    sevlv,
                                    5,
                                    5,
                                    0.030,
                                    0.150,
                                    tMinE_,
                                    tMaxE_,
                                    ((verbosity_ / 10000) % 10 > 0));
         e15x15P = spr::eECALmatrix(isoCell,
                                    barrelRecHitsHandle,
                                    endcapRecHitsHandle,
                                    *theEcalChStatus,
                                    geo,
                                    caloTopology,
                                    sevlv,
                                    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::LogVerbatim("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);
           double eByh = ((e11x11P.second) ? (e11x11P.first / std::max(h3x3, 0.001)) : 0.0);
           bool notAnElec = ((vetoEcal_ && e11x11P.second) ? ((e11x11P.first < cutEcal_) || (eByh < cutRatio_)) : true);
           if ((verbosity_ / 10) % 10 > 0)
             edm::LogVerbatim("IsoTrack") << "Tracks Reaching Hcal maxNearHcalP7x7/h5x5/h7x7 " << maxNearHcalP7x7 << "/"
                                          << h5x5 << "/" << h7x7 << " eByh " << eByh << " notAnElec " << notAnElec;
           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 && notAnElec) {
             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);
                   }
                   if (p1 < 1000) {
                     h_h_pNew[7]->Fill(p1);
                     ++counter7[(int)(p1)];
                   }
                 }
                 if (p1 < 1000) {
                   h_h_pNew[6]->Fill(p1);
                   ++counter6[(int)(p1)];
                 }
               }
               if (p1 < 1000) {
                 h_h_pNew[5]->Fill(p1);
                 ++counter5[(int)(p1)];
               }
             }
             if (p1 < 1000) {
               h_h_pNew[4]->Fill(p1);
               ++counter4[(int)(p1)];
             }
           }
           if (p1 < 1000) {
             h_h_pNew[3]->Fill(p1);
             ++counter3[(int)(p1)];
           }
         }
         if (p1 < 1000) {
           h_h_pNew[2]->Fill(p1);
           ++counter2[(int)(p1)];
         }
       }
       if (p1 < 1000) {
         h_h_pNew[1]->Fill(p1);
         ++counter1[(int)(p1)];
       }
     }
     h_ntrk[1]->Fill(ntrk, tr_eventWeight);
     if ((!tr_TrkPt.empty()) && doTree_)
       tree_->Fill();
     for (int i = 0; i < 1000; ++i) {
       if (counter0[i])
         h_counter[0]->Fill(i, counter0[i]);
       if (counter1[i])
         h_counter[1]->Fill(i, counter1[i]);
       if (counter2[i])
         h_counter[2]->Fill(i, counter2[i]);
       if (counter3[i])
         h_counter[3]->Fill(i, counter3[i]);
       if (counter4[i])
         h_counter[4]->Fill(i, counter4[i]);
       if (counter5[i])
         h_counter[5]->Fill(i, counter5[i]);
       if (counter6[i])
         h_counter[6]->Fill(i, counter6[i]);
       if (counter7[i])
         h_counter[7]->Fill(i, counter7[i]);
     }
   }
   firstEvent_ = false;
 }

References funct::abs(), Calorimetry_cff::bField, spr::chargeIsolationEcal(), spr::chargeIsolationHcal(), clear(), cutEcal_, cutMuon_, cutRatio_, reco::deltaR(), HLT_FULL_cff::dEta, doTree_, HLT_FULL_cff::dPhi, HGC3DClusterGenMatchSelector_cfi::dR, flavorHistoryFilter_cfi::dr, spr::eECALmatrix(), spr::eHCALmatrix(), reco::TrackBase::eta(), HLT_FULL_cff::eta1, fillEnergy(), fillIsolation(), fillTrack(), firstEvent_, genParticles2HepMC_cfi::genEventInfo, genParticles2HepMC_cfi::genParticles, get, edm::EventSetup::getData(), h_counter, h_goodPV, h_goodRun, h_h_pNew, h_HLTAccept, h_HLTAccepts, h_HLTCorr, h_nHLT, h_nHLTvsRN, h_ntrk, h_numberPV, photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, HLTNames_, mps_fire::i, LEDCalibrationChannels::ieta, iEvent, createfilelist::int, edm::HandleBase::isValid(), edm::InputTag::label(), mix_E10TeV_L13E31_BX432_cfi::Lumi, M_PI, match(), SiStripPI::max, maxTrackEta_, spr::trackSelectionParameters::minQuality, minTrackP_, reco::TrackBase::momentum(), newNames_, nGen_, nPVBin_, nRun_, convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), p1, reco::TrackBase::phi(), edm::Handle< T >::product(), spr::propagateCALO(), reco::TrackBase::pt(), HLT_FULL_cff::pt1, puWeights_, pvBin_, quality, reco::TrackBase::quality(), muon::segmentCompatibility(), selectionParameters_, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality_, tMaxE_, tMaxH_, tMinE_, tMinH_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_ew_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_magField_, tok_Muon_, tok_parts_, tok_recVtx_, tok_sevlv_, tok_topo_, 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_, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEvent_, L1TEGammaOffline_cfi::triggerNames, triggerResults, trigNames_, truncate_str(), verbosity_, vetoEcal_, vetoMuon_, and vetoTrigger_.

◆ beginJob()

void StudyCaloResponse::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 779 of file StudyCaloResponse.cc.

                                  {
   // 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] - 1),
                     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, 5000, -0.1, 49.9);
           h_energy[i][ip][ie][j]->Sumw2();
         }
       }
     }
   }
  
   for (int i = 0; i < 8; ++i) {
     sprintf(hname, "counter%d", i);
     sprintf(htit, "Counter with cut %d", i);
     h_counter[i] = fs_->make<TH1D>(hname, htit, 1000, 0, 1000);
     sprintf(hname, "h_pTNew%d", i);
     sprintf(htit, "Track momentum with cut %d", i);
     h_h_pNew[i] = fs_->make<TH1D>(hname, htit, 1000, 0, 1000);
   }
  
   // 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);
   }
 }

References doTree_, etaBin_, fs_, h_counter, h_ediff, h_ene1, h_ene2, h_energy, h_eta, h_goodPV, h_goodRun, h_h_pNew, h_HLTAccept, h_HLTCorr, h_maxNearP, h_nHLT, h_nHLTvsRN, h_ntrk, h_numberPV, h_p, h_phi, h_pt, mps_fire::i, dqmiolumiharvest::j, 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_.

◆ beginLuminosityBlock()

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

inlineprivatevirtual

Definition at line 97 of file StudyCaloResponse.cc.

97 {}

◆ beginRun()

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

overrideprivate

Definition at line 969 of file StudyCaloResponse.cc.

                                                                                 {
   char hname[100], htit[400];
   edm::LogVerbatim("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::LogVerbatim("IsoTrack") << "beginRun " << iRun.run();
   firstEvent_ = true;
   changed_ = false;
 }

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

◆ clear()

void StudyCaloResponse::clear ( void )

private

Definition at line 990 of file StudyCaloResponse.cc.

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

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(), and BeautifulSoup.Tag::setString().

◆ endLuminosityBlock()

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

inlineprivatevirtual

Definition at line 98 of file StudyCaloResponse.cc.

98 {}

◆ endRun()

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

overrideprivate

Definition at line 985 of file StudyCaloResponse.cc.

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

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

◆ fillDescriptions()

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

static

Definition at line 264 of file StudyCaloResponse.cc.

                                                                                    {
   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<edm::InputTag>("labelMuon", edm::InputTag("muons", "", "RECO"));
   desc.addUntracked<edm::InputTag>("labelTrack", edm::InputTag("generalTracks", "", "RECO"));
   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>("vetoTrigger", false);
   desc.addUntracked<bool>("doTree", false);
   desc.addUntracked<bool>("vetoMuon", true);
   desc.addUntracked<double>("cutMuon", 0.1);
   desc.addUntracked<bool>("vetoEcal", false);
   desc.addUntracked<double>("cutEcal", 2.0);
   desc.addUntracked<double>("cutRatio", 0.9);
   desc.addUntracked<std::vector<double> >("puWeights", weights);
   descriptions.add("studyCaloResponse", desc);
 }

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, ClusterTask_cfi::trig, and HLT_FULL_cff::weights.

◆ fillEnergy()

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

private

Definition at line 1025 of file StudyCaloResponse.cc.

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

References etaBin_, RemoveAddSevLevel::flag, h_energy, mps_fire::i, LEDCalibrationChannels::ieta, nEtaBin_, nPBin_, AlCaHLTBitMon_ParallelJobs::p, pBin_, and tr_eventWeight.

Referenced by analyze().

◆ fillIsolation()

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

private

Definition at line 1018 of file StudyCaloResponse.cc.

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

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

Referenced by analyze().

◆ fillTrack()

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

private

Definition at line 1011 of file StudyCaloResponse.cc.

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

References PVValHelper::eta, h_eta, h_p, h_phi, h_pt, mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, phi, DiDispStaMuonMonitor_cfi::pt, and tr_eventWeight.

Referenced by analyze().

◆ trackPID()

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

private

Definition at line 1062 of file StudyCaloResponse.cc.

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

References funct::abs(), HLT_FULL_cff::dEta, HLT_FULL_cff::dPhi, HGC3DClusterGenMatchSelector_cfi::dR, reco::TrackBase::eta(), genParticles2HepMC_cfi::genParticles, triggerObjects_cff::id, heavyIonCSV_trainingSettings::idx, M_PI, AlCaHLTBitMon_ParallelJobs::p, EgammaValidation_cff::pdgId, reco::TrackBase::phi(), and mathSSE::sqrt().

Referenced by analyze().

◆ truncate_str()

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

private

Definition at line 1050 of file StudyCaloResponse.cc.

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

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

Referenced by analyze().

Member Data Documentation

◆ changed_

bool StudyCaloResponse::changed_

private

Definition at line 124 of file StudyCaloResponse.cc.

Referenced by beginRun(), and StudyCaloResponse().

◆ cutEcal_

const double StudyCaloResponse::cutEcal_

private

Definition at line 119 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ cutMuon_

const double StudyCaloResponse::cutMuon_

private

Definition at line 119 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ cutRatio_

const double StudyCaloResponse::cutRatio_

private

Definition at line 119 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ doTree_

const bool StudyCaloResponse::doTree_

private

Definition at line 118 of file StudyCaloResponse.cc.

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

◆ etaBin_

int StudyCaloResponse::etaBin_[nEtaBin_+1]

private

Definition at line 155 of file StudyCaloResponse.cc.

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

◆ firstEvent_

bool StudyCaloResponse::firstEvent_

private

Definition at line 124 of file StudyCaloResponse.cc.

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

◆ fs_

edm::Service<TFileService> StudyCaloResponse::fs_

private

Definition at line 111 of file StudyCaloResponse.cc.

Referenced by beginJob(), and beginRun().

◆ h_counter

TH1D * StudyCaloResponse::h_counter[8]

private

Definition at line 149 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_ediff

TH1D * StudyCaloResponse::h_ediff[2]

private

Definition at line 152 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillIsolation().

◆ h_ene1

TH1D * StudyCaloResponse::h_ene1[2]

private

Definition at line 152 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillIsolation().

◆ h_ene2

TH1D * StudyCaloResponse::h_ene2[2]

private

Definition at line 152 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillIsolation().

◆ h_energy

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

private

Definition at line 153 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillEnergy().

◆ h_eta

TH1D* StudyCaloResponse::h_eta[nGen_+2]

private

Definition at line 150 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillTrack().

◆ h_goodPV

TH1I* StudyCaloResponse::h_goodPV

private

Definition at line 146 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_goodRun

TH1I * StudyCaloResponse::h_goodRun

private

Definition at line 146 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_h_pNew

TH1D * StudyCaloResponse::h_h_pNew[8]

private

Definition at line 150 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_HLTAccept

TH1I * StudyCaloResponse::h_HLTAccept

private

Definition at line 145 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_HLTAccepts

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

private

Definition at line 148 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginRun().

◆ h_HLTCorr

TH1I * StudyCaloResponse::h_HLTCorr

private

Definition at line 145 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_maxNearP

TH1D* StudyCaloResponse::h_maxNearP[2]

private

Definition at line 152 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillIsolation().

◆ h_nHLT

TH1I* StudyCaloResponse::h_nHLT

private

Definition at line 145 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_nHLTvsRN

TH2I* StudyCaloResponse::h_nHLTvsRN

private

Definition at line 147 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_ntrk

TH1I* StudyCaloResponse::h_ntrk[2]

private

Definition at line 151 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_numberPV

TH1I * StudyCaloResponse::h_numberPV

private

Definition at line 145 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ h_p

TH1D* StudyCaloResponse::h_p[nGen_+2]

private

Definition at line 149 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillTrack().

◆ h_phi

TH1D * StudyCaloResponse::h_phi[nGen_+2]

private

Definition at line 150 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillTrack().

◆ h_pt

TH1D * StudyCaloResponse::h_pt[nGen_+2]

private

Definition at line 149 of file StudyCaloResponse.cc.

Referenced by beginJob(), and fillTrack().

◆ hltConfig_

HLTConfigProvider StudyCaloResponse::hltConfig_

private

Definition at line 110 of file StudyCaloResponse.cc.

Referenced by beginRun().

◆ HLTNames_

std::vector<std::string> StudyCaloResponse::HLTNames_

private

Definition at line 123 of file StudyCaloResponse.cc.

Referenced by analyze().

◆ isItAOD_

const bool StudyCaloResponse::isItAOD_

private

Definition at line 118 of file StudyCaloResponse.cc.

Referenced by StudyCaloResponse().

◆ labelGenTrack_

const edm::InputTag StudyCaloResponse::labelGenTrack_

private

Definition at line 114 of file StudyCaloResponse.cc.

Referenced by StudyCaloResponse().

◆ labelMuon_

const edm::InputTag StudyCaloResponse::labelMuon_

private

Definition at line 114 of file StudyCaloResponse.cc.

Referenced by StudyCaloResponse().

◆ maxTrackEta_

const double StudyCaloResponse::maxTrackEta_

private

Definition at line 116 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ minTrackP_

const double StudyCaloResponse::minTrackP_

private

Definition at line 116 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ nEtaBin_

const int StudyCaloResponse::nEtaBin_ = 4

staticprivate

Definition at line 108 of file StudyCaloResponse.cc.

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

◆ newNames_

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

private

Definition at line 113 of file StudyCaloResponse.cc.

Referenced by analyze().

◆ nGen_

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

staticprivate

Definition at line 109 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ nPBin_

const int StudyCaloResponse::nPBin_ = 15

staticprivate

Definition at line 108 of file StudyCaloResponse.cc.

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

◆ nPVBin_

const int StudyCaloResponse::nPVBin_ = 4

staticprivate

Definition at line 108 of file StudyCaloResponse.cc.

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

◆ nRun_

int StudyCaloResponse::nRun_

private

Definition at line 155 of file StudyCaloResponse.cc.

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

◆ pBin_

double StudyCaloResponse::pBin_[nPBin_+1]

private

Definition at line 156 of file StudyCaloResponse.cc.

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

◆ puWeights_

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

private

Definition at line 120 of file StudyCaloResponse.cc.

Referenced by analyze().

◆ pvBin_

int StudyCaloResponse::pvBin_[nPVBin_+1]

private

Definition at line 155 of file StudyCaloResponse.cc.

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

◆ selectionParameters_

spr::trackSelectionParameters StudyCaloResponse::selectionParameters_

private

Definition at line 122 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ theTrackQuality_

const std::string StudyCaloResponse::theTrackQuality_

private

Definition at line 115 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ theTriggerResultsLabel_

const edm::InputTag StudyCaloResponse::theTriggerResultsLabel_

private

Definition at line 121 of file StudyCaloResponse.cc.

Referenced by StudyCaloResponse().

◆ tMaxE_

const double StudyCaloResponse::tMaxE_

private

Definition at line 117 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tMaxH_

const double StudyCaloResponse::tMaxH_

private

Definition at line 117 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tMinE_

const double StudyCaloResponse::tMinE_

private

Definition at line 117 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tMinH_

const double StudyCaloResponse::tMinH_

private

Definition at line 117 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_caloTopology_

edm::ESGetToken<CaloTopology, CaloTopologyRecord> StudyCaloResponse::tok_caloTopology_

private

Definition at line 139 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> StudyCaloResponse::tok_EB_

private

Definition at line 133 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_ecalChStatus_

edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> StudyCaloResponse::tok_ecalChStatus_

private

Definition at line 142 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> StudyCaloResponse::tok_EE_

private

Definition at line 134 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_ew_

edm::EDGetTokenT<GenEventInfoProduct> StudyCaloResponse::tok_ew_

private

Definition at line 136 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_genTrack_

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

private

Definition at line 131 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_geom_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> StudyCaloResponse::tok_geom_

private

Definition at line 138 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> StudyCaloResponse::tok_hbhe_

private

Definition at line 135 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_lumi

edm::EDGetTokenT<LumiDetails> StudyCaloResponse::tok_lumi

private

Definition at line 126 of file StudyCaloResponse.cc.

Referenced by StudyCaloResponse().

◆ tok_magField_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> StudyCaloResponse::tok_magField_

private

Definition at line 141 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_Muon_

edm::EDGetTokenT<reco::MuonCollection> StudyCaloResponse::tok_Muon_

private

Definition at line 130 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_parts_

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

private

Definition at line 129 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_recVtx_

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

private

Definition at line 132 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_sevlv_

edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> StudyCaloResponse::tok_sevlv_

private

Definition at line 143 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_topo_

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> StudyCaloResponse::tok_topo_

private

Definition at line 140 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_trigEvt

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

private

Definition at line 127 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tok_trigRes

edm::EDGetTokenT<edm::TriggerResults> StudyCaloResponse::tok_trigRes

private

Definition at line 128 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ tr_eventWeight

double StudyCaloResponse::tr_eventWeight

private

Definition at line 158 of file StudyCaloResponse.cc.

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

◆ tr_FE11x11P

std::vector<double> StudyCaloResponse::tr_FE11x11P

private

Definition at line 163 of file StudyCaloResponse.cc.

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

◆ tr_FE15x15P

std::vector<double> StudyCaloResponse::tr_FE15x15P

private

Definition at line 163 of file StudyCaloResponse.cc.

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

◆ tr_FE7x7P

std::vector<double> StudyCaloResponse::tr_FE7x7P

private

Definition at line 163 of file StudyCaloResponse.cc.

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

◆ tr_goodPV

int StudyCaloResponse::tr_goodPV

private

Definition at line 157 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ tr_goodRun

int StudyCaloResponse::tr_goodRun

private

Definition at line 157 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ tr_H3x3

std::vector<double> StudyCaloResponse::tr_H3x3

private

Definition at line 162 of file StudyCaloResponse.cc.

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

◆ tr_H5x5

std::vector<double> StudyCaloResponse::tr_H5x5

private

Definition at line 162 of file StudyCaloResponse.cc.

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

◆ tr_H7x7

std::vector<double> StudyCaloResponse::tr_H7x7

private

Definition at line 162 of file StudyCaloResponse.cc.

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

◆ tr_iEta

std::vector<int> StudyCaloResponse::tr_iEta

private

Definition at line 165 of file StudyCaloResponse.cc.

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

◆ tr_MaxNearHcalP7x7

std::vector<double> StudyCaloResponse::tr_MaxNearHcalP7x7

private

Definition at line 161 of file StudyCaloResponse.cc.

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

◆ tr_MaxNearP31X31

std::vector<double> StudyCaloResponse::tr_MaxNearP31X31

private

Definition at line 161 of file StudyCaloResponse.cc.

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

◆ tr_SE11x11P

std::vector<bool> StudyCaloResponse::tr_SE11x11P

private

Definition at line 164 of file StudyCaloResponse.cc.

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

◆ tr_SE15x15P

std::vector<bool> StudyCaloResponse::tr_SE15x15P

private

Definition at line 164 of file StudyCaloResponse.cc.

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

◆ tr_SE7x7P

std::vector<bool> StudyCaloResponse::tr_SE7x7P

private

Definition at line 164 of file StudyCaloResponse.cc.

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

◆ tr_TrigName

std::vector<std::string> StudyCaloResponse::tr_TrigName

private

Definition at line 159 of file StudyCaloResponse.cc.

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

◆ tr_TrkEta

std::vector<double> StudyCaloResponse::tr_TrkEta

private

Definition at line 160 of file StudyCaloResponse.cc.

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

◆ tr_TrkID

std::vector<int> StudyCaloResponse::tr_TrkID

private

Definition at line 165 of file StudyCaloResponse.cc.

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

◆ tr_TrkP

std::vector<double> StudyCaloResponse::tr_TrkP

private

Definition at line 160 of file StudyCaloResponse.cc.

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

◆ tr_TrkPhi

std::vector<double> StudyCaloResponse::tr_TrkPhi

private

Definition at line 160 of file StudyCaloResponse.cc.

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

◆ tr_TrkPt

std::vector<double> StudyCaloResponse::tr_TrkPt

private

Definition at line 160 of file StudyCaloResponse.cc.

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

◆ tree_

TTree* StudyCaloResponse::tree_

private

Definition at line 154 of file StudyCaloResponse.cc.

Referenced by analyze(), and beginJob().

◆ triggerEvent_

const edm::InputTag StudyCaloResponse::triggerEvent_

private

Definition at line 121 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ trigNames_

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

private

Definition at line 113 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ verbosity_

const int StudyCaloResponse::verbosity_

private

Definition at line 112 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ vetoEcal_

const bool StudyCaloResponse::vetoEcal_

private

Definition at line 118 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ vetoMuon_

const bool StudyCaloResponse::vetoMuon_

private

Definition at line 118 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

◆ vetoTrigger_

const bool StudyCaloResponse::vetoTrigger_

private

Definition at line 118 of file StudyCaloResponse.cc.

Referenced by analyze(), and StudyCaloResponse().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ StudyCaloResponse()

◆ ~StudyCaloResponse()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginLuminosityBlock()

◆ beginRun()

◆ clear()

◆ endLuminosityBlock()

◆ endRun()

◆ fillDescriptions()

◆ fillEnergy()

◆ fillIsolation()

◆ fillTrack()

◆ trackPID()

◆ truncate_str()

Member Data Documentation

◆ changed_

◆ cutEcal_

◆ cutMuon_

◆ cutRatio_

◆ doTree_

◆ etaBin_

◆ firstEvent_

◆ fs_

◆ h_counter

◆ h_ediff

◆ h_ene1

◆ h_ene2

◆ h_energy

◆ h_eta

◆ h_goodPV

◆ h_goodRun

◆ h_h_pNew

◆ h_HLTAccept

◆ h_HLTAccepts

◆ h_HLTCorr

◆ h_maxNearP

◆ h_nHLT

◆ h_nHLTvsRN

◆ h_ntrk

◆ h_numberPV

◆ h_p

◆ h_phi

◆ h_pt

◆ hltConfig_

◆ HLTNames_

◆ isItAOD_

◆ labelGenTrack_

◆ labelMuon_

◆ maxTrackEta_

◆ minTrackP_

◆ nEtaBin_

◆ newNames_

◆ nGen_

◆ nPBin_

◆ nPVBin_

◆ nRun_

◆ pBin_

◆ puWeights_

◆ pvBin_

◆ selectionParameters_

◆ theTrackQuality_

◆ theTriggerResultsLabel_

◆ tMaxE_

◆ tMaxH_

◆ tMinE_

◆ tMinH_

◆ tok_caloTopology_

◆ tok_EB_

◆ tok_ecalChStatus_

◆ tok_EE_