ValidIsoTrkCalib Class Reference

Inheritance diagram for ValidIsoTrkCalib:

Public Member Functions
	ValidIsoTrkCalib (const edm::ParameterSet &)
	~ValidIsoTrkCalib ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

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]
InputTag	genTracksLabel_
const CaloGeometry *	geo
InputTag	hbheLabel_
InputTag	hoLabel_
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
TrackDetectorAssociator	trackAssociator_
float	trackE
float	trackEta
InputTag	trackLabel1_
InputTag	trackLabel_
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
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 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)

Detailed Description

Definition at line 67 of file ValidIsoTrkCalib.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 198 of file ValidIsoTrkCalib.cc.

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

{

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

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

  hbheLabel_= iConfig.getParameter<edm::InputTag>("hbheInput");
  hoLabel_=iConfig.getParameter<edm::InputTag>("hoInput");
  //eLabel_=iConfig.getParameter<edm::InputTag>("eInput");
  trackLabel_ = iConfig.getParameter<edm::InputTag>("HcalIsolTrackInput");
  trackLabel1_ = 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");
  parameters_.loadParameters( parameters );
  trackAssociator_.useDefaultPropagator();

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

}

ValidIsoTrkCalib::~ValidIsoTrkCalib

(

)

Definition at line 244 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

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 255 of file ValidIsoTrkCalib.cc.

References abs, gather_cfg::cout, MaxHit_struct::depthhit, MaxHit_struct::dr, alignCSCRings::e, edm::EventID::event(), HcalObjRepresent::Fill(), relval_steps::gen(), MergeTrackCollections_cff::generalTracks, edm::EventSetup::get(), edm::Event::getByLabel(), getDistInPlaneSimple(), reco::JetExtendedAssociation::getValue(), DetId::Hcal, HcalBarrel, MaxHit_struct::hitenergy, i, edm::EventBase::id(), HcalDetId::ieta(), MaxHit_struct::ietahitm, info, HcalDetId::iphi(), MaxHit_struct::iphihitm, n, pos, edm::ESHandle< class >::product(), edm::EventID::run(), convertSQLiteXML::runNumber, mathSSE::sqrt(), reco::btau::trackEta, reco::btau::trackPhi, listHistos::trackPt, TrackDetMatchInfo::trkGlobPosAtEcal, and TrackDetMatchInfo::trkGlobPosAtHcal.

{
   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.getByLabel(genTracksLabel_, generalTracks);

  edm::Handle<reco::TrackCollection> isoProdTracks;
  iEvent.getByLabel(trackLabel1_,isoProdTracks);


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

  edm::Handle<HBHERecHitCollection> hbhe;
  iEvent.getByLabel(hbheLabel_,hbhe);
  const HBHERecHitCollection Hithbhe = *(hbhe.product());

  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  geo = pG.product();
  
  const CaloSubdetectorGeometry* gHcal = 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->size()==0) 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();

}

void ValidIsoTrkCalib::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 774 of file ValidIsoTrkCalib.cc.

References relval_steps::gen(), and convertSQLiteXML::runNumber.

{

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


}

void ValidIsoTrkCalib::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 875 of file ValidIsoTrkCalib.cc.

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

Member Data Documentation

bool ValidIsoTrkCalib::allowMissingInputs_

private

Definition at line 108 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::associationConeSize_

private

Definition at line 104 of file ValidIsoTrkCalib.cc.

string ValidIsoTrkCalib::AxB_

private

Definition at line 105 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::calibrationConeSize_

private

Definition at line 106 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::CentHitFactor

private

Definition at line 165 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::delR

private

Definition at line 170 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::dietatr

private

Definition at line 171 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::diphitr

private

Definition at line 172 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::e3x3After

private

Definition at line 177 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::e3x3Before

private

Definition at line 176 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::e5x5After

private

Definition at line 179 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::e5x5Before

private

Definition at line 178 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eAfterDepth1

private

Definition at line 160 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eAfterDepth2

private

Definition at line 162 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eBeforeDepth1

private

Definition at line 159 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eBeforeDepth2

private

Definition at line 161 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eCentHitAfter

private

Definition at line 164 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eCentHitBefore

private

Definition at line 163 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eClustAfter

private

Definition at line 151 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eClustBefore

private

Definition at line 150 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eECAL

private

Definition at line 155 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::emEnergy

private

Definition at line 134 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::energyECALmip

private

Definition at line 130 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::energyMaxIso

private

Definition at line 131 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::energyMinIso

private

