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
bool	allowMissingInputs_
double	associationConeSize_
string	AxB_
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
double	energyECALmip
double	energyMaxIso
double	energyMinIso
float	etaTrack
float	eTrack
int	eventNumber
edm::Service< TFileService >	fs
TTree *	fTree
int	gen
std::vector< edm::InputTag >	genecalLabel_
float	genEta [500]
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
double	maxPNear
int	MinNTECHitsEndcap
int	MinNTrackHitsBarrel
int	Nhits
int	NisoTrk
TH1F *	nTracks
int	numHits
float	numValidTrkHits
float	numValidTrkStrips
float	numVH
float	numVS
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_
bool	takeGenTracks_
Float_t	targetE
edm::EDGetTokenT < reco::TrackCollection >	tok_genTrack_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_
edm::EDGetTokenT < HORecHitCollection >	tok_ho_
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_recalibCorrs_
edm::EDGetTokenT < reco::TrackCollection >	tok_track1_
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 () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 64 of file ValidIsoTrkCalib.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 190 of file ValidIsoTrkCalib.cc.

References allowMissingInputs_, associationConeSize_, AxB_, calibrationConeSize_, edm::EDConsumerBase::consumesCollector(), energyECALmip, energyMaxIso, energyMinIso, edm::EDConsumerBase::esConsumes(), genhbheLabel_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), TFileService::kSharedResource, TrackAssociatorParameters::loadParameters(), maxPNear, MinNTECHitsEndcap, MinNTrackHitsBarrel, parameters_, AlCaHLTBitMon_QueryRunRegistry::string, takeGenTracks_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_ho_, tok_recalibCorrs_, tok_track1_, tok_track_, trackAssociator_, and TrackDetectorAssociator::useDefaultPropagator().

                                                                 {
  usesResource(TFileService::kSharedResource);

  //takeAllRecHits_=iConfig.getUntrackedParameter<bool>("takeAllRecHits");
  takeGenTracks_ = iConfig.getUntrackedParameter<bool>("takeGenTracks");

  tok_genTrack_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("genTracksLabel"));
  genhbheLabel_ = iConfig.getParameter<edm::InputTag>("genHBHE");
  //genhoLabel_=iConfig.getParameter<edm::InputTag>("genHO");
  //genecalLabel_=iConfig.getParameter<std::vector<edm::InputTag> >("genECAL");

  // m_hcalLabel = iConfig.getUntrackedParameter<std::string> ("hcalRecHitsLabel","hbhereco");

  tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
  tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
  //eLabel_=iConfig.getParameter<edm::InputTag>("eInput");
  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();

  associationConeSize_ = iConfig.getParameter<double>("associationConeSize");
  allowMissingInputs_ = iConfig.getUntrackedParameter<bool>("allowMissingInputs", true);
  //  outputFileName_=iConfig.getParameter<std::string>("outputFileName");
  //  calibFactorsFileName_=iConfig.getParameter<std::string>("calibFactorsFileName");

  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");

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

  // taECALCone_=iConfig.getUntrackedParameter<double>("TrackAssociatorECALCone",0.5);
  //taHCALCone_=iConfig.getUntrackedParameter<double>("TrackAssociatorHCALCone",0.6);
}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 238 of file ValidIsoTrkCalib.cc.

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

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

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

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

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

  edm::Handle<reco::TrackCollection> isoProdTracks;
  iEvent.getByToken(tok_track1_, isoProdTracks);

  edm::Handle<reco::IsolatedPixelTrackCandidateCollection> isoPixelTracks;
  //edm::Handle<reco::TrackCollection> isoPixelTracks;
  iEvent.getByToken(tok_track_, isoPixelTracks);

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

  edm::Handle<HBHERecHitCollection> hbhe;
  iEvent.getByToken(tok_hbhe_, 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("HcalIsoTrack") << "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("HcalIsoTrack") << "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("HcalIsoTrack") << "Point 0.2.1 after numofvalidhits HB: "
                                     << trit->hitPattern().numberOfValidHits();
    edm::LogVerbatim("HcalIsoTrack") << "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("HcalIsoTrack") << "Point 0.3.  Matched :: pt: " << trit->pt() << " wholeEnergy: " << trackE
                                     << "  emEnergy: " << emEnergy << "  eta: " << trackEta << " phi: " << trackPhi;
    edm::LogVerbatim("HcalIsoTrack") << "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("HcalIsoTrack") << "zero point at Ecal";
      continue;
    }
    if (xTrkHcal == 0 && yTrkHcal == 0 && zTrkHcal == 0) {
      edm::LogVerbatim("HcalIsoTrack") << "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("HcalIsoTrack") << "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("HcalIsoTrack") << "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("HcalIsoTrack") << "Hottest ieta: " << MaxHit.ietahitm << "  iphi: " << MaxHit.iphihitm
                                     << "  dr: " << MaxHit.dr;
    edm::LogVerbatim("HcalIsoTrack") << "Track   ieta: " << ietatrue << "  iphi: " << iphitrue;
    edm::LogVerbatim("HcalIsoTrack") << "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("HcalIsoTrack") << "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("HcalIsoTrack") << " 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("HcalIsoTrack") << " 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();
}

