PhotonIDSimpleAnalyzer Class Reference

#include <RecoEgamma/PhotonIdentification/test/PhotonIDSimpleAnalyzer.cc>

Inheritance diagram for PhotonIDSimpleAnalyzer:

Classes
struct	struct_recPhoton

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()
	PhotonIDSimpleAnalyzer (const edm::ParameterSet &)
	~PhotonIDSimpleAnalyzer ()
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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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 Attributes
bool	createPhotonTTree_
TH1F *	h_ebeeGap_
TH1F *	h_ebetagap_
TH1F *	h_ebphigap_
TH1F *	h_eedeeGap_
TH1F *	h_eeringGap_
TH1F *	h_hadoverem_
TH1F *	h_isoEcalRecHit_
TH1F *	h_isoHcalRecHit_
TH1F *	h_nPassEM_
TH1F *	h_nPassingPho_
TH1F *	h_nPho_
TH1F *	h_ntrk_hollow_
TH1F *	h_ntrk_solid_
TH1F *	h_photonEt_
TH1F *	h_photonEta_
TH1F *	h_photonInAnyGap_
TH1F *	h_photonPhi_
TH1F *	h_photonScEt_
TH1F *	h_photonScEta_
TH1F *	h_photonScEtaWidth_
TH1F *	h_photonScPhi_
TH1F *	h_r9_
TH1F *	h_trk_pt_hollow_
TH1F *	h_trk_pt_solid_
double	maxPhotonAbsEta_
double	maxPhotonHoverE_
double	minPhotonAbsEta_
double	minPhotonEt_
double	minPhotonR9_
std::string	outputFile_
struct_recPhoton	recPhoton
TFile *	rootFile_
TTree *	tree_PhotonAll_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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: Generate various histograms for cuts and important photon ID parameters using a data sample of photons in QCD events.

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

Description: Analyzer to make a load of histograms for the improvement of the PhotonID object

Implementation: \

Author

: J. Stilley, A. Askew May 2008

Definition at line 37 of file PhotonIDSimpleAnalyzer.h.

Constructor & Destructor Documentation

PhotonIDSimpleAnalyzer::PhotonIDSimpleAnalyzer ( const edm::ParameterSet &  ps )

explicit

Definition at line 64 of file PhotonIDSimpleAnalyzer.cc.

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

{
  // Read Parameters from configuration file

  // output filename
  outputFile_   = ps.getParameter<std::string>("outputFile");
  // Read variables that must be passed to allow a 
  //  supercluster to be placed in histograms as a photon.
  minPhotonEt_     = ps.getParameter<double>("minPhotonEt");
  minPhotonAbsEta_ = ps.getParameter<double>("minPhotonAbsEta");
  maxPhotonAbsEta_ = ps.getParameter<double>("maxPhotonAbsEta");
  minPhotonR9_     = ps.getParameter<double>("minPhotonR9");
  maxPhotonHoverE_ = ps.getParameter<double>("maxPhotonHoverE");

  // Read variable to that decidedes whether
  // a TTree of photons is created or not
  createPhotonTTree_ = ps.getParameter<bool>("createPhotonTTree");

  // open output file to store histograms
  rootFile_ = TFile::Open(outputFile_.c_str(),"RECREATE");
}

PhotonIDSimpleAnalyzer::~PhotonIDSimpleAnalyzer

(

)

Definition at line 89 of file PhotonIDSimpleAnalyzer.cc.

{

// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)

  delete rootFile_;

}

Member Function Documentation

