ValidIsoTrkCalib Class Reference

Inheritance diagram for ValidIsoTrkCalib:

Public Member Functions
	ValidIsoTrkCalib (const edm::ParameterSet &)
	~ValidIsoTrkCalib () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
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
void	endJob () 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]
InputTag	genhbheLabel_
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::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_
edm::EDGetTokenT< HORecHitCollection >	tok_ho_
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::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 66 of file ValidIsoTrkCalib.cc.

Constructor & Destructor Documentation

ValidIsoTrkCalib()

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

explicit

Definition at line 191 of file ValidIsoTrkCalib.cc.

{
  //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"));
 
  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);
}

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

~ValidIsoTrkCalib()

ValidIsoTrkCalib::~ValidIsoTrkCalib ( )

override

Definition at line 235 of file ValidIsoTrkCalib.cc.

                                    {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 241 of file ValidIsoTrkCalib.cc.

                                                                                  {
  using namespace edm;
 
  try {
    edm::ESHandle<HcalRespCorrs> recalibCorrs;
    iSetup.get<HcalRespCorrsRcd>().get("recalibrate", recalibCorrs);
    respRecalib = recalibCorrs.product();
 
    LogInfo("CalibConstants") << "  Loaded:  OK ";
 
  } catch (const cms::Exception& e) {
    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());
 
  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  geo = pG.product();
 
  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_;
 
  //cout<<"Hello World. TrackCollectionSize: "<< isoPixelTracks->size()<<endl;
 
  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();
    //cout<<"Point 0.1  isoMatch. ptnear: "<<ptNear<<endl;
 
    // CUT
    if (trit->hitPattern().numberOfValidHits() < MinNTrackHitsBarrel)
      continue;
    if (fabs(trit->eta()) > 1.47 && trit->hitPattern().numberOfValidStripTECHits() < MinNTECHitsEndcap)
      continue;
 
    //cout<<"Point 0.2.1 after numofvalidhits HB: "<<trit->hitPattern().numberOfValidHits()<<endl;
    //cout<<"Point 0.2.2 after numofvalidstrips HE: "<<trit->hitPattern().numberOfValidStripTECHits()<<endl;
 
    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();
 
    //cout<<"Point 0.3.  Matched :: pt: "<<trit->pt()<<" wholeEnergy: "<<trackE<<"  emEnergy: "<<emEnergy<<"  eta: "<<etahcal<<" phi: "<<phihcal<<endl;
    //cout<<"Point 0.4.  EM energy in cone: "<<emEnergy<<"  EtaHcal: "<<etahcal<<"  PhiHcal: "<<phihcal<<endl;
 
    TrackDetMatchInfo info =
        trackAssociator_.associate(iEvent, iSetup, trackAssociator_.getFreeTrajectoryState(iSetup, *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) {
      cout << "zero point at Ecal" << endl;
      continue;
    }
    if (xTrkHcal == 0 && yTrkHcal == 0 && zTrkHcal == 0) {
      cout << "zero point at Hcal" << endl;
      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();
    //cout <<"Point 1. Entrance to HBHECollection"<<endl;
 
    //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;
 
      /*
        float dphi = fabs(phihcal - phihit); 
        if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
        float deta = fabs(etahcal - etahit); 
        float dr = sqrt(dphi*dphi + deta*deta);
      */
 
      //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;
        }
 
        //cout<<"IN CONE ieta: "<<ietahitm<<"  iphi: "<<iphihitm<<" depthhit: "<<depthhit<<"  dr: "<<dr<<" energy: "<<enehit<<endl;
 
        //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();
 
    //cout<<"Hottest ieta: "<<MaxHit.ietahitm<<"  iphi: "<<MaxHit.iphihitm<<"  dr: "<<MaxHit.dr<<endl;
    //cout<<"Track   ieta: "<<ietatrue<<"  iphi: "<<iphitrue<<endl;
 
    //cout<<"Point 3.  MaxHit :::En "<<MaxHit.hitenergy<<"  ieta: "<<MaxHit.ietahitm<<"  iphi: "<<MaxHit.iphihitm<<endl;
 
    Bool_t passCuts = kFALSE;
    if (trackE > energyMinIso && trackE < energyMaxIso && emEnergy < energyECALmip && MaxHit.hitenergy > 0. &&
        abs(MaxHit.ietahitm) < 29)
      passCuts = kTRUE;
 
    //cout<<"Pont 0.1.1.  trackE:"<<trackE<<"  emEn: "<<emEnergy<<endl;
 
    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();
 
          //if(i<5){cout<<" calib factors: "<<factor<<"  ij "<<100*i+j<<endl;}
 
          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")
            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  ----------   ----*/
  //cout<<" generalTracks Size: "<< generalTracks->size()<<endl;
  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();
}

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), gather_cfg::cout, LEDCalibrationChannels::depth, MaxHit_struct::depthhit, MaxHit_struct::dr, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), DQMScaleToClient_cfi::factor, isolatedTracks_cfi::generalTracks, PFJetFilter_cfi::genPt, edm::EventSetup::get(), edm::get(), HcalGeometry::getClosestCell(), getDistInPlaneSimple(), reco::JetExtendedAssociation::getValue(), photonIsolationHIProducer_cfi::hbhe, DetId::Hcal, HcalBarrel, MaxHit_struct::hitenergy, mps_fire::i, L1TowerCalibrationProducer_cfi::iEta, LEDCalibrationChannels::ieta, HcalDetId::ieta(), MaxHit_struct::ietahitm, iEvent, info(), LEDCalibrationChannels::iphi, HcalDetId::iphi(), MaxHit_struct::iphihitm, muonTagProbeFilters_cff::matched, dqmiodumpmetadata::n, nTracks(), edm::ESHandle< T >::product(), convertSQLiteXML::runNumber, mathSSE::sqrt(), reco::btau::trackEta, reco::btau::trackPhi, and listHistos::trackPt.

