IsolatedTracksCone Class Reference

#include <IsolatedTracksCone.h>

Inheritance diagram for IsolatedTracksCone:

Public Member Functions
	IsolatedTracksCone (const edm::ParameterSet &)
	~IsolatedTracksCone () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Public Attributes
double	genPartEtaBins [5]
double	genPartPBins [22]

Static Public Attributes
static const int	NEtaBins = 4
static const int	NPBins = 21
static const int	NPTBins = 21

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	BuildTree ()
void	clearTrackVectors ()
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
bool	applyEcalIsolation_
int	debugEcalSimInfo_
int	debugTrks_
bool	doMC
double	drLeadJetVeto_
edm::Service< TFileService >	fs
TH1F *	hRawEta
TH1F *	hRawP
TH1F *	hRawPhi
TH1F *	hRawPt
double	leadL1JetEta
double	leadL1JetPhi
double	leadL1JetPT
double	maxNearTrackP_
double	maxTrackEta_
double	minTrackP_
int	myverbose_
int	nEVT
int	nEVT_failL1
int	nFailEta
int	nFailHighPurityQaul
int	nFailPt
int	nMissEcal
int	nMissHcal
int	nRawTRK
TTree *	ntp
int	nTRK
TrackAssociatorParameters	parameters_
bool	printTrkHitPattern_
double	ptMinLeadJet_
std::vector< double > *	t_distFromHotCell_h3x3
std::vector< double > *	t_distFromHotCell_h5x5
std::vector< double > *	t_dtFromLeadJet
std::vector< double > *	t_e3x3
std::vector< double > *	t_h3x3
std::vector< double > *	t_h5x5
std::vector< int > *	t_hltDoubleLooseIsoTau15_Trk5
std::vector< int > *	t_hltHB
std::vector< int > *	t_hltHE
std::vector< int > *	t_hltJet110
std::vector< int > *	t_hltJet140
std::vector< int > *	t_hltJet180
std::vector< int > *	t_hltJet30
std::vector< int > *	t_hltJet50
std::vector< int > *	t_hltJet80
std::vector< int > *	t_hltL1Jet15
std::vector< int > *	t_hltL1SingleEG5
std::vector< int > *	t_hltMinBiasEcal
std::vector< int > *	t_hltMinBiasHcal
std::vector< int > *	t_hltMinBiasPixel
std::vector< int > *	t_hltSingleIsoTau30_Trk5
std::vector< int > *	t_hltZeroBias
std::vector< double > *	t_hsim3x3
std::vector< double > *	t_hsim3x3CharHad
std::vector< double > *	t_hsim3x3Matched
std::vector< int > *	t_hsim3x3NCharHad
std::vector< double > *	t_hsim3x3NeutHad
std::vector< int > *	t_hsim3x3NMatched
std::vector< int > *	t_hsim3x3NNeutHad
std::vector< int > *	t_hsim3x3NPhoton
std::vector< int > *	t_hsim3x3NRest
std::vector< int > *	t_hsim3x3NTotal
std::vector< double > *	t_hsim3x3PdgMatched
std::vector< double > *	t_hsim3x3Photon
std::vector< double > *	t_hsim3x3Rest
std::vector< double > *	t_hsim3x3Total
std::vector< double > *	t_hsim5x5
std::vector< double > *	t_hsim5x5CharHad
std::vector< double > *	t_hsim5x5Matched
std::vector< int > *	t_hsim5x5NCharHad
std::vector< double > *	t_hsim5x5NeutHad
std::vector< int > *	t_hsim5x5NMatched
std::vector< int > *	t_hsim5x5NNeutHad
std::vector< int > *	t_hsim5x5NPhoton
std::vector< int > *	t_hsim5x5NRest
std::vector< int > *	t_hsim5x5NTotal
std::vector< double > *	t_hsim5x5PdgMatched
std::vector< double > *	t_hsim5x5Photon
std::vector< double > *	t_hsim5x5Rest
std::vector< double > *	t_hsim5x5Total
std::vector< int > *	t_ietaFromHotCell_h3x3
std::vector< int > *	t_ietaFromHotCell_h5x5
std::vector< unsigned int > *	t_ievt
std::vector< unsigned int > *	t_ilum
std::vector< int > *	t_iphiFromHotCell_h3x3
std::vector< int > *	t_iphiFromHotCell_h5x5
std::vector< unsigned int > *	t_irun
std::vector< double > *	t_nRH_h3x3
std::vector< double > *	t_nRH_h5x5
std::vector< double > *	t_simP
std::vector< double > *	t_trkEcalEne
std::vector< double > *	t_trkEta
std::vector< double > *	t_trkHcalEne
std::vector< double > *	t_trkNLayersCrossed
std::vector< double > *	t_trkNOuterHits
std::vector< double > *	t_trkP
std::vector< double > *	t_trkPhi
std::vector< double > *	t_trkPt
std::vector< std::vector< double > > *	t_v_cone_hmaxNearP
std::vector< std::vector< double > > *	t_v_cone_hmaxNearP_goodTrk
std::vector< std::vector< int > > *	t_v_cone_hnLayers_maxNearP
std::vector< std::vector< int > > *	t_v_cone_hnNearTRKs
std::vector< std::vector< int > > *	t_v_cone_htrkQual_maxNearP
std::vector< std::vector< double > > *	t_v_distFromHotCell
std::vector< std::vector< double > > *	t_v_eDR
std::vector< std::vector< double > > *	t_v_eMipDR
std::vector< std::vector< double > > *	t_v_hCone
std::vector< std::vector< int > > *	t_v_hlTriggers
std::vector< std::vector< double > > *	t_v_hmaxNearP
std::vector< std::vector< double > > *	t_v_hmaxNearP_goodTrk
std::vector< std::vector< int > > *	t_v_hnLayers_maxNearP
std::vector< std::vector< int > > *	t_v_hnNearTRKs
std::vector< std::vector< double > > *	t_v_hsimCone
std::vector< std::vector< double > > *	t_v_hsimInfoConeCharHad
std::vector< std::vector< double > > *	t_v_hsimInfoConeMatched
std::vector< std::vector< int > > *	t_v_hsimInfoConeNCharHad
std::vector< std::vector< double > > *	t_v_hsimInfoConeNeutHad
std::vector< std::vector< int > > *	t_v_hsimInfoConeNMatched
std::vector< std::vector< int > > *	t_v_hsimInfoConeNNeutHad
std::vector< std::vector< int > > *	t_v_hsimInfoConeNPhoton
std::vector< std::vector< int > > *	t_v_hsimInfoConeNRest
std::vector< std::vector< int > > *	t_v_hsimInfoConeNTotal
std::vector< std::vector< double > > *	t_v_hsimInfoConePdgMatched
std::vector< std::vector< double > > *	t_v_hsimInfoConePhoton
std::vector< std::vector< double > > *	t_v_hsimInfoConeRest
std::vector< std::vector< double > > *	t_v_hsimInfoConeTotal
std::vector< std::vector< int > > *	t_v_htrkQual_maxNearP
std::vector< std::vector< int > > *	t_v_ietaFromHotCell
std::vector< std::vector< int > > *	t_v_iphiFromHotCell
std::vector< std::vector< int > > *	t_v_nRecHitsCone
std::vector< std::vector< int > > *	t_v_nSimHitsCone
std::vector< std::vector< double > > *	t_v_RH_h3x3_ene
std::vector< std::vector< int > > *	t_v_RH_h3x3_ieta
std::vector< std::vector< int > > *	t_v_RH_h3x3_iphi
std::vector< std::vector< double > > *	t_v_RH_h5x5_ene
std::vector< std::vector< int > > *	t_v_RH_h5x5_ieta
std::vector< std::vector< int > > *	t_v_RH_h5x5_iphi
std::vector< std::vector< double > > *	t_v_RH_r26_ene
std::vector< std::vector< int > > *	t_v_RH_r26_ieta
std::vector< std::vector< int > > *	t_v_RH_r26_iphi
std::vector< std::vector< double > > *	t_v_RH_r44_ene
std::vector< std::vector< int > > *	t_v_RH_r44_ieta
std::vector< std::vector< int > > *	t_v_RH_r44_iphi
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< l1extra::L1JetParticleCollection >	tok_L1extCenJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection >	tok_L1extFwdJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection >	tok_L1extTauJet_
edm::EDGetTokenT< edm::SimTrackContainer >	tok_simTk_
edm::EDGetTokenT< edm::SimVertexContainer >	tok_simVtx_
edm::EDGetTokenT< edm::TriggerResults >	tok_trigger_
TrackDetectorAssociator *	trackAssociator_
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
bool	useJetTrigger_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Analysis/IsolatedTracksCone/src/IsolatedTracksCone.cc

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

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

Definition at line 94 of file IsolatedTracksCone.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 37 of file IsolatedTracksCone.cc.

References applyEcalIsolation_, edm::EDConsumerBase::consumesCollector(), gather_cfg::cout, debugEcalSimInfo_, debugTrks_, doMC, drLeadJetVeto_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), TrackAssociatorParameters::loadParameters(), maxNearTrackP_, maxTrackEta_, minTrackP_, myverbose_, metProducer_cfi::parameters, parameters_, printTrkHitPattern_, ptMinLeadJet_, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extTauJet_, tok_simTk_, tok_simVtx_, tok_trigger_, trackAssociator_, TrackDetectorAssociator::useDefaultPropagator(), and useJetTrigger_.

                                                                     :
   trackerHitAssociatorConfig_(consumesCollector()) {

  //now do what ever initialization is needed
  doMC            = iConfig.getUntrackedParameter<bool>  ("DoMC", false); 
  myverbose_      = iConfig.getUntrackedParameter<int>( "Verbosity", 5 );
  useJetTrigger_  = iConfig.getUntrackedParameter<bool>( "useJetTrigger", false);
  drLeadJetVeto_  = iConfig.getUntrackedParameter<double>( "drLeadJetVeto",  1.2 );
  ptMinLeadJet_   = iConfig.getUntrackedParameter<double>( "ptMinLeadJet",  15.0 );

  debugTrks_          = iConfig.getUntrackedParameter<int>("DebugTracks");
  printTrkHitPattern_ = iConfig.getUntrackedParameter<bool>("PrintTrkHitPattern");
  
  minTrackP_      = iConfig.getUntrackedParameter<double>( "minTrackP", 10.0);
  maxTrackEta_    = iConfig.getUntrackedParameter<double>( "maxTrackEta", 5.0);
  maxNearTrackP_  = iConfig.getUntrackedParameter<double>( "maxNearTrackP", 1.0);

  debugEcalSimInfo_   = iConfig.getUntrackedParameter<int>("DebugEcalSimInfo");

  applyEcalIsolation_ = iConfig.getUntrackedParameter<bool>("ApplyEcalIsolation");

  tok_L1extTauJet_ = consumes<l1extra::L1JetParticleCollection>(iConfig.getParameter<edm::InputTag>("L1extraTauJetSource"));
  tok_L1extCenJet_ = consumes<l1extra::L1JetParticleCollection>(iConfig.getParameter<edm::InputTag>("L1extraCenJetSource"));
  tok_L1extFwdJet_ = consumes<l1extra::L1JetParticleCollection>(iConfig.getParameter<edm::InputTag>("L1extraFwdJetSource"));

  // hard coded collection access
  tok_EB_   = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
  tok_EE_   = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
  tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
  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"));
  tok_trigger_  = consumes<edm::TriggerResults>(edm::InputTag("TriggerResults","","HLT"));

  edm::ParameterSet parameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
  edm::ConsumesCollector iC = consumesCollector();
  parameters_.loadParameters( parameters, iC );
  trackAssociator_ =  new TrackDetectorAssociator();
  trackAssociator_->useDefaultPropagator();

  if(myverbose_>=0) {
    std::cout <<"Parameters read from config file \n" 
          << "myverbose_ "          << myverbose_     << "\n"     
          << "useJetTrigger_ "      << useJetTrigger_ << "\n"
          << "drLeadJetVeto_ "      << drLeadJetVeto_ << "\n"
          << "minTrackP_ "         << minTrackP_    << "\n"
          << "maxTrackEta_ "        << maxTrackEta_   << "\n"
          << "maxNearTrackP_ "      << maxNearTrackP_ 
          << std::endl;
  }
}

IsolatedTracksCone::~IsolatedTracksCone ( )

override

Definition at line 92 of file IsolatedTracksCone.cc.

References trackAssociator_.

                                        {
  delete  trackAssociator_;
}

Member Function Documentation

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

overrideprivate

