ValidIsoTrkCalib Class Reference

Inheritance diagram for ValidIsoTrkCalib:

Public Member Functions
	ValidIsoTrkCalib (const edm::ParameterSet &)
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)

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

Private Attributes
const bool	allowMissingInputs_
const double	associationConeSize_
const std::string	AxB_
const double	calibrationConeSize_
float	CentHitFactor
float	delR
int	dietatr
int	diphitr
float	e3x3After
float	e3x3Before
float	e5x5After
float	e5x5Before
float	eAfterDepth1
float	eAfterDepth2
float	eBeforeDepth1
float	eBeforeDepth2
float	eCentHitAfter
float	eCentHitBefore
float	eClustAfter
float	eClustBefore
float	eECAL
Float_t	emEnergy
const double	energyECALmip
const double	energyMaxIso
const double	energyMinIso
float	etaTrack
float	eTrack
int	eventNumber
TTree *	fTree
int	gen
std::vector< edm::InputTag >	genecalLabel_
float	genEta [500]
const edm::InputTag	genhbheLabel_
edm::InputTag	genhoLabel_
float	genPhi [500]
float	genPt [500]
const CaloGeometry *	geo
float	HTime [100]
float	iDr
int	iEta
int	iEtaTr
int	iPhi
int	iPhiTr
int	iso
float	isoEta [500]
float	isoPhi [500]
float	isoPt [500]
float	iTime
const double	maxPNear
const int	MinNTECHitsEndcap
const int	MinNTrackHitsBarrel
int	Nhits
int	NisoTrk
TH1F *	nTracks
int	numHits
float	numValidTrkHits
float	numValidTrkStrips
float	numVH
float	numVS
std::string	outputFileName_
TrackAssociatorParameters	parameters_
float	phiTrack
int	pix
float	pixEta [500]
float	pixPhi [500]
float	pixPt [500]
float	ptNear
float	PtNearBy
float	ptrack
const HcalRespCorrs *	respRecalib
int	runNumber
float	rvert
double	taECALCone_
double	taHCALCone_
const bool	takeGenTracks_
Float_t	targetE
const edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_
const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_
const edm::EDGetTokenT< HORecHitCollection >	tok_ho_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_recalibCorrs_
const edm::EDGetTokenT< reco::TrackCollection >	tok_track1_
const edm::EDGetTokenT< reco::IsolatedPixelTrackCandidateCollection >	tok_track_
TrackDetectorAssociator	trackAssociator_
float	trackE
float	trackEta
float	trackPhi
float	trackPt
TTree *	tTree
Float_t	xTrkEcal
Float_t	xTrkHcal
Float_t	yTrkEcal
Float_t	yTrkHcal
Float_t	zTrkEcal
Float_t	zTrkHcal

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 ()
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 64 of file ValidIsoTrkCalib.cc.

Constructor & Destructor Documentation

ValidIsoTrkCalib()

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

explicit

Definition at line 175 of file ValidIsoTrkCalib.cc.

References edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::getParameter(), TFileService::kSharedResource, TrackAssociatorParameters::loadParameters(), parameters_, trackAssociator_, and TrackDetectorAssociator::useDefaultPropagator().

    : takeGenTracks_(iConfig.getUntrackedParameter<bool>("takeGenTracks")),
      genhbheLabel_(iConfig.getParameter<edm::InputTag>("genHBHE")),
      associationConeSize_(iConfig.getParameter<double>("associationConeSize")),
      allowMissingInputs_(iConfig.getUntrackedParameter<bool>("allowMissingInputs", true)),
      AxB_(iConfig.getParameter<std::string>("AxB")),
      calibrationConeSize_(iConfig.getParameter<double>("calibrationConeSize")),
      MinNTrackHitsBarrel(iConfig.getParameter<int>("MinNTrackHitsBarrel")),
      MinNTECHitsEndcap(iConfig.getParameter<int>("MinNTECHitsEndcap")),
      energyECALmip(iConfig.getParameter<double>("energyECALmip")),
      energyMinIso(iConfig.getParameter<double>("energyMinIso")),
      energyMaxIso(iConfig.getParameter<double>("energyMaxIso")),
      maxPNear(iConfig.getParameter<double>("maxPNear")),
      tok_genTrack_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("genTracksLabel"))),
      tok_hbhe_(consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"))),
      tok_ho_(consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"))),
      tok_track_(consumes<reco::IsolatedPixelTrackCandidateCollection>(
          iConfig.getParameter<edm::InputTag>("HcalIsolTrackInput"))),
      tok_track1_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("trackInput"))),
      tok_recalibCorrs_(esConsumes(edm::ESInputTag("", "recalibrate"))),
      tok_geom_(esConsumes()) {
  usesResource(TFileService::kSharedResource);

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

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 205 of file ValidIsoTrkCalib.cc.

References funct::abs(), TrackDetectorAssociator::associate(), associationConeSize_, AxB_, edm::SortedCollection< T, SORT >::begin(), TrackAssociatorParameters::bFieldToken, calibrationConeSize_, CentHitFactor, delR, LEDCalibrationChannels::depth, MaxHit_struct::depthhit, dietatr, diphitr, MaxHit_struct::dr, l1ctLayer1_cff::dr, MillePedeFileConverter_cfg::e, e3x3After, e3x3Before, e5x5After, e5x5Before, eAfterDepth1, eAfterDepth2, eBeforeDepth1, eBeforeDepth2, eCentHitAfter, eCentHitBefore, eClustAfter, eClustBefore, eECAL, emEnergy, edm::SortedCollection< T, SORT >::end(), energyMaxIso, energyMinIso, etaTrack, eTrack, eventNumber, DQMScaleToClient_cfi::factor, fTree, generalTracks_cfi::generalTracks, genEta, genPhi, genPt, geo, HcalGeometry::getClosestCell(), edm::EventSetup::getData(), getDistInPlaneSimple(), TrackDetectorAssociator::getFreeTrajectoryState(), CaloGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), photonIsolationHIProducer_cfi::hbhe, DetId::Hcal, HcalBarrel, MaxHit_struct::hitenergy, HTime, mps_fire::i, iDr, LEDCalibrationChannels::ieta, iEta, HcalDetId::ieta(), MaxHit_struct::ietahitm, iEtaTr, iEvent, info(), LEDCalibrationChannels::iphi, iPhi, HcalDetId::iphi(), MaxHit_struct::iphihitm, iPhiTr, iso, isoEta, isoPhi, isoPt, iTime, muonTagProbeFilters_cff::matched, maxPNear, MinNTECHitsEndcap, MinNTrackHitsBarrel, dqmiodumpmetadata::n, nTracks, numHits, numValidTrkHits, numValidTrkStrips, numVH, numVS, parameters_, phiTrack, pix, pixEta, pixPhi, pixPt, ptNear, PtNearBy, respRecalib, runNumber, mathSSE::sqrt(), takeGenTracks_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_recalibCorrs_, tok_track1_, tok_track_, trackAssociator_, trackE, trackEta, trackPhi, trackPt, tTree, TrackAssociatorParameters::useCalo, TrackAssociatorParameters::useEcal, TrackAssociatorParameters::useHcal, TrackAssociatorParameters::useMuon, xTrkEcal, xTrkHcal, yTrkEcal, yTrkHcal, zTrkEcal, and zTrkHcal.

                                                                                  {
  try {
    respRecalib = &iSetup.getData(tok_recalibCorrs_);

    edm::LogVerbatim("ValidIsoTrkCalib") << "  Loaded:  OK ";

  } catch (const cms::Exception& e) {
    edm::LogWarning("ValidIsoTrkCalib") << "   Not Found!! ";
  }

  const edm::Handle<reco::TrackCollection>& generalTracks = iEvent.getHandle(tok_genTrack_);

  const edm::Handle<reco::TrackCollection>& isoProdTracks = iEvent.getHandle(tok_track1_);

  const edm::Handle<reco::IsolatedPixelTrackCandidateCollection>& isoPixelTracks = iEvent.getHandle(tok_track_);
  //edm::Handle<reco::TrackCollection> isoPixelTracks;

  /*
  edm::Handle<EcalRecHitCollection> ecal;
  iEvent.getByLabel(eLabel_,ecal);
  const EcalRecHitCollection Hitecal = *(ecal.product());
  */

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

  geo = &iSetup.getData(tok_geom_);

  const HcalGeometry* gHcal = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));
  //Note: even though it says HcalBarrel, we actually get the whole Hcal detector geometry!

  // Lumi_n=iEvent.luminosityBlock();
  parameters_.useEcal = true;
  parameters_.useHcal = true;
  parameters_.useCalo = false;
  parameters_.useMuon = false;
  //parameters_.dREcal = taECALCone_;
  //parameters_.dRHcal = taHCALCone_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("ValidIsoTrkCalib") << "Hello World. TrackCollectionSize: " << isoPixelTracks->size();
