ValidationHcalIsoTrackAlCaReco Class Reference

#include <Calibration/HcalCalibAlgos/plugins/ValidationHcalIsoTrackAlCaReco.cc>

Inheritance diagram for ValidationHcalIsoTrackAlCaReco:

Public Member Functions
	ValidationHcalIsoTrackAlCaReco (const edm::ParameterSet &)
	~ValidationHcalIsoTrackAlCaReco ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) 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 ()
double	getDist (double, double, double, double)
double	getDistInCM (double eta1, double phi1, double eta2, double phi2)
std::pair< int, int >	towerIndex (double eta, double phi)

Private Attributes
DQMStore *	dbe_
std::string	folderName_
MonitorElement *	hDeposEcalInner
MonitorElement *	hDeposEcalOuter
double	heLow_
double	heUp_
MonitorElement *	hl3AbsEta
MonitorElement *	hl3eta
MonitorElement *	hl3phi
MonitorElement *	hl3Pt
edm::InputTag	hltFilterTag_
MonitorElement *	hOccupancyFull
MonitorElement *	hOccupancyHighEn
MonitorElement *	hOffAbsEta
MonitorElement *	hOffEta
MonitorElement *	hOffEtaFP
MonitorElement *	hOffL3TrackMatch
MonitorElement *	hOffL3TrackPtRat
MonitorElement *	hOffP
MonitorElement *	hOffP_0005
MonitorElement *	hOffP_0510
MonitorElement *	hOffP_1015
MonitorElement *	hOffP_1520
MonitorElement *	hOffPhi
MonitorElement *	hOffPhiFP
MonitorElement *	hPurityEta
MonitorElement *	hPurityPhi
MonitorElement *	hSimAbsEta
MonitorElement *	hSimDist
MonitorElement *	hSimEta
MonitorElement *	hSimN
MonitorElement *	hSimNE
MonitorElement *	hSimNM
MonitorElement *	hSimNN
MonitorElement *	hSimP
MonitorElement *	hSimPhi
MonitorElement *	hSimPt
MonitorElement *	hSimPtRatOff
MonitorElement *	hTracksMaxP
MonitorElement *	hTracksSumP
int	nHLTL3accepts
int	nTotal
std::string	outRootFileName_
double	pThr_
edm::InputTag	recoTrLabel_
bool	saveToFile_
edm::EDGetTokenT < reco::IsolatedPixelTrackCandidateCollection >	tok_arITr_
edm::EDGetTokenT < trigger::TriggerEvent >	tok_hlt_
edm::EDGetTokenT < edm::SimTrackContainer >	tok_simTrack_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
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::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

Description: <one line="" class="" summary>="">

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

Definition at line 73 of file ValidationHcalIsoTrackAlCaReco.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 130 of file ValidationHcalIsoTrackAlCaReco.cc.

References folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), heLow_, hltFilterTag_, nHLTL3accepts, nTotal, outRootFileName_, pThr_, recoTrLabel_, saveToFile_, AlCaHLTBitMon_QueryRunRegistry::string, tok_arITr_, tok_hlt_, and tok_simTrack_.

                                                                                             {

  tok_simTrack_ = consumes<edm::SimTrackContainer>(iConfig.getParameter<edm::InputTag>("simTracksTag"));

  folderName_ = iConfig.getParameter<std::string>("folderName");
  saveToFile_=iConfig.getParameter<bool>("saveToFile");
  outRootFileName_=iConfig.getParameter<std::string>("outputRootFileName");
  tok_hlt_ = consumes<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag>("hltTriggerEventLabel"));
  hltFilterTag_=iConfig.getParameter<edm::InputTag>("hltL3FilterLabel");
  tok_arITr_ = consumes<reco::IsolatedPixelTrackCandidateCollection>(iConfig.getParameter<edm::InputTag>("alcarecoIsoTracksLabel"));
  recoTrLabel_=iConfig.getParameter<edm::InputTag>("recoTracksLabel");
  pThr_=iConfig.getUntrackedParameter<double>("pThrL3",0);
  heLow_=iConfig.getUntrackedParameter<double>("lowerHighEnergyCut",40);
  heLow_=iConfig.getUntrackedParameter<double>("upperHighEnergyCut",60);

  nTotal=0;
  nHLTL3accepts=0;
}

