ATLAS Offline Software
Public Member Functions | Private Member Functions | Private Attributes | List of all members
TrigEgammaPrecisionPhotonHypoTool Class Reference

#include <TrigEgammaPrecisionPhotonHypoTool.h>

Inheritance diagram for TrigEgammaPrecisionPhotonHypoTool:
Collaboration diagram for TrigEgammaPrecisionPhotonHypoTool:

Public Member Functions

 TrigEgammaPrecisionPhotonHypoTool (const std::string &type, const std::string &name, const IInterface *parent)
 
virtual StatusCode initialize () override
 
virtual StatusCode decide (std::vector< ITrigEgammaPrecisionPhotonHypoTool::PhotonInfo > &input) const override
 
virtual bool decide (const ITrigEgammaPrecisionPhotonHypoTool::PhotonInfo &i) const override
 

Private Member Functions

int findCutIndex (float eta) const
 

Private Attributes

HLT::Identifier m_decisionId
 
Gaudi::Property< std::vector< float > > m_etabin { this, "EtaBins", {} , "Bins of eta" }
 selection variable for PRECISION calo selection:eta bins More...
 
Gaudi::Property< std::vector< float > > m_eTthr { this, "ETthr", {}, "ET Threshold" }
 
Gaudi::Property< float > m_detacluster { this, "dETACLUSTERthr", 0. , "" }
 
Gaudi::Property< float > m_dphicluster { this, "dPHICLUSTERthr", 0. , "" }
 
Gaudi::Property< bool > m_doNoPid { this, "DoNoPid", false , "No Pid/Isolation applied" }
 
Gaudi::Property< std::string > m_pidName {this, "PidName", "", "Pid name"}
 
ToolHandle< GenericMonitoringTool > m_monTool { this, "MonTool", "", "Monitoring tool" }
 
SG::ReadDecorHandleKey< xAOD::EventInfo > m_avgMuKey { this, "averageInteractionsPerCrossingKey", "EventInfo.averageInteractionsPerCrossing", "Decoration for Average Interaction Per Crossing" }
 

Detailed Description

Definition at line 19 of file TrigEgammaPrecisionPhotonHypoTool.h.

Constructor & Destructor Documentation

◆ TrigEgammaPrecisionPhotonHypoTool()

TrigEgammaPrecisionPhotonHypoTool::TrigEgammaPrecisionPhotonHypoTool ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Definition at line 17 of file TrigEgammaPrecisionPhotonHypoTool.cxx.

20  : base_class( type, name, parent ),
21  m_decisionId( HLT::Identifier::fromToolName( name ) ) {
22 }

Member Function Documentation

◆ decide() [1/2]

bool TrigEgammaPrecisionPhotonHypoTool::decide ( const ITrigEgammaPrecisionPhotonHypoTool::PhotonInfo &  i) const
overridevirtual

Definition at line 63 of file TrigEgammaPrecisionPhotonHypoTool.cxx.