#endif
  if (isoPixelTracks->empty())
    return;

  for (reco::TrackCollection::const_iterator trit = isoProdTracks->begin(); trit != isoProdTracks->end(); trit++) {
    reco::IsolatedPixelTrackCandidateCollection::const_iterator isoMatched = isoPixelTracks->begin();
    //reco::TrackCollection::const_iterator isoMatched=isoPixelTracks->begin();
    bool matched = false;

    //for (reco::IsolatedPixelTrackCandidateCollection::const_iterator trit = isoPixelTracks->begin(); trit!=isoPixelTracks->end(); trit++)
    for (reco::IsolatedPixelTrackCandidateCollection::const_iterator it = isoPixelTracks->begin();
         it != isoPixelTracks->end();
         it++)
    //for (reco::TrackCollection::const_iterator it = isoPixelTracks->begin(); it!=isoPixelTracks->end(); it++)
    {
      if (abs((trit->pt() - it->pt()) / it->pt()) < 0.005 && abs(trit->eta() - it->eta()) < 0.01) {
        isoMatched = it;
        matched = true;
        break;
      }
    }
    // CUT

    if (!matched)
      continue;
    if (isoMatched->maxPtPxl() > maxPNear)
      continue;

    ptNear = isoMatched->maxPtPxl();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 0.1  isoMatch. ptnear: " << ptNear;
#endif
    // CUT
    if (trit->hitPattern().numberOfValidHits() < MinNTrackHitsBarrel)
      continue;
    if (fabs(trit->eta()) > 1.47 && trit->hitPattern().numberOfValidStripTECHits() < MinNTECHitsEndcap)
      continue;

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 0.2.1 after numofvalidhits HB: "
                                         << trit->hitPattern().numberOfValidHits();
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 0.2.2 after numofvalidstrips HE: "
                                         << trit->hitPattern().numberOfValidStripTECHits();
#endif
    numVH = trit->hitPattern().numberOfValidHits();
    numVS = trit->hitPattern().numberOfValidStripTECHits();

    trackE = sqrt(trit->px() * trit->px() + trit->py() * trit->py() + trit->pz() * trit->pz() + 0.14 * 0.14);
    trackPt = trit->pt();
    trackEta = trit->eta();
    trackPhi = trit->phi();

    emEnergy = isoMatched->energyIn();

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 0.3.  Matched :: pt: " << trit->pt() << " wholeEnergy: " << trackE
                                         << "  emEnergy: " << emEnergy << "  eta: " << trackEta << " phi: " << trackPhi;
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 0.4.  EM energy in cone: " << emEnergy;
#endif
    TrackDetMatchInfo info = trackAssociator_.associate(
        iEvent,
        iSetup,
        trackAssociator_.getFreeTrajectoryState(&iSetup.getData(parameters_.bFieldToken), *trit),
        parameters_);

    //float etaecal=info.trkGlobPosAtEcal.eta();
    //float phiecal=info.trkGlobPosAtEcal.phi();
    //  float etahcal=info.trkGlobPosAtHcal.eta();
    // float phihcal=info.trkGlobPosAtHcal.phi();

    xTrkEcal = info.trkGlobPosAtEcal.x();
    yTrkEcal = info.trkGlobPosAtEcal.y();
    zTrkEcal = info.trkGlobPosAtEcal.z();

    xTrkHcal = info.trkGlobPosAtHcal.x();
    yTrkHcal = info.trkGlobPosAtHcal.y();
    zTrkHcal = info.trkGlobPosAtHcal.z();

    if (xTrkEcal == 0 && yTrkEcal == 0 && zTrkEcal == 0) {
      edm::LogVerbatim("ValidIsoTrkCalib") << "zero point at Ecal";
      continue;
    }
    if (xTrkHcal == 0 && yTrkHcal == 0 && zTrkHcal == 0) {
      edm::LogVerbatim("ValidIsoTrkCalib") << "zero point at Hcal";
      continue;
    }

    /*GlobalVector trackMomAtEcal = info.trkMomAtEcal;
      GlobalVector trackMomAtHcal = info.trkMomAtHcal;
    
      PxTrkHcal = trackMomAtHcal.x();
      PyTrkHcal = trackMomAtHcal.y();
      PzTrkHcal = trackMomAtHcal.z();
      */

    GlobalPoint gPointEcal(xTrkEcal, yTrkEcal, zTrkEcal);
    GlobalPoint gPointHcal(xTrkHcal, yTrkHcal, zTrkHcal);

    int iphitrue = -10;
    int ietatrue = 100;
    const HcalDetId tempId = gHcal->getClosestCell(gPointHcal);
    ietatrue = tempId.ieta();
    iphitrue = tempId.iphi();

    MaxHit_struct MaxHit;

    MaxHit.hitenergy = -100.;

    //container for used recHits
    std::vector<DetId> usedHits;
    //
    usedHits.clear();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 1. Entrance to HBHECollection";
#endif
    //float dddeta = 1000.;
    //float dddphi = 1000.;
    //int iphitrue = 1234;
    //int ietatrue = 1234;

    GlobalPoint gPhot;

    for (HBHERecHitCollection::const_iterator hhit = Hithbhe.begin(); hhit != Hithbhe.end(); hhit++) {
      //check that this hit was not considered before and push it into usedHits
      bool hitIsUsed = false;
      for (uint32_t i = 0; i < usedHits.size(); i++) {
        if (usedHits[i] == hhit->id())
          hitIsUsed = true;
      }
      if (hitIsUsed)
        continue;
      usedHits.push_back(hhit->id());
      //

      // rof 16.05.2008 start: include the possibility for recalibration
      float recal = 1;
      // rof end

      GlobalPoint pos = geo->getPosition(hhit->detid());
      //float phihit = pos.phi();
      //float etahit = pos.eta();

      int iphihitm = (hhit->id()).iphi();
      int ietahitm = (hhit->id()).ieta();
      int depthhit = (hhit->id()).depth();
      float enehit = hhit->energy() * recal;

      if (depthhit != 1)
        continue;

#ifdef EDM_ML_DEBUG
      float dphi = fabs(info.trkGlobPosAtHcal.phi() - pos.phi());
      if (dphi > 4. * atan(1.))
        dphi = 8. * atan(1.) - dphi;
      float deta = fabs(info.trkGlobPosAtHcal.eta() - pos.eta());
      float dr = sqrt(dphi * dphi + deta * deta);
#endif

      //double distAtHcal =  getDistInPlaneTrackDir(gPointHcal, trackMomAtHcal, pos);
      double distAtHcal = getDistInPlaneSimple(gPointHcal, pos);

      if (distAtHcal < associationConeSize_) {
        for (HBHERecHitCollection::const_iterator hhit2 = Hithbhe.begin(); hhit2 != Hithbhe.end(); hhit2++) {
          int iphihitm2 = (hhit2->id()).iphi();
          int ietahitm2 = (hhit2->id()).ieta();
          int depthhit2 = (hhit2->id()).depth();
          float enehit2 = hhit2->energy() * recal;

          if (iphihitm == iphihitm2 && ietahitm == ietahitm2 && depthhit != depthhit2)
            enehit = enehit + enehit2;
        }

#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("ValidIsoTrkCalib") << "IN CONE ieta: " << ietahitm << "  iphi: " << iphihitm
                                             << " depthhit: " << depthhit << "  dr: " << dr << " energy: " << enehit;
#endif
        //Find a Hit with Maximum Energy

        if (enehit > MaxHit.hitenergy) {
          MaxHit.hitenergy = enehit;
          MaxHit.ietahitm = (hhit->id()).ieta();
          MaxHit.iphihitm = (hhit->id()).iphi();
          MaxHit.dr = distAtHcal;
          //MaxHit.depthhit  = (hhit->id()).depth();
          MaxHit.depthhit = 1;

          //gPhot = geo->getPosition(hhit->detid());
        }
      }
    }  //end of all HBHE hits cycle

    usedHits.clear();

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ValidIsoTrkCalib") << "Hottest ieta: " << MaxHit.ietahitm << "  iphi: " << MaxHit.iphihitm
                                         << "  dr: " << MaxHit.dr;
    edm::LogVerbatim("ValidIsoTrkCalib") << "Track   ieta: " << ietatrue << "  iphi: " << iphitrue;
    edm::LogVerbatim("ValidIsoTrkCalib") << "Point 3.  MaxHit :::En " << MaxHit.hitenergy
                                         << "  ieta: " << MaxHit.ietahitm << "  iphi: " << MaxHit.iphihitm;
#endif

    Bool_t passCuts = kFALSE;
    if (trackE > energyMinIso && trackE < energyMaxIso && emEnergy < energyECALmip && MaxHit.hitenergy > 0. &&
        abs(MaxHit.ietahitm) < 29)
      passCuts = kTRUE;

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ValidIsoTrkCalib") << "Pont 0.1.1.  trackE:" << trackE << "  emEn: " << emEnergy;
#endif
    numHits = 0;

    eClustBefore = 0.0;
    eClustAfter = 0.0;
    eBeforeDepth1 = 0.0;
    eAfterDepth1 = 0.0;
    eBeforeDepth2 = 0.0;
    eAfterDepth2 = 0.0;
    CentHitFactor = 0.0;
    e3x3After = 0.0;
    e3x3Before = 0.0;
    e5x5After = 0.0;
    e5x5Before = 0.0;

    for (HBHERecHitCollection::const_iterator hhit = Hithbhe.begin(); hhit != Hithbhe.end(); hhit++) {
      //check that this hit was not considered before and push it into usedHits
      bool hitIsUsed = false;
      for (uint32_t i = 0; i < usedHits.size(); i++) {
        if (usedHits[i] == hhit->id())
          hitIsUsed = true;
      }
      if (hitIsUsed)
        continue;
      usedHits.push_back(hhit->id());

      int DIETA = 100;
      if (MaxHit.ietahitm * (hhit->id()).ieta() > 0) {
        DIETA = MaxHit.ietahitm - (hhit->id()).ieta();
      }
      if (MaxHit.ietahitm * (hhit->id()).ieta() < 0) {
        DIETA = MaxHit.ietahitm - (hhit->id()).ieta();
        DIETA = DIETA > 0 ? DIETA - 1 : DIETA + 1;
      }

      int DIPHI = abs(MaxHit.iphihitm - (hhit->id()).iphi());
      DIPHI = DIPHI > 36 ? 72 - DIPHI : DIPHI;

      int numbercell = 100;  //always collect Wide clastor!

      //if(AxB_=="5x5" || AxB_=="3x3" || AxB_=="7x7"|| AxB_=="Cone")

      //if(AxB_=="3x3") numbercell = 1;
      //if(AxB_=="5x5") numbercell = 2;
      //if(AxB_=="Cone") numbercell = 1000;

      if (abs(DIETA) <= numbercell &&
          (abs(DIPHI) <= numbercell || (abs(MaxHit.ietahitm) >= 20 && abs(DIPHI) <= numbercell + 1))) {
        const GlobalPoint pos2 = geo->getPosition(hhit->detid());

        if (passCuts && hhit->energy() > 0) {
          float factor = 0.0;
          // factor = CalibFactors[hhit->id()];
          factor = respRecalib->getValues(hhit->id())->getValue();

#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("ValidIsoTrkCalib") << " calib factors: " << factor;
#endif

          if (hhit->id().ieta() == MaxHit.ietahitm && hhit->id().iphi() == MaxHit.iphihitm)
            CentHitFactor = factor;

          if (hhit->id().ieta() == MaxHit.ietahitm && hhit->id().iphi() == MaxHit.iphihitm)
            iTime = hhit->time();

          if (AxB_ != "3x3" && AxB_ != "5x5" && AxB_ != "Cone")
            edm::LogWarning(" AxB ") << "   Not supported: " << AxB_;

          if (abs(DIETA) <= 2 && (abs(DIPHI) <= 2 || ((abs(MaxHit.ietahitm) > 20 && abs(DIPHI) <= 4) &&
                                                      !((abs(MaxHit.ietahitm) == 21 || abs(MaxHit.ietahitm) == 22) &&
                                                        abs((hhit->id()).ieta()) <= 20 && abs(DIPHI) > 2)))) {
            e5x5Before += hhit->energy();
            e5x5After += hhit->energy() * factor;
          }

          if (abs(DIETA) <= 1 && (abs(DIPHI) <= 1 || ((abs(MaxHit.ietahitm) > 20 && abs(DIPHI) <= 2) &&
                                                      !(abs(MaxHit.ietahitm) == 21 && abs((hhit->id()).ieta()) <= 20 &&
                                                        abs(DIPHI) > 1)))) {
            e3x3Before += hhit->energy();
            e3x3After += hhit->energy() * factor;
          }

          if (AxB_ == "5x5") {
            if (abs(DIETA) <= 2 && (abs(DIPHI) <= 2 || (abs(MaxHit.ietahitm) > 20 && abs(DIPHI) <= 4))) {
              if (abs(MaxHit.ietahitm) == 21 && abs((hhit->id()).ieta()) <= 20 && abs(DIPHI) > 3)
                continue;

              HTime[numHits] = hhit->time();
              numHits++;

              eClustBefore += hhit->energy();
              eClustAfter += hhit->energy() * factor;

              if ((hhit->id().depth() == 1) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29)) {
                eBeforeDepth1 += hhit->energy();
                eAfterDepth1 += hhit->energy() * factor;
              } else if ((hhit->id().depth() == 2) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29)) {
                eBeforeDepth2 += hhit->energy();
                eAfterDepth2 += hhit->energy() * factor;
              }
            }
          }  //end of 5x5

          if (AxB_ == "3x3") {
            if (abs(DIETA) <= 1 && (abs(DIPHI) <= 1 || (abs(MaxHit.ietahitm) > 20 && abs(DIPHI) <= 2))) {
              if (abs(MaxHit.ietahitm) == 21 && abs((hhit->id()).ieta()) <= 20 && abs(DIPHI) > 2)
                continue;

              HTime[numHits] = hhit->time();
              numHits++;

              eClustBefore += hhit->energy();
              eClustAfter += hhit->energy() * factor;

              if ((hhit->id().depth() == 1) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29)) {
                eBeforeDepth1 += hhit->energy();
                eAfterDepth1 += hhit->energy() * factor;
              } else if ((hhit->id().depth() == 2) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29)) {
                eBeforeDepth2 += hhit->energy();
                eAfterDepth2 += hhit->energy() * factor;
              }
            }
          }  //end of 3x3

          if (AxB_ == "Cone" && getDistInPlaneSimple(gPointHcal, pos2) < calibrationConeSize_) {
            HTime[numHits] = hhit->time();
            numHits++;

            eClustBefore += hhit->energy();
            eClustAfter += hhit->energy() * factor;

            if ((hhit->id().depth() == 1) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29)) {
              eBeforeDepth1 += hhit->energy();
              eAfterDepth1 += hhit->energy() * factor;
            } else if ((hhit->id().depth() == 2) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29)) {
              eBeforeDepth2 += hhit->energy();
              eAfterDepth2 += hhit->energy() * factor;
            }

          }  //end of Cone

        }  //end of passCuts

      }  //end of DIETA DIPHI

    }  //end of associatedcone HBHE hits cycle

    int dieta_M_P = 100;
    int diphi_M_P = 100;
    if (MaxHit.ietahitm * ietatrue > 0) {
      dieta_M_P = abs(MaxHit.ietahitm - ietatrue);
    }
    if (MaxHit.ietahitm * ietatrue < 0) {
      dieta_M_P = abs(MaxHit.ietahitm - ietatrue) - 1;
    }
    diphi_M_P = abs(MaxHit.iphihitm - iphitrue);
    diphi_M_P = diphi_M_P > 36 ? 72 - diphi_M_P : diphi_M_P;

    if (passCuts)

    {
      eventNumber = iEvent.id().event();
      runNumber = iEvent.id().run();

      eCentHitBefore = MaxHit.hitenergy;
      eCentHitAfter = MaxHit.hitenergy * CentHitFactor;
      eECAL = emEnergy;
      numValidTrkHits = numVH;
      numValidTrkStrips = numVS;
      PtNearBy = ptNear;

      eTrack = trackE;
      phiTrack = trackPhi;
      etaTrack = trackEta;

      iEta = MaxHit.ietahitm;
      iPhi = MaxHit.iphihitm;

      iEtaTr = ietatrue;
      iPhiTr = iphitrue;
      iDr = sqrt(diphi_M_P * diphi_M_P + dieta_M_P * dieta_M_P);
      delR = MaxHit.dr;
      dietatr = dieta_M_P;
      diphitr = diphi_M_P;

      fTree->Fill();
    }

  }  //end of isoProdTracks cycle

  /* ------------------   Some stuff for general tracks  ----------   ----*/
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("ValidIsoTrkCalib") << " generalTracks Size: " << generalTracks->size();
#endif
  int n = generalTracks->size();
  nTracks->Fill(n);

  if (takeGenTracks_ && iEvent.id().event() % 10 == 1) {
    gen = generalTracks->size();
    iso = isoProdTracks->size();
    pix = isoPixelTracks->size();

    genPt[0] = -33;
    genPhi[0] = -33;
    genEta[0] = -33;

    isoPt[0] = -33;
    isoPhi[0] = -33;
    isoEta[0] = -33;

    pixPt[0] = -33;
    pixPhi[0] = -33;
    pixEta[0] = -33;

    Int_t gencount = 0, isocount = 0, pixcount = 0;
    for (reco::TrackCollection::const_iterator gentr = generalTracks->begin(); gentr != generalTracks->end(); gentr++) {
      genPt[gencount] = gentr->pt();
      genPhi[gencount] = gentr->phi();
      genEta[gencount] = gentr->eta();
      gencount++;
    }

    for (reco::TrackCollection::const_iterator isotr = isoProdTracks->begin(); isotr != isoProdTracks->end(); isotr++) {
      isoPt[isocount] = isotr->pt();
      isoPhi[isocount] = isotr->phi();
      isoEta[isocount] = isotr->eta();
      isocount++;
    }

    for (reco::IsolatedPixelTrackCandidateCollection::const_iterator pixtr = isoPixelTracks->begin();
         pixtr != isoPixelTracks->end();
         pixtr++) {
      pixPt[pixcount] = pixtr->pt();
      pixPhi[pixcount] = pixtr->phi();
      pixEta[pixcount] = pixtr->eta();
      pixcount++;
    }
  }

  tTree->Fill();
}

