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

Inheritance diagram for PhotonIDSimpleAnalyzer:

Classes
struct	struct_recPhoton

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

	PhotonIDSimpleAnalyzer (const edm::ParameterSet &)

	~PhotonIDSimpleAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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>="">

Definition at line 73 of file PhotonIDSimpleAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 160 of file PhotonIDSimpleAnalyzer.cc.

References createPhotonTTree_, edm::ParameterSet::getParameter(), maxPhotonAbsEta_, maxPhotonHoverE_, minPhotonAbsEta_, minPhotonEt_, minPhotonR9_, outputFile_, rootFile_, 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 ( )

override

Definition at line 184 of file PhotonIDSimpleAnalyzer.cc.

References rootFile_.

                                                 {
   // 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
	)

override

Definition at line 245 of file PhotonIDSimpleAnalyzer.cc.

                                                                                  {
   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 )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 194 of file PhotonIDSimpleAnalyzer.cc.

References createPhotonTTree_, h_ebeeGap_, h_ebetagap_, h_ebphigap_, h_eedeeGap_, h_eeringGap_, h_hadoverem_, h_isoEcalRecHit_, h_isoHcalRecHit_, h_nPassEM_, h_nPassingPho_, h_nPho_, h_ntrk_hollow_, h_ntrk_solid_, h_photonEt_, h_photonEta_, h_photonInAnyGap_, h_photonPhi_, h_photonScEt_, h_photonScEta_, h_photonScEtaWidth_, h_photonScPhi_, h_r9_, h_trk_pt_hollow_, h_trk_pt_solid_, PhotonIDSimpleAnalyzer::struct_recPhoton::isolationEcalRecHit, Pi, recPhoton, rootFile_, and tree_PhotonAll_.

                                       {
   // 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 )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 367 of file PhotonIDSimpleAnalyzer.cc.

References DEFINE_FWK_MODULE, h_ebeeGap_, h_ebetagap_, h_ebphigap_, h_eedeeGap_, h_eeringGap_, h_hadoverem_, h_isoEcalRecHit_, h_isoHcalRecHit_, h_nPassEM_, h_nPassingPho_, h_nPho_, h_ntrk_hollow_, h_ntrk_solid_, h_photonEt_, h_photonEta_, h_photonInAnyGap_, h_photonPhi_, h_photonScEt_, h_photonScEta_, h_photonScEtaWidth_, h_photonScPhi_, h_r9_, h_trk_pt_hollow_, h_trk_pt_solid_, and rootFile_.

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