Photon selection for top analyses. More...

#include <Photon.h>

Inheritance diagram for top::Photon:

Collaboration diagram for top::Photon:

Public Member Functions
	Photon (double ptcut, double etamax, IsolationBase *isolation)
	Class to help select good photons. More...

	Photon (double ptcut, double etamax, const std::string &tightID, const std::string &looseID, IsolationBase *isolation)
	Class to help select good photons. More...

bool	passSelection (const xAOD::Photon &ph) const override
	The cuts to select good photons for your analysis should be implemented in here. More...

bool	passSelectionLoose (const xAOD::Photon &ph) const override
	The loose selection. More...

void	print (std::ostream &) const override
	Print some useful information about the photon selection. More...

Private Member Functions
bool	passSelectionNoIsolation (const xAOD::Photon &ph, const std::string &photon_selection) const
	Do all the cuts except for the isolation. More...

std::string	getConfigFile (int operatingPoint)

Private Attributes
double	m_ptcut

double	m_etamax

std::string	m_photon_selection

std::string	m_loose_photon_selection

std::unique_ptr< top::IsolationBase >	m_isolation

asg::AnaToolHandle< IAsgDeadHVCellRemovalTool >	m_deadHVTool

Detailed Description

Photon selection for top analyses.

For Release 21 specific cleaning

Definition at line 23 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/TopObjectSelectionTools/Photon.h.

Constructor & Destructor Documentation

◆ Photon() [1/2]

Photon::Photon	(	double	ptcut,
		double	etamax,
		IsolationBase *	isolation
	)

Class to help select good photons.

Parameters

ptcut	The minimum pT cut to apply to the photons.
etamax	The maximum eta cut
isolation	nullptr for un-isolated, or a new "isolation object" to apply isolation cuts

Definition at line 17 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Photon.cxx.

                                                                       :
     m_ptcut(ptcut),
     m_etamax(etamax),
     m_photon_selection("DFCommonPhotonsIsEMTight"),
     m_loose_photon_selection("DFCommonPhotonsIsEMLoose"),
     m_isolation(isolation) {
   }

◆ Photon() [2/2]

Photon::Photon	(	double	ptcut,
		double	etamax,
		const std::string &	tightID,
		const std::string &	looseID,
		IsolationBase *	isolation
	)

Class to help select good photons.

Parameters

ptcut	The minimum pT cut to apply to the photons.
etamax	The maximum eta cut
tightID	The ID level used to select tight photons.
looseID	The ID level used to select loose photons.
isolation	nullptr for un-isolated, or a new "isolation object" to apply isolation cuts

Definition at line 25 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Photon.cxx.

                                            :
     m_ptcut(ptcut),
     m_etamax(etamax),
     m_photon_selection(tightID),
     m_loose_photon_selection(looseID),
     m_isolation(isolation) {
     // Make a map of shortcuts e.g "Tight = DFCommonPhotonsIsEMTight"
     std::map<std::string, std::string> id_map;
     id_map["Tight"] = "DFCommonPhotonsIsEMTight";
     id_map["Loose"] = "DFCommonPhotonsIsEMLoose";
     id_map["None"] = "None";
     // If ID in map then set it to full name, else keep as is.
     if (id_map.find(tightID) != id_map.end()) m_photon_selection = id_map[tightID];
     if (id_map.find(looseID) != id_map.end()) m_loose_photon_selection = id_map[looseID];
  
     m_deadHVTool.setTypeAndName("AsgDeadHVCellRemovalTool/deadHVTool");
     top::check(m_deadHVTool.retrieve(), "Failed to setup Egamma DeadHVCellRemovalTool");
   }

Member Function Documentation

◆ getConfigFile()

std::string top::Photon::getConfigFile ( int operatingPoint )

private

◆ passSelection()

bool Photon::passSelection ( const xAOD::Photon & ph ) const

overridevirtual

The cuts to select good photons for your analysis should be implemented in here.

Parameters

ph	The photon to cut on (all photons in the event are passed to the tool)

Returns: True if this is a good photon, false otherwise.

Implements top::PhotonSelectionBase.

Definition at line 47 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Photon.cxx.

                                                        {
     if (!passSelectionNoIsolation(ph, m_photon_selection)) return false;
  
     if (m_isolation && !m_isolation->passSelection(ph)) return false;
  
     return true;
   }

◆ passSelectionLoose()

bool Photon::passSelectionLoose ( const xAOD::Photon & ph ) const

overridevirtual

The loose selection.

Parameters

ph

Returns

Implements top::PhotonSelectionBase.

Definition at line 55 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Photon.cxx.

                                                             {
     if (!passSelectionNoIsolation(ph, m_loose_photon_selection)) return false;
  
     if (m_isolation && !m_isolation->passSelectionLoose(ph)) return false;
  
     return true;
   }

◆ passSelectionNoIsolation()

bool Photon::passSelectionNoIsolation	(	const xAOD::Photon &	ph,
		const std::string &	photon_selection
	)		const

private

Do all the cuts except for the isolation.

Parameters

ph	The photon in question

Returns: True if the photon passes

Definition at line 63 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Photon.cxx.

                                                                                      {
     // Photon author : AuthorPhoton/AuthorAmbiguous.
     // Also, for special-case recovery of soft photons:
     // xAOD::EgammaParameters::AuthorCaloTopo35.
     if (ph.author() != xAOD::EgammaParameters::AuthorPhoton
         && ph.author() != xAOD::EgammaParameters::AuthorAmbiguous) return false;
  
     if (ph.pt() < m_ptcut) return false;
  
     // removing bad photon cluser [http://cern.ch/go/kp8F]
     if (!ph.isGoodOQ(xAOD::EgammaParameters::BADCLUSPHOTON)) return false;
  
     //  Photon selection using derivation decorations
     if (photon_selection != "None")
       if (!ph.auxdataConst<char>(photon_selection)) return false;
  
     if (!ph.caloCluster()) return false;
  
     float abs_eta = std::abs(ph.caloCluster()->etaBE(2));
     if (abs_eta > m_etamax) return false;
  
     // Remove crack region
     if (abs_eta > 1.37 && abs_eta < 1.52) return false;
  
     // Photon cleaning
     static const SG::AuxElement::ConstAccessor<char> accDFCommonPhotonsCleaning("DFCommonPhotonsCleaning");
     static const SG::AuxElement::ConstAccessor<char> accDFCommonPhotonsCleaningNoTime("DFCommonPhotonsCleaningNoTime");
  
     if (!accDFCommonPhotonsCleaning(ph)) return false;
     if (!accDFCommonPhotonsCleaningNoTime(ph)) return false;
     
     // removing photon cluster in EMEC bad HV regions
     // https://twiki.cern.ch/twiki/bin/view/AtlasProtected/EGammaIdentificationRun2#Removal_of_Electron_Photon_clust
     if (!m_deadHVTool->accept(ph)) return false;
  
  
     return true;
   }

◆ print()

void Photon::print ( std::ostream & os ) const

overridevirtual

Print some useful information about the photon selection.

Usually this goes to the log file, so you know what you ran with.

Parameters

Where the print-out should go, e.g. msg stream.

Implements top::PhotonSelectionBase.

Definition at line 103 of file PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Photon.cxx.

                                          {
     os << "Photon\n";
     os << "    * pT > " << m_ptcut << "\n";
     os << "    * |eta| < " << m_etamax << "\n";
     os << "    * Photon PID: " << m_photon_selection << "\n";
     os << "    * Loose Photon PID: " << m_loose_photon_selection << "\n";
     if (!m_isolation) os << "    * No isolation requirement\n";
     else m_isolation->print(os);
   }