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
82 PassedCuts = PassedCuts + 1;
83
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;
92 return pass;
93 }
94
96 << ": Eta = " << roiDescriptor->eta()
97 << ", Phi = " << roiDescriptor->phi() );
98
99
100 double etaRef = roiDescriptor->eta();
101 double phiRef = roiDescriptor->phi();
102
103 if ( fabs( phiRef ) >
M_PI ) phiRef -= 2*
M_PI;
104
105 auto pClus =
input.photon->caloCluster();
106
107 float absEta = fabs( pClus->eta() );
109
110 dEta = pClus->eta() - etaRef;
111
112 dPhi = fabs( pClus->phi() - phiRef );
114 ET = pClus->et();
117
119 << " roi eta=" << etaRef << " DeltaEta=" << dEta
121
124 return pass;
125 }
128 PassedCuts = PassedCuts + 1;
129
130
132 << " roi phi="<< phiRef << " DeltaPhi="<< dPhi
134
137 return pass;
138 }
141 PassedCuts = PassedCuts + 1;
142
143
144 if ( cutIndex == -1 ) {
146 return pass;
147 } else {
149 }
151 PassedCuts = PassedCuts + 1;
152
153
155 if ( ET <
m_eTthr[cutIndex] ) {
157 return pass;
158 }
159 mon_ET = ET;
160 PassedCuts = PassedCuts + 1;
162 pass = true;
163 return pass;
164 }
165
168 }
169
170
171 float avg_mu = 0;
172 SG::ReadDecorHandle<xAOD::EventInfo,float> eventInfoDecor(
m_avgMuKey);
173 if(eventInfoDecor.isPresent()) {
174 avg_mu = eventInfoDecor(0);
176 }
177 mon_mu = avg_mu;
178
181
182
183
184
186
188
189
190
192
194
196
198
199
200
202
204
205
206
208
210
212
214
227
228 if ( !pass ){
230 return pass;
231 } else {
233 }
234
236
237 return pass;
238
239}
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
bool dPhi(const xAOD::TauJet &tau, const xAOD::CaloVertexedTopoCluster &cluster, float &out)
bool dEta(const xAOD::TauJet &tau, const xAOD::CaloVertexedTopoCluster &cluster, float &out)
bool absEta(const xAOD::TauJet &tau, float &out)
@ wtots1
shower width is determined in a window detaxdphi = 0,0625 ×~0,2, corresponding typically to 20 strips...
@ e277
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
@ f1
E1/E = fraction of energy reconstructed in the first sampling, where E1 is energy in all strips belon...
@ Eratio
(emaxs1-e2tsts1)/(emaxs1+e2tsts1)
@ fracs1
shower shape in the shower core : [E(+/-3)-E(+/-1)]/E(+/-1), where E(+/-n) is the energy in ± n strip...
@ weta2
the lateral width is calculated with a window of 3x5 cells using the energy weighted sum over all cel...
@ weta1
shower width using +/-3 strips around the one with the maximal energy deposit: w3 strips = sqrt{sum(E...
1.14.0