IsolatedTracksHcalScale Class Reference

Inheritance diagram for IsolatedTracksHcalScale:

Public Member Functions
	IsolatedTracksHcalScale (const edm::ParameterSet &)
	~IsolatedTracksHcalScale () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Private Attributes
double	a_charIsoR_
double	a_coneR_
double	a_mipR_
double	a_neutIsoR_
bool	doMC_
std::string	minQuality_
int	myverbose_
int	nEventProc_
spr::trackSelectionParameters	selectionParameters_
int	t_Bunch
std::vector< double > *	t_conehmaxNearP
std::vector< double > *	t_e11x11_20Sig
std::vector< double > *	t_e15x15_20Sig
std::vector< double > *	t_eECALDR
std::vector< double > *	t_eECALDR_1
std::vector< double > *	t_eECALDR_2
std::vector< double > *	t_eHCALDR
std::vector< double > *	t_eHCALDRHB
std::vector< double > *	t_eMipDR
std::vector< double > *	t_eMipDR_1
std::vector< double > *	t_eMipDR_2
int	t_EvtNo
std::vector< double > *	t_hCone
std::vector< double > *	t_hConeHB
std::vector< double > *	t_hsim
std::vector< double > *	t_hsimInfoCharHad
std::vector< double > *	t_hsimInfoMatched
std::vector< int > *	t_hsimInfoNCharHad
std::vector< double > *	t_hsimInfoNeutHad
std::vector< int > *	t_hsimInfoNMatched
std::vector< int > *	t_hsimInfoNNeutHad
std::vector< int > *	t_hsimInfoNPhoton
std::vector< int > *	t_hsimInfoNRest
std::vector< int > *	t_hsimInfoNTotal
std::vector< double > *	t_hsimInfoPdgMatched
std::vector< double > *	t_hsimInfoPhoton
std::vector< double > *	t_hsimInfoRest
std::vector< double > *	t_hsimInfoTotal
int	t_Lumi
std::vector< int > *	t_nSimHits
int	t_nTracks
int	t_RunNo
std::vector< double > *	t_trackEta
std::vector< double > *	t_trackHcalEta
std::vector< double > *	t_trackHcalPhi
std::vector< double > *	t_trackP
std::vector< double > *	t_trackPhi
std::vector< double > *	t_trackPt
std::string	theTrackQuality_
double	tMaxE_
double	tMinE_
edm::EDGetTokenT< reco::BeamSpot >	tok_bs_
edm::EDGetTokenT < edm::PCaloHitContainer >	tok_caloEB_
edm::EDGetTokenT < edm::PCaloHitContainer >	tok_caloEE_
edm::EDGetTokenT < edm::PCaloHitContainer >	tok_caloHH_
edm::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_
edm::EDGetTokenT < EcalRecHitCollection >	tok_EB_
edm::ESGetToken < EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_
edm::EDGetTokenT < EcalRecHitCollection >	tok_EE_
edm::EDGetTokenT < reco::TrackCollection >	tok_genTrack_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_magField_
edm::EDGetTokenT < reco::VertexCollection >	tok_recVtx_
edm::ESGetToken < EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_
edm::EDGetTokenT < edm::SimTrackContainer >	tok_simTk_
edm::EDGetTokenT < edm::SimVertexContainer >	tok_simVtx_
TrackerHitAssociator::Config	trackerHitAssociatorConfig_
TTree *	tree_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 99 of file IsolatedTracksHcalScale.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 161 of file IsolatedTracksHcalScale.cc.

