HcalIsoTrkAnalyzer Class Reference

#include <Calibration/HcalCalibAlgos/src/HcalIsoTrkAnalyzer.cc>

Inheritance diagram for HcalIsoTrkAnalyzer:

Public Member Functions
double	getDistInPlaneSimple (const GlobalPoint caloPoint, const GlobalPoint rechitPoint)
	HcalIsoTrkAnalyzer (const edm::ParameterSet &)
	~HcalIsoTrkAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

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_
UInt_t	cell
Float_t	cellEnergy
TClonesArray *	cells
TRefArray *	cells3x3
TRefArray *	cellsPF
char	dirname [50]
InputTag	eLabel_
Float_t	emEnergy
double	energyECALmip
double	energyMaxIso
double	energyMinIso
Float_t	etVetoJet
UInt_t	eventNumber
float	eventWeight
const CaloGeometry *	geo
InputTag	hbheLabel_
char	hname [20]
InputTag	hoLabel_
char	htitle [80]
Int_t	iEtaHit
ofstream	input_to_L3
UInt_t	iPhiHit
std::string	m_ebInstance
std::string	m_ecalLabel
std::string	m_eeInstance
std::string	m_hcalLabel
int	m_histoFlag
std::string	m_inputTrackLabel
double	maxPNear
int	MinNTECHitsEndcap
int	MinNTrackHitsBarrel
int	nIterations
UInt_t	numberOfCells
string	outputFileName_
TrackAssociatorParameters	parameters_
Float_t	probeJetEmFrac
TLorentzVector *	probeJetP4
TFile *	rootFile
UInt_t	runNumber
double	rvert
Float_t	tagJetEmFrac
TLorentzVector *	tagJetP4
Float_t	targetE
TrackDetectorAssociator	trackAssociator_
double	trackEta
InputTag	trackLabel1_
InputTag	trackLabel_
double	trackPhi
TTree *	tree
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)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Description: see twiki for details: https://twiki.cern.ch/twiki/bin/view/CMS/IsoTrackAnalyzer

Implementation: <Notes on="" implementation>="">

Definition at line 98 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 239 of file HcalIsoTrkAnalyzer.cc.

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

{
  
  m_ecalLabel = iConfig.getUntrackedParameter<std::string> ("ecalRecHitsLabel","ecalRecHit");
  m_ebInstance = iConfig.getUntrackedParameter<std::string> ("ebRecHitsInstance","EcalRecHitsEB");
  m_eeInstance = iConfig.getUntrackedParameter<std::string> ("eeRecHitsInstance","EcalRecHitsEE");
  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");

  AxB_=iConfig.getParameter<std::string>("AxB");
  calibrationConeSize_=iConfig.getParameter<double>("calibrationConeSize");

  nIterations = iConfig.getParameter<int>("noOfIterations");
  eventWeight = iConfig.getParameter<double>("eventWeight");
  energyMinIso = iConfig.getParameter<double>("energyMinIso");
  energyMaxIso = iConfig.getParameter<double>("energyMaxIso");
  energyECALmip = iConfig.getParameter<double>("energyECALmip");
  maxPNear = iConfig.getParameter<double>("maxPNear");
  MinNTrackHitsBarrel =  iConfig.getParameter<int>("MinNTrackHitsBarrel");
  MinNTECHitsEndcap =  iConfig.getParameter<int>("MinNTECHitsEndcap");

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

}

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer

(

)

Definition at line 274 of file HcalIsoTrkAnalyzer.cc.

{
}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 280 of file HcalIsoTrkAnalyzer.cc.

References abs, edm::SortedCollection< T, SORT >::begin(), MaxHit_struct::depthhit, patCandidatesForDimuonsSequences_cff::ecal, EcalBarrel, EcalEndcap, edm::SortedCollection< T, SORT >::end(), PV3DBase< T, PVType, FrameType >::eta(), edm::EventID::event(), GeneralTracks_cfi::generalTracks, edm::EventSetup::get(), edm::Event::getByLabel(), CaloSubdetectorGeometry::getClosestCell(), getDistInPlaneSimple(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), DetId::Hcal, HcalBarrel, HcalEndcap, MaxHit_struct::hitenergy, i, edm::EventBase::id(), HcalDetId::ieta(), MaxHit_struct::ietahitm, info, HcalDetId::iphi(), MaxHit_struct::iphihitm, PV3DBase< T, PVType, FrameType >::phi(), pos, MaxHit_struct::posMax, edm::Handle< T >::product(), edm::SortedCollection< T, SORT >::push_back(), edm::EventID::run(), inputsource_file_cfi::runNumber, mathSSE::sqrt(), reco::btau::trackEta, reco::btau::trackPhi, and diffTreeTool::tree.