Definition at line 131 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::etaTrack

private

Definition at line 153 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::eTrack

private

Definition at line 152 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::eventNumber

private

Definition at line 180 of file ValidIsoTrkCalib.cc.

Referenced by Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::setFilterBranches().

edm::Service<TFileService> ValidIsoTrkCalib::fs

private

Definition at line 191 of file ValidIsoTrkCalib.cc.

TTree * ValidIsoTrkCalib::fTree

private

Definition at line 139 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::gen

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

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

private

Definition at line 95 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::genEta[500]

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::genhbheLabel_

private

Definition at line 93 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::genhoLabel_

private

Definition at line 94 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::genPhi[500]

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::genPt[500]

private

Definition at line 116 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::genTracksLabel_

private

Definition at line 92 of file ValidIsoTrkCalib.cc.

const CaloGeometry* ValidIsoTrkCalib::geo

private

Definition at line 91 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::hbheLabel_

private

Definition at line 96 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::hoLabel_

private

Definition at line 97 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::HTime[100]

private

Definition at line 175 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::iDr

private

Definition at line 170 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::iEta

private

Definition at line 166 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::iEtaTr

private

Definition at line 168 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::iPhi

private

Definition at line 167 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::iPhiTr

private

Definition at line 169 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::iso

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::isoEta[500]

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::isoPhi[500]

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::isoPt[500]

private

Definition at line 117 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::iTime

private

Definition at line 174 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::maxPNear

private

Definition at line 130 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::MinNTECHitsEndcap

private

Definition at line 129 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::MinNTrackHitsBarrel

private

Definition at line 129 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::Nhits

private

Definition at line 149 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::NisoTrk

private

Definition at line 123 of file ValidIsoTrkCalib.cc.

TH1F* ValidIsoTrkCalib::nTracks

private

Definition at line 189 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::numHits

private

Definition at line 156 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::numValidTrkHits

private

Definition at line 183 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::numValidTrkStrips

private

Definition at line 183 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::numVH

private

Definition at line 183 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::numVS

private

Definition at line 183 of file ValidIsoTrkCalib.cc.

string ValidIsoTrkCalib::outputFileName_

private

Definition at line 109 of file ValidIsoTrkCalib.cc.

TrackAssociatorParameters ValidIsoTrkCalib::parameters_

private

Definition at line 87 of file ValidIsoTrkCalib.cc.

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

float ValidIsoTrkCalib::phiTrack

private

Definition at line 154 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::pix

private

Definition at line 115 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::pixEta[500]

private

Definition at line 118 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::pixPhi[500]

private

Definition at line 118 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::pixPt[500]

private

Definition at line 118 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::ptNear

private

Definition at line 125 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::PtNearBy

private

Definition at line 182 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::ptrack

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

const HcalRespCorrs* ValidIsoTrkCalib::respRecalib

private

Definition at line 185 of file ValidIsoTrkCalib.cc.

int ValidIsoTrkCalib::runNumber

private

Definition at line 181 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::rvert

private

Definition at line 126 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::taECALCone_

private

Definition at line 88 of file ValidIsoTrkCalib.cc.

double ValidIsoTrkCalib::taHCALCone_

private

Definition at line 89 of file ValidIsoTrkCalib.cc.

bool ValidIsoTrkCalib::takeGenTracks_

private

Definition at line 113 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::targetE

private

Definition at line 135 of file ValidIsoTrkCalib.cc.

TrackDetectorAssociator ValidIsoTrkCalib::trackAssociator_

private

Definition at line 86 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::trackE

private

Definition at line 124 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::trackEta

private

Definition at line 124 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::trackLabel1_

private

Definition at line 99 of file ValidIsoTrkCalib.cc.

InputTag ValidIsoTrkCalib::trackLabel_

private

Definition at line 98 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::trackPhi

private

Definition at line 124 of file ValidIsoTrkCalib.cc.

float ValidIsoTrkCalib::trackPt

private

Definition at line 124 of file ValidIsoTrkCalib.cc.

TTree* ValidIsoTrkCalib::tTree

private

Definition at line 139 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::xTrkEcal

private

Definition at line 141 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::xTrkHcal

private

Definition at line 145 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::yTrkEcal

private

Definition at line 142 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::yTrkHcal

private

Definition at line 146 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::zTrkEcal

private

Definition at line 143 of file ValidIsoTrkCalib.cc.

Float_t ValidIsoTrkCalib::zTrkHcal

private

Definition at line 147 of file ValidIsoTrkCalib.cc.