References a_charIsoR_, a_coneR_, a_mipR_, a_neutIsoR_, doMC_, edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::InputTag, TFileService::kSharedResource, spr::trackSelectionParameters::maxChi2, spr::trackSelectionParameters::maxDpOverP, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, spr::trackSelectionParameters::minLayerCrossed, spr::trackSelectionParameters::minOuterHit, spr::trackSelectionParameters::minPt, spr::trackSelectionParameters::minQuality, myverbose_, reco::TrackBase::qualityByName(), selectionParameters_, theTrackQuality_, tMaxE_, tMinE_, tok_bs_, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_magField_, tok_recVtx_, tok_sevlv_, tok_simTk_, tok_simVtx_, and Reconstruction_hiPF_cff::trackQuality.

    : doMC_(iConfig.getUntrackedParameter<bool>("DoMC", false)),
      myverbose_(iConfig.getUntrackedParameter<int>("Verbosity", 5)),
      theTrackQuality_(iConfig.getUntrackedParameter<std::string>("TrackQuality", "highPurity")),
      a_mipR_(iConfig.getUntrackedParameter<double>("ConeRadiusMIP", 14.0)),
      a_coneR_(iConfig.getUntrackedParameter<double>("ConeRadius", 34.98)),
      tMinE_(iConfig.getUntrackedParameter<double>("TimeMinCutECAL", -500.)),
      tMaxE_(iConfig.getUntrackedParameter<double>("TimeMaxCutECAL", 500.)),
      trackerHitAssociatorConfig_(consumesCollector()) {
  usesResource(TFileService::kSharedResource);

  //now do what ever initialization is needed
  reco::TrackBase::TrackQuality trackQuality = reco::TrackBase::qualityByName(theTrackQuality_);
  selectionParameters_.minPt = iConfig.getUntrackedParameter<double>("MinTrackPt", 10.0);
  selectionParameters_.minQuality = trackQuality;
  selectionParameters_.maxDxyPV = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.2);
  selectionParameters_.maxDzPV = iConfig.getUntrackedParameter<double>("MaxDzPV", 5.0);
  selectionParameters_.maxChi2 = iConfig.getUntrackedParameter<double>("MaxChi2", 5.0);
  selectionParameters_.maxDpOverP = iConfig.getUntrackedParameter<double>("MaxDpOverP", 0.1);
  selectionParameters_.minOuterHit = iConfig.getUntrackedParameter<int>("MinOuterHit", 4);
  selectionParameters_.minLayerCrossed = iConfig.getUntrackedParameter<int>("MinLayerCrossed", 8);
  selectionParameters_.maxInMiss = iConfig.getUntrackedParameter<int>("MaxInMiss", 0);
  selectionParameters_.maxOutMiss = iConfig.getUntrackedParameter<int>("MaxOutMiss", 0);
  a_charIsoR_ = a_coneR_ + 28.9;
  a_neutIsoR_ = a_charIsoR_ * 0.726;

  tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
  tok_bs_ = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
  tok_simTk_ = consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
  tok_simVtx_ = consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
  tok_caloEB_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEB"));
  tok_caloEE_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEE"));
  tok_caloHH_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalHits"));

  if (myverbose_ >= 0) {
    edm::LogVerbatim("IsoTrack") << "Parameters read from config file \n"
                                 << " doMC " << doMC_ << "\t myverbose " << myverbose_ << "\t minPt "
                                 << selectionParameters_.minPt << "\t theTrackQuality " << theTrackQuality_
                                 << "\t minQuality " << selectionParameters_.minQuality << "\t maxDxyPV "
                                 << selectionParameters_.maxDxyPV << "\t maxDzPV " << selectionParameters_.maxDzPV
                                 << "\t maxChi2 " << selectionParameters_.maxChi2 << "\t maxDpOverP "
                                 << selectionParameters_.maxDpOverP << "\t minOuterHit "
                                 << selectionParameters_.minOuterHit << "\t minLayerCrossed "
                                 << selectionParameters_.minLayerCrossed << "\t maxInMiss "
                                 << selectionParameters_.maxInMiss << "\t maxOutMiss "
                                 << selectionParameters_.maxOutMiss << "\t a_coneR " << a_coneR_ << "\t a_charIsoR "
                                 << a_charIsoR_ << "\t a_neutIsoR " << a_neutIsoR_ << "\t a_mipR " << a_mipR_
                                 << "\t time Range (" << tMinE_ << ":" << tMaxE_ << ")";
  }

  tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
  tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
  tok_magField_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
  tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
  tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
}

IsolatedTracksHcalScale::~IsolatedTracksHcalScale ( )

inlineoverride

Definition at line 102 of file IsolatedTracksHcalScale.cc.

{}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 243 of file IsolatedTracksHcalScale.cc.