ValidationHcalIsoTrackAlCaReco::~ValidationHcalIsoTrackAlCaReco

(

)

Definition at line 149 of file ValidationHcalIsoTrackAlCaReco.cc.

{}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 152 of file ValidationHcalIsoTrackAlCaReco.cc.

References gather_cfg::cout, HLT_25ns14e33_v1_cff::distance, alignCSCRings::e, trigger::TriggerObject::eta(), MonitorElement::Fill(), edm::Event::getByToken(), getDist(), getDistInCM(), hDeposEcalInner, hDeposEcalOuter, heLow_, heUp_, hl3AbsEta, hl3eta, hl3phi, hl3Pt, hltFilterTag_, hOccupancyFull, hOccupancyHighEn, hOffAbsEta, hOffEta, hOffEtaFP, hOffL3TrackMatch, hOffL3TrackPtRat, hOffP, hOffP_0005, hOffP_0510, hOffP_1015, hOffP_1520, hOffPhi, hOffPhiFP, hSimAbsEta, hSimDist, hSimEta, hSimN, hSimNE, hSimNM, hSimNN, hSimP, hSimPhi, hSimPt, hSimPtRatOff, hTracksMaxP, hTracksSumP, prof2calltree::l, match(), HLT_25ns14e33_v1_cff::nFilt, nTotal, trigger::TriggerObject::p(), trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), pThr_, dttmaxenums::R, mathSSE::sqrt(), tok_arITr_, tok_hlt_, tok_simTrack_, and towerIndex().