{
  using namespace edm;
  using namespace std;

  vector<float> rawEnergyVec;
  vector<int> detiphi;
  vector<int> detieta;
  vector<int> i3i5;
  vector<HcalDetId> detidvec;
  float calEnergy;

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

  edm::Handle<reco::IsolatedPixelTrackCandidateCollection> trackCollection;
  iEvent.getByLabel(trackLabel_,trackCollection);
    
  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();

// rof 16.05.2008 start: include the possibility for recalibration (use "recalibrate" label for safety)
/*
  edm::ESHandle <HcalRespCorrs> recalibCorrs;
  iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
  const HcalRespCorrs* myRecalib = recalibCorrs.product();
*/
// rof end

  parameters_.useEcal = true;
  parameters_.useHcal = true;
  parameters_.useCalo = false;
  parameters_.useMuon = false;
  parameters_.dREcal = 0.5;
  parameters_.dRHcal = 0.6;  

  if (trackCollection->size()==0) return;

  for (reco::TrackCollection::const_iterator trit=generalTracks->begin(); trit!=generalTracks->end(); trit++)
    {
      reco::IsolatedPixelTrackCandidateCollection::const_iterator isoMatched=trackCollection->begin();
      bool matched=false;
      for (reco::IsolatedPixelTrackCandidateCollection::const_iterator it = trackCollection->begin(); it!=trackCollection->end(); it++)
    { 

      if (abs((trit->pt() - it->pt())/it->pt()) < 0.005 && abs(trit->eta() - it->eta()) < 0.01) 
        {
          isoMatched=it;
          matched=true;
          break;
        }
    }
      if (!matched) continue;
      
      if (trit->hitPattern().numberOfValidHits()<MinNTrackHitsBarrel) continue;
      if (fabs(trit->eta())>1.47&&trit->hitPattern().numberOfValidStripTECHits()<MinNTECHitsEndcap) continue;
      
      //container for used recHits
      std::vector<int> usedHits; 
      //      

      calEnergy = sqrt(trit->px()*trit->px()+trit->py()*trit->py()+trit->pz()*trit->pz()+0.14*0.14);

      trackEta = trit->eta();
      trackPhi = trit->phi();
      
      double corrHCAL = 1.; //another possibility for correction.  - why?


      //      cout << endl << " ISO TRACK E = "<< calEnergy << " ETA = " << trackEta<< " PHI = " << trackPhi <<  " Correction " <<  corrHCAL<< endl;
      
      rvert = sqrt(trit->vx()*trit->vx()+trit->vy()*trit->vy()+trit->vz()*trit->vz());      
      
      //Associate track with a calorimeter
      TrackDetMatchInfo info = trackAssociator_.associate(iEvent, iSetup,trackAssociator_.getFreeTrajectoryState(iSetup, *trit),parameters_);
     //*(it->track().get())
      double etaecal=info.trkGlobPosAtEcal.eta();
      double phiecal=info.trkGlobPosAtEcal.phi();

      double etahcal=info.trkGlobPosAtHcal.eta();
      //double 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();
    
    GlobalPoint gP(xTrkHcal,yTrkHcal,zTrkHcal);

    int iphitrue = -10;
    int ietatrue = 100;

    if (etahcal<1.392) 
      {
        const CaloSubdetectorGeometry* gHB = geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
        //    const GlobalPoint tempPoint(newx, newy, newz);
        //const DetId tempId = gHB->getClosestCell(tempPoint);
        const HcalDetId tempId = gHB->getClosestCell(gP);
        ietatrue = tempId.ieta();
        iphitrue = tempId.iphi();
      }

    if (etahcal>1.392 &&  etahcal<3.0) 
      {
        const CaloSubdetectorGeometry* gHE = geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap);
        const HcalDetId tempId = gHE->getClosestCell(gP);
        ietatrue = tempId.ieta();
        iphitrue = tempId.iphi();
      }

    //      eecal=info.coneEnergy(parameters_.dREcal, TrackDetMatchInfo::EcalRecHits);
//      ehcal=info.coneEnergy(parameters_.dRHcal, TrackDetMatchInfo::HcalRecHits);

      double rmin = 0.07;
      if( fabs(etaecal) > 1.47 ) rmin = 0.07*(fabs(etaecal)-0.47)*1.2;
      if( fabs(etaecal) > 2.2 ) rmin = 0.07*(fabs(etaecal)-0.47)*1.4;

      MaxHit_struct MaxHit;

      // Find Ecal RecHit with maximum energy and collect other information
      MaxHit.hitenergy=-100;
      
      double econus = 0.;
      float ecal_cluster = 0.;      

      //clear usedHits
      usedHits.clear();
      //

      for (std::vector<EcalRecHit>::const_iterator ehit=Hitecal.begin(); ehit!=Hitecal.end(); ehit++)   
    {
      //check that this hit was not considered before and push it into usedHits
      bool hitIsUsed=false;
      int hitHashedIndex=-10000; 
      if (ehit->id().subdetId()==EcalBarrel) 
        {
          EBDetId did(ehit->id());
          hitHashedIndex=did.hashedIndex();
        }
      
      if (ehit->id().subdetId()==EcalEndcap) 
        {
          EEDetId did(ehit->id());
          hitHashedIndex=did.hashedIndex();
        }
      for (uint32_t i=0; i<usedHits.size(); i++)
        {
          if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
        }
      if (hitIsUsed) continue;
      usedHits.push_back(hitHashedIndex);
      //

      if((*ehit).energy() > MaxHit.hitenergy) 
        {
          MaxHit.hitenergy = (*ehit).energy();
        }

      GlobalPoint pos = geo->getPosition((*ehit).detid());
      double phihit = pos.phi();
      double etahit = pos.eta();
     
      double dphi = fabs(phiecal - phihit); 
      if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
      double deta = fabs(etaecal - etahit); 
      double dr = sqrt(dphi*dphi + deta*deta);

      if (dr < rmin) {
        econus = econus + (*ehit).energy();
      }
    
        if (dr < 0.13) ecal_cluster += (*ehit).energy();

    }
      MaxHit.hitenergy=-100;
      
      //clear usedHits
      usedHits.clear();
      //
      
      //float dddeta = 1000.;
      //float dddphi = 1000.;
      //int iphitrue = 1234;
      //int ietatrue = 1234;

      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;
      int hitHashedIndex=hhit->id().hashed_index();
      for (uint32_t i=0; i<usedHits.size(); i++)
        {
          if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
        }
      if (hitIsUsed) continue;
      usedHits.push_back(hitHashedIndex);
      //

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

      GlobalPoint pos = geo->getPosition(hhit->detid());
      double phihit = pos.phi();
      double etahit = pos.eta();
      
      int iphihitm  = (hhit->id()).iphi();
      int ietahitm  = (hhit->id()).ieta();
      int depthhit = (hhit->id()).depth();
      double enehit = hhit->energy() * recal;
      
      if (depthhit!=1) continue;
       
      double dphi = fabs(info.trkGlobPosAtHcal.phi() - phihit); 
      if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
      double deta = fabs(info.trkGlobPosAtHcal.eta() - etahit); 
      double dr = sqrt(dphi*dphi + deta*deta);
      
      /* Commented. Another way of finding true projection is implemented
      if (deta<dddeta) {
       ietatrue = ietahitm;
       dddeta=deta;
      }
      
      if (dphi<dddphi) {
      iphitrue = iphihitm;
      dddphi=dphi;
      }*/
      
      if(dr<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();
          double enehit2 = hhit2->energy() * recal;
          
          if ( iphihitm==iphihitm2 && ietahitm==ietahitm2  && depthhit!=depthhit2){
            
            enehit = enehit+enehit2;
            
          }
          
        }
          
          if(enehit > MaxHit.hitenergy) 
        {
          MaxHit.hitenergy =  enehit;
          MaxHit.ietahitm   = (hhit->id()).ieta();
          MaxHit.iphihitm   = (hhit->id()).iphi();
          MaxHit.depthhit  = (hhit->id()).depth();
          
          MaxHit.posMax = geo->getPosition(hhit->detid());
          
        }
        }
    }
      

      Bool_t passCuts = kFALSE;
      if(calEnergy > energyMinIso && calEnergy < energyMaxIso && isoMatched->energyIn() < energyECALmip && 
     isoMatched->maxPtPxl() < maxPNear && abs(MaxHit.ietahitm)<30 && MaxHit.hitenergy > 0.){ passCuts = kTRUE; }
      

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

      //clear usedHits
      usedHits.clear();
      //

      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;
          int hitHashedIndex=hhit->id().hashed_index();
          for (uint32_t i=0; i<usedHits.size(); i++)
        {
          if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
        }
          if (hitIsUsed) continue;
          usedHits.push_back(hitHashedIndex);
          //

          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;
          /*AP DIPHI = DIPHI<-36 ? 72+DIPHI : DIPHI; */
          
          int numbercell=0;
          if(AxB_=="3x3") numbercell = 1;
          if(AxB_=="5x5") numbercell = 2;
          if(AxB_=="7x7") numbercell = 3;
          if(AxB_=="Cone") numbercell = 1000;
          
          if( abs(DIETA)<=numbercell && (abs(DIPHI)<=numbercell || ( abs(MaxHit.ietahitm)>=20 && abs(DIPHI)<=numbercell+1)) )  {
        
        // rof 16.05.2008 start: include the possibility for recalibration
        float recal = 1;
        
            int iii3i5 = 0;

            
        const GlobalPoint pos2 = geo->getPosition(hhit->detid());       
        
        if(passCuts){

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

          rawEnergyVec.push_back(hhit->energy() * recal * corrHCAL);
          detidvec.push_back(hhit->id());
          detiphi.push_back((hhit->id()).iphi());
          detieta.push_back((hhit->id()).ieta());
          i3i5.push_back(iii3i5);
          
          }
         
          if (AxB_=="Cone" && getDistInPlaneSimple(gP,pos2) < calibrationConeSize_) {
      
          rawEnergyVec.push_back(hhit->energy() * recal * corrHCAL);
          detidvec.push_back(hhit->id());
          detiphi.push_back((hhit->id()).iphi());
          detieta.push_back((hhit->id()).ieta());
          i3i5.push_back(iii3i5);

                  }
        
        }
          }
        }
    }

      if(AxB_!="3x3" && AxB_!="5x5" && AxB_!="7x7" && AxB_!="Cone") LogWarning(" AxB ")<<"   Not supported: "<< AxB_;
      
      if(passCuts){
    
    input_to_L3 << rawEnergyVec.size() << "   " << calEnergy;

    
    for (unsigned int i=0; i<rawEnergyVec.size(); i++)
      {
        input_to_L3 << "   " << rawEnergyVec.at(i) << "   " << detidvec.at(i).rawId() ;    

      }
    input_to_L3 <<endl;
    
        eventNumber = iEvent.id().event();
        runNumber = iEvent.id().run();
        iEtaHit = ietatrue;
        iPhiHit = iphitrue;
        emEnergy = isoMatched->energyIn();
//        exampleP4->SetPxPyPzE(2, 1, 1, 10);
    
    numberOfCells=rawEnergyVec.size();
    targetE = calEnergy;
        
        for (unsigned int ia=0; ia<numberOfCells; ++ia) {
      cellEnergy = rawEnergyVec.at(ia);
      cell = detidvec.at(ia).rawId();
          
      new((*cells)[ia])  TCell(cell, cellEnergy);
      

        } 
    
        tree->Fill();
        
        cells->Clear();   
    
      }
      
      rawEnergyVec.clear();
      detidvec.clear();
      detiphi.clear();
      detieta.clear();
      i3i5.clear();

      try {
    Handle<HORecHitCollection> ho;
    iEvent.getByLabel(hoLabel_,ho);
    const HORecHitCollection Hitho = *(ho.product());
    
    //clear usedHits
    usedHits.clear();
    //
    
    for(HORecHitCollection::const_iterator hoItr=Hitho.begin(); hoItr!=Hitho.end(); hoItr++)
      {

        //check that this hit was not considered before and push it into usedHits
        bool hitIsUsed=false;
        int hitHashedIndex=hoItr->id().hashed_index();
        for (uint32_t i=0; i<usedHits.size(); i++)
          {
        if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
          }
        if (hitIsUsed) continue;
        usedHits.push_back(hitHashedIndex);

        /*AP
        GlobalPoint pos = geo->getPosition(hoItr->detid());
        double phihit = pos.phi();
        double etahit = pos.eta();
        
        int iphihitm = (hoItr->id()).iphi();
        int ietahitm = (hoItr->id()).ieta();
        int depthhit = (hoItr->id()).depth();
        
        double dphi = fabs(trackPhi - phihit); 
        if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
        double deta = fabs(trackEta - etahit); 
        double dr = sqrt(dphi*dphi + deta*deta);
        */
        
      }
      } catch (cms::Exception& e) { // can't find it!
    if (!allowMissingInputs_) throw e;
      }
    }
}