References a_charIsoR_, a_coneR_, a_mipR_, a_neutIsoR_, ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationCone(), clearTreeVectors(), doMC_, spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), edm::EventID::event(), edm::Event::getByToken(), edm::EventSetup::getData(), spr::goodTrack(), HcalBarrel, edm::EventBase::id(), HcalDetId::ieta(), iEvent, HcalDetId::iphi(), edm::HandleBase::isValid(), edm::EventBase::luminosityBlock(), reco::TrackBase::momentum(), myverbose_, nEventProc_, BeamSpotPI::nTracks, reco::TrackBase::p(), edm::Handle< T >::product(), spr::propagateCALO(), reco::TrackBase::pt(), edm::EventID::run(), selectionParameters_, t_Bunch, t_conehmaxNearP, t_e11x11_20Sig, t_e15x15_20Sig, t_eECALDR, t_eECALDR_1, t_eECALDR_2, t_eHCALDR, t_eHCALDRHB, t_eMipDR, t_eMipDR_1, t_eMipDR_2, t_EvtNo, t_hCone, t_hConeHB, t_hsim, t_hsimInfoCharHad, t_hsimInfoMatched, t_hsimInfoNCharHad, t_hsimInfoNeutHad, t_hsimInfoNMatched, t_hsimInfoNNeutHad, t_hsimInfoNPhoton, t_hsimInfoNRest, t_hsimInfoNTotal, t_hsimInfoPdgMatched, t_hsimInfoPhoton, t_hsimInfoRest, t_hsimInfoTotal, t_Lumi, t_nSimHits, t_RunNo, t_trackEta, t_trackHcalEta, t_trackHcalPhi, t_trackP, t_trackPhi, t_trackPt, theTrackQuality_, tMaxE_, tMinE_, tok_bs_, tok_caloEB_, tok_caloEE_, tok_caloHH_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_magField_, tok_recVtx_, tok_sevlv_, tok_simTk_, tok_simVtx_, trackerHitAssociatorConfig_, tree_, x, y, and z.

                                                                                         {
  const CaloGeometry *geo = &iSetup.getData(tok_geom_);
  const MagneticField *bField = &iSetup.getData(tok_magField_);
  const EcalChannelStatus *theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
  const EcalSeverityLevelAlgo *sevlv = &iSetup.getData(tok_sevlv_);
  const CaloTopology *caloTopology = &iSetup.getData(tok_caloTopology_);

  clearTreeVectors();

  ++nEventProc_;

  t_RunNo = iEvent.id().run();
  t_EvtNo = iEvent.id().event();
  t_Lumi = iEvent.luminosityBlock();
  t_Bunch = iEvent.bunchCrossing();
  if (myverbose_ > 0)
    edm::LogVerbatim("IsoTrack") << nEventProc_ << " Run " << t_RunNo << " Event " << t_EvtNo << " Lumi " << t_Lumi
                                 << " Bunch " << t_Bunch;

  edm::Handle<reco::TrackCollection> trkCollection;
  iEvent.getByToken(tok_genTrack_, trkCollection);

  edm::Handle<reco::VertexCollection> recVtxs;
  iEvent.getByToken(tok_recVtx_, recVtxs);

  // Get the beamspot
  edm::Handle<reco::BeamSpot> beamSpotH;
  iEvent.getByToken(tok_bs_, beamSpotH);

  math::XYZPoint leadPV(0, 0, 0);
  if (!recVtxs->empty() && !((*recVtxs)[0].isFake())) {
    leadPV = math::XYZPoint((*recVtxs)[0].x(), (*recVtxs)[0].y(), (*recVtxs)[0].z());
  } else if (beamSpotH.isValid()) {
    leadPV = beamSpotH->position();
  }

  if (myverbose_ > 0) {
    edm::LogVerbatim("IsoTrack") << "Primary Vertex " << leadPV;
    if (beamSpotH.isValid())
      edm::LogVerbatim("IsoTrack") << "Beam Spot " << beamSpotH->position();
  }

  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, (myverbose_ > 2));
  std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;

  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);

  edm::Handle<HBHERecHitCollection> hbhe;
  iEvent.getByToken(tok_hbhe_, hbhe);
  const HBHERecHitCollection Hithbhe = *(hbhe.product());

  //get Handles to SimTracks and SimHits
  edm::Handle<edm::SimTrackContainer> SimTk;
  edm::Handle<edm::SimVertexContainer> SimVtx;

  //get Handles to PCaloHitContainers of eb/ee/hbhe
  edm::Handle<edm::PCaloHitContainer> pcaloeb;
  edm::Handle<edm::PCaloHitContainer> pcaloee;
  edm::Handle<edm::PCaloHitContainer> pcalohh;

  //associates tracker rechits/simhits to a track
  std::unique_ptr<TrackerHitAssociator> associate;

  if (doMC_) {
    iEvent.getByToken(tok_simTk_, SimTk);
    iEvent.getByToken(tok_simVtx_, SimVtx);
    iEvent.getByToken(tok_caloEB_, pcaloeb);
    iEvent.getByToken(tok_caloEE_, pcaloee);
    iEvent.getByToken(tok_caloHH_, pcalohh);
    associate = std::make_unique<TrackerHitAssociator>(iEvent, trackerHitAssociatorConfig_);
  }

  unsigned int nTracks = 0;
  for (trkDetItr = trkCaloDirections.begin(), nTracks = 0; trkDetItr != trkCaloDirections.end();
       trkDetItr++, nTracks++) {
    const reco::Track *pTrack = &(*(trkDetItr->trkItr));
    if (spr::goodTrack(pTrack, leadPV, selectionParameters_, (myverbose_ > 2)) && trkDetItr->okECAL &&
        trkDetItr->okHCAL) {
      int nRH_eMipDR = 0, nRH_eDR = 0, nNearTRKs = 0, nRecHitsCone = -99;
      double distFromHotCell = -99.0, distFromHotCell2 = -99.0;
      int ietaHotCell = -99, iphiHotCell = -99;
      int ietaHotCell2 = -99, iphiHotCell2 = -99;
      GlobalPoint gposHotCell(0., 0., 0.), gposHotCell2(0., 0., 0.);
      std::vector<DetId> coneRecHitDetIds, coneRecHitDetIds2;
      std::pair<double, bool> e11x11_20SigP, e15x15_20SigP;
      double hCone = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_coneR_,
                                     trkDetItr->directionHCAL,
                                     nRecHitsCone,
                                     coneRecHitDetIds,
                                     distFromHotCell,
                                     ietaHotCell,
                                     iphiHotCell,
                                     gposHotCell,
                                     -1);
      double hConeHB = spr::eCone_hcal(geo,
                                       hbhe,
                                       trkDetItr->pointHCAL,
                                       trkDetItr->pointECAL,
                                       a_coneR_,
                                       trkDetItr->directionHCAL,
                                       nRecHitsCone,
                                       coneRecHitDetIds,
                                       distFromHotCell,
                                       ietaHotCell,
                                       iphiHotCell,
                                       gposHotCell,
                                       (int)(HcalBarrel));
      double eHCALDR = spr::eCone_hcal(geo,
                                       hbhe,
                                       trkDetItr->pointHCAL,
                                       trkDetItr->pointECAL,
                                       a_charIsoR_,
                                       trkDetItr->directionHCAL,
                                       nRecHitsCone,
                                       coneRecHitDetIds2,
                                       distFromHotCell2,
                                       ietaHotCell2,
                                       iphiHotCell2,
                                       gposHotCell2,
                                       -1);
      double eHCALDRHB = spr::eCone_hcal(geo,
                                         hbhe,
                                         trkDetItr->pointHCAL,
                                         trkDetItr->pointECAL,
                                         a_charIsoR_,
                                         trkDetItr->directionHCAL,
                                         nRecHitsCone,
                                         coneRecHitDetIds2,
                                         distFromHotCell2,
                                         ietaHotCell2,
                                         iphiHotCell2,
                                         gposHotCell2,
                                         (int)(HcalBarrel));

      double conehmaxNearP =
          spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, (myverbose_ > 3));

      double eMipDR = spr::eCone_ecal(geo,
                                      barrelRecHitsHandle,
                                      endcapRecHitsHandle,
                                      trkDetItr->pointHCAL,
                                      trkDetItr->pointECAL,
                                      a_mipR_,
                                      trkDetItr->directionECAL,
                                      nRH_eMipDR);
      double eECALDR = spr::eCone_ecal(geo,
                                       barrelRecHitsHandle,
                                       endcapRecHitsHandle,
                                       trkDetItr->pointHCAL,
                                       trkDetItr->pointECAL,
                                       a_neutIsoR_,
                                       trkDetItr->directionECAL,
                                       nRH_eDR);
      double eMipDR_1 = spr::eCone_ecal(geo,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        trkDetItr->pointHCAL,
                                        trkDetItr->pointECAL,
                                        a_mipR_,
                                        trkDetItr->directionECAL,
                                        nRH_eMipDR,
                                        0.030,
                                        0.150);
      double eECALDR_1 = spr::eCone_ecal(geo,
                                         barrelRecHitsHandle,
                                         endcapRecHitsHandle,
                                         trkDetItr->pointHCAL,
                                         trkDetItr->pointECAL,
                                         a_neutIsoR_,
                                         trkDetItr->directionECAL,
                                         nRH_eDR,
                                         0.030,
                                         0.150);
      double eMipDR_2 = spr::eCone_ecal(geo,
                                        barrelRecHitsHandle,
                                        endcapRecHitsHandle,
                                        trkDetItr->pointHCAL,
                                        trkDetItr->pointECAL,
                                        a_mipR_,
                                        trkDetItr->directionECAL,
                                        nRH_eMipDR,
                                        0.060,
                                        0.300);
      double eECALDR_2 = spr::eCone_ecal(geo,
                                         barrelRecHitsHandle,
                                         endcapRecHitsHandle,
                                         trkDetItr->pointHCAL,
                                         trkDetItr->pointECAL,
                                         a_neutIsoR_,
                                         trkDetItr->directionECAL,
                                         nRH_eDR,
                                         0.060,
                                         0.300);

      HcalDetId closestCell = (HcalDetId)(trkDetItr->detIdHCAL);

      e11x11_20SigP = spr::eECALmatrix(trkDetItr->detIdECAL,
                                       barrelRecHitsHandle,
                                       endcapRecHitsHandle,
                                       *theEcalChStatus,
                                       geo,
                                       caloTopology,
                                       sevlv,
                                       5,
                                       5,
                                       0.060,
                                       0.300,
                                       tMinE_,
                                       tMaxE_);
      e15x15_20SigP = spr::eECALmatrix(trkDetItr->detIdECAL,
                                       barrelRecHitsHandle,
                                       endcapRecHitsHandle,
                                       *theEcalChStatus,
                                       geo,
                                       caloTopology,
                                       sevlv,
                                       7,
                                       7,
                                       0.060,
                                       0.300,
                                       tMinE_,
                                       tMaxE_);

      // Fill the tree Branches here
      t_trackP->push_back(pTrack->p());
      t_trackPt->push_back(pTrack->pt());
      t_trackEta->push_back(pTrack->momentum().eta());
      t_trackPhi->push_back(pTrack->momentum().phi());
      t_trackHcalEta->push_back(closestCell.ieta());
      t_trackHcalPhi->push_back(closestCell.iphi());
      t_hCone->push_back(hCone);
      t_conehmaxNearP->push_back(conehmaxNearP);
      t_eMipDR->push_back(eMipDR);
      t_eECALDR->push_back(eECALDR);
      t_eHCALDR->push_back(eHCALDR);
      t_e11x11_20Sig->push_back(e11x11_20SigP.first);
      t_e15x15_20Sig->push_back(e15x15_20SigP.first);
      t_eMipDR_1->push_back(eMipDR_1);
      t_eECALDR_1->push_back(eECALDR_1);
      t_eMipDR_2->push_back(eMipDR_2);
      t_eECALDR_2->push_back(eECALDR_2);
      t_hConeHB->push_back(hConeHB);
      t_eHCALDRHB->push_back(eHCALDRHB);

      if (myverbose_ > 0) {
        edm::LogVerbatim("IsoTrack") << "Track p " << pTrack->p() << " pt " << pTrack->pt() << " eta "
                                     << pTrack->momentum().eta() << " phi " << pTrack->momentum().phi()
                                     << " ieta/iphi (" << closestCell.ieta() << ", " << closestCell.iphi()
                                     << ") Energy in cone " << hCone << " Charge Isolation " << conehmaxNearP
                                     << " eMIP (" << eMipDR << ", " << eMipDR_1 << ", " << eMipDR_2 << ")"
                                     << " Neutral isolation (ECAL) (" << eECALDR - eMipDR << ", "
                                     << eECALDR_1 - eMipDR_1 << ", " << eECALDR_2 - eMipDR_2 << ") (ECAL NxN) "
                                     << e15x15_20SigP.first - e11x11_20SigP.first << " (HCAL) " << eHCALDR - hCone;
      }

      if (doMC_) {
        int nSimHits = -999;
        double hsim;
        std::map<std::string, double> hsimInfo;
        std::vector<int> multiplicity;
        hsim = spr::eCone_hcal(
            geo, pcalohh, trkDetItr->pointHCAL, trkDetItr->pointECAL, a_coneR_, trkDetItr->directionHCAL, nSimHits);
        hsimInfo = spr::eHCALSimInfoCone(iEvent,
                                         pcalohh,
                                         SimTk,
                                         SimVtx,
                                         pTrack,
                                         *associate,
                                         geo,
                                         trkDetItr->pointHCAL,
                                         trkDetItr->pointECAL,
                                         a_coneR_,
                                         trkDetItr->directionHCAL,
                                         multiplicity);

        t_hsimInfoMatched->push_back(hsimInfo["eMatched"]);
        t_hsimInfoRest->push_back(hsimInfo["eRest"]);
        t_hsimInfoPhoton->push_back(hsimInfo["eGamma"]);
        t_hsimInfoNeutHad->push_back(hsimInfo["eNeutralHad"]);
        t_hsimInfoCharHad->push_back(hsimInfo["eChargedHad"]);
        t_hsimInfoPdgMatched->push_back(hsimInfo["pdgMatched"]);
        t_hsimInfoTotal->push_back(hsimInfo["eTotal"]);

        t_hsimInfoNMatched->push_back(multiplicity.at(0));
        t_hsimInfoNTotal->push_back(multiplicity.at(1));
        t_hsimInfoNNeutHad->push_back(multiplicity.at(2));
        t_hsimInfoNCharHad->push_back(multiplicity.at(3));
        t_hsimInfoNPhoton->push_back(multiplicity.at(4));
        t_hsimInfoNRest->push_back(multiplicity.at(5));

        t_hsim->push_back(hsim);
        t_nSimHits->push_back(nSimHits);

        if (myverbose_ > 0) {
          edm::LogVerbatim("IsoTrack") << "Matched (E) " << hsimInfo["eMatched"] << " (N) " << multiplicity.at(0)
                                       << " Rest (E) " << hsimInfo["eRest"] << " (N) " << multiplicity.at(1)
                                       << " Gamma (E) " << hsimInfo["eGamma"] << " (N) " << multiplicity.at(2)
                                       << " Neutral Had (E) " << hsimInfo["eNeutralHad"] << " (N) "
                                       << multiplicity.at(3) << " Charged Had (E) " << hsimInfo["eChargedHad"]
                                       << " (N) " << multiplicity.at(4) << " Total (E) " << hsimInfo["eTotal"]
                                       << " (N) " << multiplicity.at(5) << " PDG " << hsimInfo["pdgMatched"]
                                       << " Total E " << hsim << " NHit " << nSimHits;
        }
      }
    }
  }

  tree_->Fill();
}