beginJob()

void ValidIsoTrkCalib::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 696 of file ValidIsoTrkCalib.cc.

References delR, dietatr, diphitr, e3x3After, e3x3Before, e5x5After, e5x5Before, eAfterDepth1, eAfterDepth2, eBeforeDepth1, eBeforeDepth2, eCentHitAfter, eCentHitBefore, eClustAfter, eClustBefore, eECAL, etaTrack, eTrack, eventNumber, compareTotals::fs, fTree, genEta, genPhi, genPt, HTime, iDr, iEta, iEtaTr, iPhi, iPhiTr, iso, isoEta, isoPhi, isoPt, iTime, nTracks, numHits, numValidTrkHits, numValidTrkStrips, phiTrack, pix, pixEta, pixPhi, pixPt, PtNearBy, runNumber, takeGenTracks_, tTree, xTrkEcal, xTrkHcal, yTrkEcal, yTrkHcal, zTrkEcal, and zTrkHcal.

                                {
  edm::Service<TFileService> fs;

  nTracks = fs->make<TH1F>("nTracks", "general;number of general tracks", 11, -0.5, 10.5);

  tTree = fs->make<TTree>("tTree", "Tree for gen info");

  fTree = fs->make<TTree>("fTree", "Tree for IsoTrack Calibration");

  fTree->Branch("eventNumber", &eventNumber, "eventNumber/I");
  fTree->Branch("runNumber", &runNumber, "runNumber/I");

  fTree->Branch("eClustBefore", &eClustBefore, "eClustBefore/F");
  fTree->Branch("eClustAfter", &eClustAfter, "eClustAfter/F");
  fTree->Branch("eTrack", &eTrack, "eTrack/F");
  fTree->Branch("etaTrack", &etaTrack, "etaTrack/F");
  fTree->Branch("phiTrack", &phiTrack, "phiTrack/F");

  fTree->Branch("numHits", &numHits, "numHits/I");
  fTree->Branch("eECAL", &eECAL, "eECAL/F");
  fTree->Branch("PtNearBy", &PtNearBy, "PtNearBy/F");
  fTree->Branch("numValidTrkHits", &numValidTrkHits, "numValidTrkHits/F");
  fTree->Branch("numValidTrkStrips", &numValidTrkStrips, "numValidTrkStrips/F");

  fTree->Branch("eBeforeDepth1", &eBeforeDepth1, "eBeforeDepth1/F");
  fTree->Branch("eBeforeDepth2", &eBeforeDepth2, "eBeforeDepth2/F");
  fTree->Branch("eAfterDepth1", &eAfterDepth1, "eAfterDepth1/F");
  fTree->Branch("eAfterDepth2", &eAfterDepth2, "eAfterDepth2/F");

  fTree->Branch("e3x3Before", &e3x3Before, "e3x3Before/F");
  fTree->Branch("e3x3After", &e3x3After, "e3x3After/F");
  fTree->Branch("e5x5Before", &e5x5Before, "e5x5Before/F");
  fTree->Branch("e5x5After", &e5x5After, "e5x5After/F");

  fTree->Branch("eCentHitAfter", &eCentHitAfter, "eCentHitAfter/F");
  fTree->Branch("eCentHitBefore", &eCentHitBefore, "eCentHitBefore/F");
  fTree->Branch("iEta", &iEta, "iEta/I");
  fTree->Branch("iPhi", &iPhi, "iPhi/I");

  fTree->Branch("iEtaTr", &iEtaTr, "iEtaTr/I");
  fTree->Branch("iPhiTr", &iPhiTr, "iPhiTr/I");
  fTree->Branch("dietatr", &dietatr, "dietatr/I");
  fTree->Branch("diphitr", &diphitr, "diphitr/I");
  fTree->Branch("iDr", &iDr, "iDr/F");
  fTree->Branch("delR", &delR, "delR/F");

  fTree->Branch("iTime", &iTime, "iTime/F");
  fTree->Branch("HTime", HTime, "HTime[numHits]/F");

  fTree->Branch("xTrkEcal", &xTrkEcal, "xTrkEcal/F");
  fTree->Branch("yTrkEcal", &yTrkEcal, "yTrkEcal/F");
  fTree->Branch("zTrkEcal", &zTrkEcal, "zTrkEcal/F");
  fTree->Branch("xTrkHcal", &xTrkHcal, "xTrkHcal/F");
  fTree->Branch("yTrkHcal", &yTrkHcal, "yTrkHcal/F");
  fTree->Branch("zTrkHcal", &zTrkHcal, "zTrkHcal/F");

  if (takeGenTracks_) {
    tTree->Branch("gen", &gen, "gen/I");
    tTree->Branch("iso", &iso, "iso/I");
    tTree->Branch("pix", &pix, "pix/I");
    tTree->Branch("genPt", genPt, "genPt[gen]/F");
    tTree->Branch("genPhi", genPhi, "genPhi[gen]/F");
    tTree->Branch("genEta", genEta, "genEta[gen]/F");

    tTree->Branch("isoPt", isoPt, "isoPt[iso]/F");
    tTree->Branch("isoPhi", isoPhi, "isoPhi[iso]/F");
    tTree->Branch("isoEta", isoEta, "isoEta[iso]/F");

    tTree->Branch("pixPt", pixPt, "pixPt[pix]/F");
    tTree->Branch("pixPhi", pixPhi, "pixPhi[pix]/F");
    tTree->Branch("pixEta", pixEta, "pixEta[pix]/F");
  }
}