Definition at line 96 of file IsolatedTracksCone.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), TrackDetectorAssociator::associate(), clearTrackVectors(), spr::coneChargeIsolation(), gather_cfg::cout, debugTrks_, DeltaPhi(), doMC, DetId::Ecal, EcalBarrel, EcalEndcap, spr::eCaloSimInfo(), spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), spr::eHCALmatrix(), spr::eHCALSimInfo(), reco::TrackBase::eta(), edm::EventID::event(), plotBeamSpotDB::first, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getClosestCell(), spr::getDistInPlaneTrackDir(), TrackDetectorAssociator::getFreeTrajectoryState(), edm::Event::getLuminosityBlock(), CaloGeometry::getSubdetectorGeometry(), photonIsolationHIProducer_cfi::hbhe, DetId::Hcal, HcalBarrel, HcalEndcap, reco::TrackBase::hitPattern(), hRawEta, hRawP, hRawPhi, hRawPt, mps_fire::i, edm::EventBase::id(), createfilelist::int, edm::HandleBase::isValid(), leadL1JetEta, leadL1JetPhi, leadL1JetPT, edm::LuminosityBlockBase::luminosityBlock(), spr::matchedSimTrack(), maxTrackEta_, minTrackP_, reco::TrackBase::momentum(), nEVT, nEVT_failL1, nFailEta, nFailHighPurityQaul, nFailPt, nMissEcal, nMissHcal, nRawTRK, ntp, nTRK, reco::TrackBase::p(), p1, parameters_, edm::Handle< T >::product(), edm::ESHandle< T >::product(), reco::TrackBase::pt(), hiDetachedQuadStep_cff::pt1, ptMinLeadJet_, reco::TrackBase::quality(), reco::TrackBase::qualityByName(), DetId::rawId(), edm::EventID::run(), edm::HLTGlobalStatus::size(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, t_distFromHotCell_h3x3, t_distFromHotCell_h5x5, t_dtFromLeadJet, t_e3x3, t_h3x3, t_h5x5, t_hltDoubleLooseIsoTau15_Trk5, t_hltHB, t_hltHE, t_hltJet110, t_hltJet140, t_hltJet180, t_hltJet30, t_hltJet50, t_hltJet80, t_hltL1Jet15, t_hltL1SingleEG5, t_hltMinBiasEcal, t_hltMinBiasHcal, t_hltMinBiasPixel, t_hltSingleIsoTau30_Trk5, t_hltZeroBias, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NCharHad, t_hsim3x3NeutHad, t_hsim3x3NMatched, t_hsim3x3NNeutHad, t_hsim3x3NPhoton, t_hsim3x3NRest, t_hsim3x3NTotal, t_hsim3x3PdgMatched, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim3x3Total, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NCharHad, t_hsim5x5NeutHad, t_hsim5x5NMatched, t_hsim5x5NNeutHad, t_hsim5x5NPhoton, t_hsim5x5NRest, t_hsim5x5NTotal, t_hsim5x5PdgMatched, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim5x5Total, t_ietaFromHotCell_h3x3, t_ietaFromHotCell_h5x5, t_ievt, t_ilum, t_iphiFromHotCell_h3x3, t_iphiFromHotCell_h5x5, t_irun, t_nRH_h3x3, t_nRH_h5x5, t_simP, t_trkEcalEne, t_trkEta, t_trkHcalEne, t_trkNLayersCrossed, t_trkNOuterHits, t_trkP, t_trkPhi, t_trkPt, t_v_cone_hmaxNearP, t_v_cone_hmaxNearP_goodTrk, t_v_cone_hnLayers_maxNearP, t_v_cone_hnNearTRKs, t_v_cone_htrkQual_maxNearP, t_v_distFromHotCell, t_v_eDR, t_v_eMipDR, t_v_hCone, t_v_hlTriggers, t_v_hmaxNearP, t_v_hmaxNearP_goodTrk, t_v_hnLayers_maxNearP, t_v_hnNearTRKs, t_v_hsimCone, t_v_hsimInfoConeCharHad, t_v_hsimInfoConeMatched, t_v_hsimInfoConeNCharHad, t_v_hsimInfoConeNeutHad, t_v_hsimInfoConeNMatched, t_v_hsimInfoConeNNeutHad, t_v_hsimInfoConeNPhoton, t_v_hsimInfoConeNRest, t_v_hsimInfoConeNTotal, t_v_hsimInfoConePdgMatched, t_v_hsimInfoConePhoton, t_v_hsimInfoConeRest, t_v_hsimInfoConeTotal, t_v_htrkQual_maxNearP, t_v_ietaFromHotCell, t_v_iphiFromHotCell, t_v_nRecHitsCone, t_v_nSimHitsCone, t_v_RH_h3x3_ene, t_v_RH_h3x3_ieta, t_v_RH_h3x3_iphi, t_v_RH_h5x5_ene, t_v_RH_h5x5_ieta, t_v_RH_h5x5_iphi, t_v_RH_r26_ene, t_v_RH_r26_ieta, t_v_RH_r26_iphi, t_v_RH_r44_ene, t_v_RH_r44_ieta, t_v_RH_r44_iphi, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extTauJet_, tok_simTk_, tok_simVtx_, tok_trigger_, trackAssociator_, trackerHitAssociatorConfig_, reco::HitPattern::trackerLayersWithMeasurement(), edm::TriggerNames::triggerName(), edm::Event::triggerNames(), triggerResults, TrackDetMatchInfo::trkGlobPosAtEcal, TrackDetMatchInfo::trkGlobPosAtHcal, TrackDetMatchInfo::trkMomAtEcal, TrackDetMatchInfo::trkMomAtHcal, and useJetTrigger_.

                                              {
  
  unsigned int irun = (unsigned int)iEvent.id().run();
  unsigned int ilum = (unsigned int)iEvent.getLuminosityBlock().luminosityBlock();
  unsigned int ievt = (unsigned int)iEvent.id().event();

  
  
  clearTrackVectors();
  
  // check the L1 objects
  bool   L1Pass = false;
  leadL1JetPT=-999, leadL1JetEta=-999,  leadL1JetPhi=-999;
  if( !useJetTrigger_ ) {
    L1Pass = true;
  } else {
    edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
    iEvent.getByToken(tok_L1extTauJet_,l1TauHandle);
    l1extra::L1JetParticleCollection::const_iterator itr;
    for(itr = l1TauHandle->begin(); itr != l1TauHandle->end(); ++itr ) 
    {
      if( itr->pt()>leadL1JetPT ) {
    leadL1JetPT  = itr->pt();
    leadL1JetEta = itr->eta();
    leadL1JetPhi = itr->phi();
      }
    }
    edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
    iEvent.getByToken(tok_L1extCenJet_,l1CenJetHandle);
    for( itr = l1CenJetHandle->begin();  itr != l1CenJetHandle->end(); ++itr ) 
    {
      if( itr->pt()>leadL1JetPT ) {
    leadL1JetPT  = itr->pt();
    leadL1JetEta = itr->eta();
    leadL1JetPhi = itr->phi();
      }
    }
    edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
    iEvent.getByToken(tok_L1extFwdJet_,l1FwdJetHandle);
    for( itr = l1FwdJetHandle->begin();  itr != l1FwdJetHandle->end(); ++itr ) 
    {
      if( itr->pt()>leadL1JetPT ) {
    leadL1JetPT  = itr->pt();
    leadL1JetEta = itr->eta();
    leadL1JetPhi = itr->phi();
      }
    }
    if(leadL1JetPT>ptMinLeadJet_) 
    {
      L1Pass = true;
    }
    
  }
  

  // Break now if L1Pass is false
  if (!L1Pass) {
    nEVT_failL1++;
  //     //    std::cout << "L1Pass is false : " << L1Pass << std::endl;
  //     return;  
  }
  
  // Get the collection handles
  

  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  const CaloGeometry* geo = pG.product();
  const CaloSubdetectorGeometry* gEB = 
    geo->getSubdetectorGeometry(DetId::Ecal,EcalBarrel);
  const CaloSubdetectorGeometry* gEE = 
    geo->getSubdetectorGeometry(DetId::Ecal,EcalEndcap);
  const CaloSubdetectorGeometry* gHB = 
    geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
  const CaloSubdetectorGeometry* gHE = 
    geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap);
  
  edm::ESHandle<CaloTopology> theCaloTopology;
  iSetup.get<CaloTopologyRecord>().get(theCaloTopology); 
  const CaloTopology *caloTopology = theCaloTopology.product();
  //  const CaloSubdetectorTopology* theEBTopology   = theCaloTopology->getSubdetectorTopology(DetId::Ecal,EcalBarrel);
  //  const CaloSubdetectorTopology* theEETopology   = theCaloTopology->getSubdetectorTopology(DetId::Ecal,EcalEndcap);
  
  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();
  
  edm::Handle<HBHERecHitCollection> hbhe;
  iEvent.getByToken(tok_hbhe_,hbhe);
  const HBHERecHitCollection Hithbhe = *(hbhe.product());

  edm::Handle<reco::TrackCollection> trkCollection;
  iEvent.getByToken(tok_genTrack_, trkCollection);
  reco::TrackCollection::const_iterator trkItr;
  if(debugTrks_>1){
    std::cout << "Track Collection: " << std::endl;
    std::cout << "Number of Tracks " << trkCollection->size() << std::endl;
  }
  std::string theTrackQuality = "highPurity";
  reco::TrackBase::TrackQuality trackQuality_=
    reco::TrackBase::qualityByName(theTrackQuality);
  
  //get Handles to SimTracks and SimHits
  edm::Handle<edm::SimTrackContainer> SimTk;
  if (doMC) iEvent.getByToken(tok_simTk_,SimTk);
  edm::SimTrackContainer::const_iterator simTrkItr;

  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);
  
  
  
  // Get HLT_IsoTrackHB/HE Information
    
  edm::Handle<edm::TriggerResults> triggerResults;
  iEvent.getByToken( tok_trigger_, triggerResults);

  

  std::vector<int> v_hlTriggers;
  int hltHB(-99);
  int hltHE(-99);
  int hltL1Jet15                (-99);
  int hltJet30              (-99);
  int hltJet50              (-99);
  int hltJet80              (-99);
  int hltJet110             (-99);
  int hltJet140             (-99);
  int hltJet180             (-99);
  int hltL1SingleEG5        (-99);
  int hltZeroBias       (-99);
  int hltMinBiasHcal        (-99);
  int hltMinBiasEcal        (-99);
  int hltMinBiasPixel           (-99);
  int hltSingleIsoTau30_Trk5    (-99);
  int hltDoubleLooseIsoTau15_Trk5(-99);

  if (triggerResults.isValid()) {
    
    const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
    // TriggerNames class  triggerNames.init(*triggerResults);
    

    for (unsigned int i=0; i<triggerResults->size(); i++){
      //      std::cout << "triggerNames.triggerName(" << i << ") = " << triggerNames.triggerName(i) << std::endl;
      if (triggerNames.triggerName(i) == "HLT_IsoTrackHE_1E31") hltHE = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_IsoTrackHB_1E31") hltHB = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_L1Jet15"              ) hltL1Jet15                 = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_Jet30"                ) hltJet30               = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_Jet50"                ) hltJet50               = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_Jet80"                ) hltJet80               = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_Jet110"               ) hltJet110              = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_Jet140"               ) hltJet140              = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_Jet180"               ) hltJet180              = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_L1SingleEG5"      ) hltL1SingleEG5         = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_ZeroBias"             ) hltZeroBias                = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_MinBiasHcal"      ) hltMinBiasHcal         = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_MinBiasEcal"      ) hltMinBiasEcal         = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_MinBiasPixel"         ) hltMinBiasPixel            = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_SingleIsoTau30_Trk5"  ) hltSingleIsoTau30_Trk5     = triggerResults->accept(i);
      if (triggerNames.triggerName(i) == "HLT_DoubleLooseIsoTau15_Trk5" ) hltDoubleLooseIsoTau15_Trk5    = triggerResults->accept(i);
    }
  }
  
    

  
  // Primary loop over tracks
  std::unique_ptr<TrackerHitAssociator> associate;
  if (doMC) associate.reset(new TrackerHitAssociator(iEvent, trackerHitAssociatorConfig_));
  

  nTRK      = 0;
  nRawTRK   = 0;
  nFailPt   = 0;
  nFailEta  = 0;
  nFailHighPurityQaul = 0;
  nMissEcal = 0;
  nMissHcal = 0;

  for( trkItr = trkCollection->begin(); 
       trkItr != trkCollection->end(); ++trkItr)
  {

    nRawTRK++;

    const reco::Track* pTrack = &(*trkItr);

    // Check for min Pt and max Eta P

    bool trkQual  = pTrack->quality(trackQuality_);
    bool goodPt   = pTrack->p()>minTrackP_;
    bool goodEta  = std::abs(pTrack->momentum().eta())<maxTrackEta_;

    double eta1       = pTrack->momentum().eta();
    double phi1       = pTrack->momentum().phi();
    double pt1        = pTrack->pt();
    double p1         = pTrack->p();


    if (!goodEta){
      nFailEta++;
    }
    if (!goodPt){
      nFailPt++;
    }
    if (!trkQual){
      nFailHighPurityQaul++;
    }
    
    hRawPt ->Fill(pt1 );
    hRawP  ->Fill(p1  );
    hRawEta->Fill(eta1);
    hRawPhi->Fill(phi1);
      
    if( !goodEta || !goodPt || !trkQual ) continue; // Skip to next track
    
    // Find track trajectory
    
      
    const FreeTrajectoryState fts1 = 
      trackAssociator_->getFreeTrajectoryState(iSetup, *pTrack);
    
      
    TrackDetMatchInfo info1 = 
      trackAssociator_->associate(iEvent, iSetup, fts1, parameters_);
    
      

    // First confirm track makes it to Hcal

    if (info1.trkGlobPosAtHcal.x()==0 && 
    info1.trkGlobPosAtHcal.y()==0 && 
    info1.trkGlobPosAtHcal.z()==0) 
    {
      nMissHcal++;
      continue;      
    }
    
    const GlobalPoint hpoint1(info1.trkGlobPosAtHcal.x(),
                  info1.trkGlobPosAtHcal.y(),
                  info1.trkGlobPosAtHcal.z());

      

    // Get basic quantities

    const reco::HitPattern& hitp = pTrack->hitPattern();
    int nLayersCrossed = hitp.trackerLayersWithMeasurement();        
    int nOuterHits     = hitp.stripTOBLayersWithMeasurement()
      +hitp.stripTECLayersWithMeasurement() ;

    
    double simP = 0;
    if (doMC) {
      edm::SimTrackContainer::const_iterator matchedSimTrk = 
    spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
      simP = matchedSimTrk->momentum().P();
    }
    // Get Ecal Point

    const GlobalPoint point1(info1.trkGlobPosAtEcal.x(),
                 info1.trkGlobPosAtEcal.y(),
                 info1.trkGlobPosAtEcal.z());
    

    // Sanity check that track hits Ecal

    if (info1.trkGlobPosAtEcal.x()==0 && 
    info1.trkGlobPosAtEcal.y()==0 && 
    info1.trkGlobPosAtEcal.z()==0) 
    {
      std::cout << "Track doesn't reach Ecal." << std::endl;
      nMissEcal++;
      continue;
    }

    // Get Track Momentum - make sure you have latest version of TrackDetMatchInfo
    
    GlobalVector trackMomAtEcal = info1.trkMomAtEcal;
    GlobalVector trackMomAtHcal = info1.trkMomAtHcal;

    // If using Jet trigger, get distance from leading jet

    double drFromLeadJet = 999.0;
    if( useJetTrigger_ ) {
      double dphi = DeltaPhi(phi1, leadL1JetPhi);
      double deta = eta1 - leadL1JetEta;
      drFromLeadJet =  sqrt(dphi*dphi + deta*deta);
    }
    

    // Define Arrays for sizes of Charge, Neut Iso Radii and
    // Clustering Cone Radius.

    const int a_size = 7;
    double a_coneR[a_size];
    double a_charIsoR[a_size];
    double a_neutIsoR[a_size];

    a_coneR[0] = 17.49; // = area of 2x2
    a_coneR[1] = 26.23; // = area of 3x3
    a_coneR[2] = 30.61;
    a_coneR[3] = 34.98; // = area of 4x4
    a_coneR[4] = 39.35;
    a_coneR[5] = 43.72; // = area of 5x5
    a_coneR[6] = 52.46; // = area of 6x6
    
    for (int i=0; i<a_size; i++){
      a_charIsoR[i] = a_coneR[i]+28.9; // 28.9 gives 55.1 for 3x3 benchmark 
      a_neutIsoR[i] = a_charIsoR[i]*0.726; // Ecal radius = 0.726*Hcal radius
    }
    
    // Do Neutral Iso in radius on Ecal surface.
  
    // NxN cluster
    double e3x3=-999.0;
    double trkEcalEne =-999.0;
    edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
    iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);

    if(std::abs(point1.eta())<1.479) {
      const DetId isoCell = gEB->getClosestCell(point1);
      e3x3   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),1,1).first;  
      trkEcalEne   = spr::eCaloSimInfo(iEvent, geo, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate);
    } else {
      const DetId isoCell = gEE->getClosestCell(point1);
      e3x3   = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),1,1).first;
      trkEcalEne   = spr::eCaloSimInfo(iEvent, geo, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate);
    }

    // Cone cluster

    // Set up array of cone sizes for MIP cut
    const int a_mip_size = 5;
    double a_mipR[a_mip_size];
    a_mipR[0] = 3.84; // = area of 3x3 ecal
    a_mipR[1] = 14.0;
    a_mipR[2] = 19.0;
    a_mipR[3] = 24.0;
    a_mipR[4] = 9.0;  // = from standard analyzer

    std::vector<double> v_eDR;
    for (int i = 0 ; i < a_size ; i++){
      int nRH_eDR = 0;

      // Cone in ecal
      double eDR = spr::eCone_ecal(geo, 
                   barrelRecHitsHandle, 
                   endcapRecHitsHandle, 
                   hpoint1, point1, 
                   a_neutIsoR[i],  
                   trackMomAtEcal, nRH_eDR);
      v_eDR.push_back(eDR);
      
    }

    std::vector<double> v_eMipDR;
    for (int i = 0 ; i < a_mip_size ; i++){
      int nRH_eMipDR = 0;
      double eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, 
                      endcapRecHitsHandle, 
                      hpoint1, point1, 
                      a_mipR[i], trackMomAtEcal, nRH_eMipDR);
      
      v_eMipDR.push_back(eMipDR);
    }
    
        
    // Do charge isolation in radius at Hcal surface for 5 different
    // radii defined above in a_charIso
    
    std::vector<double> v_hmaxNearP_goodTrk;
    std::vector<double> v_hmaxNearP        ;
    std::vector<int>    v_hnNearTRKs       ;
    std::vector<int>    v_hnLayers_maxNearP;
    std::vector<int>    v_htrkQual_maxNearP;

    std::vector<double> v_cone_hmaxNearP_goodTrk;
    std::vector<double> v_cone_hmaxNearP        ;
    std::vector<int>    v_cone_hnNearTRKs       ;
    std::vector<int>    v_cone_hnLayers_maxNearP;
    std::vector<int>    v_cone_htrkQual_maxNearP;

    for (int i = 0 ; i < a_size ; i++){

      double hmaxNearP         = -999.0;
      int    hnNearTRKs        =  0;
      int    hnLayers_maxNearP =  0;
      int    htrkQual_maxNearP = -1; 
      double hmaxNearP_goodTrk = -999.0;

      double conehmaxNearP         = -999.0;
      int    conehnNearTRKs        =  0;
      int    conehnLayers_maxNearP =  0;
      int    conehtrkQual_maxNearP = -1; 
      double conehmaxNearP_goodTrk = -999.0;

      conehmaxNearP = spr::coneChargeIsolation(iEvent, iSetup, 
                           trkItr, trkCollection, 
                           *trackAssociator_, parameters_, 
                           theTrackQuality, conehnNearTRKs,
                           conehnLayers_maxNearP,
                           conehtrkQual_maxNearP, 
                           conehmaxNearP_goodTrk, 
                           hpoint1, trackMomAtHcal, 
                           a_charIsoR[i]); 

      v_hmaxNearP_goodTrk.push_back(hmaxNearP_goodTrk);
      v_hmaxNearP        .push_back(hmaxNearP        );
      v_hnNearTRKs       .push_back(hnNearTRKs       );
      v_hnLayers_maxNearP.push_back(hnLayers_maxNearP);
      v_htrkQual_maxNearP.push_back(htrkQual_maxNearP);

      v_cone_hmaxNearP_goodTrk.push_back(conehmaxNearP_goodTrk);
      v_cone_hmaxNearP        .push_back(conehmaxNearP        );
      v_cone_hnNearTRKs       .push_back(conehnNearTRKs       );
      v_cone_hnLayers_maxNearP.push_back(conehnLayers_maxNearP);
      v_cone_htrkQual_maxNearP.push_back(conehtrkQual_maxNearP);
      
    }
    
    double h3x3=-999.0, h5x5=-999.0;
    double hsim3x3=-999.0, hsim5x5=-999.0, trkHcalEne=-999.0;
    std::map<std::string, double> hsimInfo3x3, hsimInfo5x5;
    double distFromHotCell_h3x3 = -99.;
    int ietaFromHotCell_h3x3    = -99;
    int iphiFromHotCell_h3x3    = -99;
    double distFromHotCell_h5x5 = -99.;
    int ietaFromHotCell_h5x5    = -99;
    int iphiFromHotCell_h5x5    = -99;

    GlobalPoint gPosHotCell_h3x3(0.,0.,0.);
    GlobalPoint gPosHotCell_h5x5(0.,0.,0.);

    int nRH_h3x3(0), nRH_h5x5(0);

    // Hcal Energy Clustering
    
    // Get closetcell for ietaFromHotCell and iphiFromHotCell
    DetId ClosestCell;
    if( std::abs(pTrack->eta())<1.392 ) {
      ClosestCell = gHB->getClosestCell(hpoint1);
    } else {
      ClosestCell = gHE->getClosestCell(hpoint1);
    }
    // Transform into HcalDetId so that I can get ieta, iphi later.
    HcalDetId ClosestCell_HcalDetId(ClosestCell.rawId());

    // Using NxN Cluster
    std::vector<int>    v_RH_h3x3_ieta;
    std::vector<int>    v_RH_h3x3_iphi;
    std::vector<double> v_RH_h3x3_ene;
    std::vector<int>    v_RH_h5x5_ieta;
    std::vector<int>    v_RH_h5x5_iphi;
    std::vector<double> v_RH_h5x5_ene;
    
    
    h3x3 = spr::eHCALmatrix(geo, theHBHETopology, ClosestCell, hbhe,1,1,
                nRH_h3x3, v_RH_h3x3_ieta, v_RH_h3x3_iphi, v_RH_h3x3_ene, 
                gPosHotCell_h3x3);
    distFromHotCell_h3x3 = spr::getDistInPlaneTrackDir(hpoint1, trackMomAtHcal, gPosHotCell_h3x3);
    
    h5x5 = spr::eHCALmatrix(geo, theHBHETopology, ClosestCell, hbhe,2,2,
                nRH_h5x5, v_RH_h5x5_ieta, v_RH_h5x5_iphi,  v_RH_h5x5_ene, 
                gPosHotCell_h5x5);
    distFromHotCell_h5x5 = spr::getDistInPlaneTrackDir(hpoint1, trackMomAtHcal, gPosHotCell_h5x5);
    

    //     double heta = (double)hpoint1.eta();
    //     double hphi = (double)hpoint1.phi();
    std::vector<int> multiplicity3x3;
    std::vector<int> multiplicity5x5;
    if (doMC) {
      hsim3x3   = spr::eHCALmatrix(theHBHETopology, ClosestCell, 
                   pcalohh,1,1);
      hsim5x5   = spr::eHCALmatrix(theHBHETopology, ClosestCell,
                   pcalohh,2,2);

      hsimInfo3x3 = spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, pcalohh, SimTk, SimVtx, pTrack, *associate, 1,1, multiplicity3x3);
      hsimInfo5x5 = spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, pcalohh, SimTk, SimVtx, pTrack, *associate, 2,2, multiplicity5x5);
    
      // Get energy from all simhits in hcal associated with iso track
      trkHcalEne   = spr::eCaloSimInfo(iEvent, geo, pcalohh, SimTk, SimVtx, pTrack, *associate);
    }

    // Finally for cones of varying radii.
    std::vector<double> v_hsimInfoConeMatched;
    std::vector<double> v_hsimInfoConeRest   ;
    std::vector<double> v_hsimInfoConePhoton  ;
    std::vector<double> v_hsimInfoConeNeutHad;
    std::vector<double> v_hsimInfoConeCharHad;
    std::vector<double> v_hsimInfoConePdgMatched;
    std::vector<double> v_hsimInfoConeTotal  ;

    std::vector<int> v_hsimInfoConeNMatched;
    std::vector<int> v_hsimInfoConeNTotal  ;
    std::vector<int> v_hsimInfoConeNNeutHad;
    std::vector<int> v_hsimInfoConeNCharHad;
    std::vector<int> v_hsimInfoConeNPhoton ;
    std::vector<int> v_hsimInfoConeNRest   ;

    std::vector<double> v_hsimCone           ;
    std::vector<double> v_hCone              ;
    
    std::vector<int>    v_nRecHitsCone       ;
    std::vector<int>    v_nSimHitsCone       ;

    std::vector<double> v_distFromHotCell;
    std::vector<int> v_ietaFromHotCell;
    std::vector<int> v_iphiFromHotCell;
    GlobalPoint gposHotCell(0.,0.,0.);
    
    
    std::vector<int>    v_RH_r26_ieta;
    std::vector<int>    v_RH_r26_iphi;
    std::vector<double> v_RH_r26_ene;
    std::vector<int>    v_RH_r44_ieta;
    std::vector<int>    v_RH_r44_iphi;
    std::vector<double> v_RH_r44_ene;
    
    
    
    for (int i = 0 ; i < a_size ; i++){

      
      std::map<std::string, double> hsimInfoCone;
      double hsimCone = -999.0, hCone = -999.0;
      double distFromHotCell = -99.0;
      int ietaFromHotCell = -99;
      int iphiFromHotCell = -99;
      int ietaHotCell = -99;
      int iphiHotCell = -99;
      int nRecHitsCone = -999;
      int nSimHitsCone = -999;

      std::vector<int> multiplicityCone;
      std::vector<DetId> coneRecHitDetIds;
      if (doMC) 
    hsimCone = spr::eCone_hcal(geo, pcalohh, hpoint1, point1, 
                   a_coneR[i], trackMomAtHcal, nSimHitsCone);
      
      // If needed, get ieta and iphi of rechits for cones of 23.25
      // and for hitmap for debugging
      bool makeHitmaps = false;
      if (a_coneR[i] == 26.23 && makeHitmaps)
      {
    
    hCone = spr::eCone_hcal(geo, hbhe, hpoint1, point1, 
                a_coneR[i], trackMomAtHcal,nRecHitsCone,
                v_RH_r26_ieta, v_RH_r26_iphi,  v_RH_r26_ene, 
                coneRecHitDetIds, distFromHotCell, 
                ietaHotCell, iphiHotCell, gposHotCell);
      } 
      else if (a_coneR[i] == 43.72 && makeHitmaps)
      {
    
    hCone = spr::eCone_hcal(geo, hbhe, hpoint1, point1, 
                a_coneR[i], trackMomAtHcal,nRecHitsCone,
                v_RH_r44_ieta, v_RH_r44_iphi,  v_RH_r44_ene, 
                coneRecHitDetIds, distFromHotCell, 
                ietaHotCell, iphiHotCell, gposHotCell);
      } 
      else 
      {
    
    hCone = spr::eCone_hcal(geo, hbhe, hpoint1, point1, 
                a_coneR[i], trackMomAtHcal, nRecHitsCone, 
                coneRecHitDetIds, distFromHotCell, 
                ietaHotCell, iphiHotCell, gposHotCell);
      }

      
      
      if (ietaHotCell != 99){
    ietaFromHotCell = ietaHotCell-ClosestCell_HcalDetId.ieta();
    iphiFromHotCell = iphiHotCell-ClosestCell_HcalDetId.iphi();
      }
      
      // SimHits NOT matched to RecHits
      if (doMC) {
    hsimInfoCone = spr::eHCALSimInfoCone(iEvent,pcalohh, SimTk, SimVtx, pTrack, *associate, geo, hpoint1, point1, a_coneR[i], trackMomAtHcal, multiplicityCone);      
      
      
      
    // SimHits matched to RecHits
    //       hsimInfoCone = spr::eHCALSimInfoCone(iEvent,pcalohh, SimTk, SimVtx, 
    //                             pTrack, *associate, 
    //                             geo, hpoint1, point1, 
    //                             a_coneR[i], trackMomAtHcal, 
    //                             multiplicityCone,
    //                             coneRecHitDetIds);      
      
    v_hsimInfoConeMatched   .push_back(hsimInfoCone["eMatched"   ]);
    v_hsimInfoConeRest      .push_back(hsimInfoCone["eRest"      ]);
    v_hsimInfoConePhoton    .push_back(hsimInfoCone["eGamma"     ]);
    v_hsimInfoConeNeutHad   .push_back(hsimInfoCone["eNeutralHad"]);
    v_hsimInfoConeCharHad   .push_back(hsimInfoCone["eChargedHad"]);
    v_hsimInfoConePdgMatched.push_back(hsimInfoCone["pdgMatched" ]);
    v_hsimInfoConeTotal     .push_back(hsimInfoCone["eTotal"     ]);

    v_hsimInfoConeNMatched  .push_back(multiplicityCone.at(0));

    v_hsimInfoConeNTotal      .push_back(multiplicityCone.at(1));
    v_hsimInfoConeNNeutHad    .push_back(multiplicityCone.at(2));
    v_hsimInfoConeNCharHad    .push_back(multiplicityCone.at(3));
    v_hsimInfoConeNPhoton     .push_back(multiplicityCone.at(4));
    v_hsimInfoConeNRest       .push_back(multiplicityCone.at(5));

    v_hsimCone                .push_back(hsimCone                   );
    v_nSimHitsCone            .push_back(nSimHitsCone               );
      }
      v_hCone                   .push_back(hCone                      );
      v_nRecHitsCone            .push_back(nRecHitsCone               );
      
      v_distFromHotCell         .push_back(distFromHotCell            );
      v_ietaFromHotCell         .push_back(ietaFromHotCell            );
      v_iphiFromHotCell         .push_back(iphiFromHotCell            );
      
      
    }
    
 
    // Fill Vectors that go into root file

    t_v_hnNearTRKs        ->push_back(v_hnNearTRKs                 );
    t_v_hnLayers_maxNearP ->push_back(v_hnLayers_maxNearP          );
    t_v_htrkQual_maxNearP ->push_back(v_htrkQual_maxNearP          );
    t_v_hmaxNearP_goodTrk ->push_back(v_hmaxNearP_goodTrk          ); 
    t_v_hmaxNearP         ->push_back(v_hmaxNearP                  );    

    t_v_cone_hnNearTRKs        ->push_back(v_cone_hnNearTRKs       );
    t_v_cone_hnLayers_maxNearP ->push_back(v_cone_hnLayers_maxNearP);
    t_v_cone_htrkQual_maxNearP ->push_back(v_cone_htrkQual_maxNearP);
    t_v_cone_hmaxNearP_goodTrk ->push_back(v_cone_hmaxNearP_goodTrk); 
    t_v_cone_hmaxNearP         ->push_back(v_cone_hmaxNearP        );    

    //    t_hScale            ->push_back(hScale                     );
    t_trkNOuterHits     ->push_back(nOuterHits                 );
    t_trkNLayersCrossed ->push_back(nLayersCrossed             );
    t_dtFromLeadJet     ->push_back(drFromLeadJet              );
    t_trkP              ->push_back(p1                         );
    t_trkPt             ->push_back(pt1                        );
    t_trkEta            ->push_back(eta1                       );
    t_trkPhi            ->push_back(phi1                       );

    t_e3x3              ->push_back(e3x3                       );
    t_v_eDR             ->push_back(v_eDR                      );
    t_v_eMipDR          ->push_back(v_eMipDR                   );

    t_h3x3              ->push_back(h3x3                       );
    t_h5x5              ->push_back(h5x5                       );
    t_nRH_h3x3          ->push_back(nRH_h3x3                   );
    t_nRH_h5x5          ->push_back(nRH_h5x5                   );
    
    t_v_RH_h3x3_ieta    ->push_back(v_RH_h3x3_ieta);
    t_v_RH_h3x3_iphi    ->push_back(v_RH_h3x3_iphi);
    t_v_RH_h3x3_ene     ->push_back(v_RH_h3x3_ene);
    t_v_RH_h5x5_ieta    ->push_back(v_RH_h5x5_ieta);
    t_v_RH_h5x5_iphi    ->push_back(v_RH_h5x5_iphi);
    t_v_RH_h5x5_ene     ->push_back(v_RH_h5x5_ene);

    if (doMC) {
      t_simP              ->push_back(simP                       );
      t_hsim3x3           ->push_back(hsim3x3                    );
      t_hsim5x5           ->push_back(hsim5x5                    );

      t_hsim3x3Matched    ->push_back(hsimInfo3x3["eMatched"]    );
      t_hsim5x5Matched    ->push_back(hsimInfo5x5["eMatched"]    );
      t_hsim3x3Rest       ->push_back(hsimInfo3x3["eRest"]       );
      t_hsim5x5Rest       ->push_back(hsimInfo5x5["eRest"]       );
      t_hsim3x3Photon     ->push_back(hsimInfo3x3["eGamma"]      );
      t_hsim5x5Photon     ->push_back(hsimInfo5x5["eGamma"]      );
      t_hsim3x3NeutHad    ->push_back(hsimInfo3x3["eNeutralHad"] );
      t_hsim5x5NeutHad    ->push_back(hsimInfo5x5["eNeutralHad"] );
      t_hsim3x3CharHad    ->push_back(hsimInfo3x3["eChargedHad"] );
      t_hsim5x5CharHad    ->push_back(hsimInfo5x5["eChargedHad"] );
      t_hsim3x3Total      ->push_back(hsimInfo3x3["eTotal"]      );
      t_hsim5x5Total      ->push_back(hsimInfo5x5["eTotal"]      );
      t_hsim3x3PdgMatched ->push_back(hsimInfo3x3["pdgMatched"]  );
      t_hsim5x5PdgMatched ->push_back(hsimInfo5x5["pdgMatched"]  );

      t_hsim3x3NMatched   ->push_back(multiplicity3x3.at(0));
      t_hsim3x3NTotal     ->push_back(multiplicity3x3.at(1));
      t_hsim3x3NNeutHad   ->push_back(multiplicity3x3.at(2));
      t_hsim3x3NCharHad   ->push_back(multiplicity3x3.at(3));
      t_hsim3x3NPhoton    ->push_back(multiplicity3x3.at(4));
      t_hsim3x3NRest      ->push_back(multiplicity3x3.at(5));

      t_hsim5x5NMatched   ->push_back(multiplicity5x5.at(0));
      t_hsim5x5NTotal     ->push_back(multiplicity5x5.at(1));
      t_hsim5x5NNeutHad   ->push_back(multiplicity5x5.at(2));
      t_hsim5x5NCharHad   ->push_back(multiplicity5x5.at(3));
      t_hsim5x5NPhoton    ->push_back(multiplicity5x5.at(4));
      t_hsim5x5NRest      ->push_back(multiplicity5x5.at(5));
    }
    
    t_distFromHotCell_h3x3->push_back(distFromHotCell_h3x3);
    t_ietaFromHotCell_h3x3->push_back(ietaFromHotCell_h3x3);
    t_iphiFromHotCell_h3x3->push_back(iphiFromHotCell_h3x3);
    t_distFromHotCell_h5x5->push_back(distFromHotCell_h5x5);
    t_ietaFromHotCell_h5x5->push_back(ietaFromHotCell_h5x5);
    t_iphiFromHotCell_h5x5->push_back(iphiFromHotCell_h5x5);
    
    if (doMC) {
      t_trkHcalEne                ->push_back(trkHcalEne                 );
      t_trkEcalEne                ->push_back(trkEcalEne                 );

      t_v_hsimInfoConeMatched     ->push_back(v_hsimInfoConeMatched   );
      t_v_hsimInfoConeRest        ->push_back(v_hsimInfoConeRest      );
      t_v_hsimInfoConePhoton      ->push_back(v_hsimInfoConePhoton   );
      t_v_hsimInfoConeNeutHad     ->push_back(v_hsimInfoConeNeutHad   );
      t_v_hsimInfoConeCharHad     ->push_back(v_hsimInfoConeCharHad   );
      t_v_hsimInfoConePdgMatched  ->push_back(v_hsimInfoConePdgMatched);  
      t_v_hsimInfoConeTotal       ->push_back(v_hsimInfoConeTotal     );

      t_v_hsimInfoConeNMatched    ->push_back(v_hsimInfoConeNMatched  );
      t_v_hsimInfoConeNTotal      ->push_back(v_hsimInfoConeNTotal    );
      t_v_hsimInfoConeNNeutHad    ->push_back(v_hsimInfoConeNNeutHad  );
      t_v_hsimInfoConeNCharHad    ->push_back(v_hsimInfoConeNCharHad  );
      t_v_hsimInfoConeNPhoton     ->push_back(v_hsimInfoConeNPhoton    );
      t_v_hsimInfoConeNRest       ->push_back(v_hsimInfoConeNRest     );  

      t_v_hsimCone    ->push_back(v_hsimCone    );
      t_v_hCone       ->push_back(v_hCone       );
      t_v_nRecHitsCone->push_back(v_nRecHitsCone);
      t_v_nSimHitsCone->push_back(v_nSimHitsCone);
    }

    
    t_v_distFromHotCell->push_back(v_distFromHotCell);
    t_v_ietaFromHotCell->push_back(v_ietaFromHotCell);
    t_v_iphiFromHotCell->push_back(v_iphiFromHotCell);
    
    t_v_RH_r26_ieta    ->push_back(v_RH_r26_ieta);
    t_v_RH_r26_iphi    ->push_back(v_RH_r26_iphi);
    t_v_RH_r26_ene     ->push_back(v_RH_r26_ene);
    t_v_RH_r44_ieta    ->push_back(v_RH_r44_ieta);
    t_v_RH_r44_iphi    ->push_back(v_RH_r44_iphi);
    t_v_RH_r44_ene     ->push_back(v_RH_r44_ene);

    

    t_v_hlTriggers    ->push_back(v_hlTriggers);
    t_hltHB ->push_back(hltHB);
    t_hltHE ->push_back(hltHE);
    t_hltL1Jet15                 ->push_back(hltL1Jet15             );
    t_hltJet30               ->push_back(hltJet30               );
    t_hltJet50               ->push_back(hltJet50               );
    t_hltJet80               ->push_back(hltJet80               );
    t_hltJet110              ->push_back(hltJet110              );
    t_hltJet140              ->push_back(hltJet140              );
    t_hltJet180              ->push_back(hltJet180              );
    t_hltL1SingleEG5         ->push_back(hltL1SingleEG5     );
    t_hltZeroBias        ->push_back(hltZeroBias        );
    t_hltMinBiasHcal         ->push_back(hltMinBiasHcal     );
    t_hltMinBiasEcal         ->push_back(hltMinBiasEcal     );
    t_hltMinBiasPixel            ->push_back(hltMinBiasPixel            );
    t_hltSingleIsoTau30_Trk5     ->push_back(hltSingleIsoTau30_Trk5 );
    t_hltDoubleLooseIsoTau15_Trk5->push_back(hltDoubleLooseIsoTau15_Trk5);

    t_irun->push_back(irun);
    t_ievt->push_back(ievt);
    t_ilum->push_back(ilum);
    
    nTRK++;
    
    
  } // Loop over track collection
  
    //  std::cout << "nEVT= " << nEVT << std::endl;
  
  ntp->Fill();
  nEVT++;
  
}