void ValidIsoTrkCalib::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 733 of file ValidIsoTrkCalib.cc.

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

                                {
#ifdef EDM_ML_DEBUG
//   if (!ReadCalibFactors(calibFactorsFileName_.c_str() ))
//     edm::LogVerbatim("HcalIsoTrack")<<"Cant read file with cailib coefficients!! ---";
#endif
  //  rootFile = new TFile(outputFileName_.c_str(),"RECREATE");

  //@@@@@@@@@@@@@
  //TFileDirectory ValDir = fs->mkdir("Validation");

  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

bool ValidIsoTrkCalib::allowMissingInputs_

private

Definition at line 105 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

double ValidIsoTrkCalib::associationConeSize_

private

Definition at line 101 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

string ValidIsoTrkCalib::AxB_

private

Definition at line 102 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

double ValidIsoTrkCalib::calibrationConeSize_

private

Definition at line 103 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

float ValidIsoTrkCalib::CentHitFactor

private

Definition at line 160 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::delR

private

Definition at line 165 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::dietatr

private

Definition at line 166 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::diphitr

private

Definition at line 167 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::e3x3After

private

Definition at line 172 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::e3x3Before

private

Definition at line 171 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::e5x5After

private

Definition at line 174 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::e5x5Before

private

Definition at line 173 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eAfterDepth1

private

Definition at line 155 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eAfterDepth2

private

Definition at line 157 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eBeforeDepth1

private

Definition at line 154 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eBeforeDepth2

private

Definition at line 156 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eCentHitAfter

private

Definition at line 159 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eCentHitBefore

private

Definition at line 158 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eClustAfter

private

Definition at line 146 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eClustBefore

private

Definition at line 145 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eECAL

private

Definition at line 150 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::emEnergy

private

Definition at line 129 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

double ValidIsoTrkCalib::energyECALmip

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

double ValidIsoTrkCalib::energyMaxIso

private

Definition at line 127 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

double ValidIsoTrkCalib::energyMinIso

private

Definition at line 127 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

float ValidIsoTrkCalib::etaTrack

private

Definition at line 148 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::eTrack

private

Definition at line 147 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::eventNumber

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

edm::Service<TFileService> ValidIsoTrkCalib::fs

private

Definition at line 186 of file ValidIsoTrkCalib.cc.

Referenced by beginJob().

TTree * ValidIsoTrkCalib::fTree

private

Definition at line 134 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::gen

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 87 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::genEta[500]

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

edm::InputTag ValidIsoTrkCalib::genhbheLabel_

private

Definition at line 85 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

edm::InputTag ValidIsoTrkCalib::genhoLabel_

private

Definition at line 86 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::genPhi[500]

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::genPt[500]

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

const CaloGeometry* ValidIsoTrkCalib::geo

private

Definition at line 83 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::HTime[100]

private

Definition at line 170 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::iDr

private

Definition at line 165 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::iEta

private

Definition at line 161 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::iEtaTr

private

Definition at line 163 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::iPhi

private

Definition at line 162 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::iPhiTr

private

Definition at line 164 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::iso

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::isoEta[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::isoPhi[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::isoPt[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::iTime

private

Definition at line 169 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

double ValidIsoTrkCalib::maxPNear

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

int ValidIsoTrkCalib::MinNTECHitsEndcap

private

Definition at line 125 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

int ValidIsoTrkCalib::MinNTrackHitsBarrel

private

Definition at line 125 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

int ValidIsoTrkCalib::Nhits

private

Definition at line 144 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::NisoTrk

private

Definition at line 119 of file ValidIsoTrkCalib.cc.

TH1F* ValidIsoTrkCalib::nTracks

private

Definition at line 184 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::numHits

private

Definition at line 151 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::numValidTrkHits

private

Definition at line 178 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::numValidTrkStrips

private

Definition at line 178 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::numVH

private

Definition at line 178 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::numVS

private

Definition at line 178 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

string ValidIsoTrkCalib::outputFileName_

private

Definition at line 106 of file ValidIsoTrkCalib.cc.

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().

float ValidIsoTrkCalib::phiTrack

private

Definition at line 149 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

int ValidIsoTrkCalib::pix

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::pixEta[500]

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::pixPhi[500]

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::pixPt[500]

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::ptNear

private

Definition at line 121 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::PtNearBy

private

Definition at line 177 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

float ValidIsoTrkCalib::ptrack

private

Definition at line 122 of file ValidIsoTrkCalib.cc.

const HcalRespCorrs* ValidIsoTrkCalib::respRecalib

private

Definition at line 180 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

int ValidIsoTrkCalib::runNumber

private

Definition at line 176 of file ValidIsoTrkCalib.cc.

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

float ValidIsoTrkCalib::rvert

private

Definition at line 122 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::taECALCone_

private

Definition at line 80 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::taHCALCone_

private

Definition at line 81 of file ValidIsoTrkCalib.cc.

bool ValidIsoTrkCalib::takeGenTracks_

private

Definition at line 110 of file ValidIsoTrkCalib.cc.

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

Float_t ValidIsoTrkCalib::targetE

private

Definition at line 130 of file ValidIsoTrkCalib.cc.

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

private

Definition at line 89 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

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

private

Definition at line 96 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

edm::EDGetTokenT<HBHERecHitCollection> ValidIsoTrkCalib::tok_hbhe_

private

Definition at line 90 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

edm::EDGetTokenT<HORecHitCollection> ValidIsoTrkCalib::tok_ho_

private

Definition at line 91 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

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

private

Definition at line 95 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

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

private

Definition at line 93 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

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

private

Definition at line 92 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

TrackDetectorAssociator ValidIsoTrkCalib::trackAssociator_

private

Definition at line 78 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

float ValidIsoTrkCalib::trackE

private

Definition at line 120 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::trackEta

private

Definition at line 120 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::trackPhi

private

Definition at line 120 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

float ValidIsoTrkCalib::trackPt

private

Definition at line 120 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

TTree* ValidIsoTrkCalib::tTree

private

Definition at line 134 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::xTrkEcal

private

Definition at line 136 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::xTrkHcal

private

Definition at line 140 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::yTrkEcal

private

Definition at line 137 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::yTrkHcal

private

Definition at line 141 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::zTrkEcal

private

Definition at line 138 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

Float_t ValidIsoTrkCalib::zTrkHcal

private

Definition at line 142 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().