{
  nTotal++;

  edm::Handle<trigger::TriggerEvent> trEv;
  iEvent.getByToken(tok_hlt_,trEv);
  
  edm::Handle<reco::IsolatedPixelTrackCandidateCollection> recoIsoTracks;
  iEvent.getByToken(tok_arITr_,recoIsoTracks);

  const trigger::TriggerObjectCollection& TOCol(trEv->getObjects());
  
  trigger::Keys KEYS;
  const trigger::size_type nFilt(trEv->sizeFilters());
  for (trigger::size_type iFilt=0; iFilt!=nFilt; iFilt++) 
    {
      if (trEv->filterTag(iFilt)==hltFilterTag_) 
    {
      KEYS=trEv->filterKeys(iFilt);
    }
    }
  
  trigger::size_type nReg=KEYS.size();
  
  std::vector<double> trigEta;
  std::vector<double> trigPhi;
  bool trig=false;

  //checks with IsoTrack trigger results
  for (trigger::size_type iReg=0; iReg<nReg; iReg++)
    {
      const trigger::TriggerObject& TObj(TOCol[KEYS[iReg]]);
      if (TObj.p()<pThr_) continue;
      hl3eta->Fill(TObj.eta(),1);
      hl3AbsEta->Fill(fabs(TObj.eta()),1);
      hl3phi->Fill(TObj.phi(),1);

      if (recoIsoTracks->size()>0)
    {
      double minRecoL3dist=1000;
      reco::IsolatedPixelTrackCandidateCollection::const_iterator mrtr;
      for (reco::IsolatedPixelTrackCandidateCollection::const_iterator rtrit=recoIsoTracks->begin(); rtrit!=recoIsoTracks->end(); rtrit++)
        {
          double R=getDist(rtrit->eta(),rtrit->phi(),TObj.eta(),TObj.phi()); 
          if (R<minRecoL3dist) 
        {
          mrtr=rtrit;
          minRecoL3dist=R;
        }
        }
      hOffL3TrackMatch->Fill(minRecoL3dist,1);
      hOffL3TrackPtRat->Fill(TObj.pt()/mrtr->pt(),1);
    }
      
      hl3Pt->Fill(TObj.pt(),1);
      trig=true; 
      trigEta.push_back(TObj.eta());
      trigPhi.push_back(TObj.phi());
    }

  //general distributions
  for (reco::IsolatedPixelTrackCandidateCollection::const_iterator itr=recoIsoTracks->begin(); itr!=recoIsoTracks->end(); itr++)
    {
      bool match=false;
      for (unsigned int l=0; l<trigEta.size(); l++)
    {
      if (getDist(itr->eta(),itr->phi(),trigEta[l],trigPhi[l])<0.2) match=true;
    }
      if (match&&trig)
    {   
      hOffEtaFP->Fill(itr->eta(),1);
      hOffPhiFP->Fill(itr->phi(),1);
    }
      
      hOffEta->Fill(itr->eta(),1);
      hOffPhi->Fill(itr->phi(),1);

      hOffAbsEta->Fill(fabs(itr->eta()),1);
   
      hDeposEcalInner->Fill(itr->energyIn(),1);
      hDeposEcalOuter->Fill(itr->energyOut(),1);

      hTracksSumP->Fill(itr->sumPtPxl(),1);
      hTracksMaxP->Fill(itr->maxPtPxl(),1);

      if (fabs(itr->eta())<0.5) hOffP_0005->Fill(itr->p(),1);
      if (fabs(itr->eta())>0.5&&fabs(itr->eta())<1.0) hOffP_0510->Fill(itr->p(),1);
      if (fabs(itr->eta())>1.0&&fabs(itr->eta())<1.5) hOffP_1015->Fill(itr->p(),1);
      if (fabs(itr->eta())<1.5&&fabs(itr->eta())<2.0) hOffP_1520->Fill(itr->p(),1);

      hOffP->Fill(itr->p(),1);

      std::pair<int,int> TI=towerIndex(itr->eta(),itr->phi());
      hOccupancyFull->Fill(TI.first,TI.second,1);
      if (itr->p()>heLow_&&itr->p()<heUp_) hOccupancyHighEn->Fill(TI.first,TI.second,1);
    }
    
// Sergey +

   std::cout <<  std::endl << "  End / Start " << std::endl;

   edm::Handle<edm::SimTrackContainer> simTracks;
   iEvent.getByToken<edm::SimTrackContainer>(tok_simTrack_, simTracks);

  for (reco::IsolatedPixelTrackCandidateCollection::const_iterator bll=recoIsoTracks->begin(); bll!=recoIsoTracks->end(); bll++)
    {         

    std::cout<<"ISO Pt " << bll->pt() <<  " P     "  << bll->p() << " Eta "<< bll->eta() << " Phi "<< bll->phi()<< std::endl;   
    
    double distanceMin = 1.;
    double SimPtMatched = 1.;
    double SimPhiMatched = 1.;
    double SimEtaMatched = 1.;
    double SimDistMatched = 1.;
    double SimPMatched = 1.;
    double neuen = 0.;
    double neuenm = 0.;
    int neun = 0;

   for(edm::SimTrackContainer::const_iterator tracksCI = simTracks->begin(); 
       tracksCI != simTracks->end(); tracksCI++){
        
    int partIndex = tracksCI->genpartIndex();
        if (tracksCI->momentum().eta() > (bll->eta()-0.1) 
       && tracksCI->momentum().eta() < (bll->eta()+0.1)
       && tracksCI->momentum().phi() > (bll->phi()-0.1) 
       && tracksCI->momentum().phi() < (bll->phi()+0.1) 
//     && tracksCI->momentum().e() > (0.5*bll->p())
//     && tracksCI->momentum().e() < (2.*bll->p())
       && tracksCI->momentum().e() > 2.
       && fabs(tracksCI->charge()) == 1 && partIndex >0) 
       {
       
    double distance=getDist(tracksCI->momentum().eta(),tracksCI->momentum().phi(),bll->eta(),bll->phi());
    double distanceCM=getDistInCM(tracksCI->momentum().eta(),tracksCI->momentum().phi(),bll->eta(),bll->phi());
    
    if (distanceMin > distance) {
      distanceMin = distance;
      SimPtMatched = tracksCI->momentum().pt();
      SimPhiMatched = tracksCI->momentum().phi();
      SimEtaMatched = tracksCI->momentum().eta();
      SimDistMatched = distance;
      SimPMatched = sqrt(tracksCI->momentum().pt()*tracksCI->momentum().pt() + tracksCI->momentum().pz()*tracksCI->momentum().pz());
    }

        std::cout<<"    Pt "<<tracksCI->momentum().pt() 
    << " Energy " << tracksCI->momentum().e() 
    << " Eta "<< tracksCI->momentum().eta() 
    << " Phi "<< tracksCI->momentum().phi() 
    << " Ind "  << partIndex
    << " Cha "  << tracksCI->charge() 
    << " Dis "  << distance
    << " DCM "  << distanceCM
    << std::endl;       
    
       }
       
        
        if (
          tracksCI->momentum().eta() > (bll->eta()-0.5) 
       && tracksCI->momentum().eta() < (bll->eta()+0.5)
       && tracksCI->momentum().phi() > (bll->phi()-0.5) 
       && tracksCI->momentum().phi() < (bll->phi()+0.5) 
//     && tracksCI->momentum().e() > 2.
       && tracksCI->charge() == 0 && partIndex >0) 
       {    

    double distance=getDist(tracksCI->momentum().eta(),tracksCI->momentum().phi(),bll->eta(),bll->phi());
    double distanceCM=getDistInCM(tracksCI->momentum().eta(),tracksCI->momentum().phi(),bll->eta(),bll->phi());

        std::cout<<"NEU Pt "<<tracksCI->momentum().pt() 
    << " Energy " << tracksCI->momentum().e() 
    << " Eta "<< tracksCI->momentum().eta() 
    << " Phi "<< tracksCI->momentum().phi() 
    << " Ind "  << partIndex
    << " Cha "  << tracksCI->charge() 
    << " Dis "  << distance
    << " DCM "  << distanceCM
    << std::endl;
    
     if (distanceCM < 40.){
    
      neuen = neuen + tracksCI->momentum().e();
      neun = neun + 1; 
      if (neuenm < tracksCI->momentum().e()) neuenm = tracksCI->momentum().e();
    
     } 
    
    } 
       
       }

          hSimNN->Fill(neun,1);
          hSimNE->Fill(neuen,1);
          hSimNM->Fill(neuenm,1);

    if (distanceMin < 0.1) {

          hSimPt->Fill(SimPtMatched,1);
          hSimPhi->Fill(SimPhiMatched,1);
          hSimEta->Fill(SimEtaMatched,1);
          hSimAbsEta->Fill(fabs(SimEtaMatched),1);
      hSimDist->Fill(SimDistMatched,1);
      hSimPtRatOff->Fill(SimPtMatched/bll->pt(),1);
          hSimP->Fill(SimPMatched,1);
          hSimN->Fill(1,1);

          std::cout<<"S    Pt "<<  SimPtMatched
      << std::endl;
    }
    
    if (distanceMin > 0.1) {
          hSimN->Fill(0,1);
    }
    
    
    
    
  }

// Sergey -      
    
    
    
}

