CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

EgammaObjects Class Reference

#include <EgammaObjects.h>

Inheritance diagram for EgammaObjects:
edm::EDAnalyzer

List of all members.

Public Member Functions

virtual void analyze (const edm::Event &, const edm::EventSetup &)
virtual void beginJob ()
 EgammaObjects (const edm::ParameterSet &)
virtual void endJob ()
 ~EgammaObjects ()

Private Member Functions

void analyzeElectrons (const edm::Event &, const edm::EventSetup &)
void analyzePhotons (const edm::Event &, const edm::EventSetup &)
void applyLabels ()
void createBookedHistoObjects ()
void createTempHistoObjects ()
float ecalEta (float EtaParticle, float Zvertex, float plane_Radius)
double findRecoMass (reco::Photon pOne, reco::Photon pTwo)
double findRecoMass (reco::GsfElectron eOne, reco::GsfElectron eTwo)
void fitHistos ()
void getDeltaResHistosViaSlicing ()
void getEfficiencyHistosViaDividing ()
void loadCMSSWObjects (const edm::ParameterSet &ps)
void loadHistoParameters (const edm::ParameterSet &ps)
void saveHistos ()
void setDrawOptions ()

Private Attributes

TH2D * _TEMP_scatterPlot_deltaEtaVsE_
TH2D * _TEMP_scatterPlot_deltaEtaVsEt_
TH2D * _TEMP_scatterPlot_deltaEtaVsEta_
TH2D * _TEMP_scatterPlot_deltaEtaVsPhi_
TH2D * _TEMP_scatterPlot_deltaPhiVsE_
TH2D * _TEMP_scatterPlot_deltaPhiVsEt_
TH2D * _TEMP_scatterPlot_deltaPhiVsEta_
TH2D * _TEMP_scatterPlot_deltaPhiVsPhi_
TH2D * _TEMP_scatterPlot_EOverTruthVsE_
TH2D * _TEMP_scatterPlot_EOverTruthVsEt_
TH2D * _TEMP_scatterPlot_EOverTruthVsEta_
TH2D * _TEMP_scatterPlot_EOverTruthVsPhi_
TH2D * _TEMP_scatterPlot_EtOverTruthVsE_
TH2D * _TEMP_scatterPlot_EtOverTruthVsEt_
TH2D * _TEMP_scatterPlot_EtOverTruthVsEta_
TH2D * _TEMP_scatterPlot_EtOverTruthVsPhi_
int EtCut
TH1D * hist_All_recoMass_
TH1D * hist_BarrelOnly_recoMass_
int hist_bins_deltaEta_
int hist_bins_deltaPhi_
int hist_bins_E_
int hist_bins_EOverTruth_
int hist_bins_Et_
int hist_bins_Eta_
int hist_bins_EtaOverTruth_
int hist_bins_EtOverTruth_
int hist_bins_Phi_
int hist_bins_PhiOverTruth_
int hist_bins_recoMass_
TH1D * hist_deltaEtaVsE_
TH1D * hist_deltaEtaVsEt_
TH1D * hist_deltaEtaVsEta_
TH1D * hist_deltaEtaVsPhi_
TH1D * hist_deltaPhiVsE_
TH1D * hist_deltaPhiVsEt_
TH1D * hist_deltaPhiVsEta_
TH1D * hist_deltaPhiVsPhi_
TH1D * hist_E_
TH1D * hist_EEfficiency_
TH1D * hist_EndcapOnly_recoMass_
TH1D * hist_ENumRecoOverNumTrue_
TH1D * hist_EOverTruth_
TH1D * hist_EOverTruthVsE_
TH1D * hist_EOverTruthVsEt_
TH1D * hist_EOverTruthVsEta_
TH1D * hist_EOverTruthVsPhi_
TH1D * hist_Et_
TH1D * hist_Eta_
TH1D * hist_EtaEfficiency_
TH1D * hist_EtaNumRecoOverNumTrue_
TH1D * hist_EtaOverTruth_
TH1D * hist_EtEfficiency_
TH1D * hist_EtNumRecoOverNumTrue_
TH1D * hist_EtOverTruth_
TH1D * hist_EtOverTruthVsE_
TH1D * hist_EtOverTruthVsEt_
TH1D * hist_EtOverTruthVsEta_
TH1D * hist_EtOverTruthVsPhi_
double hist_max_deltaEta_
double hist_max_deltaPhi_
double hist_max_E_
double hist_max_EOverTruth_
double hist_max_Et_
double hist_max_Eta_
double hist_max_EtaOverTruth_
double hist_max_EtOverTruth_
double hist_max_Phi_
double hist_max_PhiOverTruth_
double hist_max_recoMass_
double hist_min_deltaEta_
double hist_min_deltaPhi_
double hist_min_E_
double hist_min_EOverTruth_
double hist_min_Et_
double hist_min_Eta_
double hist_min_EtaOverTruth_
double hist_min_EtOverTruth_
double hist_min_Phi_
double hist_min_PhiOverTruth_
double hist_min_recoMass_
TH1D * hist_Mixed_recoMass_
TH1D * hist_Phi_
TH1D * hist_PhiEfficiency_
TH1D * hist_PhiNumRecoOverNumTrue_
TH1D * hist_PhiOverTruth_
TH1D * hist_recoMass_withBackgroud_10EtCut_
TH1D * hist_recoMass_withBackgroud_20EtCut_
TH1D * hist_recoMass_withBackgroud_5EtCut_
TH1D * hist_recoMass_withBackgroud_NoEtCut_
TH1D * hist_resolutionEtaVsE_
TH1D * hist_resolutionEtaVsEt_
TH1D * hist_resolutionEtaVsEta_
TH1D * hist_resolutionEtaVsPhi_
TH1D * hist_resolutionEtVsE_
TH1D * hist_resolutionEtVsEt_
TH1D * hist_resolutionEtVsEta_
TH1D * hist_resolutionEtVsPhi_
TH1D * hist_resolutionEVsE_
TH1D * hist_resolutionEVsEt_
TH1D * hist_resolutionEVsEta_
TH1D * hist_resolutionEVsPhi_
TH1D * hist_resolutionPhiVsE_
TH1D * hist_resolutionPhiVsEt_
TH1D * hist_resolutionPhiVsEta_
TH1D * hist_resolutionPhiVsPhi_
edm::InputTag MCTruthCollection_
int particleID
std::string particleString
edm::InputTag RecoCollection_
TFile * rootFile_

Detailed Description

Definition at line 21 of file EgammaObjects.h.


Constructor & Destructor Documentation

EgammaObjects::EgammaObjects ( const edm::ParameterSet ps) [explicit]

Definition at line 17 of file EgammaObjects.cc.

References _TEMP_scatterPlot_deltaEtaVsE_, _TEMP_scatterPlot_deltaEtaVsEt_, _TEMP_scatterPlot_deltaEtaVsEta_, _TEMP_scatterPlot_deltaEtaVsPhi_, _TEMP_scatterPlot_deltaPhiVsE_, _TEMP_scatterPlot_deltaPhiVsEt_, _TEMP_scatterPlot_deltaPhiVsEta_, _TEMP_scatterPlot_deltaPhiVsPhi_, _TEMP_scatterPlot_EOverTruthVsE_, _TEMP_scatterPlot_EOverTruthVsEt_, _TEMP_scatterPlot_EOverTruthVsEta_, _TEMP_scatterPlot_EOverTruthVsPhi_, _TEMP_scatterPlot_EtOverTruthVsE_, _TEMP_scatterPlot_EtOverTruthVsEt_, _TEMP_scatterPlot_EtOverTruthVsEta_, _TEMP_scatterPlot_EtOverTruthVsPhi_, EtCut, edm::ParameterSet::getParameter(), hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_deltaEtaVsE_, hist_deltaEtaVsEt_, hist_deltaEtaVsEta_, hist_deltaEtaVsPhi_, hist_deltaPhiVsE_, hist_deltaPhiVsEt_, hist_deltaPhiVsEta_, hist_deltaPhiVsPhi_, hist_EndcapOnly_recoMass_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, hist_resolutionEtaVsE_, hist_resolutionEtaVsEt_, hist_resolutionEtaVsEta_, hist_resolutionEtaVsPhi_, hist_resolutionPhiVsE_, hist_resolutionPhiVsEt_, hist_resolutionPhiVsEta_, hist_resolutionPhiVsPhi_, loadCMSSWObjects(), loadHistoParameters(), particleID, particleString, and rootFile_.

{
  particleID = ps.getParameter<int>("particleID");
  EtCut = ps.getParameter<int>("EtCut");

  if( particleID == 22 ) particleString = "Photon";
  else if( particleID == 11 ) particleString = "Electron";
  else
    throw(std::runtime_error("\n\nEgammaObjects: Only 11 = Photon and 22 = Electron are acceptable parictleIDs! Exiting...\n\n"));

  loadCMSSWObjects(ps);
  loadHistoParameters(ps);

  rootFile_ = TFile::Open(ps.getParameter<std::string>("outputFile").c_str(),"RECREATE");

  hist_EtaEfficiency_ = 0 ;
  hist_EtaNumRecoOverNumTrue_ = 0 ;
  hist_deltaEtaVsEt_ = 0 ;
  hist_deltaEtaVsE_ = 0 ;
  hist_deltaEtaVsEta_ = 0 ;
  hist_deltaEtaVsPhi_ = 0 ;
  hist_resolutionEtaVsEt_ = 0 ;
  hist_resolutionEtaVsE_ = 0 ;
  hist_resolutionEtaVsEta_ = 0 ;
  hist_resolutionEtaVsPhi_ = 0 ;

  hist_Phi_ = 0 ;
  hist_PhiOverTruth_ = 0 ;
  hist_PhiEfficiency_ = 0 ;
  hist_PhiNumRecoOverNumTrue_ = 0 ;
  hist_deltaPhiVsEt_ = 0 ;
  hist_deltaPhiVsE_ = 0 ;
  hist_deltaPhiVsEta_ = 0 ;
  hist_deltaPhiVsPhi_ = 0 ;
  hist_resolutionPhiVsEt_ = 0 ;
  hist_resolutionPhiVsE_ = 0 ;
  hist_resolutionPhiVsEta_ = 0 ;
  hist_resolutionPhiVsPhi_ = 0 ;

  hist_All_recoMass_ = 0 ;
  hist_BarrelOnly_recoMass_ = 0 ;
  hist_EndcapOnly_recoMass_ = 0 ;
  hist_Mixed_recoMass_ = 0 ;

  hist_recoMass_withBackgroud_NoEtCut_ = 0 ;
  hist_recoMass_withBackgroud_5EtCut_ = 0 ;
  hist_recoMass_withBackgroud_10EtCut_ = 0 ;
  hist_recoMass_withBackgroud_20EtCut_ = 0 ;

  _TEMP_scatterPlot_EtOverTruthVsEt_ = 0 ;
  _TEMP_scatterPlot_EtOverTruthVsE_ = 0 ;
  _TEMP_scatterPlot_EtOverTruthVsEta_ = 0 ;
  _TEMP_scatterPlot_EtOverTruthVsPhi_ = 0 ;

  _TEMP_scatterPlot_EOverTruthVsEt_ = 0 ;
  _TEMP_scatterPlot_EOverTruthVsE_ = 0 ;
  _TEMP_scatterPlot_EOverTruthVsEta_ = 0 ;
  _TEMP_scatterPlot_EOverTruthVsPhi_ = 0 ;

  _TEMP_scatterPlot_deltaEtaVsEt_ = 0 ;
  _TEMP_scatterPlot_deltaEtaVsE_ = 0 ;
  _TEMP_scatterPlot_deltaEtaVsEta_ = 0 ;
  _TEMP_scatterPlot_deltaEtaVsPhi_ = 0 ;

  _TEMP_scatterPlot_deltaPhiVsEt_ = 0 ;
  _TEMP_scatterPlot_deltaPhiVsE_ = 0 ;
  _TEMP_scatterPlot_deltaPhiVsEta_ = 0 ;
  _TEMP_scatterPlot_deltaPhiVsPhi_ = 0 ;

}
EgammaObjects::~EgammaObjects ( )