void HcalIsoTrkAnalyzer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 736 of file HcalIsoTrkAnalyzer.cc.

References BuildWebpage::rootFile, inputsource_file_cfi::runNumber, and diffTreeTool::tree.

{

  //  MyL3Algo = new MinL3AlgoUniv<HcalDetId>(eventWeight);

    input_to_L3.open("input_to_l3.txt");

    rootFile = new TFile("rootFile.root", "RECREATE");
    tree = new TTree("hcalCalibTree", "Tree for IsoTrack Calibration");  
    cells = new TClonesArray("TCell", 10000);
//    cells3x3 = new TRefArray;
//    cellsPF = new TRefArray;
//    exampleP4 = new TLorentzVector();
    tagJetP4 = new TLorentzVector(); // dijet
    probeJetP4 = new TLorentzVector();  // dijet   
    
    tree->Branch("cells", &cells, 64000, 0);
    tree->Branch("targetE", &targetE, "targetE/F");
    tree->Branch("emEnergy", &emEnergy, "emEnergy/F");

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

//    tree->Branch("exampleP4", "TLorentzVector", &exampleP4);
//    tree->Branch("cells3x3", &cells3x3, 64000, 0);
//    tree->Branch("cellsPF", &cellsPF, 64000, 0);
    tree->Branch("iEtaHit", &iEtaHit, "iEtaHit/I");
    tree->Branch("iPhiHit", &iPhiHit, "iPhiHit/i");
    tree->Branch("eventNumber", &eventNumber, "eventNumber/i");
    tree->Branch("runNumber", &runNumber, "runNumber/i");

    tree->Branch("tagJetP4", "TLorentzVector", &tagJetP4); // dijet
    tree->Branch("probeJetP4", "TLorentzVector", &probeJetP4); // dijet
    tree->Branch("etVetoJet", &etVetoJet, "etVetoJet/F"); // dijet
    tree->Branch("tagJetEmFrac", &tagJetEmFrac,"tagJetEmFrac/F"); // dijet
    tree->Branch("probeJetEmFrac", &probeJetEmFrac,"probeJetEmFrac/F"); // dijet

}

void HcalIsoTrkAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 781 of file HcalIsoTrkAnalyzer.cc.

References BuildWebpage::rootFile.

                           {

    input_to_L3.close();

    rootFile->Write();
    rootFile->Close();

    if (cells) delete cells;
//    if (cells3x3) delete cells3x3;
//    if (cellsPF) delete cellsPF;
//    if (exampleP4) delete exampleP4;
    if (tagJetP4) delete tagJetP4; // dijet
    if (probeJetP4) delete probeJetP4; // dijet

}

double HcalIsoTrkAnalyzer::getDistInPlaneSimple	(	const GlobalPoint	caloPoint,
		const GlobalPoint	rechitPoint
	)

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

References Vector3DBase< T, FrameTag >::dot(), PV3DBase< T, PVType, FrameType >::mag(), Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
  
  // Simplified version of getDistInPlane
  // Assume track direction is origin -> point of hcal intersection
  
  const GlobalVector caloIntersectVector(caloPoint.x(), 
                     caloPoint.y(), 
                     caloPoint.z());

  const GlobalVector caloIntersectUnitVector = caloIntersectVector.unit();
  
  const GlobalVector rechitVector(rechitPoint.x(),
                  rechitPoint.y(),
                  rechitPoint.z());

  const GlobalVector rechitUnitVector = rechitVector.unit();

  double dotprod = caloIntersectUnitVector.dot(rechitUnitVector);
  double rechitdist = caloIntersectVector.mag()/dotprod;
  
  
  const GlobalVector effectiveRechitVector = rechitdist*rechitUnitVector;
  const GlobalPoint effectiveRechitPoint(effectiveRechitVector.x(),
                     effectiveRechitVector.y(),
                     effectiveRechitVector.z());
  
  
  GlobalVector distance_vector = effectiveRechitPoint-caloPoint;
  
  if (dotprod > 0.)
  {
    return distance_vector.mag();
  }
  else
  {
    return 999999.;
    
  }

    
}