void ValidationHcalIsoTrackAlCaReco::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 378 of file ValidationHcalIsoTrackAlCaReco.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe_, folderName_, MonitorElement::getTH2F(), hDeposEcalInner, hDeposEcalOuter, hl3AbsEta, hl3eta, hl3phi, hl3Pt, hOccupancyFull, hOccupancyHighEn, hOffAbsEta, hOffEta, hOffEtaFP, hOffL3TrackMatch, hOffL3TrackPtRat, hOffP, hOffP_0005, hOffP_0510, hOffP_1015, hOffP_1520, hOffPhi, hOffPhiFP, hSimAbsEta, hSimDist, hSimEta, hSimN, hSimNE, hSimNM, hSimNN, hSimP, hSimPhi, hSimPt, hSimPtRatOff, hTracksMaxP, hTracksSumP, cppFunctionSkipper::operator, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

{
  dbe_ = edm::Service<DQMStore>().operator->();
  dbe_->setCurrentFolder(folderName_);

  hl3Pt=dbe_->book1D("hl3Pt","pT of hlt L3 objects",1000,0,1000);
  hl3Pt->setAxisTitle("pT(GeV)",1);
  hl3eta=dbe_->book1D("hl3eta","eta of hlt L3 objects",16,-2,2);
  hl3eta->setAxisTitle("eta",1);
  hl3AbsEta=dbe_->book1D("hl3AbsEta","|eta| of hlt L3 objects",8,0,2);
  hl3AbsEta->setAxisTitle("eta",1);
  hl3phi=dbe_->book1D("hl3phi","phi of hlt L3 objects",16,-3.2,3.2);
  hl3phi->setAxisTitle("phi",1);

  hOffEta=dbe_->book1D("hOffEta","eta of alcareco objects",100,-2,2);
  hOffEta->setAxisTitle("eta",1);
  hOffPhi=dbe_->book1D("hOffPhi","phi of alcareco objects",100,-3.2,3.2);
  hOffPhi->setAxisTitle("phi",1);
  hOffP=dbe_->book1D("hOffP","p of alcareco objects",1000,0,1000);
  hOffP->setAxisTitle("E(GeV)",1);
  hOffP_0005=dbe_->book1D("hOffP_0005","p of alcareco objects, |eta|<0.5",1000,0,1000);
  hOffP_0005->setAxisTitle("E(GeV)",1);
  hOffP_0510=dbe_->book1D("hOffP_0510","p of alcareco objects, 0.5<|eta|<1.0",1000,0,1000);
  hOffP_0510->setAxisTitle("E(GeV)",1);
  hOffP_1015=dbe_->book1D("hOffP_1015","p of alcareco objects, 1.0<|eta|<1.5",1000,0,1000);
  hOffP_1015->setAxisTitle("E(GeV)",1);
  hOffP_1520=dbe_->book1D("hOffP_1520","p of alcareco objects, 1.5<|eta|<2.0",1000,0,1000);
  hOffP_1520->setAxisTitle("E(GeV)",1);
  hOffEtaFP=dbe_->book1D("hOffEtaFP","eta of alcareco objects, FP",16,-2,2);
  hOffEtaFP->setAxisTitle("eta",1);
  hOffAbsEta=dbe_->book1D("hOffAbsEta","|eta| of alcareco objects",8,0,2);
  hOffAbsEta->setAxisTitle("|eta|",1);
  hOffPhiFP=dbe_->book1D("hOffPhiFP","phi of alcareco objects, FP",16,-3.2,3.2);
  hOffPhiFP->setAxisTitle("phi",1);
  hTracksSumP=dbe_->book1D("hTracksSumP","summary p of tracks in the isolation cone",100,0,20);
  hTracksSumP->setAxisTitle("E(GeV)");
  hTracksMaxP=dbe_->book1D("hTracksMaxP","maximum p among tracks in the isolation cone",100,0,20);
  hTracksMaxP->setAxisTitle("E(GeV)");

  hDeposEcalInner=dbe_->book1D("hDeposEcalInner","ecal energy deposition in inner cone around track",1000,0,1000);
  hDeposEcalInner->setAxisTitle("E(GeV)");
  hDeposEcalOuter=dbe_->book1D("hDeposEcalOuter","ecal energy deposition in outer cone around track",1000,0,1000);
  hDeposEcalOuter->setAxisTitle("E(GeV)");

  hOccupancyFull=dbe_->book2D("hOccupancyFull","number of tracks per tower, full energy range",48,-25,25,73,0,73);
  hOccupancyFull->setAxisTitle("ieta",1);
  hOccupancyFull->setAxisTitle("iphi",2);
  hOccupancyFull->getTH2F()->SetOption("colz");
  hOccupancyFull->getTH2F()->SetStats(kFALSE);
  hOccupancyHighEn=dbe_->book2D("hOccupancyHighEn","number of tracks per tower, high energy tracks",48,-25,25,73,0,73);
  hOccupancyHighEn->setAxisTitle("ieta",1);
  hOccupancyHighEn->setAxisTitle("iphi",2);
  hOccupancyHighEn->getTH2F()->SetOption("colz");
  hOccupancyHighEn->getTH2F()->SetStats(kFALSE);

  hOffL3TrackMatch=dbe_->book1D("hOffL3TrackMatch","Distance from L3 object to offline track",200,0,0.5);
  hOffL3TrackMatch->setAxisTitle("R(eta,phi)",1);
  hOffL3TrackPtRat=dbe_->book1D("hOffL3TrackPtRat","Ratio of pT: L3/offline",100,0,10);
  hOffL3TrackPtRat->setAxisTitle("ratio L3/offline",1);


// Sergey +

  hSimPt=dbe_->book1D("hSimPt","pT of matched SimTrack",1000,0,1000);
  hSimPt->setAxisTitle("pT(GeV)",1);

  hSimPhi=dbe_->book1D("hSimPhi","Phi of matched SimTrack",100,-3.2,3.2);
  hSimPhi->setAxisTitle("phi",1);
  
  hSimEta=dbe_->book1D("hSimEta","Eta of matched SimTrack",100,-2.,2.);
  hSimEta->setAxisTitle("eta",1);
  
  hSimAbsEta=dbe_->book1D("hSimAbsEta","|eta| of matched SimTrack",8,0.,2.);
  hSimAbsEta->setAxisTitle("|eta|",1);

  hSimDist=dbe_->book1D("hSimDist","Distance from matched SimTrack to Offline Track",200,0,0.1);
  hSimDist->setAxisTitle("R(eta,phi)",1);

  hSimPtRatOff=dbe_->book1D("hSimPtRatOff","pT Sim / pT Offline",100,0,10);
  hSimPtRatOff->setAxisTitle("pT Sim / pT Offline",1);

  hSimP=dbe_->book1D("hSimP","p of matched SimTrack",1000,0,1000);
  hSimP->setAxisTitle("p(GeV)",1);
  
  hSimN=dbe_->book1D("hSimN","Number matched",2,0,2);
  hSimN->setAxisTitle("Offline/SimTRack - Matched or Not",1);

  hSimNN=dbe_->book1D("hSimNN","Number of the neutral particles in cone on ECAL",100,0,100);
  hSimNN->setAxisTitle("Number",1);

  hSimNE=dbe_->book1D("hSimNE","Total energy of the neutral particles in cone on ECAL",100,0,100);
  hSimNE->setAxisTitle("Energy",1);

  hSimNM=dbe_->book1D("hSimNM","Maximum energy of the neutral particles in cone on ECAL",100,0,100);
  hSimNM->setAxisTitle("Energy",1);
  
// Sergey -


}