void IsolatedTracksCone::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 924 of file IsolatedTracksCone.cc.

References BuildTree(), doMC, genPartEtaBins, genPartPBins, mps_fire::i, nEVT, nEVT_failL1, nTRK, t_distFromHotCell_h3x3, t_distFromHotCell_h5x5, t_dtFromLeadJet, t_e3x3, t_h3x3, t_h5x5, t_hltDoubleLooseIsoTau15_Trk5, t_hltHB, t_hltHE, t_hltJet110, t_hltJet140, t_hltJet180, t_hltJet30, t_hltJet50, t_hltJet80, t_hltL1Jet15, t_hltL1SingleEG5, t_hltMinBiasEcal, t_hltMinBiasHcal, t_hltMinBiasPixel, t_hltSingleIsoTau30_Trk5, t_hltZeroBias, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NCharHad, t_hsim3x3NeutHad, t_hsim3x3NMatched, t_hsim3x3NNeutHad, t_hsim3x3NPhoton, t_hsim3x3NRest, t_hsim3x3NTotal, t_hsim3x3PdgMatched, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim3x3Total, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NCharHad, t_hsim5x5NeutHad, t_hsim5x5NMatched, t_hsim5x5NNeutHad, t_hsim5x5NPhoton, t_hsim5x5NRest, t_hsim5x5NTotal, t_hsim5x5PdgMatched, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim5x5Total, t_ietaFromHotCell_h3x3, t_ietaFromHotCell_h5x5, t_ievt, t_ilum, t_iphiFromHotCell_h3x3, t_iphiFromHotCell_h5x5, t_irun, t_nRH_h3x3, t_nRH_h5x5, t_simP, t_trkEcalEne, t_trkEta, t_trkHcalEne, t_trkNLayersCrossed, t_trkNOuterHits, t_trkP, t_trkPhi, t_trkPt, t_v_cone_hmaxNearP, t_v_cone_hmaxNearP_goodTrk, t_v_cone_hnLayers_maxNearP, t_v_cone_hnNearTRKs, t_v_cone_htrkQual_maxNearP, t_v_distFromHotCell, t_v_eDR, t_v_eMipDR, t_v_hCone, t_v_hlTriggers, t_v_hmaxNearP, t_v_hmaxNearP_goodTrk, t_v_hnLayers_maxNearP, t_v_hnNearTRKs, t_v_hsimCone, t_v_hsimInfoConeCharHad, t_v_hsimInfoConeMatched, t_v_hsimInfoConeNCharHad, t_v_hsimInfoConeNeutHad, t_v_hsimInfoConeNMatched, t_v_hsimInfoConeNNeutHad, t_v_hsimInfoConeNPhoton, t_v_hsimInfoConeNRest, t_v_hsimInfoConeNTotal, t_v_hsimInfoConePdgMatched, t_v_hsimInfoConePhoton, t_v_hsimInfoConeRest, t_v_hsimInfoConeTotal, t_v_htrkQual_maxNearP, t_v_ietaFromHotCell, t_v_iphiFromHotCell, t_v_nRecHitsCone, t_v_nSimHitsCone, t_v_RH_h3x3_ene, t_v_RH_h3x3_ieta, t_v_RH_h3x3_iphi, t_v_RH_h5x5_ene, t_v_RH_h5x5_ieta, t_v_RH_h5x5_iphi, t_v_RH_r26_ene, t_v_RH_r26_ieta, t_v_RH_r26_iphi, t_v_RH_r44_ene, t_v_RH_r44_ieta, and t_v_RH_r44_iphi.

                                  {

  //   hbScale = 120.0;
  //   heScale = 135.0;
  nEVT=0;
  nEVT_failL1=0;
  nTRK=0;
  
  double tempgen_TH[22] = { 0.0,  1.0,  2.0,  3.0,  4.0,  
                5.0,  6.0,  7.0,  8.0,  9.0, 
                10.0, 12.0, 15.0, 20.0, 25.0, 
                30.0, 40.0, 50.0, 60.0, 70.0, 80.0, 100};
  for(int i=0; i<22; i++)  genPartPBins[i]  = tempgen_TH[i];


  double tempgen_Eta[5] = {0.0, 0.5, 1.1, 1.7, 2.0};
  for(int i=0; i<5; i++) genPartEtaBins[i] = tempgen_Eta[i];

  t_v_hnNearTRKs           = new std::vector<std::vector<int> >   (); 
  t_v_hnLayers_maxNearP    = new std::vector<std::vector<int> >   (); 
  t_v_htrkQual_maxNearP    = new std::vector<std::vector<int> >  (); 
  t_v_hmaxNearP_goodTrk    = new std::vector<std::vector<double> >();
  t_v_hmaxNearP            = new std::vector<std::vector<double> >();
                                    
  t_v_cone_hnNearTRKs           = new std::vector<std::vector<int> >   (); 
  t_v_cone_hnLayers_maxNearP    = new std::vector<std::vector<int> >   (); 
  t_v_cone_htrkQual_maxNearP    = new std::vector<std::vector<int> >  (); 
  t_v_cone_hmaxNearP_goodTrk    = new std::vector<std::vector<double> >();
  t_v_cone_hmaxNearP            = new std::vector<std::vector<double> >();
                                    
  //  t_hScale             = new std::vector<double>();
  t_trkNOuterHits      = new std::vector<double>();
  t_trkNLayersCrossed  = new std::vector<double>();
  t_dtFromLeadJet      = new std::vector<double>();
  t_trkP               = new std::vector<double>();
  t_trkPt              = new std::vector<double>();
  t_trkEta             = new std::vector<double>();
  t_trkPhi             = new std::vector<double>();

  t_e3x3               = new std::vector<double>();
  t_v_eDR              = new std::vector<std::vector<double> >();
  t_v_eMipDR           = new std::vector<std::vector<double> >();

  t_h3x3               = new std::vector<double>();
  t_h5x5               = new std::vector<double>();

  t_nRH_h3x3           = new std::vector<double>();
  t_nRH_h5x5           = new std::vector<double>();

  if (doMC) {
    t_simP               = new std::vector<double>();
    t_hsim3x3            = new std::vector<double>();
    t_hsim5x5            = new std::vector<double>();

    t_hsim3x3Matched     = new std::vector<double>();
    t_hsim5x5Matched     = new std::vector<double>();
    t_hsim3x3Rest        = new std::vector<double>();
    t_hsim5x5Rest        = new std::vector<double>();
    t_hsim3x3Photon      = new std::vector<double>();
    t_hsim5x5Photon      = new std::vector<double>();
    t_hsim3x3NeutHad     = new std::vector<double>();
    t_hsim5x5NeutHad     = new std::vector<double>();
    t_hsim3x3CharHad     = new std::vector<double>();
    t_hsim5x5CharHad     = new std::vector<double>();
    t_hsim3x3PdgMatched  = new std::vector<double>();
    t_hsim5x5PdgMatched  = new std::vector<double>();
    t_hsim3x3Total       = new std::vector<double>();
    t_hsim5x5Total       = new std::vector<double>();

    t_hsim3x3NMatched  = new std::vector<int>();
    t_hsim3x3NTotal    = new std::vector<int>();
    t_hsim3x3NNeutHad  = new std::vector<int>();
    t_hsim3x3NCharHad  = new std::vector<int>();
    t_hsim3x3NPhoton   = new std::vector<int>();
    t_hsim3x3NRest     = new std::vector<int>();

    t_hsim5x5NMatched  = new std::vector<int>();
    t_hsim5x5NTotal    = new std::vector<int>();
    t_hsim5x5NNeutHad  = new std::vector<int>();
    t_hsim5x5NCharHad  = new std::vector<int>();
    t_hsim5x5NPhoton   = new std::vector<int>();
    t_hsim5x5NRest     = new std::vector<int>();

    t_trkHcalEne         = new std::vector<double>();
    t_trkEcalEne         = new std::vector<double>();
  }

  t_distFromHotCell_h3x3  = new std::vector<double>();
  t_ietaFromHotCell_h3x3  = new std::vector<int>();
  t_iphiFromHotCell_h3x3  = new std::vector<int>();
  t_distFromHotCell_h5x5  = new std::vector<double>();
  t_ietaFromHotCell_h5x5  = new std::vector<int>();
  t_iphiFromHotCell_h5x5  = new std::vector<int>();

  if (doMC) {
    t_v_hsimInfoConeMatched   = new std::vector<std::vector<double> >();
    t_v_hsimInfoConeRest      = new std::vector<std::vector<double> >();
    t_v_hsimInfoConePhoton     = new std::vector<std::vector<double> >();
    t_v_hsimInfoConeNeutHad   = new std::vector<std::vector<double> >();
    t_v_hsimInfoConeCharHad   = new std::vector<std::vector<double> >();
    t_v_hsimInfoConePdgMatched= new std::vector<std::vector<double> >();
    t_v_hsimInfoConeTotal     = new std::vector<std::vector<double> >();

    t_v_hsimInfoConeNMatched  = new std::vector<std::vector<int> >();
    t_v_hsimInfoConeNTotal    = new std::vector<std::vector<int> >();
    t_v_hsimInfoConeNNeutHad  = new std::vector<std::vector<int> >();
    t_v_hsimInfoConeNCharHad  = new std::vector<std::vector<int> >();
    t_v_hsimInfoConeNPhoton    = new std::vector<std::vector<int> >();
    t_v_hsimInfoConeNRest     = new std::vector<std::vector<int> >();

    t_v_hsimCone    = new std::vector<std::vector<double> >();
  }

  t_v_hCone       = new std::vector<std::vector<double> >();
  t_v_nRecHitsCone= new std::vector<std::vector<int> >();
  t_v_nSimHitsCone= new std::vector<std::vector<int> >();

  t_v_distFromHotCell= new std::vector<std::vector<double> >();
  t_v_ietaFromHotCell= new std::vector<std::vector<int> >();
  t_v_iphiFromHotCell= new std::vector<std::vector<int> >();

  t_v_RH_h3x3_ieta = new std::vector<std::vector<int> >();
  t_v_RH_h3x3_iphi = new std::vector<std::vector<int> >();
  t_v_RH_h3x3_ene  = new std::vector<std::vector<double> >();
  t_v_RH_h5x5_ieta = new std::vector<std::vector<int> >();
  t_v_RH_h5x5_iphi = new std::vector<std::vector<int> >();
  t_v_RH_h5x5_ene  = new std::vector<std::vector<double> >();
  t_v_RH_r26_ieta  = new std::vector<std::vector<int> >();
  t_v_RH_r26_iphi  = new std::vector<std::vector<int> >();
  t_v_RH_r26_ene   = new std::vector<std::vector<double> >();
  t_v_RH_r44_ieta  = new std::vector<std::vector<int> >();
  t_v_RH_r44_iphi  = new std::vector<std::vector<int> >();
  t_v_RH_r44_ene   = new std::vector<std::vector<double> >();


  t_v_hlTriggers    = new std::vector<std::vector<int> >();

  t_hltHE                         = new std::vector<int>();
  t_hltHB                         = new std::vector<int>();
  t_hltL1Jet15                = new std::vector<int>();
  t_hltJet30                  = new std::vector<int>();
  t_hltJet50                  = new std::vector<int>();
  t_hltJet80                  = new std::vector<int>();
  t_hltJet110                 = new std::vector<int>();
  t_hltJet140                 = new std::vector<int>();
  t_hltJet180                 = new std::vector<int>();
  t_hltL1SingleEG5            = new std::vector<int>();
  t_hltZeroBias               = new std::vector<int>();
  t_hltMinBiasHcal            = new std::vector<int>();
  t_hltMinBiasEcal            = new std::vector<int>();
  t_hltMinBiasPixel           = new std::vector<int>();
  t_hltSingleIsoTau30_Trk5        = new std::vector<int>();
  t_hltDoubleLooseIsoTau15_Trk5   = new std::vector<int>();
  

  t_irun = new std::vector<unsigned int>();
  t_ievt = new std::vector<unsigned int>();
  t_ilum = new std::vector<unsigned int>();

  BuildTree();
}

