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

Inheritance diagram for IsolatedTracksNxN:

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

	~IsolatedTracksNxN () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

void	beginJob () override

void	bookHistograms ()

void	clearTreeVectors ()

void	endJob () override

void	printTrack (const reco::Track *pTrack)

Private Attributes
const bool	debugL1Info_

const int	debugTrks_

const bool	doMC_

double	genPartEtaBins [NEtaBins+1]

double	genPartPBins [NPBins+1]

TH1I *	h_L1AlgoNames

TH1F *	h_maxNearP15x15 [NPBins][NEtaBins]

TH1F *	h_maxNearP21x21 [NPBins][NEtaBins]

TH1F *	h_maxNearP25x25 [NPBins][NEtaBins]

TH1F *	h_maxNearP31x31 [NPBins][NEtaBins]

TH1F *	h_nTracks

TH1F *	h_PVTracksWt

TH1F *	h_recEta_0

TH1F *	h_recEta_1

TH1F *	h_recEta_2

TH2F *	h_recEtaP_0

TH2F *	h_recEtaP_1

TH2F *	h_recEtaP_2

TH2F *	h_recEtaPt_0

TH2F *	h_recEtaPt_1

TH2F *	h_recEtaPt_2

TH1F *	h_recP_0

TH1F *	h_recP_1

TH1F *	h_recP_2

TH1F *	h_recPhi_0

TH1F *	h_recPhi_1

TH1F *	h_recPhi_2

TH1F *	h_recPt_0

TH1F *	h_recPt_1

TH1F *	h_recPt_2

bool	initL1_

std::map< std::pair< unsigned int, std::string >, int >	l1AlgoMap_

const bool	L1TriggerAlgoInfo_

std::unique_ptr< L1GtUtils >	m_l1GtUtils

std::vector< unsigned int >	m_triggerMaskAlgoTrig

const double	maxTrackEta_

const double	minTrackP_

const int	myverbose_

int	nbad_

int	nEventProc_

const bool	printTrkHitPattern_

const double	pvTracksPtMin_

int	t_Bunch

std::vector< double > *	t_e11x11

std::vector< double > *	t_e11x11_10Sig

std::vector< double > *	t_e11x11_15Sig

std::vector< double > *	t_e11x11_20Sig

std::vector< double > *	t_e11x11_25Sig

std::vector< double > *	t_e11x11_30Sig

std::vector< double > *	t_e15x15

std::vector< double > *	t_e15x15_10Sig

std::vector< double > *	t_e15x15_15Sig

std::vector< double > *	t_e15x15_20Sig

std::vector< double > *	t_e15x15_25Sig

std::vector< double > *	t_e15x15_30Sig

std::vector< double > *	t_e7x7

std::vector< double > *	t_e7x7_10Sig

std::vector< double > *	t_e7x7_15Sig

std::vector< double > *	t_e7x7_20Sig

std::vector< double > *	t_e7x7_25Sig

std::vector< double > *	t_e7x7_30Sig

std::vector< double > *	t_e9x9

std::vector< double > *	t_e9x9_10Sig

std::vector< double > *	t_e9x9_15Sig

std::vector< double > *	t_e9x9_20Sig

std::vector< double > *	t_e9x9_25Sig

std::vector< double > *	t_e9x9_30Sig

std::vector< int > *	t_ecalSpike11x11

std::vector< double > *	t_esim11x11

std::vector< double > *	t_esim11x11CharHad

std::vector< double > *	t_esim11x11Matched

std::vector< double > *	t_esim11x11NeutHad

std::vector< double > *	t_esim11x11Photon

std::vector< double > *	t_esim11x11Rest

std::vector< double > *	t_esim15x15

std::vector< double > *	t_esim15x15CharHad

std::vector< double > *	t_esim15x15Matched

std::vector< double > *	t_esim15x15NeutHad

std::vector< double > *	t_esim15x15Photon

std::vector< double > *	t_esim15x15Rest

std::vector< double > *	t_esim7x7

std::vector< double > *	t_esim7x7CharHad

std::vector< double > *	t_esim7x7Matched

std::vector< double > *	t_esim7x7NeutHad

std::vector< double > *	t_esim7x7Photon

std::vector< double > *	t_esim7x7Rest

std::vector< double > *	t_esim9x9

std::vector< double > *	t_esim9x9CharHad

std::vector< double > *	t_esim9x9Matched

std::vector< double > *	t_esim9x9NeutHad

std::vector< double > *	t_esim9x9Photon

std::vector< double > *	t_esim9x9Rest

std::vector< double > *	t_esimPdgId

int	t_EvtNo

std::vector< double > *	t_h3x3

std::vector< double > *	t_h3x3Sig

std::vector< double > *	t_h5x5

std::vector< double > *	t_h5x5Sig

std::vector< double > *	t_h7x7

std::vector< double > *	t_h7x7Sig

std::vector< double > *	t_hsim3x3

std::vector< double > *	t_hsim3x3CharHad

std::vector< double > *	t_hsim3x3Matched

std::vector< double > *	t_hsim3x3NeutHad

std::vector< double > *	t_hsim3x3Photon

std::vector< double > *	t_hsim3x3Rest

std::vector< double > *	t_hsim5x5

std::vector< double > *	t_hsim5x5CharHad

std::vector< double > *	t_hsim5x5Matched

std::vector< double > *	t_hsim5x5NeutHad

std::vector< double > *	t_hsim5x5Photon

std::vector< double > *	t_hsim5x5Rest

std::vector< double > *	t_hsim7x7

std::vector< double > *	t_hsim7x7CharHad

std::vector< double > *	t_hsim7x7Matched

std::vector< double > *	t_hsim7x7NeutHad

std::vector< double > *	t_hsim7x7Photon

std::vector< double > *	t_hsim7x7Rest

std::vector< int > *	t_infoHcal

std::vector< double > *	t_jetEta

std::vector< double > *	t_jetPhi

std::vector< double > *	t_jetPt

std::vector< std::string > *	t_L1AlgoNames

std::vector< double > *	t_L1CenJetEta

std::vector< double > *	t_L1CenJetPhi

std::vector< double > *	t_L1CenJetPt

int	t_L1Decision [128]

std::vector< double > *	t_L1FwdJetEta

std::vector< double > *	t_L1FwdJetPhi

std::vector< double > *	t_L1FwdJetPt

std::vector< double > *	t_L1IsoEMEta

std::vector< double > *	t_L1IsoEMPhi

std::vector< double > *	t_L1IsoEMPt

std::vector< double > *	t_L1METEta

std::vector< double > *	t_L1METPhi

std::vector< double > *	t_L1METPt

std::vector< double > *	t_L1MuonEta

std::vector< double > *	t_L1MuonPhi

std::vector< double > *	t_L1MuonPt

std::vector< double > *	t_L1NonIsoEMEta

std::vector< double > *	t_L1NonIsoEMPhi

std::vector< double > *	t_L1NonIsoEMPt

std::vector< int > *	t_L1PreScale

std::vector< double > *	t_L1TauJetEta

std::vector< double > *	t_L1TauJetPhi

std::vector< double > *	t_L1TauJetPt

int	t_Lumi

std::vector< double > *	t_maxNearHcalP3x3

std::vector< double > *	t_maxNearHcalP5x5

std::vector< double > *	t_maxNearHcalP7x7

std::vector< double > *	t_maxNearP21x21

std::vector< double > *	t_maxNearP31x31

std::vector< int > *	t_NLayersCrossed

int	t_nTracks

std::vector< double > *	t_nTrksJetCalo

std::vector< double > *	t_nTrksJetVtx

std::vector< int > *	t_PVisValid

std::vector< int > *	t_PVndof

std::vector< int > *	t_PVNTracks

std::vector< int > *	t_PVNTracksHP

std::vector< int > *	t_PVNTracksHPWt

std::vector< int > *	t_PVNTracksWt

std::vector< double > *	t_PVTracksSumPt

std::vector< double > *	t_PVTracksSumPtHP

std::vector< double > *	t_PVTracksSumPtHPWt

std::vector< double > *	t_PVTracksSumPtWt

std::vector< double > *	t_PVx

std::vector< double > *	t_PVy

std::vector< double > *	t_PVz

int	t_RunNo

std::vector< double > *	t_simTrackP

std::vector< double > *	t_trackChiSq

std::vector< double > *	t_trackChiSqAll

std::vector< double > *	t_trackDxy

std::vector< double > *	t_trackDxyAll

std::vector< double > *	t_trackDxyBS

std::vector< double > *	t_trackDxyPV

std::vector< double > *	t_trackDxyPVAll

std::vector< double > *	t_trackDz

std::vector< double > *	t_trackDzAll

std::vector< double > *	t_trackDzBS

std::vector< double > *	t_trackDzPV

std::vector< double > *	t_trackDzPVAll

std::vector< double > *	t_trackEcalEta

std::vector< double > *	t_trackEcalPhi

std::vector< double > *	t_trackEta

std::vector< double > *	t_trackEtaAll

std::vector< double > *	t_trackHcalEta

std::vector< double > *	t_trackHcalPhi

std::vector< int > *	t_trackHitInMeasTEC

std::vector< int > *	t_trackHitInMeasTIB

std::vector< int > *	t_trackHitInMeasTID

std::vector< int > *	t_trackHitInMeasTOB

std::vector< int > *	t_trackHitInMissTEC

std::vector< int > *	t_trackHitInMissTIB

std::vector< int > *	t_trackHitInMissTIBTID

std::vector< int > *	t_trackHitInMissTID

std::vector< int > *	t_trackHitInMissTOB

std::vector< int > *	t_trackHitOutMeasTEC

std::vector< int > *	t_trackHitOutMeasTIB

std::vector< int > *	t_trackHitOutMeasTID

std::vector< int > *	t_trackHitOutMeasTOB

std::vector< int > *	t_trackHitOutMissTEC

std::vector< int > *	t_trackHitOutMissTIB

std::vector< int > *	t_trackHitOutMissTID

std::vector< int > *	t_trackHitOutMissTOB

std::vector< int > *	t_trackHitOutMissTOBTEC

std::vector< int > *	t_trackHitsTEC

std::vector< int > *	t_trackHitsTOB

std::vector< double > *	t_trackL

std::vector< int > *	t_trackNOuterHits

std::vector< double > *	t_trackOutPosOutHitDr

std::vector< double > *	t_trackP

std::vector< double > *	t_trackPAll

std::vector< double > *	t_trackPdgIdAll

std::vector< double > *	t_trackPhi

std::vector< double > *	t_trackPhiAll

std::vector< double > *	t_trackPt

std::vector< double > *	t_trackPtAll

std::vector< int > *	t_trackPVIdx

std::vector< double > *	t_trkEcalEne

std::vector< double > *	t_trkHcalEne

const double	tMaxE_

const double	tMaxH_

const double	tMinE_

const double	tMinH_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_caloEB_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_caloEE_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_caloHH_

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	tok_htmap_

edm::EDGetTokenT< reco::CaloJetCollection >	tok_jets_

edm::EDGetTokenT< l1extra::L1JetParticleCollection >	tok_L1extCenJet_

edm::EDGetTokenT< l1extra::L1JetParticleCollection >	tok_L1extFwdJet_

edm::EDGetTokenT< l1extra::L1EmParticleCollection >	tok_L1extIsoEm_

edm::EDGetTokenT< l1extra::L1MuonParticleCollection >	tok_L1extMu_

edm::EDGetTokenT< l1extra::L1EmParticleCollection >	tok_L1extNoIsoEm_

edm::EDGetTokenT< l1extra::L1JetParticleCollection >	tok_L1extTauJet_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_magField_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_

edm::EDGetTokenT< edm::SimTrackContainer >	tok_simTk_

edm::EDGetTokenT< edm::SimVertexContainer >	tok_simVtx_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_topo_

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

TTree *	tree_

const bool	writeAllTracks_

Static Private Attributes
static constexpr size_t	NEtaBins = 3

static constexpr size_t	nL1BitsMax = 128

static constexpr size_t	NPBins = 15

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

template<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 142 of file IsolatedTracksNxN.cc.

Constructor & Destructor Documentation

◆ IsolatedTracksNxN()

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

explicit

Definition at line 307 of file IsolatedTracksNxN.cc.

References edm::EDConsumerBase::consumesCollector(), debugL1Info_, doMC_, L1GtUtils::Event, edm::ParameterSet::getParameter(), HLT_2022v14_cff::InputTag, TFileService::kSharedResource, L1TriggerAlgoInfo_, m_l1GtUtils, maxTrackEta_, minTrackP_, myverbose_, tMaxE_, tMaxH_, tMinE_, tMinH_, tok_bs_, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_htmap_, tok_jets_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extIsoEm_, tok_L1extMu_, tok_L1extNoIsoEm_, tok_L1extTauJet_, tok_magField_, tok_recVtx_, tok_sevlv_, tok_simTk_, tok_simVtx_, tok_topo_, and useL1GtTriggerMenuLite().

     : trackerHitAssociatorConfig_(consumesCollector()),
       doMC_(iConfig.getUntrackedParameter<bool>("doMC", false)),
       writeAllTracks_(iConfig.getUntrackedParameter<bool>("writeAllTracks", false)),
       myverbose_(iConfig.getUntrackedParameter<int>("verbosity", 5)),
       pvTracksPtMin_(iConfig.getUntrackedParameter<double>("pvTracksPtMin", 1.0)),
       debugTrks_(iConfig.getUntrackedParameter<int>("debugTracks", 0)),
       printTrkHitPattern_(iConfig.getUntrackedParameter<bool>("printTrkHitPattern", false)),
       minTrackP_(iConfig.getUntrackedParameter<double>("minTrackP", 1.0)),
       maxTrackEta_(iConfig.getUntrackedParameter<double>("maxTrackEta", 5.0)),
       debugL1Info_(iConfig.getUntrackedParameter<bool>("debugL1Info", false)),
       L1TriggerAlgoInfo_(iConfig.getUntrackedParameter<bool>("l1TriggerAlgoInfo", false)),
       tMinE_(iConfig.getUntrackedParameter<double>("timeMinCutECAL", -500.)),
       tMaxE_(iConfig.getUntrackedParameter<double>("timeMaxCutECAL", 500.)),
       tMinH_(iConfig.getUntrackedParameter<double>("timeMinCutHCAL", -500.)),
       tMaxH_(iConfig.getUntrackedParameter<double>("timeMaxCutHCAL", 500.)),
       t_L1AlgoNames(nullptr),
       t_L1PreScale(nullptr),
       t_L1CenJetPt(nullptr),
       t_L1CenJetEta(nullptr),
       t_L1CenJetPhi(nullptr),
       t_L1FwdJetPt(nullptr),
       t_L1FwdJetEta(nullptr),
       t_L1FwdJetPhi(nullptr),
       t_L1TauJetPt(nullptr),
       t_L1TauJetEta(nullptr),
       t_L1TauJetPhi(nullptr),
       t_L1MuonPt(nullptr),
       t_L1MuonEta(nullptr),
       t_L1MuonPhi(nullptr),
       t_L1IsoEMPt(nullptr),
       t_L1IsoEMEta(nullptr),
       t_L1IsoEMPhi(nullptr),
       t_L1NonIsoEMPt(nullptr),
       t_L1NonIsoEMEta(nullptr),
       t_L1NonIsoEMPhi(nullptr),
       t_L1METPt(nullptr),
       t_L1METEta(nullptr),
       t_L1METPhi(nullptr),
       t_PVx(nullptr),
       t_PVy(nullptr),
       t_PVz(nullptr),
       t_PVTracksSumPt(nullptr),
       t_PVTracksSumPtWt(nullptr),
       t_PVTracksSumPtHP(nullptr),
       t_PVTracksSumPtHPWt(nullptr),
       t_PVisValid(nullptr),
       t_PVNTracks(nullptr),
       t_PVNTracksWt(nullptr),
       t_PVndof(nullptr),
       t_PVNTracksHP(nullptr),
       t_PVNTracksHPWt(nullptr),
       t_jetPt(nullptr),
       t_jetEta(nullptr),
       t_jetPhi(nullptr),
       t_nTrksJetCalo(nullptr),
       t_nTrksJetVtx(nullptr),
       t_trackPAll(nullptr),
       t_trackEtaAll(nullptr),
       t_trackPhiAll(nullptr),
       t_trackPdgIdAll(nullptr),
       t_trackPtAll(nullptr),
       t_trackDxyAll(nullptr),
       t_trackDzAll(nullptr),
       t_trackDxyPVAll(nullptr),
       t_trackDzPVAll(nullptr),
       t_trackChiSqAll(nullptr),
       t_trackP(nullptr),
       t_trackPt(nullptr),
       t_trackEta(nullptr),
       t_trackPhi(nullptr),
       t_trackEcalEta(nullptr),
       t_trackEcalPhi(nullptr),
       t_trackHcalEta(nullptr),
       t_trackHcalPhi(nullptr),
       t_trackDxy(nullptr),
       t_trackDxyBS(nullptr),
       t_trackDz(nullptr),
       t_trackDzBS(nullptr),
       t_trackDxyPV(nullptr),
       t_trackDzPV(nullptr),
       t_trackChiSq(nullptr),
       t_trackPVIdx(nullptr),
       t_NLayersCrossed(nullptr),
       t_trackNOuterHits(nullptr),
       t_trackHitsTOB(nullptr),
       t_trackHitsTEC(nullptr),
       t_trackHitInMissTOB(nullptr),
       t_trackHitInMissTEC(nullptr),
       t_trackHitInMissTIB(nullptr),
       t_trackHitInMissTID(nullptr),
       t_trackHitInMissTIBTID(nullptr),
       t_trackHitOutMissTOB(nullptr),
       t_trackHitOutMissTEC(nullptr),
       t_trackHitOutMissTIB(nullptr),
       t_trackHitOutMissTID(nullptr),
       t_trackHitOutMissTOBTEC(nullptr),
       t_trackHitInMeasTOB(nullptr),
       t_trackHitInMeasTEC(nullptr),
       t_trackHitInMeasTIB(nullptr),
       t_trackHitInMeasTID(nullptr),
       t_trackHitOutMeasTOB(nullptr),
       t_trackHitOutMeasTEC(nullptr),
       t_trackHitOutMeasTIB(nullptr),
       t_trackHitOutMeasTID(nullptr),
       t_trackOutPosOutHitDr(nullptr),
       t_trackL(nullptr),
       t_maxNearP31x31(nullptr),
       t_maxNearP21x21(nullptr),
       t_ecalSpike11x11(nullptr),
       t_e7x7(nullptr),
       t_e9x9(nullptr),
       t_e11x11(nullptr),
       t_e15x15(nullptr),
       t_e7x7_10Sig(nullptr),
       t_e9x9_10Sig(nullptr),
       t_e11x11_10Sig(nullptr),
       t_e15x15_10Sig(nullptr),
       t_e7x7_15Sig(nullptr),
       t_e9x9_15Sig(nullptr),
       t_e11x11_15Sig(nullptr),
       t_e15x15_15Sig(nullptr),
       t_e7x7_20Sig(nullptr),
       t_e9x9_20Sig(nullptr),
       t_e11x11_20Sig(nullptr),
       t_e15x15_20Sig(nullptr),
       t_e7x7_25Sig(nullptr),
       t_e9x9_25Sig(nullptr),
       t_e11x11_25Sig(nullptr),
       t_e15x15_25Sig(nullptr),
       t_e7x7_30Sig(nullptr),
       t_e9x9_30Sig(nullptr),
       t_e11x11_30Sig(nullptr),
       t_e15x15_30Sig(nullptr),
       t_esimPdgId(nullptr),
       t_simTrackP(nullptr),
       t_trkEcalEne(nullptr),
       t_esim7x7(nullptr),
       t_esim9x9(nullptr),
       t_esim11x11(nullptr),
       t_esim15x15(nullptr),
       t_esim7x7Matched(nullptr),
       t_esim9x9Matched(nullptr),
       t_esim11x11Matched(nullptr),
       t_esim15x15Matched(nullptr),
       t_esim7x7Rest(nullptr),
       t_esim9x9Rest(nullptr),
       t_esim11x11Rest(nullptr),
       t_esim15x15Rest(nullptr),
       t_esim7x7Photon(nullptr),
       t_esim9x9Photon(nullptr),
       t_esim11x11Photon(nullptr),
       t_esim15x15Photon(nullptr),
       t_esim7x7NeutHad(nullptr),
       t_esim9x9NeutHad(nullptr),
       t_esim11x11NeutHad(nullptr),
       t_esim15x15NeutHad(nullptr),
       t_esim7x7CharHad(nullptr),
       t_esim9x9CharHad(nullptr),
       t_esim11x11CharHad(nullptr),
       t_esim15x15CharHad(nullptr),
       t_maxNearHcalP3x3(nullptr),
       t_maxNearHcalP5x5(nullptr),
       t_maxNearHcalP7x7(nullptr),
       t_h3x3(nullptr),
       t_h5x5(nullptr),
       t_h7x7(nullptr),
       t_h3x3Sig(nullptr),
       t_h5x5Sig(nullptr),
       t_h7x7Sig(nullptr),
       t_infoHcal(nullptr),
       t_trkHcalEne(nullptr),
       t_hsim3x3(nullptr),
       t_hsim5x5(nullptr),
       t_hsim7x7(nullptr),
       t_hsim3x3Matched(nullptr),
       t_hsim5x5Matched(nullptr),
       t_hsim7x7Matched(nullptr),
       t_hsim3x3Rest(nullptr),
       t_hsim5x5Rest(nullptr),
       t_hsim7x7Rest(nullptr),
       t_hsim3x3Photon(nullptr),
       t_hsim5x5Photon(nullptr),
       t_hsim7x7Photon(nullptr),
       t_hsim3x3NeutHad(nullptr),
       t_hsim5x5NeutHad(nullptr),
       t_hsim7x7NeutHad(nullptr),
       t_hsim3x3CharHad(nullptr),
       t_hsim5x5CharHad(nullptr),
       t_hsim7x7CharHad(nullptr) {
   if (L1TriggerAlgoInfo_) {
     m_l1GtUtils = std::make_unique<L1GtUtils>(
         iConfig, consumesCollector(), useL1GtTriggerMenuLite, *this, L1GtUtils::UseEventSetupIn::Event);
   }
 
   usesResource(TFileService::kSharedResource);
 
   //now do what ever initialization is needed
 
   edm::InputTag L1extraTauJetSource_ = iConfig.getParameter<edm::InputTag>("l1extraTauJetSource");
   edm::InputTag L1extraCenJetSource_ = iConfig.getParameter<edm::InputTag>("l1extraCenJetSource");
   edm::InputTag L1extraFwdJetSource_ = iConfig.getParameter<edm::InputTag>("l1extraFwdJetSource");
   edm::InputTag L1extraMuonSource_ = iConfig.getParameter<edm::InputTag>("l1extraMuonSource");
   edm::InputTag L1extraIsoEmSource_ = iConfig.getParameter<edm::InputTag>("l1extraIsoEmSource");
   edm::InputTag L1extraNonIsoEmSource_ = iConfig.getParameter<edm::InputTag>("l1extraNonIsoEmSource");
   edm::InputTag L1GTReadoutRcdSource_ = iConfig.getParameter<edm::InputTag>("l1GTReadoutRcdSource");
   edm::InputTag L1GTObjectMapRcdSource_ = iConfig.getParameter<edm::InputTag>("l1GTObjectMapRcdSource");
   edm::InputTag JetSrc_ = iConfig.getParameter<edm::InputTag>("jetSource");
   edm::InputTag JetExtender_ = iConfig.getParameter<edm::InputTag>("jetExtender");
   edm::InputTag HBHERecHitSource_ = iConfig.getParameter<edm::InputTag>("hbheRecHitSource");
 
   // define tokens for access
   tok_L1extTauJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraTauJetSource_);
   tok_L1extCenJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraCenJetSource_);
   tok_L1extFwdJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraFwdJetSource_);
   tok_L1extMu_ = consumes<l1extra::L1MuonParticleCollection>(L1extraMuonSource_);
   tok_L1extIsoEm_ = consumes<l1extra::L1EmParticleCollection>(L1extraIsoEmSource_);
   tok_L1extNoIsoEm_ = consumes<l1extra::L1EmParticleCollection>(L1extraNonIsoEmSource_);
   tok_jets_ = consumes<reco::CaloJetCollection>(JetSrc_);
   tok_hbhe_ = consumes<HBHERecHitCollection>(HBHERecHitSource_);
 
   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_simTk_ = consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
   tok_simVtx_ = consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
   tok_caloEB_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEB"));
   tok_caloEE_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEE"));
   tok_caloHH_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalHits"));
 
   if (myverbose_ >= 0) {
     edm::LogVerbatim("IsoTrack") << "Parameters read from config file \n"
                                  << " doMC        " << doMC_ << "\t myverbose " << myverbose_ << "\t minTrackP "
                                  << minTrackP_ << "\t maxTrackEta " << maxTrackEta_ << "\t tMinE " << tMinE_
                                  << "\t tMaxE " << tMaxE_ << "\t tMinH " << tMinH_ << "\t tMaxH " << tMaxH_
                                  << "\n debugL1Info " << debugL1Info_ << "\t L1TriggerAlgoInfo " << L1TriggerAlgoInfo_
                                  << "\n";
   }
 
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
   tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
   tok_topo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
   tok_magField_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
   tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
   tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
   tok_htmap_ = esConsumes<EcalTrigTowerConstituentsMap, IdealGeometryRecord>();
 }