Definition at line 141 of file EgammaObjects.cc.

References rootFile_.

{
  delete rootFile_;
}

Member Function Documentation

void EgammaObjects::analyze ( const edm::Event evt,
const edm::EventSetup es 
) [virtual]

Implements edm::EDAnalyzer.

Definition at line 295 of file EgammaObjects.cc.

References analyzeElectrons(), analyzePhotons(), and particleID.

{
  if( particleID == 22 ) analyzePhotons(evt, es);
  else if( particleID == 11 ) analyzeElectrons(evt, es);
}
void EgammaObjects::analyzeElectrons ( const edm::Event evt,
const edm::EventSetup es 
) [private]

Definition at line 463 of file EgammaObjects.cc.

References _TEMP_scatterPlot_deltaEtaVsE_, _TEMP_scatterPlot_deltaEtaVsEt_, _TEMP_scatterPlot_deltaEtaVsEta_, _TEMP_scatterPlot_deltaEtaVsPhi_, _TEMP_scatterPlot_deltaPhiVsE_, _TEMP_scatterPlot_deltaPhiVsEt_, _TEMP_scatterPlot_deltaPhiVsEta_, _TEMP_scatterPlot_deltaPhiVsPhi_, _TEMP_scatterPlot_EOverTruthVsE_, _TEMP_scatterPlot_EOverTruthVsEt_, _TEMP_scatterPlot_EOverTruthVsEta_, _TEMP_scatterPlot_EOverTruthVsPhi_, _TEMP_scatterPlot_EtOverTruthVsE_, _TEMP_scatterPlot_EtOverTruthVsEt_, _TEMP_scatterPlot_EtOverTruthVsEta_, _TEMP_scatterPlot_EtOverTruthVsPhi_, abs, SiPixelRawToDigiRegional_cfi::deltaPhi, deltaR(), HI_PhotonSkim_cff::electrons, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), EtCut, findRecoMass(), MCTruth::genEvent, edm::Event::getByLabel(), hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, edm::HandleBase::isValid(), edm::InputTag::label(), MCTruthCollection_, Geom::pi(), pi, funct::pow(), edm::Handle< T >::product(), RecoCollection_, mathSSE::sqrt(), and reco::GsfElectron::superCluster().

Referenced by analyze().

{
  edm::Handle<reco::GsfElectronCollection> pElectrons;
  evt.getByLabel(RecoCollection_, pElectrons);
  if (!pElectrons.isValid()) {
    edm::LogError("DOEPlotsProducerElectrons") << "Error! can't get collection with label " << RecoCollection_.label();
  }

  const reco::GsfElectronCollection* electrons = pElectrons.product();
  std::vector<reco::GsfElectron> electronsMCMatched;

  for(reco::GsfElectronCollection::const_iterator aClus = electrons->begin(); aClus != electrons->end(); aClus++)
  {
    if(aClus->et() >= EtCut)
    {
      hist_Et_->Fill(aClus->et());
      hist_E_->Fill(aClus->energy());
      hist_Eta_->Fill(aClus->eta());
      hist_Phi_->Fill(aClus->phi());
     }
  }

  for(int firstElectron = 0, numElectrons = electrons->size(); firstElectron < numElectrons - 1; firstElectron++)
    for(int secondElectron = firstElectron + 1; secondElectron < numElectrons; secondElectron++)
    {
      reco::GsfElectron eOne = (*electrons)[firstElectron];
      reco::GsfElectron eTwo = (*electrons)[secondElectron];

      double recoMass = findRecoMass(eOne, eTwo);

      hist_recoMass_withBackgroud_NoEtCut_->Fill(recoMass);

      if(eOne.et() >= 5 && eTwo.et() >= 5)
        hist_recoMass_withBackgroud_5EtCut_->Fill(recoMass);

      if(eOne.et() >= 10 && eTwo.et() >= 10)
        hist_recoMass_withBackgroud_10EtCut_->Fill(recoMass);

      if(eOne.et() >= 20 && eTwo.et() >= 20)
        hist_recoMass_withBackgroud_20EtCut_->Fill(recoMass);
    }

  edm::Handle<edm::HepMCProduct> pMCTruth ;
  evt.getByLabel(MCTruthCollection_, pMCTruth);
  if (!pMCTruth.isValid()) {
    edm::LogError("DOEPlotsProducerElectrons") << "Error! can't get collection with label " << MCTruthCollection_.label();
  }

  const HepMC::GenEvent* genEvent = pMCTruth->GetEvent();
  for(HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
      currentParticle != genEvent->particles_end(); currentParticle++ )
  {
    if(abs((*currentParticle)->pdg_id())==11 && (*currentParticle)->status()==1
      && (*currentParticle)->momentum().e()/cosh((*currentParticle)->momentum().eta()) >= EtCut)
    {
                        double phiTrue = (*currentParticle)->momentum().phi();
                        double etaTrue = (*currentParticle)->momentum().eta();
      double eTrue  = (*currentParticle)->momentum().e();
                        double etTrue  = (*currentParticle)->momentum().e()/cosh(etaTrue);

      double etaCurrent, etaFound = -999;
      double phiCurrent, phiFound = -999;
      double etCurrent,  etFound  = -999;
      double eCurrent,   eFound  = -999;

      reco::GsfElectron bestMatchElectron;

      double closestParticleDistance = 999;

      for(reco::GsfElectronCollection::const_iterator aClus = electrons->begin(); aClus != electrons->end(); aClus++)
      {
        if(aClus->et() > EtCut)
        {
          etaCurrent = aClus->eta();
          phiCurrent = aClus->phi();
          etCurrent  = aClus->et();
          eCurrent  = aClus->energy();

          double deltaPhi = phiCurrent-phiTrue;
          if(deltaPhi > Geom::pi()) deltaPhi -= 2.*Geom::pi();
          if(deltaPhi < -Geom::pi()) deltaPhi += 2.*Geom::pi();
          double deltaR = std::sqrt(std::pow(etaCurrent-etaTrue,2)+std::pow(deltaPhi,2));

              if(deltaR < closestParticleDistance)
          {
            etFound  = etCurrent;
            eFound   = eCurrent;
            etaFound = etaCurrent;
            phiFound = phiCurrent;
            closestParticleDistance = deltaR;
            bestMatchElectron = *aClus;
          }
        }
      }

      if(closestParticleDistance < 0.05 && etFound/etTrue > .5 && etFound/etTrue < 1.5)
      {
        hist_EtOverTruth_->Fill(etFound/etTrue);
        hist_EOverTruth_->Fill(eFound/eTrue);
        hist_EtaOverTruth_->Fill(etaFound/etaTrue);
        hist_PhiOverTruth_->Fill(phiFound/phiTrue);

        hist_EtEfficiency_->Fill(etTrue);
        hist_EEfficiency_->Fill(eTrue);
        hist_EtaEfficiency_->Fill(etaTrue);
        hist_PhiEfficiency_->Fill(phiTrue);

        double deltaPhi = phiFound-phiTrue;
        if(deltaPhi > Geom::pi()) deltaPhi -= 2.*Geom::pi();
        if(deltaPhi < -Geom::pi()) deltaPhi += 2.*Geom::pi();

        _TEMP_scatterPlot_EtOverTruthVsEt_->Fill(etFound,etFound/etTrue);
        _TEMP_scatterPlot_EtOverTruthVsE_->Fill(eFound,etFound/etTrue);
        _TEMP_scatterPlot_EtOverTruthVsEta_->Fill(etaFound,etFound/etTrue);
        _TEMP_scatterPlot_EtOverTruthVsPhi_->Fill(phiFound,etFound/etTrue);

        _TEMP_scatterPlot_EOverTruthVsEt_->Fill(etFound,eFound/eTrue);
        _TEMP_scatterPlot_EOverTruthVsE_->Fill(eFound,eFound/eTrue);
        _TEMP_scatterPlot_EOverTruthVsEta_->Fill(etaFound,eFound/eTrue);
        _TEMP_scatterPlot_EOverTruthVsPhi_->Fill(phiFound,eFound/eTrue);

        _TEMP_scatterPlot_deltaEtaVsEt_->Fill(etFound,etaFound-etaTrue);
        _TEMP_scatterPlot_deltaEtaVsE_->Fill(eFound,etaFound-etaTrue);
        _TEMP_scatterPlot_deltaEtaVsEta_->Fill(etaFound,etaFound-etaTrue);
        _TEMP_scatterPlot_deltaEtaVsPhi_->Fill(phiFound,etaFound-etaTrue);

        _TEMP_scatterPlot_deltaPhiVsEt_->Fill(etFound,deltaPhi);
        _TEMP_scatterPlot_deltaPhiVsE_->Fill(eFound,deltaPhi);
        _TEMP_scatterPlot_deltaPhiVsEta_->Fill(etaFound,deltaPhi);
        _TEMP_scatterPlot_deltaPhiVsPhi_->Fill(phiFound,deltaPhi);

        electronsMCMatched.push_back(bestMatchElectron);
      }

      hist_EtNumRecoOverNumTrue_->Fill(etTrue);
      hist_ENumRecoOverNumTrue_->Fill(eTrue);
      hist_EtaNumRecoOverNumTrue_->Fill(etaTrue);
      hist_PhiNumRecoOverNumTrue_->Fill(phiTrue);
    }
  }

  if(electronsMCMatched.size() == 2)
  {
    reco::GsfElectron eOne = electronsMCMatched[0];
    reco::GsfElectron eTwo = electronsMCMatched[1];

    double recoMass = findRecoMass(eOne, eTwo);

    hist_All_recoMass_->Fill(recoMass);

    if(eOne.superCluster()->seed()->algo() == 1 && eTwo.superCluster()->seed()->algo() == 1)
      hist_BarrelOnly_recoMass_->Fill(recoMass);
    else if(eOne.superCluster()->seed()->algo() == 0 && eTwo.superCluster()->seed()->algo() == 0)
      hist_EndcapOnly_recoMass_->Fill(recoMass);
    else
      hist_Mixed_recoMass_->Fill(recoMass);
  }
}
void EgammaObjects::analyzePhotons ( const edm::Event evt,
const edm::EventSetup es 
) [private]

Definition at line 301 of file EgammaObjects.cc.