void ValidationHcalIsoTrackAlCaReco::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 479 of file ValidationHcalIsoTrackAlCaReco.cc.

References dbe_, outRootFileName_, DQMStore::save(), and saveToFile_.

                                            {

if(dbe_) 
  {
    if (saveToFile_) dbe_->save(outRootFileName_);
  }
}

double ValidationHcalIsoTrackAlCaReco::getDist ( double  eta1, double  phi1, double  eta2, double  phi2  )

private

Definition at line 73 of file ValidationHcalIsoTrackAlCaReco.cc.

References funct::pow(), and mathSSE::sqrt().

Referenced by analyze().

{
  double dphi = fabs(phi1 - phi2); 
  if(dphi>acos(-1)) dphi = 2*acos(-1)-dphi;
  double dr = sqrt(dphi*dphi + pow(eta1-eta2,2));
  return dr;
}

double ValidationHcalIsoTrackAlCaReco::getDistInCM ( double  eta1, double  phi1, double  eta2, double  phi2  )

private

Definition at line 83 of file ValidationHcalIsoTrackAlCaReco.cc.

References funct::cos(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, create_public_lumi_plots::exp, funct::sin(), and funct::tan().

Referenced by analyze().

{
  double dR, Rec;
  double theta1=2*atan(exp(-eta1));
  double theta2=2*atan(exp(-eta2));
  if (fabs(eta1)<1.479) Rec=129;
  else Rec=275;
  //|vect| times tg of acos(scalar product)
  dR=fabs((Rec/sin(theta1))*tan(acos(sin(theta1)*sin(theta2)*(sin(phi1)*sin(phi2)+cos(phi1)*cos(phi2))+cos(theta1)*cos(theta2))));
  return dR;
}

std::pair< int, int > ValidationHcalIsoTrackAlCaReco::towerIndex ( double  eta, double  phi  )

private

Definition at line 97 of file ValidationHcalIsoTrackAlCaReco.cc.

References delta, and i.

Referenced by analyze().

{
  int ieta=0;
  int iphi=0;
  for (int i=1; i<21; i++)
    {
      if (fabs(eta)<(i*0.087)&&fabs(eta)>(i-1)*0.087) ieta=int(fabs(eta)/eta)*i;
    }
  if (fabs(eta)>1.740&&fabs(eta)<1.830) ieta=int(fabs(eta)/eta)*21;
  if (fabs(eta)>1.830&&fabs(eta)<1.930) ieta=int(fabs(eta)/eta)*22;
  if (fabs(eta)>1.930&&fabs(eta)<2.043) ieta=int(fabs(eta)/eta)*23;

  double delta=phi+0.174532925;
  if (delta<0) delta=delta+2*acos(-1);
  if (fabs(eta)<1.740) 
    {
      for (int i=0; i<72; i++)
    {
      if (delta<(i+1)*0.087266462&&delta>i*0.087266462) iphi=i;
    }
    }
  else 
    {
      for (int i=0; i<36; i++)
    {
      if (delta<2*(i+1)*0.087266462&&delta>2*i*0.087266462) iphi=2*i;
    }
    }

  return std::pair<int,int>(ieta,iphi);
}

Member Data Documentation

DQMStore* ValidationHcalIsoTrackAlCaReco::dbe_

private

Definition at line 81 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by beginJob(), and endJob().

std::string ValidationHcalIsoTrackAlCaReco::folderName_

private

Definition at line 87 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by beginJob(), and ValidationHcalIsoTrackAlCaReco().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hDeposEcalInner

private

Definition at line 118 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hDeposEcalOuter

private

Definition at line 119 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

double ValidationHcalIsoTrackAlCaReco::heLow_

private

Definition at line 98 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().

double ValidationHcalIsoTrackAlCaReco::heUp_

private

Definition at line 99 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hl3AbsEta

private

Definition at line 103 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hl3eta

private

Definition at line 102 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hl3phi

private

Definition at line 104 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hl3Pt

private

Definition at line 101 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

edm::InputTag ValidationHcalIsoTrackAlCaReco::hltFilterTag_

private

Definition at line 90 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOccupancyFull

private

Definition at line 128 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOccupancyHighEn

private

Definition at line 129 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffAbsEta

private

Definition at line 122 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffEta

private

Definition at line 125 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffEtaFP

private

Definition at line 121 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffL3TrackMatch

private

Definition at line 105 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffL3TrackPtRat

private

Definition at line 106 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffP

private

Definition at line 113 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffP_0005

private

Definition at line 108 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffP_0510

private

Definition at line 109 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffP_1015

private

Definition at line 110 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffP_1520

private

Definition at line 111 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffPhi

private

Definition at line 126 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hOffPhiFP

private

Definition at line 123 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hPurityEta

private

Definition at line 131 of file ValidationHcalIsoTrackAlCaReco.h.

MonitorElement* ValidationHcalIsoTrackAlCaReco::hPurityPhi

private

Definition at line 132 of file ValidationHcalIsoTrackAlCaReco.h.

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimAbsEta

private

Definition at line 139 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimDist

private

Definition at line 140 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimEta

private

Definition at line 138 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimN

private

Definition at line 143 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimNE

private

Definition at line 145 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimNM

private

Definition at line 146 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimNN

private

Definition at line 144 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimP

private

Definition at line 142 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimPhi

private

Definition at line 137 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimPt

private

Definition at line 136 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hSimPtRatOff

private

Definition at line 141 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hTracksMaxP

private

Definition at line 116 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

MonitorElement* ValidationHcalIsoTrackAlCaReco::hTracksSumP

private

Definition at line 115 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and beginJob().

int ValidationHcalIsoTrackAlCaReco::nHLTL3accepts

private

Definition at line 151 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by ValidationHcalIsoTrackAlCaReco().

int ValidationHcalIsoTrackAlCaReco::nTotal

private

Definition at line 150 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().

std::string ValidationHcalIsoTrackAlCaReco::outRootFileName_

private

Definition at line 89 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by endJob(), and ValidationHcalIsoTrackAlCaReco().

double ValidationHcalIsoTrackAlCaReco::pThr_

private

Definition at line 97 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().

edm::InputTag ValidationHcalIsoTrackAlCaReco::recoTrLabel_

private

Definition at line 91 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by ValidationHcalIsoTrackAlCaReco().

bool ValidationHcalIsoTrackAlCaReco::saveToFile_

private

Definition at line 88 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by endJob(), and ValidationHcalIsoTrackAlCaReco().

edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> ValidationHcalIsoTrackAlCaReco::tok_arITr_

private

Definition at line 94 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().

edm::EDGetTokenT<trigger::TriggerEvent> ValidationHcalIsoTrackAlCaReco::tok_hlt_

private

Definition at line 93 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().

edm::EDGetTokenT<edm::SimTrackContainer> ValidationHcalIsoTrackAlCaReco::tok_simTrack_

private

Definition at line 95 of file ValidationHcalIsoTrackAlCaReco.h.

Referenced by analyze(), and ValidationHcalIsoTrackAlCaReco().