64 {
65 
66  bool pass = false;
67 
68  auto mon_ET = Monitored::Scalar( "Et_em", -1.0 );
69  auto mon_dEta = Monitored::Scalar( "dEta", -1.0 );
70  auto mon_dPhi = Monitored::Scalar( "dPhi", -1.0 );
71  auto mon_etaBin = Monitored::Scalar( "EtaBin", -1.0 );
72  auto mon_Eta = Monitored::Scalar( "Eta", -99. );
73  auto mon_Phi = Monitored::Scalar( "Phi", -99. );
74  auto mon_mu = Monitored::Scalar("mu", -1.);
75 
76  auto PassedCuts = Monitored::Scalar<int>( "CutCounter", -1 );
77  auto monitorIt = Monitored::Group( m_monTool, mon_ET, mon_dEta, mon_dPhi,
78  mon_etaBin, mon_Eta, mon_Phi, mon_mu,
79  PassedCuts );
80 
81  // when leaving scope it will ship data to monTool
82  PassedCuts = PassedCuts + 1; //got called (data in place)
83 
84  float ET(0), dEta(0), dPhi(0), eta(0), phi(0);
85 
86  auto roiDescriptor = input.roi;
87 
88  if ( fabs( roiDescriptor->eta() ) > 2.6 ) {
89  ATH_MSG_DEBUG( "REJECT The photon had eta coordinates beyond the EM fiducial volume : "
90  << roiDescriptor->eta() << "; stop the chain now" );
91  pass=false; // special case
92  return pass;
93  }
94 
95  ATH_MSG_DEBUG( "; RoI ID = " << roiDescriptor->roiId()
96  << ": Eta = " << roiDescriptor->eta()
97  << ", Phi = " << roiDescriptor->phi() );
98 
99  // fill local variables for RoI reference position
100  double etaRef = roiDescriptor->eta();
101  double phiRef = roiDescriptor->phi();
102  // correct phi the to right range ( probably not needed anymore )
103  if ( fabs( phiRef ) > M_PI ) phiRef -= 2*M_PI; // correct phi if outside range
104 
105  auto pClus = input.photon->caloCluster();
106 
107  float absEta = fabs( pClus->eta() );
108  const int cutIndex = findCutIndex( absEta );
109 
110  dEta = pClus->eta() - etaRef;
111  // Deal with angle diferences greater than Pi
112  dPhi = fabs( pClus->phi() - phiRef );
113  dPhi = ( dPhi < M_PI ? dPhi : 2*M_PI - dPhi ); // TB why only <
114  ET = pClus->et();
115  eta = pClus->eta();
116  phi = pClus->phi();
117  // apply cuts: DeltaEta( clus-ROI )
118  ATH_MSG_DEBUG( "Photon : eta=" << pClus->eta()
119  << " roi eta=" << etaRef << " DeltaEta=" << dEta
120  << " cut: <" << m_detacluster );
121 
122  if ( fabs( pClus->eta() - etaRef ) > m_detacluster ) {
123  ATH_MSG_DEBUG("REJECT Photon a cut failed");
124  return pass;
125  }
126  mon_Eta = eta;
127  mon_dEta = dEta;
128  PassedCuts = PassedCuts + 1; //Deta
129 
130  // DeltaPhi( clus-ROI )
131  ATH_MSG_DEBUG( ": phi=" << pClus->phi()
132  << " roi phi="<< phiRef << " DeltaPhi="<< dPhi
133  << " cut: <" << m_dphicluster );
134 
135  if( dPhi > m_dphicluster ) {
136  ATH_MSG_DEBUG("REJECT Clsuter dPhi cut failed");
137  return pass;
138  }
139  mon_Phi = phi;
140  mon_dPhi = dPhi;
141  PassedCuts = PassedCuts + 1; //DPhi
142 
143  // eta range
144  if ( cutIndex == -1 ) { // VD
145  ATH_MSG_DEBUG( "Photon : " << absEta << " outside eta range " << m_etabin[m_etabin.size()-1] );
146  return pass;
147  } else {
148  ATH_MSG_DEBUG( "eta bin used for cuts " << cutIndex );
149  }
150  mon_etaBin = m_etabin[cutIndex];
151  PassedCuts = PassedCuts + 1; // passed eta cut
152 
153  // ET_em
154  ATH_MSG_DEBUG( "Photon: ET_em=" << ET << " cut: >" << m_eTthr[cutIndex] );
155  if ( ET < m_eTthr[cutIndex] ) {
156  ATH_MSG_DEBUG("REJECT et cut failed");
157  return pass;
158  }
159  mon_ET = ET;
160  PassedCuts = PassedCuts + 1; // ET_em
161  if(m_doNoPid){
162  pass = true;
163  return pass;
164  }
165 
166  if(input.pidDecorator.count(m_pidName)){
167  pass = input.pidDecorator.at(m_pidName);
168  }
169 
170  // get average luminosity information to calculate LH
171  float avg_mu = 0;
172  SG::ReadDecorHandle<xAOD::EventInfo,float> eventInfoDecor(m_avgMuKey);
173  if(eventInfoDecor.isPresent()) {
174  avg_mu = eventInfoDecor(0);
175  ATH_MSG_DEBUG("Average mu " << avg_mu);
176  }
177  mon_mu = avg_mu;
178 
179  float Rhad1(0), Rhad(0), Reta(0), Rphi(0), e277(0), weta2c(0), //emax2(0),
180  Eratio(0), DeltaE(0), f1(0), weta1c(0), wtot(0), fracm(0);
181 
182 
183  // variables based on HCAL
184  // transverse energy in 1st scintillator of hadronic calorimeter/ET
185  input.photon->showerShapeValue(Rhad1, xAOD::EgammaParameters::Rhad1);
186  // transverse energy in hadronic calorimeter/ET
187  input.photon->showerShapeValue(Rhad, xAOD::EgammaParameters::Rhad);
188 
189  // variables based on S2 of EM CAL
190  // E(7*7) in 2nd sampling
191  input.photon->showerShapeValue(e277, xAOD::EgammaParameters::e277);
192  // E(3*7)/E(7*7) in 2nd sampling
193  input.photon->showerShapeValue(Reta, xAOD::EgammaParameters::Reta);
194  // E(3*3)/E(3*7) in 2nd sampling
195  input.photon->showerShapeValue(Rphi, xAOD::EgammaParameters::Rphi);
196  // shower width in 2nd sampling
197  input.photon->showerShapeValue(weta2c, xAOD::EgammaParameters::weta2);
198 
199  // variables based on S1 of EM CAL
200  // fraction of energy reconstructed in the 1st sampling
201  input.photon->showerShapeValue(f1, xAOD::EgammaParameters::f1);
202  // shower width in 3 strips in 1st sampling
203  input.photon->showerShapeValue(weta1c, xAOD::EgammaParameters::weta1);
204  // E of 2nd max between max and min in strips [NOT USED]
205  // eg->showerShapeValue(emax2, xAOD::EgammaParameters::e2tsts1);
206  // (E of 1st max in strips-E of 2nd max)/(E of 1st max+E of 2nd max)
207  input.photon->showerShapeValue(Eratio, xAOD::EgammaParameters::Eratio);
208  // E(2nd max)-E(min) in strips
209  input.photon->showerShapeValue(DeltaE, xAOD::EgammaParameters::DeltaE);
210  // total shower width in 1st sampling
211  input.photon->showerShapeValue(wtot, xAOD::EgammaParameters::wtots1);
212  // E(+/-3)-E(+/-1)/E(+/-1)
213  input.photon->showerShapeValue(fracm, xAOD::EgammaParameters::fracs1);
214 
215  ATH_MSG_DEBUG( " Rhad = " << Rhad ) ;
216  ATH_MSG_DEBUG( " Rhad1 = " << Rhad1 ) ;
217  ATH_MSG_DEBUG( " e277 = " << e277 ) ;
218  ATH_MSG_DEBUG( " Reta = " << Reta ) ;
219  ATH_MSG_DEBUG( " Rphi = " << Rphi ) ;
220  ATH_MSG_DEBUG( " weta2c = " << weta2c ) ;
221  ATH_MSG_DEBUG( " f1 = " << f1 ) ;
222  ATH_MSG_DEBUG( " weta1c = " << weta1c ) ;
223  ATH_MSG_DEBUG( " Eratio = " << Eratio ) ;
224  ATH_MSG_DEBUG( " DeltaE = " << DeltaE ) ;
225  ATH_MSG_DEBUG( " wtot = " << wtot ) ;
226  ATH_MSG_DEBUG( " fracm = " << fracm ) ;
227 
228  if ( !pass ){
229  ATH_MSG_DEBUG("REJECT isEM failed");
230  return pass;
231  } else {
232  ATH_MSG_DEBUG("ACCEPT isEM passed");
233  }
234  // Monitor showershapes
235  ATH_MSG_DEBUG( "pass = " << pass );
236 
237  return pass;
238 
239 }

