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
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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::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

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
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 = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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...
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 61 of file ValidIsoTrkCalib.cc.

Constructor & Destructor Documentation

ValidIsoTrkCalib()

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

explicit

Definition at line 187 of file ValidIsoTrkCalib.cc.

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

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

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 235 of file ValidIsoTrkCalib.cc.

                                                                                  {
  try {
    respRecalib = &iSetup.getData(tok_recalibCorrs_);
 
    edm::LogInfo("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_;
 
  //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.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) {
      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")
            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  ----------   ----*/
  //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(), TrackDetectorAssociator::associate(), associationConeSize_, AxB_, edm::SortedCollection< T, SORT >::begin(), TrackAssociatorParameters::bFieldToken, calibrationConeSize_, CentHitFactor, gather_cfg::cout, delR, LEDCalibrationChannels::depth, MaxHit_struct::depthhit, dietatr, diphitr, MaxHit_struct::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, isolatedTracks_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, HcalDetId::ieta(), iEta, MaxHit_struct::ietahitm, iEtaTr, iEvent, info(), LEDCalibrationChannels::iphi, HcalDetId::iphi(), 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.

beginJob()

void ValidIsoTrkCalib::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 712 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();
  //
  //   edm::LogInfo("CalibConstants")<<"  Loaded:  OK ";
  //
  // }catch(const cms::Exception & e) {
  //   edm::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 delR, dietatr, diphitr, e3x3After, e3x3Before, e5x5After, e5x5Before, eAfterDepth1, eAfterDepth2, eBeforeDepth1, eBeforeDepth2, eCentHitAfter, eCentHitBefore, eClustAfter, eClustBefore, eECAL, etaTrack, eTrack, eventNumber, fs, fTree, 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.

Member Data Documentation

allowMissingInputs_

bool ValidIsoTrkCalib::allowMissingInputs_

private

Definition at line 102 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

associationConeSize_

double ValidIsoTrkCalib::associationConeSize_

private

Definition at line 98 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

AxB_

string ValidIsoTrkCalib::AxB_

private

Definition at line 99 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

calibrationConeSize_

double ValidIsoTrkCalib::calibrationConeSize_

private

Definition at line 100 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

CentHitFactor

float ValidIsoTrkCalib::CentHitFactor

private

Definition at line 157 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

delR

float ValidIsoTrkCalib::delR

private

Definition at line 162 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

dietatr

int ValidIsoTrkCalib::dietatr

private

Definition at line 163 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

diphitr

int ValidIsoTrkCalib::diphitr

private

Definition at line 164 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e3x3After

float ValidIsoTrkCalib::e3x3After

private

Definition at line 169 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e3x3Before

float ValidIsoTrkCalib::e3x3Before

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e5x5After

float ValidIsoTrkCalib::e5x5After

private

Definition at line 171 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

e5x5Before

float ValidIsoTrkCalib::e5x5Before

private

Definition at line 170 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eAfterDepth1

float ValidIsoTrkCalib::eAfterDepth1

private

Definition at line 152 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eAfterDepth2

float ValidIsoTrkCalib::eAfterDepth2

private

Definition at line 154 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eBeforeDepth1

float ValidIsoTrkCalib::eBeforeDepth1

private

Definition at line 151 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eBeforeDepth2

float ValidIsoTrkCalib::eBeforeDepth2

private

Definition at line 153 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eCentHitAfter

float ValidIsoTrkCalib::eCentHitAfter

private

Definition at line 156 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eCentHitBefore

float ValidIsoTrkCalib::eCentHitBefore

private

Definition at line 155 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eClustAfter

float ValidIsoTrkCalib::eClustAfter

private

Definition at line 143 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eClustBefore

float ValidIsoTrkCalib::eClustBefore

private

Definition at line 142 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eECAL

float ValidIsoTrkCalib::eECAL

private

Definition at line 147 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

emEnergy

Float_t ValidIsoTrkCalib::emEnergy

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

energyECALmip

double ValidIsoTrkCalib::energyECALmip

private

Definition at line 123 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

energyMaxIso

double ValidIsoTrkCalib::energyMaxIso

private

Definition at line 124 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

energyMinIso

double ValidIsoTrkCalib::energyMinIso

private

Definition at line 124 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

etaTrack

float ValidIsoTrkCalib::etaTrack

private

Definition at line 145 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eTrack

float ValidIsoTrkCalib::eTrack

private

Definition at line 144 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

eventNumber

int ValidIsoTrkCalib::eventNumber

private

Definition at line 172 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

fs

edm::Service<TFileService> ValidIsoTrkCalib::fs

private

Definition at line 183 of file ValidIsoTrkCalib.cc.

Referenced by beginJob().

fTree

TTree * ValidIsoTrkCalib::fTree

private

Definition at line 131 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

gen

int ValidIsoTrkCalib::gen

private

Definition at line 109 of file ValidIsoTrkCalib.cc.

genecalLabel_

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

private

Definition at line 84 of file ValidIsoTrkCalib.cc.

genEta

float ValidIsoTrkCalib::genEta[500]

private

Definition at line 110 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

genhbheLabel_

edm::InputTag ValidIsoTrkCalib::genhbheLabel_

private

Definition at line 82 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

genhoLabel_

edm::InputTag ValidIsoTrkCalib::genhoLabel_

private

Definition at line 83 of file ValidIsoTrkCalib.cc.

genPhi

float ValidIsoTrkCalib::genPhi[500]

private

Definition at line 110 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

genPt

float ValidIsoTrkCalib::genPt[500]

private

Definition at line 110 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

geo

const CaloGeometry* ValidIsoTrkCalib::geo

private

Definition at line 80 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

HTime

float ValidIsoTrkCalib::HTime[100]

private

Definition at line 167 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iDr

float ValidIsoTrkCalib::iDr

private

Definition at line 162 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iEta

int ValidIsoTrkCalib::iEta

private

Definition at line 158 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iEtaTr

int ValidIsoTrkCalib::iEtaTr

private

Definition at line 160 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iPhi

int ValidIsoTrkCalib::iPhi

private

Definition at line 159 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iPhiTr

int ValidIsoTrkCalib::iPhiTr

private

Definition at line 161 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iso

int ValidIsoTrkCalib::iso

private

Definition at line 109 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

isoEta

float ValidIsoTrkCalib::isoEta[500]

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

isoPhi

float ValidIsoTrkCalib::isoPhi[500]

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

isoPt

float ValidIsoTrkCalib::isoPt[500]

private

Definition at line 111 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

iTime

float ValidIsoTrkCalib::iTime

private

Definition at line 166 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

maxPNear

double ValidIsoTrkCalib::maxPNear

private

Definition at line 123 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

MinNTECHitsEndcap

int ValidIsoTrkCalib::MinNTECHitsEndcap

private

Definition at line 122 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

MinNTrackHitsBarrel

int ValidIsoTrkCalib::MinNTrackHitsBarrel

private

Definition at line 122 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

Nhits

int ValidIsoTrkCalib::Nhits

private

Definition at line 141 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 181 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numHits

int ValidIsoTrkCalib::numHits

private

Definition at line 148 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numValidTrkHits

float ValidIsoTrkCalib::numValidTrkHits

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numValidTrkStrips

float ValidIsoTrkCalib::numValidTrkStrips

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

numVH

float ValidIsoTrkCalib::numVH

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

numVS

float ValidIsoTrkCalib::numVS

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

outputFileName_

string ValidIsoTrkCalib::outputFileName_

private

Definition at line 103 of file ValidIsoTrkCalib.cc.

parameters_

TrackAssociatorParameters ValidIsoTrkCalib::parameters_

private

Definition at line 76 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 146 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pix

int ValidIsoTrkCalib::pix

private

Definition at line 109 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pixEta

float ValidIsoTrkCalib::pixEta[500]

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pixPhi

float ValidIsoTrkCalib::pixPhi[500]

private

Definition at line 112 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

pixPt

float ValidIsoTrkCalib::pixPt[500]

private

Definition at line 112 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 174 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 177 of file ValidIsoTrkCalib.cc.

Referenced by analyze().

runNumber

int ValidIsoTrkCalib::runNumber

private

Definition at line 173 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 77 of file ValidIsoTrkCalib.cc.

taHCALCone_

double ValidIsoTrkCalib::taHCALCone_

private

Definition at line 78 of file ValidIsoTrkCalib.cc.

takeGenTracks_

bool ValidIsoTrkCalib::takeGenTracks_

private

Definition at line 107 of file ValidIsoTrkCalib.cc.

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

targetE

Float_t ValidIsoTrkCalib::targetE

private

Definition at line 127 of file ValidIsoTrkCalib.cc.

tok_genTrack_

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

private

Definition at line 86 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

tok_geom_

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

private

Definition at line 93 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> ValidIsoTrkCalib::tok_hbhe_

private

Definition at line 87 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

tok_ho_

edm::EDGetTokenT<HORecHitCollection> ValidIsoTrkCalib::tok_ho_

private

Definition at line 88 of file ValidIsoTrkCalib.cc.

Referenced by ValidIsoTrkCalib().

tok_recalibCorrs_

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

private

Definition at line 92 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

tok_track1_

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

private

Definition at line 90 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

tok_track_

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

private

Definition at line 89 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and ValidIsoTrkCalib().

trackAssociator_

TrackDetectorAssociator ValidIsoTrkCalib::trackAssociator_

private

Definition at line 75 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 131 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

xTrkEcal

Float_t ValidIsoTrkCalib::xTrkEcal

private

Definition at line 133 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

xTrkHcal

Float_t ValidIsoTrkCalib::xTrkHcal

private

Definition at line 137 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

yTrkEcal

Float_t ValidIsoTrkCalib::yTrkEcal

private

Definition at line 134 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

yTrkHcal

Float_t ValidIsoTrkCalib::yTrkHcal

private

Definition at line 138 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

zTrkEcal

Float_t ValidIsoTrkCalib::zTrkEcal

private

Definition at line 135 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().

zTrkHcal

Float_t ValidIsoTrkCalib::zTrkHcal

private

Definition at line 139 of file ValidIsoTrkCalib.cc.

Referenced by analyze(), and beginJob().