beginJob()

void ValidIsoTrkCalib::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 721 of file ValidIsoTrkCalib.cc.

                                {
  // if(!ReadCalibFactors(calibFactorsFileName_.c_str() )) {cout<<"Cant read file with cailib coefficients!! ---"<<endl;}
 
  // try{
  //   edm::ESHandle <HcalRespCorrs> recalibCorrs;
  //   iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
  //   respRecalib = recalibCorrs.product();
  //
  //   LogInfo("CalibConstants")<<"  Loaded:  OK ";
  //
  // }catch(const cms::Exception & e) {
  //   LogWarning("CalibConstants")<<"   Not Found!! ";
  // }
 
  //  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");
  }
}

References PFJetFilter_cfi::genPt, L1TowerCalibrationProducer_cfi::iEta, nTracks(), and convertSQLiteXML::runNumber.

endJob()

void ValidIsoTrkCalib::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 812 of file ValidIsoTrkCalib.cc.

                              {
  //  rootFile->Write();
  //rootFile->Close();
}

Referenced by o2olib.O2ORunMgr::executeJob().

Member Data Documentation

allowMissingInputs_

bool ValidIsoTrkCalib::allowMissingInputs_

private

Definition at line 106 of file ValidIsoTrkCalib.cc.

associationConeSize_

double ValidIsoTrkCalib::associationConeSize_

private

Definition at line 102 of file ValidIsoTrkCalib.cc.

AxB_

string ValidIsoTrkCalib::AxB_

private

Definition at line 103 of file ValidIsoTrkCalib.cc.

calibrationConeSize_

double ValidIsoTrkCalib::calibrationConeSize_

private

Definition at line 104 of file ValidIsoTrkCalib.cc.

CentHitFactor

float ValidIsoTrkCalib::CentHitFactor

private