void IsolatedTracksHcalScale::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 561 of file IsolatedTracksHcalScale.cc.

References doMC_, TFileService::make(), nEventProc_, t_Bunch, t_conehmaxNearP, t_e11x11_20Sig, t_e15x15_20Sig, t_eECALDR, t_eECALDR_1, t_eECALDR_2, t_eHCALDR, t_eHCALDRHB, t_eMipDR, t_eMipDR_1, t_eMipDR_2, t_hCone, t_hConeHB, t_hsim, t_hsimInfoCharHad, t_hsimInfoMatched, t_hsimInfoNCharHad, t_hsimInfoNeutHad, t_hsimInfoNMatched, t_hsimInfoNNeutHad, t_hsimInfoNPhoton, t_hsimInfoNRest, t_hsimInfoNTotal, t_hsimInfoPdgMatched, t_hsimInfoPhoton, t_hsimInfoRest, t_hsimInfoTotal, t_Lumi, t_nSimHits, t_RunNo, t_trackEta, t_trackHcalEta, t_trackHcalPhi, t_trackP, t_trackPhi, t_trackPt, and tree_.

                                       {
  nEventProc_ = 0;
  edm::Service<TFileService> fs;

  tree_ = fs->make<TTree>("tree", "tree");
  tree_->SetAutoSave(10000);

  tree_->Branch("t_RunNo", &t_RunNo, "t_RunNo/I");
  tree_->Branch("t_Lumi", &t_Lumi, "t_Lumi/I");
  tree_->Branch("t_Bunch", &t_Bunch, "t_Bunch/I");

  t_trackP = new std::vector<double>();
  t_trackPt = new std::vector<double>();
  t_trackEta = new std::vector<double>();
  t_trackPhi = new std::vector<double>();
  t_trackHcalEta = new std::vector<double>();
  t_trackHcalPhi = new std::vector<double>();
  t_hCone = new std::vector<double>();
  t_conehmaxNearP = new std::vector<double>();
  t_eMipDR = new std::vector<double>();
  t_eECALDR = new std::vector<double>();
  t_eHCALDR = new std::vector<double>();
  t_e11x11_20Sig = new std::vector<double>();
  t_e15x15_20Sig = new std::vector<double>();
  t_eMipDR_1 = new std::vector<double>();
  t_eECALDR_1 = new std::vector<double>();
  t_eMipDR_2 = new std::vector<double>();
  t_eECALDR_2 = new std::vector<double>();
  t_hConeHB = new std::vector<double>();
  t_eHCALDRHB = new std::vector<double>();

  tree_->Branch("t_trackP", "std::vector<double>", &t_trackP);
  tree_->Branch("t_trackPt", "std::vector<double>", &t_trackPt);
  tree_->Branch("t_trackEta", "std::vector<double>", &t_trackEta);
  tree_->Branch("t_trackPhi", "std::vector<double>", &t_trackPhi);
  tree_->Branch("t_trackHcalEta", "std::vector<double>", &t_trackHcalEta);
  tree_->Branch("t_trackHcalPhi", "std::vector<double>", &t_trackHcalPhi);
  tree_->Branch("t_hCone", "std::vector<double>", &t_hCone);
  tree_->Branch("t_conehmaxNearP", "std::vector<double>", &t_conehmaxNearP);
  tree_->Branch("t_eMipDR", "std::vector<double>", &t_eMipDR);
  tree_->Branch("t_eECALDR", "std::vector<double>", &t_eECALDR);
  tree_->Branch("t_eHCALDR", "std::vector<double>", &t_eHCALDR);
  tree_->Branch("t_e11x11_20Sig", "std::vector<double>", &t_e11x11_20Sig);
  tree_->Branch("t_e15x15_20Sig", "std::vector<double>", &t_e15x15_20Sig);
  tree_->Branch("t_eMipDR_1", "std::vector<double>", &t_eMipDR_1);
  tree_->Branch("t_eECALDR_1", "std::vector<double>", &t_eECALDR_1);
  tree_->Branch("t_eMipDR_2", "std::vector<double>", &t_eMipDR_2);
  tree_->Branch("t_eECALDR_2", "std::vector<double>", &t_eECALDR_2);
  tree_->Branch("t_hConeHB", "std::vector<double>", &t_hConeHB);
  tree_->Branch("t_eHCALDRHB", "std::vector<double>", &t_eHCALDRHB);

  if (doMC_) {
    t_hsimInfoMatched = new std::vector<double>();
    t_hsimInfoRest = new std::vector<double>();
    t_hsimInfoPhoton = new std::vector<double>();
    t_hsimInfoNeutHad = new std::vector<double>();
    t_hsimInfoCharHad = new std::vector<double>();
    t_hsimInfoPdgMatched = new std::vector<double>();
    t_hsimInfoTotal = new std::vector<double>();
    t_hsimInfoNMatched = new std::vector<int>();
    t_hsimInfoNTotal = new std::vector<int>();
    t_hsimInfoNNeutHad = new std::vector<int>();
    t_hsimInfoNCharHad = new std::vector<int>();
    t_hsimInfoNPhoton = new std::vector<int>();
    t_hsimInfoNRest = new std::vector<int>();
    t_hsim = new std::vector<double>();
    t_nSimHits = new std::vector<int>();

    tree_->Branch("t_hsimInfoMatched", "std::vector<double>", &t_hsimInfoMatched);
    tree_->Branch("t_hsimInfoRest", "std::vector<double>", &t_hsimInfoRest);
    tree_->Branch("t_hsimInfoPhoton", "std::vector<double>", &t_hsimInfoPhoton);
    tree_->Branch("t_hsimInfoNeutHad", "std::vector<double>", &t_hsimInfoNeutHad);
    tree_->Branch("t_hsimInfoCharHad", "std::vector<double>", &t_hsimInfoCharHad);
    tree_->Branch("t_hsimInfoPdgMatched", "std::vector<double>", &t_hsimInfoPdgMatched);
    tree_->Branch("t_hsimInfoTotal", "std::vector<double>", &t_hsimInfoTotal);
    tree_->Branch("t_hsimInfoNMatched", "std::vector<int>", &t_hsimInfoNMatched);
    tree_->Branch("t_hsimInfoNTotal", "std::vector<int>", &t_hsimInfoNTotal);
    tree_->Branch("t_hsimInfoNNeutHad", "std::vector<int>", &t_hsimInfoNNeutHad);
    tree_->Branch("t_hsimInfoNCharHad", "std::vector<int>", &t_hsimInfoNCharHad);
    tree_->Branch("t_hsimInfoNPhoton", "std::vector<int>", &t_hsimInfoNPhoton);
    tree_->Branch("t_hsimInfoNRest", "std::vector<int>", &t_hsimInfoNRest);
    tree_->Branch("t_hsim", "std::vector<double>", &t_hsim);
    tree_->Branch("t_nSimHits", "std::vector<int>", &t_nSimHits);
  }
}