Member Data Documentation

allowMissingInputs_

const bool ValidIsoTrkCalib::allowMissingInputs_

private

Definition at line 88 of file ValidIsoTrkCalib.cc.

associationConeSize_

const double ValidIsoTrkCalib::associationConeSize_

private

Definition at line 87 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

AxB_

const std::string ValidIsoTrkCalib::AxB_

private

Definition at line 89 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

calibrationConeSize_

const double ValidIsoTrkCalib::calibrationConeSize_

private

Definition at line 90 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

CentHitFactor

float ValidIsoTrkCalib::CentHitFactor

private

Definition at line 150 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

delR

float ValidIsoTrkCalib::delR

private

Definition at line 155 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

dietatr

int ValidIsoTrkCalib::dietatr

private

Definition at line 156 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

diphitr

int ValidIsoTrkCalib::diphitr

private

Definition at line 157 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e3x3After

float ValidIsoTrkCalib::e3x3After

private

Definition at line 162 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e3x3Before

float ValidIsoTrkCalib::e3x3Before

private

Definition at line 161 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e5x5After

float ValidIsoTrkCalib::e5x5After

private

Definition at line 164 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e5x5Before

float ValidIsoTrkCalib::e5x5Before

private

Definition at line 163 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eAfterDepth1

float ValidIsoTrkCalib::eAfterDepth1