Member Data Documentation

bool HcalIsoTrkAnalyzer::allowMissingInputs_

private

Definition at line 132 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::associationConeSize_

private

Definition at line 130 of file HcalIsoTrkAnalyzer.cc.

string HcalIsoTrkAnalyzer::AxB_

private

Definition at line 131 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::calibrationConeSize_

private

Definition at line 130 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::cell

private

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::cellEnergy

private

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

TClonesArray* HcalIsoTrkAnalyzer::cells

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

TRefArray* HcalIsoTrkAnalyzer::cells3x3

private

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

TRefArray* HcalIsoTrkAnalyzer::cellsPF

private

Definition at line 168 of file HcalIsoTrkAnalyzer.cc.

char HcalIsoTrkAnalyzer::dirname[50]

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

InputTag HcalIsoTrkAnalyzer::eLabel_

private

Definition at line 119 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::emEnergy

private

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::energyECALmip

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::energyMaxIso

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::energyMinIso

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::etVetoJet

private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::eventNumber

private

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

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

float HcalIsoTrkAnalyzer::eventWeight

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

const CaloGeometry* HcalIsoTrkAnalyzer::geo

private

Definition at line 116 of file HcalIsoTrkAnalyzer.cc.

InputTag HcalIsoTrkAnalyzer::hbheLabel_