void IsolatedTracksHcalScale::clearTreeVectors ( )

private

Definition at line 650 of file IsolatedTracksHcalScale.cc.

References doMC_, t_conehmaxNearP, t_e11x11_20Sig, t_e15x15_20Sig, t_eECALDR, t_eECALDR_1, t_eECALDR_2, t_eHCALDR, t_eHCALDRHB, t_eMipDR, t_eMipDR_1, t_eMipDR_2, t_hCone, t_hConeHB, t_hsim, t_hsimInfoCharHad, t_hsimInfoMatched, t_hsimInfoNCharHad, t_hsimInfoNeutHad, t_hsimInfoNMatched, t_hsimInfoNNeutHad, t_hsimInfoNPhoton, t_hsimInfoNRest, t_hsimInfoNTotal, t_hsimInfoPdgMatched, t_hsimInfoPhoton, t_hsimInfoRest, t_hsimInfoTotal, t_nSimHits, t_trackEta, t_trackHcalEta, t_trackHcalPhi, t_trackP, t_trackPhi, and t_trackPt.

Referenced by analyze().

                                               {
  t_trackP->clear();
  t_trackPt->clear();
  t_trackEta->clear();
  t_trackPhi->clear();
  t_trackHcalEta->clear();
  t_trackHcalPhi->clear();
  t_hCone->clear();
  t_conehmaxNearP->clear();
  t_eMipDR->clear();
  t_eECALDR->clear();
  t_eHCALDR->clear();
  t_e11x11_20Sig->clear();
  t_e15x15_20Sig->clear();
  t_eMipDR_1->clear();
  t_eECALDR_1->clear();
  t_eMipDR_2->clear();
  t_eECALDR_2->clear();
  t_hConeHB->clear();
  t_eHCALDRHB->clear();

  if (doMC_) {
    t_hsimInfoMatched->clear();
    t_hsimInfoRest->clear();
    t_hsimInfoPhoton->clear();
    t_hsimInfoNeutHad->clear();
    t_hsimInfoCharHad->clear();
    t_hsimInfoPdgMatched->clear();
    t_hsimInfoTotal->clear();
    t_hsimInfoNMatched->clear();
    t_hsimInfoNTotal->clear();
    t_hsimInfoNNeutHad->clear();
    t_hsimInfoNCharHad->clear();
    t_hsimInfoNPhoton->clear();
    t_hsimInfoNRest->clear();
    t_hsim->clear();
    t_nSimHits->clear();
  }
}