private

Definition at line 145 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eAfterDepth2

float ValidIsoTrkCalib::eAfterDepth2

private

Definition at line 147 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eBeforeDepth1

float ValidIsoTrkCalib::eBeforeDepth1

private

Definition at line 144 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eBeforeDepth2

float ValidIsoTrkCalib::eBeforeDepth2

private

Definition at line 146 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eCentHitAfter

float ValidIsoTrkCalib::eCentHitAfter

private

Definition at line 149 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eCentHitBefore

float ValidIsoTrkCalib::eCentHitBefore

private

Definition at line 148 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eClustAfter

float ValidIsoTrkCalib::eClustAfter

private

Definition at line 137 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eClustBefore

float ValidIsoTrkCalib::eClustBefore

private

Definition at line 136 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eECAL

float ValidIsoTrkCalib::eECAL

private

Definition at line 141 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

emEnergy

Float_t ValidIsoTrkCalib::emEnergy

private

Definition at line 122 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

energyECALmip

const double ValidIsoTrkCalib::energyECALmip

private

Definition at line 93 of file ValidIsoTrkCalib.cc.

energyMaxIso

const double ValidIsoTrkCalib::energyMaxIso

private

Definition at line 95 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

energyMinIso

const double ValidIsoTrkCalib::energyMinIso