void IsolatedTracksCone::BuildTree ( )

private

Definition at line 1235 of file IsolatedTracksCone.cc.

References doMC, fs, hRawEta, hRawP, hRawPhi, hRawPt, leadL1JetEta, leadL1JetPhi, leadL1JetPT, TFileService::make(), nEVT, nFailEta, nFailHighPurityQaul, nFailPt, nMissEcal, nMissHcal, nRawTRK, ntp, nTRK, t_distFromHotCell_h3x3, t_distFromHotCell_h5x5, t_dtFromLeadJet, t_e3x3, t_h3x3, t_h5x5, t_hltDoubleLooseIsoTau15_Trk5, t_hltHB, t_hltHE, t_hltJet110, t_hltJet140, t_hltJet180, t_hltJet30, t_hltJet50, t_hltJet80, t_hltL1Jet15, t_hltL1SingleEG5, t_hltMinBiasEcal, t_hltMinBiasHcal, t_hltMinBiasPixel, t_hltSingleIsoTau30_Trk5, t_hltZeroBias, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NCharHad, t_hsim3x3NeutHad, t_hsim3x3NMatched, t_hsim3x3NNeutHad, t_hsim3x3NPhoton, t_hsim3x3NRest, t_hsim3x3NTotal, t_hsim3x3PdgMatched, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim3x3Total, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NCharHad, t_hsim5x5NeutHad, t_hsim5x5NMatched, t_hsim5x5NNeutHad, t_hsim5x5NPhoton, t_hsim5x5NRest, t_hsim5x5NTotal, t_hsim5x5PdgMatched, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim5x5Total, t_ietaFromHotCell_h3x3, t_ietaFromHotCell_h5x5, t_ievt, t_ilum, t_iphiFromHotCell_h3x3, t_iphiFromHotCell_h5x5, t_irun, t_nRH_h3x3, t_nRH_h5x5, t_simP, t_trkEcalEne, t_trkEta, t_trkHcalEne, t_trkNLayersCrossed, t_trkNOuterHits, t_trkP, t_trkPhi, t_trkPt, t_v_cone_hmaxNearP, t_v_cone_hmaxNearP_goodTrk, t_v_cone_hnLayers_maxNearP, t_v_cone_hnNearTRKs, t_v_cone_htrkQual_maxNearP, t_v_distFromHotCell, t_v_eDR, t_v_eMipDR, t_v_hCone, t_v_hlTriggers, t_v_hmaxNearP, t_v_hmaxNearP_goodTrk, t_v_hnLayers_maxNearP, t_v_hnNearTRKs, t_v_hsimCone, t_v_hsimInfoConeCharHad, t_v_hsimInfoConeMatched, t_v_hsimInfoConeNCharHad, t_v_hsimInfoConeNeutHad, t_v_hsimInfoConeNMatched, t_v_hsimInfoConeNNeutHad, t_v_hsimInfoConeNPhoton, t_v_hsimInfoConeNRest, t_v_hsimInfoConeNTotal, t_v_hsimInfoConePdgMatched, t_v_hsimInfoConePhoton, t_v_hsimInfoConeRest, t_v_hsimInfoConeTotal, t_v_htrkQual_maxNearP, t_v_ietaFromHotCell, t_v_iphiFromHotCell, t_v_nRecHitsCone, t_v_nSimHitsCone, t_v_RH_h3x3_ene, t_v_RH_h3x3_ieta, t_v_RH_h3x3_iphi, t_v_RH_h5x5_ene, t_v_RH_h5x5_ieta, t_v_RH_h5x5_iphi, t_v_RH_r26_ene, t_v_RH_r26_ieta, t_v_RH_r26_iphi, t_v_RH_r44_ene, t_v_RH_r44_ieta, and t_v_RH_r44_iphi.