void IsolatedTracksHcalScale::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 648 of file IsolatedTracksHcalScale.cc.

References nEventProc_.

Referenced by o2olib.O2ORunMgr::executeJob().

{ edm::LogVerbatim("IsoTrack") << "Number of Events Processed " << nEventProc_; }

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

static

Definition at line 222 of file IsolatedTracksHcalScale.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::addUntracked(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                         {
  edm::ParameterSetDescription desc;
  desc.addUntracked<bool>("doMC", false);
  desc.addUntracked<int>("Verbosity", 0);
  desc.addUntracked<std::string>("TrackQuality", "highPurity");
  desc.addUntracked<double>("MinTrackPt", 10.0);
  desc.addUntracked<double>("MaxDxyPV", 0.02);
  desc.addUntracked<double>("MaxDzPV", 0.02);
  desc.addUntracked<double>("MaxChi2", 5.0);
  desc.addUntracked<double>("MaxDpOverP", 0.1);
  desc.addUntracked<int>("MinOuterHit", 4);
  desc.addUntracked<int>("MinLayerCrossed", 8);
  desc.addUntracked<int>("MaxInMiss", 0);
  desc.addUntracked<int>("MaxOutMiss", 0);
  desc.addUntracked<double>("ConeRadius", 34.98);
  desc.addUntracked<double>("ConeRadiusMIP", 14.0);
  desc.addUntracked<double>("TimeMinCutECAL", -500.0);
  desc.addUntracked<double>("TimeMaxCutECAL", 500.0);
  descriptions.add("isolatedTracksHcalScale", desc);
}

Member Data Documentation

double IsolatedTracksHcalScale::a_charIsoR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

double IsolatedTracksHcalScale::a_coneR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

double IsolatedTracksHcalScale::a_mipR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

double IsolatedTracksHcalScale::a_neutIsoR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

bool IsolatedTracksHcalScale::doMC_

private

Definition at line 114 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), beginJob(), clearTreeVectors(), and IsolatedTracksHcalScale().

