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 noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 99 of file IsolatedTracksHcalScale.cc.

Constructor & Destructor Documentation

IsolatedTracksHcalScale()

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(), ProducerED_cfi::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 StandaloneTrackMonitor_cfi::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::~IsolatedTracksHcalScale ( )

inlineoverride

Definition at line 102 of file IsolatedTracksHcalScale.cc.

{}

Member Function Documentation

analyze()

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_, Calorimetry_cff::bField, spr::chargeIsolationCone(), clearTreeVectors(), doMC_, spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), edm::EventSetup::getData(), spr::goodTrack(), photonIsolationHIProducer_cfi::hbhe, HcalBarrel, HcalDetId::ieta(), iEvent, HcalDetId::iphi(), edm::HandleBase::isValid(), reco::TrackBase::momentum(), myverbose_, nEventProc_, beamSpotPI::nTracks, reco::TrackBase::p(), reco::BeamSpot::position(), spr::propagateCALO(), reco::TrackBase::pt(), 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();
}

beginJob()

void IsolatedTracksHcalScale::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 561 of file IsolatedTracksHcalScale.cc.

References doMC_, compareTotals::fs, 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);
  }
}

clearTreeVectors()

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

endJob()

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

fillDescriptions()

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

static

Definition at line 222 of file IsolatedTracksHcalScale.cc.

References edm::ConfigurationDescriptions::add(), 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

a_charIsoR_

double IsolatedTracksHcalScale::a_charIsoR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

a_coneR_

double IsolatedTracksHcalScale::a_coneR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

a_mipR_

double IsolatedTracksHcalScale::a_mipR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

a_neutIsoR_

double IsolatedTracksHcalScale::a_neutIsoR_

private

Definition at line 118 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

doMC_

bool IsolatedTracksHcalScale::doMC_

private

Definition at line 114 of file IsolatedTracksHcalScale.cc.

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

minQuality_

std::string IsolatedTracksHcalScale::minQuality_

private

Definition at line 116 of file IsolatedTracksHcalScale.cc.

myverbose_

int IsolatedTracksHcalScale::myverbose_

private

Definition at line 115 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

nEventProc_

int IsolatedTracksHcalScale::nEventProc_

private

Definition at line 142 of file IsolatedTracksHcalScale.cc.

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

selectionParameters_

spr::trackSelectionParameters IsolatedTracksHcalScale::selectionParameters_

private

Definition at line 117 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

t_Bunch

int IsolatedTracksHcalScale::t_Bunch

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().

t_conehmaxNearP

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

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

t_e11x11_20Sig

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

private

Definition at line 150 of file IsolatedTracksHcalScale.cc.

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

t_e15x15_20Sig

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

private

Definition at line 150 of file IsolatedTracksHcalScale.cc.

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

t_eECALDR

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

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

t_eECALDR_1

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

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

t_eECALDR_2

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

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

t_eHCALDR

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

private

Definition at line 148 of file IsolatedTracksHcalScale.cc.

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

t_eHCALDRHB

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

private

Definition at line 152 of file IsolatedTracksHcalScale.cc.

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

t_eMipDR

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

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

t_eMipDR_1

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

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

t_eMipDR_2

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

private

Definition at line 151 of file IsolatedTracksHcalScale.cc.

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

t_EvtNo

int IsolatedTracksHcalScale::t_EvtNo

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze().

t_hCone

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

private

Definition at line 149 of file IsolatedTracksHcalScale.cc.

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

t_hConeHB

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

private

Definition at line 152 of file IsolatedTracksHcalScale.cc.

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

t_hsim

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

private

Definition at line 155 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoCharHad

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

private

Definition at line 154 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoMatched

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

private

Definition at line 153 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNCharHad

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

private

Definition at line 157 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNeutHad

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

private

Definition at line 154 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNMatched

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

private

Definition at line 156 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNNeutHad

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

private

Definition at line 156 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNPhoton

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

private

Definition at line 157 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNRest

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

private

Definition at line 157 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoNTotal

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

private

Definition at line 156 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoPdgMatched

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

private

Definition at line 154 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoPhoton

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

private

Definition at line 153 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoRest

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

private

Definition at line 153 of file IsolatedTracksHcalScale.cc.

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

t_hsimInfoTotal

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

private

Definition at line 155 of file IsolatedTracksHcalScale.cc.

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

t_Lumi

int IsolatedTracksHcalScale::t_Lumi

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().

t_nSimHits

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

private

Definition at line 158 of file IsolatedTracksHcalScale.cc.

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

t_nTracks

int IsolatedTracksHcalScale::t_nTracks

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

t_RunNo

int IsolatedTracksHcalScale::t_RunNo

private

Definition at line 146 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().

t_trackEta

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

t_trackHcalEta

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

private

Definition at line 148 of file IsolatedTracksHcalScale.cc.

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

t_trackHcalPhi

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

private

Definition at line 148 of file IsolatedTracksHcalScale.cc.

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

t_trackP

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

t_trackPhi

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

t_trackPt

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

private

Definition at line 147 of file IsolatedTracksHcalScale.cc.

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

theTrackQuality_

std::string IsolatedTracksHcalScale::theTrackQuality_

private

Definition at line 116 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tMaxE_

double IsolatedTracksHcalScale::tMaxE_

private

Definition at line 119 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tMinE_

double IsolatedTracksHcalScale::tMinE_

private

Definition at line 119 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_bs_

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

private

Definition at line 125 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_caloEB_

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

private

Definition at line 132 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_caloEE_

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

private

Definition at line 133 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_caloHH_

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

private

Definition at line 134 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_caloTopology_

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

private

Definition at line 137 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksHcalScale::tok_EB_

private

Definition at line 126 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_ecalChStatus_

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

private

Definition at line 139 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> IsolatedTracksHcalScale::tok_EE_

private

Definition at line 127 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_genTrack_

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

private

Definition at line 123 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_geom_

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

private

Definition at line 136 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> IsolatedTracksHcalScale::tok_hbhe_

private

Definition at line 128 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_magField_

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

private

Definition at line 138 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_recVtx_

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

private

Definition at line 124 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_sevlv_

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

private

Definition at line 140 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_simTk_

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

private

Definition at line 130 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

tok_simVtx_

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

private

Definition at line 131 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and IsolatedTracksHcalScale().

trackerHitAssociatorConfig_

TrackerHitAssociator::Config IsolatedTracksHcalScale::trackerHitAssociatorConfig_

private

Definition at line 121 of file IsolatedTracksHcalScale.cc.

Referenced by analyze().

tree_

TTree* IsolatedTracksHcalScale::tree_

private

Definition at line 144 of file IsolatedTracksHcalScale.cc.

Referenced by analyze(), and beginJob().