private

Definition at line 117 of file HcalIsoTrkAnalyzer.cc.

char HcalIsoTrkAnalyzer::hname[20]

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

InputTag HcalIsoTrkAnalyzer::hoLabel_

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

char HcalIsoTrkAnalyzer::htitle[80]

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Int_t HcalIsoTrkAnalyzer::iEtaHit

private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

ofstream HcalIsoTrkAnalyzer::input_to_L3

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::iPhiHit

private

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_ebInstance

private

Definition at line 125 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_ecalLabel

private

Definition at line 124 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_eeInstance

private

Definition at line 126 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_hcalLabel

private

Definition at line 127 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::m_histoFlag

private

Definition at line 128 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_inputTrackLabel

private

Definition at line 123 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::maxPNear

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::MinNTECHitsEndcap

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::MinNTrackHitsBarrel

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::nIterations

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::numberOfCells

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

string HcalIsoTrkAnalyzer::outputFileName_

private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

TrackAssociatorParameters HcalIsoTrkAnalyzer::parameters_

private

Definition at line 114 of file HcalIsoTrkAnalyzer.cc.

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

Float_t HcalIsoTrkAnalyzer::probeJetEmFrac

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

TLorentzVector* HcalIsoTrkAnalyzer::probeJetP4

private

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

TFile* HcalIsoTrkAnalyzer::rootFile

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::runNumber

private

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::rvert

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::tagJetEmFrac

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

TLorentzVector* HcalIsoTrkAnalyzer::tagJetP4

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::targetE

private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

TrackDetectorAssociator HcalIsoTrkAnalyzer::trackAssociator_

private

Definition at line 113 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::trackEta

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

InputTag HcalIsoTrkAnalyzer::trackLabel1_

private

Definition at line 121 of file HcalIsoTrkAnalyzer.cc.

InputTag HcalIsoTrkAnalyzer::trackLabel_

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::trackPhi

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by svgfig.SVG::__str__(), and python.cmstools.EventTree::SetAlias().

Float_t HcalIsoTrkAnalyzer::xTrkEcal

private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::xTrkHcal

private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::yTrkEcal

private

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::yTrkHcal

private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::zTrkEcal

private

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::zTrkHcal

private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.