std::string IsolatedTracksHcalScale::minQuality_

private

Definition at line 116 of file IsolatedTracksHcalScale.cc.

int IsolatedTracksHcalScale::myverbose_

private

Definition at line 115 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

int IsolatedTracksHcalScale::nEventProc_

private

Definition at line 142 of file IsolatedTracksHcalScale.cc.

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

spr::trackSelectionParameters IsolatedTracksHcalScale::selectionParameters_

private

Definition at line 117 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

int IsolatedTracksHcalScale::t_Bunch

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsolatedTracksHcalScale::t_conehmaxNearP

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 150 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 150 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eECALDR

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eECALDR_1

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eECALDR_2

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eHCALDR

private

Definition at line 148 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eHCALDRHB

private

Definition at line 152 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eMipDR

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

std::vector<double>* IsolatedTracksHcalScale::t_eMipDR_1

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_eMipDR_2

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

int IsolatedTracksHcalScale::t_EvtNo

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze().

std::vector<double>* IsolatedTracksHcalScale::t_hCone

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

std::vector<double>* IsolatedTracksHcalScale::t_hConeHB

private

Definition at line 152 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_hsim

private

Definition at line 155 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_hsimInfoCharHad

private

Definition at line 154 of file IsolatedTracksHcalScale.cc.

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

