IsolatedTracksNxN Class Reference

#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
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
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 ()
double	deltaPhi (double v1, double v2)
double	deltaR (double eta1, double phi1, double eta2, double phi2)
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_
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::EDGetTokenT< EcalRecHitCollection >	tok_EB_
edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_
edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_
edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_
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::EDGetTokenT< reco::VertexCollection >	tok_recVtx_
edm::EDGetTokenT< edm::SimTrackContainer >	tok_simTk_
edm::EDGetTokenT< edm::SimVertexContainer >	tok_simVtx_
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
TTree *	tree_
const bool	writeAllTracks_

Static Private Attributes
static size_t	NEtaBins = 3
static size_t	nL1BitsMax =128
static 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)
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<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: <one line="" class="" summary>="">

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

Definition at line 139 of file IsolatedTracksNxN.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 305 of file IsolatedTracksNxN.cc.

References debugL1Info_, doMC_, edm::ParameterSet::getParameter(), TFileService::kSharedResource, L1TriggerAlgoInfo_, maxTrackEta_, minTrackP_, myverbose_, tMaxE_, tMaxH_, tMinE_, tMinH_, tok_bs_, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_jets_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extIsoEm_, tok_L1extMu_, tok_L1extNoIsoEm_, tok_L1extTauJet_, tok_recVtx_, tok_simTk_, and tok_simVtx_.

                                                                   :
  m_l1GtUtils(iConfig, consumesCollector(), useL1GtTriggerMenuLite, *this),
  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) {

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

IsolatedTracksNxN::~IsolatedTracksNxN ( )

override

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

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

overrideprivate

Definition at line 650 of file IsolatedTracksNxN.cc.

References funct::abs(), L1GtUtils::availableL1Configuration(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationEcal(), spr::chargeIsolationHcal(), clearTreeVectors(), debugL1Info_, debugTrks_, TauDecayModes::dec, L1GtUtils::decision(), deltaR(), diffTreeTool::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(), spr::caloSimInfo::eTotal, edm::EventID::event(), plotBeamSpotDB::first, genPartEtaBins, genPartPBins, edm::EventSetup::get(), edm::RefToBase< T >::get(), edm::Event::getByToken(), L1GtUtils::getL1GtRunCache(), 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, edm::EventBase::id(), cuy::ii, initL1_, edm::HandleBase::isValid(), fwrapper::jets, l1AlgoMap_, L1TriggerAlgoInfo_, edm::EventBase::luminosityBlock(), 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, nTracks(), reco::Track::outerPosition(), reco::TrackBase::p(), p1, spr::caloSimInfo::pdgMatched, phi, reco::BeamSpot::position(), L1GtUtils::prescaleFactor(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), spr::propagateCALO(), spr::propagateTrackerEnd(), EnergyCorrector::pt, reco::TrackBase::pt(), hiDetachedQuadStep_cff::pt1, L1GtUtils::ptrL1TriggerMenuEventSetup(), pvTracksPtMin_, reco::TrackBase::quality(), reco::TrackBase::qualityByName(), edm::EventID::run(), 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_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_jets_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extIsoEm_, tok_L1extMu_, tok_L1extNoIsoEm_, tok_L1extTauJet_, tok_recVtx_, tok_simTk_, tok_simVtx_, trackerHitAssociatorConfig_, reco::HitPattern::trackerLayersWithMeasurement(), tree_, useL1EventSetup(), useL1GtTriggerMenuLite(), writeAllTracks_, x, y, and z.

                                                                                   {
  
  bool haveIsoTrack=false;  

  edm::ESHandle<MagneticField> bFieldH;
  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
  const MagneticField* bField = bFieldH.product();
  
  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
  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  const CaloGeometry* geo = pG.product();
  
  edm::ESHandle<CaloTopology> theCaloTopology;
  iSetup.get<CaloTopologyRecord>().get(theCaloTopology); 
  const CaloTopology *caloTopology = theCaloTopology.product();
  
  edm::ESHandle<HcalTopology> htopo;
  iSetup.get<HcalRecNumberingRecord>().get(htopo);
  const HcalTopology* theHBHETopology = htopo.product();
  
  edm::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
  edm::ESHandle<EcalChannelStatus> ecalChStatus;
  iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
  const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();

  // Retrieve trigger tower map
  edm::ESHandle<EcalTrigTowerConstituentsMap> hTtmap;
  iSetup.get<IdealGeometryRecord>().get(hTtmap);
  const EcalTrigTowerConstituentsMap& ttMap = *hTtmap;

  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.reset(new 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();
    }

    edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
    iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);

    // 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.product(),3,3,   -100.0, -100.0, tMinE_,tMaxE_);
    e9x9P         = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4,   -100.0, -100.0, tMinE_,tMaxE_);
    e11x11P       = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5,   -100.0, -100.0, tMinE_,tMaxE_);
    e15x15P       = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7,   -100.0, -100.0, tMinE_,tMaxE_);
    
    e7x7_10SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3,   0.030,  0.150, tMinE_,tMaxE_);
    e9x9_10SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4,   0.030,  0.150, tMinE_,tMaxE_);
    e11x11_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5,   0.030,  0.150, tMinE_,tMaxE_);
    e15x15_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7,   0.030,  0.150, tMinE_,tMaxE_);

    e7x7_15SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),ttMap, 3,3, 0.20,0.45, tMinE_,tMaxE_);
    e9x9_15SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),ttMap, 4,4, 0.20,0.45, tMinE_,tMaxE_);
    e11x11_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(), ttMap, 5,5, 0.20,0.45, tMinE_,tMaxE_);
    e15x15_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),ttMap, 7,7, 0.20,0.45, tMinE_,tMaxE_, false);
    
    e7x7_20SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3,   0.060,  0.300, tMinE_,tMaxE_);
    e9x9_20SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4,   0.060,  0.300, tMinE_,tMaxE_);
    e11x11_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5,   0.060,  0.300, tMinE_,tMaxE_);
    e15x15_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7,   0.060,  0.300, tMinE_,tMaxE_);
    
    e7x7_25SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3,   0.075,  0.375, tMinE_,tMaxE_);
    e9x9_25SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4,   0.075,  0.375, tMinE_,tMaxE_);
    e11x11_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5,   0.075,  0.375, tMinE_,tMaxE_);
    e15x15_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7,   0.075,  0.375, tMinE_,tMaxE_);
    
    e7x7_30SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3,   0.090,  0.450, tMinE_,tMaxE_);
    e9x9_30SigP   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4,   0.090,  0.450, tMinE_,tMaxE_);
    e11x11_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5,   0.090,  0.450, tMinE_,tMaxE_);
    e15x15_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),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 (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();
}

void IsolatedTracksNxN::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1590 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();
}

void IsolatedTracksNxN::bookHistograms ( )

private

Definition at line 1831 of file IsolatedTracksNxN.cc.

References dir, 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, mps_fire::i, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), 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");

}