◆ ~IsolatedTracksNxN()

IsolatedTracksNxN::~IsolatedTracksNxN ( )

override

Definition at line 558 of file IsolatedTracksNxN.cc.

References t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1AlgoNames, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1METEta, t_L1METPhi, t_L1METPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1PreScale, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP21x21, t_maxNearP31x31, t_NLayersCrossed, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, and t_trkHcalEne.

                                       {
   delete t_PVx;
   delete t_PVy;
   delete t_PVz;
   delete t_PVisValid;
   delete t_PVndof;
   delete t_PVNTracks;
   delete t_PVNTracksWt;
   delete t_PVTracksSumPt;
   delete t_PVTracksSumPtWt;
   delete t_PVNTracksHP;
   delete t_PVNTracksHPWt;
   delete t_PVTracksSumPtHP;
   delete t_PVTracksSumPtHPWt;
   delete t_L1AlgoNames;
   delete t_L1PreScale;
   delete t_L1CenJetPt;
   delete t_L1CenJetEta;
   delete t_L1CenJetPhi;
   delete t_L1FwdJetPt;
   delete t_L1FwdJetEta;
   delete t_L1FwdJetPhi;
   delete t_L1TauJetPt;
   delete t_L1TauJetEta;
   delete t_L1TauJetPhi;
   delete t_L1MuonPt;
   delete t_L1MuonEta;
   delete t_L1MuonPhi;
   delete t_L1IsoEMPt;
   delete t_L1IsoEMEta;
   delete t_L1IsoEMPhi;
   delete t_L1NonIsoEMPt;
   delete t_L1NonIsoEMEta;
   delete t_L1NonIsoEMPhi;
   delete t_L1METPt;
   delete t_L1METEta;
   delete t_L1METPhi;
   delete t_jetPt;
   delete t_jetEta;
   delete t_jetPhi;
   delete t_nTrksJetCalo;
   delete t_nTrksJetVtx;
   delete t_trackPAll;
   delete t_trackEtaAll;
   delete t_trackPhiAll;
   delete t_trackPdgIdAll;
   delete t_trackPtAll;
   delete t_trackDxyAll;
   delete t_trackDzAll;
   delete t_trackDxyPVAll;
   delete t_trackDzPVAll;
   delete t_trackChiSqAll;
   delete t_trackP;
   delete t_trackPt;
   delete t_trackEta;
   delete t_trackPhi;
   delete t_trackEcalEta;
   delete t_trackEcalPhi;
   delete t_trackHcalEta;
   delete t_trackHcalPhi;
   delete t_trackNOuterHits;
   delete t_NLayersCrossed;
   delete t_trackDxy;
   delete t_trackDxyBS;
   delete t_trackDz;
   delete t_trackDzBS;
   delete t_trackDxyPV;
   delete t_trackDzPV;
   delete t_trackPVIdx;
   delete t_trackChiSq;
   delete t_trackHitsTOB;
   delete t_trackHitsTEC;
   delete t_trackHitInMissTOB;
   delete t_trackHitInMissTEC;
   delete t_trackHitInMissTIB;
   delete t_trackHitInMissTID;
   delete t_trackHitInMissTIBTID;
   delete t_trackHitOutMissTOB;
   delete t_trackHitOutMissTEC;
   delete t_trackHitOutMissTIB;
   delete t_trackHitOutMissTID;
   delete t_trackHitOutMissTOBTEC;
   delete t_trackHitInMeasTOB;
   delete t_trackHitInMeasTEC;
   delete t_trackHitInMeasTIB;
   delete t_trackHitInMeasTID;
   delete t_trackHitOutMeasTOB;
   delete t_trackHitOutMeasTEC;
   delete t_trackHitOutMeasTIB;
   delete t_trackHitOutMeasTID;
   delete t_trackOutPosOutHitDr;
   delete t_trackL;
   delete t_maxNearP31x31;
   delete t_maxNearP21x21;
   delete t_ecalSpike11x11;
   delete t_e7x7;
   delete t_e9x9;
   delete t_e11x11;
   delete t_e15x15;
   delete t_e7x7_10Sig;
   delete t_e9x9_10Sig;
   delete t_e11x11_10Sig;
   delete t_e15x15_10Sig;
   delete t_e7x7_15Sig;
   delete t_e9x9_15Sig;
   delete t_e11x11_15Sig;
   delete t_e15x15_15Sig;
   delete t_e7x7_20Sig;
   delete t_e9x9_20Sig;
   delete t_e11x11_20Sig;
   delete t_e15x15_20Sig;
   delete t_e7x7_25Sig;
   delete t_e9x9_25Sig;
   delete t_e11x11_25Sig;
   delete t_e15x15_25Sig;
   delete t_e7x7_30Sig;
   delete t_e9x9_30Sig;
   delete t_e11x11_30Sig;
   delete t_e15x15_30Sig;
   delete t_esim7x7;
   delete t_esim9x9;
   delete t_esim11x11;
   delete t_esim15x15;
   delete t_esim7x7Matched;
   delete t_esim9x9Matched;
   delete t_esim11x11Matched;
   delete t_esim15x15Matched;
   delete t_esim7x7Rest;
   delete t_esim9x9Rest;
   delete t_esim11x11Rest;
   delete t_esim15x15Rest;
   delete t_esim7x7Photon;
   delete t_esim9x9Photon;
   delete t_esim11x11Photon;
   delete t_esim15x15Photon;
   delete t_esim7x7NeutHad;
   delete t_esim9x9NeutHad;
   delete t_esim11x11NeutHad;
   delete t_esim15x15NeutHad;
   delete t_esim7x7CharHad;
   delete t_esim9x9CharHad;
   delete t_esim11x11CharHad;
   delete t_esim15x15CharHad;
   delete t_trkEcalEne;
   delete t_simTrackP;
   delete t_esimPdgId;
   delete t_maxNearHcalP3x3;
   delete t_maxNearHcalP5x5;
   delete t_maxNearHcalP7x7;
   delete t_h3x3;
   delete t_h5x5;
   delete t_h7x7;
   delete t_h3x3Sig;
   delete t_h5x5Sig;
   delete t_h7x7Sig;
   delete t_infoHcal;
   delete t_trkHcalEne;
   delete t_hsim3x3;
   delete t_hsim5x5;
   delete t_hsim7x7;
   delete t_hsim3x3Matched;
   delete t_hsim5x5Matched;
   delete t_hsim7x7Matched;
   delete t_hsim3x3Rest;
   delete t_hsim5x5Rest;
   delete t_hsim7x7Rest;
   delete t_hsim3x3Photon;
   delete t_hsim5x5Photon;
   delete t_hsim7x7Photon;
   delete t_hsim3x3NeutHad;
   delete t_hsim5x5NeutHad;
   delete t_hsim7x7NeutHad;
   delete t_hsim3x3CharHad;
   delete t_hsim5x5CharHad;
   delete t_hsim7x7CharHad;
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 766 of file IsolatedTracksNxN.cc.

References funct::abs(), Calorimetry_cff::bField, spr::chargeIsolationEcal(), spr::chargeIsolationHcal(), clearTreeVectors(), debugL1Info_, debugTrks_, TauDecayModes::dec, reco::deltaR(), change_name::diff, doMC_, reco::TrackBase::dxy(), reco::TrackBase::dz(), EcalBarrel, spr::eCaloSimInfo(), spr::caloSimInfo::eChargedHad, spr::eECALmatrix(), spr::eECALSimInfo(), spr::caloSimInfo::eGamma, spr::eHCALmatrix(), spr::eHCALmatrixCell(), spr::eHCALSimInfo(), spr::caloSimInfo::eMatched, spr::caloSimInfo::eNeutralHad, spr::caloSimInfo::eRest, PVValHelper::eta, reco::TrackBase::eta(), HLT_2022v14_cff::eta1, spr::caloSimInfo::eTotal, first, genPartEtaBins, genPartPBins, edm::RefToBase< T >::get(), edm::EventSetup::getData(), L1GtTriggerMenu::gtAlgorithmMap(), L1GtTriggerMenu::gtTriggerMenuName(), h_L1AlgoNames, h_maxNearP15x15, h_maxNearP21x21, h_maxNearP25x25, h_maxNearP31x31, h_nTracks, h_PVTracksWt, h_recEta_0, h_recEta_1, h_recEta_2, h_recEtaP_0, h_recEtaP_1, h_recEtaP_2, h_recEtaPt_0, h_recEtaPt_1, h_recEtaPt_2, h_recP_0, h_recP_1, h_recP_2, h_recPhi_0, h_recPhi_1, h_recPhi_2, h_recPt_0, h_recPt_1, h_recPt_2, photonIsolationHIProducer_cfi::hbhe, reco::TrackBase::hitPattern(), mps_fire::i, LEDCalibrationChannels::ieta, iEvent, cuy::ii, initL1_, sistrip::SpyUtilities::isValid(), edm::HandleBase::isValid(), gpuVertexFinder::iv, PDWG_EXODelayedJetMET_cff::jets, l1AlgoMap_, L1TriggerAlgoInfo_, m_l1GtUtils, spr::matchedSimTrack(), maxTrackEta_, minTrackP_, reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::MISSING_OUTER_HITS, reco::TrackBase::momentum(), myverbose_, nbad_, ndof, NEtaBins, nEventProc_, reco::TrackBase::normalizedChi2(), NPBins, BeamSpotPI::nTracks, reco::Track::outerPosition(), reco::TrackBase::p(), LaserDQM_cfg::p1, spr::caloSimInfo::pdgMatched, phi, reco::BeamSpot::position(), spr::propagateCALO(), spr::propagateTrackerEnd(), DiDispStaMuonMonitor_cfi::pt, reco::TrackBase::pt(), HLT_2022v14_cff::pt1, pvTracksPtMin_, reco::TrackBase::quality(), reco::TrackBase::qualityByName(), AlCaHLTBitMon_QueryRunRegistry::string, reco::HitPattern::stripTECLayersWithMeasurement(), reco::HitPattern::stripTECLayersWithoutMeasurement(), reco::HitPattern::stripTIBLayersWithMeasurement(), reco::HitPattern::stripTIBLayersWithoutMeasurement(), reco::HitPattern::stripTIDLayersWithMeasurement(), reco::HitPattern::stripTIDLayersWithoutMeasurement(), reco::HitPattern::stripTOBLayersWithMeasurement(), reco::HitPattern::stripTOBLayersWithoutMeasurement(), t_Bunch, t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_EvtNo, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1AlgoNames, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1Decision, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1PreScale, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_Lumi, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP21x21, t_maxNearP31x31, t_NLayersCrossed, t_nTracks, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_RunNo, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, t_trkHcalEne, tMaxE_, tMaxH_, tMinE_, tMinH_, tok_bs_, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_htmap_, tok_jets_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extIsoEm_, tok_L1extMu_, tok_L1extNoIsoEm_, tok_L1extTauJet_, tok_magField_, tok_recVtx_, tok_sevlv_, tok_simTk_, tok_simVtx_, tok_topo_, trackerHitAssociatorConfig_, reco::HitPattern::trackerLayersWithMeasurement(), tree_, useL1EventSetup(), useL1GtTriggerMenuLite(), writeAllTracks_, x, y, and z.

                                                                                    {
   bool haveIsoTrack = false;
 
   const MagneticField *bField = &iSetup.getData(tok_magField_);
 
   clearTreeVectors();
 
   t_RunNo = iEvent.id().run();
   t_EvtNo = iEvent.id().event();
   t_Lumi = iEvent.luminosityBlock();
   t_Bunch = iEvent.bunchCrossing();
 
   ++nEventProc_;
 
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   if (debugTrks_ > 1) {
     edm::LogVerbatim("IsoTrack") << "Track Collection: ";
     edm::LogVerbatim("IsoTrack") << "Number of Tracks " << trkCollection->size();
   }
   std::string theTrackQuality = "highPurity";
   reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
 
   //===================== save L1 Trigger information =======================
   if (L1TriggerAlgoInfo_) {
     m_l1GtUtils->getL1GtRunCache(iEvent, iSetup, useL1EventSetup, useL1GtTriggerMenuLite);
 
     int iErrorCode = -1;
     int l1ConfCode = -1;
     const bool l1Conf = m_l1GtUtils->availableL1Configuration(iErrorCode, l1ConfCode);
     if (!l1Conf) {
       edm::LogVerbatim("IsoTrack")
           << "\nL1 configuration code:" << l1ConfCode << "\nNo valid L1 trigger configuration available."
           << "\nSee text above for error code interpretation"
           << "\nNo return here, in order to test each method, protected against configuration error.";
     }
 
     const L1GtTriggerMenu *m_l1GtMenu = m_l1GtUtils->ptrL1TriggerMenuEventSetup(iErrorCode);
     const AlgorithmMap &algorithmMap = m_l1GtMenu->gtAlgorithmMap();
     const std::string &menuName = m_l1GtMenu->gtTriggerMenuName();
 
     if (!initL1_) {
       initL1_ = true;
       edm::LogVerbatim("IsoTrack") << "menuName " << menuName;
       for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); itAlgo++) {
         std::string algName = itAlgo->first;
         int algBitNumber = (itAlgo->second).algoBitNumber();
         l1AlgoMap_.insert(std::pair<std::pair<unsigned int, std::string>, int>(
             std::pair<unsigned int, std::string>(algBitNumber, algName), 0));
       }
       std::map<std::pair<unsigned int, std::string>, int>::iterator itr;
       for (itr = l1AlgoMap_.begin(); itr != l1AlgoMap_.end(); itr++) {
         edm::LogVerbatim("IsoTrack") << " ********** " << (itr->first).first << " " << (itr->first).second << " "
                                      << itr->second;
       }
     }
 
     std::vector<int> algbits;
     for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); itAlgo++) {
       std::string algName = itAlgo->first;
       int algBitNumber = (itAlgo->second).algoBitNumber();
       bool decision = m_l1GtUtils->decision(iEvent, itAlgo->first, iErrorCode);
       int preScale = m_l1GtUtils->prescaleFactor(iEvent, itAlgo->first, iErrorCode);
 
       // save the algo names which fired
       if (decision) {
         l1AlgoMap_[std::pair<unsigned int, std::string>(algBitNumber, algName)] += 1;
         h_L1AlgoNames->Fill(algBitNumber);
         t_L1AlgoNames->push_back(itAlgo->first);
         t_L1PreScale->push_back(preScale);
         t_L1Decision[algBitNumber] = 1;
         algbits.push_back(algBitNumber);
       }
     }
 
     if (debugL1Info_) {
       for (unsigned int ii = 0; ii < t_L1AlgoNames->size(); ii++) {
         edm::LogVerbatim("IsoTrack") << ii << " " << (*t_L1AlgoNames)[ii] << " " << (*t_L1PreScale)[ii] << " "
                                      << algbits[ii];
       }
       for (int i = 0; i < 128; ++i) {
         edm::LogVerbatim("IsoTrack") << "L1Decision: " << i << ":" << t_L1Decision[i];
       }
     }
 
     // L1Taus
     edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
     iEvent.getByToken(tok_L1extTauJet_, l1TauHandle);
     l1extra::L1JetParticleCollection::const_iterator itr;
     int iL1Obj = 0;
     for (itr = l1TauHandle->begin(), iL1Obj = 0; itr != l1TauHandle->end(); ++itr, iL1Obj++) {
       if (iL1Obj < 1) {
         t_L1TauJetPt->push_back(itr->pt());
         t_L1TauJetEta->push_back(itr->eta());
         t_L1TauJetPhi->push_back(itr->phi());
       }
       if (debugL1Info_) {
         edm::LogVerbatim("IsoTrack") << "tauJ p/pt  " << itr->momentum() << " " << itr->pt() << "  eta/phi "
                                      << itr->eta() << " " << itr->phi();
       }
     }
 
     // L1 Central Jets
     edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
     iEvent.getByToken(tok_L1extCenJet_, l1CenJetHandle);
     for (itr = l1CenJetHandle->begin(), iL1Obj = 0; itr != l1CenJetHandle->end(); ++itr, iL1Obj++) {
       if (iL1Obj < 1) {
         t_L1CenJetPt->push_back(itr->pt());
         t_L1CenJetEta->push_back(itr->eta());
         t_L1CenJetPhi->push_back(itr->phi());
       }
       if (debugL1Info_) {
         edm::LogVerbatim("IsoTrack") << "cenJ p/pt     " << itr->momentum() << " " << itr->pt() << "  eta/phi "
                                      << itr->eta() << " " << itr->phi();
       }
     }
 
     // L1 Forward Jets
     edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
     iEvent.getByToken(tok_L1extFwdJet_, l1FwdJetHandle);
     for (itr = l1FwdJetHandle->begin(), iL1Obj = 0; itr != l1FwdJetHandle->end(); ++itr, iL1Obj++) {
       if (iL1Obj < 1) {
         t_L1FwdJetPt->push_back(itr->pt());
         t_L1FwdJetEta->push_back(itr->eta());
         t_L1FwdJetPhi->push_back(itr->phi());
       }
       if (debugL1Info_) {
         edm::LogVerbatim("IsoTrack") << "fwdJ p/pt     " << itr->momentum() << " " << itr->pt() << "  eta/phi "
                                      << itr->eta() << " " << itr->phi();
       }
     }
 
     // L1 Isolated EM onjects
     l1extra::L1EmParticleCollection::const_iterator itrEm;
     edm::Handle<l1extra::L1EmParticleCollection> l1IsoEmHandle;
     iEvent.getByToken(tok_L1extIsoEm_, l1IsoEmHandle);
     for (itrEm = l1IsoEmHandle->begin(), iL1Obj = 0; itrEm != l1IsoEmHandle->end(); ++itrEm, iL1Obj++) {
       if (iL1Obj < 1) {
         t_L1IsoEMPt->push_back(itrEm->pt());
         t_L1IsoEMEta->push_back(itrEm->eta());
         t_L1IsoEMPhi->push_back(itrEm->phi());
       }
       if (debugL1Info_) {
         edm::LogVerbatim("IsoTrack") << "isoEm p/pt    " << itrEm->momentum() << " " << itrEm->pt() << "  eta/phi "
                                      << itrEm->eta() << " " << itrEm->phi();
       }
     }
 
     // L1 Non-Isolated EM onjects
     edm::Handle<l1extra::L1EmParticleCollection> l1NonIsoEmHandle;
     iEvent.getByToken(tok_L1extNoIsoEm_, l1NonIsoEmHandle);
     for (itrEm = l1NonIsoEmHandle->begin(), iL1Obj = 0; itrEm != l1NonIsoEmHandle->end(); ++itrEm, iL1Obj++) {
       if (iL1Obj < 1) {
         t_L1NonIsoEMPt->push_back(itrEm->pt());
         t_L1NonIsoEMEta->push_back(itrEm->eta());
         t_L1NonIsoEMPhi->push_back(itrEm->phi());
       }
       if (debugL1Info_) {
         edm::LogVerbatim("IsoTrack") << "nonIsoEm p/pt " << itrEm->momentum() << " " << itrEm->pt() << "  eta/phi "
                                      << itrEm->eta() << " " << itrEm->phi();
       }
     }
 
     // L1 Muons
     l1extra::L1MuonParticleCollection::const_iterator itrMu;
     edm::Handle<l1extra::L1MuonParticleCollection> l1MuHandle;
     iEvent.getByToken(tok_L1extMu_, l1MuHandle);
     for (itrMu = l1MuHandle->begin(), iL1Obj = 0; itrMu != l1MuHandle->end(); ++itrMu, iL1Obj++) {
       if (iL1Obj < 1) {
         t_L1MuonPt->push_back(itrMu->pt());
         t_L1MuonEta->push_back(itrMu->eta());
         t_L1MuonPhi->push_back(itrMu->phi());
       }
       if (debugL1Info_) {
         edm::LogVerbatim("IsoTrack") << "l1muon p/pt   " << itrMu->momentum() << " " << itrMu->pt() << "  eta/phi "
                                      << itrMu->eta() << " " << itrMu->phi();
       }
     }
   }
 
   //============== store the information about all the Non-Fake vertices ===============
 
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);
 
   std::vector<reco::Track> svTracks;
   math::XYZPoint leadPV(0, 0, 0);
   double sumPtMax = -1.0;
   for (unsigned int ind = 0; ind < recVtxs->size(); ind++) {
     if (!((*recVtxs)[ind].isFake())) {
       double vtxTrkSumPt = 0.0, vtxTrkSumPtWt = 0.0;
       int vtxTrkNWt = 0;
       double vtxTrkSumPtHP = 0.0, vtxTrkSumPtHPWt = 0.0;
       int vtxTrkNHP = 0, vtxTrkNHPWt = 0;
 
       reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
 
       for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack != (*recVtxs)[ind].tracks_end(); vtxTrack++) {
         if ((*vtxTrack)->pt() < pvTracksPtMin_)
           continue;
 
         bool trkQuality = (*vtxTrack)->quality(trackQuality_);
 
         vtxTrkSumPt += (*vtxTrack)->pt();
 
         vtxTrkSumPt += (*vtxTrack)->pt();
         if (trkQuality) {
           vtxTrkSumPtHP += (*vtxTrack)->pt();
           vtxTrkNHP++;
         }
 
         double weight = (*recVtxs)[ind].trackWeight(*vtxTrack);
         h_PVTracksWt->Fill(weight);
         if (weight > 0.5) {
           vtxTrkSumPtWt += (*vtxTrack)->pt();
           vtxTrkNWt++;
           if (trkQuality) {
             vtxTrkSumPtHPWt += (*vtxTrack)->pt();
             vtxTrkNHPWt++;
           }
         }
       }
 
       if (vtxTrkSumPt > sumPtMax) {
         sumPtMax = vtxTrkSumPt;
         leadPV = math::XYZPoint((*recVtxs)[ind].x(), (*recVtxs)[ind].y(), (*recVtxs)[ind].z());
       }
 
       t_PVx->push_back((*recVtxs)[ind].x());
       t_PVy->push_back((*recVtxs)[ind].y());
       t_PVz->push_back((*recVtxs)[ind].z());
       t_PVisValid->push_back((*recVtxs)[ind].isValid());
       t_PVNTracks->push_back((*recVtxs)[ind].tracksSize());
       t_PVndof->push_back((*recVtxs)[ind].ndof());
       t_PVNTracksWt->push_back(vtxTrkNWt);
       t_PVTracksSumPt->push_back(vtxTrkSumPt);
       t_PVTracksSumPtWt->push_back(vtxTrkSumPtWt);
 
       t_PVNTracksHP->push_back(vtxTrkNHP);
       t_PVNTracksHPWt->push_back(vtxTrkNHPWt);
       t_PVTracksSumPtHP->push_back(vtxTrkSumPtHP);
       t_PVTracksSumPtHPWt->push_back(vtxTrkSumPtHPWt);
 
       if (myverbose_ == 4) {
         edm::LogVerbatim("IsoTrack") << "PV " << ind << " isValid " << (*recVtxs)[ind].isValid() << " isFake "
                                      << (*recVtxs)[ind].isFake() << " hasRefittedTracks() " << ind << " "
                                      << (*recVtxs)[ind].hasRefittedTracks() << " refittedTrksSize "
                                      << (*recVtxs)[ind].refittedTracks().size() << "  tracksSize() "
                                      << (*recVtxs)[ind].tracksSize() << " sumPt " << vtxTrkSumPt;
       }
     }  // if vtx is not Fake
   }    // loop over PVs
   //===================================================================================
 
   // Get the beamspot
   edm::Handle<reco::BeamSpot> beamSpotH;
   iEvent.getByToken(tok_bs_, beamSpotH);
   math::XYZPoint bspot;
   bspot = (beamSpotH.isValid()) ? beamSpotH->position() : math::XYZPoint(0, 0, 0);
 
   //=====================================================================
 
   edm::Handle<reco::CaloJetCollection> jets;
   iEvent.getByToken(tok_jets_, jets);
   //  edm::Handle<reco::JetExtendedAssociation::Container> jetExtender;
   //  iEvent.getByLabel(JetExtender_,jetExtender);
 
   for (unsigned int ijet = 0; ijet < (*jets).size(); ijet++) {
     t_jetPt->push_back((*jets)[ijet].pt());
     t_jetEta->push_back((*jets)[ijet].eta());
     t_jetPhi->push_back((*jets)[ijet].phi());
     t_nTrksJetVtx->push_back(-1.0);
     t_nTrksJetCalo->push_back(-1.0);
   }
 
   //===================================================================================
 
   // get handles to calogeometry and calotopology
   const CaloGeometry *geo = &iSetup.getData(tok_geom_);
   const CaloTopology *caloTopology = &iSetup.getData(tok_caloTopology_);
   const HcalTopology *theHBHETopology = &iSetup.getData(tok_topo_);
 
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
 
   // Retrieve the good/bad ECAL channels from the DB
   const EcalChannelStatus *theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
   const EcalSeverityLevelAlgo *sevlv = &iSetup.getData(tok_sevlv_);
 
   // Retrieve trigger tower map
   const EcalTrigTowerConstituentsMap &ttMap = iSetup.getData(tok_htmap_);
 
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
   if (!hbhe.isValid()) {
     ++nbad_;
     if (nbad_ < 10)
       edm::LogVerbatim("IsoTrack") << "No HBHE rechit collection";
     return;
   }
   const HBHERecHitCollection Hithbhe = *(hbhe.product());
 
   //get Handles to SimTracks and SimHits
   edm::Handle<edm::SimTrackContainer> SimTk;
   if (doMC_)
     iEvent.getByToken(tok_simTk_, SimTk);
 
   edm::Handle<edm::SimVertexContainer> SimVtx;
   if (doMC_)
     iEvent.getByToken(tok_simVtx_, SimVtx);
 
   //get Handles to PCaloHitContainers of eb/ee/hbhe
   edm::Handle<edm::PCaloHitContainer> pcaloeb;
   if (doMC_)
     iEvent.getByToken(tok_caloEB_, pcaloeb);
 
   edm::Handle<edm::PCaloHitContainer> pcaloee;
   if (doMC_)
     iEvent.getByToken(tok_caloEE_, pcaloee);
 
   edm::Handle<edm::PCaloHitContainer> pcalohh;
   if (doMC_)
     iEvent.getByToken(tok_caloHH_, pcalohh);
 
   //associates tracker rechits/simhits to a track
   std::unique_ptr<TrackerHitAssociator> associate;
   if (doMC_)
     associate = std::make_unique<TrackerHitAssociator>(iEvent, trackerHitAssociatorConfig_);
 
   //===================================================================================
 
   h_nTracks->Fill(trkCollection->size());
 
   int nTracks = 0;
 
   t_nTracks = trkCollection->size();
 
   // get the list of DetIds closest to the impact point of track on surface calorimeters
   std::vector<spr::propagatedTrackID> trkCaloDets;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDets, false);
   std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
 
   if (myverbose_ > 2) {
     for (trkDetItr = trkCaloDets.begin(); trkDetItr != trkCaloDets.end(); trkDetItr++) {
       edm::LogVerbatim("IsoTrack") << trkDetItr->trkItr->p() << " " << trkDetItr->trkItr->eta() << " "
                                    << trkDetItr->okECAL << " " << trkDetItr->okHCAL;
       if (trkDetItr->okECAL) {
         if (trkDetItr->detIdECAL.subdetId() == EcalBarrel)
           edm::LogVerbatim("IsoTrack") << (EBDetId)trkDetItr->detIdECAL;
         else
           edm::LogVerbatim("IsoTrack") << (EEDetId)trkDetItr->detIdECAL;
       }
       if (trkDetItr->okHCAL)
         edm::LogVerbatim("IsoTrack") << (HcalDetId)trkDetItr->detIdHCAL;
     }
   }
 
   int nvtxTracks = 0;
   for (trkDetItr = trkCaloDets.begin(), nTracks = 0; trkDetItr != trkCaloDets.end(); trkDetItr++, nTracks++) {
     const reco::Track *pTrack = &(*(trkDetItr->trkItr));
 
     // find vertex index the track is associated with
     int pVtxTkId = -1;
     for (unsigned int ind = 0; ind < recVtxs->size(); ind++) {
       if (!((*recVtxs)[ind].isFake())) {
         reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
         for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack != (*recVtxs)[ind].tracks_end(); vtxTrack++) {
           const edm::RefToBase<reco::Track> pvtxTrack = (*vtxTrack);
           if (pTrack == pvtxTrack.get()) {
             pVtxTkId = ind;
             break;
             if (myverbose_ > 2) {
               if (pTrack->pt() > 1.0) {
                 edm::LogVerbatim("IsoTrack") << "Debug the track association with vertex ";
                 edm::LogVerbatim("IsoTrack") << pTrack << " " << pvtxTrack.get();
                 edm::LogVerbatim("IsoTrack")
                     << " trkVtxIndex " << nvtxTracks << " vtx " << ind << " pt " << pTrack->pt() << " eta "
                     << pTrack->eta() << " " << pTrack->pt() - pvtxTrack->pt() << " "
                     << pTrack->eta() - pvtxTrack->eta();
                 nvtxTracks++;
               }
             }
           }
         }
       }
     }
 
     const reco::HitPattern &hitp = pTrack->hitPattern();
     int nLayersCrossed = hitp.trackerLayersWithMeasurement();
     int nOuterHits = hitp.stripTOBLayersWithMeasurement() + hitp.stripTECLayersWithMeasurement();
 
     bool ifGood = pTrack->quality(trackQuality_);
     double pt1 = pTrack->pt();
     double p1 = pTrack->p();
     double eta1 = pTrack->momentum().eta();
     double phi1 = pTrack->momentum().phi();
     double etaEcal1 = trkDetItr->etaECAL;
     double phiEcal1 = trkDetItr->phiECAL;
     double etaHcal1 = trkDetItr->etaHCAL;
     double phiHcal1 = trkDetItr->phiHCAL;
     double dxy1 = pTrack->dxy();
     double dz1 = pTrack->dz();
     double chisq1 = pTrack->normalizedChi2();
     double dxybs1 = beamSpotH.isValid() ? pTrack->dxy(bspot) : pTrack->dxy();
     double dzbs1 = beamSpotH.isValid() ? pTrack->dz(bspot) : pTrack->dz();
     double dxypv1 = pTrack->dxy();
     double dzpv1 = pTrack->dz();
     if (pVtxTkId >= 0) {
       math::XYZPoint thisTkPV =
           math::XYZPoint((*recVtxs)[pVtxTkId].x(), (*recVtxs)[pVtxTkId].y(), (*recVtxs)[pVtxTkId].z());
       dxypv1 = pTrack->dxy(thisTkPV);
       dzpv1 = pTrack->dz(thisTkPV);
     }
 
     h_recEtaPt_0->Fill(eta1, pt1);
     h_recEtaP_0->Fill(eta1, p1);
     h_recPt_0->Fill(pt1);
     h_recP_0->Fill(p1);
     h_recEta_0->Fill(eta1);
     h_recPhi_0->Fill(phi1);
 
     if (ifGood && nLayersCrossed > 7) {
       h_recEtaPt_1->Fill(eta1, pt1);
       h_recEtaP_1->Fill(eta1, p1);
       h_recPt_1->Fill(pt1);
       h_recP_1->Fill(p1);
       h_recEta_1->Fill(eta1);
       h_recPhi_1->Fill(phi1);
     }
 
     if (!ifGood)
       continue;
 
     if (writeAllTracks_ && p1 > 2.0 && nLayersCrossed > 7) {
       t_trackPAll->push_back(p1);
       t_trackEtaAll->push_back(eta1);
       t_trackPhiAll->push_back(phi1);
       t_trackPtAll->push_back(pt1);
       t_trackDxyAll->push_back(dxy1);
       t_trackDzAll->push_back(dz1);
       t_trackDxyPVAll->push_back(dxypv1);
       t_trackDzPVAll->push_back(dzpv1);
       t_trackChiSqAll->push_back(chisq1);
     }
     if (doMC_) {
       edm::SimTrackContainer::const_iterator matchedSimTrkAll =
           spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
       if (writeAllTracks_ && matchedSimTrkAll != SimTk->end())
         t_trackPdgIdAll->push_back(matchedSimTrkAll->type());
     }
 
     if (pt1 > minTrackP_ && std::abs(eta1) < maxTrackEta_ && trkDetItr->okECAL) {
       double maxNearP31x31 = 999.0, maxNearP25x25 = 999.0, maxNearP21x21 = 999.0, maxNearP15x15 = 999.0;
       maxNearP31x31 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15, 15);
       maxNearP25x25 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 12, 12);
       maxNearP21x21 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 10, 10);
       maxNearP15x15 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 7, 7);
 
       int iTrkEtaBin = -1, iTrkMomBin = -1;
       for (unsigned int ieta = 0; ieta < NEtaBins; ieta++) {
         if (std::abs(eta1) > genPartEtaBins[ieta] && std::abs(eta1) < genPartEtaBins[ieta + 1])
           iTrkEtaBin = ieta;
       }
       for (unsigned int ipt = 0; ipt < NPBins; ipt++) {
         if (p1 > genPartPBins[ipt] && p1 < genPartPBins[ipt + 1])
           iTrkMomBin = ipt;
       }
       if (iTrkMomBin >= 0 && iTrkEtaBin >= 0) {
         h_maxNearP31x31[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP31x31);
         h_maxNearP25x25[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP25x25);
         h_maxNearP21x21[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP21x21);
         h_maxNearP15x15[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP15x15);
       }
       if (maxNearP31x31 < 0.0 && nLayersCrossed > 7 && nOuterHits > 4) {
         h_recEtaPt_2->Fill(eta1, pt1);
         h_recEtaP_2->Fill(eta1, p1);
         h_recPt_2->Fill(pt1);
         h_recP_2->Fill(p1);
         h_recEta_2->Fill(eta1);
         h_recPhi_2->Fill(phi1);
       }
 
       // if charge isolated in ECAL, store the further quantities
       if (maxNearP31x31 < 1.0) {
         haveIsoTrack = true;
 
         // get the matching simTrack
         double simTrackP = -1;
         if (doMC_) {
           edm::SimTrackContainer::const_iterator matchedSimTrk =
               spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
           if (matchedSimTrk != SimTk->end())
             simTrackP = matchedSimTrk->momentum().P();
         }
 
         // get ECal Tranverse Profile
         std::pair<double, bool> e7x7P, e9x9P, e11x11P, e15x15P;
         std::pair<double, bool> e7x7_10SigP, e9x9_10SigP, e11x11_10SigP, e15x15_10SigP;
         std::pair<double, bool> e7x7_15SigP, e9x9_15SigP, e11x11_15SigP, e15x15_15SigP;
         std::pair<double, bool> e7x7_20SigP, e9x9_20SigP, e11x11_20SigP, e15x15_20SigP;
         std::pair<double, bool> e7x7_25SigP, e9x9_25SigP, e11x11_25SigP, e15x15_25SigP;
         std::pair<double, bool> e7x7_30SigP, e9x9_30SigP, e11x11_30SigP, e15x15_30SigP;
 
         spr::caloSimInfo simInfo3x3, simInfo5x5, simInfo7x7, simInfo9x9;
         spr::caloSimInfo simInfo11x11, simInfo13x13, simInfo15x15, simInfo21x21, simInfo25x25, simInfo31x31;
         double trkEcalEne = 0;
 
         const DetId isoCell = trkDetItr->detIdECAL;
         e7x7P = spr::eECALmatrix(isoCell,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  *theEcalChStatus,
                                  geo,
                                  caloTopology,
                                  sevlv,
                                  3,
                                  3,
                                  -100.0,
                                  -100.0,
                                  tMinE_,
                                  tMaxE_);
         e9x9P = spr::eECALmatrix(isoCell,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  *theEcalChStatus,
                                  geo,
                                  caloTopology,
                                  sevlv,
                                  4,
                                  4,
                                  -100.0,
                                  -100.0,
                                  tMinE_,
                                  tMaxE_);
         e11x11P = spr::eECALmatrix(isoCell,
                                    barrelRecHitsHandle,
                                    endcapRecHitsHandle,
                                    *theEcalChStatus,
                                    geo,
                                    caloTopology,
                                    sevlv,
                                    5,
                                    5,
                                    -100.0,
                                    -100.0,
                                    tMinE_,
                                    tMaxE_);
         e15x15P = spr::eECALmatrix(isoCell,
                                    barrelRecHitsHandle,
                                    endcapRecHitsHandle,
                                    *theEcalChStatus,
                                    geo,
                                    caloTopology,
                                    sevlv,
                                    7,
                                    7,
                                    -100.0,
                                    -100.0,
                                    tMinE_,
                                    tMaxE_);
 
         e7x7_10SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        3,
                                        3,
                                        0.030,
                                        0.150,
                                        tMinE_,
                                        tMaxE_);
         e9x9_10SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        4,
                                        4,
                                        0.030,
                                        0.150,
                                        tMinE_,
                                        tMaxE_);
         e11x11_10SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          5,
                                          5,
                                          0.030,
                                          0.150,
                                          tMinE_,
                                          tMaxE_);
         e15x15_10SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          7,
                                          7,
                                          0.030,
                                          0.150,
                                          tMinE_,
                                          tMaxE_);
 
         e7x7_15SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        ttMap,
                                        3,
                                        3,
                                        0.20,
                                        0.45,
                                        tMinE_,
                                        tMaxE_);
         e9x9_15SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        ttMap,
                                        4,
                                        4,
                                        0.20,
                                        0.45,
                                        tMinE_,
                                        tMaxE_);
         e11x11_15SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          ttMap,
                                          5,
                                          5,
                                          0.20,
                                          0.45,
                                          tMinE_,
                                          tMaxE_);
         e15x15_15SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          ttMap,
                                          7,
                                          7,
                                          0.20,
                                          0.45,
                                          tMinE_,
                                          tMaxE_,
                                          false);
 
         e7x7_20SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        3,
                                        3,
                                        0.060,
                                        0.300,
                                        tMinE_,
                                        tMaxE_);
         e9x9_20SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        4,
                                        4,
                                        0.060,
                                        0.300,
                                        tMinE_,
                                        tMaxE_);
         e11x11_20SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          5,
                                          5,
                                          0.060,
                                          0.300,
                                          tMinE_,
                                          tMaxE_);
         e15x15_20SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          7,
                                          7,
                                          0.060,
                                          0.300,
                                          tMinE_,
                                          tMaxE_);
 
         e7x7_25SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        3,
                                        3,
                                        0.075,
                                        0.375,
                                        tMinE_,
                                        tMaxE_);
         e9x9_25SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        4,
                                        4,
                                        0.075,
                                        0.375,
                                        tMinE_,
                                        tMaxE_);
         e11x11_25SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          5,
                                          5,
                                          0.075,
                                          0.375,
                                          tMinE_,
                                          tMaxE_);
         e15x15_25SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          7,
                                          7,
                                          0.075,
                                          0.375,
                                          tMinE_,
                                          tMaxE_);
 
         e7x7_30SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        3,
                                        3,
                                        0.090,
                                        0.450,
                                        tMinE_,
                                        tMaxE_);
         e9x9_30SigP = spr::eECALmatrix(isoCell,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        *theEcalChStatus,
                                        geo,
                                        caloTopology,
                                        sevlv,
                                        4,
                                        4,
                                        0.090,
                                        0.450,
                                        tMinE_,
                                        tMaxE_);
         e11x11_30SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          5,
                                          5,
                                          0.090,
                                          0.450,
                                          tMinE_,
                                          tMaxE_);
         e15x15_30SigP = spr::eECALmatrix(isoCell,
                                          barrelRecHitsHandle,
                                          endcapRecHitsHandle,
                                          *theEcalChStatus,
                                          geo,
                                          caloTopology,
                                          sevlv,
                                          7,
                                          7,
                                          0.090,
                                          0.450,
                                          tMinE_,
                                          tMaxE_);
         if (myverbose_ == 2) {
           edm::LogVerbatim("IsoTrack") << "clean  ecal rechit ";
           edm::LogVerbatim("IsoTrack") << "e7x7 " << e7x7P.first << " e9x9 " << e9x9P.first << " e11x11 "
                                        << e11x11P.first << " e15x15 " << e15x15P.first;
           edm::LogVerbatim("IsoTrack") << "e7x7_10Sig " << e7x7_10SigP.first << " e11x11_10Sig " << e11x11_10SigP.first
                                        << " e15x15_10Sig " << e15x15_10SigP.first;
         }
 
         if (doMC_) {
           // check the energy from SimHits
           spr::eECALSimInfo(
               iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 1, 1, simInfo3x3);
           spr::eECALSimInfo(
               iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 2, 2, simInfo5x5);
           spr::eECALSimInfo(
               iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 3, 3, simInfo7x7);
           spr::eECALSimInfo(
               iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 4, 4, simInfo9x9);
           spr::eECALSimInfo(iEvent,
                             isoCell,
                             geo,
                             caloTopology,
                             pcaloeb,
                             pcaloee,
                             SimTk,
                             SimVtx,
                             pTrack,
                             *associate,
                             5,
                             5,
                             simInfo11x11);
           spr::eECALSimInfo(iEvent,
                             isoCell,
                             geo,
                             caloTopology,
                             pcaloeb,
                             pcaloee,
                             SimTk,
                             SimVtx,
                             pTrack,
                             *associate,
                             6,
                             6,
                             simInfo13x13);
           spr::eECALSimInfo(iEvent,
                             isoCell,
                             geo,
                             caloTopology,
                             pcaloeb,
                             pcaloee,
                             SimTk,
                             SimVtx,
                             pTrack,
                             *associate,
                             7,
                             7,
                             simInfo15x15,
                             150.0,
                             false);
           spr::eECALSimInfo(iEvent,
                             isoCell,
                             geo,
                             caloTopology,
                             pcaloeb,
                             pcaloee,
                             SimTk,
                             SimVtx,
                             pTrack,
                             *associate,
                             10,
                             10,
                             simInfo21x21);
           spr::eECALSimInfo(iEvent,
                             isoCell,
                             geo,
                             caloTopology,
                             pcaloeb,
                             pcaloee,
                             SimTk,
                             SimVtx,
                             pTrack,
                             *associate,
                             12,
                             12,
                             simInfo25x25);
           spr::eECALSimInfo(iEvent,
                             isoCell,
                             geo,
                             caloTopology,
                             pcaloeb,
                             pcaloee,
                             SimTk,
                             SimVtx,
                             pTrack,
                             *associate,
                             15,
                             15,
                             simInfo31x31);
 
           trkEcalEne = spr::eCaloSimInfo(iEvent, geo, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate);
           if (myverbose_ == 1) {
             edm::LogVerbatim("IsoTrack") << "Track momentum " << pt1;
             edm::LogVerbatim("IsoTrack") << "ecal siminfo ";
             edm::LogVerbatim("IsoTrack") << "simInfo3x3: eTotal " << simInfo3x3.eTotal << " eMatched "
                                          << simInfo3x3.eMatched << " eRest " << simInfo3x3.eRest << " eGamma "
                                          << simInfo3x3.eGamma << " eNeutralHad " << simInfo3x3.eNeutralHad
                                          << " eChargedHad " << simInfo3x3.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "simInfo5x5: eTotal " << simInfo5x5.eTotal << " eMatched "
                                          << simInfo5x5.eMatched << " eRest " << simInfo5x5.eRest << " eGamma "
                                          << simInfo5x5.eGamma << " eNeutralHad " << simInfo5x5.eNeutralHad
                                          << " eChargedHad " << simInfo5x5.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "simInfo7x7: eTotal " << simInfo7x7.eTotal << " eMatched "
                                          << simInfo7x7.eMatched << " eRest " << simInfo7x7.eRest << " eGamma "
                                          << simInfo7x7.eGamma << " eNeutralHad " << simInfo7x7.eNeutralHad
                                          << " eChargedHad " << simInfo7x7.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "simInfo9x9: eTotal " << simInfo9x9.eTotal << " eMatched "
                                          << simInfo9x9.eMatched << " eRest " << simInfo9x9.eRest << " eGamma "
                                          << simInfo9x9.eGamma << " eNeutralHad " << simInfo9x9.eNeutralHad
                                          << " eChargedHad " << simInfo9x9.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "simInfo11x11: eTotal " << simInfo11x11.eTotal << " eMatched "
                                          << simInfo11x11.eMatched << " eRest " << simInfo11x11.eRest << " eGamma "
                                          << simInfo11x11.eGamma << " eNeutralHad " << simInfo11x11.eNeutralHad
                                          << " eChargedHad " << simInfo11x11.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "simInfo15x15: eTotal " << simInfo15x15.eTotal << " eMatched "
                                          << simInfo15x15.eMatched << " eRest " << simInfo15x15.eRest << " eGamma "
                                          << simInfo15x15.eGamma << " eNeutralHad " << simInfo15x15.eNeutralHad
                                          << " eChargedHad " << simInfo15x15.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "simInfo31x31: eTotal " << simInfo31x31.eTotal << " eMatched "
                                          << simInfo31x31.eMatched << " eRest " << simInfo31x31.eRest << " eGamma "
                                          << simInfo31x31.eGamma << " eNeutralHad " << simInfo31x31.eNeutralHad
                                          << " eChargedHad " << simInfo31x31.eChargedHad;
             edm::LogVerbatim("IsoTrack") << "trkEcalEne" << trkEcalEne;
           }
         }
 
         // =======  Get HCAL information
         double hcalScale = 1.0;
         if (std::abs(pTrack->eta()) < 1.4) {
           hcalScale = 120.0;
         } else {
           hcalScale = 135.0;
         }
 
         double maxNearHcalP3x3 = -1, maxNearHcalP5x5 = -1, maxNearHcalP7x7 = -1;
         maxNearHcalP3x3 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 1, 1);
         maxNearHcalP5x5 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 2, 2);
         maxNearHcalP7x7 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 3, 3);
 
         double h3x3 = 0, h5x5 = 0, h7x7 = 0;
         double h3x3Sig = 0, h5x5Sig = 0, h7x7Sig = 0;
         double trkHcalEne = 0;
         spr::caloSimInfo hsimInfo3x3, hsimInfo5x5, hsimInfo7x7;
 
         if (trkDetItr->okHCAL) {
           const DetId ClosestCell(trkDetItr->detIdHCAL);
           // bool includeHO=false, bool algoNew=true, bool debug=false
           h3x3 = spr::eHCALmatrix(
               theHBHETopology, ClosestCell, hbhe, 1, 1, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_, tMaxH_);
           h5x5 = spr::eHCALmatrix(
               theHBHETopology, ClosestCell, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_, tMaxH_);
           h7x7 = spr::eHCALmatrix(
               theHBHETopology, ClosestCell, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_, tMaxH_);
           h3x3Sig = spr::eHCALmatrix(
               theHBHETopology, ClosestCell, hbhe, 1, 1, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_, tMaxH_);
           h5x5Sig = spr::eHCALmatrix(
               theHBHETopology, ClosestCell, hbhe, 2, 2, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_, tMaxH_);
           h7x7Sig = spr::eHCALmatrix(
               theHBHETopology, ClosestCell, hbhe, 3, 3, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_, tMaxH_);
 
           if (myverbose_ == 2) {
             edm::LogVerbatim("IsoTrack") << "HCAL 3x3 " << h3x3 << " " << h3x3Sig << " 5x5 " << h5x5 << " " << h5x5Sig
                                          << " 7x7 " << h7x7 << " " << h7x7Sig;
           }
 
           if (doMC_) {
             spr::eHCALSimInfo(
                 iEvent, theHBHETopology, ClosestCell, geo, pcalohh, SimTk, SimVtx, pTrack, *associate, 1, 1, hsimInfo3x3);
             spr::eHCALSimInfo(
                 iEvent, theHBHETopology, ClosestCell, geo, pcalohh, SimTk, SimVtx, pTrack, *associate, 2, 2, hsimInfo5x5);
             spr::eHCALSimInfo(iEvent,
                               theHBHETopology,
                               ClosestCell,
                               geo,
                               pcalohh,
                               SimTk,
                               SimVtx,
                               pTrack,
                               *associate,
                               3,
                               3,
                               hsimInfo7x7,
                               150.0,
                               false,
                               false);
             trkHcalEne = spr::eCaloSimInfo(iEvent, geo, pcalohh, SimTk, SimVtx, pTrack, *associate);
             if (myverbose_ == 1) {
               edm::LogVerbatim("IsoTrack") << "Hcal siminfo ";
               edm::LogVerbatim("IsoTrack")
                   << "hsimInfo3x3: eTotal " << hsimInfo3x3.eTotal << " eMatched " << hsimInfo3x3.eMatched << " eRest "
                   << hsimInfo3x3.eRest << " eGamma " << hsimInfo3x3.eGamma << " eNeutralHad " << hsimInfo3x3.eNeutralHad
                   << " eChargedHad " << hsimInfo3x3.eChargedHad;
               edm::LogVerbatim("IsoTrack")
                   << "hsimInfo5x5: eTotal " << hsimInfo5x5.eTotal << " eMatched " << hsimInfo5x5.eMatched << " eRest "
                   << hsimInfo5x5.eRest << " eGamma " << hsimInfo5x5.eGamma << " eNeutralHad " << hsimInfo5x5.eNeutralHad
                   << " eChargedHad " << hsimInfo5x5.eChargedHad;
               edm::LogVerbatim("IsoTrack")
                   << "hsimInfo7x7: eTotal " << hsimInfo7x7.eTotal << " eMatched " << hsimInfo7x7.eMatched << " eRest "
                   << hsimInfo7x7.eRest << " eGamma " << hsimInfo7x7.eGamma << " eNeutralHad " << hsimInfo7x7.eNeutralHad
                   << " eChargedHad " << hsimInfo7x7.eChargedHad;
               edm::LogVerbatim("IsoTrack") << "trkHcalEne " << trkHcalEne;
             }
           }
 
           // debug the ecal and hcal matrix
           if (myverbose_ == 4) {
             edm::LogVerbatim("IsoTrack") << "Run " << iEvent.id().run() << "  Event " << iEvent.id().event();
             std::vector<std::pair<DetId, double> > v7x7 =
                 spr::eHCALmatrixCell(theHBHETopology, ClosestCell, hbhe, 3, 3, false, false);
             double sumv = 0.0;
 
             for (unsigned int iv = 0; iv < v7x7.size(); iv++) {
               sumv += v7x7[iv].second;
             }
             edm::LogVerbatim("IsoTrack") << "h7x7 " << h7x7 << " v7x7 " << sumv << " in " << v7x7.size();
             for (unsigned int iv = 0; iv < v7x7.size(); iv++) {
               HcalDetId id = v7x7[iv].first;
               edm::LogVerbatim("IsoTrack") << " Cell " << iv << " 0x" << std::hex << v7x7[iv].first() << std::dec << " "
                                            << id << " Energy " << v7x7[iv].second;
             }
           }
         }
         if (doMC_) {
           trkHcalEne = spr::eCaloSimInfo(iEvent, geo, pcalohh, SimTk, SimVtx, pTrack, *associate);
         }
 
         // ====================================================================================================
         // get diff between track outermost hit position and the propagation point at outermost surface of tracker
         std::pair<math::XYZPoint, double> point2_TK0 = spr::propagateTrackerEnd(pTrack, bField, false);
         math::XYZPoint diff(pTrack->outerPosition().X() - point2_TK0.first.X(),
                             pTrack->outerPosition().Y() - point2_TK0.first.Y(),
                             pTrack->outerPosition().Z() - point2_TK0.first.Z());
         double trackOutPosOutHitDr = diff.R();
         double trackL = point2_TK0.second;
         if (myverbose_ == 5) {
           edm::LogVerbatim("IsoTrack") << " propagted " << point2_TK0.first << " " << point2_TK0.first.eta() << " "
                                        << point2_TK0.first.phi();
           edm::LogVerbatim("IsoTrack") << " outerPosition() " << pTrack->outerPosition() << " "
                                        << pTrack->outerPosition().eta() << " " << pTrack->outerPosition().phi();
           edm::LogVerbatim("IsoTrack") << "diff " << diff << " diffR " << diff.R() << " diffR/L "
                                        << diff.R() / point2_TK0.second;
         }
 
         for (unsigned int ind = 0; ind < recVtxs->size(); ind++) {
           if (!((*recVtxs)[ind].isFake())) {
             reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
             if (reco::deltaR(eta1, phi1, (*vtxTrack)->eta(), (*vtxTrack)->phi()) < 0.01)
               t_trackPVIdx->push_back(ind);
             else
               t_trackPVIdx->push_back(-1);
           }
         }
 
         // Fill the tree Branches here
         t_trackPVIdx->push_back(pVtxTkId);
         t_trackP->push_back(p1);
         t_trackPt->push_back(pt1);
         t_trackEta->push_back(eta1);
         t_trackPhi->push_back(phi1);
         t_trackEcalEta->push_back(etaEcal1);
         t_trackEcalPhi->push_back(phiEcal1);
         t_trackHcalEta->push_back(etaHcal1);
         t_trackHcalPhi->push_back(phiHcal1);
         t_trackDxy->push_back(dxy1);
         t_trackDz->push_back(dz1);
         t_trackDxyBS->push_back(dxybs1);
         t_trackDzBS->push_back(dzbs1);
         t_trackDxyPV->push_back(dxypv1);
         t_trackDzPV->push_back(dzpv1);
         t_trackChiSq->push_back(chisq1);
         t_trackNOuterHits->push_back(nOuterHits);
         t_NLayersCrossed->push_back(nLayersCrossed);
 
         t_trackHitsTOB->push_back(hitp.stripTOBLayersWithMeasurement());
         t_trackHitsTEC->push_back(hitp.stripTECLayersWithMeasurement());
         t_trackHitInMissTOB->push_back(hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
         t_trackHitInMissTEC->push_back(hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
         t_trackHitInMissTIB->push_back(hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
         t_trackHitInMissTID->push_back(hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
         t_trackHitInMissTIBTID->push_back(hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) +
                                           hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
 
         t_trackHitOutMissTOB->push_back(hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
         t_trackHitOutMissTEC->push_back(hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
         t_trackHitOutMissTIB->push_back(hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
         t_trackHitOutMissTID->push_back(hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
         t_trackHitOutMissTOBTEC->push_back(hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) +
                                            hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
 
         t_trackHitInMeasTOB->push_back(hitp.stripTOBLayersWithMeasurement());
         t_trackHitInMeasTEC->push_back(hitp.stripTECLayersWithMeasurement());
         t_trackHitInMeasTIB->push_back(hitp.stripTIBLayersWithMeasurement());
         t_trackHitInMeasTID->push_back(hitp.stripTIDLayersWithMeasurement());
         t_trackHitOutMeasTOB->push_back(hitp.stripTOBLayersWithMeasurement());
         t_trackHitOutMeasTEC->push_back(hitp.stripTECLayersWithMeasurement());
         t_trackHitOutMeasTIB->push_back(hitp.stripTIBLayersWithMeasurement());
         t_trackHitOutMeasTID->push_back(hitp.stripTIDLayersWithMeasurement());
         t_trackOutPosOutHitDr->push_back(trackOutPosOutHitDr);
         t_trackL->push_back(trackL);
 
         t_maxNearP31x31->push_back(maxNearP31x31);
         t_maxNearP21x21->push_back(maxNearP21x21);
 
         t_ecalSpike11x11->push_back(e11x11P.second);
         t_e7x7->push_back(e7x7P.first);
         t_e9x9->push_back(e9x9P.first);
         t_e11x11->push_back(e11x11P.first);
         t_e15x15->push_back(e15x15P.first);
 
         t_e7x7_10Sig->push_back(e7x7_10SigP.first);
         t_e9x9_10Sig->push_back(e9x9_10SigP.first);
         t_e11x11_10Sig->push_back(e11x11_10SigP.first);
         t_e15x15_10Sig->push_back(e15x15_10SigP.first);
         t_e7x7_15Sig->push_back(e7x7_15SigP.first);
         t_e9x9_15Sig->push_back(e9x9_15SigP.first);
         t_e11x11_15Sig->push_back(e11x11_15SigP.first);
         t_e15x15_15Sig->push_back(e15x15_15SigP.first);
         t_e7x7_20Sig->push_back(e7x7_20SigP.first);
         t_e9x9_20Sig->push_back(e9x9_20SigP.first);
         t_e11x11_20Sig->push_back(e11x11_20SigP.first);
         t_e15x15_20Sig->push_back(e15x15_20SigP.first);
         t_e7x7_25Sig->push_back(e7x7_25SigP.first);
         t_e9x9_25Sig->push_back(e9x9_25SigP.first);
         t_e11x11_25Sig->push_back(e11x11_25SigP.first);
         t_e15x15_25Sig->push_back(e15x15_25SigP.first);
         t_e7x7_30Sig->push_back(e7x7_30SigP.first);
         t_e9x9_30Sig->push_back(e9x9_30SigP.first);
         t_e11x11_30Sig->push_back(e11x11_30SigP.first);
         t_e15x15_30Sig->push_back(e15x15_30SigP.first);
 
         if (doMC_) {
           t_esim7x7->push_back(simInfo7x7.eTotal);
           t_esim9x9->push_back(simInfo9x9.eTotal);
           t_esim11x11->push_back(simInfo11x11.eTotal);
           t_esim15x15->push_back(simInfo15x15.eTotal);
 
           t_esim7x7Matched->push_back(simInfo7x7.eMatched);
           t_esim9x9Matched->push_back(simInfo9x9.eMatched);
           t_esim11x11Matched->push_back(simInfo11x11.eMatched);
           t_esim15x15Matched->push_back(simInfo15x15.eMatched);
 
           t_esim7x7Rest->push_back(simInfo7x7.eRest);
           t_esim9x9Rest->push_back(simInfo9x9.eRest);
           t_esim11x11Rest->push_back(simInfo11x11.eRest);
           t_esim15x15Rest->push_back(simInfo15x15.eRest);
 
           t_esim7x7Photon->push_back(simInfo7x7.eGamma);
           t_esim9x9Photon->push_back(simInfo9x9.eGamma);
           t_esim11x11Photon->push_back(simInfo11x11.eGamma);
           t_esim15x15Photon->push_back(simInfo15x15.eGamma);
 
           t_esim7x7NeutHad->push_back(simInfo7x7.eNeutralHad);
           t_esim9x9NeutHad->push_back(simInfo9x9.eNeutralHad);
           t_esim11x11NeutHad->push_back(simInfo11x11.eNeutralHad);
           t_esim15x15NeutHad->push_back(simInfo15x15.eNeutralHad);
 
           t_esim7x7CharHad->push_back(simInfo7x7.eChargedHad);
           t_esim9x9CharHad->push_back(simInfo9x9.eChargedHad);
           t_esim11x11CharHad->push_back(simInfo11x11.eChargedHad);
           t_esim15x15CharHad->push_back(simInfo15x15.eChargedHad);
 
           t_trkEcalEne->push_back(trkEcalEne);
           t_simTrackP->push_back(simTrackP);
           t_esimPdgId->push_back(simInfo11x11.pdgMatched);
         }
 
         t_maxNearHcalP3x3->push_back(maxNearHcalP3x3);
         t_maxNearHcalP5x5->push_back(maxNearHcalP5x5);
         t_maxNearHcalP7x7->push_back(maxNearHcalP7x7);
 
         t_h3x3->push_back(h3x3);
         t_h5x5->push_back(h5x5);
         t_h7x7->push_back(h7x7);
         t_h3x3Sig->push_back(h3x3Sig);
         t_h5x5Sig->push_back(h5x5Sig);
         t_h7x7Sig->push_back(h7x7Sig);
 
         t_infoHcal->push_back(trkDetItr->okHCAL);
         if (doMC_) {
           t_trkHcalEne->push_back(hcalScale * trkHcalEne);
 
           t_hsim3x3->push_back(hcalScale * hsimInfo3x3.eTotal);
           t_hsim5x5->push_back(hcalScale * hsimInfo5x5.eTotal);
           t_hsim7x7->push_back(hcalScale * hsimInfo7x7.eTotal);
 
           t_hsim3x3Matched->push_back(hcalScale * hsimInfo3x3.eMatched);
           t_hsim5x5Matched->push_back(hcalScale * hsimInfo5x5.eMatched);
           t_hsim7x7Matched->push_back(hcalScale * hsimInfo7x7.eMatched);
 
           t_hsim3x3Rest->push_back(hcalScale * hsimInfo3x3.eRest);
           t_hsim5x5Rest->push_back(hcalScale * hsimInfo5x5.eRest);
           t_hsim7x7Rest->push_back(hcalScale * hsimInfo7x7.eRest);
 
           t_hsim3x3Photon->push_back(hcalScale * hsimInfo3x3.eGamma);
           t_hsim5x5Photon->push_back(hcalScale * hsimInfo5x5.eGamma);
           t_hsim7x7Photon->push_back(hcalScale * hsimInfo7x7.eGamma);
 
           t_hsim3x3NeutHad->push_back(hcalScale * hsimInfo3x3.eNeutralHad);
           t_hsim5x5NeutHad->push_back(hcalScale * hsimInfo5x5.eNeutralHad);
           t_hsim7x7NeutHad->push_back(hcalScale * hsimInfo7x7.eNeutralHad);
 
           t_hsim3x3CharHad->push_back(hcalScale * hsimInfo3x3.eChargedHad);
           t_hsim5x5CharHad->push_back(hcalScale * hsimInfo5x5.eChargedHad);
           t_hsim7x7CharHad->push_back(hcalScale * hsimInfo7x7.eChargedHad);
         }
 
       }  // if loosely isolated track
     }    // check p1/eta
   }      // loop over track collection
 
   if (haveIsoTrack)
     tree_->Fill();
 }

◆ beginJob()

void IsolatedTracksNxN::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 2022 of file IsolatedTracksNxN.cc.

References bookHistograms(), genPartEtaBins, genPartPBins, mps_fire::i, initL1_, nbad_, NEtaBins, nEventProc_, and NPBins.

                                  {
   nEventProc_ = 0;
   nbad_ = 0;
   initL1_ = false;
   double tempgen_TH[NPBins + 1] = {
       0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 9.0, 11.0, 15.0, 20.0, 30.0, 50.0, 75.0, 100.0};
 
   for (unsigned int i = 0; i < NPBins + 1; i++)
     genPartPBins[i] = tempgen_TH[i];
 
   double tempgen_Eta[NEtaBins + 1] = {0.0, 1.131, 1.653, 2.172};
 
   for (unsigned int i = 0; i < NEtaBins + 1; i++)
     genPartEtaBins[i] = tempgen_Eta[i];
 
   bookHistograms();
 }

◆ bookHistograms()

void IsolatedTracksNxN::bookHistograms ( )

private

Definition at line 2260 of file IsolatedTracksNxN.cc.

References DeadROC_duringRun::dir, doMC_, compareTotals::fs, genPartEtaBins, genPartPBins, h_L1AlgoNames, h_maxNearP15x15, h_maxNearP21x21, h_maxNearP25x25, h_maxNearP31x31, h_nTracks, h_PVTracksWt, h_recEta_0, h_recEta_1, h_recEta_2, h_recEtaP_0, h_recEtaP_1, h_recEtaP_2, h_recEtaPt_0, h_recEtaPt_1, h_recEtaPt_2, h_recP_0, h_recP_1, h_recP_2, h_recPhi_0, h_recPhi_1, h_recPhi_2, h_recPt_0, h_recPt_1, h_recPt_2, mps_fire::i, LEDCalibrationChannels::ieta, NEtaBins, NPBins, t_Bunch, t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_EvtNo, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1AlgoNames, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1Decision, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1METEta, t_L1METPhi, t_L1METPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1PreScale, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_Lumi, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP21x21, t_maxNearP31x31, t_NLayersCrossed, t_nTracks, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_RunNo, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, t_trkHcalEne, and tree_.

Referenced by beginJob().

                                        {
   char hname[100], htit[100];
 
   edm::Service<TFileService> fs;
   TFileDirectory dir = fs->mkdir("nearMaxTrackP");
 
   for (unsigned int ieta = 0; ieta < NEtaBins; ieta++) {
     double lowEta = -5.0, highEta = 5.0;
     lowEta = genPartEtaBins[ieta];
     highEta = genPartEtaBins[ieta + 1];
 
     for (unsigned int ipt = 0; ipt < NPBins; ipt++) {
       double lowP = 0.0, highP = 300.0;
       lowP = genPartPBins[ipt];
       highP = genPartPBins[ipt + 1];
       sprintf(hname, "h_maxNearP31x31_ptBin%i_etaBin%i", ipt, ieta);
       sprintf(htit, "maxNearP in 31x31 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
       h_maxNearP31x31[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
       h_maxNearP31x31[ipt][ieta]->Sumw2();
       sprintf(hname, "h_maxNearP25x25_ptBin%i_etaBin%i", ipt, ieta);
       sprintf(htit, "maxNearP in 25x25 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
       h_maxNearP25x25[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
       h_maxNearP25x25[ipt][ieta]->Sumw2();
       sprintf(hname, "h_maxNearP21x21_ptBin%i_etaBin%i", ipt, ieta);
       sprintf(htit, "maxNearP in 21x21 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
       h_maxNearP21x21[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
       h_maxNearP21x21[ipt][ieta]->Sumw2();
       sprintf(hname, "h_maxNearP15x15_ptBin%i_etaBin%i", ipt, ieta);
       sprintf(htit, "maxNearP in 15x15 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
       h_maxNearP15x15[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
       h_maxNearP15x15[ipt][ieta]->Sumw2();
     }
   }
 
   h_L1AlgoNames = fs->make<TH1I>("h_L1AlgoNames", "h_L1AlgoNames:Bin Labels", 128, -0.5, 127.5);
 
   // Reconstructed Tracks
 
   h_PVTracksWt = fs->make<TH1F>("h_PVTracksWt", "h_PVTracksWt", 600, -0.1, 1.1);
 
   h_nTracks = fs->make<TH1F>("h_nTracks", "h_nTracks", 1000, -0.5, 999.5);
 
   sprintf(hname, "h_recEtaPt_0");
   sprintf(htit, "h_recEtaPt (all tracks Eta vs pT)");
   h_recEtaPt_0 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
 
   sprintf(hname, "h_recEtaP_0");
   sprintf(htit, "h_recEtaP (all tracks Eta vs pT)");
   h_recEtaP_0 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
 
   h_recPt_0 = fs->make<TH1F>("h_recPt_0", "Pt (all tracks)", 15, genPartPBins);
   h_recP_0 = fs->make<TH1F>("h_recP_0", "P  (all tracks)", 15, genPartPBins);
   h_recEta_0 = fs->make<TH1F>("h_recEta_0", "Eta (all tracks)", 60, -3.0, 3.0);
   h_recPhi_0 = fs->make<TH1F>("h_recPhi_0", "Phi (all tracks)", 100, -3.2, 3.2);
   //-------------------------
   sprintf(hname, "h_recEtaPt_1");
   sprintf(htit, "h_recEtaPt (all good tracks Eta vs pT)");
   h_recEtaPt_1 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
 
   sprintf(hname, "h_recEtaP_1");
   sprintf(htit, "h_recEtaP (all good tracks Eta vs pT)");
   h_recEtaP_1 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
 
   h_recPt_1 = fs->make<TH1F>("h_recPt_1", "Pt (all good tracks)", 15, genPartPBins);
   h_recP_1 = fs->make<TH1F>("h_recP_1", "P  (all good tracks)", 15, genPartPBins);
   h_recEta_1 = fs->make<TH1F>("h_recEta_1", "Eta (all good tracks)", 60, -3.0, 3.0);
   h_recPhi_1 = fs->make<TH1F>("h_recPhi_1", "Phi (all good tracks)", 100, -3.2, 3.2);
   //-------------------------
   sprintf(hname, "h_recEtaPt_2");
   sprintf(htit, "h_recEtaPt (charge isolation Eta vs pT)");
   h_recEtaPt_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
 
   sprintf(hname, "h_recEtaP_2");
   sprintf(htit, "h_recEtaP (charge isolation Eta vs pT)");
   h_recEtaP_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
 
   h_recPt_2 = fs->make<TH1F>("h_recPt_2", "Pt (charge isolation)", 15, genPartPBins);
   h_recP_2 = fs->make<TH1F>("h_recP_2", "P  (charge isolation)", 15, genPartPBins);
   h_recEta_2 = fs->make<TH1F>("h_recEta_2", "Eta (charge isolation)", 60, -3.0, 3.0);
   h_recPhi_2 = fs->make<TH1F>("h_recPhi_2", "Phi (charge isolation)", 100, -3.2, 3.2);
 
   tree_ = fs->make<TTree>("tree", "tree");
   tree_->SetAutoSave(10000);
 
   tree_->Branch("t_EvtNo", &t_EvtNo, "t_EvtNo/I");
   tree_->Branch("t_RunNo", &t_RunNo, "t_RunNo/I");
   tree_->Branch("t_Lumi", &t_Lumi, "t_Lumi/I");
   tree_->Branch("t_Bunch", &t_Bunch, "t_Bunch/I");
 
   t_PVx = new std::vector<double>();
   t_PVy = new std::vector<double>();
   t_PVz = new std::vector<double>();
   t_PVisValid = new std::vector<int>();
   t_PVndof = new std::vector<int>();
   t_PVNTracks = new std::vector<int>();
   t_PVNTracksWt = new std::vector<int>();
   t_PVTracksSumPt = new std::vector<double>();
   t_PVTracksSumPtWt = new std::vector<double>();
   t_PVNTracksHP = new std::vector<int>();
   t_PVNTracksHPWt = new std::vector<int>();
   t_PVTracksSumPtHP = new std::vector<double>();
   t_PVTracksSumPtHPWt = new std::vector<double>();
 
   tree_->Branch("PVx", "std::vector<double>", &t_PVx);
   tree_->Branch("PVy", "std::vector<double>", &t_PVy);
   tree_->Branch("PVz", "std::vector<double>", &t_PVz);
   tree_->Branch("PVisValid", "std::vector<int>", &t_PVisValid);
   tree_->Branch("PVndof", "std::vector<int>", &t_PVndof);
   tree_->Branch("PVNTracks", "std::vector<int>", &t_PVNTracks);
   tree_->Branch("PVNTracksWt", "std::vector<int>", &t_PVNTracksWt);
   tree_->Branch("t_PVTracksSumPt", "std::vector<double>", &t_PVTracksSumPt);
   tree_->Branch("t_PVTracksSumPtWt", "std::vector<double>", &t_PVTracksSumPtWt);
   tree_->Branch("PVNTracksHP", "std::vector<int>", &t_PVNTracksHP);
   tree_->Branch("PVNTracksHPWt", "std::vector<int>", &t_PVNTracksHPWt);
   tree_->Branch("t_PVTracksSumPtHP", "std::vector<double>", &t_PVTracksSumPtHP);
   tree_->Branch("t_PVTracksSumPtHPWt", "std::vector<double>", &t_PVTracksSumPtHPWt);
 
   //----- L1Trigger information
   for (int i = 0; i < 128; i++)
     t_L1Decision[i] = 0;
   t_L1AlgoNames = new std::vector<std::string>();
   t_L1PreScale = new std::vector<int>();
   t_L1CenJetPt = new std::vector<double>();
   t_L1CenJetEta = new std::vector<double>();
   t_L1CenJetPhi = new std::vector<double>();
   t_L1FwdJetPt = new std::vector<double>();
   t_L1FwdJetEta = new std::vector<double>();
   t_L1FwdJetPhi = new std::vector<double>();
   t_L1TauJetPt = new std::vector<double>();
   t_L1TauJetEta = new std::vector<double>();
   t_L1TauJetPhi = new std::vector<double>();
   t_L1MuonPt = new std::vector<double>();
   t_L1MuonEta = new std::vector<double>();
   t_L1MuonPhi = new std::vector<double>();
   t_L1IsoEMPt = new std::vector<double>();
   t_L1IsoEMEta = new std::vector<double>();
   t_L1IsoEMPhi = new std::vector<double>();
   t_L1NonIsoEMPt = new std::vector<double>();
   t_L1NonIsoEMEta = new std::vector<double>();
   t_L1NonIsoEMPhi = new std::vector<double>();
   t_L1METPt = new std::vector<double>();
   t_L1METEta = new std::vector<double>();
   t_L1METPhi = new std::vector<double>();
 
   tree_->Branch("t_L1Decision", t_L1Decision, "t_L1Decision[128]/I");
   tree_->Branch("t_L1AlgoNames", "std::vector<string>", &t_L1AlgoNames);
   tree_->Branch("t_L1PreScale", "std::vector<int>", &t_L1PreScale);
   tree_->Branch("t_L1CenJetPt", "std::vector<double>", &t_L1CenJetPt);
   tree_->Branch("t_L1CenJetEta", "std::vector<double>", &t_L1CenJetEta);
   tree_->Branch("t_L1CenJetPhi", "std::vector<double>", &t_L1CenJetPhi);
   tree_->Branch("t_L1FwdJetPt", "std::vector<double>", &t_L1FwdJetPt);
   tree_->Branch("t_L1FwdJetEta", "std::vector<double>", &t_L1FwdJetEta);
   tree_->Branch("t_L1FwdJetPhi", "std::vector<double>", &t_L1FwdJetPhi);
   tree_->Branch("t_L1TauJetPt", "std::vector<double>", &t_L1TauJetPt);
   tree_->Branch("t_L1TauJetEta", "std::vector<double>", &t_L1TauJetEta);
   tree_->Branch("t_L1TauJetPhi", "std::vector<double>", &t_L1TauJetPhi);
   tree_->Branch("t_L1MuonPt", "std::vector<double>", &t_L1MuonPt);
   tree_->Branch("t_L1MuonEta", "std::vector<double>", &t_L1MuonEta);
   tree_->Branch("t_L1MuonPhi", "std::vector<double>", &t_L1MuonPhi);
   tree_->Branch("t_L1IsoEMPt", "std::vector<double>", &t_L1IsoEMPt);
   tree_->Branch("t_L1IsoEMEta", "std::vector<double>", &t_L1IsoEMEta);
   tree_->Branch("t_L1IsoEMPhi", "std::vector<double>", &t_L1IsoEMPhi);
   tree_->Branch("t_L1NonIsoEMPt", "std::vector<double>", &t_L1NonIsoEMPt);
   tree_->Branch("t_L1NonIsoEMEta", "std::vector<double>", &t_L1NonIsoEMEta);
   tree_->Branch("t_L1NonIsoEMPhi", "std::vector<double>", &t_L1NonIsoEMPhi);
   tree_->Branch("t_L1METPt", "std::vector<double>", &t_L1METPt);
   tree_->Branch("t_L1METEta", "std::vector<double>", &t_L1METEta);
   tree_->Branch("t_L1METPhi", "std::vector<double>", &t_L1METPhi);
 
   t_jetPt = new std::vector<double>();
   t_jetEta = new std::vector<double>();
   t_jetPhi = new std::vector<double>();
   t_nTrksJetCalo = new std::vector<double>();
   t_nTrksJetVtx = new std::vector<double>();
   tree_->Branch("t_jetPt", "std::vector<double>", &t_jetPt);
   tree_->Branch("t_jetEta", "std::vector<double>", &t_jetEta);
   tree_->Branch("t_jetPhi", "std::vector<double>", &t_jetPhi);
   tree_->Branch("t_nTrksJetCalo", "std::vector<double>", &t_nTrksJetCalo);
   tree_->Branch("t_nTrksJetVtx", "std::vector<double>", &t_nTrksJetVtx);
 
   t_trackPAll = new std::vector<double>();
   t_trackEtaAll = new std::vector<double>();
   t_trackPhiAll = new std::vector<double>();
   t_trackPdgIdAll = new std::vector<double>();
   t_trackPtAll = new std::vector<double>();
   t_trackDxyAll = new std::vector<double>();
   t_trackDzAll = new std::vector<double>();
   t_trackDxyPVAll = new std::vector<double>();
   t_trackDzPVAll = new std::vector<double>();
   t_trackChiSqAll = new std::vector<double>();
   tree_->Branch("t_trackPAll", "std::vector<double>", &t_trackPAll);
   tree_->Branch("t_trackPhiAll", "std::vector<double>", &t_trackPhiAll);
   tree_->Branch("t_trackEtaAll", "std::vector<double>", &t_trackEtaAll);
   tree_->Branch("t_trackPtAll", "std::vector<double>", &t_trackPtAll);
   tree_->Branch("t_trackDxyAll", "std::vector<double>", &t_trackDxyAll);
   tree_->Branch("t_trackDzAll", "std::vector<double>", &t_trackDzAll);
   tree_->Branch("t_trackDxyPVAll", "std::vector<double>", &t_trackDxyPVAll);
   tree_->Branch("t_trackDzPVAll", "std::vector<double>", &t_trackDzPVAll);
   tree_->Branch("t_trackChiSqAll", "std::vector<double>", &t_trackChiSqAll);
   //tree_->Branch("t_trackPdgIdAll",     "std::vector<double>", &t_trackPdgIdAll);
 
   t_trackP = new std::vector<double>();
   t_trackPt = new std::vector<double>();
   t_trackEta = new std::vector<double>();
   t_trackPhi = new std::vector<double>();
   t_trackEcalEta = new std::vector<double>();
   t_trackEcalPhi = new std::vector<double>();
   t_trackHcalEta = new std::vector<double>();
   t_trackHcalPhi = new std::vector<double>();
   t_trackNOuterHits = new std::vector<int>();
   t_NLayersCrossed = new std::vector<int>();
   t_trackDxy = new std::vector<double>();
   t_trackDxyBS = new std::vector<double>();
   t_trackDz = new std::vector<double>();
   t_trackDzBS = new std::vector<double>();
   t_trackDxyPV = new std::vector<double>();
   t_trackDzPV = new std::vector<double>();
   t_trackPVIdx = new std::vector<int>();
   t_trackChiSq = new std::vector<double>();
   t_trackHitsTOB = new std::vector<int>();
   t_trackHitsTEC = new std::vector<int>();
   t_trackHitInMissTOB = new std::vector<int>();
   t_trackHitInMissTEC = new std::vector<int>();
   t_trackHitInMissTIB = new std::vector<int>();
   t_trackHitInMissTID = new std::vector<int>();
   t_trackHitOutMissTOB = new std::vector<int>();
   t_trackHitOutMissTEC = new std::vector<int>();
   t_trackHitOutMissTIB = new std::vector<int>();
   t_trackHitOutMissTID = new std::vector<int>();
   t_trackHitInMissTIBTID = new std::vector<int>();
   t_trackHitOutMissTOB = new std::vector<int>();
   t_trackHitOutMissTEC = new std::vector<int>();
   t_trackHitOutMissTIB = new std::vector<int>();
   t_trackHitOutMissTID = new std::vector<int>();
   t_trackHitOutMissTOBTEC = new std::vector<int>();
   t_trackHitInMeasTOB = new std::vector<int>();
   t_trackHitInMeasTEC = new std::vector<int>();
   t_trackHitInMeasTIB = new std::vector<int>();
   t_trackHitInMeasTID = new std::vector<int>();
   t_trackHitOutMeasTOB = new std::vector<int>();
   t_trackHitOutMeasTEC = new std::vector<int>();
   t_trackHitOutMeasTIB = new std::vector<int>();
   t_trackHitOutMeasTID = new std::vector<int>();
   t_trackOutPosOutHitDr = new std::vector<double>();
   t_trackL = new std::vector<double>();
 
   tree_->Branch("t_trackP", "std::vector<double>", &t_trackP);
   tree_->Branch("t_trackPt", "std::vector<double>", &t_trackPt);
   tree_->Branch("t_trackEta", "std::vector<double>", &t_trackEta);
   tree_->Branch("t_trackPhi", "std::vector<double>", &t_trackPhi);
   tree_->Branch("t_trackEcalEta", "std::vector<double>", &t_trackEcalEta);
   tree_->Branch("t_trackEcalPhi", "std::vector<double>", &t_trackEcalPhi);
   tree_->Branch("t_trackHcalEta", "std::vector<double>", &t_trackHcalEta);
   tree_->Branch("t_trackHcalPhi", "std::vector<double>", &t_trackHcalPhi);
 
   tree_->Branch("t_trackNOuterHits", "std::vector<int>", &t_trackNOuterHits);
   tree_->Branch("t_NLayersCrossed", "std::vector<int>", &t_NLayersCrossed);
   tree_->Branch("t_trackHitsTOB", "std::vector<int>", &t_trackHitsTOB);
   tree_->Branch("t_trackHitsTEC", "std::vector<int>", &t_trackHitsTEC);
   tree_->Branch("t_trackHitInMissTOB", "std::vector<int>", &t_trackHitInMissTOB);
   tree_->Branch("t_trackHitInMissTEC", "std::vector<int>", &t_trackHitInMissTEC);
   tree_->Branch("t_trackHitInMissTIB", "std::vector<int>", &t_trackHitInMissTIB);
   tree_->Branch("t_trackHitInMissTID", "std::vector<int>", &t_trackHitInMissTID);
   tree_->Branch("t_trackHitInMissTIBTID", "std::vector<int>", &t_trackHitInMissTIBTID);
   tree_->Branch("t_trackHitOutMissTOB", "std::vector<int>", &t_trackHitOutMissTOB);
   tree_->Branch("t_trackHitOutMissTEC", "std::vector<int>", &t_trackHitOutMissTEC);
   tree_->Branch("t_trackHitOutMissTIB", "std::vector<int>", &t_trackHitOutMissTIB);
   tree_->Branch("t_trackHitOutMissTID", "std::vector<int>", &t_trackHitOutMissTID);
   tree_->Branch("t_trackHitOutMissTOBTEC", "std::vector<int>", &t_trackHitOutMissTOBTEC);
   tree_->Branch("t_trackHitInMeasTOB", "std::vector<int>", &t_trackHitInMeasTOB);
   tree_->Branch("t_trackHitInMeasTEC", "std::vector<int>", &t_trackHitInMeasTEC);
   tree_->Branch("t_trackHitInMeasTIB", "std::vector<int>", &t_trackHitInMeasTIB);
   tree_->Branch("t_trackHitInMeasTID", "std::vector<int>", &t_trackHitInMeasTID);
   tree_->Branch("t_trackHitOutMeasTOB", "std::vector<int>", &t_trackHitOutMeasTOB);
   tree_->Branch("t_trackHitOutMeasTEC", "std::vector<int>", &t_trackHitOutMeasTEC);
   tree_->Branch("t_trackHitOutMeasTIB", "std::vector<int>", &t_trackHitOutMeasTIB);
   tree_->Branch("t_trackHitOutMeasTID", "std::vector<int>", &t_trackHitOutMeasTID);
   tree_->Branch("t_trackOutPosOutHitDr", "std::vector<double>", &t_trackOutPosOutHitDr);
   tree_->Branch("t_trackL", "std::vector<double>", &t_trackL);
 
   tree_->Branch("t_trackDxy", "std::vector<double>", &t_trackDxy);
   tree_->Branch("t_trackDxyBS", "std::vector<double>", &t_trackDxyBS);
   tree_->Branch("t_trackDz", "std::vector<double>", &t_trackDz);
   tree_->Branch("t_trackDzBS", "std::vector<double>", &t_trackDzBS);
   tree_->Branch("t_trackDxyPV", "std::vector<double>", &t_trackDxyPV);
   tree_->Branch("t_trackDzPV", "std::vector<double>", &t_trackDzPV);
   tree_->Branch("t_trackChiSq", "std::vector<double>", &t_trackChiSq);
   tree_->Branch("t_trackPVIdx", "std::vector<int>", &t_trackPVIdx);
 
   t_maxNearP31x31 = new std::vector<double>();
   t_maxNearP21x21 = new std::vector<double>();
 
   tree_->Branch("t_maxNearP31x31", "std::vector<double>", &t_maxNearP31x31);
   tree_->Branch("t_maxNearP21x21", "std::vector<double>", &t_maxNearP21x21);
 
   t_ecalSpike11x11 = new std::vector<int>();
   t_e7x7 = new std::vector<double>();
   t_e9x9 = new std::vector<double>();
   t_e11x11 = new std::vector<double>();
   t_e15x15 = new std::vector<double>();
 
   tree_->Branch("t_ecalSpike11x11", "std::vector<int>", &t_ecalSpike11x11);
   tree_->Branch("t_e7x7", "std::vector<double>", &t_e7x7);
   tree_->Branch("t_e9x9", "std::vector<double>", &t_e9x9);
   tree_->Branch("t_e11x11", "std::vector<double>", &t_e11x11);
   tree_->Branch("t_e15x15", "std::vector<double>", &t_e15x15);
 
   t_e7x7_10Sig = new std::vector<double>();
   t_e9x9_10Sig = new std::vector<double>();
   t_e11x11_10Sig = new std::vector<double>();
   t_e15x15_10Sig = new std::vector<double>();
   t_e7x7_15Sig = new std::vector<double>();
   t_e9x9_15Sig = new std::vector<double>();
   t_e11x11_15Sig = new std::vector<double>();
   t_e15x15_15Sig = new std::vector<double>();
   t_e7x7_20Sig = new std::vector<double>();
   t_e9x9_20Sig = new std::vector<double>();
   t_e11x11_20Sig = new std::vector<double>();
   t_e15x15_20Sig = new std::vector<double>();
   t_e7x7_25Sig = new std::vector<double>();
   t_e9x9_25Sig = new std::vector<double>();
   t_e11x11_25Sig = new std::vector<double>();
   t_e15x15_25Sig = new std::vector<double>();
   t_e7x7_30Sig = new std::vector<double>();
   t_e9x9_30Sig = new std::vector<double>();
   t_e11x11_30Sig = new std::vector<double>();
   t_e15x15_30Sig = new std::vector<double>();
 
   tree_->Branch("t_e7x7_10Sig", "std::vector<double>", &t_e7x7_10Sig);
   tree_->Branch("t_e9x9_10Sig", "std::vector<double>", &t_e9x9_10Sig);
   tree_->Branch("t_e11x11_10Sig", "std::vector<double>", &t_e11x11_10Sig);
   tree_->Branch("t_e15x15_10Sig", "std::vector<double>", &t_e15x15_10Sig);
   tree_->Branch("t_e7x7_15Sig", "std::vector<double>", &t_e7x7_15Sig);
   tree_->Branch("t_e9x9_15Sig", "std::vector<double>", &t_e9x9_15Sig);
   tree_->Branch("t_e11x11_15Sig", "std::vector<double>", &t_e11x11_15Sig);
   tree_->Branch("t_e15x15_15Sig", "std::vector<double>", &t_e15x15_15Sig);
   tree_->Branch("t_e7x7_20Sig", "std::vector<double>", &t_e7x7_20Sig);
   tree_->Branch("t_e9x9_20Sig", "std::vector<double>", &t_e9x9_20Sig);
   tree_->Branch("t_e11x11_20Sig", "std::vector<double>", &t_e11x11_20Sig);
   tree_->Branch("t_e15x15_20Sig", "std::vector<double>", &t_e15x15_20Sig);
   tree_->Branch("t_e7x7_25Sig", "std::vector<double>", &t_e7x7_25Sig);
   tree_->Branch("t_e9x9_25Sig", "std::vector<double>", &t_e9x9_25Sig);
   tree_->Branch("t_e11x11_25Sig", "std::vector<double>", &t_e11x11_25Sig);
   tree_->Branch("t_e15x15_25Sig", "std::vector<double>", &t_e15x15_25Sig);
   tree_->Branch("t_e7x7_30Sig", "std::vector<double>", &t_e7x7_30Sig);
   tree_->Branch("t_e9x9_30Sig", "std::vector<double>", &t_e9x9_30Sig);
   tree_->Branch("t_e11x11_30Sig", "std::vector<double>", &t_e11x11_30Sig);
   tree_->Branch("t_e15x15_30Sig", "std::vector<double>", &t_e15x15_30Sig);
 
   if (doMC_) {
     t_esim7x7 = new std::vector<double>();
     t_esim9x9 = new std::vector<double>();
     t_esim11x11 = new std::vector<double>();
     t_esim15x15 = new std::vector<double>();
 
     t_esim7x7Matched = new std::vector<double>();
     t_esim9x9Matched = new std::vector<double>();
     t_esim11x11Matched = new std::vector<double>();
     t_esim15x15Matched = new std::vector<double>();
 
     t_esim7x7Rest = new std::vector<double>();
     t_esim9x9Rest = new std::vector<double>();
     t_esim11x11Rest = new std::vector<double>();
     t_esim15x15Rest = new std::vector<double>();
 
     t_esim7x7Photon = new std::vector<double>();
     t_esim9x9Photon = new std::vector<double>();
     t_esim11x11Photon = new std::vector<double>();
     t_esim15x15Photon = new std::vector<double>();
 
     t_esim7x7NeutHad = new std::vector<double>();
     t_esim9x9NeutHad = new std::vector<double>();
     t_esim11x11NeutHad = new std::vector<double>();
     t_esim15x15NeutHad = new std::vector<double>();
 
     t_esim7x7CharHad = new std::vector<double>();
     t_esim9x9CharHad = new std::vector<double>();
     t_esim11x11CharHad = new std::vector<double>();
     t_esim15x15CharHad = new std::vector<double>();
 
     t_trkEcalEne = new std::vector<double>();
     t_simTrackP = new std::vector<double>();
     t_esimPdgId = new std::vector<double>();
 
     tree_->Branch("t_esim7x7", "std::vector<double>", &t_esim7x7);
     tree_->Branch("t_esim9x9", "std::vector<double>", &t_esim9x9);
     tree_->Branch("t_esim11x11", "std::vector<double>", &t_esim11x11);
     tree_->Branch("t_esim15x15", "std::vector<double>", &t_esim15x15);
 
     tree_->Branch("t_esim7x7Matched", "std::vector<double>", &t_esim7x7Matched);
     tree_->Branch("t_esim9x9Matched", "std::vector<double>", &t_esim9x9Matched);
     tree_->Branch("t_esim11x11Matched", "std::vector<double>", &t_esim11x11Matched);
     tree_->Branch("t_esim15x15Matched", "std::vector<double>", &t_esim15x15Matched);
 
     tree_->Branch("t_esim7x7Rest", "std::vector<double>", &t_esim7x7Rest);
     tree_->Branch("t_esim9x9Rest", "std::vector<double>", &t_esim9x9Rest);
     tree_->Branch("t_esim11x11Rest", "std::vector<double>", &t_esim11x11Rest);
     tree_->Branch("t_esim15x15Rest", "std::vector<double>", &t_esim15x15Rest);
 
     tree_->Branch("t_esim7x7Photon", "std::vector<double>", &t_esim7x7Photon);
     tree_->Branch("t_esim9x9Photon", "std::vector<double>", &t_esim9x9Photon);
     tree_->Branch("t_esim11x11Photon", "std::vector<double>", &t_esim11x11Photon);
     tree_->Branch("t_esim15x15Photon", "std::vector<double>", &t_esim15x15Photon);
 
     tree_->Branch("t_esim7x7NeutHad", "std::vector<double>", &t_esim7x7NeutHad);
     tree_->Branch("t_esim9x9NeutHad", "std::vector<double>", &t_esim9x9NeutHad);
     tree_->Branch("t_esim11x11NeutHad", "std::vector<double>", &t_esim11x11NeutHad);
     tree_->Branch("t_esim15x15NeutHad", "std::vector<double>", &t_esim15x15NeutHad);
 
     tree_->Branch("t_esim7x7CharHad", "std::vector<double>", &t_esim7x7CharHad);
     tree_->Branch("t_esim9x9CharHad", "std::vector<double>", &t_esim9x9CharHad);
     tree_->Branch("t_esim11x11CharHad", "std::vector<double>", &t_esim11x11CharHad);
     tree_->Branch("t_esim15x15CharHad", "std::vector<double>", &t_esim15x15CharHad);
 
     tree_->Branch("t_trkEcalEne", "std::vector<double>", &t_trkEcalEne);
     tree_->Branch("t_simTrackP", "std::vector<double>", &t_simTrackP);
     tree_->Branch("t_esimPdgId", "std::vector<double>", &t_esimPdgId);
   }
 
   t_maxNearHcalP3x3 = new std::vector<double>();
   t_maxNearHcalP5x5 = new std::vector<double>();
   t_maxNearHcalP7x7 = new std::vector<double>();
   t_h3x3 = new std::vector<double>();
   t_h5x5 = new std::vector<double>();
   t_h7x7 = new std::vector<double>();
   t_h3x3Sig = new std::vector<double>();
   t_h5x5Sig = new std::vector<double>();
   t_h7x7Sig = new std::vector<double>();
   t_infoHcal = new std::vector<int>();
 
   if (doMC_) {
     t_trkHcalEne = new std::vector<double>();
     t_hsim3x3 = new std::vector<double>();
     t_hsim5x5 = new std::vector<double>();
     t_hsim7x7 = new std::vector<double>();
     t_hsim3x3Matched = new std::vector<double>();
     t_hsim5x5Matched = new std::vector<double>();
     t_hsim7x7Matched = new std::vector<double>();
     t_hsim3x3Rest = new std::vector<double>();
     t_hsim5x5Rest = new std::vector<double>();
     t_hsim7x7Rest = new std::vector<double>();
     t_hsim3x3Photon = new std::vector<double>();
     t_hsim5x5Photon = new std::vector<double>();
     t_hsim7x7Photon = new std::vector<double>();
     t_hsim3x3NeutHad = new std::vector<double>();
     t_hsim5x5NeutHad = new std::vector<double>();
     t_hsim7x7NeutHad = new std::vector<double>();
     t_hsim3x3CharHad = new std::vector<double>();
     t_hsim5x5CharHad = new std::vector<double>();
     t_hsim7x7CharHad = new std::vector<double>();
   }
 
   tree_->Branch("t_maxNearHcalP3x3", "std::vector<double>", &t_maxNearHcalP3x3);
   tree_->Branch("t_maxNearHcalP5x5", "std::vector<double>", &t_maxNearHcalP5x5);
   tree_->Branch("t_maxNearHcalP7x7", "std::vector<double>", &t_maxNearHcalP7x7);
   tree_->Branch("t_h3x3", "std::vector<double>", &t_h3x3);
   tree_->Branch("t_h5x5", "std::vector<double>", &t_h5x5);
   tree_->Branch("t_h7x7", "std::vector<double>", &t_h7x7);
   tree_->Branch("t_h3x3Sig", "std::vector<double>", &t_h3x3Sig);
   tree_->Branch("t_h5x5Sig", "std::vector<double>", &t_h5x5Sig);
   tree_->Branch("t_h7x7Sig", "std::vector<double>", &t_h7x7Sig);
   tree_->Branch("t_infoHcal", "std::vector<int>", &t_infoHcal);
 
   if (doMC_) {
     tree_->Branch("t_trkHcalEne", "std::vector<double>", &t_trkHcalEne);
     tree_->Branch("t_hsim3x3", "std::vector<double>", &t_hsim3x3);
     tree_->Branch("t_hsim5x5", "std::vector<double>", &t_hsim5x5);
     tree_->Branch("t_hsim7x7", "std::vector<double>", &t_hsim7x7);
     tree_->Branch("t_hsim3x3Matched", "std::vector<double>", &t_hsim3x3Matched);
     tree_->Branch("t_hsim5x5Matched", "std::vector<double>", &t_hsim5x5Matched);
     tree_->Branch("t_hsim7x7Matched", "std::vector<double>", &t_hsim7x7Matched);
     tree_->Branch("t_hsim3x3Rest", "std::vector<double>", &t_hsim3x3Rest);
     tree_->Branch("t_hsim5x5Rest", "std::vector<double>", &t_hsim5x5Rest);
     tree_->Branch("t_hsim7x7Rest", "std::vector<double>", &t_hsim7x7Rest);
     tree_->Branch("t_hsim3x3Photon", "std::vector<double>", &t_hsim3x3Photon);
     tree_->Branch("t_hsim5x5Photon", "std::vector<double>", &t_hsim5x5Photon);
     tree_->Branch("t_hsim7x7Photon", "std::vector<double>", &t_hsim7x7Photon);
     tree_->Branch("t_hsim3x3NeutHad", "std::vector<double>", &t_hsim3x3NeutHad);
     tree_->Branch("t_hsim5x5NeutHad", "std::vector<double>", &t_hsim5x5NeutHad);
     tree_->Branch("t_hsim7x7NeutHad", "std::vector<double>", &t_hsim7x7NeutHad);
     tree_->Branch("t_hsim3x3CharHad", "std::vector<double>", &t_hsim3x3CharHad);
     tree_->Branch("t_hsim5x5CharHad", "std::vector<double>", &t_hsim5x5CharHad);
     tree_->Branch("t_hsim7x7CharHad", "std::vector<double>", &t_hsim7x7CharHad);
   }
   tree_->Branch("t_nTracks", &t_nTracks, "t_nTracks/I");
 }

◆ clearTreeVectors()

void IsolatedTracksNxN::clearTreeVectors ( )

private

Definition at line 2057 of file IsolatedTracksNxN.cc.

References doMC_, mps_fire::i, t_e11x11, t_e11x11_10Sig, t_e11x11_15Sig, t_e11x11_20Sig, t_e11x11_25Sig, t_e11x11_30Sig, t_e15x15, t_e15x15_10Sig, t_e15x15_15Sig, t_e15x15_20Sig, t_e15x15_25Sig, t_e15x15_30Sig, t_e7x7, t_e7x7_10Sig, t_e7x7_15Sig, t_e7x7_20Sig, t_e7x7_25Sig, t_e7x7_30Sig, t_e9x9, t_e9x9_10Sig, t_e9x9_15Sig, t_e9x9_20Sig, t_e9x9_25Sig, t_e9x9_30Sig, t_ecalSpike11x11, t_esim11x11, t_esim11x11CharHad, t_esim11x11Matched, t_esim11x11NeutHad, t_esim11x11Photon, t_esim11x11Rest, t_esim15x15, t_esim15x15CharHad, t_esim15x15Matched, t_esim15x15NeutHad, t_esim15x15Photon, t_esim15x15Rest, t_esim7x7, t_esim7x7CharHad, t_esim7x7Matched, t_esim7x7NeutHad, t_esim7x7Photon, t_esim7x7Rest, t_esim9x9, t_esim9x9CharHad, t_esim9x9Matched, t_esim9x9NeutHad, t_esim9x9Photon, t_esim9x9Rest, t_esimPdgId, t_h3x3, t_h3x3Sig, t_h5x5, t_h5x5Sig, t_h7x7, t_h7x7Sig, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NeutHad, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NeutHad, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim7x7, t_hsim7x7CharHad, t_hsim7x7Matched, t_hsim7x7NeutHad, t_hsim7x7Photon, t_hsim7x7Rest, t_infoHcal, t_jetEta, t_jetPhi, t_jetPt, t_L1AlgoNames, t_L1CenJetEta, t_L1CenJetPhi, t_L1CenJetPt, t_L1Decision, t_L1FwdJetEta, t_L1FwdJetPhi, t_L1FwdJetPt, t_L1IsoEMEta, t_L1IsoEMPhi, t_L1IsoEMPt, t_L1METEta, t_L1METPhi, t_L1METPt, t_L1MuonEta, t_L1MuonPhi, t_L1MuonPt, t_L1NonIsoEMEta, t_L1NonIsoEMPhi, t_L1NonIsoEMPt, t_L1PreScale, t_L1TauJetEta, t_L1TauJetPhi, t_L1TauJetPt, t_maxNearHcalP3x3, t_maxNearHcalP5x5, t_maxNearHcalP7x7, t_maxNearP21x21, t_maxNearP31x31, t_NLayersCrossed, t_nTrksJetCalo, t_nTrksJetVtx, t_PVisValid, t_PVndof, t_PVNTracks, t_PVNTracksHP, t_PVNTracksHPWt, t_PVNTracksWt, t_PVTracksSumPt, t_PVTracksSumPtHP, t_PVTracksSumPtHPWt, t_PVTracksSumPtWt, t_PVx, t_PVy, t_PVz, t_simTrackP, t_trackChiSq, t_trackChiSqAll, t_trackDxy, t_trackDxyAll, t_trackDxyBS, t_trackDxyPV, t_trackDxyPVAll, t_trackDz, t_trackDzAll, t_trackDzBS, t_trackDzPV, t_trackDzPVAll, t_trackEcalEta, t_trackEcalPhi, t_trackEta, t_trackEtaAll, t_trackHcalEta, t_trackHcalPhi, t_trackHitInMeasTEC, t_trackHitInMeasTIB, t_trackHitInMeasTID, t_trackHitInMeasTOB, t_trackHitInMissTEC, t_trackHitInMissTIB, t_trackHitInMissTIBTID, t_trackHitInMissTID, t_trackHitInMissTOB, t_trackHitOutMeasTEC, t_trackHitOutMeasTIB, t_trackHitOutMeasTID, t_trackHitOutMeasTOB, t_trackHitOutMissTEC, t_trackHitOutMissTIB, t_trackHitOutMissTID, t_trackHitOutMissTOB, t_trackHitOutMissTOBTEC, t_trackHitsTEC, t_trackHitsTOB, t_trackL, t_trackNOuterHits, t_trackOutPosOutHitDr, t_trackP, t_trackPAll, t_trackPdgIdAll, t_trackPhi, t_trackPhiAll, t_trackPt, t_trackPtAll, t_trackPVIdx, t_trkEcalEne, and t_trkHcalEne.

Referenced by analyze().

                                          {
   t_PVx->clear();
   t_PVy->clear();
   t_PVz->clear();
   t_PVisValid->clear();
   t_PVndof->clear();
   t_PVNTracks->clear();
   t_PVNTracksWt->clear();
   t_PVTracksSumPt->clear();
   t_PVTracksSumPtWt->clear();
   t_PVNTracksHP->clear();
   t_PVNTracksHPWt->clear();
   t_PVTracksSumPtHP->clear();
   t_PVTracksSumPtHPWt->clear();
 
   for (int i = 0; i < 128; i++)
     t_L1Decision[i] = 0;
   t_L1AlgoNames->clear();
   t_L1PreScale->clear();
 
   t_L1CenJetPt->clear();
   t_L1CenJetEta->clear();
   t_L1CenJetPhi->clear();
   t_L1FwdJetPt->clear();
   t_L1FwdJetEta->clear();
   t_L1FwdJetPhi->clear();
   t_L1TauJetPt->clear();
   t_L1TauJetEta->clear();
   t_L1TauJetPhi->clear();
   t_L1MuonPt->clear();
   t_L1MuonEta->clear();
   t_L1MuonPhi->clear();
   t_L1IsoEMPt->clear();
   t_L1IsoEMEta->clear();
   t_L1IsoEMPhi->clear();
   t_L1NonIsoEMPt->clear();
   t_L1NonIsoEMEta->clear();
   t_L1NonIsoEMPhi->clear();
   t_L1METPt->clear();
   t_L1METEta->clear();
   t_L1METPhi->clear();
 
   t_jetPt->clear();
   t_jetEta->clear();
   t_jetPhi->clear();
   t_nTrksJetCalo->clear();
   t_nTrksJetVtx->clear();
 
   t_trackPAll->clear();
   t_trackEtaAll->clear();
   t_trackPhiAll->clear();
   t_trackPdgIdAll->clear();
   t_trackPtAll->clear();
   t_trackDxyAll->clear();
   t_trackDzAll->clear();
   t_trackDxyPVAll->clear();
   t_trackDzPVAll->clear();
   t_trackChiSqAll->clear();
 
   t_trackP->clear();
   t_trackPt->clear();
   t_trackEta->clear();
   t_trackPhi->clear();
   t_trackEcalEta->clear();
   t_trackEcalPhi->clear();
   t_trackHcalEta->clear();
   t_trackHcalPhi->clear();
   t_NLayersCrossed->clear();
   t_trackNOuterHits->clear();
   t_trackDxy->clear();
   t_trackDxyBS->clear();
   t_trackDz->clear();
   t_trackDzBS->clear();
   t_trackDxyPV->clear();
   t_trackDzPV->clear();
   t_trackChiSq->clear();
   t_trackPVIdx->clear();
   t_trackHitsTOB->clear();
   t_trackHitsTEC->clear();
   t_trackHitInMissTOB->clear();
   t_trackHitInMissTEC->clear();
   t_trackHitInMissTIB->clear();
   t_trackHitInMissTID->clear();
   t_trackHitInMissTIBTID->clear();
   t_trackHitOutMissTOB->clear();
   t_trackHitOutMissTEC->clear();
   t_trackHitOutMissTIB->clear();
   t_trackHitOutMissTID->clear();
   t_trackHitOutMissTOBTEC->clear();
   t_trackHitInMeasTOB->clear();
   t_trackHitInMeasTEC->clear();
   t_trackHitInMeasTIB->clear();
   t_trackHitInMeasTID->clear();
   t_trackHitOutMeasTOB->clear();
   t_trackHitOutMeasTEC->clear();
   t_trackHitOutMeasTIB->clear();
   t_trackHitOutMeasTID->clear();
   t_trackOutPosOutHitDr->clear();
   t_trackL->clear();
 
   t_maxNearP31x31->clear();
   t_maxNearP21x21->clear();
 
   t_ecalSpike11x11->clear();
   t_e7x7->clear();
   t_e9x9->clear();
   t_e11x11->clear();
   t_e15x15->clear();
 
   t_e7x7_10Sig->clear();
   t_e9x9_10Sig->clear();
   t_e11x11_10Sig->clear();
   t_e15x15_10Sig->clear();
   t_e7x7_15Sig->clear();
   t_e9x9_15Sig->clear();
   t_e11x11_15Sig->clear();
   t_e15x15_15Sig->clear();
   t_e7x7_20Sig->clear();
   t_e9x9_20Sig->clear();
   t_e11x11_20Sig->clear();
   t_e15x15_20Sig->clear();
   t_e7x7_25Sig->clear();
   t_e9x9_25Sig->clear();
   t_e11x11_25Sig->clear();
   t_e15x15_25Sig->clear();
   t_e7x7_30Sig->clear();
   t_e9x9_30Sig->clear();
   t_e11x11_30Sig->clear();
   t_e15x15_30Sig->clear();
 
   if (doMC_) {
     t_simTrackP->clear();
     t_esimPdgId->clear();
     t_trkEcalEne->clear();
 
     t_esim7x7->clear();
     t_esim9x9->clear();
     t_esim11x11->clear();
     t_esim15x15->clear();
 
     t_esim7x7Matched->clear();
     t_esim9x9Matched->clear();
     t_esim11x11Matched->clear();
     t_esim15x15Matched->clear();
 
     t_esim7x7Rest->clear();
     t_esim9x9Rest->clear();
     t_esim11x11Rest->clear();
     t_esim15x15Rest->clear();
 
     t_esim7x7Photon->clear();
     t_esim9x9Photon->clear();
     t_esim11x11Photon->clear();
     t_esim15x15Photon->clear();
 
     t_esim7x7NeutHad->clear();
     t_esim9x9NeutHad->clear();
     t_esim11x11NeutHad->clear();
     t_esim15x15NeutHad->clear();
 
     t_esim7x7CharHad->clear();
     t_esim9x9CharHad->clear();
     t_esim11x11CharHad->clear();
     t_esim15x15CharHad->clear();
   }
 
   t_maxNearHcalP3x3->clear();
   t_maxNearHcalP5x5->clear();
   t_maxNearHcalP7x7->clear();
 
   t_h3x3->clear();
   t_h5x5->clear();
   t_h7x7->clear();
   t_h3x3Sig->clear();
   t_h5x5Sig->clear();
   t_h7x7Sig->clear();
 
   t_infoHcal->clear();
 
   if (doMC_) {
     t_trkHcalEne->clear();
 
     t_hsim3x3->clear();
     t_hsim5x5->clear();
     t_hsim7x7->clear();
     t_hsim3x3Matched->clear();
     t_hsim5x5Matched->clear();
     t_hsim7x7Matched->clear();
     t_hsim3x3Rest->clear();
     t_hsim5x5Rest->clear();
     t_hsim7x7Rest->clear();
     t_hsim3x3Photon->clear();
     t_hsim5x5Photon->clear();
     t_hsim7x7Photon->clear();
     t_hsim3x3NeutHad->clear();
     t_hsim5x5NeutHad->clear();
     t_hsim7x7NeutHad->clear();
     t_hsim3x3CharHad->clear();
     t_hsim5x5CharHad->clear();
     t_hsim7x7CharHad->clear();
   }
 }

◆ endJob()

void IsolatedTracksNxN::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 2041 of file IsolatedTracksNxN.cc.

References first, h_L1AlgoNames, l1AlgoMap_, L1TriggerAlgoInfo_, Skims_PA_cff::name, and nEventProc_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                {
   if (L1TriggerAlgoInfo_) {
     std::map<std::pair<unsigned int, std::string>, int>::iterator itr;
     for (itr = l1AlgoMap_.begin(); itr != l1AlgoMap_.end(); itr++) {
       edm::LogVerbatim("IsoTrack") << " ****endjob**** " << (itr->first).first << " " << (itr->first).second << " "
                                    << itr->second;
       int ibin = (itr->first).first;
       TString name((itr->first).second);
       h_L1AlgoNames->GetXaxis()->SetBinLabel(ibin + 1, name);
     }
     edm::LogVerbatim("IsoTrack") << "Number of Events Processed " << nEventProc_;
   }
 }

◆ fillDescriptions()

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

static

Definition at line 735 of file IsolatedTracksNxN.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_2022v14_cff::InputTag.

                                                                                    {
   edm::ParameterSetDescription desc;
   desc.addUntracked<bool>("doMC", false);
   desc.addUntracked<bool>("writeAllTracks", false);
   desc.addUntracked<int>("verbosity", 1);
   desc.addUntracked<double>("pvTracksPtMin", 0.200);
   desc.addUntracked<int>("debugTracks", 0);
   desc.addUntracked<bool>("printTrkHitPattern", true);
   desc.addUntracked<double>("minTrackP", 1.0);
   desc.addUntracked<double>("maxTrackEta", 2.6);
   desc.addUntracked<bool>("debugL1Info", false);
   desc.addUntracked<bool>("l1TriggerAlgoInfo", false);
   desc.add<edm::InputTag>("l1extraTauJetSource", edm::InputTag("l1extraParticles", "Tau"));
   desc.add<edm::InputTag>("l1extraCenJetSource", edm::InputTag("l1extraParticles", "Central"));
   desc.add<edm::InputTag>("l1extraFwdJetSource", edm::InputTag("l1extraParticles", "Forward"));
   desc.add<edm::InputTag>("l1extraMuonSource", edm::InputTag("l1extraParticles"));
   desc.add<edm::InputTag>("l1extraIsoEmSource", edm::InputTag("l1extraParticles", "Isolated"));
   desc.add<edm::InputTag>("l1extraNonIsoEmSource", edm::InputTag("l1extraParticles", "NonIsolated"));
   desc.add<edm::InputTag>("l1GTReadoutRcdSource", edm::InputTag("gtDigis"));
   desc.add<edm::InputTag>("l1GTObjectMapRcdSource", edm::InputTag("hltL1GtObjectMap"));
   desc.add<edm::InputTag>("jetSource", edm::InputTag("iterativeCone5CaloJets"));
   desc.add<edm::InputTag>("jetExtender", edm::InputTag("iterativeCone5JetExtender"));
   desc.add<edm::InputTag>("hbheRecHitSource", edm::InputTag("hbhereco"));
   desc.addUntracked<double>("maxNearTrackPT", 1.0);
   desc.addUntracked<double>("timeMinCutECAL", -500.0);
   desc.addUntracked<double>("timeMaxCutECAL", 500.0);
   desc.addUntracked<double>("timeMinCutHCAL", -500.0);
   desc.addUntracked<double>("timeMaxCutHCAL", 500.0);
   descriptions.add("isolatedTracksNxN", desc);
 }

◆ printTrack()

void IsolatedTracksNxN::printTrack ( const reco::Track * pTrack )

private

Definition at line 2747 of file IsolatedTracksNxN.cc.

Referenced by ntuplePrintersDiff.TrackPrinter::diff(), and ntuplePrintersDiff.TrackPrinter::printTrackAndMatchedTrackingParticles().

                                                           {
   std::string theTrackQuality = "highPurity";
   reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
 
   edm::LogVerbatim("IsoTrack") << " Reference Point " << pTrack->referencePoint() << "\n TrackMmentum "
                                << pTrack->momentum() << " (pt,eta,phi)(" << pTrack->pt() << "," << pTrack->eta() << ","
                                << pTrack->phi() << ")"
                                << " p " << pTrack->p() << "\n Normalized chi2 " << pTrack->normalizedChi2()
                                << "  charge " << pTrack->charge() << " qoverp() " << pTrack->qoverp() << "+-"
                                << pTrack->qoverpError() << " d0 " << pTrack->d0() << "\n NValidHits "
                                << pTrack->numberOfValidHits() << "  NLostHits " << pTrack->numberOfLostHits()
                                << " TrackQuality " << pTrack->qualityName(trackQuality_) << " "
                                << pTrack->quality(trackQuality_);
 
   if (printTrkHitPattern_) {
     const reco::HitPattern &p = pTrack->hitPattern();
 
     std::ostringstream st1;
     st1 << "default ";
     for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
       p.printHitPattern(reco::HitPattern::TRACK_HITS, i, st1);
     }
     edm::LogVerbatim("IsoTrack") << st1.str();
     std::ostringstream st2;
     st2 << "trackerMissingInnerHits() ";
     for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS); i++) {
       p.printHitPattern(reco::HitPattern::MISSING_INNER_HITS, i, st2);
     }
     edm::LogVerbatim("IsoTrack") << st2.str();
     std::ostringstream st3;
     st3 << "trackerMissingOuterHits() ";
     for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::MISSING_OUTER_HITS); i++) {
       p.printHitPattern(reco::HitPattern::MISSING_OUTER_HITS, i, st3);
     }
     edm::LogVerbatim("IsoTrack") << st3.str();
 
     edm::LogVerbatim("IsoTrack") << "\n \t trackerLayersWithMeasurement() " << p.trackerLayersWithMeasurement()
                                  << "\n \t pixelLayersWithMeasurement() " << p.pixelLayersWithMeasurement()
                                  << "\n \t stripLayersWithMeasurement() " << p.stripLayersWithMeasurement()
                                  << "\n \t pixelBarrelLayersWithMeasurement() " << p.pixelBarrelLayersWithMeasurement()
                                  << "\n \t pixelEndcapLayersWithMeasurement() " << p.pixelEndcapLayersWithMeasurement()
                                  << "\n \t stripTIBLayersWithMeasurement() " << p.stripTIBLayersWithMeasurement()
                                  << "\n \t stripTIDLayersWithMeasurement() " << p.stripTIDLayersWithMeasurement()
                                  << "\n \t stripTOBLayersWithMeasurement() " << p.stripTOBLayersWithMeasurement()
                                  << "\n \t stripTECLayersWithMeasurement() " << p.stripTECLayersWithMeasurement();
   }
 }

Member Data Documentation

◆ debugL1Info_

const bool IsolatedTracksNxN::debugL1Info_

private

Definition at line 175 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ debugTrks_

const int IsolatedTracksNxN::debugTrks_

private

Definition at line 172 of file IsolatedTracksNxN.cc.

Referenced by analyze().

◆ doMC_

const bool IsolatedTracksNxN::doMC_

private

Definition at line 169 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and IsolatedTracksNxN().

◆ genPartEtaBins

double IsolatedTracksNxN::genPartEtaBins[NEtaBins+1]

private

Definition at line 212 of file IsolatedTracksNxN.cc.

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

◆ genPartPBins

double IsolatedTracksNxN::genPartPBins[NPBins+1]

private

Definition at line 212 of file IsolatedTracksNxN.cc.

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

◆ h_L1AlgoNames

TH1I* IsolatedTracksNxN::h_L1AlgoNames

private

Definition at line 216 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), and endJob().

◆ h_maxNearP15x15

TH1F* IsolatedTracksNxN::h_maxNearP15x15[NPBins][NEtaBins]

private

Definition at line 214 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_maxNearP21x21

TH1F * IsolatedTracksNxN::h_maxNearP21x21[NPBins][NEtaBins]

private

Definition at line 214 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_maxNearP25x25

TH1F * IsolatedTracksNxN::h_maxNearP25x25[NPBins][NEtaBins]

private

Definition at line 214 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_maxNearP31x31

TH1F * IsolatedTracksNxN::h_maxNearP31x31[NPBins][NEtaBins]

private

Definition at line 214 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_nTracks

TH1F* IsolatedTracksNxN::h_nTracks

private

Definition at line 218 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_PVTracksWt

TH1F* IsolatedTracksNxN::h_PVTracksWt

private

Definition at line 217 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEta_0

TH1F * IsolatedTracksNxN::h_recEta_0

private

Definition at line 219 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEta_1

TH1F * IsolatedTracksNxN::h_recEta_1

private

Definition at line 221 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEta_2

TH1F * IsolatedTracksNxN::h_recEta_2

private

Definition at line 223 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEtaP_0

TH2F * IsolatedTracksNxN::h_recEtaP_0

private

Definition at line 220 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEtaP_1

TH2F * IsolatedTracksNxN::h_recEtaP_1

private

Definition at line 222 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEtaP_2

TH2F * IsolatedTracksNxN::h_recEtaP_2

private

Definition at line 224 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEtaPt_0

TH2F* IsolatedTracksNxN::h_recEtaPt_0

private

Definition at line 220 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEtaPt_1

TH2F* IsolatedTracksNxN::h_recEtaPt_1

private

Definition at line 222 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recEtaPt_2

TH2F* IsolatedTracksNxN::h_recEtaPt_2

private

Definition at line 224 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recP_0

TH1F * IsolatedTracksNxN::h_recP_0

private

Definition at line 219 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recP_1

TH1F * IsolatedTracksNxN::h_recP_1

private

Definition at line 221 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recP_2

TH1F * IsolatedTracksNxN::h_recP_2

private

Definition at line 223 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recPhi_0

TH1F * IsolatedTracksNxN::h_recPhi_0

private

Definition at line 219 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recPhi_1

TH1F * IsolatedTracksNxN::h_recPhi_1

private

Definition at line 221 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recPhi_2

TH1F * IsolatedTracksNxN::h_recPhi_2

private

Definition at line 223 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recPt_0

TH1F* IsolatedTracksNxN::h_recPt_0

private

Definition at line 219 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recPt_1

TH1F* IsolatedTracksNxN::h_recPt_1

private

Definition at line 221 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ h_recPt_2

TH1F* IsolatedTracksNxN::h_recPt_2

private

Definition at line 223 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ initL1_

bool IsolatedTracksNxN::initL1_

private

Definition at line 177 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and beginJob().

◆ l1AlgoMap_

std::map<std::pair<unsigned int, std::string>, int> IsolatedTracksNxN::l1AlgoMap_

private

Definition at line 166 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and endJob().

◆ L1TriggerAlgoInfo_

const bool IsolatedTracksNxN::L1TriggerAlgoInfo_

private

Definition at line 175 of file IsolatedTracksNxN.cc.

Referenced by analyze(), endJob(), and IsolatedTracksNxN().

◆ m_l1GtUtils

std::unique_ptr<L1GtUtils> IsolatedTracksNxN::m_l1GtUtils

private

Definition at line 161 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ m_triggerMaskAlgoTrig

std::vector<unsigned int> IsolatedTracksNxN::m_triggerMaskAlgoTrig

private

Definition at line 167 of file IsolatedTracksNxN.cc.

◆ maxTrackEta_

const double IsolatedTracksNxN::maxTrackEta_

private

Definition at line 174 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ minTrackP_

const double IsolatedTracksNxN::minTrackP_

private

Definition at line 174 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ myverbose_

const int IsolatedTracksNxN::myverbose_

private

Definition at line 170 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ nbad_

int IsolatedTracksNxN::nbad_

private

Definition at line 178 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and beginJob().

◆ NEtaBins

constexpr size_t IsolatedTracksNxN::NEtaBins = 3

staticprivate

Definition at line 211 of file IsolatedTracksNxN.cc.

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

◆ nEventProc_

int IsolatedTracksNxN::nEventProc_

private

Definition at line 178 of file IsolatedTracksNxN.cc.

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

◆ nL1BitsMax

constexpr size_t IsolatedTracksNxN::nL1BitsMax = 128

staticprivate

Definition at line 162 of file IsolatedTracksNxN.cc.

◆ NPBins

constexpr size_t IsolatedTracksNxN::NPBins = 15

staticprivate

Definition at line 210 of file IsolatedTracksNxN.cc.

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

◆ printTrkHitPattern_

const bool IsolatedTracksNxN::printTrkHitPattern_

private

Definition at line 173 of file IsolatedTracksNxN.cc.

Referenced by printTrack().

◆ pvTracksPtMin_

const double IsolatedTracksNxN::pvTracksPtMin_

private

Definition at line 171 of file IsolatedTracksNxN.cc.

Referenced by analyze().

◆ t_Bunch

int IsolatedTracksNxN::t_Bunch

private

Definition at line 229 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ t_e11x11

std::vector<double> * IsolatedTracksNxN::t_e11x11

private

Definition at line 275 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e11x11_10Sig

std::vector<double> * IsolatedTracksNxN::t_e11x11_10Sig

private

Definition at line 276 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e11x11_15Sig

std::vector<double> * IsolatedTracksNxN::t_e11x11_15Sig

private

Definition at line 277 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e11x11_20Sig

std::vector<double> * IsolatedTracksNxN::t_e11x11_20Sig

private

Definition at line 278 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e11x11_25Sig

std::vector<double> * IsolatedTracksNxN::t_e11x11_25Sig

private

Definition at line 279 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e11x11_30Sig

std::vector<double> * IsolatedTracksNxN::t_e11x11_30Sig

private

Definition at line 280 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e15x15

std::vector<double> * IsolatedTracksNxN::t_e15x15

private

Definition at line 275 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e15x15_10Sig

std::vector<double> * IsolatedTracksNxN::t_e15x15_10Sig

private

Definition at line 276 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e15x15_15Sig

std::vector<double> * IsolatedTracksNxN::t_e15x15_15Sig

private

Definition at line 277 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e15x15_20Sig

std::vector<double> * IsolatedTracksNxN::t_e15x15_20Sig

private

Definition at line 278 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e15x15_25Sig

std::vector<double> * IsolatedTracksNxN::t_e15x15_25Sig

private

Definition at line 279 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e15x15_30Sig

std::vector<double> * IsolatedTracksNxN::t_e15x15_30Sig

private

Definition at line 280 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e7x7

std::vector<double>* IsolatedTracksNxN::t_e7x7

private

Definition at line 275 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e7x7_10Sig

std::vector<double>* IsolatedTracksNxN::t_e7x7_10Sig

private

Definition at line 276 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e7x7_15Sig

std::vector<double>* IsolatedTracksNxN::t_e7x7_15Sig

private

Definition at line 277 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e7x7_20Sig

std::vector<double>* IsolatedTracksNxN::t_e7x7_20Sig

private

Definition at line 278 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e7x7_25Sig

std::vector<double>* IsolatedTracksNxN::t_e7x7_25Sig

private

Definition at line 279 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e7x7_30Sig

std::vector<double>* IsolatedTracksNxN::t_e7x7_30Sig

private

Definition at line 280 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e9x9

std::vector<double> * IsolatedTracksNxN::t_e9x9

private

Definition at line 275 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e9x9_10Sig

std::vector<double> * IsolatedTracksNxN::t_e9x9_10Sig

private

Definition at line 276 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e9x9_15Sig

std::vector<double> * IsolatedTracksNxN::t_e9x9_15Sig

private

Definition at line 277 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e9x9_20Sig

std::vector<double> * IsolatedTracksNxN::t_e9x9_20Sig

private

Definition at line 278 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e9x9_25Sig

std::vector<double> * IsolatedTracksNxN::t_e9x9_25Sig

private

Definition at line 279 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_e9x9_30Sig

std::vector<double> * IsolatedTracksNxN::t_e9x9_30Sig

private

Definition at line 280 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_ecalSpike11x11

std::vector<int>* IsolatedTracksNxN::t_ecalSpike11x11

private

Definition at line 273 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim11x11

std::vector<double> * IsolatedTracksNxN::t_esim11x11

private

Definition at line 283 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim11x11CharHad

std::vector<double> * IsolatedTracksNxN::t_esim11x11CharHad

private

Definition at line 288 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim11x11Matched

std::vector<double> * IsolatedTracksNxN::t_esim11x11Matched

private

Definition at line 284 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim11x11NeutHad

std::vector<double> * IsolatedTracksNxN::t_esim11x11NeutHad

private

Definition at line 287 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim11x11Photon

std::vector<double> * IsolatedTracksNxN::t_esim11x11Photon

private

Definition at line 286 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim11x11Rest

std::vector<double> * IsolatedTracksNxN::t_esim11x11Rest

private

Definition at line 285 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim15x15

std::vector<double> * IsolatedTracksNxN::t_esim15x15

private

Definition at line 283 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim15x15CharHad

std::vector<double> * IsolatedTracksNxN::t_esim15x15CharHad

private

Definition at line 288 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim15x15Matched

std::vector<double> * IsolatedTracksNxN::t_esim15x15Matched

private

Definition at line 284 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim15x15NeutHad

std::vector<double> * IsolatedTracksNxN::t_esim15x15NeutHad

private

Definition at line 287 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim15x15Photon

std::vector<double> * IsolatedTracksNxN::t_esim15x15Photon

private

Definition at line 286 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim15x15Rest

std::vector<double> * IsolatedTracksNxN::t_esim15x15Rest

private

Definition at line 285 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim7x7

std::vector<double>* IsolatedTracksNxN::t_esim7x7

private

Definition at line 283 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim7x7CharHad

std::vector<double>* IsolatedTracksNxN::t_esim7x7CharHad

private

Definition at line 288 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim7x7Matched

std::vector<double>* IsolatedTracksNxN::t_esim7x7Matched

private

Definition at line 284 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim7x7NeutHad

std::vector<double>* IsolatedTracksNxN::t_esim7x7NeutHad

private

Definition at line 287 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim7x7Photon

std::vector<double>* IsolatedTracksNxN::t_esim7x7Photon

private

Definition at line 286 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim7x7Rest

std::vector<double>* IsolatedTracksNxN::t_esim7x7Rest

private

Definition at line 285 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim9x9

std::vector<double> * IsolatedTracksNxN::t_esim9x9

private

Definition at line 283 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim9x9CharHad

std::vector<double> * IsolatedTracksNxN::t_esim9x9CharHad

private

Definition at line 288 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim9x9Matched

std::vector<double> * IsolatedTracksNxN::t_esim9x9Matched

private

Definition at line 284 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim9x9NeutHad

std::vector<double> * IsolatedTracksNxN::t_esim9x9NeutHad

private

Definition at line 287 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim9x9Photon

std::vector<double> * IsolatedTracksNxN::t_esim9x9Photon

private

Definition at line 286 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esim9x9Rest

std::vector<double> * IsolatedTracksNxN::t_esim9x9Rest

private

Definition at line 285 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_esimPdgId

std::vector<double>* IsolatedTracksNxN::t_esimPdgId

private

Definition at line 282 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_EvtNo

int IsolatedTracksNxN::t_EvtNo

private

Definition at line 229 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ t_h3x3

std::vector<double>* IsolatedTracksNxN::t_h3x3

private

Definition at line 291 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_h3x3Sig

std::vector<double>* IsolatedTracksNxN::t_h3x3Sig

private

Definition at line 292 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_h5x5

std::vector<double> * IsolatedTracksNxN::t_h5x5

private

Definition at line 291 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_h5x5Sig

std::vector<double> * IsolatedTracksNxN::t_h5x5Sig

private

Definition at line 292 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_h7x7

std::vector<double> * IsolatedTracksNxN::t_h7x7

private

Definition at line 291 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_h7x7Sig

std::vector<double> * IsolatedTracksNxN::t_h7x7Sig

private

Definition at line 292 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim3x3

std::vector<double>* IsolatedTracksNxN::t_hsim3x3

private

Definition at line 296 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim3x3CharHad

std::vector<double>* IsolatedTracksNxN::t_hsim3x3CharHad

private

Definition at line 301 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim3x3Matched

std::vector<double>* IsolatedTracksNxN::t_hsim3x3Matched

private

Definition at line 297 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim3x3NeutHad

std::vector<double>* IsolatedTracksNxN::t_hsim3x3NeutHad

private

Definition at line 300 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim3x3Photon

std::vector<double>* IsolatedTracksNxN::t_hsim3x3Photon

private

Definition at line 299 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim3x3Rest

std::vector<double>* IsolatedTracksNxN::t_hsim3x3Rest

private

Definition at line 298 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim5x5

std::vector<double> * IsolatedTracksNxN::t_hsim5x5

private

Definition at line 296 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim5x5CharHad

std::vector<double> * IsolatedTracksNxN::t_hsim5x5CharHad

private

Definition at line 301 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim5x5Matched

std::vector<double> * IsolatedTracksNxN::t_hsim5x5Matched

private

Definition at line 297 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim5x5NeutHad

std::vector<double> * IsolatedTracksNxN::t_hsim5x5NeutHad

private

Definition at line 300 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim5x5Photon

std::vector<double> * IsolatedTracksNxN::t_hsim5x5Photon

private

Definition at line 299 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim5x5Rest

std::vector<double> * IsolatedTracksNxN::t_hsim5x5Rest

private

Definition at line 298 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim7x7

std::vector<double> * IsolatedTracksNxN::t_hsim7x7

private

Definition at line 296 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim7x7CharHad

std::vector<double> * IsolatedTracksNxN::t_hsim7x7CharHad

private

Definition at line 301 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim7x7Matched

std::vector<double> * IsolatedTracksNxN::t_hsim7x7Matched

private

Definition at line 297 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim7x7NeutHad

std::vector<double> * IsolatedTracksNxN::t_hsim7x7NeutHad

private

Definition at line 300 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim7x7Photon

std::vector<double> * IsolatedTracksNxN::t_hsim7x7Photon

private

Definition at line 299 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_hsim7x7Rest

std::vector<double> * IsolatedTracksNxN::t_hsim7x7Rest

private

Definition at line 298 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_infoHcal

std::vector<int>* IsolatedTracksNxN::t_infoHcal

private

Definition at line 293 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_jetEta

std::vector<double> * IsolatedTracksNxN::t_jetEta

private

Definition at line 247 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_jetPhi

std::vector<double> * IsolatedTracksNxN::t_jetPhi

private

Definition at line 247 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_jetPt

std::vector<double>* IsolatedTracksNxN::t_jetPt

private

Definition at line 247 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1AlgoNames

std::vector<std::string>* IsolatedTracksNxN::t_L1AlgoNames

private

Definition at line 230 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1CenJetEta

std::vector<double> * IsolatedTracksNxN::t_L1CenJetEta

private

Definition at line 234 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1CenJetPhi

std::vector<double> * IsolatedTracksNxN::t_L1CenJetPhi

private

Definition at line 234 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1CenJetPt

std::vector<double>* IsolatedTracksNxN::t_L1CenJetPt

private

Definition at line 234 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1Decision

int IsolatedTracksNxN::t_L1Decision[128]

private

Definition at line 232 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), and clearTreeVectors().

◆ t_L1FwdJetEta

std::vector<double> * IsolatedTracksNxN::t_L1FwdJetEta

private

Definition at line 235 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1FwdJetPhi

std::vector<double> * IsolatedTracksNxN::t_L1FwdJetPhi

private

Definition at line 235 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1FwdJetPt

std::vector<double>* IsolatedTracksNxN::t_L1FwdJetPt

private

Definition at line 235 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1IsoEMEta

std::vector<double> * IsolatedTracksNxN::t_L1IsoEMEta

private

Definition at line 238 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1IsoEMPhi

std::vector<double> * IsolatedTracksNxN::t_L1IsoEMPhi

private

Definition at line 238 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1IsoEMPt

std::vector<double>* IsolatedTracksNxN::t_L1IsoEMPt

private

Definition at line 238 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1METEta

std::vector<double> * IsolatedTracksNxN::t_L1METEta

private

Definition at line 240 of file IsolatedTracksNxN.cc.

Referenced by bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1METPhi

std::vector<double> * IsolatedTracksNxN::t_L1METPhi

private

Definition at line 240 of file IsolatedTracksNxN.cc.

Referenced by bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1METPt

std::vector<double>* IsolatedTracksNxN::t_L1METPt

private

Definition at line 240 of file IsolatedTracksNxN.cc.

Referenced by bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1MuonEta

std::vector<double> * IsolatedTracksNxN::t_L1MuonEta

private

Definition at line 237 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1MuonPhi

std::vector<double> * IsolatedTracksNxN::t_L1MuonPhi

private

Definition at line 237 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1MuonPt

std::vector<double>* IsolatedTracksNxN::t_L1MuonPt

private

Definition at line 237 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1NonIsoEMEta

std::vector<double> * IsolatedTracksNxN::t_L1NonIsoEMEta

private

Definition at line 239 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1NonIsoEMPhi

std::vector<double> * IsolatedTracksNxN::t_L1NonIsoEMPhi

private

Definition at line 239 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1NonIsoEMPt

std::vector<double>* IsolatedTracksNxN::t_L1NonIsoEMPt

private

Definition at line 239 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1PreScale

std::vector<int>* IsolatedTracksNxN::t_L1PreScale

private

Definition at line 231 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1TauJetEta

std::vector<double> * IsolatedTracksNxN::t_L1TauJetEta

private

Definition at line 236 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1TauJetPhi

std::vector<double> * IsolatedTracksNxN::t_L1TauJetPhi

private

Definition at line 236 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_L1TauJetPt

std::vector<double>* IsolatedTracksNxN::t_L1TauJetPt

private

Definition at line 236 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_Lumi

int IsolatedTracksNxN::t_Lumi

private

Definition at line 229 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ t_maxNearHcalP3x3

std::vector<double>* IsolatedTracksNxN::t_maxNearHcalP3x3

private

Definition at line 290 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_maxNearHcalP5x5

std::vector<double> * IsolatedTracksNxN::t_maxNearHcalP5x5

private

Definition at line 290 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_maxNearHcalP7x7

std::vector<double> * IsolatedTracksNxN::t_maxNearHcalP7x7

private

Definition at line 290 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_maxNearP21x21

std::vector<double>* IsolatedTracksNxN::t_maxNearP21x21

private

Definition at line 272 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_maxNearP31x31

std::vector<double>* IsolatedTracksNxN::t_maxNearP31x31

private

Definition at line 271 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_NLayersCrossed

std::vector<int>* IsolatedTracksNxN::t_NLayersCrossed

private

Definition at line 261 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_nTracks

int IsolatedTracksNxN::t_nTracks

private

Definition at line 228 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ t_nTrksJetCalo

std::vector<double>* IsolatedTracksNxN::t_nTrksJetCalo

private

Definition at line 248 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_nTrksJetVtx

std::vector<double> * IsolatedTracksNxN::t_nTrksJetVtx

private

Definition at line 248 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVisValid

std::vector<int>* IsolatedTracksNxN::t_PVisValid

private

Definition at line 244 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVndof

std::vector<int> * IsolatedTracksNxN::t_PVndof

private

Definition at line 244 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVNTracks

std::vector<int> * IsolatedTracksNxN::t_PVNTracks

private

Definition at line 244 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVNTracksHP

std::vector<int>* IsolatedTracksNxN::t_PVNTracksHP

private

Definition at line 245 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVNTracksHPWt

std::vector<int> * IsolatedTracksNxN::t_PVNTracksHPWt

private

Definition at line 245 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVNTracksWt

std::vector<int> * IsolatedTracksNxN::t_PVNTracksWt

private

Definition at line 244 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVTracksSumPt

std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPt

private

Definition at line 242 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVTracksSumPtHP

std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPtHP

private

Definition at line 243 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVTracksSumPtHPWt

std::vector<double> * IsolatedTracksNxN::t_PVTracksSumPtHPWt

private

Definition at line 243 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVTracksSumPtWt

std::vector<double>* IsolatedTracksNxN::t_PVTracksSumPtWt

private

Definition at line 243 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVx

std::vector<double>* IsolatedTracksNxN::t_PVx

private

Definition at line 242 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVy

std::vector<double> * IsolatedTracksNxN::t_PVy

private

Definition at line 242 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_PVz

std::vector<double> * IsolatedTracksNxN::t_PVz

private

Definition at line 242 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_RunNo

int IsolatedTracksNxN::t_RunNo

private

Definition at line 229 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ t_simTrackP

std::vector<double> * IsolatedTracksNxN::t_simTrackP

private

Definition at line 282 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackChiSq

std::vector<double>* IsolatedTracksNxN::t_trackChiSq

private

Definition at line 258 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackChiSqAll

std::vector<double> * IsolatedTracksNxN::t_trackChiSqAll

private

Definition at line 252 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDxy

std::vector<double>* IsolatedTracksNxN::t_trackDxy

private

Definition at line 256 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDxyAll

std::vector<double>* IsolatedTracksNxN::t_trackDxyAll

private

Definition at line 252 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDxyBS

std::vector<double> * IsolatedTracksNxN::t_trackDxyBS

private

Definition at line 256 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDxyPV

std::vector<double>* IsolatedTracksNxN::t_trackDxyPV

private

Definition at line 257 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDxyPVAll

std::vector<double> * IsolatedTracksNxN::t_trackDxyPVAll

private

Definition at line 252 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDz

std::vector<double> * IsolatedTracksNxN::t_trackDz

private

Definition at line 256 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDzAll

std::vector<double> * IsolatedTracksNxN::t_trackDzAll

private

Definition at line 252 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDzBS

std::vector<double> * IsolatedTracksNxN::t_trackDzBS

private

Definition at line 256 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDzPV

std::vector<double> * IsolatedTracksNxN::t_trackDzPV

private

Definition at line 257 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackDzPVAll

std::vector<double> * IsolatedTracksNxN::t_trackDzPVAll

private

Definition at line 252 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackEcalEta

std::vector<double>* IsolatedTracksNxN::t_trackEcalEta

private

Definition at line 255 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackEcalPhi

std::vector<double> * IsolatedTracksNxN::t_trackEcalPhi

private

Definition at line 255 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackEta

std::vector<double> * IsolatedTracksNxN::t_trackEta

private

Definition at line 254 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackEtaAll

std::vector<double> * IsolatedTracksNxN::t_trackEtaAll

private

Definition at line 250 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHcalEta

std::vector<double> * IsolatedTracksNxN::t_trackHcalEta

private

Definition at line 255 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHcalPhi

std::vector<double> * IsolatedTracksNxN::t_trackHcalPhi

private

Definition at line 255 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMeasTEC

std::vector<int> * IsolatedTracksNxN::t_trackHitInMeasTEC

private

Definition at line 267 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMeasTIB

std::vector<int> * IsolatedTracksNxN::t_trackHitInMeasTIB

private

Definition at line 267 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMeasTID

std::vector<int> * IsolatedTracksNxN::t_trackHitInMeasTID

private

Definition at line 267 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMeasTOB

std::vector<int>* IsolatedTracksNxN::t_trackHitInMeasTOB

private

Definition at line 267 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMissTEC

std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTEC

private

Definition at line 263 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMissTIB

std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTIB

private

Definition at line 263 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMissTIBTID

std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTIBTID

private

Definition at line 263 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMissTID

std::vector<int> * IsolatedTracksNxN::t_trackHitInMissTID

private

Definition at line 263 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitInMissTOB

std::vector<int>* IsolatedTracksNxN::t_trackHitInMissTOB

private

Definition at line 263 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMeasTEC

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMeasTEC

private

Definition at line 268 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMeasTIB

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMeasTIB

private

Definition at line 268 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMeasTID

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMeasTID

private

Definition at line 268 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMeasTOB

std::vector<int>* IsolatedTracksNxN::t_trackHitOutMeasTOB

private

Definition at line 268 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMissTEC

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTEC

private

Definition at line 265 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMissTIB

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTIB

private

Definition at line 265 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMissTID

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTID

private

Definition at line 265 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMissTOB

std::vector<int>* IsolatedTracksNxN::t_trackHitOutMissTOB

private

Definition at line 265 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitOutMissTOBTEC

std::vector<int> * IsolatedTracksNxN::t_trackHitOutMissTOBTEC

private

Definition at line 265 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitsTEC

std::vector<int> * IsolatedTracksNxN::t_trackHitsTEC

private

Definition at line 262 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackHitsTOB

std::vector<int>* IsolatedTracksNxN::t_trackHitsTOB

private

Definition at line 262 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackL

std::vector<double> * IsolatedTracksNxN::t_trackL

private

Definition at line 269 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackNOuterHits

std::vector<int> * IsolatedTracksNxN::t_trackNOuterHits

private

Definition at line 261 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackOutPosOutHitDr

std::vector<double>* IsolatedTracksNxN::t_trackOutPosOutHitDr

private

Definition at line 269 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackP

std::vector<double>* IsolatedTracksNxN::t_trackP

private

Definition at line 254 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPAll

std::vector<double>* IsolatedTracksNxN::t_trackPAll

private

Definition at line 250 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPdgIdAll

std::vector<double> * IsolatedTracksNxN::t_trackPdgIdAll

private

Definition at line 250 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPhi

std::vector<double> * IsolatedTracksNxN::t_trackPhi

private

Definition at line 254 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPhiAll

std::vector<double> * IsolatedTracksNxN::t_trackPhiAll

private

Definition at line 250 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPt

std::vector<double> * IsolatedTracksNxN::t_trackPt

private

Definition at line 254 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPtAll

std::vector<double>* IsolatedTracksNxN::t_trackPtAll

private

Definition at line 251 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trackPVIdx

std::vector<int>* IsolatedTracksNxN::t_trackPVIdx

private

Definition at line 259 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trkEcalEne

std::vector<double> * IsolatedTracksNxN::t_trkEcalEne

private

Definition at line 282 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ t_trkHcalEne

std::vector<double>* IsolatedTracksNxN::t_trkHcalEne

private

Definition at line 295 of file IsolatedTracksNxN.cc.

Referenced by analyze(), bookHistograms(), clearTreeVectors(), and ~IsolatedTracksNxN().

◆ tMaxE_

const double IsolatedTracksNxN::tMaxE_

private

Definition at line 176 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tMaxH_

const double IsolatedTracksNxN::tMaxH_

private

Definition at line 176 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tMinE_

const double IsolatedTracksNxN::tMinE_

private

Definition at line 176 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tMinH_

const double IsolatedTracksNxN::tMinH_

private

Definition at line 176 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_bs_

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

private

Definition at line 193 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_caloEB_

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksNxN::tok_caloEB_

private

Definition at line 199 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_caloEE_

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksNxN::tok_caloEE_

private

Definition at line 200 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_caloHH_

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksNxN::tok_caloHH_

private

Definition at line 201 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_caloTopology_

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

private

Definition at line 204 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksNxN::tok_EB_

private

Definition at line 195 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_ecalChStatus_

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

private

Definition at line 207 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksNxN::tok_EE_

private

Definition at line 196 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_genTrack_

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

private

Definition at line 191 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_geom_

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

private

Definition at line 203 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> IsolatedTracksNxN::tok_hbhe_

private

Definition at line 189 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_htmap_

edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> IsolatedTracksNxN::tok_htmap_

private

Definition at line 209 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_jets_

edm::EDGetTokenT<reco::CaloJetCollection> IsolatedTracksNxN::tok_jets_

private

Definition at line 188 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_L1extCenJet_

edm::EDGetTokenT<l1extra::L1JetParticleCollection> IsolatedTracksNxN::tok_L1extCenJet_

private

Definition at line 181 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_L1extFwdJet_

edm::EDGetTokenT<l1extra::L1JetParticleCollection> IsolatedTracksNxN::tok_L1extFwdJet_

private

Definition at line 182 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_L1extIsoEm_

edm::EDGetTokenT<l1extra::L1EmParticleCollection> IsolatedTracksNxN::tok_L1extIsoEm_

private

Definition at line 185 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_L1extMu_

edm::EDGetTokenT<l1extra::L1MuonParticleCollection> IsolatedTracksNxN::tok_L1extMu_

private

Definition at line 184 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_L1extNoIsoEm_

edm::EDGetTokenT<l1extra::L1EmParticleCollection> IsolatedTracksNxN::tok_L1extNoIsoEm_

private

Definition at line 186 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_L1extTauJet_

edm::EDGetTokenT<l1extra::L1JetParticleCollection> IsolatedTracksNxN::tok_L1extTauJet_

private

Definition at line 180 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_magField_

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

private

Definition at line 206 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_recVtx_

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

private

Definition at line 192 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_sevlv_

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

private

Definition at line 208 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_simTk_

edm::EDGetTokenT<edm::SimTrackContainer> IsolatedTracksNxN::tok_simTk_

private

Definition at line 197 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_simVtx_

edm::EDGetTokenT<edm::SimVertexContainer> IsolatedTracksNxN::tok_simVtx_

private

Definition at line 198 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ tok_topo_

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

private

Definition at line 205 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and IsolatedTracksNxN().

◆ trackerHitAssociatorConfig_

TrackerHitAssociator::Config IsolatedTracksNxN::trackerHitAssociatorConfig_

private

Definition at line 163 of file IsolatedTracksNxN.cc.

Referenced by analyze().

◆ tree_

TTree* IsolatedTracksNxN::tree_

private

Definition at line 226 of file IsolatedTracksNxN.cc.

Referenced by analyze(), and bookHistograms().

◆ writeAllTracks_

const bool IsolatedTracksNxN::writeAllTracks_

private

Definition at line 169 of file IsolatedTracksNxN.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ IsolatedTracksNxN()

◆ ~IsolatedTracksNxN()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ bookHistograms()

◆ clearTreeVectors()

◆ endJob()

◆ fillDescriptions()

◆ printTrack()

Member Data Documentation

◆ debugL1Info_

◆ debugTrks_

◆ doMC_

◆ genPartEtaBins

◆ genPartPBins

◆ h_L1AlgoNames

◆ h_maxNearP15x15

◆ h_maxNearP21x21

◆ h_maxNearP25x25

◆ h_maxNearP31x31

◆ h_nTracks

◆ h_PVTracksWt

◆ h_recEta_0

◆ h_recEta_1

◆ h_recEta_2

◆ h_recEtaP_0

◆ h_recEtaP_1

◆ h_recEtaP_2

◆ h_recEtaPt_0

◆ h_recEtaPt_1

◆ h_recEtaPt_2

◆ h_recP_0

◆ h_recP_1

◆ h_recP_2

◆ h_recPhi_0

◆ h_recPhi_1

◆ h_recPhi_2

◆ h_recPt_0

◆ h_recPt_1

◆ h_recPt_2

◆ initL1_

◆ l1AlgoMap_

◆ L1TriggerAlgoInfo_

◆ m_l1GtUtils

◆ m_triggerMaskAlgoTrig

◆ maxTrackEta_

◆ minTrackP_

◆ myverbose_

◆ nbad_

◆ NEtaBins

◆ nEventProc_

◆ nL1BitsMax

◆ NPBins

◆ printTrkHitPattern_

◆ pvTracksPtMin_

◆ t_Bunch

◆ t_e11x11

◆ t_e11x11_10Sig

◆ t_e11x11_15Sig

◆ t_e11x11_20Sig

◆ t_e11x11_25Sig

◆ t_e11x11_30Sig

◆ t_e15x15

◆ t_e15x15_10Sig

◆ t_e15x15_15Sig

◆ t_e15x15_20Sig

◆ t_e15x15_25Sig

◆ t_e15x15_30Sig

◆ t_e7x7

◆ t_e7x7_10Sig

◆ t_e7x7_15Sig