References _TEMP_scatterPlot_deltaEtaVsE_, _TEMP_scatterPlot_deltaEtaVsEt_, _TEMP_scatterPlot_deltaEtaVsEta_, _TEMP_scatterPlot_deltaEtaVsPhi_, _TEMP_scatterPlot_deltaPhiVsE_, _TEMP_scatterPlot_deltaPhiVsEt_, _TEMP_scatterPlot_deltaPhiVsEta_, _TEMP_scatterPlot_deltaPhiVsPhi_, _TEMP_scatterPlot_EOverTruthVsE_, _TEMP_scatterPlot_EOverTruthVsEt_, _TEMP_scatterPlot_EOverTruthVsEta_, _TEMP_scatterPlot_EOverTruthVsPhi_, _TEMP_scatterPlot_EtOverTruthVsE_, _TEMP_scatterPlot_EtOverTruthVsEt_, _TEMP_scatterPlot_EtOverTruthVsEta_, _TEMP_scatterPlot_EtOverTruthVsPhi_, abs, SiPixelRawToDigiRegional_cfi::deltaPhi, deltaR(), ecalEta(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), EtCut, findRecoMass(), MCTruth::genEvent, edm::Event::getByLabel(), hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, edm::HandleBase::isValid(), edm::InputTag::label(), MCTruthCollection_, interactiveExample::photons, Geom::pi(), pi, funct::pow(), edm::Handle< T >::product(), RecoCollection_, mathSSE::sqrt(), and reco::Photon::superCluster().

Referenced by analyze().

{
  edm::Handle<reco::PhotonCollection> pPhotons;
  evt.getByLabel(RecoCollection_, pPhotons);
  if (!pPhotons.isValid()) {
    edm::LogError("EgammaObjects") << "Error! can't get collection with label " << RecoCollection_.label();
  }

  const reco::PhotonCollection* photons = pPhotons.product();
  std::vector<reco::Photon> photonsMCMatched;

  for(reco::PhotonCollection::const_iterator aClus = photons->begin(); aClus != photons->end(); aClus++)
  {
    if(aClus->et() >= EtCut)
    {
      hist_Et_->Fill(aClus->et());
      hist_E_->Fill(aClus->energy());
      hist_Eta_->Fill(aClus->eta());
      hist_Phi_->Fill(aClus->phi());
     }
  }

  for(int firstPhoton = 0, numPhotons = photons->size(); firstPhoton < numPhotons - 1; firstPhoton++)
    for(int secondPhoton = firstPhoton + 1; secondPhoton < numPhotons; secondPhoton++)
    {
      reco::Photon pOne = (*photons)[firstPhoton];
      reco::Photon pTwo = (*photons)[secondPhoton];

      double recoMass = findRecoMass(pOne, pTwo);

      hist_recoMass_withBackgroud_NoEtCut_->Fill(recoMass);

      if(pOne.et() >= 5 && pTwo.et() >= 5)
        hist_recoMass_withBackgroud_5EtCut_->Fill(recoMass);

      if(pOne.et() >= 10 && pTwo.et() >= 10)
        hist_recoMass_withBackgroud_10EtCut_->Fill(recoMass);

      if(pOne.et() >= 20 && pTwo.et() >= 20)
        hist_recoMass_withBackgroud_20EtCut_->Fill(recoMass);
    }

  edm::Handle<edm::HepMCProduct> pMCTruth ;
  evt.getByLabel(MCTruthCollection_, pMCTruth);
  if (!pMCTruth.isValid()) {
    edm::LogError("EgammaObjects") << "Error! can't get collection with label " << MCTruthCollection_.label();
  }

  const HepMC::GenEvent* genEvent = pMCTruth->GetEvent();

  for(HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
      currentParticle != genEvent->particles_end(); currentParticle++ )
  {
    if(abs((*currentParticle)->pdg_id())==22 && (*currentParticle)->status()==1
      && (*currentParticle)->momentum().e()/cosh(ecalEta((*currentParticle)->momentum().eta(), (*currentParticle)->production_vertex()->position().z()/10.,
      (*currentParticle)->production_vertex()->position().perp()/10.)) >= EtCut)
    {
                        HepMC::FourVector vtx = (*currentParticle)->production_vertex()->position();
                        double phiTrue = (*currentParticle)->momentum().phi();
                        double etaTrue = ecalEta((*currentParticle)->momentum().eta(), vtx.z()/10., vtx.perp()/10.);
      double eTrue  = (*currentParticle)->momentum().e();
                        double etTrue  = (*currentParticle)->momentum().e()/cosh(etaTrue);

      double etaCurrent, etaFound = -999;
      double phiCurrent, phiFound = -999;
      double etCurrent,  etFound  = -999;
      double eCurrent,   eFound  = -999;

      reco::Photon bestMatchPhoton;

      double closestParticleDistance = 999;

      for(reco::PhotonCollection::const_iterator aClus = photons->begin(); aClus != photons->end(); aClus++)
      {
        if(aClus->et() > EtCut)
        {
          etaCurrent = aClus->eta();
          phiCurrent = aClus->phi();
          etCurrent  = aClus->et();
          eCurrent  = aClus->energy();

          double deltaPhi = phiCurrent-phiTrue;
          if(deltaPhi > Geom::pi()) deltaPhi -= 2.*Geom::pi();
          if(deltaPhi < -Geom::pi()) deltaPhi += 2.*Geom::pi();
          double deltaR = std::sqrt(std::pow(etaCurrent-etaTrue,2)+std::pow(deltaPhi,2));

              if(deltaR < closestParticleDistance)
          {
            etFound  = etCurrent;
            eFound   = eCurrent;
            etaFound = etaCurrent;
            phiFound = phiCurrent;
            closestParticleDistance = deltaR;
            bestMatchPhoton = *aClus;
          }
        }
      }

      if(closestParticleDistance < 0.05 && etFound/etTrue > .5 && etFound/etTrue < 1.5)
      {
        hist_EtOverTruth_->Fill(etFound/etTrue);
        hist_EOverTruth_->Fill(eFound/eTrue);
        hist_EtaOverTruth_->Fill(etaFound/etaTrue);
        hist_PhiOverTruth_->Fill(phiFound/phiTrue);

        hist_EtEfficiency_->Fill(etTrue);
        hist_EEfficiency_->Fill(eTrue);
        hist_EtaEfficiency_->Fill(etaTrue);
        hist_PhiEfficiency_->Fill(phiTrue);

        double deltaPhi = phiFound-phiTrue;
        if(deltaPhi > Geom::pi()) deltaPhi -= 2.*Geom::pi();
        if(deltaPhi < -Geom::pi()) deltaPhi += 2.*Geom::pi();

        _TEMP_scatterPlot_EtOverTruthVsEt_->Fill(etFound,etFound/etTrue);
        _TEMP_scatterPlot_EtOverTruthVsE_->Fill(eFound,etFound/etTrue);
        _TEMP_scatterPlot_EtOverTruthVsEta_->Fill(etaFound,etFound/etTrue);
        _TEMP_scatterPlot_EtOverTruthVsPhi_->Fill(phiFound,etFound/etTrue);

        _TEMP_scatterPlot_EOverTruthVsEt_->Fill(etFound,eFound/eTrue);
        _TEMP_scatterPlot_EOverTruthVsE_->Fill(eFound,eFound/eTrue);
        _TEMP_scatterPlot_EOverTruthVsEta_->Fill(etaFound,eFound/eTrue);
        _TEMP_scatterPlot_EOverTruthVsPhi_->Fill(phiFound,eFound/eTrue);

        _TEMP_scatterPlot_deltaEtaVsEt_->Fill(etFound,etaFound-etaTrue);
        _TEMP_scatterPlot_deltaEtaVsE_->Fill(eFound,etaFound-etaTrue);
        _TEMP_scatterPlot_deltaEtaVsEta_->Fill(etaFound,etaFound-etaTrue);
        _TEMP_scatterPlot_deltaEtaVsPhi_->Fill(phiFound,etaFound-etaTrue);

        _TEMP_scatterPlot_deltaPhiVsEt_->Fill(etFound,deltaPhi);
        _TEMP_scatterPlot_deltaPhiVsE_->Fill(eFound,deltaPhi);
        _TEMP_scatterPlot_deltaPhiVsEta_->Fill(etaFound,deltaPhi);
        _TEMP_scatterPlot_deltaPhiVsPhi_->Fill(phiFound,deltaPhi);

        photonsMCMatched.push_back(bestMatchPhoton);
      }

      hist_EtNumRecoOverNumTrue_->Fill(etTrue);
      hist_ENumRecoOverNumTrue_->Fill(eTrue);
      hist_EtaNumRecoOverNumTrue_->Fill(etaTrue);
      hist_PhiNumRecoOverNumTrue_->Fill(phiTrue);
    }
  }

  if(photonsMCMatched.size() == 2)
  {
    reco::Photon pOne = photonsMCMatched[0];
    reco::Photon pTwo = photonsMCMatched[1];

    double recoMass = findRecoMass(pOne, pTwo);

    hist_All_recoMass_->Fill(recoMass);

    if(pOne.superCluster()->seed()->algo() == 1 && pTwo.superCluster()->seed()->algo() == 1)
      hist_BarrelOnly_recoMass_->Fill(recoMass);
    else if(pOne.superCluster()->seed()->algo() == 0 && pTwo.superCluster()->seed()->algo() == 0)
      hist_EndcapOnly_recoMass_->Fill(recoMass);
    else
      hist_Mixed_recoMass_->Fill(recoMass);
  }
}
void EgammaObjects::applyLabels ( ) [private]

Definition at line 857 of file EgammaObjects.cc.

References hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_deltaEtaVsE_, hist_deltaEtaVsEt_, hist_deltaEtaVsEta_, hist_deltaEtaVsPhi_, hist_deltaPhiVsE_, hist_deltaPhiVsEt_, hist_deltaPhiVsEta_, hist_deltaPhiVsPhi_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_EOverTruthVsE_, hist_EOverTruthVsEt_, hist_EOverTruthVsEta_, hist_EOverTruthVsPhi_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_EtOverTruthVsE_, hist_EtOverTruthVsEt_, hist_EtOverTruthVsEta_, hist_EtOverTruthVsPhi_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, hist_resolutionEtaVsE_, hist_resolutionEtaVsEt_, hist_resolutionEtaVsEta_, hist_resolutionEtaVsPhi_, hist_resolutionEtVsE_, hist_resolutionEtVsEt_, hist_resolutionEtVsEta_, hist_resolutionEtVsPhi_, hist_resolutionEVsE_, hist_resolutionEVsEt_, hist_resolutionEVsEta_, hist_resolutionEVsPhi_, hist_resolutionPhiVsE_, hist_resolutionPhiVsEt_, hist_resolutionPhiVsEta_, hist_resolutionPhiVsPhi_, and particleID.

Referenced by endJob().