private

Definition at line 94 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

etaTrack

float ValidIsoTrkCalib::etaTrack

private

Definition at line 139 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eTrack

float ValidIsoTrkCalib::eTrack

private

Definition at line 138 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eventNumber

int ValidIsoTrkCalib::eventNumber

private

Definition at line 165 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

fTree

TTree * ValidIsoTrkCalib::fTree

private

Definition at line 125 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

gen

int ValidIsoTrkCalib::gen

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

genecalLabel_

std::vector<edm::InputTag> ValidIsoTrkCalib::genecalLabel_

private

Definition at line 108 of file ValidIsoTrkCalib.cc.

genEta

float ValidIsoTrkCalib::genEta[500]

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

genhbheLabel_

const edm::InputTag ValidIsoTrkCalib::genhbheLabel_

private

Definition at line 86 of file ValidIsoTrkCalib.cc.

genhoLabel_

edm::InputTag ValidIsoTrkCalib::genhoLabel_

private

Definition at line 107 of file ValidIsoTrkCalib.cc.

genPhi

float ValidIsoTrkCalib::genPhi[500]

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

genPt

float ValidIsoTrkCalib::genPt[500]

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

geo

const CaloGeometry* ValidIsoTrkCalib::geo

private

Definition at line 83 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

HTime

float ValidIsoTrkCalib::HTime[100]

private

Definition at line 160 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iDr

float ValidIsoTrkCalib::iDr

private

Definition at line 155 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iEta

int ValidIsoTrkCalib::iEta

private

Definition at line 151 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iEtaTr

int ValidIsoTrkCalib::iEtaTr

private

Definition at line 153 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iPhi

int ValidIsoTrkCalib::iPhi

private

Definition at line 152 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iPhiTr

int ValidIsoTrkCalib::iPhiTr

private

Definition at line 154 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iso

int ValidIsoTrkCalib::iso

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

isoEta

float ValidIsoTrkCalib::isoEta[500]

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

isoPhi

