#include <IsoTrig/IsoTrig/src/IsoTrig.cc>

Inheritance diagram for IsoTrig:

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

	~IsoTrig () 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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > 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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, 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::ESRecordsToProductResolverIndices 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 (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

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

void	chgIsolation (double &etaTriggered, double &phiTriggered, edm::Handle< reco::TrackCollection > &trkCollection, const edm::Event &theEvent)

void	clearChgIsolnTreeVectors ()

void	clearMipCutTreeVectors ()

double	dEta (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

double	dinvPt (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

double	dP (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

double	dPhi (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

double	dPt (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

double	dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

void	endJob () override

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

std::pair< double, double >	etaPhiTrigger ()

void	fillCuts (int, double, double, double, math::XYZTLorentzVector &, int, bool)

void	fillDifferences (int, math::XYZTLorentzVector &, math::XYZTLorentzVector &, bool)

void	fillEnergy (int, int, double, double, math::XYZTLorentzVector &)

void	fillHist (int, math::XYZTLorentzVector &)

double	getDistInCM (double eta1, double phi1, double eta2, double phi2)

std::pair< double, double >	GetEtaPhiAtEcal (double etaIP, double phiIP, double pT, int charge, double vtxZ)

void	getGoodTracks (const edm::Event &, edm::Handle< reco::TrackCollection > &)

void	pushChgIsolnTreeVecs (math::XYZTLorentzVector &Pixcand, math::XYZTLorentzVector &Trkcand, std::vector< double > &PixMaxP, double &TrkMaxP, bool &selTk)

void	pushMipCutTreeVecs (math::XYZTLorentzVector &NFcand, math::XYZTLorentzVector &Trkcand, double &EmipNFcand, double &EmipTrkcand, double &mindR, double &mindP1, std::vector< bool > &Flags, double hCone)

void	studyIsolation (edm::Handle< reco::TrackCollection > &, std::vector< reco::TrackCollection::const_iterator > &)

void	studyMipCut (edm::Handle< reco::TrackCollection > &trkCollection, edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &L2cands)

void	studyTiming (const edm::Event &theEvent)

void	studyTrigger (edm::Handle< reco::TrackCollection > &, std::vector< reco::TrackCollection::const_iterator > &)

void	StudyTrkEbyP (edm::Handle< reco::TrackCollection > &trkCollection)

Private Attributes
const double	a_charIsoR_

const double	a_coneR_

const double	a_mipR_

const double	a_neutIsoR_

const double	a_neutR1_

const double	a_neutR2_

edm::Handle< EcalRecHitCollection >	barrelRecHitsHandle_

edm::Handle< reco::BeamSpot >	beamSpotH_

const MagneticField *	bField_

double	bfVal_

bool	changed_

TTree *	ChgIsolnTree_

const double	cutCharge_

const double	cutMip_

const double	cutNeutral_

const bool	doChgIsolTree_

const bool	doL2L3_

const bool	doMipCutTree_

const bool	doStudyIsol_

const bool	doTiming_

const bool	doTrkResTree_

const double	dr_L1_

edm::Handle< EcalRecHitCollection >	endcapRecHitsHandle_

edm::Service< TFileService >	fs_

TH1I *	g_Accepts

TH1D *	g_Pre

TH1D *	g_PreHLT

TH1D *	g_PreL1

const CaloGeometry *	geo_

TH1D *	h_dEta [9]

TH1D *	h_dEtaL1 [2]

TH1D *	h_dinvPt [9]

TH1D *	h_dP [9]

TH1D *	h_dPhi [9]

TH1D *	h_dPhiL1 [2]

TH1D *	h_dPt [9]

TH1D *	h_dRL1 [2]

TH1D *	h_eCalo [5][6][48]

TH1D *	h_eHcal [5][6][48]

TH1D *	h_eMaxNearP [2]

TH1D *	h_eMip [2]

TH1D *	h_eNeutIso [2]

TH1D *	h_EnIn

TH1D *	h_EnOut

TH1D *	h_eta [20]

TH1D *	h_etaCalibTracks [5][2][2]

TH1D *	h_etaMipTracks [5][2][2]

TH1I *	h_Filters

TH1I *	h_HLT

TH1D *	h_L1ObjEnergy

TH1D *	h_mindR [9]

TH2D *	h_MipEnMatch

TH2D *	h_MipEnNoMatch

TH1I *	h_nHLT

TH1I *	h_nL3Objs

TH1D *	h_p [20]

TH1D *	h_phi [20]

TH1I *	h_Pre

TH1I *	h_PreHLT

TH1D *	h_PreHLTwt

TH1I *	h_PreL1

TH1D *	h_PreL1wt

TH1D *	h_pt [20]

edm::Handle< HBHERecHitCollection >	hbhe_

HLTPrescaleProvider	hltPrescaleProvider_

const edm::InputTag	l1CandTag_

const edm::InputTag	l2CandTag_

math::XYZPoint	leadPV_

const int	maxRunNo_

const double	minPTrackValue_

const int	minRunNo_

TTree *	MipCutTree_

const edm::InputTag	pixCandTag_

const std::vector< double >	pixelIsolationConeSizeAtEC_

std::vector< reco::TrackRef >	pixelTrackRefsHB_

std::vector< reco::TrackRef >	pixelTrackRefsHE_

const std::vector< edm::InputTag >	pixelTracksSources_

double	pLimits_ [6]

const double	prelimCone_

const std::string	processName_

double	rEB_

edm::Handle< reco::VertexCollection >	recVtxs_

spr::trackSelectionParameters	selectionParameters_

std::vector< double > *	t_ECone

std::vector< double > *	t_NFcandEmip

std::vector< double > *	t_NFcandEta

std::vector< double > *	t_NFcandP

std::vector< double > *	t_NFcandPhi

std::vector< double > *	t_NFcandPt

std::vector< double > *	t_NFTrkcandEmip

std::vector< double > *	t_NFTrkcandEta

std::vector< double > *	t_NFTrkcandP

std::vector< double > *	t_NFTrkcandPhi

std::vector< double > *	t_NFTrkcandPt

std::vector< bool > *	t_NFTrkChgIsoFlag

std::vector< double > *	t_NFTrkMinDP1

std::vector< double > *	t_NFTrkMinDR

std::vector< bool > *	t_NFTrkMipFlag

std::vector< bool > *	t_NFTrkMissFlag

std::vector< bool > *	t_NFTrkNeuIsoFlag

std::vector< bool > *	t_NFTrkPropFlag

std::vector< bool > *	t_NFTrkPVFlag

std::vector< bool > *	t_NFTrkqltyFlag

std::vector< bool > *	t_NFTrkselTkFlag

std::vector< int > *	t_nGoodTk

std::vector< int > *	t_nPixCand

std::vector< int > *	t_nPixSeed

std::vector< double > *	t_PixcandEta

std::vector< std::vector< double > > *	t_PixcandMaxP

std::vector< double > *	t_PixcandP

std::vector< double > *	t_PixcandPhi

std::vector< double > *	t_PixcandPt

std::vector< double > *	t_PixTrkcandEta

std::vector< double > *	t_PixTrkcandMaxP

std::vector< double > *	t_PixTrkcandP

std::vector< double > *	t_PixTrkcandPhi

std::vector< double > *	t_PixTrkcandPt

std::vector< bool > *	t_PixTrkcandselTk

std::vector< double > *	t_timeL2Prod

std::vector< double > *	t_TrkhCone

std::vector< bool > *	t_TrkMissFlag

std::vector< bool > *	t_TrkNuIsolFlag

std::vector< double > *	t_TrkP

std::vector< bool > *	t_TrkPVFlag

std::vector< bool > *	t_TrkqltyFlag

std::vector< bool > *	t_TrkselTkFlag

const double	tauUnbiasCone_

std::string	theTrackQuality_

TTree *	TimingTree_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	tok_hlt_

edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	tok_l1cand_

edm::EDGetTokenT< reco::IsolatedPixelTrackCandidateCollection >	tok_l2cand_

edm::EDGetTokenT< LumiDetails >	tok_lumi_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_magField_

edm::EDGetTokenT< reco::IsolatedPixelTrackCandidateCollection >	tok_pixtk_

std::vector< edm::EDGetTokenT< reco::TrackCollection > >	tok_pixtks_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT< SeedingLayerSetsHits >	tok_SeedingLayerHB_

edm::EDGetTokenT< SeedingLayerSetsHits >	tok_SeedingLayerHE_

edm::EDGetTokenT< SiPixelRecHitCollection >	tok_SiPixelRecHits_

edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_

edm::EDGetTokenT< reco::VertexCollection >	tok_verthb_

edm::EDGetTokenT< reco::VertexCollection >	tok_verthe_

std::unordered_map< unsigned int, unsigned int >	trigList_

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

std::unordered_map< unsigned int, const std::pair< double, double > >	trigPreList_

TTree *	TrkResTree_

std::vector< math::XYZTLorentzVector >	vec_ [3]

const int	verbosity_

const double	vtxCutIsol_

const double	vtxCutSeed_

double	zEE_

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)

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

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

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

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

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

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

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 96 of file IsoTrig.cc.

Constructor & Destructor Documentation

◆ IsoTrig()

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

explicit

Definition at line 272 of file IsoTrig.cc.

     : hltPrescaleProvider_(iConfig, consumesCollector(), *this),
       trigNames_(iConfig.getUntrackedParameter<std::vector<std::string>>("Triggers")),
       pixCandTag_(iConfig.getUntrackedParameter<edm::InputTag>("pixCandTag")),
       l1CandTag_(iConfig.getUntrackedParameter<edm::InputTag>("l1CandTag")),
       l2CandTag_(iConfig.getUntrackedParameter<edm::InputTag>("l2CandTag")),
       pixelTracksSources_(iConfig.getUntrackedParameter<std::vector<edm::InputTag>>("pixelTracksSources")),
       doL2L3_(iConfig.getUntrackedParameter<bool>("doL2L3", true)),
       doTiming_(iConfig.getUntrackedParameter<bool>("doTimingTree", true)),
       doMipCutTree_(iConfig.getUntrackedParameter<bool>("doMipCutTree", true)),
       doTrkResTree_(iConfig.getUntrackedParameter<bool>("doTrkResTree", true)),
       doChgIsolTree_(iConfig.getUntrackedParameter<bool>("doChgIsolTree", true)),
       doStudyIsol_(iConfig.getUntrackedParameter<bool>("doStudyIsol", true)),
       verbosity_(iConfig.getUntrackedParameter<int>("verbosity", 0)),
       pixelIsolationConeSizeAtEC_(iConfig.getUntrackedParameter<std::vector<double>>("pixelIsolationConeSizeAtEC")),
       minPTrackValue_(iConfig.getUntrackedParameter<double>("minPTrackValue")),
       vtxCutSeed_(iConfig.getUntrackedParameter<double>("vertexCutSeed")),
       vtxCutIsol_(iConfig.getUntrackedParameter<double>("vertexCutIsol")),
       tauUnbiasCone_(iConfig.getUntrackedParameter<double>("tauUnbiasCone")),
       prelimCone_(iConfig.getUntrackedParameter<double>("prelimCone")),
       theTrackQuality_(iConfig.getUntrackedParameter<std::string>("trackQuality", "highPurity")),
       processName_(iConfig.getUntrackedParameter<std::string>("processName", "HLT")),
       dr_L1_(iConfig.getUntrackedParameter<double>("isolationL1", 1.0)),
       a_coneR_(iConfig.getUntrackedParameter<double>("coneRadius", 34.98)),
       a_charIsoR_(a_coneR_ + 28.9),
       a_neutIsoR_(a_charIsoR_ * 0.726),
       a_mipR_(iConfig.getUntrackedParameter<double>("coneRadiusMIP", 14.0)),
       a_neutR1_(iConfig.getUntrackedParameter<double>("coneRadiusNeut1", 21.0)),
       a_neutR2_(iConfig.getUntrackedParameter<double>("coneRadiusNeut2", 29.0)),
       cutMip_(iConfig.getUntrackedParameter<double>("cutMIP", 1.0)),
       cutCharge_(iConfig.getUntrackedParameter<double>("chargeIsolation", 2.0)),
       cutNeutral_(iConfig.getUntrackedParameter<double>("neutralIsolation", 2.0)),
       minRunNo_(iConfig.getUntrackedParameter<int>("minRun")),
       maxRunNo_(iConfig.getUntrackedParameter<int>("maxRun")),
       changed_(false),
       t_timeL2Prod(nullptr),
       t_nPixCand(nullptr),
       t_nPixSeed(nullptr),
       t_nGoodTk(nullptr),
       t_TrkhCone(nullptr),
       t_TrkP(nullptr),
       t_TrkselTkFlag(nullptr),
       t_TrkqltyFlag(nullptr),
       t_TrkMissFlag(nullptr),
       t_TrkPVFlag(nullptr),
       t_TrkNuIsolFlag(nullptr),
       t_PixcandP(nullptr),
       t_PixcandPt(nullptr),
       t_PixcandEta(nullptr),
       t_PixcandPhi(nullptr),
       t_PixcandMaxP(nullptr),
       t_PixTrkcandP(nullptr),
       t_PixTrkcandPt(nullptr),
       t_PixTrkcandEta(nullptr),
       t_PixTrkcandPhi(nullptr),
       t_PixTrkcandMaxP(nullptr),
       t_PixTrkcandselTk(nullptr),
       t_NFcandP(nullptr),
       t_NFcandPt(nullptr),
       t_NFcandEta(nullptr),
       t_NFcandPhi(nullptr),
       t_NFcandEmip(nullptr),
       t_NFTrkcandP(nullptr),
       t_NFTrkcandPt(nullptr),
       t_NFTrkcandEta(nullptr),
       t_NFTrkcandPhi(nullptr),
       t_NFTrkcandEmip(nullptr),
       t_NFTrkMinDR(nullptr),
       t_NFTrkMinDP1(nullptr),
       t_NFTrkselTkFlag(nullptr),
       t_NFTrkqltyFlag(nullptr),
       t_NFTrkMissFlag(nullptr),
       t_NFTrkPVFlag(nullptr),
       t_NFTrkPropFlag(nullptr),
       t_NFTrkChgIsoFlag(nullptr),
       t_NFTrkNeuIsoFlag(nullptr),
       t_NFTrkMipFlag(nullptr),
       t_ECone(nullptr) {
   usesResource(TFileService::kSharedResource);
 
   //now do whatever initialization is needed
   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"));
   edm::InputTag triggerEvent_("hltTriggerSummaryAOD", "", processName_);
   tok_trigEvt_ = consumes<trigger::TriggerEvent>(triggerEvent_);
   edm::InputTag theTriggerResultsLabel("TriggerResults", "", processName_);
   tok_trigRes_ = consumes<edm::TriggerResults>(theTriggerResultsLabel);
   tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
   tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
   tok_bs_ = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
   tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
   tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
   tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
   tok_pixtk_ = consumes<reco::IsolatedPixelTrackCandidateCollection>(pixCandTag_);
   tok_l1cand_ = consumes<trigger::TriggerFilterObjectWithRefs>(l1CandTag_);
   tok_l2cand_ = consumes<reco::IsolatedPixelTrackCandidateCollection>(l2CandTag_);
   if (doTiming_) {
     tok_verthb_ = consumes<reco::VertexCollection>(edm::InputTag("hltHITPixelVerticesHB"));
     tok_verthe_ = consumes<reco::VertexCollection>(edm::InputTag("hltHITPixelVerticesHB"));
     tok_hlt_ = consumes<trigger::TriggerFilterObjectWithRefs>(edm::InputTag("hltL1sL1SingleJet68"));
     tok_SeedingLayerHB_ = consumes<SeedingLayerSetsHits>(edm::InputTag("hltPixelLayerTripletsHITHB"));
     tok_SeedingLayerHE_ = consumes<SeedingLayerSetsHits>(edm::InputTag("hltPixelLayerTripletsHITHE"));
     tok_SiPixelRecHits_ = consumes<SiPixelRecHitCollection>(edm::InputTag("hltSiPixelRecHits"));
   }
   if (doChgIsolTree_) {
     for (unsigned int k = 0; k < pixelTracksSources_.size(); ++k) {
       //      edm::InputTag  pix (pixelTracksSources_[k],"",processName_);
       //      tok_pixtks_.push_back(consumes<reco::TrackCollection>(pix));
       tok_pixtks_.push_back(consumes<reco::TrackCollection>(pixelTracksSources_[k]));
     }
   }
   if (verbosity_ >= 0) {
     edm::LogVerbatim("IsoTrack") << "Parameters read from config file \n"
                                  << "\t minPt " << selectionParameters_.minPt << "\t theTrackQuality "
                                  << theTrackQuality_ << "\t minQuality " << selectionParameters_.minQuality
                                  << "\t maxDxyPV " << selectionParameters_.maxDxyPV << "\t maxDzPV "
                                  << selectionParameters_.maxDzPV << "\t maxChi2 " << selectionParameters_.maxChi2
                                  << "\t maxDpOverP " << selectionParameters_.maxDpOverP << "\t minOuterHit "
                                  << selectionParameters_.minOuterHit << "\t minLayerCrossed "
                                  << selectionParameters_.minLayerCrossed << "\t maxInMiss "
                                  << selectionParameters_.maxInMiss << "\t maxOutMiss "
                                  << selectionParameters_.maxOutMiss << "\t a_coneR " << a_coneR_ << "\t a_charIsoR "
                                  << a_charIsoR_ << "\t a_neutIsoR " << a_neutIsoR_ << "\t a_mipR " << a_mipR_
                                  << "\t a_neutR " << a_neutR1_ << ":" << a_neutR2_ << "\t cuts (MIP " << cutMip_
                                  << " : Charge " << cutCharge_ << " : Neutral " << cutNeutral_ << ")";
     edm::LogVerbatim("IsoTrack") << "Charge Isolation parameters:"
                                  << "\t minPTrackValue " << minPTrackValue_ << "\t vtxCutSeed " << vtxCutSeed_
                                  << "\t vtxCutIsol " << vtxCutIsol_ << "\t tauUnbiasCone " << tauUnbiasCone_
                                  << "\t prelimCone " << prelimCone_ << "\t pixelIsolationConeSizeAtEC";
     for (unsigned int k = 0; k < pixelIsolationConeSizeAtEC_.size(); ++k)
       edm::LogVerbatim("IsoTrack") << "[" << k << "] " << pixelIsolationConeSizeAtEC_[k];
   }
   double pl[] = {20, 30, 40, 60, 80, 120};
   for (int i = 0; i < 6; ++i)
     pLimits_[i] = pl[i];
   rEB_ = 123.8;
   zEE_ = 317.0;
 
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
   tok_magField_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
 }

◆ ~IsoTrig()

IsoTrig::~IsoTrig ( )

override

Definition at line 426 of file IsoTrig.cc.

References t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, t_NFTrkselTkFlag, t_nGoodTk, t_nPixCand, t_nPixSeed, t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, t_PixTrkcandselTk, t_timeL2Prod, t_TrkhCone, t_TrkMissFlag, t_TrkNuIsolFlag, t_TrkP, t_TrkPVFlag, t_TrkqltyFlag, and t_TrkselTkFlag.

                   {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   if (t_timeL2Prod)
     delete t_timeL2Prod;
   if (t_nPixCand)
     delete t_nPixCand;
   if (t_nPixSeed)
     delete t_nPixSeed;
   if (t_nGoodTk)
     delete t_nGoodTk;
   if (t_TrkhCone)
     delete t_TrkhCone;
   if (t_TrkP)
     delete t_TrkP;
   if (t_TrkselTkFlag)
     delete t_TrkselTkFlag;
   if (t_TrkqltyFlag)
     delete t_TrkqltyFlag;
   if (t_TrkMissFlag)
     delete t_TrkMissFlag;
   if (t_TrkPVFlag)
     delete t_TrkPVFlag;
   if (t_TrkNuIsolFlag)
     delete t_TrkNuIsolFlag;
   if (t_PixcandP)
     delete t_PixcandP;
   if (t_PixcandPt)
     delete t_PixcandPt;
   if (t_PixcandEta)
     delete t_PixcandEta;
   if (t_PixcandPhi)
     delete t_PixcandPhi;
   if (t_PixcandMaxP)
     delete t_PixcandMaxP;
   if (t_PixTrkcandP)
     delete t_PixTrkcandP;
   if (t_PixTrkcandPt)
     delete t_PixTrkcandPt;
   if (t_PixTrkcandEta)
     delete t_PixTrkcandEta;
   if (t_PixTrkcandPhi)
     delete t_PixTrkcandPhi;
   if (t_PixTrkcandMaxP)
     delete t_PixTrkcandMaxP;
   if (t_PixTrkcandselTk)
     delete t_PixTrkcandselTk;
   if (t_NFcandP)
     delete t_NFcandP;
   if (t_NFcandPt)
     delete t_NFcandPt;
   if (t_NFcandEta)
     delete t_NFcandEta;
   if (t_NFcandPhi)
     delete t_NFcandPhi;
   if (t_NFcandEmip)
     delete t_NFcandEmip;
   if (t_NFTrkcandP)
     delete t_NFTrkcandP;
   if (t_NFTrkcandPt)
     delete t_NFTrkcandPt;
   if (t_NFTrkcandEta)
     delete t_NFTrkcandEta;
   if (t_NFTrkcandPhi)
     delete t_NFTrkcandPhi;
   if (t_NFTrkcandEmip)
     delete t_NFTrkcandEmip;
   if (t_NFTrkMinDR)
     delete t_NFTrkMinDR;
   if (t_NFTrkMinDP1)
     delete t_NFTrkMinDP1;
   if (t_NFTrkselTkFlag)
     delete t_NFTrkselTkFlag;
   if (t_NFTrkqltyFlag)
     delete t_NFTrkqltyFlag;
   if (t_NFTrkMissFlag)
     delete t_NFTrkMissFlag;
   if (t_NFTrkPVFlag)
     delete t_NFTrkPVFlag;
   if (t_NFTrkPropFlag)
     delete t_NFTrkPropFlag;
   if (t_NFTrkChgIsoFlag)
     delete t_NFTrkChgIsoFlag;
   if (t_NFTrkNeuIsoFlag)
     delete t_NFTrkNeuIsoFlag;
   if (t_NFTrkMipFlag)
     delete t_NFTrkMipFlag;
   if (t_ECone)
     delete t_ECone;
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 565 of file IsoTrig.cc.

References barrelRecHitsHandle_, beamSpotH_, ProducerSetup_cfi::BField, bField_, bfVal_, changed_, chgIsolation(), l1ct::clear(), doChgIsolTree_, doL2L3_, doMipCutTree_, doStudyIsol_, doTiming_, doTrkResTree_, endcapRecHitsHandle_, trigger::TriggerObject::energy(), trigger::TriggerObject::eta(), fastmath::etaphi(), etaPhiTrigger(), spr::find(), geo_, edm::EventSetup::getData(), getGoodTracks(), h_Filters, h_HLT, h_L1ObjEnergy, h_nHLT, h_nL3Objs, h_Pre, h_PreHLT, h_PreHLTwt, h_PreL1, h_PreL1wt, hbhe_, ValidationMatrix::hlt, HltComparatorCreateWorkflow::hltConfig, HLTPrescaleProvider::hltConfigProvider(), hltPrescaleProvider_, mps_fire::i, trigger::TriggerObject::id(), iEvent, recoMuon::in, MagneticField::inTesla(), edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, LumiDetails::kOCC1, label, leadPV_, spclusmultlumicorr_cfi::lumiProducer, LumiDetails::lumiValue(), trigger::TriggerObject::mass(), LogMessageMonitor_cff::modules, dqmiodumpmetadata::n, trigger::TriggerObject::phi(), pixelTrackRefsHB_, pixelTrackRefsHE_, pixelTracksSources_, reco::BeamSpot::position(), hltrates_dqm_sourceclient-live_cfg::prescales, HLTPrescaleProvider::prescaleValues(), edm::Handle< T >::product(), trigger::TriggerObject::pt(), trigger::TriggerObject::px(), trigger::TriggerObject::py(), trigger::TriggerObject::pz(), recVtxs_, findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, studyIsolation(), studyMipCut(), studyTiming(), studyTrigger(), StudyTrkEbyP(), tok_bs_, tok_EB_, tok_EE_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_l1cand_, tok_l2cand_, tok_lumi_, tok_magField_, tok_pixtk_, tok_pixtks_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, PDWG_DiPhoton_SD_cff::triggerEvent, L1TEGammaOffline_cfi::triggerNames, triggerResults, trigList_, trigNames_, trigPreList_, vec_, verbosity_, x, y, and z.

                                                                          {
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "Event starts====================================";
 
   int RunNo = iEvent.id().run();
 
   HLTConfigProvider const &hltConfig = hltPrescaleProvider_.hltConfigProvider();
 
   bField_ = &iSetup.getData(tok_magField_);
   geo_ = &iSetup.getData(tok_geom_);
   GlobalVector BField = bField_->inTesla(GlobalPoint(0, 0, 0));
   bfVal_ = BField.mag();
 
   trigger::TriggerEvent triggerEvent;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
   if (!triggerEventHandle.isValid()) {
     edm::LogWarning("IsoTrack") << "Error! Can't get the product hltTriggerSummaryAOD";
 
   } else {
     triggerEvent = *(triggerEventHandle.product());
   }
   const trigger::TriggerObjectCollection &TOC(triggerEvent.getObjects());
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(tok_trigRes_, triggerResults);
 
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
 
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle_);
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle_);
 
   iEvent.getByToken(tok_hbhe_, hbhe_);
 
   iEvent.getByToken(tok_recVtx_, recVtxs_);
   iEvent.getByToken(tok_bs_, beamSpotH_);
   if (!recVtxs_->empty() && !((*recVtxs_)[0].isFake())) {
     leadPV_ = math::XYZPoint((*recVtxs_)[0].x(), (*recVtxs_)[0].y(), (*recVtxs_)[0].z());
   } else if (beamSpotH_.isValid()) {
     leadPV_ = beamSpotH_->position();
   }
 
   if ((verbosity_ / 100) % 10 > 0) {
     edm::LogVerbatim("IsoTrack") << "Primary Vertex " << leadPV_;
     if (beamSpotH_.isValid())
       edm::LogVerbatim("IsoTrack") << "Beam Spot " << beamSpotH_->position();
   }
   pixelTrackRefsHE_.clear();
   pixelTrackRefsHB_.clear();
   for (unsigned int iPix = 0; iPix < pixelTracksSources_.size(); iPix++) {
     edm::Handle<reco::TrackCollection> iPixCol;
     iEvent.getByToken(tok_pixtks_[iPix], iPixCol);
     if (iPixCol.isValid()) {
       for (reco::TrackCollection::const_iterator pit = iPixCol->begin(); pit != iPixCol->end(); pit++) {
         if (iPix == 0)
           pixelTrackRefsHB_.push_back(reco::TrackRef(iPixCol, pit - iPixCol->begin()));
         pixelTrackRefsHE_.push_back(reco::TrackRef(iPixCol, pit - iPixCol->begin()));
       }
     }
   }
   if (doTiming_)
     getGoodTracks(iEvent, trkCollection);
 
   for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
     std::string FilterNames[7] = {"hltL1sL1SingleJet68",
                                   "hltIsolPixelTrackL2FilterHE",
                                   "ecalIsolPixelTrackFilterHE",
                                   "hltIsolPixelTrackL3FilterHE",
                                   "hltIsolPixelTrackL2FilterHB",
                                   "ecalIsolPixelTrackFilterHB",
                                   "hltIsolPixelTrackL3FilterHB"};
     std::string label = triggerEvent.filterTag(ifilter).label();
     for (int i = 0; i < 7; i++) {
       if (label == FilterNames[i])
         h_Filters->Fill(i);
     }
   }
   edm::InputTag lumiProducer("LumiProducer", "", "RECO");
   edm::Handle<LumiDetails> Lumid;
   iEvent.getLuminosityBlock().getByToken(tok_lumi_, Lumid);
   float mybxlumi = -1;
   if (Lumid.isValid())
     mybxlumi = Lumid->lumiValue(LumiDetails::kOCC1, iEvent.bunchCrossing()) * 6.37;
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "RunNo " << RunNo << " EvtNo " << iEvent.id().event() << " Lumi "
                                  << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing() << " mybxlumi "
                                  << mybxlumi;
   if (!triggerResults.isValid()) {
     edm::LogWarning("IsoTrack") << "Error! Can't get the product triggerResults";
     //      std::shared_ptr<cms::Exception> const & error = triggerResults.whyFailed();
     //      edm::LogWarning(error->category()) << error->what();
   } else {
     std::vector<std::string> modules;
     h_nHLT->Fill(triggerResults->size());
     const edm::TriggerNames &triggerNames = iEvent.triggerNames(*triggerResults);
 
     const std::vector<std::string> &triggerNames_ = triggerNames.triggerNames();
     if (verbosity_ % 10 > 1)
       edm::LogVerbatim("IsoTrack") << "number of HLTs " << triggerNames_.size();
     double preL1(-1), preHLT(-1), prescale(-1);
     int hlt(-1);
     for (unsigned int i = 0; i < triggerResults->size(); i++) {
       unsigned int triggerindx = hltConfig.triggerIndex(triggerNames_[i]);
       const std::vector<std::string> &moduleLabels(hltConfig.moduleLabels(triggerindx));
 
       for (unsigned int in = 0; in < trigNames_.size(); ++in) {
         //    if (triggerNames_[i].find(trigNames_[in].c_str())!=std::string::npos || triggerNames_[i]==" ") {
         if (triggerNames_[i].find(trigNames_[in]) != std::string::npos) {
           if (verbosity_ % 10 > 0)
             edm::LogVerbatim("IsoTrack") << "trigger that i want " << triggerNames_[i] << " accept "
                                          << triggerResults->accept(i);
           hlt = triggerResults->accept(i);
           h_HLT->Fill(hlt);
           //        if (hlt>0 || triggerNames_[i]==" ") {
           if (hlt > 0) {
             edm::Handle<reco::IsolatedPixelTrackCandidateCollection> Pixcands;
             iEvent.getByToken(tok_pixtk_, Pixcands);
             edm::Handle<trigger::TriggerFilterObjectWithRefs> L1cands;
             iEvent.getByToken(tok_l1cand_, L1cands);
 
             auto const prescales = hltPrescaleProvider_.prescaleValues<double>(iEvent, iSetup, triggerNames_[i]);
             preL1 = prescales.first;
             preHLT = prescales.second;
             prescale = preL1 * preHLT;
             if (verbosity_ % 10 > 0)
               edm::LogVerbatim("IsoTrack")
                   << triggerNames_[i] << " accept " << hlt << " preL1 " << preL1 << " preHLT " << preHLT;
             for (int iv = 0; iv < 3; ++iv)
               vec_[iv].clear();
             if (trigList_.find(RunNo) != trigList_.end()) {
               trigList_[RunNo] += 1;
             } else {
               trigList_.insert({RunNo, 1});
               trigPreList_.insert({RunNo, prescales});
             }
             //loop over all trigger filters in event (i.e. filters passed)
             for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
               std::vector<int> Keys;
               std::string label = triggerEvent.filterTag(ifilter).label();
               //loop over keys to objects passing this filter
               for (unsigned int imodule = 0; imodule < moduleLabels.size(); imodule++) {
                 if (label.find(moduleLabels[imodule]) != std::string::npos) {
                   if (verbosity_ % 10 > 0)
                     edm::LogVerbatim("IsoTrack") << "FILTERNAME " << label;
                   for (unsigned int ifiltrKey = 0; ifiltrKey < triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
                     Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
                     const trigger::TriggerObject &TO(TOC[Keys[ifiltrKey]]);
                     math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
                     if (label.find("L2Filter") != std::string::npos) {
                       vec_[1].push_back(v4);
                     } else if (label.find("L3Filter") != std::string::npos) {
                       vec_[2].push_back(v4);
                     } else {
                       vec_[0].push_back(v4);
                       h_L1ObjEnergy->Fill(TO.energy());
                     }
                     if (verbosity_ % 10 > 0)
                       edm::LogVerbatim("IsoTrack") << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta()
                                                    << " phi " << TO.phi() << " mass " << TO.mass() << " Id " << TO.id();
                   }
                 }
               }
             }
             std::vector<reco::TrackCollection::const_iterator> goodTks;
             if (doL2L3_) {
               h_nL3Objs->Fill(vec_[2].size());
               studyTrigger(trkCollection, goodTks);
             } else {
               if (trkCollection.isValid()) {
                 reco::TrackCollection::const_iterator trkItr;
                 for (trkItr = trkCollection->begin(); trkItr != trkCollection->end(); trkItr++)
                   goodTks.push_back(trkItr);
               }
             }
             // Now study isolation etc
             if (doStudyIsol_)
               studyIsolation(trkCollection, goodTks);
             if (doTrkResTree_)
               StudyTrkEbyP(trkCollection);
 
             std::pair<double, double> etaphi = etaPhiTrigger();
             edm::Handle<reco::IsolatedPixelTrackCandidateCollection> L2cands;
             iEvent.getByToken(tok_l2cand_, L2cands);
             if (!L2cands.isValid()) {
               if (verbosity_ % 10 > 0)
                 edm::LogVerbatim("IsoTrack") << " trigCand is not valid ";
             } else {
               if (doMipCutTree_)
                 studyMipCut(trkCollection, L2cands);
             }
             if (!pixelTracksSources_.empty())
               if (doChgIsolTree_ && !pixelTrackRefsHE_.empty())
                 chgIsolation(etaphi.first, etaphi.second, trkCollection, iEvent);
           }
           break;
         }
       }
     }
     h_PreL1->Fill(preL1);
     h_PreHLT->Fill(preHLT);
     h_Pre->Fill(prescale);
     h_PreL1wt->Fill(preL1, mybxlumi);
     h_PreHLTwt->Fill(preHLT, mybxlumi);
 
     // check if trigger names in (new) config
     //      edm::LogVerbatim("IsoTrack") << "changed " << changed_;
     if (changed_) {
       changed_ = false;
       if ((verbosity_ / 10) % 10 > 1) {
         edm::LogVerbatim("IsoTrack") << "New trigger menu found !!!";
         const unsigned int n(hltConfig.size());
         for (unsigned itrig = 0; itrig < triggerNames_.size(); itrig++) {
           unsigned int triggerindx = hltConfig.triggerIndex(triggerNames_[itrig]);
           if (triggerindx >= n)
             edm::LogVerbatim("IsoTrack") << triggerNames_[itrig] << " " << triggerindx << " does not exist in"
                                          << " the current menu";
           else
             edm::LogVerbatim("IsoTrack") << triggerNames_[itrig] << " " << triggerindx << " exists";
         }
       }
     }
   }
   if (doTiming_)
     studyTiming(iEvent);
 }

◆ beginJob()

void IsoTrig::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 881 of file IsoTrig.cc.

References ChgIsolnTree_, doChgIsolTree_, doL2L3_, doMipCutTree_, doTiming_, doTrkResTree_, fs_, h_dEta, h_dEtaL1, h_dinvPt, h_dP, h_dPhi, h_dPhiL1, h_dPt, h_dRL1, h_eCalo, h_eHcal, h_eMaxNearP, h_eMip, h_eNeutIso, h_EnIn, h_EnOut, h_eta, h_etaCalibTracks, h_etaMipTracks, h_Filters, h_HLT, h_L1ObjEnergy, h_mindR, h_MipEnMatch, h_MipEnNoMatch, h_nHLT, h_nL3Objs, h_p, h_phi, h_Pre, h_PreHLT, h_PreHLTwt, h_PreL1, h_PreL1wt, h_pt, mps_fire::i, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, GCP_tree_cfg::levels, TFileService::make(), MipCutTree_, pLimits_, AlCaHLTBitMon_QueryRunRegistry::string, t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, t_NFTrkselTkFlag, t_nGoodTk, t_nPixCand, t_nPixSeed, t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, t_PixTrkcandselTk, t_timeL2Prod, t_TrkhCone, t_TrkMissFlag, t_TrkNuIsolFlag, t_TrkP, t_TrkPVFlag, t_TrkqltyFlag, t_TrkselTkFlag, TimingTree_, and TrkResTree_.

                        {
   char hname[100], htit[100];
   std::string levels[20] = {"L1",        "L2",          "L3",        "Reco",      "RecoMatch",     "RecoNoMatch",
                             "L2Match",   "L2NoMatch",   "L3Match",   "L3NoMatch", "HLTTrk",        "HLTGoodTrk",
                             "HLTIsoTrk", "HLTMip",      "HLTSelect", "nonHLTTrk", "nonHLTGoodTrk", "nonHLTIsoTrk",
                             "nonHLTMip", "nonHLTSelect"};
   if (doTiming_) {
     TimingTree_ = fs_->make<TTree>("TimingTree", "TimingTree");
     t_timeL2Prod = new std::vector<double>();
     t_nPixCand = new std::vector<int>();
     t_nPixSeed = new std::vector<int>();
     t_nGoodTk = new std::vector<int>();
 
     TimingTree_->Branch("t_timeL2Prod", "std::vector<double>", &t_timeL2Prod);
     TimingTree_->Branch("t_nPixCand", "std::vector<int>", &t_nPixCand);
     TimingTree_->Branch("t_nPixSeed", "std::vector<int>", &t_nPixSeed);
     TimingTree_->Branch("t_nGoodTk", "std::vector<int>", &t_nGoodTk);
   }
   if (doTrkResTree_) {
     TrkResTree_ = fs_->make<TTree>("TrkRestree", "TrkResTree");
     t_TrkhCone = new std::vector<double>();
     t_TrkP = new std::vector<double>();
     t_TrkselTkFlag = new std::vector<bool>();
     t_TrkqltyFlag = new std::vector<bool>();
     t_TrkMissFlag = new std::vector<bool>();
     t_TrkPVFlag = new std::vector<bool>();
     t_TrkNuIsolFlag = new std::vector<bool>();
 
     TrkResTree_->Branch("t_TrkhCone", "std::vector<double>", &t_TrkhCone);
     TrkResTree_->Branch("t_TrkP", "std::vector<double>", &t_TrkP);
     TrkResTree_->Branch("t_TrkselTkFlag", "std::vector<bool>", &t_TrkselTkFlag);
     TrkResTree_->Branch("t_TrkqltyFlag", "std::vector<bool>", &t_TrkqltyFlag);
     TrkResTree_->Branch("t_TrkMissFlag", "std::vector<bool>", &t_TrkMissFlag);
     TrkResTree_->Branch("t_TrkPVFlag", "std::vector<bool>", &t_TrkPVFlag);
     TrkResTree_->Branch("t_TrkNuIsolFlag", "std::vector<bool>", &t_TrkNuIsolFlag);
   }
   if (doChgIsolTree_) {
     ChgIsolnTree_ = fs_->make<TTree>("ChgIsolnTree", "ChgIsolnTree");
     t_PixcandP = new std::vector<double>();
     t_PixcandPt = new std::vector<double>();
     t_PixcandEta = new std::vector<double>();
     t_PixcandPhi = new std::vector<double>();
     t_PixcandMaxP = new std::vector<std::vector<double>>();
     t_PixTrkcandP = new std::vector<double>();
     t_PixTrkcandPt = new std::vector<double>();
     t_PixTrkcandEta = new std::vector<double>();
     t_PixTrkcandPhi = new std::vector<double>();
     t_PixTrkcandMaxP = new std::vector<double>();
     t_PixTrkcandselTk = new std::vector<bool>();
 
     ChgIsolnTree_->Branch("t_PixcandP", "std::vector<double>", &t_PixcandP);
     ChgIsolnTree_->Branch("t_PixcandPt", "std::vector<double>", &t_PixcandPt);
     ChgIsolnTree_->Branch("t_PixcandEta", "std::vector<double>", &t_PixcandEta);
     ChgIsolnTree_->Branch("t_PixcandPhi", "std::vector<double>", &t_PixcandPhi);
     ChgIsolnTree_->Branch("t_PixcandMaxP", "std::vector<std::vector<double> >", &t_PixcandMaxP);
     ChgIsolnTree_->Branch("t_PixTrkcandP", "std::vector<double>", &t_PixTrkcandP);
     ChgIsolnTree_->Branch("t_PixTrkcandPt", "std::vector<double>", &t_PixTrkcandPt);
     ChgIsolnTree_->Branch("t_PixTrkcandEta", "std::vector<double>", &t_PixTrkcandEta);
     ChgIsolnTree_->Branch("t_PixTrkcandPhi", "std::vector<double>", &t_PixTrkcandPhi);
     ChgIsolnTree_->Branch("t_PixTrkcandMaxP", "std::vector<double>", &t_PixTrkcandMaxP);
     ChgIsolnTree_->Branch("t_PixTrkcandselTk", "std::vector<bool>", &t_PixTrkcandselTk);
   }
   if (doMipCutTree_) {
     MipCutTree_ = fs_->make<TTree>("MipCutTree", "MipCutTree");
     t_NFcandP = new std::vector<double>();
     t_NFcandPt = new std::vector<double>();
     t_NFcandEta = new std::vector<double>();
     t_NFcandPhi = new std::vector<double>();
     t_NFcandEmip = new std::vector<double>();
     t_NFTrkcandP = new std::vector<double>();
     t_NFTrkcandPt = new std::vector<double>();
     t_NFTrkcandEta = new std::vector<double>();
     t_NFTrkcandPhi = new std::vector<double>();
     t_NFTrkcandEmip = new std::vector<double>();
     t_NFTrkMinDR = new std::vector<double>();
     t_NFTrkMinDP1 = new std::vector<double>();
     t_NFTrkselTkFlag = new std::vector<bool>();
     t_NFTrkqltyFlag = new std::vector<bool>();
     t_NFTrkMissFlag = new std::vector<bool>();
     t_NFTrkPVFlag = new std::vector<bool>();
     t_NFTrkPropFlag = new std::vector<bool>();
     t_NFTrkChgIsoFlag = new std::vector<bool>();
     t_NFTrkNeuIsoFlag = new std::vector<bool>();
     t_NFTrkMipFlag = new std::vector<bool>();
     t_ECone = new std::vector<double>();
 
     MipCutTree_->Branch("t_NFcandP", "std::vector<double>", &t_NFcandP);
     MipCutTree_->Branch("t_NFcandPt", "std::vector<double>", &t_NFcandPt);
     MipCutTree_->Branch("t_NFcandEta", "std::vector<double>", &t_NFcandEta);
     MipCutTree_->Branch("t_NFcandPhi", "std::vector<double>", &t_NFcandPhi);
     MipCutTree_->Branch("t_NFcandEmip", "std::vector<double>", &t_NFcandEmip);
     MipCutTree_->Branch("t_NFTrkcandP", "std::vector<double>", &t_NFTrkcandP);
     MipCutTree_->Branch("t_NFTrkcandPt", "std::vector<double>", &t_NFTrkcandPt);
     MipCutTree_->Branch("t_NFTrkcandEta", "std::vector<double>", &t_NFTrkcandEta);
     MipCutTree_->Branch("t_NFTrkcandPhi", "std::vector<double>", &t_NFTrkcandPhi);
     MipCutTree_->Branch("t_NFTrkcandEmip", "std::vector<double>", &t_NFTrkcandEmip);
     MipCutTree_->Branch("t_NFTrkMinDR", "std::vector<double>", &t_NFTrkMinDR);
     MipCutTree_->Branch("t_NFTrkMinDP1", "std::vector<double>", &t_NFTrkMinDP1);
     MipCutTree_->Branch("t_NFTrkselTkFlag", "std::vector<bool>", &t_NFTrkselTkFlag);
     MipCutTree_->Branch("t_NFTrkqltyFlag", "std::vector<bool>", &t_NFTrkqltyFlag);
     MipCutTree_->Branch("t_NFTrkMissFlag", "std::vector<bool>", &t_NFTrkMissFlag);
     MipCutTree_->Branch("t_NFTrkPVFlag", "std::vector<bool>", &t_NFTrkPVFlag);
     MipCutTree_->Branch("t_NFTrkPropFlag", "std::vector<bool>", &t_NFTrkPropFlag);
     MipCutTree_->Branch("t_NFTrkChgIsoFlag", "std::vector<bool>", &t_NFTrkChgIsoFlag);
     MipCutTree_->Branch("t_NFTrkNeuIsoFlag", "std::vector<bool>", &t_NFTrkNeuIsoFlag);
     MipCutTree_->Branch("t_NFTrkMipFlag", "std::vector<bool>", &t_NFTrkMipFlag);
     MipCutTree_->Branch("t_ECone", "std::vector<double>", &t_ECone);
   }
   h_Filters = fs_->make<TH1I>("h_Filters", "Filter Accepts", 10, 0, 10);
   std::string FilterNames[7] = {"hltL1sL1SingleJet68",
                                 "hltIsolPixelTrackL2FilterHE",
                                 "ecalIsolPixelTrackFilterHE",
                                 "hltIsolPixelTrackL3FilterHE",
                                 "hltIsolPixelTrackL2FilterHB",
                                 "ecalIsolPixelTrackFilterHB",
                                 "hltIsolPixelTrackL3FilterHB"};
   for (int i = 0; i < 7; i++) {
     h_Filters->GetXaxis()->SetBinLabel(i + 1, FilterNames[i].c_str());
   }
 
   h_nHLT = fs_->make<TH1I>("h_nHLT", "Size of trigger Names", 1000, 1, 1000);
   h_HLT = fs_->make<TH1I>("h_HLT", "HLT accept", 3, -1, 2);
   h_PreL1 = fs_->make<TH1I>("h_PreL1", "L1 Prescale", 500, 0, 500);
   h_PreHLT = fs_->make<TH1I>("h_PreHLT", "HLT Prescale", 50, 0, 50);
   h_Pre = fs_->make<TH1I>("h_Pre", "Prescale", 3000, 0, 3000);
 
   h_PreL1wt = fs_->make<TH1D>("h_PreL1wt", "Weighted L1 Prescale", 500, 0, 500);
   h_PreHLTwt = fs_->make<TH1D>("h_PreHLTwt", "Weighted HLT Prescale", 500, 0, 100);
   h_L1ObjEnergy = fs_->make<TH1D>("h_L1ObjEnergy", "Energy of L1Object", 500, 0.0, 500.0);
 
   h_EnIn = fs_->make<TH1D>("h_EnInEcal", "EnergyIn Ecal", 200, 0.0, 20.0);
   h_EnOut = fs_->make<TH1D>("h_EnOutEcal", "EnergyOut Ecal", 200, 0.0, 20.0);
   h_MipEnMatch =
       fs_->make<TH2D>("h_MipEnMatch", "MipEn: HLT level vs Reco Level (Matched)", 200, 0.0, 20.0, 200, 0.0, 20.0);
   h_MipEnNoMatch = fs_->make<TH2D>(
       "h_MipEnNoMatch", "MipEn: HLT level vs Reco Level (No Match Found)", 200, 0.0, 20.0, 200, 0.0, 20.0);
 
   if (doL2L3_) {
     h_nL3Objs = fs_->make<TH1I>("h_nL3Objs", "Number of L3 objects", 10, 0, 10);
 
     std::string pairs[9] = {"L2L3",
                             "L2L3Match",
                             "L2L3NoMatch",
                             "L3Reco",
                             "L3RecoMatch",
                             "L3RecoNoMatch",
                             "NewFilterReco",
                             "NewFilterRecoMatch",
                             "NewFilterRecoNoMatch"};
     for (int ipair = 0; ipair < 9; ipair++) {
       sprintf(hname, "h_dEta%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta#eta for %s", pairs[ipair].c_str());
       h_dEta[ipair] = fs_->make<TH1D>(hname, htit, 200, -10.0, 10.0);
       h_dEta[ipair]->GetXaxis()->SetTitle("d#eta");
 
       sprintf(hname, "h_dPhi%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta#phi for %s", pairs[ipair].c_str());
       h_dPhi[ipair] = fs_->make<TH1D>(hname, htit, 140, -7.0, 7.0);
       h_dPhi[ipair]->GetXaxis()->SetTitle("d#phi");
 
       sprintf(hname, "h_dPt%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta dp_{T} for %s objects", pairs[ipair].c_str());
       h_dPt[ipair] = fs_->make<TH1D>(hname, htit, 400, -200.0, 200.0);
       h_dPt[ipair]->GetXaxis()->SetTitle("dp_{T} (GeV)");
 
       sprintf(hname, "h_dP%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta p for %s objects", pairs[ipair].c_str());
       h_dP[ipair] = fs_->make<TH1D>(hname, htit, 400, -200.0, 200.0);
       h_dP[ipair]->GetXaxis()->SetTitle("dP (GeV)");
 
       sprintf(hname, "h_dinvPt%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta (1/p_{T}) for %s objects", pairs[ipair].c_str());
       h_dinvPt[ipair] = fs_->make<TH1D>(hname, htit, 500, -0.4, 0.1);
       h_dinvPt[ipair]->GetXaxis()->SetTitle("d(1/p_{T})");
       sprintf(hname, "h_mindR%s", pairs[ipair].c_str());
       sprintf(htit, "min(#Delta R) for %s objects", pairs[ipair].c_str());
       h_mindR[ipair] = fs_->make<TH1D>(hname, htit, 500, 0.0, 1.0);
       h_mindR[ipair]->GetXaxis()->SetTitle("dR");
     }
 
     for (int lvl = 0; lvl < 2; lvl++) {
       sprintf(hname, "h_dEtaL1%s", levels[lvl + 1].c_str());
       sprintf(htit, "#Delta#eta for L1 and %s objects", levels[lvl + 1].c_str());
       h_dEtaL1[lvl] = fs_->make<TH1D>(hname, htit, 400, -10.0, 10.0);
 
       sprintf(hname, "h_dPhiL1%s", levels[lvl + 1].c_str());
       sprintf(htit, "#Delta#phi for L1 and %s objects", levels[lvl + 1].c_str());
       h_dPhiL1[lvl] = fs_->make<TH1D>(hname, htit, 280, -7.0, 7.0);
 
       sprintf(hname, "h_dRL1%s", levels[lvl + 1].c_str());
       sprintf(htit, "#Delta R for L1 and %s objects", levels[lvl + 1].c_str());
       h_dRL1[lvl] = fs_->make<TH1D>(hname, htit, 100, 0.0, 10.0);
     }
   }
 
   int levmin = (doL2L3_ ? 0 : 10);
   for (int ilevel = levmin; ilevel < 20; ilevel++) {
     sprintf(hname, "h_p%s", levels[ilevel].c_str());
     sprintf(htit, "p for %s objects", levels[ilevel].c_str());
     h_p[ilevel] = fs_->make<TH1D>(hname, htit, 100, 0.0, 500.0);
     h_p[ilevel]->GetXaxis()->SetTitle("p (GeV)");
 
     sprintf(hname, "h_pt%s", levels[ilevel].c_str());
     sprintf(htit, "p_{T} for %s objects", levels[ilevel].c_str());
     h_pt[ilevel] = fs_->make<TH1D>(hname, htit, 100, 0.0, 500.0);
     h_pt[ilevel]->GetXaxis()->SetTitle("p_{T} (GeV)");
 
     sprintf(hname, "h_eta%s", levels[ilevel].c_str());
     sprintf(htit, "#eta for %s objects", levels[ilevel].c_str());
     h_eta[ilevel] = fs_->make<TH1D>(hname, htit, 100, -5.0, 5.0);
     h_eta[ilevel]->GetXaxis()->SetTitle("#eta");
 
     sprintf(hname, "h_phi%s", levels[ilevel].c_str());
     sprintf(htit, "#phi for %s objects", levels[ilevel].c_str());
     h_phi[ilevel] = fs_->make<TH1D>(hname, htit, 70, -3.5, 3.50);
     h_phi[ilevel]->GetXaxis()->SetTitle("#phi");
   }
 
   std::string cuts[2] = {"HLTMatched", "HLTNotMatched"};
   std::string cuts2[2] = {"All", "Away from L1"};
   for (int icut = 0; icut < 2; icut++) {
     sprintf(hname, "h_eMip%s", cuts[icut].c_str());
     sprintf(htit, "eMip for %s tracks", cuts[icut].c_str());
     h_eMip[icut] = fs_->make<TH1D>(hname, htit, 200, 0.0, 10.0);
     h_eMip[icut]->GetXaxis()->SetTitle("E_{Mip} (GeV)");
 
     sprintf(hname, "h_eMaxNearP%s", cuts[icut].c_str());
     sprintf(htit, "eMaxNearP for %s tracks", cuts[icut].c_str());
     h_eMaxNearP[icut] = fs_->make<TH1D>(hname, htit, 240, -2.0, 10.0);
     h_eMaxNearP[icut]->GetXaxis()->SetTitle("E_{MaxNearP} (GeV)");
 
     sprintf(hname, "h_eNeutIso%s", cuts[icut].c_str());
     sprintf(htit, "eNeutIso for %s ", cuts[icut].c_str());
     h_eNeutIso[icut] = fs_->make<TH1D>(hname, htit, 200, 0.0, 10.0);
     h_eNeutIso[icut]->GetXaxis()->SetTitle("E_{NeutIso} (GeV)");
 
     for (int kcut = 0; kcut < 2; ++kcut) {
       for (int lim = 0; lim < 5; ++lim) {
         sprintf(hname, "h_etaCalibTracks%sCut%dLim%d", cuts[icut].c_str(), kcut, lim);
         sprintf(htit,
                 "#eta for %s isolated MIP tracks (%4.1f < p < %5.1f Gev/c %s)",
                 cuts[icut].c_str(),
                 pLimits_[lim],
                 pLimits_[lim + 1],
                 cuts2[kcut].c_str());
         h_etaCalibTracks[lim][icut][kcut] = fs_->make<TH1D>(hname, htit, 60, -30.0, 30.0);
         h_etaCalibTracks[lim][icut][kcut]->GetXaxis()->SetTitle("i#eta");
 
         sprintf(hname, "h_etaMipTracks%sCut%dLim%d", cuts[icut].c_str(), kcut, lim);
         sprintf(htit,
                 "#eta for %s charge isolated MIP tracks (%4.1f < p < %5.1f Gev/c %s)",
                 cuts[icut].c_str(),
                 pLimits_[lim],
                 pLimits_[lim + 1],
                 cuts2[kcut].c_str());
         h_etaMipTracks[lim][icut][kcut] = fs_->make<TH1D>(hname, htit, 60, -30.0, 30.0);
         h_etaMipTracks[lim][icut][kcut]->GetXaxis()->SetTitle("i#eta");
       }
     }
   }
 
   std::string ecut1[3] = {"all", "HLTMatched", "HLTNotMatched"};
   std::string ecut2[2] = {"without", "with"};
   int etac[48] = {-1,  -2,  -3,  -4,  -5,  -6,  -7,  -8,  -9, -10, -11, -12, -13, -14, -15, -16,
                   -17, -18, -19, -20, -21, -22, -23, -24, 1,  2,   3,   4,   5,   6,   7,   8,
                   9,   10,  11,  12,  13,  14,  15,  16,  17, 18,  19,  20,  21,  22,  23,  24};
   for (int icut = 0; icut < 6; icut++) {
     //    int i1 = (icut>3 ? 1 : 0);
     int i1 = (icut > 2 ? 1 : 0);
     int i2 = icut - i1 * 3;
     for (int kcut = 0; kcut < 48; kcut++) {
       for (int lim = 0; lim < 5; ++lim) {
         sprintf(hname, "h_eta%dEnHcal%s%s%d", etac[kcut], ecut1[i2].c_str(), ecut2[i1].c_str(), lim);
         sprintf(htit,
                 "HCAL energy for #eta=%d for %s tracks (p=%4.1f:%5.1f Gev) %s neutral isolation",
                 etac[kcut],
                 ecut1[i2].c_str(),
                 pLimits_[lim],
                 pLimits_[lim + 1],
                 ecut2[i1].c_str());
         h_eHcal[lim][icut][kcut] = fs_->make<TH1D>(hname, htit, 750, 0.0, 150.0);
         h_eHcal[lim][icut][kcut]->GetXaxis()->SetTitle("Energy (GeV)");
         sprintf(hname, "h_eta%dEnCalo%s%s%d", etac[kcut], ecut1[i2].c_str(), ecut2[i1].c_str(), lim);
         sprintf(htit,
                 "Calorimter energy for #eta=%d for %s tracks (p=%4.1f:%5.1f Gev) %s neutral isolation",
                 etac[kcut],
                 ecut1[i2].c_str(),
                 pLimits_[lim],
                 pLimits_[lim + 1],
                 ecut2[i1].c_str());
         h_eCalo[lim][icut][kcut] = fs_->make<TH1D>(hname, htit, 750, 0.0, 150.0);
         h_eCalo[lim][icut][kcut]->GetXaxis()->SetTitle("Energy (GeV)");
       }
     }
   }
 }

◆ beginRun()

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

overrideprivate

Definition at line 1201 of file IsoTrig.cc.

References changed_, hltPrescaleProvider_, HLTPrescaleProvider::init(), processName_, and edm::RunBase::run().

                                                                       {
   edm::LogWarning("IsoTrack") << "Run " << iRun.run() << " hltconfig.init "
                               << hltPrescaleProvider_.init(iRun, iSetup, processName_, changed_);
 }

◆ chgIsolation()

void IsoTrig::chgIsolation	(	double &	etaTriggered,
		double &	phiTriggered,
		edm::Handle< reco::TrackCollection > &	trkCollection,
		const edm::Event &	theEvent
	)

private

Definition at line 1825 of file IsoTrig.cc.

References a_charIsoR_, bField_, ALCARECOTkAlJpsiMuMu_cff::charge, spr::chargeIsolationCone(), ChgIsolnTree_, clearChgIsolnTreeVectors(), PbPb_ZMuSkimMuonDPG_cff::deltaR, l1ctLayer1_cff::dr, dR(), PVValHelper::dxy, PVValHelper::dz, PVValHelper::eta, geo_, getDistInCM(), GetEtaPhiAtEcal(), spr::goodTrack(), mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, dqmiolumiharvest::j, dqmdumpme::k, MainPageGenerator::l, leadPV_, RecoDisplacedMuonValidator_cfi::maxP, minPTrackValue_, beamSpotPI::nTracks, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), phi, pixelIsolationConeSizeAtEC_, pixelTrackRefsHE_, position, prelimCone_, spr::propagateCALO(), DiDispStaMuonMonitor_cfi::pt, pushChgIsolnTreeVecs(), multPhiCorr_741_25nsDY_cfi::px, reco::TrackBase::px(), multPhiCorr_741_25nsDY_cfi::py, reco::TrackBase::py(), reco::TrackBase::pz(), dttmaxenums::R, recVtxs_, edm::second(), selectionParameters_, tauUnbiasCone_, theTrackQuality_, verbosity_, vtxCutIsol_, vtxCutSeed_, and z.

Referenced by analyze().

                                                      {
   clearChgIsolnTreeVectors();
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "Inside chgIsolation() with eta/phi Triggered: " << etaTriggered << "/"
                                  << phiTriggered;
   std::vector<double> maxP;
 
   std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
   std::vector<spr::propagatedTrackDirection> trkCaloDirections1;
   spr::propagateCALO(trkCollection, geo_, bField_, theTrackQuality_, trkCaloDirections1, ((verbosity_ / 100) % 10 > 2));
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "Propagated TrkCollection";
   for (unsigned int k = 0; k < pixelIsolationConeSizeAtEC_.size(); ++k)
     maxP.push_back(0);
   unsigned i = pixelTrackRefsHE_.size();
   std::vector<std::pair<unsigned int, std::pair<double, double>>> VecSeedsatEC;
   //loop to select isolated tracks
   for (unsigned iS = 0; iS < pixelTrackRefsHE_.size(); iS++) {
     if (pixelTrackRefsHE_[iS]->p() > minPTrackValue_) {
       bool vtxMatch = false;
       //associate to vertex (in Z)
       unsigned int ivSel = recVtxs_->size();
       double minDZ = 100;
       for (unsigned int iv = 0; iv < recVtxs_->size(); ++iv) {
         if (fabs(pixelTrackRefsHE_[iS]->dz((*recVtxs_)[iv].position())) < minDZ) {
           minDZ = fabs(pixelTrackRefsHE_[iS]->dz((*recVtxs_)[iv].position()));
           ivSel = iv;
         }
       }
       //cut on dYX:
       if (ivSel == recVtxs_->size()) {
         vtxMatch = true;
       } else if (fabs(pixelTrackRefsHE_[iS]->dxy((*recVtxs_)[ivSel].position())) < vtxCutSeed_) {
         vtxMatch = true;
       }
       //select tracks not matched to triggered L1 jet
       double R = deltaR(etaTriggered, phiTriggered, pixelTrackRefsHE_[iS]->eta(), pixelTrackRefsHE_[iS]->phi());
       if (R > tauUnbiasCone_ && vtxMatch) {
         //propagate seed track to ECAL surface:
         std::pair<double, double> seedCooAtEC;
         // in case vertex is found:
         if (minDZ != 100)
           seedCooAtEC = GetEtaPhiAtEcal(pixelTrackRefsHE_[iS]->eta(),
                                         pixelTrackRefsHE_[iS]->phi(),
                                         pixelTrackRefsHE_[iS]->pt(),
                                         pixelTrackRefsHE_[iS]->charge(),
                                         (*recVtxs_)[ivSel].z());
         //in case vertex is not found:
         else
           seedCooAtEC = GetEtaPhiAtEcal(pixelTrackRefsHE_[iS]->eta(),
                                         pixelTrackRefsHE_[iS]->phi(),
                                         pixelTrackRefsHE_[iS]->pt(),
                                         pixelTrackRefsHE_[iS]->charge(),
                                         0);
         VecSeedsatEC.push_back(std::make_pair(iS, seedCooAtEC));
       }
     }
   }
   for (unsigned int l = 0; l < VecSeedsatEC.size(); l++) {
     unsigned int iSeed = VecSeedsatEC[l].first;
     math::XYZTLorentzVector v1(pixelTrackRefsHE_[iSeed]->px(),
                                pixelTrackRefsHE_[iSeed]->py(),
                                pixelTrackRefsHE_[iSeed]->pz(),
                                pixelTrackRefsHE_[iSeed]->p());
 
     for (unsigned int j = 0; j < VecSeedsatEC.size(); j++) {
       unsigned int iSurr = VecSeedsatEC[j].first;
       if (iSeed != iSurr) {
         //define preliminary cone around seed track impact point from which tracks will be extrapolated:
         //  edm::Ref<reco::IsolatedPixelTrackCandidateCollection> cand2ref =
         //    edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(L2cands, iSurr);
         if (deltaR(pixelTrackRefsHE_[iSeed]->eta(),
                    pixelTrackRefsHE_[iSeed]->phi(),
                    pixelTrackRefsHE_[iSurr]->eta(),
                    pixelTrackRefsHE_[iSurr]->phi()) < prelimCone_) {
           unsigned int ivSel = recVtxs_->size();
           double minDZ2 = 100;
           for (unsigned int iv = 0; iv < recVtxs_->size(); ++iv) {
             if (fabs(pixelTrackRefsHE_[iSurr]->dz((*recVtxs_)[iv].position())) < minDZ2) {
               minDZ2 = fabs(pixelTrackRefsHE_[iSurr]->dz((*recVtxs_)[iv].position()));
               ivSel = iv;
             }
           }
           //cut ot dXY:
           if (minDZ2 == 100 || fabs(pixelTrackRefsHE_[iSurr]->dxy((*recVtxs_)[ivSel].position())) < vtxCutIsol_) {
             //calculate distance at ECAL surface and update isolation:
             double dist = getDistInCM(VecSeedsatEC[i].second.first,
                                       VecSeedsatEC[i].second.second,
                                       VecSeedsatEC[j].second.first,
                                       VecSeedsatEC[j].second.second);
             for (unsigned int k = 0; k < pixelIsolationConeSizeAtEC_.size(); ++k) {
               if (dist < pixelIsolationConeSizeAtEC_[k]) {
                 if (pixelTrackRefsHE_[iSurr]->p() > maxP[k])
                   maxP[k] = pixelTrackRefsHE_[iSurr]->p();
               }
             }
           }
         }
       }
     }
 
     double conehmaxNearP = -1;
     bool selectTk = false;
     double mindR = 999.9;
     int nTracks = 0;
     math::XYZTLorentzVector mindRvec;
     for (trkDetItr = trkCaloDirections1.begin(); trkDetItr != trkCaloDirections1.end(); trkDetItr++, nTracks++) {
       int nNearTRKs = 0;
       const reco::Track *pTrack = &(*(trkDetItr->trkItr));
       math::XYZTLorentzVector v2(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
       double dr = dR(v1, v2);
       if (dr < mindR) {
         selectTk = spr::goodTrack(pTrack, leadPV_, selectionParameters_, ((verbosity_ / 100) % 10 > 1));
         conehmaxNearP = spr::chargeIsolationCone(
             nTracks, trkCaloDirections1, a_charIsoR_, nNearTRKs, ((verbosity_ / 100) % 10 > 1));
         mindR = dr;
         mindRvec = v2;
       }
     }
     pushChgIsolnTreeVecs(v1, mindRvec, maxP, conehmaxNearP, selectTk);
   }
   ChgIsolnTree_->Fill();
 }

◆ clearChgIsolnTreeVectors()

void IsoTrig::clearChgIsolnTreeVectors ( )

private

Definition at line 792 of file IsoTrig.cc.

References l1ct::clear(), mps_fire::i, t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, and t_PixTrkcandselTk.

Referenced by chgIsolation().

                                        {
   t_PixcandP->clear();
   t_PixcandPt->clear();
   t_PixcandEta->clear();
   t_PixcandPhi->clear();
   for (unsigned int i = 0; i < t_PixcandMaxP->size(); i++)
     t_PixcandMaxP[i].clear();
   t_PixcandMaxP->clear();
   t_PixTrkcandP->clear();
   t_PixTrkcandPt->clear();
   t_PixTrkcandEta->clear();
   t_PixTrkcandPhi->clear();
   t_PixTrkcandMaxP->clear();
   t_PixTrkcandselTk->clear();
 }

◆ clearMipCutTreeVectors()

void IsoTrig::clearMipCutTreeVectors ( )

private

Definition at line 808 of file IsoTrig.cc.

References t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, and t_NFTrkselTkFlag.

Referenced by studyMipCut().

                                      {
   t_NFcandP->clear();
   t_NFcandPt->clear();
   t_NFcandEta->clear();
   t_NFcandPhi->clear();
   t_NFcandEmip->clear();
   t_NFTrkcandP->clear();
   t_NFTrkcandPt->clear();
   t_NFTrkcandEta->clear();
   t_NFTrkcandPhi->clear();
   t_NFTrkcandEmip->clear();
   t_NFTrkMinDR->clear();
   t_NFTrkMinDP1->clear();
   t_NFTrkselTkFlag->clear();
   t_NFTrkqltyFlag->clear();
   t_NFTrkMissFlag->clear();
   t_NFTrkPVFlag->clear();
   t_NFTrkPropFlag->clear();
   t_NFTrkChgIsoFlag->clear();
   t_NFTrkNeuIsoFlag->clear();
   t_NFTrkMipFlag->clear();
   t_ECone->clear();
 }

◆ dEta()

double IsoTrig::dEta	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 2100 of file IsoTrig.cc.

Referenced by dR(), fillDifferences(), and studyTrigger().

2100 { return (vec1.eta() - vec2.eta()); }

vec1

std::vector< double > vec1

Definition: HCALResponse.h:15

vec2

std::vector< vec1 > vec2

Definition: HCALResponse.h:16

◆ dinvPt()

double IsoTrig::dinvPt	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 2129 of file IsoTrig.cc.

Referenced by fillDifferences().

                                                                                {
   return ((1 / vec1.pt()) - (1 / vec2.pt()));
 }

◆ dP()

double IsoTrig::dP	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 2125 of file IsoTrig.cc.

References funct::abs().

Referenced by fillDifferences().

                                                                            {
   return (std::abs(vec1.r() - vec2.r()));
 }

◆ dPhi()

double IsoTrig::dPhi	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 2102 of file IsoTrig.cc.

References M_PI.

Referenced by dR(), fillDifferences(), and studyTrigger().

                                                                              {
   double phi1 = vec1.phi();
   if (phi1 < 0)
     phi1 += 2.0 * M_PI;
   double phi2 = vec2.phi();
   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;
   return dphi;
 }

◆ dPt()

double IsoTrig::dPt	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 2123 of file IsoTrig.cc.

Referenced by fillDifferences().

2123 { return (vec1.pt() - vec2.pt()); }

vec1

std::vector< double > vec1

Definition: HCALResponse.h:15

vec2

std::vector< vec1 > vec2

Definition: HCALResponse.h:16

◆ dR()

double IsoTrig::dR	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 2117 of file IsoTrig.cc.

References dEta(), dPhi(), and mathSSE::sqrt().

Referenced by chgIsolation(), fillDifferences(), getGoodTracks(), studyIsolation(), studyMipCut(), and studyTrigger().

                                                                            {
   double deta = dEta(vec1, vec2);
   double dphi = dPhi(vec1, vec2);
   return std::sqrt(deta * deta + dphi * dphi);
 }

◆ endJob()

void IsoTrig::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1179 of file IsoTrig.cc.

References fs_, g_Accepts, g_Pre, g_PreHLT, g_PreL1, mps_fire::i, TFileService::make(), maxRunNo_, minRunNo_, dqmiodumpmetadata::n, DeadROC_duringRun::runNum, trigList_, trigNames_, and trigPreList_.

Referenced by o2olib.O2ORunMgr::executeJob().

                      {
   edm::LogWarning("IsoTrack") << trigNames_.size() << "Triggers were run. RunNo vs HLT accepts for";
   for (unsigned int i = 0; i < trigNames_.size(); ++i)
     edm::LogWarning("IsoTrack") << "[" << i << "]: " << trigNames_[i];
   unsigned int n = maxRunNo_ - minRunNo_ + 1;
   g_Pre = fs_->make<TH1D>("h_PrevsRN", "PreScale Vs Run Number", n, minRunNo_, maxRunNo_);
   g_PreL1 = fs_->make<TH1D>("h_PreL1vsRN", "L1 PreScale Vs Run Number", n, minRunNo_, maxRunNo_);
   g_PreHLT = fs_->make<TH1D>("h_PreHLTvsRN", "HLT PreScale Vs Run Number", n, minRunNo_, maxRunNo_);
   g_Accepts = fs_->make<TH1I>("h_HLTAcceptsvsRN", "HLT Accepts Vs Run Number", n, minRunNo_, maxRunNo_);
 
   for (auto const &[runNum, nAccept] : trigList_) {
     auto const &triggerPrescales = trigPreList_[runNum];
     auto const preL1 = triggerPrescales.first;
     auto const preHLT = triggerPrescales.second;
     edm::LogVerbatim("IsoTrack") << runNum << " " << nAccept << " " << preL1 << " " << preHLT;
     g_Accepts->Fill(runNum, nAccept);
     g_PreL1->Fill(runNum, preL1);
     g_PreHLT->Fill(runNum, preHLT);
     g_Pre->Fill(runNum, preL1 * preHLT);
   }
 }

◆ endRun()

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

inlineoverrideprivate

Definition at line 108 of file IsoTrig.cc.

108 {}

◆ etaPhiTrigger()

std::pair< double, double > IsoTrig::etaPhiTrigger ( )

private

Definition at line 2133 of file IsoTrig.cc.

References PVValHelper::eta, dqmdumpme::k, phi, DiDispStaMuonMonitor_cfi::pt, muonTiming_cfi::ptmax, and vec_.

Referenced by analyze().

                                                {
   double eta(0), phi(0), ptmax(0);
   for (unsigned int k = 0; k < vec_[0].size(); ++k) {
     if (k == 0) {
       eta = vec_[0][k].eta();
       phi = vec_[0][k].phi();
       ptmax = vec_[0][k].pt();
     } else if (vec_[0][k].pt() > ptmax) {
       eta = vec_[0][k].eta();
       phi = vec_[0][k].phi();
       ptmax = vec_[0][k].pt();
     }
   }
   return std::pair<double, double>(eta, phi);
 }

◆ fillCuts()

void IsoTrig::fillCuts	(	int	indx,
		double	eMipDR,
		double	conehmaxNearP,
		double	e_inCone,
		math::XYZTLorentzVector &	vec,
		int	ieta,
		bool	cut
	)

private

Definition at line 2063 of file IsoTrig.cc.

References TkAlMuonSelectors_cfi::cut, cutCharge_, cutMip_, cutNeutral_, h_eMaxNearP, h_eMip, h_eNeutIso, h_etaCalibTracks, h_etaMipTracks, hcalRecHitTable_cff::ieta, and pLimits_.

Referenced by studyIsolation().

                                                                                                                     {
   h_eMip[indx]->Fill(eMipDR);
   h_eMaxNearP[indx]->Fill(conehmaxNearP);
   h_eNeutIso[indx]->Fill(e_inCone);
   if ((conehmaxNearP < cutCharge_) && (eMipDR < cutMip_)) {
     for (int lim = 0; lim < 5; ++lim) {
       if ((vec.r() > pLimits_[lim]) && (vec.r() <= pLimits_[lim + 1])) {
         h_etaMipTracks[lim][indx][0]->Fill((double)(ieta));
         if (cut)
           h_etaMipTracks[lim][indx][1]->Fill((double)(ieta));
         if (e_inCone < cutNeutral_) {
           h_etaCalibTracks[lim][indx][0]->Fill((double)(ieta));
           if (cut)
             h_etaCalibTracks[lim][indx][1]->Fill((double)(ieta));
         }
       }
     }
   }
 }

◆ fillDescriptions()

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

static

Definition at line 517 of file IsoTrig.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, and triggerMatcherToHLTDebug_cfi::tags.

                                                                          {
   std::vector<std::string> triggers = {"HLT_IsoTrackHB"};
   std::vector<edm::InputTag> tags = {edm::InputTag("hltHITPixelTracksHB"), edm::InputTag("hltHITPixelTracksHE")};
   std::vector<double> cones = {35.0, 40.0, 45.0, 50.0, 55.0, 60.0, 63.9, 70.0};
   edm::ParameterSetDescription desc;
   desc.addUntracked<std::vector<std::string>>("Triggers", triggers);
   desc.addUntracked<edm::InputTag>("pixCandTag", edm::InputTag(" "));
   desc.addUntracked<edm::InputTag>("l1CandTag", edm::InputTag("hltL1sV0SingleJet60"));
   desc.addUntracked<edm::InputTag>("l2CandTag", edm::InputTag("isolEcalPixelTrackProd"));
   desc.addUntracked<bool>("doL2L3", false);
   desc.addUntracked<bool>("doTimingTree", false);
   desc.addUntracked<bool>("doMipCutTree", false);
   desc.addUntracked<bool>("doTrkResTree", true);
   desc.addUntracked<bool>("doChgIsolTree", false);
   desc.addUntracked<bool>("doStudyIsol", false);
   desc.addUntracked<int>("verbosity", 0);
   desc.addUntracked<std::string>("processName", "HLT");
   desc.addUntracked<std::string>("trackQuality", "highPurity");
   desc.addUntracked<double>("minTrackPt", 10.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>("isolationL1", 1.0);
   desc.addUntracked<double>("coneRadius", 34.98);
   desc.addUntracked<double>("coneRadiusMIP", 14.0);
   desc.addUntracked<double>("coneRadiusNeut1", 21.0);
   desc.addUntracked<double>("coneRadiusNeut2", 29.0);
   desc.addUntracked<double>("cutMIP", 1.0);
   desc.addUntracked<double>("chargeIsolation", 2.0);
   desc.addUntracked<double>("neutralIsolation", 2.0);
   desc.addUntracked<int>("minRun", 190456);
   desc.addUntracked<int>("maxRun", 203002);
   desc.addUntracked<std::vector<edm::InputTag>>("pixelTracksSources", tags);
   desc.addUntracked<std::vector<double>>("pixelIsolationConeSizeAtEC", cones);
   desc.addUntracked<double>("minPTrackValue", 0.0);
   desc.addUntracked<double>("vertexCutSeed", 101.0);
   desc.addUntracked<double>("vertexCutIsol", 101.0);
   desc.addUntracked<double>("tauUnbiasCone", 1.2);
   desc.addUntracked<double>("prelimCone", 1.0);
   desc.add<unsigned int>("stageL1Trigger", 1);
   descriptions.add("isoTrigDefault", desc);
 }

◆ fillDifferences()

void IsoTrig::fillDifferences	(	int	indx,
		math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2,
		bool	debug
	)

private

Definition at line 2045 of file IsoTrig.cc.

References debug, dEta(), dinvPt(), dP(), Calorimetry_cff::dp, dPhi(), dPt(), l1ctLayer1_cff::dr, dR(), h_dEta, h_dinvPt, h_dP, h_dPhi, h_dPt, and h_mindR.

Referenced by studyMipCut(), and studyTrigger().

                                                                                                             {
   double dr = dR(vec1, vec2);
   double deta = dEta(vec1, vec2);
   double dphi = dPhi(vec1, vec2);
   double dpt = dPt(vec1, vec2);
   double dp = dP(vec1, vec2);
   double dinvpt = dinvPt(vec1, vec2);
   h_dEta[indx]->Fill(deta);
   h_dPhi[indx]->Fill(dphi);
   h_dPt[indx]->Fill(dpt);
   h_dP[indx]->Fill(dp);
   h_dinvPt[indx]->Fill(dinvpt);
   h_mindR[indx]->Fill(dr);
   if (debug)
     edm::LogVerbatim("IsoTrack") << "mindR for index " << indx << " is " << dr << " deta " << deta << " dphi " << dphi
                                  << " dpt " << dpt << " dinvpt " << dinvpt;
 }

◆ fillEnergy()

void IsoTrig::fillEnergy	(	int	indx,
		int	ieta,
		double	hCone,
		double	eMipDR,
		math::XYZTLorentzVector &	vec
	)

private

Definition at line 2084 of file IsoTrig.cc.

References h_eCalo, h_eHcal, hcalRecHitTable_cff::ieta, GetRecoTauVFromDQM_MC_cff::kk, and pLimits_.

Referenced by studyIsolation().

                                                                                                     {
   int kk(-1);
   if (ieta > 0 && ieta < 25)
     kk = 23 + ieta;
   else if (ieta > -25 && ieta < 0)
     kk = -(ieta + 1);
   if (kk >= 0 && eMipDR > 0.01 && hCone > 1.0) {
     for (int lim = 0; lim < 5; ++lim) {
       if ((vec.r() > pLimits_[lim]) && (vec.r() <= pLimits_[lim + 1])) {
         h_eHcal[lim][indx][kk]->Fill(hCone);
         h_eCalo[lim][indx][kk]->Fill(hCone + eMipDR);
       }
     }
   }
 }

◆ fillHist()

void IsoTrig::fillHist	(	int	indx,
		math::XYZTLorentzVector &	vec
	)

private

Definition at line 2038 of file IsoTrig.cc.

References h_eta, h_p, h_phi, and h_pt.

Referenced by studyIsolation(), and studyTrigger().

                                                            {
   h_p[indx]->Fill(vec.r());
   h_pt[indx]->Fill(vec.pt());
   h_eta[indx]->Fill(vec.eta());
   h_phi[indx]->Fill(vec.phi());
 }

◆ getDistInCM()

double IsoTrig::getDistInCM	(	double	eta1,
		double	phi1,
		double	eta2,
		double	phi2
	)

private

Definition at line 2212 of file IsoTrig.cc.

References angle(), funct::cos(), HLT_2024v13_cff::eta1, HLT_2024v13_cff::eta2, JetChargeProducer_cfi::exp, M_PI, rEB_, funct::sin(), funct::tan(), and zEE_.

Referenced by chgIsolation().

                                                                               {
   double Rec;
   double theta1 = 2 * atan(exp(-eta1));
   double theta2 = 2 * atan(exp(-eta2));
   if (fabs(eta1) < 1.479)
     Rec = rEB_;  //radius of ECAL barrel
   else if (fabs(eta1) > 1.479 && fabs(eta1) < 7.0)
     Rec = tan(theta1) * zEE_;  //distance from IP to ECAL endcap
   else
     return 1000;
 
   //|vect| times tg of acos(scalar product)
   double angle =
       acos((sin(theta1) * sin(theta2) * (sin(phi1) * sin(phi2) + cos(phi1) * cos(phi2)) + cos(theta1) * cos(theta2)));
   if (angle < 0.5 * M_PI)
     return fabs((Rec / sin(theta1)) * tan(angle));
   else
     return 1000;
 }

◆ GetEtaPhiAtEcal()

std::pair< double, double > IsoTrig::GetEtaPhiAtEcal	(	double	etaIP,
		double	phiIP,
		double	pT,
		int	charge,
		double	vtxZ
	)

private

Definition at line 2149 of file IsoTrig.cc.

References simBeamSpotPI::alpha, bfVal_, ALCARECOTkAlJpsiMuMu_cff::charge, SiPixelRawToDigiRegional_cfi::deltaPhi, JetChargeProducer_cfi::exp, dqm-mbProfile::log, M_PI, pv::pT, rEB_, funct::sin(), funct::tan(), theta(), z, and zEE_.

Referenced by chgIsolation().

                                                                                                                {
   double deltaPhi = 0;
   double etaEC = 100;
   double phiEC = 100;
 
   double Rcurv = 9999999;
   if (bfVal_ != 0)
     Rcurv = pT * 33.3 * 100 / (bfVal_ * 10);  //r(m)=pT(GeV)*33.3/B(kG)
 
   double ecDist = zEE_;
   double ecRad = rEB_;  //radius of ECAL barrel (cm)
   double theta = 2 * atan(exp(-etaIP));
   double zNew = 0;
   if (theta > 0.5 * M_PI)
     theta = M_PI - theta;
   if (fabs(etaIP) < 1.479) {
     if ((0.5 * ecRad / Rcurv) > 1) {
       etaEC = 10000;
       deltaPhi = 0;
     } else {
       deltaPhi = -charge * asin(0.5 * ecRad / Rcurv);
       double alpha1 = 2 * asin(0.5 * ecRad / Rcurv);
       double z = ecRad / tan(theta);
       if (etaIP > 0)
         zNew = z * (Rcurv * alpha1) / ecRad + vtxZ;  //new z-coordinate of track
       else
         zNew = -z * (Rcurv * alpha1) / ecRad + vtxZ;  //new z-coordinate of track
       double zAbs = fabs(zNew);
       if (zAbs < ecDist) {
         etaEC = -log(tan(0.5 * atan(ecRad / zAbs)));
         deltaPhi = -charge * asin(0.5 * ecRad / Rcurv);
       }
       if (zAbs > ecDist) {
         zAbs = (fabs(etaIP) / etaIP) * ecDist;
         double Zflight = fabs(zAbs - vtxZ);
         double alpha = (Zflight * ecRad) / (z * Rcurv);
         double Rec = 2 * Rcurv * sin(alpha / 2);
         deltaPhi = -charge * alpha / 2;
         etaEC = -log(tan(0.5 * atan(Rec / ecDist)));
       }
     }
   } else {
     zNew = (fabs(etaIP) / etaIP) * ecDist;
     double Zflight = fabs(zNew - vtxZ);
     double Rvirt = fabs(Zflight * tan(theta));
     double Rec = 2 * Rcurv * sin(Rvirt / (2 * Rcurv));
     deltaPhi = -(charge) * (Rvirt / (2 * Rcurv));
     etaEC = -log(tan(0.5 * atan(Rec / ecDist)));
   }
 
   if (zNew < 0)
     etaEC = -etaEC;
   phiEC = phiIP + deltaPhi;
 
   if (phiEC < -M_PI)
     phiEC += 2 * M_PI;
   if (phiEC > M_PI)
     phiEC -= 2 * M_PI;
 
   std::pair<double, double> retVal(etaEC, phiEC);
   return retVal;
 }

◆ getGoodTracks()

void IsoTrig::getGoodTracks	(	const edm::Event &	iEvent,
		edm::Handle< reco::TrackCollection > &	trkCollection
	)

private

Definition at line 1952 of file IsoTrig.cc.

References a_charIsoR_, a_mipR_, barrelRecHitsHandle_, bField_, spr::chargeIsolationCone(), funct::cos(), dR(), spr::eCone_ecal(), endcapRecHitsHandle_, geo_, trigger::TriggerRefsCollections::getObjects(), spr::goodTrack(), iEvent, cuy::ii, edm::HandleBase::isValid(), leadPV_, beamSpotPI::nTracks, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), spr::propagateCALO(), DiDispStaMuonMonitor_cfi::pt, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters_, funct::sin(), t_nGoodTk, nnet::tanh(), theTrackQuality_, tok_l1cand_, trigger::TriggerL1CenJet, trigger::TriggerL1ForJet, trigger::TriggerL1TauJet, and verbosity_.

Referenced by analyze().

                                                                                                  {
   t_nGoodTk->clear();
   std::vector<int> nGood(4, 0);
   if (trkCollection.isValid()) {
     std::vector<spr::propagatedTrackDirection> trkCaloDirections;
     spr::propagateCALO(
         trkCollection, geo_, bField_, theTrackQuality_, trkCaloDirections, ((verbosity_ / 100) % 10 > 2));
 
     // get the trigger jet
     edm::Handle<trigger::TriggerFilterObjectWithRefs> l1trigobj;
     iEvent.getByToken(tok_l1cand_, l1trigobj);
 
     std::vector<edm::Ref<l1extra::L1JetParticleCollection>> l1tauobjref;
     l1trigobj->getObjects(trigger::TriggerL1TauJet, l1tauobjref);
     std::vector<edm::Ref<l1extra::L1JetParticleCollection>> l1jetobjref;
     l1trigobj->getObjects(trigger::TriggerL1CenJet, l1jetobjref);
     std::vector<edm::Ref<l1extra::L1JetParticleCollection>> l1forjetobjref;
     l1trigobj->getObjects(trigger::TriggerL1ForJet, l1forjetobjref);
 
     double ptTriggered(-10), etaTriggered(-100), phiTriggered(-100);
     for (unsigned int p = 0; p < l1tauobjref.size(); p++) {
       if (l1tauobjref[p]->pt() > ptTriggered) {
         ptTriggered = l1tauobjref[p]->pt();
         phiTriggered = l1tauobjref[p]->phi();
         etaTriggered = l1tauobjref[p]->eta();
       }
     }
     for (unsigned int p = 0; p < l1jetobjref.size(); p++) {
       if (l1jetobjref[p]->pt() > ptTriggered) {
         ptTriggered = l1jetobjref[p]->pt();
         phiTriggered = l1jetobjref[p]->phi();
         etaTriggered = l1jetobjref[p]->eta();
       }
     }
     for (unsigned int p = 0; p < l1forjetobjref.size(); p++) {
       if (l1forjetobjref[p]->pt() > ptTriggered) {
         ptTriggered = l1forjetobjref[p]->pt();
         phiTriggered = l1forjetobjref[p]->phi();
         etaTriggered = l1forjetobjref[p]->eta();
       }
     }
     double pTriggered = ptTriggered * cosh(etaTriggered);
     math::XYZTLorentzVector pTrigger(
         ptTriggered * cos(phiTriggered), ptTriggered * sin(phiTriggered), pTriggered * tanh(etaTriggered), pTriggered);
 
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
     unsigned int nTracks(0);
     for (trkDetItr = trkCaloDirections.begin(); trkDetItr != trkCaloDirections.end(); trkDetItr++, nTracks++) {
       const reco::Track *pTrack = &(*(trkDetItr->trkItr));
       math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
       bool selectTk = spr::goodTrack(pTrack, leadPV_, selectionParameters_, ((verbosity_ / 100) % 10 > 1));
       double mindR = dR(v4, pTrigger);
       if ((verbosity_ / 100) % 10 > 1)
         edm::LogVerbatim("IsoTrack") << "Track ECAL " << trkDetItr->okECAL << " HCAL " << trkDetItr->okHCAL << " Flag "
                                      << selectTk;
       if (selectTk && trkDetItr->okECAL && trkDetItr->okHCAL && mindR > 1.0) {
         int nRH_eMipDR(0), nNearTRKs(0);
         double eMipDR = spr::eCone_ecal(geo_,
                                         barrelRecHitsHandle_,
                                         endcapRecHitsHandle_,
                                         trkDetItr->pointHCAL,
                                         trkDetItr->pointECAL,
                                         a_mipR_,
                                         trkDetItr->directionECAL,
                                         nRH_eMipDR);
         double conehmaxNearP =
             spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, ((verbosity_ / 100) % 10 > 1));
         if (conehmaxNearP < 2.0 && eMipDR < 1.0) {
           if (pTrack->p() >= 20 && pTrack->p() < 30) {
             ++nGood[0];
           } else if (pTrack->p() >= 30 && pTrack->p() < 40) {
             ++nGood[1];
           } else if (pTrack->p() >= 40 && pTrack->p() < 60) {
             ++nGood[2];
           } else if (pTrack->p() >= 60 && pTrack->p() < 100) {
             ++nGood[3];
           }
         }
       }
     }
   }
 
   for (unsigned int ii = 0; ii < nGood.size(); ++ii)
     t_nGoodTk->push_back(nGood[ii]);
 }

◆ pushChgIsolnTreeVecs()

void IsoTrig::pushChgIsolnTreeVecs	(	math::XYZTLorentzVector &	Pixcand,
		math::XYZTLorentzVector &	Trkcand,
		std::vector< double > &	PixMaxP,
		double &	TrkMaxP,
		bool &	selTk
	)

private

Definition at line 832 of file IsoTrig.cc.

References t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, and t_PixTrkcandselTk.

Referenced by chgIsolation().

                                                 {
   t_PixcandP->push_back(Pixcand.r());
   t_PixcandPt->push_back(Pixcand.pt());
   t_PixcandEta->push_back(Pixcand.eta());
   t_PixcandPhi->push_back(Pixcand.phi());
   t_PixcandMaxP->push_back(PixMaxP);
   t_PixTrkcandP->push_back(Trkcand.r());
   t_PixTrkcandPt->push_back(Trkcand.pt());
   t_PixTrkcandEta->push_back(Trkcand.eta());
   t_PixTrkcandPhi->push_back(Trkcand.phi());
   t_PixTrkcandMaxP->push_back(TrkMaxP);
   t_PixTrkcandselTk->push_back(selTk);
 }

◆ pushMipCutTreeVecs()

void IsoTrig::pushMipCutTreeVecs	(	math::XYZTLorentzVector &	NFcand,
		math::XYZTLorentzVector &	Trkcand,
		double &	EmipNFcand,
		double &	EmipTrkcand,
		double &	mindR,
		double &	mindP1,
		std::vector< bool > &	Flags,
		double	hCone
	)

private

Definition at line 850 of file IsoTrig.cc.

References t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, and t_NFTrkselTkFlag.

Referenced by studyMipCut().

                                                {
   t_NFcandP->push_back(NFcand.r());
   t_NFcandPt->push_back(NFcand.pt());
   t_NFcandEta->push_back(NFcand.eta());
   t_NFcandPhi->push_back(NFcand.phi());
   t_NFcandEmip->push_back(EmipNFcand);
   t_NFTrkcandP->push_back(Trkcand.r());
   t_NFTrkcandPt->push_back(Trkcand.pt());
   t_NFTrkcandEta->push_back(Trkcand.eta());
   t_NFTrkcandPhi->push_back(Trkcand.phi());
   t_NFTrkcandEmip->push_back(EmipTrkcand);
   t_NFTrkMinDR->push_back(mindR);
   t_NFTrkMinDP1->push_back(mindP1);
   t_NFTrkselTkFlag->push_back(Flags[0]);
   t_NFTrkqltyFlag->push_back(Flags[1]);
   t_NFTrkMissFlag->push_back(Flags[2]);
   t_NFTrkPVFlag->push_back(Flags[3]);
   t_NFTrkPropFlag->push_back(Flags[4]);
   t_NFTrkChgIsoFlag->push_back(Flags[5]);
   t_NFTrkNeuIsoFlag->push_back(Flags[6]);
   t_NFTrkMipFlag->push_back(Flags[7]);
   t_ECone->push_back(hCone);
 }

◆ studyIsolation()

void IsoTrig::studyIsolation	(	edm::Handle< reco::TrackCollection > &	trkCollection,
		std::vector< reco::TrackCollection::const_iterator > &	goodTks
	)

private

Definition at line 1688 of file IsoTrig.cc.

References a_charIsoR_, a_coneR_, a_mipR_, a_neutR1_, a_neutR2_, barrelRecHitsHandle_, bField_, spr::chargeIsolationCone(), cutCharge_, cutMip_, cutNeutral_, PbPb_ZMuSkimMuonDPG_cff::deltaR, hcalRecHitTable_cff::detId, doL2L3_, dR(), dr_L1_, StorageManager_cfg::e1, spr::eCone_ecal(), spr::eCone_hcal(), endcapRecHitsHandle_, fillCuts(), fillEnergy(), fillHist(), spr::find(), geo_, spr::goodTrack(), hbhe_, hit::id, hcalRecHitTable_cff::ieta, edm::HandleBase::isValid(), dqmdumpme::k, leadPV_, beamSpotPI::nTracks, reco::TrackBase::p(), spr::propagateCALO(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters_, theTrackQuality_, vec_, and verbosity_.

Referenced by analyze().

                                                                                     {
   if (trkCollection.isValid()) {
     std::vector<spr::propagatedTrackDirection> trkCaloDirections;
     spr::propagateCALO(
         trkCollection, geo_, bField_, theTrackQuality_, trkCaloDirections, ((verbosity_ / 100) % 10 > 2));
 
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
     if ((verbosity_ / 1000) % 10 > 1) {
       edm::LogVerbatim("IsoTrack") << "n of barrelRecHitsHandle " << barrelRecHitsHandle_->size();
       for (EcalRecHitCollection::const_iterator hit = barrelRecHitsHandle_->begin(); hit != barrelRecHitsHandle_->end();
            ++hit) {
         edm::LogVerbatim("IsoTrack") << "hit : detid(ieta,iphi) " << (EBDetId)(hit->id()) << " time " << hit->time()
                                      << " energy " << hit->energy();
       }
       edm::LogVerbatim("IsoTrack") << "n of endcapRecHitsHandle " << endcapRecHitsHandle_->size();
       for (EcalRecHitCollection::const_iterator hit = endcapRecHitsHandle_->begin(); hit != endcapRecHitsHandle_->end();
            ++hit) {
         edm::LogVerbatim("IsoTrack") << "hit : detid(ieta,iphi) " << (EEDetId)(hit->id()) << " time " << hit->time()
                                      << " energy " << hit->energy();
       }
       edm::LogVerbatim("IsoTrack") << "n of hbhe " << hbhe_->size();
       for (HBHERecHitCollection::const_iterator hit = hbhe_->begin(); hit != hbhe_->end(); ++hit) {
         edm::LogVerbatim("IsoTrack") << "hit : detid(ieta,iphi) " << hit->id() << " time " << hit->time() << " energy "
                                      << hit->energy();
       }
     }
     unsigned int nTracks = 0;
     int ieta = 999;
     for (trkDetItr = trkCaloDirections.begin(); trkDetItr != trkCaloDirections.end(); trkDetItr++, nTracks++) {
       bool l3Track = (std::find(goodTks.begin(), goodTks.end(), trkDetItr->trkItr) != goodTks.end());
       const reco::Track *pTrack = &(*(trkDetItr->trkItr));
       math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
       bool selectTk = spr::goodTrack(pTrack, leadPV_, selectionParameters_, ((verbosity_ / 100) % 10 > 1));
       double eMipDR = 9999., e_inCone = 0, conehmaxNearP = 0, mindR = 999.9, hCone = 0;
       if (trkDetItr->okHCAL) {
         HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
         ieta = detId.ieta();
       }
       for (unsigned k = 0; k < vec_[0].size(); ++k) {
         double deltaR = dR(v4, vec_[0][k]);
         if (deltaR < mindR)
           mindR = deltaR;
       }
       if ((verbosity_ / 100) % 10 > 1)
         edm::LogVerbatim("IsoTrack") << "Track ECAL " << trkDetItr->okECAL << " HCAL " << trkDetItr->okHCAL << " Flag "
                                      << selectTk;
       if (selectTk && trkDetItr->okECAL && trkDetItr->okHCAL) {
         int nRH_eMipDR = 0, nNearTRKs = 0;
         double e1 = spr::eCone_ecal(geo_,
                                     barrelRecHitsHandle_,
                                     endcapRecHitsHandle_,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_neutR1_,
                                     trkDetItr->directionECAL,
                                     nRH_eMipDR);
         double e2 = spr::eCone_ecal(geo_,
                                     barrelRecHitsHandle_,
                                     endcapRecHitsHandle_,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_neutR2_,
                                     trkDetItr->directionECAL,
                                     nRH_eMipDR);
         eMipDR = spr::eCone_ecal(geo_,
                                  barrelRecHitsHandle_,
                                  endcapRecHitsHandle_,
                                  trkDetItr->pointHCAL,
                                  trkDetItr->pointECAL,
                                  a_mipR_,
                                  trkDetItr->directionECAL,
                                  nRH_eMipDR);
         conehmaxNearP =
             spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, ((verbosity_ / 100) % 10 > 1));
         e_inCone = e2 - e1;
         double distFromHotCell = -99.0;
         int nRecHitsCone = -99, ietaHotCell = -99, iphiHotCell = -99;
         GlobalPoint gposHotCell(0., 0., 0.);
         std::vector<DetId> coneRecHitDetIds;
         hCone = spr::eCone_hcal(geo_,
                                 hbhe_,
                                 trkDetItr->pointHCAL,
                                 trkDetItr->pointECAL,
                                 a_coneR_,
                                 trkDetItr->directionHCAL,
                                 nRecHitsCone,
                                 coneRecHitDetIds,
                                 distFromHotCell,
                                 ietaHotCell,
                                 iphiHotCell,
                                 gposHotCell,
                                 -1);
       }
       if (l3Track) {
         fillHist(10, v4);
         if (selectTk) {
           fillHist(11, v4);
           fillCuts(0, eMipDR, conehmaxNearP, e_inCone, v4, ieta, (mindR > dr_L1_));
           if (conehmaxNearP < cutCharge_) {
             fillHist(12, v4);
             if (eMipDR < cutMip_) {
               fillHist(13, v4);
               fillEnergy(1, ieta, hCone, eMipDR, v4);
               fillEnergy(0, ieta, hCone, eMipDR, v4);
               if (e_inCone < cutNeutral_) {
                 fillHist(14, v4);
                 fillEnergy(4, ieta, hCone, eMipDR, v4);
                 fillEnergy(3, ieta, hCone, eMipDR, v4);
               }
             }
           }
         }
       } else if (doL2L3_) {
         fillHist(15, v4);
         if (selectTk) {
           fillHist(16, v4);
           fillCuts(1, eMipDR, conehmaxNearP, e_inCone, v4, ieta, (mindR > dr_L1_));
           if (conehmaxNearP < cutCharge_) {
             fillHist(17, v4);
             if (eMipDR < cutMip_) {
               fillHist(18, v4);
               fillEnergy(2, ieta, hCone, eMipDR, v4);
               fillEnergy(0, ieta, hCone, eMipDR, v4);
               if (e_inCone < cutNeutral_) {
                 fillHist(19, v4);
                 fillEnergy(5, ieta, hCone, eMipDR, v4);
                 fillEnergy(3, ieta, hCone, eMipDR, v4);
               }
             }
           }
         }
       }
     }
   }
 }

◆ studyMipCut()

void IsoTrig::studyMipCut	(	edm::Handle< reco::TrackCollection > &	trkCollection,
		edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &	L2cands
	)

private

Definition at line 1441 of file IsoTrig.cc.

References a_charIsoR_, a_coneR_, a_mipR_, a_neutR1_, a_neutR2_, funct::abs(), barrelRecHitsHandle_, bField_, spr::chargeIsolationCone(), clearMipCutTreeVectors(), cutCharge_, cutMip_, cutNeutral_, l1ctLayer1_cff::dr, dR(), StorageManager_cfg::e1, spr::eCone_ecal(), spr::eCone_hcal(), endcapRecHitsHandle_, reco::TrackBase::eta(), fillDifferences(), geo_, spr::goodTrack(), h_EnIn, h_EnOut, h_MipEnMatch, h_MipEnNoMatch, hbhe_, mps_fire::i, edm::HandleBase::isValid(), leadPV_, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, MipCutTree_, beamSpotPI::nTracks, reco::TrackBase::p(), reco::TrackBase::phi(), spr::propagateCALO(), reco::TrackBase::pt(), pushMipCutTreeVecs(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters_, theTrackQuality_, and verbosity_.

Referenced by analyze().

                                                                                          {
   clearMipCutTreeVectors();
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "inside studymipcut";
   if (!trkCollection.isValid()) {
     edm::LogWarning("IsoTrack") << "trkCollection.isValid is false";
   } else {
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
     std::vector<spr::propagatedTrackDirection> trkCaloDirections1;
     spr::propagateCALO(
         trkCollection, geo_, bField_, theTrackQuality_, trkCaloDirections1, ((verbosity_ / 100) % 10 > 2));
     if (verbosity_ % 10 > 0)
       edm::LogVerbatim("IsoTrack") << "Number of L2cands:" << L2cands->size() << " to be matched to something out of "
                                    << trkCaloDirections1.size() << " reco tracks";
     for (unsigned int i = 0; i < L2cands->size(); i++) {
       edm::Ref<reco::IsolatedPixelTrackCandidateCollection> candref =
           edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(L2cands, i);
       double enIn = candref->energyIn();
       h_EnIn->Fill(candref->energyIn());
       h_EnOut->Fill(candref->energyOut());
       math::XYZTLorentzVector v1(
           candref->track()->px(), candref->track()->py(), candref->track()->pz(), candref->track()->p());
       if (verbosity_ % 10 > 1)
         edm::LogVerbatim("IsoTrack") << "HLT Level Candidate eta/phi/pt/enIn:" << candref->track()->eta() << "/"
                                      << candref->track()->phi() << "/" << candref->track()->pt() << "/"
                                      << candref->energyIn();
       math::XYZTLorentzVector mindRvec;
       double mindR = 999.9, mindP1 = 999.9, eMipDR = 0.0;
       std::vector<bool> selFlags;
       unsigned int nTracks = 0;
       double conehmaxNearP = 0, hCone = 0;
       for (trkDetItr = trkCaloDirections1.begin(); trkDetItr != trkCaloDirections1.end(); trkDetItr++, nTracks++) {
         const reco::Track *pTrack = &(*(trkDetItr->trkItr));
         math::XYZTLorentzVector v2(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
         double dr = dR(v1, v2);
         double dp1 = std::abs(1. / v1.r() - 1. / v2.r());
         if (verbosity_ % 1000 > 0)
           edm::LogVerbatim("IsoTrack") << "This recotrack(eta/phi/pt) " << pTrack->eta() << "/" << pTrack->phi() << "/"
                                        << pTrack->pt() << " has dr/dp= " << dr << "/" << dp1;
         if (dr < mindR) {
           mindR = dr;
           mindP1 = dp1;
           mindRvec = v2;
           int nRH_eMipDR = 0, nNearTRKs = 0;
           eMipDR = spr::eCone_ecal(geo_,
                                    barrelRecHitsHandle_,
                                    endcapRecHitsHandle_,
                                    trkDetItr->pointHCAL,
                                    trkDetItr->pointECAL,
                                    a_mipR_,
                                    trkDetItr->directionECAL,
                                    nRH_eMipDR);
           bool selectTk = spr::goodTrack(pTrack, leadPV_, selectionParameters_, ((verbosity_ / 100) % 10 > 1));
           spr::trackSelectionParameters oneCutParameters = selectionParameters_;
           oneCutParameters.maxDxyPV = 10;
           oneCutParameters.maxDzPV = 100;
           oneCutParameters.maxInMiss = 2;
           oneCutParameters.maxOutMiss = 2;
           bool qltyFlag = spr::goodTrack(pTrack, leadPV_, oneCutParameters, ((verbosity_ / 100) % 10 > 1));
           oneCutParameters = selectionParameters_;
           oneCutParameters.maxDxyPV = 10;
           oneCutParameters.maxDzPV = 100;
           bool qltyMissFlag = spr::goodTrack(pTrack, leadPV_, oneCutParameters, ((verbosity_ / 100) % 10 > 1));
           oneCutParameters = selectionParameters_;
           oneCutParameters.maxInMiss = 2;
           oneCutParameters.maxOutMiss = 2;
           bool qltyPVFlag = spr::goodTrack(pTrack, leadPV_, oneCutParameters, ((verbosity_ / 100) % 10 > 1));
           conehmaxNearP = spr::chargeIsolationCone(
               nTracks, trkCaloDirections1, a_charIsoR_, nNearTRKs, ((verbosity_ / 100) % 10 > 1));
           double e1 = spr::eCone_ecal(geo_,
                                       barrelRecHitsHandle_,
                                       endcapRecHitsHandle_,
                                       trkDetItr->pointHCAL,
                                       trkDetItr->pointECAL,
                                       a_neutR1_,
                                       trkDetItr->directionECAL,
                                       nRH_eMipDR);
           double e2 = spr::eCone_ecal(geo_,
                                       barrelRecHitsHandle_,
                                       endcapRecHitsHandle_,
                                       trkDetItr->pointHCAL,
                                       trkDetItr->pointECAL,
                                       a_neutR2_,
                                       trkDetItr->directionECAL,
                                       nRH_eMipDR);
           double e_inCone = e2 - e1;
           bool chgIsolFlag = (conehmaxNearP < cutCharge_);
           bool mipFlag = (eMipDR < cutMip_);
           bool neuIsolFlag = (e_inCone < cutNeutral_);
           bool trkpropFlag = ((trkDetItr->okECAL) && (trkDetItr->okHCAL));
           selFlags.clear();
           selFlags.push_back(selectTk);
           selFlags.push_back(qltyFlag);
           selFlags.push_back(qltyMissFlag);
           selFlags.push_back(qltyPVFlag);
           selFlags.push_back(trkpropFlag);
           selFlags.push_back(chgIsolFlag);
           selFlags.push_back(neuIsolFlag);
           selFlags.push_back(mipFlag);
           double distFromHotCell = -99.0;
           int nRecHitsCone = -99, ietaHotCell = -99, iphiHotCell = -99;
           GlobalPoint gposHotCell(0., 0., 0.);
           std::vector<DetId> coneRecHitDetIds;
           hCone = spr::eCone_hcal(geo_,
                                   hbhe_,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_coneR_,
                                   trkDetItr->directionHCAL,
                                   nRecHitsCone,
                                   coneRecHitDetIds,
                                   distFromHotCell,
                                   ietaHotCell,
                                   iphiHotCell,
                                   gposHotCell,
                                   -1);
         }
       }
       pushMipCutTreeVecs(v1, mindRvec, enIn, eMipDR, mindR, mindP1, selFlags, hCone);
       fillDifferences(6, v1, mindRvec, (verbosity_ % 10 > 0));
       if (mindR >= 0.05) {
         fillDifferences(8, v1, mindRvec, (verbosity_ % 10 > 0));
         h_MipEnNoMatch->Fill(candref->energyIn(), eMipDR);
       } else {
         fillDifferences(7, v1, mindRvec, (verbosity_ % 10 > 0));
         h_MipEnMatch->Fill(candref->energyIn(), eMipDR);
       }
     }
   }
   MipCutTree_->Fill();
 }

◆ studyTiming()

void IsoTrig::studyTiming ( const edm::Event & theEvent )

private

Definition at line 1330 of file IsoTrig.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, PVValHelper::dxy, PVValHelper::dz, PVValHelper::eta, edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), AlCaHLTBitMon_ParallelJobs::p, phi, pixelTrackRefsHB_, pixelTrackRefsHE_, edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, dttmaxenums::R, SeedingLayerSetsHits::size(), edmNew::DetSetVector< T >::size(), t_nPixCand, t_nPixSeed, t_timeL2Prod, tauUnbiasCone_, TimingTree_, tok_l1cand_, tok_SeedingLayerHB_, tok_SeedingLayerHE_, tok_SiPixelRecHits_, tok_verthb_, tok_verthe_, trigger::TriggerL1CenJet, trigger::TriggerL1ForJet, trigger::TriggerL1TauJet, verbosity_, and vtxCutSeed_.

Referenced by analyze().

                                                   {
   t_timeL2Prod->clear();
   t_nPixCand->clear();
   t_nPixSeed->clear();
 
   if (verbosity_ % 10 > 0) {
     edm::Handle<SeedingLayerSetsHits> hblayers, helayers;
     theEvent.getByToken(tok_SeedingLayerHB_, hblayers);
     theEvent.getByToken(tok_SeedingLayerHE_, helayers);
     const SeedingLayerSetsHits *layershb = hblayers.product();
     const SeedingLayerSetsHits *layershe = helayers.product();
     edm::LogVerbatim("IsoTrack") << "size of Seeding TripletLayers hb/he " << layershb->size() << "/"
                                  << layershe->size();
     edm::Handle<SiPixelRecHitCollection> rchts;
     theEvent.getByToken(tok_SiPixelRecHits_, rchts);
     const SiPixelRecHitCollection *rechits = rchts.product();
     edm::LogVerbatim("IsoTrack") << "size of SiPixelRechits " << rechits->size();
   }
   int nCandHB = pixelTrackRefsHB_.size(), nCandHE = pixelTrackRefsHE_.size();
   int nSeedHB = 0, nSeedHE = 0;
 
   if (nCandHE > 0) {
     edm::Handle<reco::VertexCollection> pVertHB, pVertHE;
     theEvent.getByToken(tok_verthb_, pVertHB);
     theEvent.getByToken(tok_verthe_, pVertHE);
     edm::Handle<trigger::TriggerFilterObjectWithRefs> l1trigobj;
     theEvent.getByToken(tok_l1cand_, l1trigobj);
 
     std::vector<edm::Ref<l1extra::L1JetParticleCollection>> l1tauobjref;
     std::vector<edm::Ref<l1extra::L1JetParticleCollection>> l1jetobjref;
     std::vector<edm::Ref<l1extra::L1JetParticleCollection>> l1forjetobjref;
 
     l1trigobj->getObjects(trigger::TriggerL1TauJet, l1tauobjref);
     l1trigobj->getObjects(trigger::TriggerL1CenJet, l1jetobjref);
     l1trigobj->getObjects(trigger::TriggerL1ForJet, l1forjetobjref);
 
     double ptTriggered = -10;
     double etaTriggered = -100;
     double phiTriggered = -100;
     for (unsigned int p = 0; p < l1tauobjref.size(); p++) {
       if (l1tauobjref[p]->pt() > ptTriggered) {
         ptTriggered = l1tauobjref[p]->pt();
         phiTriggered = l1tauobjref[p]->phi();
         etaTriggered = l1tauobjref[p]->eta();
       }
     }
     for (unsigned int p = 0; p < l1jetobjref.size(); p++) {
       if (l1jetobjref[p]->pt() > ptTriggered) {
         ptTriggered = l1jetobjref[p]->pt();
         phiTriggered = l1jetobjref[p]->phi();
         etaTriggered = l1jetobjref[p]->eta();
       }
     }
     for (unsigned int p = 0; p < l1forjetobjref.size(); p++) {
       if (l1forjetobjref[p]->pt() > ptTriggered) {
         ptTriggered = l1forjetobjref[p]->pt();
         phiTriggered = l1forjetobjref[p]->phi();
         etaTriggered = l1forjetobjref[p]->eta();
       }
     }
     for (unsigned iS = 0; iS < pixelTrackRefsHE_.size(); iS++) {
       reco::VertexCollection::const_iterator vitSel;
       double minDZ = 100;
       bool vtxMatch;
       for (reco::VertexCollection::const_iterator vit = pVertHE->begin(); vit != pVertHE->end(); vit++) {
         if (fabs(pixelTrackRefsHE_[iS]->dz(vit->position())) < minDZ) {
           minDZ = fabs(pixelTrackRefsHE_[iS]->dz(vit->position()));
           vitSel = vit;
         }
       }
       //cut on dYX:
       if ((fabs(pixelTrackRefsHE_[iS]->dxy(vitSel->position())) < vtxCutSeed_) || (minDZ == 100))
         vtxMatch = true;
 
       //select tracks not matched to triggered L1 jet
       double R = deltaR(etaTriggered, phiTriggered, pixelTrackRefsHE_[iS]->eta(), pixelTrackRefsHE_[iS]->phi());
       if (R > tauUnbiasCone_ && vtxMatch)
         nSeedHE++;
     }
     for (unsigned iS = 0; iS < pixelTrackRefsHB_.size(); iS++) {
       reco::VertexCollection::const_iterator vitSel;
       double minDZ = 100;
       bool vtxMatch(false);
       for (reco::VertexCollection::const_iterator vit = pVertHB->begin(); vit != pVertHB->end(); vit++) {
         if (fabs(pixelTrackRefsHB_[iS]->dz(vit->position())) < minDZ) {
           minDZ = fabs(pixelTrackRefsHB_[iS]->dz(vit->position()));
           vitSel = vit;
         }
       }
       //cut on dYX:
       if ((fabs(pixelTrackRefsHB_[iS]->dxy(vitSel->position())) < 101.0) || (minDZ == 100))
         vtxMatch = true;
 
       //select tracks not matched to triggered L1 jet
       double R = deltaR(etaTriggered, phiTriggered, pixelTrackRefsHB_[iS]->eta(), pixelTrackRefsHB_[iS]->phi());
       if (R > 1.2 && vtxMatch)
         nSeedHB++;
     }
 
     if (verbosity_ % 10 > 0) {
       edm::LogVerbatim("IsoTrack") << "(HB/HE) nCand: " << nCandHB << "/" << nCandHE << "nSeed: " << nSeedHB << "/"
                                    << nSeedHE;
     }
   }
   t_nPixSeed->push_back(nSeedHB);
   t_nPixSeed->push_back(nSeedHE);
   t_nPixCand->push_back(nCandHB);
   t_nPixCand->push_back(nCandHE);
 
   TimingTree_->Fill();
 }

◆ studyTrigger()

void IsoTrig::studyTrigger	(	edm::Handle< reco::TrackCollection > &	trkCollection,
		std::vector< reco::TrackCollection::const_iterator > &	goodTks
	)

private

Definition at line 1574 of file IsoTrig.cc.

References funct::abs(), PbPb_ZMuSkimMuonDPG_cff::deltaR, dEta(), Calorimetry_cff::dp, dPhi(), l1ctLayer1_cff::dr, dR(), fillDifferences(), fillHist(), h_dEtaL1, h_dPhiL1, h_dRL1, mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmdumpme::k, mathSSE::sqrt(), vec_, and verbosity_.

Referenced by analyze().

                                                                                   {
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("IsoTrack") << "Inside StudyTrigger";
   for (int j = 0; j < 3; j++) {
     for (unsigned int k = 0; k < vec_[j].size(); k++) {
       if (verbosity_ % 10 > 0)
         edm::LogVerbatim("IsoTrack") << "vec[" << j << "][" << k << "] pt " << vec_[j][k].pt() << " eta "
                                      << vec_[j][k].eta() << " phi " << vec_[j][k].phi();
       fillHist(j, vec_[j][k]);
     }
   }
 
   double deta, dphi, dr;
   for (int lvl = 1; lvl < 3; lvl++) {
     for (unsigned int i = 0; i < vec_[lvl].size(); i++) {
       deta = dEta(vec_[0][0], vec_[lvl][i]);
       dphi = dPhi(vec_[0][0], vec_[lvl][i]);
       dr = dR(vec_[0][0], vec_[lvl][i]);
       if (verbosity_ % 10 > 1)
         edm::LogVerbatim("IsoTrack") << "lvl " << lvl << " i " << i << " deta " << deta << " dphi " << dphi << " dR "
                                      << dr;
       h_dEtaL1[lvl - 1]->Fill(deta);
       h_dPhiL1[lvl - 1]->Fill(dphi);
       h_dRL1[lvl - 1]->Fill(std::sqrt(dr));
     }
   }
 
   math::XYZTLorentzVector mindRvec;
   double mindR;
   for (unsigned int k = 0; k < vec_[2].size(); ++k) {
     mindR = 999.9;
     if (verbosity_ % 10 > 1)
       edm::LogVerbatim("IsoTrack") << "L3obj: pt " << vec_[2][k].pt() << " eta " << vec_[2][k].eta() << " phi "
                                    << vec_[2][k].phi();
     for (unsigned int j = 0; j < vec_[1].size(); j++) {
       dr = dR(vec_[2][k], vec_[1][j]);
       if (dr < mindR) {
         mindR = dr;
         mindRvec = vec_[1][j];
       }
     }
     fillDifferences(0, vec_[2][k], mindRvec, (verbosity_ % 10 > 0));
     if (mindR < 0.03) {
       fillDifferences(1, vec_[2][k], mindRvec, (verbosity_ % 10 > 0));
       fillHist(6, mindRvec);
       fillHist(8, vec_[2][k]);
     } else {
       fillDifferences(2, vec_[2][k], mindRvec, (verbosity_ % 10 > 0));
       fillHist(7, mindRvec);
       fillHist(9, vec_[2][k]);
     }
 
     mindR = 999.9;
     if (verbosity_ % 10 > 0)
       edm::LogVerbatim("IsoTrack") << "Now Matching L3 track with reco: L3 Track (eta, phi) " << vec_[2][k].eta() << ":"
                                    << vec_[2][k].phi() << " L2 Track " << mindRvec.eta() << ":" << mindRvec.phi()
                                    << " dR " << mindR;
     reco::TrackCollection::const_iterator goodTk = trkCollection->end();
     if (trkCollection.isValid()) {
       double mindP(9999.9);
       reco::TrackCollection::const_iterator trkItr;
       for (trkItr = trkCollection->begin(); trkItr != trkCollection->end(); trkItr++) {
         math::XYZTLorentzVector v4(trkItr->px(), trkItr->py(), trkItr->pz(), trkItr->p());
         double deltaR = dR(v4, vec_[2][k]);
         double dp = std::abs(v4.r() / vec_[2][k].r() - 1.0);
         if (deltaR < mindR) {
           mindR = deltaR;
           mindP = dp;
           mindRvec = v4;
           goodTk = trkItr;
         }
         if ((verbosity_ / 10) % 10 > 1 && deltaR < 1.0)
           edm::LogVerbatim("IsoTrack") << "reco track: pt " << v4.pt() << " eta " << v4.eta() << " phi " << v4.phi()
                                        << " DR " << deltaR;
       }
       if (verbosity_ % 10 > 0)
         edm::LogVerbatim("IsoTrack") << "Now Matching at Reco level in step 1 DR: " << mindR << ":" << mindP
                                      << " eta:phi " << mindRvec.eta() << ":" << mindRvec.phi();
       if (mindR < 0.03 && mindP > 0.1) {
         for (trkItr = trkCollection->begin(); trkItr != trkCollection->end(); trkItr++) {
           math::XYZTLorentzVector v4(trkItr->px(), trkItr->py(), trkItr->pz(), trkItr->p());
           double deltaR = dR(v4, vec_[2][k]);
           double dp = std::abs(v4.r() / vec_[2][k].r() - 1.0);
           if (dp < mindP && deltaR < 0.03) {
             mindR = deltaR;
             mindP = dp;
             mindRvec = v4;
             goodTk = trkItr;
           }
         }
         if (verbosity_ % 10 > 0)
           edm::LogVerbatim("IsoTrack") << "Now Matching at Reco level in step 2 DR: " << mindR << ":" << mindP
                                        << " eta:phi " << mindRvec.eta() << ":" << mindRvec.phi();
       }
       fillDifferences(3, vec_[2][k], mindRvec, (verbosity_ % 10 > 0));
       fillHist(3, mindRvec);
       if (mindR < 0.03) {
         fillDifferences(4, vec_[2][k], mindRvec, (verbosity_ % 10 > 0));
         fillHist(4, mindRvec);
       } else {
         fillDifferences(5, vec_[2][k], mindRvec, (verbosity_ % 10 > 0));
         fillHist(5, mindRvec);
       }
       if (goodTk != trkCollection->end())
         goodTks.push_back(goodTk);
     }
   }
 }

◆ StudyTrkEbyP()

void IsoTrig::StudyTrkEbyP ( edm::Handle< reco::TrackCollection > & trkCollection )

private

Definition at line 1206 of file IsoTrig.cc.

References a_charIsoR_, a_coneR_, a_mipR_, a_neutR1_, a_neutR2_, barrelRecHitsHandle_, bField_, spr::chargeIsolationCone(), cutCharge_, cutMip_, cutNeutral_, StorageManager_cfg::e1, spr::eCone_ecal(), spr::eCone_hcal(), endcapRecHitsHandle_, geo_, spr::goodTrack(), hbhe_, edm::HandleBase::isValid(), leadPV_, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, beamSpotPI::nTracks, reco::TrackBase::p(), spr::propagateCALO(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters_, t_TrkhCone, t_TrkMissFlag, t_TrkNuIsolFlag, t_TrkP, t_TrkPVFlag, t_TrkqltyFlag, t_TrkselTkFlag, theTrackQuality_, TrkResTree_, and verbosity_.

Referenced by analyze().

                                                                         {
   t_TrkselTkFlag->clear();
   t_TrkqltyFlag->clear();
   t_TrkMissFlag->clear();
   t_TrkPVFlag->clear();
   t_TrkNuIsolFlag->clear();
   t_TrkhCone->clear();
   t_TrkP->clear();
 
   if (!trkCollection.isValid()) {
     edm::LogVerbatim("IsoTrack") << "trkCollection.isValid is false";
   } else {
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
     std::vector<spr::propagatedTrackDirection> trkCaloDirections1;
     spr::propagateCALO(
         trkCollection, geo_, bField_, theTrackQuality_, trkCaloDirections1, ((verbosity_ / 100) % 10 > 2));
     unsigned int nTracks = 0;
     int nRH_eMipDR = 0, nNearTRKs = 0;
     std::vector<bool> selFlags;
     for (trkDetItr = trkCaloDirections1.begin(); trkDetItr != trkCaloDirections1.end(); trkDetItr++, nTracks++) {
       double conehmaxNearP = 0, hCone = 0, eMipDR = 0.0;
       const reco::Track *pTrack = &(*(trkDetItr->trkItr));
       if (verbosity_ % 10 > 0)
         edm::LogVerbatim("IsoTrack") << "track no. " << nTracks << " p(): " << pTrack->p();
       if (pTrack->p() > 20) {
         math::XYZTLorentzVector v2(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
         eMipDR = spr::eCone_ecal(geo_,
                                  barrelRecHitsHandle_,
                                  endcapRecHitsHandle_,
                                  trkDetItr->pointHCAL,
                                  trkDetItr->pointECAL,
                                  a_mipR_,
                                  trkDetItr->directionECAL,
                                  nRH_eMipDR);
         bool selectTk = spr::goodTrack(pTrack, leadPV_, selectionParameters_, ((verbosity_ / 100) % 10 > 1));
         spr::trackSelectionParameters oneCutParameters = selectionParameters_;
         oneCutParameters.maxDxyPV = 10;
         oneCutParameters.maxDzPV = 100;
         oneCutParameters.maxInMiss = 2;
         oneCutParameters.maxOutMiss = 2;
         bool qltyFlag = spr::goodTrack(pTrack, leadPV_, oneCutParameters, ((verbosity_ / 100) % 10 > 1));
         oneCutParameters = selectionParameters_;
         oneCutParameters.maxDxyPV = 10;
         oneCutParameters.maxDzPV = 100;
         bool qltyMissFlag = spr::goodTrack(pTrack, leadPV_, oneCutParameters, ((verbosity_ / 100) % 10 > 1));
         oneCutParameters = selectionParameters_;
         oneCutParameters.maxInMiss = 2;
         oneCutParameters.maxOutMiss = 2;
         bool qltyPVFlag = spr::goodTrack(pTrack, leadPV_, oneCutParameters, ((verbosity_ / 100) % 10 > 1));
         if ((verbosity_ / 1000) % 10 > 1)
           edm::LogVerbatim("IsoTrack") << "sel " << selectTk << "ntracks " << nTracks << " a_charIsoR " << a_charIsoR_
                                        << " nNearTRKs " << nNearTRKs;
         conehmaxNearP = spr::chargeIsolationCone(
             nTracks, trkCaloDirections1, a_charIsoR_, nNearTRKs, ((verbosity_ / 100) % 10 > 1));
         if ((verbosity_ / 1000) % 10 > 1)
           edm::LogVerbatim("IsoTrack") << "coneh " << conehmaxNearP << "ok " << trkDetItr->okECAL << " "
                                        << trkDetItr->okHCAL;
         double e1 = spr::eCone_ecal(geo_,
                                     barrelRecHitsHandle_,
                                     endcapRecHitsHandle_,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_neutR1_,
                                     trkDetItr->directionECAL,
                                     nRH_eMipDR);
         double e2 = spr::eCone_ecal(geo_,
                                     barrelRecHitsHandle_,
                                     endcapRecHitsHandle_,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_neutR2_,
                                     trkDetItr->directionECAL,
                                     nRH_eMipDR);
         double e_inCone = e2 - e1;
         bool chgIsolFlag = (conehmaxNearP < cutCharge_);
         bool mipFlag = (eMipDR < cutMip_);
         bool neuIsolFlag = (e_inCone < cutNeutral_);
         bool trkpropFlag = ((trkDetItr->okECAL) && (trkDetItr->okHCAL));
         selFlags.clear();
         selFlags.push_back(selectTk);
         selFlags.push_back(qltyFlag);
         selFlags.push_back(qltyMissFlag);
         selFlags.push_back(qltyPVFlag);
         if (verbosity_ % 10 > 0)
           edm::LogVerbatim("IsoTrack") << "emip: " << eMipDR << "<" << cutMip_ << "(" << mipFlag << ")"
                                        << " ; ok: " << trkDetItr->okECAL << "/" << trkDetItr->okHCAL
                                        << " ; chgiso: " << conehmaxNearP << "<" << cutCharge_ << "(" << chgIsolFlag
                                        << ")";
 
         if (chgIsolFlag && mipFlag && trkpropFlag) {
           double distFromHotCell = -99.0;
           int nRecHitsCone = -99, ietaHotCell = -99, iphiHotCell = -99;
           GlobalPoint gposHotCell(0., 0., 0.);
           std::vector<DetId> coneRecHitDetIds;
           hCone = spr::eCone_hcal(geo_,
                                   hbhe_,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_coneR_,
                                   trkDetItr->directionHCAL,
                                   nRecHitsCone,
                                   coneRecHitDetIds,
                                   distFromHotCell,
                                   ietaHotCell,
                                   iphiHotCell,
                                   gposHotCell,
                                   -1);
           // push vectors into the Tree
           t_TrkselTkFlag->push_back(selFlags[0]);
           t_TrkqltyFlag->push_back(selFlags[1]);
           t_TrkMissFlag->push_back(selFlags[2]);
           t_TrkPVFlag->push_back(selFlags[3]);
           t_TrkNuIsolFlag->push_back(neuIsolFlag);
           t_TrkhCone->push_back(hCone);
           t_TrkP->push_back(pTrack->p());
         }
       }
     }
     if (verbosity_ % 10 > 0)
       edm::LogVerbatim("IsoTrack") << "Filling " << t_TrkP->size() << " tracks in TrkRestree out of " << nTracks;
   }
   TrkResTree_->Fill();
 }

Member Data Documentation

◆ a_charIsoR_

const double IsoTrig::a_charIsoR_

private

Definition at line 165 of file IsoTrig.cc.

Referenced by chgIsolation(), getGoodTracks(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ a_coneR_

const double IsoTrig::a_coneR_

private

Definition at line 165 of file IsoTrig.cc.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ a_mipR_

const double IsoTrig::a_mipR_

private

Definition at line 166 of file IsoTrig.cc.

Referenced by getGoodTracks(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ a_neutIsoR_

const double IsoTrig::a_neutIsoR_

private

Definition at line 165 of file IsoTrig.cc.

Referenced by IsoTrig().

◆ a_neutR1_

const double IsoTrig::a_neutR1_

private

Definition at line 166 of file IsoTrig.cc.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ a_neutR2_

const double IsoTrig::a_neutR2_

private

Definition at line 166 of file IsoTrig.cc.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ barrelRecHitsHandle_

edm::Handle<EcalRecHitCollection> IsoTrig::barrelRecHitsHandle_

private

Definition at line 193 of file IsoTrig.cc.

Referenced by analyze(), getGoodTracks(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ beamSpotH_

edm::Handle<reco::BeamSpot> IsoTrig::beamSpotH_

private

Definition at line 195 of file IsoTrig.cc.

Referenced by analyze().

◆ bField_

const MagneticField* IsoTrig::bField_

private

Definition at line 198 of file IsoTrig.cc.

Referenced by analyze(), chgIsolation(), getGoodTracks(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ bfVal_

double IsoTrig::bfVal_

private

Definition at line 163 of file IsoTrig.cc.

Referenced by analyze(), and GetEtaPhiAtEcal().

◆ changed_

bool IsoTrig::changed_

private

Definition at line 204 of file IsoTrig.cc.

Referenced by analyze(), and beginRun().

◆ ChgIsolnTree_

TTree * IsoTrig::ChgIsolnTree_

private

Definition at line 207 of file IsoTrig.cc.

Referenced by beginJob(), and chgIsolation().

◆ cutCharge_

const double IsoTrig::cutCharge_

private

Definition at line 167 of file IsoTrig.cc.

Referenced by fillCuts(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ cutMip_

const double IsoTrig::cutMip_

private

Definition at line 166 of file IsoTrig.cc.

Referenced by fillCuts(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ cutNeutral_

const double IsoTrig::cutNeutral_

private

Definition at line 167 of file IsoTrig.cc.

Referenced by fillCuts(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ doChgIsolTree_

const bool IsoTrig::doChgIsolTree_

private

Definition at line 156 of file IsoTrig.cc.

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

◆ doL2L3_

const bool IsoTrig::doL2L3_

private

Definition at line 155 of file IsoTrig.cc.

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

◆ doMipCutTree_

const bool IsoTrig::doMipCutTree_

private

Definition at line 155 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ doStudyIsol_

const bool IsoTrig::doStudyIsol_

private

Definition at line 156 of file IsoTrig.cc.

Referenced by analyze().

◆ doTiming_

const bool IsoTrig::doTiming_

private

Definition at line 155 of file IsoTrig.cc.

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

◆ doTrkResTree_

const bool IsoTrig::doTrkResTree_

private

Definition at line 156 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ dr_L1_

const double IsoTrig::dr_L1_

private

Definition at line 165 of file IsoTrig.cc.

Referenced by studyIsolation().

◆ endcapRecHitsHandle_

edm::Handle<EcalRecHitCollection> IsoTrig::endcapRecHitsHandle_

private

Definition at line 194 of file IsoTrig.cc.

Referenced by analyze(), getGoodTracks(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ fs_

edm::Service<TFileService> IsoTrig::fs_

private

Definition at line 206 of file IsoTrig.cc.

Referenced by beginJob(), and endJob().

◆ g_Accepts

TH1I* IsoTrig::g_Accepts

private

Definition at line 268 of file IsoTrig.cc.

Referenced by endJob().

◆ g_Pre

TH1D* IsoTrig::g_Pre

private

Definition at line 267 of file IsoTrig.cc.

Referenced by endJob().

◆ g_PreHLT

TH1D * IsoTrig::g_PreHLT

private

Definition at line 267 of file IsoTrig.cc.

Referenced by endJob().

◆ g_PreL1

TH1D * IsoTrig::g_PreL1

private

Definition at line 267 of file IsoTrig.cc.

Referenced by endJob().

◆ geo_

const CaloGeometry* IsoTrig::geo_

private

Definition at line 199 of file IsoTrig.cc.

Referenced by analyze(), chgIsolation(), getGoodTracks(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ h_dEta

TH1D* IsoTrig::h_dEta[9]

private

Definition at line 262 of file IsoTrig.cc.

Referenced by beginJob(), and fillDifferences().

◆ h_dEtaL1

TH1D* IsoTrig::h_dEtaL1[2]

private

Definition at line 261 of file IsoTrig.cc.

Referenced by beginJob(), and studyTrigger().

◆ h_dinvPt

TH1D* IsoTrig::h_dinvPt[9]

private

Definition at line 263 of file IsoTrig.cc.

Referenced by beginJob(), and fillDifferences().

◆ h_dP

TH1D * IsoTrig::h_dP[9]

private

Definition at line 262 of file IsoTrig.cc.

Referenced by beginJob(), and fillDifferences().

◆ h_dPhi

TH1D * IsoTrig::h_dPhi[9]

private

Definition at line 262 of file IsoTrig.cc.

Referenced by beginJob(), and fillDifferences().

◆ h_dPhiL1

TH1D * IsoTrig::h_dPhiL1[2]

private

Definition at line 261 of file IsoTrig.cc.

Referenced by beginJob(), and studyTrigger().

◆ h_dPt

TH1D * IsoTrig::h_dPt[9]

private

Definition at line 262 of file IsoTrig.cc.

Referenced by beginJob(), and fillDifferences().

◆ h_dRL1

TH1D * IsoTrig::h_dRL1[2]

private

Definition at line 261 of file IsoTrig.cc.

Referenced by beginJob(), and studyTrigger().

◆ h_eCalo

TH1D * IsoTrig::h_eCalo[5][6][48]

private

Definition at line 266 of file IsoTrig.cc.

Referenced by beginJob(), and fillEnergy().

◆ h_eHcal

TH1D* IsoTrig::h_eHcal[5][6][48]

private

Definition at line 266 of file IsoTrig.cc.

Referenced by beginJob(), and fillEnergy().

◆ h_eMaxNearP

TH1D* IsoTrig::h_eMaxNearP[2]

private

Definition at line 264 of file IsoTrig.cc.

Referenced by beginJob(), and fillCuts().

◆ h_eMip

TH1D * IsoTrig::h_eMip[2]

private

Definition at line 263 of file IsoTrig.cc.

Referenced by beginJob(), and fillCuts().

◆ h_eNeutIso

TH1D * IsoTrig::h_eNeutIso[2]

private

Definition at line 264 of file IsoTrig.cc.

Referenced by beginJob(), and fillCuts().

◆ h_EnIn

TH1D* IsoTrig::h_EnIn

private

Definition at line 255 of file IsoTrig.cc.

Referenced by beginJob(), and studyMipCut().

◆ h_EnOut

TH1D * IsoTrig::h_EnOut

private

Definition at line 255 of file IsoTrig.cc.

Referenced by beginJob(), and studyMipCut().

◆ h_eta

TH1D * IsoTrig::h_eta[20]

private

Definition at line 260 of file IsoTrig.cc.

Referenced by beginJob(), and fillHist().

◆ h_etaCalibTracks

TH1D* IsoTrig::h_etaCalibTracks[5][2][2]

private

Definition at line 265 of file IsoTrig.cc.

Referenced by beginJob(), and fillCuts().

◆ h_etaMipTracks

TH1D * IsoTrig::h_etaMipTracks[5][2][2]

private

Definition at line 265 of file IsoTrig.cc.

Referenced by beginJob(), and fillCuts().

◆ h_Filters

TH1I * IsoTrig::h_Filters

private

Definition at line 258 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_HLT

TH1I * IsoTrig::h_HLT

private

Definition at line 257 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_L1ObjEnergy

TH1D * IsoTrig::h_L1ObjEnergy

private

Definition at line 259 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_mindR

TH1D * IsoTrig::h_mindR[9]

private

Definition at line 263 of file IsoTrig.cc.

Referenced by beginJob(), and fillDifferences().

◆ h_MipEnMatch

TH2D* IsoTrig::h_MipEnMatch

private

Definition at line 256 of file IsoTrig.cc.

Referenced by beginJob(), and studyMipCut().

◆ h_MipEnNoMatch

TH2D * IsoTrig::h_MipEnNoMatch

private

Definition at line 256 of file IsoTrig.cc.

Referenced by beginJob(), and studyMipCut().

◆ h_nHLT

TH1I* IsoTrig::h_nHLT

private

Definition at line 257 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_nL3Objs

TH1I * IsoTrig::h_nL3Objs

private

Definition at line 258 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_p

TH1D* IsoTrig::h_p[20]

private

Definition at line 260 of file IsoTrig.cc.

Referenced by beginJob(), and fillHist().

◆ h_phi

TH1D * IsoTrig::h_phi[20]

private

Definition at line 260 of file IsoTrig.cc.

Referenced by beginJob(), and fillHist().

◆ h_Pre

TH1I* IsoTrig::h_Pre

private

Definition at line 258 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_PreHLT

TH1I * IsoTrig::h_PreHLT

private

Definition at line 257 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_PreHLTwt

TH1D * IsoTrig::h_PreHLTwt

private

Definition at line 259 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_PreL1

TH1I * IsoTrig::h_PreL1

private

Definition at line 257 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_PreL1wt

TH1D* IsoTrig::h_PreL1wt

private

Definition at line 259 of file IsoTrig.cc.

Referenced by analyze(), and beginJob().

◆ h_pt

TH1D * IsoTrig::h_pt[20]

private

Definition at line 260 of file IsoTrig.cc.

Referenced by beginJob(), and fillHist().

◆ hbhe_

edm::Handle<HBHERecHitCollection> IsoTrig::hbhe_

private

Definition at line 192 of file IsoTrig.cc.

Referenced by analyze(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ hltPrescaleProvider_

HLTPrescaleProvider IsoTrig::hltPrescaleProvider_

private

Definition at line 151 of file IsoTrig.cc.

Referenced by analyze(), and beginRun().

◆ l1CandTag_

const edm::InputTag IsoTrig::l1CandTag_

private

Definition at line 153 of file IsoTrig.cc.

Referenced by IsoTrig().

◆ l2CandTag_

const edm::InputTag IsoTrig::l2CandTag_

private

Definition at line 153 of file IsoTrig.cc.

Referenced by IsoTrig().

◆ leadPV_

math::XYZPoint IsoTrig::leadPV_

private

Definition at line 200 of file IsoTrig.cc.

Referenced by analyze(), chgIsolation(), getGoodTracks(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ maxRunNo_

const int IsoTrig::maxRunNo_

private

Definition at line 168 of file IsoTrig.cc.

Referenced by endJob().

◆ minPTrackValue_

const double IsoTrig::minPTrackValue_

private

Definition at line 159 of file IsoTrig.cc.

Referenced by chgIsolation(), and IsoTrig().

◆ minRunNo_

const int IsoTrig::minRunNo_

private

Definition at line 168 of file IsoTrig.cc.

Referenced by endJob().

◆ MipCutTree_

TTree* IsoTrig::MipCutTree_

private

Definition at line 207 of file IsoTrig.cc.

Referenced by beginJob(), and studyMipCut().

◆ pixCandTag_

const edm::InputTag IsoTrig::pixCandTag_

private

Definition at line 153 of file IsoTrig.cc.

Referenced by IsoTrig().

◆ pixelIsolationConeSizeAtEC_

const std::vector<double> IsoTrig::pixelIsolationConeSizeAtEC_

private

Definition at line 158 of file IsoTrig.cc.

Referenced by chgIsolation(), and IsoTrig().

◆ pixelTrackRefsHB_

std::vector<reco::TrackRef> IsoTrig::pixelTrackRefsHB_

private

Definition at line 191 of file IsoTrig.cc.

Referenced by analyze(), and studyTiming().

◆ pixelTrackRefsHE_

std::vector<reco::TrackRef> IsoTrig::pixelTrackRefsHE_

private

Definition at line 191 of file IsoTrig.cc.

Referenced by analyze(), chgIsolation(), and studyTiming().

◆ pixelTracksSources_

const std::vector<edm::InputTag> IsoTrig::pixelTracksSources_

private

Definition at line 154 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ pLimits_

double IsoTrig::pLimits_[6]

private

Definition at line 205 of file IsoTrig.cc.

Referenced by beginJob(), fillCuts(), fillEnergy(), and IsoTrig().

◆ prelimCone_

const double IsoTrig::prelimCone_

private

Definition at line 160 of file IsoTrig.cc.

Referenced by chgIsolation(), and IsoTrig().

◆ processName_

const std::string IsoTrig::processName_

private

Definition at line 162 of file IsoTrig.cc.

Referenced by beginRun(), and IsoTrig().

◆ rEB_

double IsoTrig::rEB_

private

Definition at line 163 of file IsoTrig.cc.

Referenced by getDistInCM(), GetEtaPhiAtEcal(), and IsoTrig().

◆ recVtxs_

edm::Handle<reco::VertexCollection> IsoTrig::recVtxs_

private

Definition at line 196 of file IsoTrig.cc.

Referenced by analyze(), and chgIsolation().

◆ selectionParameters_

spr::trackSelectionParameters IsoTrig::selectionParameters_

private

Definition at line 164 of file IsoTrig.cc.

Referenced by chgIsolation(), getGoodTracks(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ t_ECone

std::vector<double>* IsoTrig::t_ECone

private

Definition at line 253 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFcandEmip

std::vector<double>* IsoTrig::t_NFcandEmip

private

Definition at line 237 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFcandEta

std::vector<double>* IsoTrig::t_NFcandEta

private

Definition at line 235 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFcandP

std::vector<double>* IsoTrig::t_NFcandP

private

Definition at line 233 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFcandPhi

std::vector<double>* IsoTrig::t_NFcandPhi

private

Definition at line 236 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFcandPt

std::vector<double>* IsoTrig::t_NFcandPt

private

Definition at line 234 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkcandEmip

std::vector<double>* IsoTrig::t_NFTrkcandEmip

private

Definition at line 242 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkcandEta

std::vector<double>* IsoTrig::t_NFTrkcandEta

private

Definition at line 240 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkcandP

std::vector<double>* IsoTrig::t_NFTrkcandP

private

Definition at line 238 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkcandPhi

std::vector<double>* IsoTrig::t_NFTrkcandPhi

private

Definition at line 241 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkcandPt

std::vector<double>* IsoTrig::t_NFTrkcandPt

private

Definition at line 239 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkChgIsoFlag

std::vector<bool>* IsoTrig::t_NFTrkChgIsoFlag

private

Definition at line 250 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkMinDP1

std::vector<double>* IsoTrig::t_NFTrkMinDP1

private

Definition at line 244 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkMinDR

std::vector<double>* IsoTrig::t_NFTrkMinDR

private

Definition at line 243 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkMipFlag

std::vector<bool>* IsoTrig::t_NFTrkMipFlag

private

Definition at line 252 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkMissFlag

std::vector<bool>* IsoTrig::t_NFTrkMissFlag

private

Definition at line 247 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkNeuIsoFlag

std::vector<bool>* IsoTrig::t_NFTrkNeuIsoFlag

private

Definition at line 251 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkPropFlag

std::vector<bool>* IsoTrig::t_NFTrkPropFlag

private

Definition at line 249 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkPVFlag

std::vector<bool>* IsoTrig::t_NFTrkPVFlag

private

Definition at line 248 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkqltyFlag

std::vector<bool>* IsoTrig::t_NFTrkqltyFlag

private

Definition at line 246 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_NFTrkselTkFlag

std::vector<bool>* IsoTrig::t_NFTrkselTkFlag

private

Definition at line 245 of file IsoTrig.cc.

Referenced by beginJob(), clearMipCutTreeVectors(), pushMipCutTreeVecs(), and ~IsoTrig().

◆ t_nGoodTk

std::vector<int>* IsoTrig::t_nGoodTk

private

Definition at line 211 of file IsoTrig.cc.

Referenced by beginJob(), getGoodTracks(), and ~IsoTrig().

◆ t_nPixCand

std::vector<int>* IsoTrig::t_nPixCand

private

Definition at line 209 of file IsoTrig.cc.

Referenced by beginJob(), studyTiming(), and ~IsoTrig().

◆ t_nPixSeed

std::vector<int>* IsoTrig::t_nPixSeed

private

Definition at line 210 of file IsoTrig.cc.

Referenced by beginJob(), studyTiming(), and ~IsoTrig().

◆ t_PixcandEta

std::vector<double>* IsoTrig::t_PixcandEta

private

Definition at line 223 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixcandMaxP

std::vector<std::vector<double> >* IsoTrig::t_PixcandMaxP

private

Definition at line 225 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixcandP

std::vector<double>* IsoTrig::t_PixcandP

private

Definition at line 221 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixcandPhi

std::vector<double>* IsoTrig::t_PixcandPhi

private

Definition at line 224 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixcandPt

std::vector<double>* IsoTrig::t_PixcandPt

private

Definition at line 222 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixTrkcandEta

std::vector<double>* IsoTrig::t_PixTrkcandEta

private

Definition at line 228 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixTrkcandMaxP

std::vector<double>* IsoTrig::t_PixTrkcandMaxP

private

Definition at line 230 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixTrkcandP

std::vector<double>* IsoTrig::t_PixTrkcandP

private

Definition at line 226 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixTrkcandPhi

std::vector<double>* IsoTrig::t_PixTrkcandPhi

private

Definition at line 229 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixTrkcandPt

std::vector<double>* IsoTrig::t_PixTrkcandPt

private

Definition at line 227 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_PixTrkcandselTk

std::vector<bool>* IsoTrig::t_PixTrkcandselTk

private

Definition at line 231 of file IsoTrig.cc.

Referenced by beginJob(), clearChgIsolnTreeVectors(), pushChgIsolnTreeVecs(), and ~IsoTrig().

◆ t_timeL2Prod

std::vector<double>* IsoTrig::t_timeL2Prod

private

Definition at line 208 of file IsoTrig.cc.

Referenced by beginJob(), studyTiming(), and ~IsoTrig().

◆ t_TrkhCone

std::vector<double>* IsoTrig::t_TrkhCone

private

Definition at line 213 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ t_TrkMissFlag

std::vector<bool>* IsoTrig::t_TrkMissFlag

private

Definition at line 217 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ t_TrkNuIsolFlag

std::vector<bool>* IsoTrig::t_TrkNuIsolFlag

private

Definition at line 219 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ t_TrkP

std::vector<double>* IsoTrig::t_TrkP

private

Definition at line 214 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ t_TrkPVFlag

std::vector<bool>* IsoTrig::t_TrkPVFlag

private

Definition at line 218 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ t_TrkqltyFlag

std::vector<bool>* IsoTrig::t_TrkqltyFlag

private

Definition at line 216 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ t_TrkselTkFlag

std::vector<bool>* IsoTrig::t_TrkselTkFlag

private

Definition at line 215 of file IsoTrig.cc.

Referenced by beginJob(), StudyTrkEbyP(), and ~IsoTrig().

◆ tauUnbiasCone_

const double IsoTrig::tauUnbiasCone_

private

Definition at line 160 of file IsoTrig.cc.

Referenced by chgIsolation(), IsoTrig(), and studyTiming().

◆ theTrackQuality_

std::string IsoTrig::theTrackQuality_

private

Definition at line 161 of file IsoTrig.cc.

Referenced by chgIsolation(), getGoodTracks(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

◆ TimingTree_

TTree * IsoTrig::TimingTree_

private

Definition at line 207 of file IsoTrig.cc.

Referenced by beginJob(), and studyTiming().

◆ tok_bs_

edm::EDGetTokenT<reco::BeamSpot> IsoTrig::tok_bs_

private

Definition at line 175 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> IsoTrig::tok_EB_

private

Definition at line 176 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> IsoTrig::tok_EE_

private

Definition at line 177 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_genTrack_

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

private

Definition at line 173 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_geom_

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

private

Definition at line 188 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> IsoTrig::tok_hbhe_

private

Definition at line 178 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_hlt_

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> IsoTrig::tok_hlt_

private

Definition at line 172 of file IsoTrig.cc.

Referenced by IsoTrig().

◆ tok_l1cand_

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> IsoTrig::tok_l1cand_

private

Definition at line 184 of file IsoTrig.cc.

Referenced by analyze(), getGoodTracks(), IsoTrig(), and studyTiming().

◆ tok_l2cand_

edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> IsoTrig::tok_l2cand_

private

Definition at line 185 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_lumi_

edm::EDGetTokenT<LumiDetails> IsoTrig::tok_lumi_

private

Definition at line 169 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_magField_

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

private

Definition at line 189 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_pixtk_

edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> IsoTrig::tok_pixtk_

private

Definition at line 183 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_pixtks_

std::vector<edm::EDGetTokenT<reco::TrackCollection> > IsoTrig::tok_pixtks_

private

Definition at line 186 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_recVtx_

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

private

Definition at line 174 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_SeedingLayerHB_

edm::EDGetTokenT<SeedingLayerSetsHits> IsoTrig::tok_SeedingLayerHB_

private

Definition at line 180 of file IsoTrig.cc.

Referenced by IsoTrig(), and studyTiming().

◆ tok_SeedingLayerHE_

edm::EDGetTokenT<SeedingLayerSetsHits> IsoTrig::tok_SeedingLayerHE_

private

Definition at line 181 of file IsoTrig.cc.

Referenced by IsoTrig(), and studyTiming().

◆ tok_SiPixelRecHits_

edm::EDGetTokenT<SiPixelRecHitCollection> IsoTrig::tok_SiPixelRecHits_

private

Definition at line 182 of file IsoTrig.cc.

Referenced by IsoTrig(), and studyTiming().

◆ tok_trigEvt_

edm::EDGetTokenT<trigger::TriggerEvent> IsoTrig::tok_trigEvt_

private

Definition at line 170 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_trigRes_

edm::EDGetTokenT<edm::TriggerResults> IsoTrig::tok_trigRes_

private

Definition at line 171 of file IsoTrig.cc.

Referenced by analyze(), and IsoTrig().

◆ tok_verthb_

edm::EDGetTokenT<reco::VertexCollection> IsoTrig::tok_verthb_

private

Definition at line 179 of file IsoTrig.cc.

Referenced by IsoTrig(), and studyTiming().

◆ tok_verthe_

edm::EDGetTokenT<reco::VertexCollection> IsoTrig::tok_verthe_

private

Definition at line 179 of file IsoTrig.cc.

Referenced by IsoTrig(), and studyTiming().

◆ trigList_

std::unordered_map<unsigned int, unsigned int> IsoTrig::trigList_

private

Definition at line 202 of file IsoTrig.cc.

Referenced by analyze(), and endJob().

◆ trigNames_

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

private

Definition at line 152 of file IsoTrig.cc.

Referenced by analyze(), and endJob().

◆ trigPreList_

std::unordered_map<unsigned int, const std::pair<double, double> > IsoTrig::trigPreList_

private

Definition at line 203 of file IsoTrig.cc.

Referenced by analyze(), and endJob().

◆ TrkResTree_

TTree * IsoTrig::TrkResTree_

private

Definition at line 207 of file IsoTrig.cc.

Referenced by beginJob(), and StudyTrkEbyP().

◆ vec_

std::vector<math::XYZTLorentzVector> IsoTrig::vec_[3]

private

Definition at line 269 of file IsoTrig.cc.

Referenced by analyze(), etaPhiTrigger(), studyIsolation(), and studyTrigger().

◆ verbosity_

const int IsoTrig::verbosity_

private

Definition at line 157 of file IsoTrig.cc.

Referenced by analyze(), chgIsolation(), getGoodTracks(), IsoTrig(), studyIsolation(), studyMipCut(), studyTiming(), studyTrigger(), and StudyTrkEbyP().

◆ vtxCutIsol_

const double IsoTrig::vtxCutIsol_

private

Definition at line 159 of file IsoTrig.cc.

Referenced by chgIsolation(), and IsoTrig().

◆ vtxCutSeed_

const double IsoTrig::vtxCutSeed_

private

Definition at line 159 of file IsoTrig.cc.

Referenced by chgIsolation(), IsoTrig(), and studyTiming().

◆ zEE_

double IsoTrig::zEE_

private

Definition at line 163 of file IsoTrig.cc.

Referenced by getDistInCM(), GetEtaPhiAtEcal(), and IsoTrig().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ IsoTrig()

◆ ~IsoTrig()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ chgIsolation()

◆ clearChgIsolnTreeVectors()

◆ clearMipCutTreeVectors()

◆ dEta()

◆ dinvPt()

◆ dP()

◆ dPhi()

◆ dPt()

◆ dR()

◆ endJob()

◆ endRun()

◆ etaPhiTrigger()

◆ fillCuts()

◆ fillDescriptions()

◆ fillDifferences()

◆ fillEnergy()

◆ fillHist()

◆ getDistInCM()

◆ GetEtaPhiAtEcal()

◆ getGoodTracks()

◆ pushChgIsolnTreeVecs()

◆ pushMipCutTreeVecs()

◆ studyIsolation()

◆ studyMipCut()

◆ studyTiming()

◆ studyTrigger()

◆ StudyTrkEbyP()

Member Data Documentation

◆ a_charIsoR_

◆ a_coneR_

◆ a_mipR_

◆ a_neutIsoR_

◆ a_neutR1_

◆ a_neutR2_

◆ barrelRecHitsHandle_

◆ beamSpotH_

◆ bField_

◆ bfVal_

◆ changed_

◆ ChgIsolnTree_

◆ cutCharge_

◆ cutMip_

◆ cutNeutral_

◆ doChgIsolTree_

◆ doL2L3_

◆ doMipCutTree_

◆ doStudyIsol_

◆ doTiming_

◆ doTrkResTree_

◆ dr_L1_

◆ endcapRecHitsHandle_

◆ fs_

◆ g_Accepts

◆ g_Pre

◆ g_PreHLT

◆ g_PreL1

◆ geo_

◆ h_dEta

◆ h_dEtaL1

◆ h_dinvPt

◆ h_dP

◆ h_dPhi

◆ h_dPhiL1

◆ h_dPt

◆ h_dRL1

◆ h_eCalo

◆ h_eHcal