{
  hist_Et_->GetXaxis()->SetTitle("Et (GeV)");
  hist_Et_->GetYaxis()->SetTitle("# per Et Bin");
  hist_EtOverTruth_->GetXaxis()->SetTitle("Reco Et/True Et");
  hist_EtOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
  hist_EtEfficiency_->GetXaxis()->SetTitle("Et (GeV)");
  hist_EtEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per Et Bin");
  hist_EtNumRecoOverNumTrue_->GetXaxis()->SetTitle("Et (GeV)");
  hist_EtNumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per Et Bin");
  hist_EtOverTruthVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_EtOverTruthVsEt_->GetYaxis()->SetTitle("Reco Et/True Et per Et Bin");
  hist_EtOverTruthVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_EtOverTruthVsE_->GetYaxis()->SetTitle("Reco Et/True Et per E Bin");
  hist_EtOverTruthVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_EtOverTruthVsEta_->GetYaxis()->SetTitle("Reco Et/True Et per Eta Bin");
  hist_EtOverTruthVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_EtOverTruthVsPhi_->GetYaxis()->SetTitle("Reco Et/True Et per Phi Bin");
  hist_resolutionEtVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_resolutionEtVsEt_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per Et Bin");
  hist_resolutionEtVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_resolutionEtVsE_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per E Bin");
  hist_resolutionEtVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_resolutionEtVsEta_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per Eta Bin");
  hist_resolutionEtVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_resolutionEtVsPhi_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per Phi Bin");

  hist_E_->GetXaxis()->SetTitle("E (GeV)");
  hist_E_->GetYaxis()->SetTitle("# per E Bin");
  hist_EOverTruth_->GetXaxis()->SetTitle("Reco E/True E");
  hist_EOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
  hist_EEfficiency_->GetXaxis()->SetTitle("E (GeV)");
  hist_EEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per E Bin");
  hist_ENumRecoOverNumTrue_->GetXaxis()->SetTitle("E (GeV)");
  hist_ENumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per E Bin");
  hist_EOverTruthVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_EOverTruthVsEt_->GetYaxis()->SetTitle("Reco E/True E per Et Bin");
  hist_EOverTruthVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_EOverTruthVsE_->GetYaxis()->SetTitle("Reco E/True E per E Bin");
  hist_EOverTruthVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_EOverTruthVsEta_->GetYaxis()->SetTitle("Reco E/True E per Eta Bin");
  hist_EOverTruthVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_EOverTruthVsPhi_->GetYaxis()->SetTitle("Reco E/True E per Phi Bin");
  hist_resolutionEVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_resolutionEVsEt_->GetYaxis()->SetTitle("#sigma of Reco E/True E per Et Bin");
  hist_resolutionEVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_resolutionEVsE_->GetYaxis()->SetTitle("#sigma of Reco E/True E per E Bin");
  hist_resolutionEVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_resolutionEVsEta_->GetYaxis()->SetTitle("#sigma of Reco E/True E per Eta Bin");
  hist_resolutionEVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_resolutionEVsPhi_->GetYaxis()->SetTitle("#sigma of Reco E/True E per Phi Bin");

  hist_Eta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_Eta_->GetYaxis()->SetTitle("# per Eta Bin");
  hist_EtaOverTruth_->GetXaxis()->SetTitle("Reco Eta/True Eta");
  hist_EtaOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
  hist_EtaEfficiency_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_EtaEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per Eta Bin");
  hist_EtaNumRecoOverNumTrue_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_EtaNumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per Eta Bin");
  hist_deltaEtaVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_deltaEtaVsEt_->GetYaxis()->SetTitle("Reco Eta - True Eta per Et Bin");
  hist_deltaEtaVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_deltaEtaVsE_->GetYaxis()->SetTitle("Reco Eta - True Eta per E Bin");
  hist_deltaEtaVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_deltaEtaVsEta_->GetYaxis()->SetTitle("Reco Eta - True Eta per Eta Bin");
  hist_deltaEtaVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_deltaEtaVsPhi_->GetYaxis()->SetTitle("Reco Eta - True Eta per Phi Bin");
  hist_resolutionEtaVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_resolutionEtaVsEt_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per Et Bin");
  hist_resolutionEtaVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_resolutionEtaVsE_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per E Bin");
  hist_resolutionEtaVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_resolutionEtaVsEta_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per Eta Bin");
  hist_resolutionEtaVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_resolutionEtaVsPhi_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per Phi Bin");

  hist_Phi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_Phi_->GetYaxis()->SetTitle("# per Phi Bin");
  hist_PhiOverTruth_->GetXaxis()->SetTitle("Reco Phi/True Phi");
  hist_PhiOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
  hist_PhiEfficiency_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_PhiEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per Phi Bin");
  hist_PhiNumRecoOverNumTrue_->GetXaxis()->SetTitle("#Phi (Radians)");
  hist_PhiNumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per Phi Bin");
  hist_deltaPhiVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_deltaPhiVsEt_->GetYaxis()->SetTitle("Reco Phi - True Phi per Et Bin");
  hist_deltaPhiVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_deltaPhiVsE_->GetYaxis()->SetTitle("Reco Phi - True Phi per E Bin");
  hist_deltaPhiVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_deltaPhiVsEta_->GetYaxis()->SetTitle("Reco Phi - True Phi per Eta Bin");
  hist_deltaPhiVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_deltaPhiVsPhi_->GetYaxis()->SetTitle("Reco Phi - True Phi per Phi Bin");
  hist_resolutionPhiVsEt_->GetXaxis()->SetTitle("Et (GeV)");
  hist_resolutionPhiVsEt_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per Et Bin");
  hist_resolutionPhiVsE_->GetXaxis()->SetTitle("E (GeV)");
  hist_resolutionPhiVsE_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per E Bin");
  hist_resolutionPhiVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
  hist_resolutionPhiVsEta_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per Eta Bin");
  hist_resolutionPhiVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
  hist_resolutionPhiVsPhi_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per Phi Bin");

  std::string recoParticleName;

  if( particleID == 22 ) recoParticleName = "Reco Higgs";
  else if( particleID == 11 ) recoParticleName = "Reco Z";

  hist_All_recoMass_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_All_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_BarrelOnly_recoMass_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_BarrelOnly_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_EndcapOnly_recoMass_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_EndcapOnly_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_Mixed_recoMass_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_Mixed_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_recoMass_withBackgroud_NoEtCut_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_recoMass_withBackgroud_NoEtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_recoMass_withBackgroud_5EtCut_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_recoMass_withBackgroud_5EtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_recoMass_withBackgroud_10EtCut_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_recoMass_withBackgroud_10EtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
  hist_recoMass_withBackgroud_20EtCut_->GetXaxis()->SetTitle((recoParticleName+" Mass (GeV)").c_str());
  hist_recoMass_withBackgroud_20EtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
}
void EgammaObjects::beginJob ( void  ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 146 of file EgammaObjects.cc.

References createBookedHistoObjects(), and createTempHistoObjects().

{
  TH1::SetDefaultSumw2(true);

  createBookedHistoObjects();
  createTempHistoObjects();
}
void EgammaObjects::createBookedHistoObjects ( ) [private]

Definition at line 154 of file EgammaObjects.cc.

References hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_bins_E_, hist_bins_EOverTruth_, hist_bins_Et_, hist_bins_Eta_, hist_bins_EtaOverTruth_, hist_bins_EtOverTruth_, hist_bins_Phi_, hist_bins_PhiOverTruth_, hist_bins_recoMass_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_max_E_, hist_max_EOverTruth_, hist_max_Et_, hist_max_Eta_, hist_max_EtaOverTruth_, hist_max_EtOverTruth_, hist_max_Phi_, hist_max_PhiOverTruth_, hist_max_recoMass_, hist_min_E_, hist_min_EOverTruth_, hist_min_Et_, hist_min_Eta_, hist_min_EtaOverTruth_, hist_min_EtOverTruth_, hist_min_Phi_, hist_min_PhiOverTruth_, hist_min_recoMass_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, particleID, and particleString.

Referenced by beginJob().

{
  hist_Et_
    = new TH1D("hist_Et_",("Et Distribution of "+particleString).c_str(),
      hist_bins_Et_,hist_min_Et_,hist_max_Et_);
  hist_EtOverTruth_
    = new TH1D("hist_EtOverTruth_",("Reco Et over True Et of "+particleString).c_str(),
      hist_bins_EtOverTruth_,hist_min_EtOverTruth_,hist_max_EtOverTruth_);
  hist_EtEfficiency_
    = new TH1D("hist_EtEfficiency_",("# of True "+particleString+" Reconstructed over # of True "+particleString+" VS Et of "+particleString).c_str(),
      hist_bins_Et_,hist_min_Et_,hist_max_Et_);
  hist_EtNumRecoOverNumTrue_
    = new TH1D("hist_EtNumRecoOverNumTrue_",("# of Reco "+particleString+" over # of True "+particleString+" VS Et of "+particleString).c_str(),
      hist_bins_Et_,hist_min_Et_,hist_max_Et_);

  hist_E_
    = new TH1D("hist_E_",("E Distribution of "+particleString).c_str(),
      hist_bins_E_,hist_min_E_,hist_max_E_);
  hist_EOverTruth_
    = new TH1D("hist_EOverTruth_",("Reco E over True E of "+particleString).c_str(),
      hist_bins_EOverTruth_,hist_min_EOverTruth_,hist_max_EOverTruth_);
  hist_EEfficiency_
    = new TH1D("hist_EEfficiency_",("# of True "+particleString+" Reconstructed over # of True "+particleString+" VS E of "+particleString).c_str(),
      hist_bins_E_,hist_min_E_,hist_max_E_);
  hist_ENumRecoOverNumTrue_
    = new TH1D("hist_ENumRecoOverNumTrue_",("# of Reco "+particleString+" over # of True "+particleString+" VS E of "+particleString).c_str(),
      hist_bins_E_,hist_min_E_,hist_max_E_);

  hist_Eta_
    = new TH1D("hist_Eta_",("Eta Distribution of "+particleString).c_str(),
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_);
  hist_EtaOverTruth_
    = new TH1D("hist_EtaOverTruth_",("Reco Eta over True Eta of "+particleString).c_str(),
      hist_bins_EtaOverTruth_,hist_min_EtaOverTruth_,hist_max_EtaOverTruth_);
  hist_EtaEfficiency_
    = new TH1D("hist_EtaEfficiency_",("# of True "+particleString+" Reconstructed over # of True "+particleString+" VS Eta of "+particleString).c_str(),
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_);
  hist_EtaNumRecoOverNumTrue_
    = new TH1D("hist_EtaNumRecoOverNumTrue_",("# of Reco "+particleString+" over # of True "+particleString+" VS Eta of "+particleString).c_str(),
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_);

  hist_Phi_
    = new TH1D("hist_Phi_",("Phi Distribution of "+particleString).c_str(),
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_);
  hist_PhiOverTruth_
    = new TH1D("hist_PhiOverTruth_",("Reco Phi over True Phi of "+particleString).c_str(),
      hist_bins_PhiOverTruth_,hist_min_PhiOverTruth_,hist_max_PhiOverTruth_);
  hist_PhiEfficiency_
    = new TH1D("hist_PhiEfficiency_",("# of True "+particleString+" Reconstructed over # of True "+particleString+" VS Phi of "+particleString).c_str(),
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_);
   hist_PhiNumRecoOverNumTrue_
    = new TH1D("hist_PhiNumRecoOverNumTrue_",("# of Reco "+particleString+" over # of True "+particleString+" VS Phi of "+particleString).c_str(),
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_);

  std::string recoParticleName;

  if( particleID == 22 ) recoParticleName = "Higgs";
  else if( particleID == 11 ) recoParticleName = "Z";

  hist_All_recoMass_
    = new TH1D("hist_All_recoMass_",(recoParticleName+" Mass from "+particleString+" in All Regions").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
  hist_BarrelOnly_recoMass_
    = new TH1D("hist_BarrelOnly_recoMass_",(recoParticleName+" Mass from "+particleString+" in Barrel").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
  hist_EndcapOnly_recoMass_
    = new TH1D("hist_EndcapOnly_recoMass_",(recoParticleName+" Mass from "+particleString+" in EndCap").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
  hist_Mixed_recoMass_
    = new TH1D("hist_Mixed_recoMass_",(recoParticleName+" Mass from "+particleString+" in Split Detectors").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);

  hist_recoMass_withBackgroud_NoEtCut_
    = new TH1D("hist_recoMass_withBackgroud_NoEtCut_",(recoParticleName+" Mass from "+particleString+" with Background, No Et Cut").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
  hist_recoMass_withBackgroud_5EtCut_
    = new TH1D("hist_recoMass_withBackgroud_5EtCut_",(recoParticleName+" Mass from "+particleString+" with Background, 5 Et Cut").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
  hist_recoMass_withBackgroud_10EtCut_
    = new TH1D("hist_recoMass_withBackgroud_10EtCut_",(recoParticleName+" Mass from "+particleString+" with Background, 10 Et Cut").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
  hist_recoMass_withBackgroud_20EtCut_
    = new TH1D("hist_recoMass_withBackgroud_20EtCut_",(recoParticleName+" Mass from "+particleString+" with Background, 20 Et Cut").c_str(),
      hist_bins_recoMass_,hist_min_recoMass_,hist_max_recoMass_);
}
void EgammaObjects::createTempHistoObjects ( ) [private]

Definition at line 240 of file EgammaObjects.cc.

References _TEMP_scatterPlot_deltaEtaVsE_, _TEMP_scatterPlot_deltaEtaVsEt_, _TEMP_scatterPlot_deltaEtaVsEta_, _TEMP_scatterPlot_deltaEtaVsPhi_, _TEMP_scatterPlot_deltaPhiVsE_, _TEMP_scatterPlot_deltaPhiVsEt_, _TEMP_scatterPlot_deltaPhiVsEta_, _TEMP_scatterPlot_deltaPhiVsPhi_, _TEMP_scatterPlot_EOverTruthVsE_, _TEMP_scatterPlot_EOverTruthVsEt_, _TEMP_scatterPlot_EOverTruthVsEta_, _TEMP_scatterPlot_EOverTruthVsPhi_, _TEMP_scatterPlot_EtOverTruthVsE_, _TEMP_scatterPlot_EtOverTruthVsEt_, _TEMP_scatterPlot_EtOverTruthVsEta_, _TEMP_scatterPlot_EtOverTruthVsPhi_, hist_bins_deltaEta_, hist_bins_deltaPhi_, hist_bins_E_, hist_bins_EOverTruth_, hist_bins_Et_, hist_bins_Eta_, hist_bins_EtOverTruth_, hist_bins_Phi_, hist_max_deltaEta_, hist_max_deltaPhi_, hist_max_E_, hist_max_EOverTruth_, hist_max_Et_, hist_max_Eta_, hist_max_EtOverTruth_, hist_max_Phi_, hist_min_deltaEta_, hist_min_deltaPhi_, hist_min_E_, hist_min_EOverTruth_, hist_min_Et_, hist_min_Eta_, hist_min_EtOverTruth_, and hist_min_Phi_.

Referenced by beginJob().

{
  _TEMP_scatterPlot_EtOverTruthVsEt_
    = new TH2D("_TEMP_scatterPlot_EtOverTruthVsEt_","_TEMP_scatterPlot_EtOverTruthVsEt_",
      hist_bins_Et_,hist_min_Et_,hist_max_Et_,hist_bins_EtOverTruth_,hist_min_EtOverTruth_,hist_max_EtOverTruth_);
  _TEMP_scatterPlot_EtOverTruthVsE_
    = new TH2D("_TEMP_scatterPlot_EtOverTruthVsE_","_TEMP_scatterPlot_EtOverTruthVsE_",
      hist_bins_E_,hist_min_E_,hist_max_E_,hist_bins_EtOverTruth_,hist_min_EtOverTruth_,hist_max_EtOverTruth_);
  _TEMP_scatterPlot_EtOverTruthVsEta_
    = new TH2D("_TEMP_scatterPlot_EtOverTruthVsEta_","_TEMP_scatterPlot_EtOverTruthVsEta_",
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_,hist_bins_EtOverTruth_,hist_min_EtOverTruth_,hist_max_EtOverTruth_);
  _TEMP_scatterPlot_EtOverTruthVsPhi_
    = new TH2D("_TEMP_scatterPlot_EtOverTruthVsPhi_","_TEMP_scatterPlot_EtOverTruthVsPhi_",
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_,hist_bins_EtOverTruth_,hist_min_EtOverTruth_,hist_max_EtOverTruth_);

  _TEMP_scatterPlot_EOverTruthVsEt_
    = new TH2D("_TEMP_scatterPlot_EOverTruthVsEt_","_TEMP_scatterPlot_EOverTruthVsEt_",
      hist_bins_Et_,hist_min_Et_,hist_max_Et_,hist_bins_EOverTruth_,hist_min_EOverTruth_,hist_max_EOverTruth_);
  _TEMP_scatterPlot_EOverTruthVsE_
    = new TH2D("_TEMP_scatterPlot_EOverTruthVsE_","_TEMP_scatterPlot_EOverTruthVsE_",
      hist_bins_E_,hist_min_E_,hist_max_E_,hist_bins_EOverTruth_,hist_min_EOverTruth_,hist_max_EOverTruth_);
  _TEMP_scatterPlot_EOverTruthVsEta_
    = new TH2D("_TEMP_scatterPlot_EOverTruthVsEta_","_TEMP_scatterPlot_EOverTruthVsEta_",
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_,hist_bins_EOverTruth_,hist_min_EOverTruth_,hist_max_EOverTruth_);
  _TEMP_scatterPlot_EOverTruthVsPhi_
    = new TH2D("_TEMP_scatterPlot_EOverTruthVsPhi_","_TEMP_scatterPlot_EOverTruthVsPhi_",
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_,hist_bins_EOverTruth_,hist_min_EOverTruth_,hist_max_EOverTruth_);

  _TEMP_scatterPlot_deltaEtaVsEt_
    = new TH2D("_TEMP_scatterPlot_deltaEtaVsEt_","_TEMP_scatterPlot_deltaEtaVsEt_",
      hist_bins_Et_,hist_min_Et_,hist_max_Et_,hist_bins_deltaEta_,hist_min_deltaEta_,hist_max_deltaEta_);
  _TEMP_scatterPlot_deltaEtaVsE_
    = new TH2D("_TEMP_scatterPlot_deltaEtaVsE_","_TEMP_scatterPlot_deltaEtaVsE_",
      hist_bins_E_,hist_min_E_,hist_max_E_,hist_bins_deltaEta_,hist_min_deltaEta_,hist_max_deltaEta_);
  _TEMP_scatterPlot_deltaEtaVsEta_
    = new TH2D("_TEMP_scatterPlot_deltaEtaVsEta_","_TEMP_scatterPlot_deltaEtaVsEta_",
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_,hist_bins_deltaEta_,hist_min_deltaEta_,hist_max_deltaEta_);
  _TEMP_scatterPlot_deltaEtaVsPhi_
    = new TH2D("_TEMP_scatterPlot_deltaEtaVsPhi_","_TEMP_scatterPlot_deltaEtaVsPhi_",
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_,hist_bins_deltaEta_,hist_min_deltaEta_,hist_max_deltaEta_);

  _TEMP_scatterPlot_deltaPhiVsEt_
    = new TH2D("_TEMP_scatterPlot_deltaPhiVsEt_","_TEMP_scatterPlot_deltaPhiVsEt_",
      hist_bins_Et_,hist_min_Et_,hist_max_Et_,hist_bins_deltaPhi_,hist_min_deltaPhi_,hist_max_deltaPhi_);
  _TEMP_scatterPlot_deltaPhiVsE_
    = new TH2D("_TEMP_scatterPlot_deltaPhiVsE_","_TEMP_scatterPlot_deltaPhiVsE_",
      hist_bins_E_,hist_min_E_,hist_max_E_,hist_bins_deltaPhi_,hist_min_deltaPhi_,hist_max_deltaPhi_);
  _TEMP_scatterPlot_deltaPhiVsEta_
    = new TH2D("_TEMP_scatterPlot_deltaPhiVsEta_","_TEMP_scatterPlot_deltaPhiVsEta_",
      hist_bins_Eta_,hist_min_Eta_,hist_max_Eta_,hist_bins_deltaPhi_,hist_min_deltaPhi_,hist_max_deltaPhi_);
  _TEMP_scatterPlot_deltaPhiVsPhi_
    = new TH2D("_TEMP_scatterPlot_deltaPhiVsPhi_","_TEMP_scatterPlot_deltaPhiVsPhi_",
      hist_bins_Phi_,hist_min_Phi_,hist_max_Phi_,hist_bins_deltaPhi_,hist_min_deltaPhi_,hist_max_deltaPhi_);
}
float EgammaObjects::ecalEta ( float  EtaParticle,
float  Zvertex,
float  plane_Radius 
) [private]

Definition at line 644 of file EgammaObjects.cc.

References abs, ETA, etaBarrelEndcap, funct::log(), pi, R_ECAL, funct::tan(), and Z_Endcap.

Referenced by analyzePhotons().

{
        const float R_ECAL           = 136.5;
        const float Z_Endcap         = 328.0;
        const float etaBarrelEndcap  = 1.479;

        if(EtaParticle != 0.)
        {
                float Theta = 0.0  ;
                float ZEcal = (R_ECAL-plane_Radius)*sinh(EtaParticle)+Zvertex;

                if(ZEcal != 0.0) Theta = atan(R_ECAL/ZEcal);
                if(Theta<0.0) Theta = Theta+Geom::pi() ;

                float ETA = - log(tan(0.5*Theta));

                if( std::abs(ETA) > etaBarrelEndcap )
                {
                        float Zend = Z_Endcap ;
                        if(EtaParticle<0.0 )  Zend = -Zend ;
                        float Zlen = Zend - Zvertex ;
                        float RR = Zlen/sinh(EtaParticle);
                        Theta = atan((RR+plane_Radius)/Zend);
                        if(Theta<0.0) Theta = Theta+Geom::pi() ;
                        ETA = - log(tan(0.5*Theta));
                }

                return ETA;
        }
        else
        {
                edm::LogWarning("")  << "[EgammaObjects::ecalEta] Warning: Eta equals to zero, not correcting" ;
                return EtaParticle;
        }
}
void EgammaObjects::endJob ( void  ) [virtual]
double EgammaObjects::findRecoMass ( reco::Photon  pOne,
reco::Photon  pTwo 
) [private]

Definition at line 622 of file EgammaObjects.cc.

References funct::cos(), relval_parameters_module::energy, mathSSE::sqrt(), and reco::Photon::superCluster().

Referenced by analyzeElectrons(), and analyzePhotons().

{
  double cosTheta
    = (cos(pOne.superCluster()->phi() - pTwo.superCluster()->phi()) + sinh(pOne.superCluster()->eta()) * sinh(pTwo.superCluster()->eta())) /
      (cosh(pOne.superCluster()->eta()) * cosh(pTwo.superCluster()->eta()));

  double recoMass = sqrt(2 * (pOne.superCluster())->energy() * (pTwo.superCluster())->energy() * (1 - cosTheta));

  return recoMass;
}
double EgammaObjects::findRecoMass ( reco::GsfElectron  eOne,
reco::GsfElectron  eTwo 
) [private]

Definition at line 633 of file EgammaObjects.cc.

References reco::GsfElectron::caloEnergy(), reco::GsfElectron::caloPosition(), funct::cos(), and mathSSE::sqrt().

{
  double cosTheta
    = (cos(eOne.caloPosition().phi() - eTwo.caloPosition().phi()) + sinh(eOne.caloPosition().eta()) * sinh(eTwo.caloPosition().eta())) /
      (cosh(eOne.caloPosition().eta()) * cosh(eTwo.caloPosition().eta()));

  double recoMass = sqrt(2 * eOne.caloEnergy() * eTwo.caloEnergy() * (1 - cosTheta));

  return recoMass;
}
void EgammaObjects::fitHistos ( ) [private]

Definition at line 812 of file EgammaObjects.cc.

References hist_EOverTruth_, hist_EtaOverTruth_, hist_EtOverTruth_, and hist_PhiOverTruth_.

Referenced by endJob().

{
  hist_EtOverTruth_->Fit("gaus","QEM");
//  hist_EtNumRecoOverNumTrue_->Fit("pol1","QEM");

  hist_EOverTruth_->Fit("gaus","QEM");
//  hist_ENumRecoOverNumTrue_->Fit("pol1","QEM");

  hist_EtaOverTruth_->Fit("gaus","QEM");
//  hist_EtaNumRecoOverNumTrue_->Fit("pol1","QEM");

  hist_PhiOverTruth_->Fit("gaus","QEM");
//  hist_PhiNumRecoOverNumTrue_->Fit("pol1","QEM");

  /*
  hist_EtOverTruthVsEt_->Fit("pol1","QEM");
  hist_EtOverTruthVsEta_->Fit("pol1","QEM");
  hist_EtOverTruthVsPhi_->Fit("pol1","QEM");
  hist_resolutionEtVsEt_->Fit("pol1","QEM");
  hist_resolutionEtVsEta_->Fit("pol1","QEM");
  hist_resolutionEtVsPhi_->Fit("pol1","QEM");

  hist_EOverTruthVsEt_->Fit("pol1","QEM");
  hist_EOverTruthVsEta_->Fit("pol1","QEM");
  hist_EOverTruthVsPhi_->Fit("pol1","QEM");
  hist_resolutionEVsEt_->Fit("pol1","QEM");
  hist_resolutionEVsEta_->Fit("pol1","QEM");
  hist_resolutionEVsPhi_->Fit("pol1","QEM");

  hist_deltaEtaVsEt_->Fit("pol1","QEM");
  hist_deltaEtaVsEta_->Fit("pol1","QEM");
  hist_deltaEtaVsPhi_->Fit("pol1","QEM");
  hist_resolutionEtaVsEt_->Fit("pol1","QEM");
  hist_resolutionEtaVsEta_->Fit("pol1","QEM");
  hist_resolutionEtaVsPhi_->Fit("pol1","QEM");

  hist_deltaPhiVsEt_->Fit("pol1","QEM");
  hist_deltaPhiVsEta_->Fit("pol1","QEM");
  hist_deltaPhiVsPhi_->Fit("pol1","QEM");
  hist_resolutionPhiVsEt_->Fit("pol1","QEM");
  hist_resolutionPhiVsEta_->Fit("pol1","QEM");
  hist_resolutionPhiVsPhi_->Fit("pol1","QEM");
  */
}
void EgammaObjects::getDeltaResHistosViaSlicing ( ) [private]

Definition at line 696 of file EgammaObjects.cc.

References _TEMP_scatterPlot_deltaEtaVsE_, _TEMP_scatterPlot_deltaEtaVsEt_, _TEMP_scatterPlot_deltaEtaVsEta_, _TEMP_scatterPlot_deltaEtaVsPhi_, _TEMP_scatterPlot_deltaPhiVsE_, _TEMP_scatterPlot_deltaPhiVsEt_, _TEMP_scatterPlot_deltaPhiVsEta_, _TEMP_scatterPlot_deltaPhiVsPhi_, _TEMP_scatterPlot_EOverTruthVsE_, _TEMP_scatterPlot_EOverTruthVsEt_, _TEMP_scatterPlot_EOverTruthVsEta_, _TEMP_scatterPlot_EOverTruthVsPhi_, _TEMP_scatterPlot_EtOverTruthVsE_, _TEMP_scatterPlot_EtOverTruthVsEt_, _TEMP_scatterPlot_EtOverTruthVsEta_, _TEMP_scatterPlot_EtOverTruthVsPhi_, hist_bins_E_, hist_bins_Et_, hist_bins_Eta_, hist_bins_Phi_, hist_deltaEtaVsE_, hist_deltaEtaVsEt_, hist_deltaEtaVsEta_, hist_deltaEtaVsPhi_, hist_deltaPhiVsE_, hist_deltaPhiVsEt_, hist_deltaPhiVsEta_, hist_deltaPhiVsPhi_, hist_EOverTruthVsE_, hist_EOverTruthVsEt_, hist_EOverTruthVsEta_, hist_EOverTruthVsPhi_, hist_EtOverTruthVsE_, hist_EtOverTruthVsEt_, hist_EtOverTruthVsEta_, hist_EtOverTruthVsPhi_, hist_resolutionEtaVsE_, hist_resolutionEtaVsEt_, hist_resolutionEtaVsEta_, hist_resolutionEtaVsPhi_, hist_resolutionEtVsE_, hist_resolutionEtVsEt_, hist_resolutionEtVsEta_, hist_resolutionEtVsPhi_, hist_resolutionEVsE_, hist_resolutionEVsEt_, hist_resolutionEVsEta_, hist_resolutionEVsPhi_, hist_resolutionPhiVsE_, hist_resolutionPhiVsEt_, hist_resolutionPhiVsEta_, hist_resolutionPhiVsPhi_, and particleString.

Referenced by endJob().

{
  _TEMP_scatterPlot_EtOverTruthVsEt_->FitSlicesY(0,1,hist_bins_Et_,10,"QRG3");
  _TEMP_scatterPlot_EtOverTruthVsE_->FitSlicesY(0,1,hist_bins_E_,10,"QRG3");
  _TEMP_scatterPlot_EtOverTruthVsEta_->FitSlicesY(0,1,hist_bins_Eta_,10,"QRG2");
  _TEMP_scatterPlot_EtOverTruthVsPhi_->FitSlicesY(0,1,hist_bins_Phi_,10,"QRG2");

  _TEMP_scatterPlot_EOverTruthVsEt_->FitSlicesY(0,1,hist_bins_Et_,10,"QRG3");
  _TEMP_scatterPlot_EOverTruthVsE_->FitSlicesY(0,1,hist_bins_E_,10,"QRG3");
  _TEMP_scatterPlot_EOverTruthVsEta_->FitSlicesY(0,1,hist_bins_Eta_,10,"QRG2");
  _TEMP_scatterPlot_EOverTruthVsPhi_->FitSlicesY(0,1,hist_bins_Phi_,10,"QRG2");

  _TEMP_scatterPlot_deltaEtaVsEt_->FitSlicesY(0,1,hist_bins_Et_,10,"QRG3");
  _TEMP_scatterPlot_deltaEtaVsE_->FitSlicesY(0,1,hist_bins_E_,10,"QRG3");
  _TEMP_scatterPlot_deltaEtaVsEta_->FitSlicesY(0,1,hist_bins_Eta_,10,"QRG2");
  _TEMP_scatterPlot_deltaEtaVsPhi_->FitSlicesY(0,1,hist_bins_Phi_,10,"QRG2");

  _TEMP_scatterPlot_deltaPhiVsEt_->FitSlicesY(0,1,hist_bins_Et_,10,"QRG3");
  _TEMP_scatterPlot_deltaPhiVsE_->FitSlicesY(0,1,hist_bins_E_,10,"QRG3");
  _TEMP_scatterPlot_deltaPhiVsEta_->FitSlicesY(0,1,hist_bins_Eta_,10,"QRG2");
  _TEMP_scatterPlot_deltaPhiVsPhi_->FitSlicesY(0,1,hist_bins_Phi_,10,"QRG2");

  hist_EtOverTruthVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEt__1");
  hist_EtOverTruthVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsE__1");
  hist_EtOverTruthVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEta__1");
  hist_EtOverTruthVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsPhi__1");

  hist_EOverTruthVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEt__1");
  hist_EOverTruthVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsE__1");
  hist_EOverTruthVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEta__1");
  hist_EOverTruthVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsPhi__1");

  hist_deltaEtaVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEt__1");
  hist_deltaEtaVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsE__1");
  hist_deltaEtaVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEta__1");
  hist_deltaEtaVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsPhi__1");

  hist_deltaPhiVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEt__1");
  hist_deltaPhiVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsE__1");
  hist_deltaPhiVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEta__1");
  hist_deltaPhiVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsPhi__1");

  hist_EtOverTruthVsEt_->SetNameTitle("hist_EtOverTruthVsEt_",("Reco Et over True Et VS Et of "+particleString).c_str());
  hist_EtOverTruthVsE_->SetNameTitle("hist_EtOverTruthVsE_",("Reco Et over True Et VS E of "+particleString).c_str());
  hist_EtOverTruthVsEta_->SetNameTitle("hist_EtOverTruthVsEta_",("Reco Et over True Et VS Eta of "+particleString).c_str());
  hist_EtOverTruthVsPhi_->SetNameTitle("hist_EtOverTruthVsPhi_",("Reco Et over True Et VS Phi of "+particleString).c_str());

  hist_EOverTruthVsEt_->SetNameTitle("hist_EOverTruthVsEt_",("Reco E over True E VS Et of "+particleString).c_str());
  hist_EOverTruthVsE_->SetNameTitle("hist_EOverTruthVsE_",("Reco E over True E VS E of "+particleString).c_str());
  hist_EOverTruthVsEta_->SetNameTitle("hist_EOverTruthVsEta_",("Reco E over True E VS Eta of "+particleString).c_str());
  hist_EOverTruthVsPhi_->SetNameTitle("hist_EOverTruthVsPhi_",("Reco E over True E VS Phi of "+particleString).c_str());

  hist_deltaEtaVsEt_->SetNameTitle("hist_deltaEtaVsEt_",("delta Eta VS Et of "+particleString).c_str());
  hist_deltaEtaVsE_->SetNameTitle("hist_deltaEtaVsE_",("delta Eta VS E of "+particleString).c_str());
  hist_deltaEtaVsEta_->SetNameTitle("hist_deltaEtaVsEta_",("delta Eta VS Eta of "+particleString).c_str());
  hist_deltaEtaVsPhi_->SetNameTitle("hist_deltaEtaVsPhi_",("delta Eta VS Phi of "+particleString).c_str());

  hist_deltaPhiVsEt_->SetNameTitle("hist_deltaPhiVsEt_",("delta Phi VS Et of "+particleString).c_str());
  hist_deltaPhiVsE_->SetNameTitle("hist_deltaPhiVsE_",("delta Phi VS E of "+particleString).c_str());
  hist_deltaPhiVsEta_->SetNameTitle("hist_deltaPhiVsEta_",("delta Phi VS Eta of "+particleString).c_str());
  hist_deltaPhiVsPhi_->SetNameTitle("hist_deltaPhiVsPhi_",("delta Phi VS Phi of "+particleString).c_str());

  hist_resolutionEtVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEt__2");
  hist_resolutionEtVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsE__2");
  hist_resolutionEtVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEta__2");
  hist_resolutionEtVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsPhi__2");

  hist_resolutionEVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEt__2");
  hist_resolutionEVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsE__2");
  hist_resolutionEVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEta__2");
  hist_resolutionEVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsPhi__2");

  hist_resolutionEtaVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEt__2");
  hist_resolutionEtaVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsE__2");
  hist_resolutionEtaVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEta__2");
  hist_resolutionEtaVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsPhi__2");

  hist_resolutionPhiVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEt__2");
  hist_resolutionPhiVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsE__2");
  hist_resolutionPhiVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEta__2");
  hist_resolutionPhiVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsPhi__2");

  hist_resolutionEtVsEt_->SetNameTitle("hist_resolutionEtVsEt_",("#sigma of Reco Et over True Et VS Et of "+particleString).c_str());
  hist_resolutionEtVsE_->SetNameTitle("hist_resolutionEtVsE_",("#sigma of Reco Et over True Et VS E of "+particleString).c_str());
  hist_resolutionEtVsEta_->SetNameTitle("hist_resolutionEtVsEta_",("#sigma of Reco Et over True Et VS Eta of "+particleString).c_str());
  hist_resolutionEtVsPhi_->SetNameTitle("hist_resolutionEtVsPhi_",("#sigma of Reco Et over True Et VS Phi of "+particleString).c_str());

  hist_resolutionEVsEt_->SetNameTitle("hist_resolutionEVsEt_",("#sigma of Reco E over True E VS Et of "+particleString).c_str());
  hist_resolutionEVsE_->SetNameTitle("hist_resolutionEVsE_",("#sigma of Reco E over True E VS E of "+particleString).c_str());
  hist_resolutionEVsEta_->SetNameTitle("hist_resolutionEVsEta_",("#sigma of Reco E over True E VS Eta of "+particleString).c_str());
  hist_resolutionEVsPhi_->SetNameTitle("hist_resolutionEVsPhi_",("#sigma of Reco E over True E VS Phi of "+particleString).c_str());

  hist_resolutionEtaVsEt_->SetNameTitle("hist_resolutionEtaVsEt_",("#sigma of delta Eta VS Et of "+particleString).c_str());
  hist_resolutionEtaVsE_->SetNameTitle("hist_resolutionEtaVsE_",("#sigma of delta Eta VS E of "+particleString).c_str());
  hist_resolutionEtaVsEta_->SetNameTitle("hist_resolutionEtaVsEta_",("#sigma of delta Eta VS Eta of "+particleString).c_str());
  hist_resolutionEtaVsPhi_->SetNameTitle("hist_resolutionEtaVsPhi_",("#sigma of delta Eta VS Phi of "+particleString).c_str());

  hist_resolutionPhiVsEt_->SetNameTitle("hist_resolutionPhiVsEt_",("#sigma of delta Phi VS Et of "+particleString).c_str());
  hist_resolutionPhiVsE_->SetNameTitle("hist_resolutionPhiVsE_",("#sigma of delta Phi VS E of "+particleString).c_str());
  hist_resolutionPhiVsEta_->SetNameTitle("hist_resolutionPhiVsEta_",("#sigma of delta Phi VS Eta of "+particleString).c_str());
  hist_resolutionPhiVsPhi_->SetNameTitle("hist_resolutionPhiVsPhi_",("#sigma of delta Phi VS Phi of "+particleString).c_str());
}
void EgammaObjects::getEfficiencyHistosViaDividing ( ) [private]
void EgammaObjects::loadCMSSWObjects ( const edm::ParameterSet ps) [private]

Definition at line 88 of file EgammaObjects.cc.

References edm::ParameterSet::getParameter(), MCTruthCollection_, and RecoCollection_.

Referenced by EgammaObjects().

{
  MCTruthCollection_ = ps.getParameter<edm::InputTag>("MCTruthCollection");
  RecoCollection_  = ps.getParameter<edm::InputTag>("RecoCollection");
}
void EgammaObjects::loadHistoParameters ( const edm::ParameterSet ps) [private]

Definition at line 94 of file EgammaObjects.cc.

References edm::ParameterSet::getParameter(), hist_bins_deltaEta_, hist_bins_deltaPhi_, hist_bins_E_, hist_bins_EOverTruth_, hist_bins_Et_, hist_bins_Eta_, hist_bins_EtaOverTruth_, hist_bins_EtOverTruth_, hist_bins_Phi_, hist_bins_PhiOverTruth_, hist_bins_recoMass_, hist_max_deltaEta_, hist_max_deltaPhi_, hist_max_E_, hist_max_EOverTruth_, hist_max_Et_, hist_max_Eta_, hist_max_EtaOverTruth_, hist_max_EtOverTruth_, hist_max_Phi_, hist_max_PhiOverTruth_, hist_max_recoMass_, hist_min_deltaEta_, hist_min_deltaPhi_, hist_min_E_, hist_min_EOverTruth_, hist_min_Et_, hist_min_Eta_, hist_min_EtaOverTruth_, hist_min_EtOverTruth_, hist_min_Phi_, hist_min_PhiOverTruth_, and hist_min_recoMass_.

Referenced by EgammaObjects().

{
  hist_min_Et_  = ps.getParameter<double>("hist_min_Et");
  hist_max_Et_  = ps.getParameter<double>("hist_max_Et");
  hist_bins_Et_ = ps.getParameter<int>   ("hist_bins_Et");

  hist_min_E_  = ps.getParameter<double>("hist_min_E");
  hist_max_E_  = ps.getParameter<double>("hist_max_E");
  hist_bins_E_ = ps.getParameter<int>   ("hist_bins_E");

  hist_min_Eta_  = ps.getParameter<double>("hist_min_Eta");
  hist_max_Eta_  = ps.getParameter<double>("hist_max_Eta");
  hist_bins_Eta_ = ps.getParameter<int>   ("hist_bins_Eta");

  hist_min_Phi_  = ps.getParameter<double>("hist_min_Phi");
  hist_max_Phi_  = ps.getParameter<double>("hist_max_Phi");
  hist_bins_Phi_ = ps.getParameter<int>   ("hist_bins_Phi");

  hist_min_EtOverTruth_  = ps.getParameter<double>("hist_min_EtOverTruth");
  hist_max_EtOverTruth_  = ps.getParameter<double>("hist_max_EtOverTruth");
  hist_bins_EtOverTruth_ = ps.getParameter<int>   ("hist_bins_EtOverTruth");

  hist_min_EOverTruth_  = ps.getParameter<double>("hist_min_EOverTruth");
  hist_max_EOverTruth_  = ps.getParameter<double>("hist_max_EOverTruth");
  hist_bins_EOverTruth_ = ps.getParameter<int>   ("hist_bins_EOverTruth");

  hist_min_EtaOverTruth_  = ps.getParameter<double>("hist_min_EtaOverTruth");
  hist_max_EtaOverTruth_  = ps.getParameter<double>("hist_max_EtaOverTruth");
  hist_bins_EtaOverTruth_ = ps.getParameter<int>   ("hist_bins_EtaOverTruth");

  hist_min_PhiOverTruth_  = ps.getParameter<double>("hist_min_PhiOverTruth");
  hist_max_PhiOverTruth_  = ps.getParameter<double>("hist_max_PhiOverTruth");
  hist_bins_PhiOverTruth_ = ps.getParameter<int>   ("hist_bins_PhiOverTruth");

  hist_min_deltaEta_  = ps.getParameter<double>("hist_min_deltaEta");
  hist_max_deltaEta_  = ps.getParameter<double>("hist_max_deltaEta");
  hist_bins_deltaEta_ = ps.getParameter<int>   ("hist_bins_deltaEta");

  hist_min_deltaPhi_  = ps.getParameter<double>("hist_min_deltaPhi");
  hist_max_deltaPhi_  = ps.getParameter<double>("hist_max_deltaPhi");
  hist_bins_deltaPhi_ = ps.getParameter<int>   ("hist_bins_deltaPhi");

  hist_min_recoMass_ = ps.getParameter<double>("hist_min_recoMass");
  hist_max_recoMass_ = ps.getParameter<double>("hist_max_recoMass");
  hist_bins_recoMass_ = ps.getParameter<int>   ("hist_bins_recoMass");
}
void EgammaObjects::saveHistos ( ) [private]

Definition at line 1046 of file EgammaObjects.cc.

References _TEMP_scatterPlot_deltaEtaVsE_, _TEMP_scatterPlot_deltaEtaVsEt_, _TEMP_scatterPlot_deltaEtaVsEta_, _TEMP_scatterPlot_deltaEtaVsPhi_, _TEMP_scatterPlot_deltaPhiVsE_, _TEMP_scatterPlot_deltaPhiVsEt_, _TEMP_scatterPlot_deltaPhiVsEta_, _TEMP_scatterPlot_deltaPhiVsPhi_, _TEMP_scatterPlot_EOverTruthVsE_, _TEMP_scatterPlot_EOverTruthVsEt_, _TEMP_scatterPlot_EOverTruthVsEta_, _TEMP_scatterPlot_EOverTruthVsPhi_, _TEMP_scatterPlot_EtOverTruthVsE_, _TEMP_scatterPlot_EtOverTruthVsEt_, _TEMP_scatterPlot_EtOverTruthVsEta_, _TEMP_scatterPlot_EtOverTruthVsPhi_, hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_deltaEtaVsE_, hist_deltaEtaVsEt_, hist_deltaEtaVsEta_, hist_deltaEtaVsPhi_, hist_deltaPhiVsE_, hist_deltaPhiVsEt_, hist_deltaPhiVsEta_, hist_deltaPhiVsPhi_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_EOverTruthVsE_, hist_EOverTruthVsEt_, hist_EOverTruthVsEta_, hist_EOverTruthVsPhi_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_EtOverTruthVsE_, hist_EtOverTruthVsEt_, hist_EtOverTruthVsEta_, hist_EtOverTruthVsPhi_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, hist_resolutionEtaVsE_, hist_resolutionEtaVsEt_, hist_resolutionEtaVsEta_, hist_resolutionEtaVsPhi_, hist_resolutionEtVsE_, hist_resolutionEtVsEt_, hist_resolutionEtVsEta_, hist_resolutionEtVsPhi_, hist_resolutionEVsE_, hist_resolutionEVsEt_, hist_resolutionEVsEta_, hist_resolutionEVsPhi_, hist_resolutionPhiVsE_, hist_resolutionPhiVsEt_, hist_resolutionPhiVsEta_, hist_resolutionPhiVsPhi_, particleID, particleString, and rootFile_.

Referenced by endJob().

{
  rootFile_->cd();
  rootFile_->GetDirectory(particleString.c_str())->mkdir("ET");
  rootFile_->cd(("/"+particleString+"/ET").c_str());

  hist_Et_->Write();
  hist_EtOverTruth_->Write();
  hist_EtEfficiency_->Write();
  hist_EtNumRecoOverNumTrue_->Write();
  hist_EtOverTruthVsEt_->Write();
  hist_EtOverTruthVsE_->Write();
  hist_EtOverTruthVsEta_->Write();
  hist_EtOverTruthVsPhi_->Write();
  hist_resolutionEtVsEt_->Write();
  hist_resolutionEtVsE_->Write();
  hist_resolutionEtVsEta_->Write();
  hist_resolutionEtVsPhi_->Write();

  rootFile_->cd();
  rootFile_->GetDirectory(particleString.c_str())->mkdir("E");
  rootFile_->cd(("/"+particleString+"/E").c_str());

  hist_E_->Write();
  hist_EOverTruth_->Write();
  hist_EEfficiency_->Write();
  hist_ENumRecoOverNumTrue_->Write();
  hist_EOverTruthVsEt_->Write();
  hist_EOverTruthVsE_->Write();
  hist_EOverTruthVsEta_->Write();
  hist_EOverTruthVsPhi_->Write();
  hist_resolutionEVsEt_->Write();
  hist_resolutionEVsE_->Write();
  hist_resolutionEVsEta_->Write();
  hist_resolutionEVsPhi_->Write();

  rootFile_->cd();
  rootFile_->GetDirectory(particleString.c_str())->mkdir("Eta");
  rootFile_->cd(("/"+particleString+"/Eta").c_str());

  hist_Eta_->Write();
  hist_EtaOverTruth_->Write();
  hist_EtaEfficiency_->Write();
  hist_EtaNumRecoOverNumTrue_->Write();
  hist_deltaEtaVsEt_->Write();
  hist_deltaEtaVsE_->Write();
  hist_deltaEtaVsEta_->Write();
  hist_deltaEtaVsPhi_->Write();
  hist_resolutionEtaVsEt_->Write();
  hist_resolutionEtaVsE_->Write();
  hist_resolutionEtaVsEta_->Write();
  hist_resolutionEtaVsPhi_->Write();

  rootFile_->cd();
  rootFile_->GetDirectory(particleString.c_str())->mkdir("Phi");
  rootFile_->cd(("/"+particleString+"/Phi").c_str());

  hist_Phi_->Write();
  hist_PhiOverTruth_->Write();
  hist_PhiEfficiency_->Write();
  hist_PhiNumRecoOverNumTrue_->Write();
  hist_deltaPhiVsEt_->Write();
  hist_deltaPhiVsE_->Write();
  hist_deltaPhiVsEta_->Write();
  hist_deltaPhiVsPhi_->Write();
  hist_resolutionPhiVsEt_->Write();
  hist_resolutionPhiVsE_->Write();
  hist_resolutionPhiVsEta_->Write();
  hist_resolutionPhiVsPhi_->Write();

  std::string recoParticleName;

  if( particleID == 22 ) recoParticleName = "HiggsRecoMass";
  else if( particleID == 11 ) recoParticleName = "ZRecoMass";

  rootFile_->cd();
  rootFile_->GetDirectory(particleString.c_str())->mkdir(recoParticleName.c_str());
  rootFile_->cd(("/"+particleString+"/"+recoParticleName).c_str());

  hist_All_recoMass_->Write();
  hist_BarrelOnly_recoMass_->Write();
  hist_EndcapOnly_recoMass_->Write();
  hist_Mixed_recoMass_->Write();
  hist_recoMass_withBackgroud_NoEtCut_->Write();
  hist_recoMass_withBackgroud_5EtCut_->Write();
  hist_recoMass_withBackgroud_10EtCut_->Write();
  hist_recoMass_withBackgroud_20EtCut_->Write();

  rootFile_->cd();
  rootFile_->GetDirectory(particleString.c_str())->mkdir("_TempScatterPlots");
  rootFile_->cd(("/"+particleString+"/_TempScatterPlots").c_str());

  _TEMP_scatterPlot_EtOverTruthVsEt_->Write();
  _TEMP_scatterPlot_EtOverTruthVsE_->Write();
  _TEMP_scatterPlot_EtOverTruthVsEta_->Write();
  _TEMP_scatterPlot_EtOverTruthVsPhi_->Write();

  _TEMP_scatterPlot_EOverTruthVsEt_->Write();
  _TEMP_scatterPlot_EOverTruthVsE_->Write();
  _TEMP_scatterPlot_EOverTruthVsEta_->Write();
  _TEMP_scatterPlot_EOverTruthVsPhi_->Write();

  _TEMP_scatterPlot_deltaEtaVsEt_->Write();
  _TEMP_scatterPlot_deltaEtaVsE_->Write();
  _TEMP_scatterPlot_deltaEtaVsEta_->Write();
  _TEMP_scatterPlot_deltaEtaVsPhi_->Write();

  _TEMP_scatterPlot_deltaPhiVsEt_->Write();
  _TEMP_scatterPlot_deltaPhiVsE_->Write();
  _TEMP_scatterPlot_deltaPhiVsEta_->Write();
  _TEMP_scatterPlot_deltaPhiVsPhi_->Write();

  rootFile_->cd();
}
void EgammaObjects::setDrawOptions ( ) [private]

Definition at line 982 of file EgammaObjects.cc.

References hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_deltaEtaVsE_, hist_deltaEtaVsEt_, hist_deltaEtaVsEta_, hist_deltaEtaVsPhi_, hist_deltaPhiVsE_, hist_deltaPhiVsEt_, hist_deltaPhiVsEta_, hist_deltaPhiVsPhi_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_EOverTruthVsE_, hist_EOverTruthVsEt_, hist_EOverTruthVsEta_, hist_EOverTruthVsPhi_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_EtOverTruthVsE_, hist_EtOverTruthVsEt_, hist_EtOverTruthVsEta_, hist_EtOverTruthVsPhi_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, hist_resolutionEtaVsE_, hist_resolutionEtaVsEt_, hist_resolutionEtaVsEta_, hist_resolutionEtaVsPhi_, hist_resolutionEtVsE_, hist_resolutionEtVsEt_, hist_resolutionEtVsEta_, hist_resolutionEtVsPhi_, hist_resolutionEVsE_, hist_resolutionEVsEt_, hist_resolutionEVsEta_, hist_resolutionEVsPhi_, hist_resolutionPhiVsE_, hist_resolutionPhiVsEt_, hist_resolutionPhiVsEta_, and hist_resolutionPhiVsPhi_.

Referenced by endJob().

{
  hist_Et_->SetOption("e");
  hist_EtOverTruth_->SetOption("e");
  hist_EtEfficiency_->SetOption("e");
  hist_EtNumRecoOverNumTrue_->SetOption("e");
  hist_EtOverTruthVsEt_->SetOption("e");
  hist_EtOverTruthVsE_->SetOption("e");
  hist_EtOverTruthVsEta_->SetOption("e");
  hist_EtOverTruthVsPhi_->SetOption("e");
  hist_resolutionEtVsEt_->SetOption("e");
  hist_resolutionEtVsE_->SetOption("e");
  hist_resolutionEtVsEta_->SetOption("e");
  hist_resolutionEtVsPhi_->SetOption("e");

  hist_E_->SetOption("e");
  hist_EOverTruth_->SetOption("e");
  hist_EEfficiency_->SetOption("e");
  hist_ENumRecoOverNumTrue_->SetOption("e");
  hist_EOverTruthVsEt_->SetOption("e");
  hist_EOverTruthVsE_->SetOption("e");
  hist_EOverTruthVsEta_->SetOption("e");
  hist_EOverTruthVsPhi_->SetOption("e");
  hist_resolutionEVsEt_->SetOption("e");
  hist_resolutionEVsE_->SetOption("e");
  hist_resolutionEVsEta_->SetOption("e");
  hist_resolutionEVsPhi_->SetOption("e");

  hist_Eta_->SetOption("e");
  hist_EtaOverTruth_->SetOption("e");
  hist_EtaEfficiency_->SetOption("e");
  hist_EtaNumRecoOverNumTrue_->SetOption("e");
  hist_deltaEtaVsEt_->SetOption("e");
  hist_deltaEtaVsE_->SetOption("e");
  hist_deltaEtaVsEta_->SetOption("e");
  hist_deltaEtaVsPhi_->SetOption("e");
  hist_resolutionEtaVsEt_->SetOption("e");
  hist_resolutionEtaVsE_->SetOption("e");
  hist_resolutionEtaVsEta_->SetOption("e");
  hist_resolutionEtaVsPhi_->SetOption("e");

  hist_Phi_->SetOption("e");
  hist_PhiOverTruth_->SetOption("e");
  hist_PhiEfficiency_->SetOption("e");
  hist_PhiNumRecoOverNumTrue_->SetOption("e");
  hist_deltaPhiVsEt_->SetOption("e");
  hist_deltaPhiVsE_->SetOption("e");
  hist_deltaPhiVsEta_->SetOption("e");
  hist_deltaPhiVsPhi_->SetOption("e");
  hist_resolutionPhiVsEt_->SetOption("e");
  hist_resolutionPhiVsE_->SetOption("e");
  hist_resolutionPhiVsEta_->SetOption("e");
  hist_resolutionPhiVsPhi_->SetOption("e");

  hist_All_recoMass_->SetOption("e");
  hist_BarrelOnly_recoMass_->SetOption("e");
  hist_EndcapOnly_recoMass_->SetOption("e");
  hist_Mixed_recoMass_->SetOption("e");
  hist_recoMass_withBackgroud_NoEtCut_->SetOption("e");
  hist_recoMass_withBackgroud_5EtCut_->SetOption("e");
  hist_recoMass_withBackgroud_10EtCut_->SetOption("e");
  hist_recoMass_withBackgroud_20EtCut_->SetOption("e");
}

Member Data Documentation

int EgammaObjects::EtCut [private]

Definition at line 39 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), and EgammaObjects().

Definition at line 75 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

Definition at line 79 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

Definition at line 67 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

Definition at line 71 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

Definition at line 83 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

TH1D* EgammaObjects::hist_E_ [private]
TH1D* EgammaObjects::hist_Et_ [private]
TH1D* EgammaObjects::hist_Eta_ [private]

Definition at line 74 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

Definition at line 78 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

double EgammaObjects::hist_max_E_ [private]
double EgammaObjects::hist_max_Et_ [private]
double EgammaObjects::hist_max_Eta_ [private]

Definition at line 66 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

double EgammaObjects::hist_max_Phi_ [private]

Definition at line 70 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

Definition at line 82 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

Definition at line 73 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

Definition at line 77 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

double EgammaObjects::hist_min_E_ [private]
double EgammaObjects::hist_min_Et_ [private]
double EgammaObjects::hist_min_Eta_ [private]

Definition at line 65 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

double EgammaObjects::hist_min_Phi_ [private]

Definition at line 69 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

Definition at line 81 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

TH1D* EgammaObjects::hist_Phi_ [private]

Definition at line 33 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), and loadCMSSWObjects().

std::string EgammaObjects::particleString [private]

Definition at line 34 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), and loadCMSSWObjects().

TFile* EgammaObjects::rootFile_ [private]

Definition at line 31 of file EgammaObjects.h.

Referenced by EgammaObjects(), endJob(), saveHistos(), and ~EgammaObjects().