Definition at line 161 of file ValidIsoTrkCalib.cc.

delR

float ValidIsoTrkCalib::delR

private

Definition at line 166 of file ValidIsoTrkCalib.cc.

dietatr

int ValidIsoTrkCalib::dietatr

private

Definition at line 167 of file ValidIsoTrkCalib.cc.

diphitr

int ValidIsoTrkCalib::diphitr

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

e3x3After

float ValidIsoTrkCalib::e3x3After

private

Definition at line 173 of file ValidIsoTrkCalib.cc.

e3x3Before

float ValidIsoTrkCalib::e3x3Before

private

Definition at line 172 of file ValidIsoTrkCalib.cc.

e5x5After

float ValidIsoTrkCalib::e5x5After

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

e5x5Before

float ValidIsoTrkCalib::e5x5Before

private

Definition at line 174 of file ValidIsoTrkCalib.cc.

eAfterDepth1

float ValidIsoTrkCalib::eAfterDepth1

private

Definition at line 156 of file ValidIsoTrkCalib.cc.

eAfterDepth2

float ValidIsoTrkCalib::eAfterDepth2

private

Definition at line 158 of file ValidIsoTrkCalib.cc.

eBeforeDepth1

float ValidIsoTrkCalib::eBeforeDepth1

private

Definition at line 155 of file ValidIsoTrkCalib.cc.

eBeforeDepth2

float ValidIsoTrkCalib::eBeforeDepth2

private

Definition at line 157 of file ValidIsoTrkCalib.cc.

eCentHitAfter

float ValidIsoTrkCalib::eCentHitAfter

private

Definition at line 160 of file ValidIsoTrkCalib.cc.

eCentHitBefore

float ValidIsoTrkCalib::eCentHitBefore

private

Definition at line 159 of file ValidIsoTrkCalib.cc.

eClustAfter

float ValidIsoTrkCalib::eClustAfter

private

Definition at line 147 of file ValidIsoTrkCalib.cc.

eClustBefore

float ValidIsoTrkCalib::eClustBefore

private

Definition at line 146 of file ValidIsoTrkCalib.cc.

eECAL

float ValidIsoTrkCalib::eECAL

private

Definition at line 151 of file ValidIsoTrkCalib.cc.

emEnergy

Float_t ValidIsoTrkCalib::emEnergy

private

Definition at line 130 of file ValidIsoTrkCalib.cc.

energyECALmip

double ValidIsoTrkCalib::energyECALmip

private

Definition at line 127 of file ValidIsoTrkCalib.cc.

energyMaxIso

double ValidIsoTrkCalib::energyMaxIso

private

Definition at line 128 of file ValidIsoTrkCalib.cc.

energyMinIso

double ValidIsoTrkCalib::energyMinIso

private

Definition at line 128 of file ValidIsoTrkCalib.cc.

etaTrack

float ValidIsoTrkCalib::etaTrack

private

Definition at line 149 of file ValidIsoTrkCalib.cc.

eTrack

float ValidIsoTrkCalib::eTrack

private

Definition at line 148 of file ValidIsoTrkCalib.cc.

eventNumber

int ValidIsoTrkCalib::eventNumber

private

Definition at line 176 of file ValidIsoTrkCalib.cc.

fs

edm::Service<TFileService> ValidIsoTrkCalib::fs

private

Definition at line 187 of file ValidIsoTrkCalib.cc.

fTree

TTree * ValidIsoTrkCalib::fTree

private

Definition at line 135 of file ValidIsoTrkCalib.cc.

gen

int ValidIsoTrkCalib::gen

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

genecalLabel_

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

private

Definition at line 91 of file ValidIsoTrkCalib.cc.

genEta