void PhotonIDSimpleAnalyzer::analyze ( const edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Implements edm::EDAnalyzer.

Definition at line 153 of file PhotonIDSimpleAnalyzer.cc.

References HcalObjRepresent::Fill(), edm::Event::getByLabel(), and interactiveExample::photons.

{
  
  using namespace std;
  using namespace edm;
  
  // grab photons
  Handle<reco::PhotonCollection> photonColl;
  evt.getByLabel("photons", "", photonColl);

  Handle<edm::ValueMap<Bool_t> > loosePhotonQual;
  evt.getByLabel("PhotonIDProd", "PhotonCutBasedIDLoose", loosePhotonQual);
  Handle<edm::ValueMap<Bool_t> > looseEMQual;
  evt.getByLabel("PhotonIDProd","PhotonCutBasedIDLooseEM",looseEMQual);
  // grab PhotonId objects  
//   Handle<reco::PhotonIDAssociationCollection> photonIDMapColl;
//   evt.getByLabel("PhotonIDProd", "PhotonAssociatedID", photonIDMapColl);
  
  // create reference to the object types we are interested in
  const reco::PhotonCollection *photons = photonColl.product();  
  const edm::ValueMap<Bool_t> *phoMap = loosePhotonQual.product();
  const edm::ValueMap<Bool_t> *lEMMap = looseEMQual.product();
  int photonCounter = 0;
  int idxpho=0;
  reco::PhotonCollection::const_iterator pho;
  for (pho = (*photons).begin(); pho!= (*photons).end(); pho++){   
    
    edm::Ref<reco::PhotonCollection> photonref(photonColl, idxpho);
    //reco::PhotonIDAssociationCollection::const_iterator photonIter = phoMap->find(photonref);
    //const reco::PhotonIDRef &phtn = photonIter->val;
    //const reco::PhotonRef &pho = photonIter->key;

    float photonEt       = pho->et();
    float superClusterEt = (pho->superCluster()->energy())/(cosh(pho->superCluster()->position().eta()));
    Bool_t LoosePhotonQu = (*phoMap)[photonref];
    h_nPassingPho_->Fill(LoosePhotonQu);
    Bool_t LooseEMQu = (*lEMMap)[photonref];
    h_nPassEM_->Fill(LooseEMQu);
    // Only store photon candidates (SuperClusters) that pass some simple cuts
    bool passCuts = (              photonEt > minPhotonEt_     ) &&
                    (      fabs(pho->eta()) > minPhotonAbsEta_ ) &&
                    (      fabs(pho->eta()) < maxPhotonAbsEta_ ) &&
                    (             pho->r9() > minPhotonR9_     ) &&
                    ( pho->hadronicOverEm() < maxPhotonHoverE_ ) ;

    if ( passCuts )
    {
      //                fill histograms                    //
      // PhotonID Variables
      h_isoEcalRecHit_->Fill(pho->ecalRecHitSumEtConeDR04());
      h_isoHcalRecHit_->Fill(pho->hcalTowerSumEtConeDR04());
      h_trk_pt_solid_ ->Fill(pho->trkSumPtSolidConeDR04());
      h_trk_pt_hollow_->Fill(pho->trkSumPtHollowConeDR04());
      h_ntrk_solid_->   Fill(pho->nTrkSolidConeDR04());
      h_ntrk_hollow_->  Fill(pho->nTrkHollowConeDR04());
      h_ebetagap_->        Fill(pho->isEBEtaGap());
      h_ebphigap_->        Fill(pho->isEBPhiGap());
      h_eeringGap_->        Fill(pho->isEERingGap()); 
      h_eedeeGap_->        Fill(pho->isEEDeeGap()); 
      h_ebeeGap_->      Fill(pho->isEBEEGap());
      h_r9_->           Fill(pho->r9());

      // Photon Variables
      h_photonEt_->  Fill(photonEt);
      h_photonEta_-> Fill(pho->eta());
      h_photonPhi_-> Fill(pho->phi());
      h_hadoverem_-> Fill(pho->hadronicOverEm());

      // Photon's SuperCluster Variables
      // eta is with respect to detector (not physics) vertex,
      // thus Et and eta are different from photon.
      h_photonScEt_->      Fill(superClusterEt);
      h_photonScEta_->     Fill(pho->superCluster()->position().eta());
      h_photonScPhi_->     Fill(pho->superCluster()->position().phi());
      h_photonScEtaWidth_->Fill(pho->superCluster()->etaWidth());
      
      // It passed photon cuts, mark it

      //                fill TTree (optional)              //
      if ( createPhotonTTree_ ) {
    recPhoton.isolationEcalRecHit    = pho->ecalRecHitSumEtConeDR04();
    recPhoton.isolationHcalRecHit    = pho->hcalTowerSumEtConeDR04();
    recPhoton.isolationSolidTrkCone  = pho->trkSumPtSolidConeDR04();
    recPhoton.isolationHollowTrkCone = pho->trkSumPtHollowConeDR04();
    recPhoton.nTrkSolidCone          = pho->nTrkSolidConeDR04();
    recPhoton.nTrkHollowCone         = pho->nTrkHollowConeDR04();
    recPhoton.isEBEtaGap                = pho->isEBEtaGap();
    recPhoton.isEBPhiGap                = pho->isEBPhiGap();
    recPhoton.isEERingGap                = pho->isEERingGap();
    recPhoton.isEEDeeGap                = pho->isEEDeeGap();
    recPhoton.isEBEEGap              = pho->isEBEEGap();
        recPhoton.r9                     = pho->r9();
        recPhoton.et                     = pho->et();
        recPhoton.eta                    = pho->eta();
        recPhoton.phi                    = pho->phi();
        recPhoton.hadronicOverEm         = pho->hadronicOverEm();

        // Fill the tree (this records all the recPhoton.* since
        // tree_PhotonAll_ was set to point at that.
        tree_PhotonAll_->Fill();
      }

      // Record whether it was near any module gap.
      // Very convoluted at the moment.
      bool inAnyGap = pho->isEBEEGap() || (pho->isEB()&&pho->isEBEtaGap()) ||(pho->isEB()&&pho->isEBPhiGap())  || (pho->isEE()&&pho->isEERingGap()) || (pho->isEE()&&pho->isEEDeeGap()) ;
      if (inAnyGap) {
        h_photonInAnyGap_->Fill(1.0);
      } else {
        h_photonInAnyGap_->Fill(0.0);
      }

      photonCounter++;
    } 
    else
    {
      // This didn't pass photon cuts, mark it
   
    }
    idxpho++;
  } // End Loop over photons
  h_nPho_->Fill(photonCounter);
 
}

void PhotonIDSimpleAnalyzer::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 103 of file PhotonIDSimpleAnalyzer.cc.

References Pi.

{

  // go to *OUR* rootfile
  rootFile_->cd();

  // Book Histograms
  // PhotonID Histograms
  h_isoEcalRecHit_ = new TH1F("photonEcalIso",          "Ecal Rec Hit Isolation", 300, 0, 300);
  h_isoHcalRecHit_ = new TH1F("photonHcalIso",          "Hcal Rec Hit Isolation", 300, 0, 300);
  h_trk_pt_solid_  = new TH1F("photonTrackSolidIso",    "Sum of track pT in a cone of #DeltaR" , 300, 0, 300);
  h_trk_pt_hollow_ = new TH1F("photonTrackHollowIso",   "Sum of track pT in a hollow cone" ,     300, 0, 300);
  h_ntrk_solid_    = new TH1F("photonTrackCountSolid",  "Number of tracks in a cone of #DeltaR", 100, 0, 100);
  h_ntrk_hollow_   = new TH1F("photonTrackCountHollow", "Number of tracks in a hollow cone",     100, 0, 100);
  h_ebetagap_         = new TH1F("photonInEBEtagap",          "Ecal Barrel eta gap flag",  2, -0.5, 1.5);
  h_ebphigap_         = new TH1F("photonInEBEtagap",          "Ecal Barrel phi gap flag",  2, -0.5, 1.5);
  h_eeringGap_         = new TH1F("photonInEERinggap",          "Ecal Endcap ring gap flag",  2, -0.5, 1.5);
  h_eedeeGap_         = new TH1F("photonInEEDeegap",          "Ecal Endcap dee gap flag",  2, -0.5, 1.5);
  h_ebeeGap_       = new TH1F("photonInEEgap",          "Ecal Barrel/Endcap gap flag",  2, -0.5, 1.5);
  h_r9_            = new TH1F("photonR9",               "R9 = E(3x3) / E(SuperCluster)", 300, 0, 3);

  // Photon Histograms
  h_photonEt_      = new TH1F("photonEt",     "Photon E_{T}",  200,  0, 200);
  h_photonEta_     = new TH1F("photonEta",    "Photon #eta",   800, -4,   4);
  h_photonPhi_     = new TH1F("photonPhi",    "Photon #phi",   628, -1.*TMath::Pi(), TMath::Pi());
  h_hadoverem_     = new TH1F("photonHoverE", "Hadronic over EM", 200, 0, 1);

  // Photon's SuperCluster Histograms
  h_photonScEt_       = new TH1F("photonScEt",  "Photon SuperCluster E_{T}", 200,  0, 200);
  h_photonScEta_      = new TH1F("photonScEta", "Photon #eta",               800, -4,   4);
  h_photonScPhi_      = new TH1F("photonScPhi", "Photon #phi",628, -1.*TMath::Pi(), TMath::Pi());
  h_photonScEtaWidth_ = new TH1F("photonScEtaWidth","#eta-width",            100,  0,  .1);

  // Composite or Other Histograms
  h_photonInAnyGap_   = new TH1F("photonInAnyGap",     "Photon in any gap flag",  2, -0.5, 1.5);
  h_nPassingPho_      = new TH1F("photonLoosePhoton", "Total number photons (0=NotPassing, 1=Passing)", 2, -0.5, 1.5);
  h_nPassEM_          = new TH1F("photonLooseEM", "Total number photons (0=NotPassing, 1=Passing)", 2, -0.5, 1.5);
  h_nPho_             = new TH1F("photonCount",        "Number of photons passing cuts in event",  10,  0,  10);

  // Create a TTree of photons if set to 'True' in config file
  if ( createPhotonTTree_ ) {
    tree_PhotonAll_     = new TTree("TreePhotonAll", "Reconstructed Photon");
    tree_PhotonAll_->Branch("recPhoton", &recPhoton.isolationEcalRecHit, "isolationEcalRecHit/F:isolationHcalRecHit:isolationSolidTrkCone:isolationHollowTrkCone:nTrkSolidCone:nTrkHollowCone:isEBGap:isEEGap:isEBEEGap:r9:et:eta:phi:hadronicOverEm");
  }
}

void PhotonIDSimpleAnalyzer::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 285 of file PhotonIDSimpleAnalyzer.cc.

{

  // go to *OUR* root file and store histograms
  rootFile_->cd();

  // PhotonID Histograms
  h_isoEcalRecHit_->Write();
  h_isoHcalRecHit_->Write();
  h_trk_pt_solid_-> Write();
  h_trk_pt_hollow_->Write();
  h_ntrk_solid_->   Write();
  h_ntrk_hollow_->  Write();
  h_ebetagap_->     Write();
  h_ebphigap_->     Write();
  h_eeringGap_->     Write();
  h_eedeeGap_->     Write();
  h_ebeeGap_->   Write();
  h_r9_->        Write();

  // Photon Histograms
  h_photonEt_->  Write();
  h_photonEta_-> Write();
  h_photonPhi_-> Write();
  h_hadoverem_-> Write();

  // Photon's SuperCluster Histograms
  h_photonScEt_->      Write();
  h_photonScEta_->     Write();
  h_photonScPhi_->     Write();
  h_photonScEtaWidth_->Write();

  // Composite or Other Histograms
  h_photonInAnyGap_->Write();
  h_nPassingPho_->   Write();
  h_nPassEM_->       Write();
  h_nPho_->          Write();

  // Write the root file (really writes the TTree)
  rootFile_->Write();
  rootFile_->Close();

}

Member Data Documentation

bool PhotonIDSimpleAnalyzer::createPhotonTTree_

private

Definition at line 54 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_ebeeGap_

private

Definition at line 95 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_ebetagap_

private

Definition at line 91 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_ebphigap_

private

Definition at line 92 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_eedeeGap_

private

Definition at line 94 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_eeringGap_

private

Definition at line 93 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_hadoverem_

private

Definition at line 102 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_isoEcalRecHit_

private

Definition at line 85 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_isoHcalRecHit_

private

Definition at line 86 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_nPassEM_

private

Definition at line 113 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_nPassingPho_

private

Definition at line 112 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_nPho_

private

Definition at line 114 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_ntrk_hollow_

private

Definition at line 90 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_ntrk_solid_

private

Definition at line 89 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonEt_

private

Definition at line 99 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonEta_

private

Definition at line 100 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonInAnyGap_

private

Definition at line 111 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonPhi_

private

Definition at line 101 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonScEt_

private

Definition at line 105 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonScEta_

private

Definition at line 106 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonScEtaWidth_

private

Definition at line 108 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_photonScPhi_

private

Definition at line 107 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_r9_

private

Definition at line 96 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_trk_pt_hollow_

private

Definition at line 88 of file PhotonIDSimpleAnalyzer.h.

TH1F* PhotonIDSimpleAnalyzer::h_trk_pt_solid_

private

Definition at line 87 of file PhotonIDSimpleAnalyzer.h.

double PhotonIDSimpleAnalyzer::maxPhotonAbsEta_

private

Definition at line 51 of file PhotonIDSimpleAnalyzer.h.

double PhotonIDSimpleAnalyzer::maxPhotonHoverE_

private

Definition at line 53 of file PhotonIDSimpleAnalyzer.h.

double PhotonIDSimpleAnalyzer::minPhotonAbsEta_

private

Definition at line 50 of file PhotonIDSimpleAnalyzer.h.

double PhotonIDSimpleAnalyzer::minPhotonEt_

private

Definition at line 49 of file PhotonIDSimpleAnalyzer.h.

double PhotonIDSimpleAnalyzer::minPhotonR9_

private

Definition at line 52 of file PhotonIDSimpleAnalyzer.h.

std::string PhotonIDSimpleAnalyzer::outputFile_

private

Definition at line 48 of file PhotonIDSimpleAnalyzer.h.

struct_recPhoton PhotonIDSimpleAnalyzer::recPhoton

private

Definition at line 77 of file PhotonIDSimpleAnalyzer.h.

TFile* PhotonIDSimpleAnalyzer::rootFile_

private

Definition at line 80 of file PhotonIDSimpleAnalyzer.h.

TTree* PhotonIDSimpleAnalyzer::tree_PhotonAll_

private

Definition at line 117 of file PhotonIDSimpleAnalyzer.h.