Referenced by beginJob().

                                  {


  hRawPt  = fs->make<TH1F>("hRawPt ", "hRawPt ",  100,  0.0, 100.0);
  hRawP   = fs->make<TH1F>("hRawP  ", "hRawP  ",  100,  0.0, 100.0);
  hRawEta = fs->make<TH1F>("hRawEta", "hRawEta",   15,  0.0,   3.0);
  hRawPhi = fs->make<TH1F>("hRawPhi", "hRawPhi",  100, -3.2,   3.2);

  ntp = fs->make<TTree>("ntp", "ntp");

  
  // Counters
  ntp->Branch("nEVT"        , &nEVT        , "nEVT/I"        );
  ntp->Branch("leadL1JetPT" , &leadL1JetPT , "leadL1JetPT/D" );
  ntp->Branch("leadL1JetEta", &leadL1JetEta, "leadL1JetEta/D");
  ntp->Branch("leadL1JetPhi", &leadL1JetPhi, "leadL1JetPhi/D");
  ntp->Branch("nTRK",         &nTRK,         "nTRK/I");
  ntp->Branch("nRawTRK"            , &nRawTRK            ,"nRawTRK/I"            );
  ntp->Branch("nFailHighPurityQaul", &nFailHighPurityQaul,"nFailHighPurityQaul/I");
  ntp->Branch("nFailPt"            , &nFailPt            ,"nFailPt/I"            );
  ntp->Branch("nFailEta"           , &nFailEta           ,"nFailEta/I"           );
  ntp->Branch("nMissEcal"          , &nMissEcal          ,"nMissEcal/I"          );
  ntp->Branch("nMissHcal"          , &nMissHcal          ,"nMissHcal/I"          );

  ntp->Branch("hnNearTRKs"          ,"vector<vector<int> >   ",&t_v_hnNearTRKs          );
  ntp->Branch("hnLayers_maxNearP"   ,"vector<vector<int> >   ",&t_v_hnLayers_maxNearP   );
  ntp->Branch("htrkQual_maxNearP"   ,"vector<vector<int> >   ",&t_v_htrkQual_maxNearP   );
  ntp->Branch("hmaxNearP_goodTrk"   ,"vector<vector<double> >",&t_v_hmaxNearP_goodTrk   );
  ntp->Branch("hmaxNearP"           ,"vector<vector<double> >",&t_v_hmaxNearP           );

  ntp->Branch("cone_hnNearTRKs"       ,"vector<vector<int> >   ",&t_v_cone_hnNearTRKs       );
  ntp->Branch("cone_hnLayers_maxNearP","vector<vector<int> >   ",&t_v_cone_hnLayers_maxNearP);
  ntp->Branch("cone_htrkQual_maxNearP","vector<vector<int> >   ",&t_v_cone_htrkQual_maxNearP);
  ntp->Branch("cone_hmaxNearP_goodTrk","vector<vector<double> >",&t_v_cone_hmaxNearP_goodTrk);
  ntp->Branch("cone_hmaxNearP"        ,"vector<vector<double> >",&t_v_cone_hmaxNearP        );
                                        
  //  ntp->Branch("hScale"           , "vector<double>", &t_hScale           );
  ntp->Branch("trkNOuterHits"    , "vector<double>", &t_trkNOuterHits    );
  ntp->Branch("trkNLayersCrossed", "vector<double>", &t_trkNLayersCrossed);
  ntp->Branch("drFromLeadJet"    , "vector<double>", &t_dtFromLeadJet    );
  ntp->Branch("trkP"             , "vector<double>", &t_trkP             );
  ntp->Branch("trkPt"            , "vector<double>", &t_trkPt            );
  ntp->Branch("trkEta"           , "vector<double>", &t_trkEta           );
  ntp->Branch("trkPhi"           , "vector<double>", &t_trkPhi           );
  ntp->Branch("e3x3"             , "vector<double>", &t_e3x3             );

  ntp->Branch("e3x3"          , "vector<double>"         , &t_e3x3 );
  ntp->Branch("v_eDR"         , "vector<vector<double> >", &t_v_eDR);
  ntp->Branch("v_eMipDR"      , "vector<vector<double> >", &t_v_eMipDR);

  ntp->Branch("h3x3"             , "vector<double>", &t_h3x3             );
  ntp->Branch("h5x5"             , "vector<double>", &t_h5x5             );
  ntp->Branch("nRH_h3x3"         , "vector<double>", &t_nRH_h3x3         );
  ntp->Branch("nRH_h5x5"         , "vector<double>", &t_nRH_h5x5         );

  if (doMC) {
    ntp->Branch("simP"             , "vector<double>", &t_simP             );
    ntp->Branch("hsim3x3"          , "vector<double>", &t_hsim3x3          );
    ntp->Branch("hsim5x5"          , "vector<double>", &t_hsim5x5          );

    ntp->Branch("hsim3x3Matched"   , "vector<double>", &t_hsim3x3Matched   );
    ntp->Branch("hsim5x5Matched"   , "vector<double>", &t_hsim5x5Matched   );
    ntp->Branch("hsim3x3Rest"      , "vector<double>", &t_hsim3x3Rest      );
    ntp->Branch("hsim5x5Rest"      , "vector<double>", &t_hsim5x5Rest      );
    ntp->Branch("hsim3x3Photon"    , "vector<double>", &t_hsim3x3Photon    );
    ntp->Branch("hsim5x5Photon"    , "vector<double>", &t_hsim5x5Photon    );
    ntp->Branch("hsim3x3NeutHad"   , "vector<double>", &t_hsim3x3NeutHad   );
    ntp->Branch("hsim5x5NeutHad"   , "vector<double>", &t_hsim5x5NeutHad   );
    ntp->Branch("hsim3x3CharHad"   , "vector<double>", &t_hsim3x3CharHad   );
    ntp->Branch("hsim5x5CharHad"   , "vector<double>", &t_hsim5x5CharHad   );
    ntp->Branch("hsim3x3PdgMatched", "vector<double>", &t_hsim3x3PdgMatched);
    ntp->Branch("hsim5x5PdgMatched", "vector<double>", &t_hsim5x5PdgMatched);
    ntp->Branch("hsim3x3Total"     , "vector<double>", &t_hsim3x3Total     );
    ntp->Branch("hsim5x5Total"     , "vector<double>", &t_hsim5x5Total     );

    ntp->Branch("hsim3x3NMatched"   , "vector<int>", &t_hsim3x3NMatched   );
    ntp->Branch("hsim3x3NRest"      , "vector<int>", &t_hsim3x3NRest      );
    ntp->Branch("hsim3x3NPhoton"    , "vector<int>", &t_hsim3x3NPhoton    );
    ntp->Branch("hsim3x3NNeutHad"   , "vector<int>", &t_hsim3x3NNeutHad   );
    ntp->Branch("hsim3x3NCharHad"   , "vector<int>", &t_hsim3x3NCharHad   );
    ntp->Branch("hsim3x3NTotal"     , "vector<int>", &t_hsim3x3NTotal     );

    ntp->Branch("hsim5x5NMatched"   , "vector<int>", &t_hsim5x5NMatched   );
    ntp->Branch("hsim5x5NRest"      , "vector<int>", &t_hsim5x5NRest      );
    ntp->Branch("hsim5x5NPhoton"    , "vector<int>", &t_hsim5x5NPhoton    );
    ntp->Branch("hsim5x5NNeutHad"   , "vector<int>", &t_hsim5x5NNeutHad   );
    ntp->Branch("hsim5x5NCharHad"   , "vector<int>", &t_hsim5x5NCharHad   );
    ntp->Branch("hsim5x5NTotal"     , "vector<int>", &t_hsim5x5NTotal     );

    ntp->Branch("trkHcalEne"       , "vector<double>", &t_trkHcalEne       );
    ntp->Branch("trkEcalEne"       , "vector<double>", &t_trkEcalEne       );
  }

  ntp->Branch("distFromHotCell_h3x3", "vector<double>", &t_distFromHotCell_h3x3);
  ntp->Branch("ietaFromHotCell_h3x3", "vector<int>", &t_ietaFromHotCell_h3x3);
  ntp->Branch("iphiFromHotCell_h3x3", "vector<int>", &t_iphiFromHotCell_h3x3);
  ntp->Branch("distFromHotCell_h5x5", "vector<double>", &t_distFromHotCell_h5x5);
  ntp->Branch("ietaFromHotCell_h5x5", "vector<int>", &t_ietaFromHotCell_h5x5);
  ntp->Branch("iphiFromHotCell_h5x5", "vector<int>", &t_iphiFromHotCell_h5x5);

  if (doMC) {
    ntp->Branch("v_hsimInfoConeMatched"   ,"vector<vector<double> >",&t_v_hsimInfoConeMatched   );
    ntp->Branch("v_hsimInfoConeRest"      ,"vector<vector<double> >",&t_v_hsimInfoConeRest      );
    ntp->Branch("v_hsimInfoConePhoton"     ,"vector<vector<double> >",&t_v_hsimInfoConePhoton     );
    ntp->Branch("v_hsimInfoConeNeutHad"   ,"vector<vector<double> >",&t_v_hsimInfoConeNeutHad   );
    ntp->Branch("v_hsimInfoConeCharHad"   ,"vector<vector<double> >",&t_v_hsimInfoConeCharHad   );
    ntp->Branch("v_hsimInfoConePdgMatched","vector<vector<double> >",&t_v_hsimInfoConePdgMatched);
    ntp->Branch("v_hsimInfoConeTotal"     ,"vector<vector<double> >",&t_v_hsimInfoConeTotal     );

    ntp->Branch("v_hsimInfoConeNMatched"  ,"vector<vector<int> >"   ,&t_v_hsimInfoConeNMatched   );
    ntp->Branch("v_hsimInfoConeNRest"     ,"vector<vector<int> >"   ,&t_v_hsimInfoConeNRest      );
    ntp->Branch("v_hsimInfoConeNPhoton"    ,"vector<vector<int> >"   ,&t_v_hsimInfoConeNPhoton     );
    ntp->Branch("v_hsimInfoConeNNeutHad"  ,"vector<vector<int> >"   ,&t_v_hsimInfoConeNNeutHad   );
    ntp->Branch("v_hsimInfoConeNCharHad"  ,"vector<vector<int> >"   ,&t_v_hsimInfoConeNCharHad   );
    ntp->Branch("v_hsimInfoConeNTotal"    ,"vector<vector<int> >"   ,&t_v_hsimInfoConeNTotal     );

    ntp->Branch("v_hsimCone"     ,"vector<vector<double> >",&t_v_hsimCone    );
  }

  ntp->Branch("v_hCone"        ,"vector<vector<double> >",&t_v_hCone       );
  ntp->Branch("v_nRecHitsCone" ,"vector<vector<int> >"   ,&t_v_nRecHitsCone);
  ntp->Branch("v_nSimHitsCone" ,"vector<vector<int> >"   ,&t_v_nSimHitsCone);

  ntp->Branch("v_distFromHotCell"        ,"vector<vector<double> >",&t_v_distFromHotCell       );
  ntp->Branch("v_ietaFromHotCell"        ,"vector<vector<int> >",&t_v_ietaFromHotCell       );
  ntp->Branch("v_iphiFromHotCell"        ,"vector<vector<int> >",&t_v_iphiFromHotCell       );

  ntp->Branch("v_RH_h3x3_ieta" ,"vector<vector<int> >"   ,&t_v_RH_h3x3_ieta);
  ntp->Branch("v_RH_h3x3_iphi" ,"vector<vector<int> >"   ,&t_v_RH_h3x3_iphi);
  ntp->Branch("v_RH_h3x3_ene"  ,"vector<vector<double> >",&t_v_RH_h3x3_ene );
  ntp->Branch("v_RH_h5x5_ieta" ,"vector<vector<int> >"   ,&t_v_RH_h5x5_ieta);
  ntp->Branch("v_RH_h5x5_iphi" ,"vector<vector<int> >"   ,&t_v_RH_h5x5_iphi);
  ntp->Branch("v_RH_h5x5_ene"  ,"vector<vector<double> >",&t_v_RH_h5x5_ene );
  ntp->Branch("v_RH_r26_ieta"  ,"vector<vector<int> >"   ,&t_v_RH_r26_ieta );
  ntp->Branch("v_RH_r26_iphi"  ,"vector<vector<int> >"   ,&t_v_RH_r26_iphi );
  ntp->Branch("v_RH_r26_ene"   ,"vector<vector<double> >",&t_v_RH_r26_ene  );
  ntp->Branch("v_RH_r44_ieta"  ,"vector<vector<int> >"   ,&t_v_RH_r44_ieta );
  ntp->Branch("v_RH_r44_iphi"  ,"vector<vector<int> >"   ,&t_v_RH_r44_iphi );
  ntp->Branch("v_RH_r44_ene"   ,"vector<vector<double> >",&t_v_RH_r44_ene  );

  ntp->Branch("v_hlTriggers" ,"vector<vector<int> >",&t_v_hlTriggers);
  ntp->Branch("v_hltHB" ,"vector<int>",&t_hltHB);
  ntp->Branch("v_hltHE" ,"vector<int>",&t_hltHE);
  ntp->Branch("v_hltL1Jet15" ,"vector<int>",&t_hltL1Jet15            );
  ntp->Branch("v_hltJet30" ,"vector<int>",&t_hltJet30            );
  ntp->Branch("v_hltJet50" ,"vector<int>",&t_hltJet50            );
  ntp->Branch("v_hltJet80" ,"vector<int>",&t_hltJet80            );
  ntp->Branch("v_hltJet110" ,"vector<int>",&t_hltJet110          );
  ntp->Branch("v_hltJet140" ,"vector<int>",&t_hltJet140          );
  ntp->Branch("v_hltJet180" ,"vector<int>",&t_hltJet180          );
  ntp->Branch("v_hltL1SingleEG5" ,"vector<int>",&t_hltL1SingleEG5        );
  ntp->Branch("v_hltZeroBias" ,"vector<int>",&t_hltZeroBias          );
  ntp->Branch("v_hltMinBiasHcal" ,"vector<int>",&t_hltMinBiasHcal        );
  ntp->Branch("v_hltMinBiasEcal" ,"vector<int>",&t_hltMinBiasEcal        );
  ntp->Branch("v_hltMinBiasPixel" ,"vector<int>",&t_hltMinBiasPixel      );
  ntp->Branch("v_hltSingleIsoTau30_Trk5" ,"vector<int>",&t_hltSingleIsoTau30_Trk5     );
  ntp->Branch("v_hltDoubleLooseIsoTau15_Trk5" ,"vector<int>",&t_hltDoubleLooseIsoTau15_Trk5);

  ntp->Branch("irun" ,"vector<unsigned int>", &t_irun);
  ntp->Branch("ievt" ,"vector<unsigned int>", &t_ievt);
  ntp->Branch("ilum" ,"vector<unsigned int>", &t_ilum);
   
}