◆ decide() [2/2]

virtual StatusCode TrigEgammaPrecisionPhotonHypoTool::decide ( std::vector< ITrigEgammaPrecisionPhotonHypoTool::PhotonInfo > &  input) const
overridevirtual

◆ findCutIndex()

int TrigEgammaPrecisionPhotonHypoTool::findCutIndex ( float  eta) const
private

Definition at line 241 of file TrigEgammaPrecisionPhotonHypoTool.cxx.

241  {
242  const float absEta = std::abs(eta);
243 
244  auto binIterator = std::adjacent_find( m_etabin.begin(), m_etabin.end(), [=](float left, float right){ return left < absEta and absEta < right; } );
245  if ( binIterator == m_etabin.end() ) {
246  return -1;
247  }
248  return binIterator - m_etabin.begin();
249 }

◆ initialize()

StatusCode TrigEgammaPrecisionPhotonHypoTool::initialize ( )
overridevirtual

Definition at line 25 of file TrigEgammaPrecisionPhotonHypoTool.cxx.

26 {
27  ATH_MSG_DEBUG( "Initialization completed successfully" );
28  ATH_MSG_DEBUG( "EtaBins = " << m_etabin );
29  ATH_MSG_DEBUG( "ETthr = " << m_eTthr );
30  ATH_MSG_DEBUG( "dPHICLUSTERthr = " << m_dphicluster );
31  ATH_MSG_DEBUG( "dETACLUSTERthr = " << m_detacluster );
32  ATH_MSG_DEBUG( "DoNoPid = " << m_doNoPid );
33 
34  if ( m_etabin.empty() ) {
35  ATH_MSG_ERROR( " There are no cuts set (EtaBins property is an empty list)" );
36  return StatusCode::FAILURE;
37  }
38 
39  // Now we try to retrieve the ElectronPhotonSelectorTools that we will use to apply the photon Identification. This is a *must*
40 
41 
42  // Retrieving Luminosity info
43  ATH_MSG_DEBUG( "Retrieving luminosityCondData..." );
44  ATH_CHECK( m_avgMuKey.initialize() );
45 
46  unsigned int nEtaBin = m_etabin.size();
47 
48 #define CHECK_SIZE( __n) if ( m_##__n.size() != (nEtaBin - 1) ) \
49  { ATH_MSG_DEBUG(" __n size is " << m_##__n.size() << " but needs to be " << (nEtaBin - 1) ); return StatusCode::FAILURE; }
50 
51  CHECK_SIZE( eTthr );
52 #undef CHECK_SIZE
53 
54  ATH_MSG_DEBUG( "Tool configured for chain/id: " << m_decisionId );
55 
56  if ( not m_monTool.name().empty() )
57  CHECK( m_monTool.retrieve() );
58 
59  return StatusCode::SUCCESS;
60 }

Member Data Documentation

◆ m_avgMuKey

SG::ReadDecorHandleKey<xAOD::EventInfo> TrigEgammaPrecisionPhotonHypoTool::m_avgMuKey { this, "averageInteractionsPerCrossingKey", "EventInfo.averageInteractionsPerCrossing", "Decoration for Average Interaction Per Crossing" }
private

Definition at line 47 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_decisionId

HLT::Identifier TrigEgammaPrecisionPhotonHypoTool::m_decisionId
private

Definition at line 30 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_detacluster

Gaudi::Property< float > TrigEgammaPrecisionPhotonHypoTool::m_detacluster { this, "dETACLUSTERthr", 0. , "" }
private

Definition at line 35 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_doNoPid

Gaudi::Property< bool > TrigEgammaPrecisionPhotonHypoTool::m_doNoPid { this, "DoNoPid", false , "No Pid/Isolation applied" }
private

Definition at line 37 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_dphicluster

Gaudi::Property< float > TrigEgammaPrecisionPhotonHypoTool::m_dphicluster { this, "dPHICLUSTERthr", 0. , "" }
private

Definition at line 36 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_etabin

Gaudi::Property< std::vector<float> > TrigEgammaPrecisionPhotonHypoTool::m_etabin { this, "EtaBins", {} , "Bins of eta" }
private

selection variable for PRECISION calo selection:eta bins

Definition at line 33 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_eTthr

Gaudi::Property< std::vector<float> > TrigEgammaPrecisionPhotonHypoTool::m_eTthr { this, "ETthr", {}, "ET Threshold" }
private

Definition at line 34 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_monTool

ToolHandle< GenericMonitoringTool > TrigEgammaPrecisionPhotonHypoTool::m_monTool { this, "MonTool", "", "Monitoring tool" }
private

Definition at line 44 of file TrigEgammaPrecisionPhotonHypoTool.h.

◆ m_pidName

Gaudi::Property< std::string > TrigEgammaPrecisionPhotonHypoTool::m_pidName {this, "PidName", "", "Pid name"}
private

Definition at line 41 of file TrigEgammaPrecisionPhotonHypoTool.h.

The documentation for this class was generated from the following files:
TrigEgammaPrecisionPhotonHypoTool::m_detacluster
Gaudi::Property< float > m_detacluster
Definition: TrigEgammaPrecisionPhotonHypoTool.h:35
TrigDefs::Group
Group
Properties of a chain group.
Definition: GroupProperties.h:13
CHECK_SIZE
#define CHECK_SIZE(__n)
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:67
xAOD::EgammaParameters::Reta
@ Reta
e237/e277
Definition: EgammaEnums.h:154
TrigEgammaPrecisionPhotonHypoTool::m_doNoPid
Gaudi::Property< bool > m_doNoPid
Definition: TrigEgammaPrecisionPhotonHypoTool.h:37
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:83
TrigEgammaPrecisionPhotonHypoTool::m_monTool
ToolHandle< GenericMonitoringTool > m_monTool
Definition: TrigEgammaPrecisionPhotonHypoTool.h:44
TrigEgammaPrecisionPhotonHypoTool::m_eTthr
Gaudi::Property< std::vector< float > > m_eTthr
Definition: TrigEgammaPrecisionPhotonHypoTool.h:34
M_PI
#define M_PI
Definition: ActiveFraction.h:11
xAOD::EgammaParameters::Rphi
@ Rphi
e233/e237
Definition: EgammaEnums.h:156
TrigEgammaPrecisionPhotonHypoTool::m_decisionId
HLT::Identifier m_decisionId
Definition: TrigEgammaPrecisionPhotonHypoTool.h:30
TauClusterVars::dPhi
bool dPhi(const xAOD::TauJet &tau, const xAOD::CaloVertexedTopoCluster &cluster, float &out)
Definition: ConstituentLoaderTauCluster.cxx:119
TrigEgammaPrecisionPhotonHypoTool::m_etabin
Gaudi::Property< std::vector< float > > m_etabin
selection variable for PRECISION calo selection:eta bins
Definition: TrigEgammaPrecisionPhotonHypoTool.h:33
xAOD::EgammaParameters::wtots1
@ wtots1
shower width is determined in a window detaxdphi = 0,0625 ×~0,2, corresponding typically to 20 strips...
Definition: EgammaEnums.h:140
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
TauClusterVars::dEta
bool dEta(const xAOD::TauJet &tau, const xAOD::CaloVertexedTopoCluster &cluster, float &out)
Definition: ConstituentLoaderTauCluster.cxx:114
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
SG::ReadDecorHandle
Handle class for reading a decoration on an object.
Definition: StoreGate/StoreGate/ReadDecorHandle.h:94
xAOD::EgammaParameters::f1
@ f1
E1/E = fraction of energy reconstructed in the first sampling, where E1 is energy in all strips belon...
Definition: EgammaEnums.h:52
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
TrigEgammaPrecisionPhotonHypoTool::m_dphicluster
Gaudi::Property< float > m_dphicluster
Definition: TrigEgammaPrecisionPhotonHypoTool.h:36
PlotPulseshapeFromCool.input
input
Definition: PlotPulseshapeFromCool.py:106
test_pyathena.parent
parent
Definition: test_pyathena.py:15
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:422
xAOD::EgammaParameters::Rhad1
@ Rhad1
ethad1/et
Definition: EgammaEnums.h:162
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
HLT::Identifier::fromToolName
static HLT::Identifier fromToolName(const std::string &tname)
Definition: HLTIdentifier.cxx:21
TrigEgammaPrecisionPhotonHypoTool::findCutIndex
int findCutIndex(float eta) const
Definition: TrigEgammaPrecisionPhotonHypoTool.cxx:241
TauScalarVars::absEta
bool absEta(const xAOD::TauJet &tau, float &out)
Definition: TauGNNDataLoader.cxx:118
xAOD::EgammaParameters::e277
@ e277
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
Definition: EgammaEnums.h:80
TrigEgammaPrecisionPhotonHypoTool::m_avgMuKey
SG::ReadDecorHandleKey< xAOD::EventInfo > m_avgMuKey
Definition: TrigEgammaPrecisionPhotonHypoTool.h:47
xAOD::EgammaParameters::weta1
@ weta1
shower width using +/-3 strips around the one with the maximal energy deposit: w3 strips = sqrt{sum(E...
Definition: EgammaEnums.h:97
xAOD::EgammaParameters::Eratio
@ Eratio
(emaxs1-e2tsts1)/(emaxs1+e2tsts1)
Definition: EgammaEnums.h:158
TrigEgammaPrecisionPhotonHypoTool::m_pidName
Gaudi::Property< std::string > m_pidName
Definition: TrigEgammaPrecisionPhotonHypoTool.h:41
Monitored::Scalar
Declare a monitored scalar variable.
Definition: MonitoredScalar.h:34
xAOD::EgammaParameters::Rhad
@ Rhad
ethad/et
Definition: EgammaEnums.h:160
xAOD::EgammaParameters::DeltaE
@ DeltaE
e2tsts1-emins1
Definition: EgammaEnums.h:164
xAOD::EgammaParameters::fracs1
@ fracs1
shower shape in the shower core : [E(+/-3)-E(+/-1)]/E(+/-1), where E(+/-n) is the energy in ± n strip...
Definition: EgammaEnums.h:111
read_hist_ntuple.f1
f1
Definition: read_hist_ntuple.py:4
xAOD::EgammaParameters::weta2
@ weta2
the lateral width is calculated with a window of 3x5 cells using the energy weighted sum over all cel...
Definition: EgammaEnums.h:103
Generated on Sat Nov 1 2025 21:52:21 for ATLAS Offline Software by   doxygen 1.8.18