float ValidIsoTrkCalib::genEta[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

genhbheLabel_

InputTag ValidIsoTrkCalib::genhbheLabel_

private

Definition at line 89 of file ValidIsoTrkCalib.cc.

genhoLabel_

InputTag ValidIsoTrkCalib::genhoLabel_

private

Definition at line 90 of file ValidIsoTrkCalib.cc.

genPhi

float ValidIsoTrkCalib::genPhi[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

genPt

float ValidIsoTrkCalib::genPt[500]

private

Definition at line 114 of file ValidIsoTrkCalib.cc.

geo

const CaloGeometry* ValidIsoTrkCalib::geo

private

Definition at line 87 of file ValidIsoTrkCalib.cc.

HTime

float ValidIsoTrkCalib::HTime[100]

private

Definition at line 171 of file ValidIsoTrkCalib.cc.

iDr

float ValidIsoTrkCalib::iDr

private

Definition at line 166 of file ValidIsoTrkCalib.cc.

iEta

int ValidIsoTrkCalib::iEta

private

Definition at line 162 of file ValidIsoTrkCalib.cc.

iEtaTr

int ValidIsoTrkCalib::iEtaTr

private

Definition at line 164 of file ValidIsoTrkCalib.cc.

iPhi

int ValidIsoTrkCalib::iPhi

private

Definition at line 163 of file ValidIsoTrkCalib.cc.

iPhiTr

int ValidIsoTrkCalib::iPhiTr

private

Definition at line 165 of file ValidIsoTrkCalib.cc.

iso

int ValidIsoTrkCalib::iso

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

isoEta

float ValidIsoTrkCalib::isoEta[500]

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

isoPhi

float ValidIsoTrkCalib::isoPhi[500]

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

isoPt

float ValidIsoTrkCalib::isoPt[500]

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

iTime

float ValidIsoTrkCalib::iTime

private

Definition at line 170 of file ValidIsoTrkCalib.cc.

maxPNear

double ValidIsoTrkCalib::maxPNear

private

Definition at line 127 of file ValidIsoTrkCalib.cc.

MinNTECHitsEndcap

int ValidIsoTrkCalib::MinNTECHitsEndcap

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

MinNTrackHitsBarrel

int ValidIsoTrkCalib::MinNTrackHitsBarrel

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

Nhits

int ValidIsoTrkCalib::Nhits

private

Definition at line 145 of file ValidIsoTrkCalib.cc.

NisoTrk

int ValidIsoTrkCalib::NisoTrk

private

Definition at line 120 of file ValidIsoTrkCalib.cc.

nTracks

TH1F* ValidIsoTrkCalib::nTracks

private

Definition at line 185 of file ValidIsoTrkCalib.cc.

numHits

int ValidIsoTrkCalib::numHits

private

Definition at line 152 of file ValidIsoTrkCalib.cc.

numValidTrkHits

float ValidIsoTrkCalib::numValidTrkHits

private

Definition at line 179 of file ValidIsoTrkCalib.cc.

numValidTrkStrips

float ValidIsoTrkCalib::numValidTrkStrips

private

Definition at line 179 of file ValidIsoTrkCalib.cc.

numVH

float ValidIsoTrkCalib::numVH

private

Definition at line 179 of file ValidIsoTrkCalib.cc.

numVS

float ValidIsoTrkCalib::numVS

private

Definition at line 179 of file ValidIsoTrkCalib.cc.

outputFileName_

string ValidIsoTrkCalib::outputFileName_

private

Definition at line 107 of file ValidIsoTrkCalib.cc.

parameters_

TrackAssociatorParameters ValidIsoTrkCalib::parameters_

private

Definition at line 83 of file ValidIsoTrkCalib.cc.

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

phiTrack

float ValidIsoTrkCalib::phiTrack

private

Definition at line 150 of file ValidIsoTrkCalib.cc.

pix

int ValidIsoTrkCalib::pix

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

pixEta

float ValidIsoTrkCalib::pixEta[500]

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

pixPhi

float ValidIsoTrkCalib::pixPhi[500]

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

pixPt

float ValidIsoTrkCalib::pixPt[500]

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

ptNear

float ValidIsoTrkCalib::ptNear

private

Definition at line 122 of file ValidIsoTrkCalib.cc.

PtNearBy

float ValidIsoTrkCalib::PtNearBy

private

Definition at line 178 of file ValidIsoTrkCalib.cc.

ptrack

float ValidIsoTrkCalib::ptrack

private

Definition at line 123 of file ValidIsoTrkCalib.cc.

respRecalib

const HcalRespCorrs* ValidIsoTrkCalib::respRecalib

private

Definition at line 181 of file ValidIsoTrkCalib.cc.

runNumber

int ValidIsoTrkCalib::runNumber

private

Definition at line 177 of file ValidIsoTrkCalib.cc.

Referenced by TH2PolyOfflineMaps.TH2PolyOfflineMaps::PrintTrackerMaps().

rvert

float ValidIsoTrkCalib::rvert

private

Definition at line 123 of file ValidIsoTrkCalib.cc.

taECALCone_

double ValidIsoTrkCalib::taECALCone_

private

Definition at line 84 of file ValidIsoTrkCalib.cc.

taHCALCone_

double ValidIsoTrkCalib::taHCALCone_

private

Definition at line 85 of file ValidIsoTrkCalib.cc.

takeGenTracks_

bool ValidIsoTrkCalib::takeGenTracks_

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

targetE

Float_t ValidIsoTrkCalib::targetE

private

Definition at line 131 of file ValidIsoTrkCalib.cc.

tok_genTrack_

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

private

Definition at line 93 of file ValidIsoTrkCalib.cc.

tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> ValidIsoTrkCalib::tok_hbhe_

private

Definition at line 94 of file ValidIsoTrkCalib.cc.

tok_ho_

edm::EDGetTokenT<HORecHitCollection> ValidIsoTrkCalib::tok_ho_

private

Definition at line 95 of file ValidIsoTrkCalib.cc.

tok_track1_

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

private

Definition at line 97 of file ValidIsoTrkCalib.cc.

tok_track_

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

private

Definition at line 96 of file ValidIsoTrkCalib.cc.

trackAssociator_

TrackDetectorAssociator ValidIsoTrkCalib::trackAssociator_

private

Definition at line 82 of file ValidIsoTrkCalib.cc.

trackE

float ValidIsoTrkCalib::trackE

private

Definition at line 121 of file ValidIsoTrkCalib.cc.

trackEta

float ValidIsoTrkCalib::trackEta

private

Definition at line 121 of file ValidIsoTrkCalib.cc.

trackPhi

float ValidIsoTrkCalib::trackPhi

private

Definition at line 121 of file ValidIsoTrkCalib.cc.

trackPt

float ValidIsoTrkCalib::trackPt

private

Definition at line 121 of file ValidIsoTrkCalib.cc.

tTree

TTree* ValidIsoTrkCalib::tTree

private

Definition at line 135 of file ValidIsoTrkCalib.cc.

xTrkEcal

Float_t ValidIsoTrkCalib::xTrkEcal

private

Definition at line 137 of file ValidIsoTrkCalib.cc.

xTrkHcal

Float_t ValidIsoTrkCalib::xTrkHcal

private

Definition at line 141 of file ValidIsoTrkCalib.cc.

yTrkEcal

Float_t ValidIsoTrkCalib::yTrkEcal

private

Definition at line 138 of file ValidIsoTrkCalib.cc.

yTrkHcal

Float_t ValidIsoTrkCalib::yTrkHcal

private

Definition at line 142 of file ValidIsoTrkCalib.cc.

zTrkEcal

Float_t ValidIsoTrkCalib::zTrkEcal

private

Definition at line 139 of file ValidIsoTrkCalib.cc.

zTrkHcal

Float_t ValidIsoTrkCalib::zTrkHcal

private

Definition at line 143 of file ValidIsoTrkCalib.cc.