void IsolatedTracksCone::clearTrackVectors ( )

private

Definition at line 1086 of file IsolatedTracksCone.cc.

References doMC, t_distFromHotCell_h3x3, t_distFromHotCell_h5x5, t_dtFromLeadJet, t_e3x3, t_h3x3, t_h5x5, t_hltDoubleLooseIsoTau15_Trk5, t_hltHB, t_hltHE, t_hltJet110, t_hltJet140, t_hltJet180, t_hltJet30, t_hltJet50, t_hltJet80, t_hltL1Jet15, t_hltL1SingleEG5, t_hltMinBiasEcal, t_hltMinBiasHcal, t_hltMinBiasPixel, t_hltSingleIsoTau30_Trk5, t_hltZeroBias, t_hsim3x3, t_hsim3x3CharHad, t_hsim3x3Matched, t_hsim3x3NCharHad, t_hsim3x3NeutHad, t_hsim3x3NMatched, t_hsim3x3NNeutHad, t_hsim3x3NPhoton, t_hsim3x3NRest, t_hsim3x3NTotal, t_hsim3x3PdgMatched, t_hsim3x3Photon, t_hsim3x3Rest, t_hsim3x3Total, t_hsim5x5, t_hsim5x5CharHad, t_hsim5x5Matched, t_hsim5x5NCharHad, t_hsim5x5NeutHad, t_hsim5x5NMatched, t_hsim5x5NNeutHad, t_hsim5x5NPhoton, t_hsim5x5NRest, t_hsim5x5NTotal, t_hsim5x5PdgMatched, t_hsim5x5Photon, t_hsim5x5Rest, t_hsim5x5Total, t_ietaFromHotCell_h3x3, t_ietaFromHotCell_h5x5, t_ievt, t_ilum, t_iphiFromHotCell_h3x3, t_iphiFromHotCell_h5x5, t_irun, t_nRH_h3x3, t_nRH_h5x5, t_simP, t_trkEcalEne, t_trkEta, t_trkHcalEne, t_trkNLayersCrossed, t_trkNOuterHits, t_trkP, t_trkPhi, t_trkPt, t_v_cone_hmaxNearP, t_v_cone_hmaxNearP_goodTrk, t_v_cone_hnLayers_maxNearP, t_v_cone_hnNearTRKs, t_v_cone_htrkQual_maxNearP, t_v_distFromHotCell, t_v_eDR, t_v_eMipDR, t_v_hCone, t_v_hlTriggers, t_v_hmaxNearP, t_v_hmaxNearP_goodTrk, t_v_hnLayers_maxNearP, t_v_hnNearTRKs, t_v_hsimCone, t_v_hsimInfoConeCharHad, t_v_hsimInfoConeMatched, t_v_hsimInfoConeNCharHad, t_v_hsimInfoConeNeutHad, t_v_hsimInfoConeNMatched, t_v_hsimInfoConeNNeutHad, t_v_hsimInfoConeNPhoton, t_v_hsimInfoConeNRest, t_v_hsimInfoConeNTotal, t_v_hsimInfoConePdgMatched, t_v_hsimInfoConePhoton, t_v_hsimInfoConeRest, t_v_hsimInfoConeTotal, t_v_htrkQual_maxNearP, t_v_ietaFromHotCell, t_v_iphiFromHotCell, t_v_nRecHitsCone, t_v_nSimHitsCone, t_v_RH_h3x3_ene, t_v_RH_h3x3_ieta, t_v_RH_h3x3_iphi, t_v_RH_h5x5_ene, t_v_RH_h5x5_ieta, t_v_RH_h5x5_iphi, t_v_RH_r26_ene, t_v_RH_r26_ieta, t_v_RH_r26_iphi, t_v_RH_r44_ene, t_v_RH_r44_ieta, and t_v_RH_r44_iphi.

