IsolatedTracksCone Class Reference

#include <IsolatedTracksCone.h>

Inheritance diagram for IsolatedTracksCone:

Public Member Functions
	IsolatedTracksCone (const edm::ParameterSet &)
	~IsolatedTracksCone ()
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
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob (const edm::EventSetup &)
void	BuildTree ()
void	clearTrackVectors ()
double	DeltaPhi (double v1, double v2)
double	DeltaR (double eta1, double phi1, double eta2, double phi2)
virtual void	endJob ()
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_
bool	useJetTrigger_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
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 91 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_, Parameters::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_.

                                                                     {

  //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

(

)

Definition at line 107 of file IsolatedTracksCone.cc.

References trackAssociator_.

                                        {
  delete  trackAssociator_;
}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 111 of file IsolatedTracksCone.cc.

References funct::abs(), 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(), first, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getClosestCell(), spr::getDistInPlaneTrackDir(), TrackDetectorAssociator::getFreeTrajectoryState(), edm::Event::getLuminosityBlock(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, HcalEndcap, reco::TrackBase::hitPattern(), hRawEta, hRawP, hRawPhi, hRawPt, i, edm::EventBase::id(), 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< class >::product(), reco::TrackBase::pt(), ptMinLeadJet_, reco::TrackBase::quality(), reco::TrackBase::qualityByName(), DetId::rawId(), edm::EventID::run(), 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_, 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<IdealGeometryRecord>().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
  TrackerHitAssociator* associate=0;
  if (doMC) associate = new TrackerHitAssociator(iEvent);
  

  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++;
  
  
  delete associate;
}

void IsolatedTracksCone::beginJob ( const edm::EventSetup &  )

privatevirtual

Definition at line 941 of file IsolatedTracksCone.cc.

References BuildTree(), doMC, genPartEtaBins, genPartPBins, 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 1252 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 1103 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 1443 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 1456 of file IsolatedTracksCone.cc.

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

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

void IsolatedTracksCone::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1244 of file IsolatedTracksCone.cc.

References gather_cfg::cout, nEVT, and nEVT_failL1.

                                {

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

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

private

Definition at line 1417 of file IsolatedTracksCone.cc.

References reco::TrackBase::charge(), gather_cfg::cout, reco::TrackBase::d0(), reco::TrackBase::eta(), reco::TrackBase::hitPattern(), i, reco::TrackBase::momentum(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfHits(), reco::TrackBase::numberOfLostHits(), reco::TrackBase::numberOfValidHits(), reco::TrackBase::p(), AlCaHLTBitMon_ParallelJobs::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(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                           {
  
  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.numberOfHits(); i++) {
      p.printHitPattern(i, std::cout);
    }
  }

}

Member Data Documentation

bool IsolatedTracksCone::applyEcalIsolation_

private

Definition at line 144 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

int IsolatedTracksCone::debugEcalSimInfo_

private

Definition at line 142 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

int IsolatedTracksCone::debugTrks_

private

Definition at line 117 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

bool IsolatedTracksCone::doMC

private

Definition at line 119 of file IsolatedTracksCone.h.

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

double IsolatedTracksCone::drLeadJetVeto_

private

Definition at line 122 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

edm::Service<TFileService> IsolatedTracksCone::fs

private

Definition at line 307 of file IsolatedTracksCone.h.

Referenced by BuildTree().

double IsolatedTracksCone::genPartEtaBins[5]

Definition at line 101 of file IsolatedTracksCone.h.

Referenced by beginJob().

double IsolatedTracksCone::genPartPBins[22]

Definition at line 101 of file IsolatedTracksCone.h.

Referenced by beginJob().

TH1F* IsolatedTracksCone::hRawEta

private

Definition at line 154 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TH1F* IsolatedTracksCone::hRawP

private

Definition at line 153 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TH1F* IsolatedTracksCone::hRawPhi

private

Definition at line 155 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TH1F* IsolatedTracksCone::hRawPt

private

Definition at line 152 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::leadL1JetEta

private

Definition at line 169 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::leadL1JetPhi

private

Definition at line 170 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::leadL1JetPT

private

Definition at line 168 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

double IsolatedTracksCone::maxNearTrackP_

private

Definition at line 140 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

double IsolatedTracksCone::maxTrackEta_

private

Definition at line 140 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

double IsolatedTracksCone::minTrackP_

private

Definition at line 140 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

int IsolatedTracksCone::myverbose_

private

Definition at line 120 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone().

const int IsolatedTracksCone::NEtaBins = 4

static

Definition at line 97 of file IsolatedTracksCone.h.

int IsolatedTracksCone::nEVT

private

Definition at line 163 of file IsolatedTracksCone.h.

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

int IsolatedTracksCone::nEVT_failL1

private

Definition at line 164 of file IsolatedTracksCone.h.

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

int IsolatedTracksCone::nFailEta

private

Definition at line 160 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nFailHighPurityQaul

private

Definition at line 158 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nFailPt

private

Definition at line 159 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nMissEcal

private

Definition at line 161 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nMissHcal

private

Definition at line 162 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

const int IsolatedTracksCone::NPBins = 21

static

Definition at line 99 of file IsolatedTracksCone.h.

const int IsolatedTracksCone::NPTBins = 21

static

Definition at line 98 of file IsolatedTracksCone.h.

int IsolatedTracksCone::nRawTRK

private

Definition at line 157 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

TTree* IsolatedTracksCone::ntp

private

Definition at line 150 of file IsolatedTracksCone.h.

Referenced by analyze(), and BuildTree().

int IsolatedTracksCone::nTRK

private

Definition at line 165 of file IsolatedTracksCone.h.

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

TrackAssociatorParameters IsolatedTracksCone::parameters_

private

Definition at line 147 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 118 of file IsolatedTracksCone.h.

Referenced by IsolatedTracksCone(), and printTrack().

double IsolatedTracksCone::ptMinLeadJet_

private

Definition at line 122 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 235 of file IsolatedTracksCone.h.

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

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

private

Definition at line 238 of file IsolatedTracksCone.h.

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

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

private

Definition at line 187 of file IsolatedTracksCone.h.

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

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

private

Definition at line 193 of file IsolatedTracksCone.h.

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

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

private

Definition at line 198 of file IsolatedTracksCone.h.

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

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

private

Definition at line 199 of file IsolatedTracksCone.h.

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

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

private

Definition at line 286 of file IsolatedTracksCone.h.

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

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

private

Definition at line 271 of file IsolatedTracksCone.h.

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

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

private

Definition at line 272 of file IsolatedTracksCone.h.

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

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

private

Definition at line 277 of file IsolatedTracksCone.h.

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

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

private

Definition at line 278 of file IsolatedTracksCone.h.

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

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

private

Definition at line 279 of file IsolatedTracksCone.h.

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

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

private

Definition at line 274 of file IsolatedTracksCone.h.

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

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

private

Definition at line 275 of file IsolatedTracksCone.h.

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

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

private

Definition at line 276 of file IsolatedTracksCone.h.

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

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

private

Definition at line 273 of file IsolatedTracksCone.h.

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

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

private

Definition at line 280 of file IsolatedTracksCone.h.

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

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

private

Definition at line 283 of file IsolatedTracksCone.h.

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

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

private

Definition at line 282 of file IsolatedTracksCone.h.

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

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

private

Definition at line 284 of file IsolatedTracksCone.h.

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

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

private

Definition at line 285 of file IsolatedTracksCone.h.

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

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

private

Definition at line 281 of file IsolatedTracksCone.h.

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

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

private

Definition at line 200 of file IsolatedTracksCone.h.

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

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

private

Definition at line 214 of file IsolatedTracksCone.h.

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

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

private

Definition at line 206 of file IsolatedTracksCone.h.

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

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

private

Definition at line 225 of file IsolatedTracksCone.h.

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

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

private

Definition at line 212 of file IsolatedTracksCone.h.

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

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

private

Definition at line 221 of file IsolatedTracksCone.h.

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

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

private

Definition at line 224 of file IsolatedTracksCone.h.

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

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

private

Definition at line 223 of file IsolatedTracksCone.h.

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

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

private

Definition at line 222 of file IsolatedTracksCone.h.

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

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

private

Definition at line 226 of file IsolatedTracksCone.h.

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

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

private

Definition at line 216 of file IsolatedTracksCone.h.

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

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

private

Definition at line 210 of file IsolatedTracksCone.h.

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

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

private

Definition at line 208 of file IsolatedTracksCone.h.

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

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

private

Definition at line 218 of file IsolatedTracksCone.h.

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

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

private

Definition at line 201 of file IsolatedTracksCone.h.

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

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

private

Definition at line 215 of file IsolatedTracksCone.h.

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

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

private

Definition at line 207 of file IsolatedTracksCone.h.

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

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

private

Definition at line 232 of file IsolatedTracksCone.h.

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

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

private

Definition at line 213 of file IsolatedTracksCone.h.

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

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

private

Definition at line 228 of file IsolatedTracksCone.h.

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

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

private

Definition at line 231 of file IsolatedTracksCone.h.

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

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

private

Definition at line 230 of file IsolatedTracksCone.h.

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

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

private

Definition at line 229 of file IsolatedTracksCone.h.

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

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

private

Definition at line 233 of file IsolatedTracksCone.h.

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

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

private

Definition at line 217 of file IsolatedTracksCone.h.

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

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

private

Definition at line 211 of file IsolatedTracksCone.h.

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

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

private

Definition at line 209 of file IsolatedTracksCone.h.

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

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

private

Definition at line 219 of file IsolatedTracksCone.h.

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

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

private

Definition at line 236 of file IsolatedTracksCone.h.

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

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

private

Definition at line 239 of file IsolatedTracksCone.h.

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

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

private

Definition at line 303 of file IsolatedTracksCone.h.

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

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

private

Definition at line 304 of file IsolatedTracksCone.h.

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

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

private

Definition at line 237 of file IsolatedTracksCone.h.

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

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

private

Definition at line 240 of file IsolatedTracksCone.h.

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

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

private

Definition at line 302 of file IsolatedTracksCone.h.

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

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

private

Definition at line 203 of file IsolatedTracksCone.h.

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

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

private

Definition at line 204 of file IsolatedTracksCone.h.

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

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

private

Definition at line 189 of file IsolatedTracksCone.h.

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

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

private

Definition at line 243 of file IsolatedTracksCone.h.

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

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

private

Definition at line 191 of file IsolatedTracksCone.h.

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

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

private

Definition at line 242 of file IsolatedTracksCone.h.

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

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

private

Definition at line 186 of file IsolatedTracksCone.h.

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

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

private

Definition at line 185 of file IsolatedTracksCone.h.

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

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

private

Definition at line 188 of file IsolatedTracksCone.h.

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

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

private

Definition at line 192 of file IsolatedTracksCone.h.

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

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

private

Definition at line 190 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 182 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 181 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 179 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 178 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 180 of file IsolatedTracksCone.h.

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

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

private

Definition at line 265 of file IsolatedTracksCone.h.

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

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

private

Definition at line 195 of file IsolatedTracksCone.h.

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

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

private

Definition at line 196 of file IsolatedTracksCone.h.

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

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

private

Definition at line 261 of file IsolatedTracksCone.h.

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

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

private

Definition at line 270 of file IsolatedTracksCone.h.

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

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

private

Definition at line 176 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 175 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 173 of file IsolatedTracksCone.h.

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

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

private

Definition at line 172 of file IsolatedTracksCone.h.

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

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

private

Definition at line 260 of file IsolatedTracksCone.h.

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

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

private

Definition at line 249 of file IsolatedTracksCone.h.

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

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

private

Definition at line 245 of file IsolatedTracksCone.h.

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

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

private

Definition at line 257 of file IsolatedTracksCone.h.

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

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

private

Definition at line 248 of file IsolatedTracksCone.h.

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

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

private

Definition at line 253 of file IsolatedTracksCone.h.

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

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

private

Definition at line 256 of file IsolatedTracksCone.h.

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

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

private

Definition at line 255 of file IsolatedTracksCone.h.

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

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

private

Definition at line 254 of file IsolatedTracksCone.h.

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

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

private

Definition at line 258 of file IsolatedTracksCone.h.

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

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

private

Definition at line 250 of file IsolatedTracksCone.h.

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

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

private

Definition at line 247 of file IsolatedTracksCone.h.

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

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

private

Definition at line 246 of file IsolatedTracksCone.h.

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

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

private

Definition at line 251 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 174 of file IsolatedTracksCone.h.

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

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

private

Definition at line 266 of file IsolatedTracksCone.h.

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

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

private

Definition at line 267 of file IsolatedTracksCone.h.

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

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

private

Definition at line 262 of file IsolatedTracksCone.h.

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

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

private

Definition at line 263 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 291 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 289 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 290 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 294 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 292 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 293 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 297 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 295 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 296 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 300 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 298 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 299 of file IsolatedTracksCone.h.

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

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

private

Definition at line 135 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 136 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 137 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksCone::tok_EB_

private

Definition at line 128 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksCone::tok_EE_

private

Definition at line 129 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 132 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

edm::EDGetTokenT<HBHERecHitCollection> IsolatedTracksCone::tok_hbhe_

private

Definition at line 131 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 125 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 126 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 124 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 133 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 134 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

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

private

Definition at line 138 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().

TrackDetectorAssociator* IsolatedTracksCone::trackAssociator_

mutableprivate

Definition at line 148 of file IsolatedTracksCone.h.

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

bool IsolatedTracksCone::useJetTrigger_

private

Definition at line 121 of file IsolatedTracksCone.h.

Referenced by analyze(), and IsolatedTracksCone().