float ValidIsoTrkCalib::isoPhi[500]

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

isoPt

float ValidIsoTrkCalib::isoPt[500]

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iTime

float ValidIsoTrkCalib::iTime

private

Definition at line 159 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

maxPNear

const double ValidIsoTrkCalib::maxPNear

private

Definition at line 96 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

MinNTECHitsEndcap

const int ValidIsoTrkCalib::MinNTECHitsEndcap

private

Definition at line 92 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

MinNTrackHitsBarrel

const int ValidIsoTrkCalib::MinNTrackHitsBarrel

private

Definition at line 91 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

Nhits

int ValidIsoTrkCalib::Nhits

private

Definition at line 135 of file ValidIsoTrkCalib.cc.

NisoTrk

int ValidIsoTrkCalib::NisoTrk

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

nTracks

TH1F* ValidIsoTrkCalib::nTracks

private

Definition at line 172 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numHits

int ValidIsoTrkCalib::numHits

private

Definition at line 142 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numValidTrkHits

float ValidIsoTrkCalib::numValidTrkHits

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numValidTrkStrips

float ValidIsoTrkCalib::numValidTrkStrips

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numVH

float ValidIsoTrkCalib::numVH

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

numVS

float ValidIsoTrkCalib::numVS

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

outputFileName_

std::string ValidIsoTrkCalib::outputFileName_

private

Definition at line 109 of file ValidIsoTrkCalib.cc.

parameters_

TrackAssociatorParameters ValidIsoTrkCalib::parameters_

private

Definition at line 79 of file ValidIsoTrkCalib.cc.

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

phiTrack

float ValidIsoTrkCalib::phiTrack

private

Definition at line 140 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pix

int ValidIsoTrkCalib::pix

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pixEta

float ValidIsoTrkCalib::pixEta[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pixPhi

float ValidIsoTrkCalib::pixPhi[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pixPt

float ValidIsoTrkCalib::pixPt[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

ptNear

float ValidIsoTrkCalib::ptNear

private

Definition at line 118 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

PtNearBy

float ValidIsoTrkCalib::PtNearBy

private

Definition at line 167 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

ptrack

float ValidIsoTrkCalib::ptrack

private

Definition at line 119 of file ValidIsoTrkCalib.cc.

respRecalib

const HcalRespCorrs* ValidIsoTrkCalib::respRecalib

private

Definition at line 170 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

runNumber

int ValidIsoTrkCalib::runNumber

private

Definition at line 166 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), beginJob(), and TH2PolyOfflineMaps.TH2PolyOfflineMaps::PrintTrackerMaps().

rvert

float ValidIsoTrkCalib::rvert

private

Definition at line 119 of file ValidIsoTrkCalib.cc.

taECALCone_

double ValidIsoTrkCalib::taECALCone_

private

Definition at line 80 of file ValidIsoTrkCalib.cc.

taHCALCone_

double ValidIsoTrkCalib::taHCALCone_

private

Definition at line 81 of file ValidIsoTrkCalib.cc.

takeGenTracks_

const bool ValidIsoTrkCalib::takeGenTracks_

private

Definition at line 85 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

targetE

Float_t ValidIsoTrkCalib::targetE

private

Definition at line 123 of file ValidIsoTrkCalib.cc.

tok_genTrack_

const edm::EDGetTokenT<reco::TrackCollection> ValidIsoTrkCalib::tok_genTrack_

private

Definition at line 98 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

tok_geom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> ValidIsoTrkCalib::tok_geom_

private

Definition at line 105 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

tok_hbhe_

const edm::EDGetTokenT<HBHERecHitCollection> ValidIsoTrkCalib::tok_hbhe_

private

Definition at line 99 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

tok_ho_

const edm::EDGetTokenT<HORecHitCollection> ValidIsoTrkCalib::tok_ho_

private

Definition at line 100 of file ValidIsoTrkCalib.cc.

tok_recalibCorrs_

const edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> ValidIsoTrkCalib::tok_recalibCorrs_

private

Definition at line 104 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

tok_track1_

const edm::EDGetTokenT<reco::TrackCollection> ValidIsoTrkCalib::tok_track1_

private

Definition at line 102 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

tok_track_

const edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> ValidIsoTrkCalib::tok_track_

private

Definition at line 101 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

trackAssociator_

TrackDetectorAssociator ValidIsoTrkCalib::trackAssociator_

private

Definition at line 78 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

trackE

float ValidIsoTrkCalib::trackE

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

trackEta

float ValidIsoTrkCalib::trackEta

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

trackPhi

float ValidIsoTrkCalib::trackPhi

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

trackPt

float ValidIsoTrkCalib::trackPt

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

tTree

TTree* ValidIsoTrkCalib::tTree

private

Definition at line 125 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

xTrkEcal

Float_t ValidIsoTrkCalib::xTrkEcal

private

Definition at line 127 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

xTrkHcal

Float_t ValidIsoTrkCalib::xTrkHcal

private

Definition at line 131 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

yTrkEcal

Float_t ValidIsoTrkCalib::yTrkEcal

private

Definition at line 128 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

yTrkHcal

Float_t ValidIsoTrkCalib::yTrkHcal

private

Definition at line 132 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

zTrkEcal

Float_t ValidIsoTrkCalib::zTrkEcal

private

Definition at line 129 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

zTrkHcal

Float_t ValidIsoTrkCalib::zTrkHcal

private

Definition at line 133 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().