Referenced by analyze().

                                           {

  t_v_hnNearTRKs          ->clear();   
  t_v_hnLayers_maxNearP   ->clear();   
  t_v_htrkQual_maxNearP   ->clear();   
  t_v_hmaxNearP_goodTrk   ->clear();
  t_v_hmaxNearP           ->clear();

  t_v_cone_hnNearTRKs          ->clear();   
  t_v_cone_hnLayers_maxNearP   ->clear();   
  t_v_cone_htrkQual_maxNearP   ->clear();   
  t_v_cone_hmaxNearP_goodTrk   ->clear();
  t_v_cone_hmaxNearP           ->clear();

  //  t_hScale             ->clear();
  t_trkNOuterHits      ->clear();
  t_trkNLayersCrossed  ->clear();
  t_dtFromLeadJet      ->clear();
  t_trkP               ->clear();
  t_trkPt              ->clear();
  t_trkEta             ->clear();
  t_trkPhi             ->clear();
  t_e3x3               ->clear();
  t_v_eDR              ->clear();
  t_v_eMipDR           ->clear();
  t_h3x3               ->clear();
  t_h5x5               ->clear();
  t_nRH_h3x3           ->clear();
  t_nRH_h5x5           ->clear();

  if (doMC) {
    t_simP               ->clear();
    t_hsim3x3            ->clear();
    t_hsim5x5            ->clear();
    t_hsim3x3Matched     ->clear();
    t_hsim5x5Matched     ->clear();
    t_hsim3x3Rest        ->clear();
    t_hsim5x5Rest        ->clear();
    t_hsim3x3Photon      ->clear();
    t_hsim5x5Photon      ->clear();
    t_hsim3x3NeutHad     ->clear();
    t_hsim5x5NeutHad     ->clear();
    t_hsim3x3CharHad     ->clear();
    t_hsim5x5CharHad     ->clear();
    t_hsim3x3PdgMatched  ->clear();
    t_hsim5x5PdgMatched  ->clear();
    t_hsim3x3Total       ->clear();
    t_hsim5x5Total       ->clear();

    t_hsim3x3NMatched    ->clear();
    t_hsim3x3NTotal      ->clear();
    t_hsim3x3NNeutHad    ->clear();
    t_hsim3x3NCharHad    ->clear();
    t_hsim3x3NPhoton      ->clear();
    t_hsim3x3NRest       ->clear();

    t_hsim5x5NMatched    ->clear();
    t_hsim5x5NTotal      ->clear();
    t_hsim5x5NNeutHad    ->clear();
    t_hsim5x5NCharHad    ->clear();
    t_hsim5x5NPhoton      ->clear();
    t_hsim5x5NRest       ->clear();

    t_trkHcalEne         ->clear();
    t_trkEcalEne         ->clear();
  }

  t_distFromHotCell_h3x3  ->clear();
  t_ietaFromHotCell_h3x3  ->clear();
  t_iphiFromHotCell_h3x3  ->clear();
  t_distFromHotCell_h5x5  ->clear();
  t_ietaFromHotCell_h5x5  ->clear();
  t_iphiFromHotCell_h5x5  ->clear();

  if (doMC) {
    t_v_hsimInfoConeMatched   ->clear();
    t_v_hsimInfoConeRest      ->clear();
    t_v_hsimInfoConePhoton     ->clear();
    t_v_hsimInfoConeNeutHad   ->clear();
    t_v_hsimInfoConeCharHad   ->clear();
    t_v_hsimInfoConePdgMatched->clear();  
    t_v_hsimInfoConeTotal     ->clear();

    t_v_hsimInfoConeNMatched   ->clear();
    t_v_hsimInfoConeNRest      ->clear();
    t_v_hsimInfoConeNPhoton     ->clear();
    t_v_hsimInfoConeNNeutHad   ->clear();
    t_v_hsimInfoConeNCharHad   ->clear();
    t_v_hsimInfoConeNTotal     ->clear();

    t_v_hsimCone    ->clear();
  }

  t_v_hCone       ->clear();
  t_v_nRecHitsCone->clear();
  t_v_nSimHitsCone->clear();

  t_v_distFromHotCell->clear();
  t_v_ietaFromHotCell->clear();
  t_v_iphiFromHotCell->clear();
  
  t_v_RH_h3x3_ieta ->clear();
  t_v_RH_h3x3_iphi ->clear();
  t_v_RH_h3x3_ene ->clear();
  t_v_RH_h5x5_ieta ->clear();
  t_v_RH_h5x5_iphi ->clear();
  t_v_RH_h5x5_ene ->clear();
  t_v_RH_r26_ieta  ->clear();
  t_v_RH_r26_iphi  ->clear();
  t_v_RH_r26_ene  ->clear();
  t_v_RH_r44_ieta  ->clear();
  t_v_RH_r44_iphi  ->clear();
  t_v_RH_r44_ene  ->clear();

  t_v_hlTriggers  ->clear();
  t_hltHB->clear();
  t_hltHE->clear();
  t_hltL1Jet15           ->clear();
  t_hltJet30             ->clear();
  t_hltJet50             ->clear();
  t_hltJet80             ->clear();
  t_hltJet110            ->clear();
  t_hltJet140            ->clear();
  t_hltJet180            ->clear();
  t_hltL1SingleEG5       ->clear();
  t_hltZeroBias          ->clear();
  t_hltMinBiasHcal       ->clear();
  t_hltMinBiasEcal       ->clear();
  t_hltMinBiasPixel      ->clear();
  t_hltSingleIsoTau30_Trk5     ->clear();
  t_hltDoubleLooseIsoTau15_Trk5->clear();
  
  t_irun->clear();
  t_ievt->clear();
  t_ilum->clear();



}

double IsolatedTracksCone::DeltaPhi ( double  v1, double  v2  )

private

Definition at line 1426 of file IsolatedTracksCone.cc.

References funct::abs(), corr, diffTreeTool::diff, and pi.