std::vector<double>* IsolatedTracksHcalScale::t_hsimInfoMatched

private

Definition at line 153 of file IsolatedTracksHcalScale.cc.

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

std::vector<int>* IsolatedTracksHcalScale::t_hsimInfoNCharHad

private

Definition at line 157 of file IsolatedTracksHcalScale.cc.

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

std::vector<double>* IsolatedTracksHcalScale::t_hsimInfoNeutHad

private

Definition at line 154 of file IsolatedTracksHcalScale.cc.

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

std::vector<int>* IsolatedTracksHcalScale::t_hsimInfoNMatched

private

Definition at line 156 of file IsolatedTracksHcalScale.cc.

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

std::vector<int> * IsolatedTracksHcalScale::t_hsimInfoNNeutHad

private

Definition at line 156 of file IsolatedTracksHcalScale.cc.

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

std::vector<int> * IsolatedTracksHcalScale::t_hsimInfoNPhoton

private

Definition at line 157 of file IsolatedTracksHcalScale.cc.

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

std::vector<int> * IsolatedTracksHcalScale::t_hsimInfoNRest

private

Definition at line 157 of file IsolatedTracksHcalScale.cc.

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

std::vector<int> * IsolatedTracksHcalScale::t_hsimInfoNTotal

private

Definition at line 156 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_hsimInfoPdgMatched

private

Definition at line 154 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_hsimInfoPhoton

private

Definition at line 153 of file IsolatedTracksHcalScale.cc.

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

std::vector<double> * IsolatedTracksHcalScale::t_hsimInfoRest

private

Definition at line 153 of file IsolatedTracksHcalScale.cc.

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

std::vector<double>* IsolatedTracksHcalScale::t_hsimInfoTotal

private

Definition at line 155 of file IsolatedTracksHcalScale.cc.

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

int IsolatedTracksHcalScale::t_Lumi

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().

std::vector<int>* IsolatedTracksHcalScale::t_nSimHits

private

Definition at line 158 of file IsolatedTracksHcalScale.cc.

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

int IsolatedTracksHcalScale::t_nTracks

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

int IsolatedTracksHcalScale::t_RunNo

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 148 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 148 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

std::string IsolatedTracksHcalScale::theTrackQuality_

private

Definition at line 116 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

double IsolatedTracksHcalScale::tMaxE_

private

Definition at line 119 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

double IsolatedTracksHcalScale::tMinE_

private

Definition at line 119 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 125 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksHcalScale::tok_caloEB_

private

Definition at line 132 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksHcalScale::tok_caloEE_

private

Definition at line 133 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<edm::PCaloHitContainer> IsolatedTracksHcalScale::tok_caloHH_

private

Definition at line 134 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 137 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksHcalScale::tok_EB_

private

Definition at line 126 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 139 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksHcalScale::tok_EE_

private

Definition at line 127 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 123 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 136 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<HBHERecHitCollection> IsolatedTracksHcalScale::tok_hbhe_

private

Definition at line 128 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 138 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 124 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

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

private

Definition at line 140 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<edm::SimTrackContainer> IsolatedTracksHcalScale::tok_simTk_

private

Definition at line 130 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

edm::EDGetTokenT<edm::SimVertexContainer> IsolatedTracksHcalScale::tok_simVtx_

private

Definition at line 131 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

TrackerHitAssociator::Config IsolatedTracksHcalScale::trackerHitAssociatorConfig_

private

Definition at line 121 of file IsolatedTracksHcalScale.cc.

Referenced by analyze().

TTree* IsolatedTracksHcalScale::tree_

private

Definition at line 144 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().