Referenced by analyze(), and DeltaR().

                                                        {
  // Computes the correctly normalized phi difference
  // v1, v2 = phi of object 1 and 2
  
  double pi    = 3.141592654;
  double twopi = 6.283185307;
  
  double diff = std::abs(v2 - v1);
  double corr = twopi - diff;
  if (diff < pi){ return diff;} else { return corr;} 
}

double IsolatedTracksCone::DeltaR ( double  eta1, double  phi1, double  eta2, double  phi2  )

private

Definition at line 1439 of file IsolatedTracksCone.cc.

References DEFINE_FWK_MODULE, DeltaPhi(), and mathSSE::sqrt().

                                            {
  double deta = eta1 - eta2;
  double dphi = DeltaPhi(phi1, phi2);
  return std::sqrt(deta*deta + dphi*dphi);
}

void IsolatedTracksCone::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1227 of file IsolatedTracksCone.cc.

References gather_cfg::cout, nEVT, and nEVT_failL1.

Referenced by o2olib.O2ORunMgr::executeJob().

                                {

  std::cout << "Number of Events Processed " << nEVT << " failed L1 "
        << nEVT_failL1 << std::endl;
  
}

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

private

Definition at line 1400 of file IsolatedTracksCone.cc.

References reco::TrackBase::charge(), gather_cfg::cout, reco::TrackBase::d0(), reco::TrackBase::eta(), reco::TrackBase::hitPattern(), mps_fire::i, reco::TrackBase::momentum(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfAllHits(), reco::TrackBase::numberOfLostHits(), reco::TrackBase::numberOfValidHits(), AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), reco::TrackBase::phi(), reco::HitPattern::printHitPattern(), printTrkHitPattern_, reco::TrackBase::pt(), reco::TrackBase::qoverp(), reco::TrackBase::qoverpError(), reco::TrackBase::quality(), reco::TrackBase::qualityByName(), reco::TrackBase::qualityName(), reco::TrackBase::referencePoint(), AlCaHLTBitMon_QueryRunRegistry::string, and reco::HitPattern::TRACK_HITS.

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

                                                           {
  
  std::string theTrackQuality = "highPurity";
  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);

  std::cout << " Reference Point " << pTrack->referencePoint() <<"\n"
        << " TrackMmentum " << pTrack->momentum()
        << " (pt,eta,phi)(" << pTrack->pt()<<","<<pTrack->eta()<<","<<pTrack->phi()<<")"
        << " p " << pTrack->p() << "\n"
        << " Normalized chi2 " << pTrack->normalizedChi2() <<"  charge " << pTrack->charge()
        << " qoverp() " << pTrack->qoverp() <<"+-" << pTrack->qoverpError()
        << " d0 " << pTrack->d0() << "\n"
        << " NValidHits " << pTrack->numberOfValidHits() << "  NLostHits " << pTrack->numberOfLostHits()
        << " TrackQuality " << pTrack->qualityName(trackQuality_) << " " << pTrack->quality(trackQuality_) 
        << std::endl;
  
  if( printTrkHitPattern_ ) {
    const reco::HitPattern& p = pTrack->hitPattern();
    
    for (int i=0; i<p.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
      p.printHitPattern(reco::HitPattern::TRACK_HITS, i, std::cout);
    }
  }

}

Member Data Documentation

bool IsolatedTracksCone::applyEcalIsolation_

private

Definition at line 150 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

int IsolatedTracksCone::debugEcalSimInfo_

private

Definition at line 148 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

int IsolatedTracksCone::debugTrks_

private

Definition at line 120 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

bool IsolatedTracksCone::doMC

private

Definition at line 122 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), clearTrackVectors(), and IsolatedTracksCone().

double IsolatedTracksCone::drLeadJetVeto_

private

Definition at line 125 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

edm::Service<TFileService> IsolatedTracksCone::fs

private

Definition at line 309 of file IsolatedTracksCone.h.

Referenced by BuildTree().

double IsolatedTracksCone::genPartEtaBins[5]

Definition at line 104 of file IsolatedTracksCone.h.

Referenced by beginJob().

double IsolatedTracksCone::genPartPBins[22]

Definition at line 104 of file IsolatedTracksCone.h.

Referenced by beginJob().

TH1F* IsolatedTracksCone::hRawEta

private

Definition at line 160 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TH1F* IsolatedTracksCone::hRawP

private

Definition at line 159 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TH1F* IsolatedTracksCone::hRawPhi

private

Definition at line 161 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TH1F* IsolatedTracksCone::hRawPt

private

Definition at line 158 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::leadL1JetEta

private

Definition at line 173 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::leadL1JetPhi

private

Definition at line 174 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::leadL1JetPT

private

Definition at line 172 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::maxNearTrackP_

private

Definition at line 146 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

double IsolatedTracksCone::maxTrackEta_

private

Definition at line 146 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

double IsolatedTracksCone::minTrackP_

private

Definition at line 146 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

int IsolatedTracksCone::myverbose_

private

Definition at line 123 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

const int IsolatedTracksCone::NEtaBins = 4

static

Definition at line 100 of file IsolatedTracksCone.h.

int IsolatedTracksCone::nEVT

private

Definition at line 169 of file IsolatedTracksCone.h.

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

int IsolatedTracksCone::nEVT_failL1

private

Definition at line 170 of file IsolatedTracksCone.h.

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

int IsolatedTracksCone::nFailEta

private

Definition at line 166 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nFailHighPurityQaul

private

Definition at line 164 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nFailPt

private

Definition at line 165 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nMissEcal

private

Definition at line 167 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nMissHcal

private

Definition at line 168 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

const int IsolatedTracksCone::NPBins = 21

static

Definition at line 102 of file IsolatedTracksCone.h.

const int IsolatedTracksCone::NPTBins = 21

static

Definition at line 101 of file IsolatedTracksCone.h.

int IsolatedTracksCone::nRawTRK

private

Definition at line 163 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TTree* IsolatedTracksCone::ntp

private

Definition at line 156 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nTRK

private

Definition at line 171 of file IsolatedTracksCone.h.

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

TrackAssociatorParameters IsolatedTracksCone::parameters_

private

Definition at line 153 of file IsolatedTracksCone.h.

Referenced by analyze(), Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Mixins._TypedParameterizable::copy(), and IsolatedTracksCone().

bool IsolatedTracksCone::printTrkHitPattern_

private

Definition at line 121 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone(), and printTrack().

double IsolatedTracksCone::ptMinLeadJet_

private

Definition at line 125 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

std::vector<double>* IsolatedTracksCone::t_distFromHotCell_h3x3

private

Definition at line 238 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_distFromHotCell_h5x5

private

Definition at line 241 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_dtFromLeadJet

private

Definition at line 190 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_e3x3

private

Definition at line 196 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 201 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 202 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltDoubleLooseIsoTau15_Trk5

private

Definition at line 289 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltHB

private

Definition at line 274 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltHE

private

Definition at line 275 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltJet110

private

Definition at line 280 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltJet140

private

Definition at line 281 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltJet180

private

Definition at line 282 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltJet30

private

Definition at line 277 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltJet50

private

Definition at line 278 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltJet80

private

Definition at line 279 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltL1Jet15

private

Definition at line 276 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltL1SingleEG5

private

Definition at line 283 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltMinBiasEcal

private

Definition at line 286 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltMinBiasHcal

private

Definition at line 285 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltMinBiasPixel

private

Definition at line 287 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltSingleIsoTau30_Trk5

private

Definition at line 288 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hltZeroBias

private

Definition at line 284 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 203 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 217 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 209 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim3x3NCharHad

private

Definition at line 228 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 215 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim3x3NMatched

private

Definition at line 224 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim3x3NNeutHad

private

Definition at line 227 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim3x3NPhoton

private

Definition at line 226 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim3x3NRest

private

Definition at line 225 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim3x3NTotal

private

Definition at line 229 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_hsim3x3PdgMatched

private

Definition at line 219 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 213 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 211 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_hsim3x3Total

private

Definition at line 221 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 204 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 218 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 210 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim5x5NCharHad

private

Definition at line 235 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 216 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim5x5NMatched

private

Definition at line 231 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim5x5NNeutHad

private

Definition at line 234 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim5x5NPhoton

private

Definition at line 233 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim5x5NRest

private

Definition at line 232 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_hsim5x5NTotal

private

Definition at line 236 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_hsim5x5PdgMatched

private

Definition at line 220 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 214 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 212 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_hsim5x5Total

private

Definition at line 222 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_ietaFromHotCell_h3x3

private

Definition at line 239 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_ietaFromHotCell_h5x5

private

Definition at line 242 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<unsigned int>* IsolatedTracksCone::t_ievt

private

Definition at line 306 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<unsigned int>* IsolatedTracksCone::t_ilum

private

Definition at line 307 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_iphiFromHotCell_h3x3

private

Definition at line 240 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<int>* IsolatedTracksCone::t_iphiFromHotCell_h5x5

private

Definition at line 243 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<unsigned int>* IsolatedTracksCone::t_irun

private

Definition at line 305 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_nRH_h3x3

private

Definition at line 206 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_nRH_h5x5

private

Definition at line 207 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_simP

private

Definition at line 192 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 246 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_trkEta

private

Definition at line 194 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

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

private

Definition at line 245 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_trkNLayersCrossed

private

Definition at line 189 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_trkNOuterHits

private

Definition at line 188 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_trkP

private

Definition at line 191 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_trkPhi

private

Definition at line 195 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<double>* IsolatedTracksCone::t_trkPt

private

Definition at line 193 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_cone_hmaxNearP

private

Definition at line 186 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_cone_hmaxNearP_goodTrk

private

Definition at line 185 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_cone_hnLayers_maxNearP

private

Definition at line 183 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_cone_hnNearTRKs

private

Definition at line 182 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_cone_htrkQual_maxNearP

private

Definition at line 184 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_distFromHotCell

private

Definition at line 268 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_eDR

private

Definition at line 198 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_eMipDR

private

Definition at line 199 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hCone

private

Definition at line 264 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hlTriggers

private

Definition at line 273 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hmaxNearP

private

Definition at line 180 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hmaxNearP_goodTrk

private

Definition at line 179 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hnLayers_maxNearP

private

Definition at line 177 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hnNearTRKs

private

Definition at line 176 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimCone

private

Definition at line 263 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConeCharHad

private

Definition at line 252 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConeMatched

private

Definition at line 248 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hsimInfoConeNCharHad

private

Definition at line 260 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConeNeutHad

private

Definition at line 251 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hsimInfoConeNMatched

private

Definition at line 256 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hsimInfoConeNNeutHad

private

Definition at line 259 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hsimInfoConeNPhoton

private

Definition at line 258 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hsimInfoConeNRest

private

Definition at line 257 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_hsimInfoConeNTotal

private

Definition at line 261 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConePdgMatched

private

Definition at line 253 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConePhoton

private

Definition at line 250 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConeRest

private

Definition at line 249 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_hsimInfoConeTotal

private

Definition at line 254 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_htrkQual_maxNearP

private

Definition at line 178 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_ietaFromHotCell

private

Definition at line 269 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_iphiFromHotCell

private

Definition at line 270 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_nRecHitsCone

private

Definition at line 265 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_nSimHitsCone

private

Definition at line 266 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_RH_h3x3_ene

private

Definition at line 294 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_h3x3_ieta

private

Definition at line 292 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_h3x3_iphi

private

Definition at line 293 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_RH_h5x5_ene

private

Definition at line 297 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_h5x5_ieta

private

Definition at line 295 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_h5x5_iphi

private

Definition at line 296 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_RH_r26_ene

private

Definition at line 300 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_r26_ieta

private

Definition at line 298 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_r26_iphi

private

Definition at line 299 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<double> >* IsolatedTracksCone::t_v_RH_r44_ene

private

Definition at line 303 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_r44_ieta

private

Definition at line 301 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

std::vector<std::vector<int> >* IsolatedTracksCone::t_v_RH_r44_iphi

private

Definition at line 302 of file IsolatedTracksCone.h.

Referenced by analyze(), beginJob(), BuildTree(), and clearTrackVectors().

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksCone::tok_caloEB_

private

Definition at line 141 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksCone::tok_caloEE_

private

Definition at line 142 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksCone::tok_caloHH_

private

Definition at line 143 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksCone::tok_EB_

private

Definition at line 133 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksCone::tok_EE_

private

Definition at line 134 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 137 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<HBHERecHitCollection> IsolatedTracksCone::tok_hbhe_

private

Definition at line 135 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 130 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 131 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 129 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<edm::SimTrackContainer> IsolatedTracksCone::tok_simTk_

private

Definition at line 138 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<edm::SimVertexContainer> IsolatedTracksCone::tok_simVtx_

private

Definition at line 139 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<edm::TriggerResults> IsolatedTracksCone::tok_trigger_

private

Definition at line 144 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

TrackDetectorAssociator* IsolatedTracksCone::trackAssociator_

mutableprivate

Definition at line 154 of file IsolatedTracksCone.h.

Referenced by analyze(), IsolatedTracksCone(), and ~IsolatedTracksCone().

TrackerHitAssociator::Config IsolatedTracksCone::trackerHitAssociatorConfig_

private

Definition at line 127 of file IsolatedTracksCone.h.

Referenced by analyze().

bool IsolatedTracksCone::useJetTrigger_

private

Definition at line 124 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().