73{
74 auto mon_ET = Monitored::Scalar( "Et_em", -1.0 );
75 auto mon_etaBin = Monitored::Scalar( "EtaBin", -1.0 );
76 auto mon_Eta = Monitored::Scalar( "Eta", -99. );
77 auto mon_Phi = Monitored::Scalar( "Phi", -99. );
78 auto mon_mu = Monitored::Scalar("mu", -1.);
79 auto mon_etcone20 = Monitored::Scalar("etcone20", -99.);
80 auto mon_topoetcone20 = Monitored::Scalar("topoetcone20", -99.);
81 auto mon_relEtCone20 = Monitored::Scalar("relEtCone20", -99.);
82 auto mon_relTopoEtCone20 = Monitored::Scalar("relTopoEtCone20", -99.);
83
84 auto mon_etcone30 = Monitored::Scalar("etcone30", -99.);
85 auto mon_topoetcone30 = Monitored::Scalar("topoetcone30", -99.);
86 auto mon_relEtCone30 = Monitored::Scalar("relEtCone30", -99.);
87 auto mon_relTopoEtCone30 = Monitored::Scalar("relTopoEtCone30", -99.);
88
89 auto mon_etcone40 = Monitored::Scalar("etcone40", -99.);
90 auto mon_topoetcone40 = Monitored::Scalar("topoetcone40", -99.);
91 auto mon_relEtCone40 = Monitored::Scalar("relEtCone40", -99.);
92 auto mon_relTopoEtCone40 = Monitored::Scalar("relTopoEtCone40", -99.);
93
94 auto PassedCuts = Monitored::Scalar<int>( "CutCounter", -1 );
95 auto monitorIt = Monitored::Group(
m_monTool,
96 mon_etaBin, mon_Eta, mon_Phi, mon_mu,
97 mon_etcone20, mon_topoetcone20, mon_relEtCone20, mon_relTopoEtCone20,
98 mon_etcone30, mon_topoetcone30, mon_relEtCone30, mon_relTopoEtCone30,
99 mon_etcone40, mon_topoetcone40, mon_relEtCone40, mon_relTopoEtCone40,
100 PassedCuts );
101
102
103 PassedCuts = PassedCuts + 1;
104
105 auto roiDescriptor =
input.roi;
106
107 if ( fabs( roiDescriptor->eta() ) > 2.6 ) {
108 ATH_MSG_DEBUG(
"REJECT The photon had eta coordinates beyond the EM fiducial volume : "
109 << roiDescriptor->eta() << "; stop the chain now" );
110 return false;
111 }
112
114 << ": Eta = " << roiDescriptor->eta()
115 << ", Phi = " << roiDescriptor->phi() );
116
117
118 const auto pClus =
input.photon->caloCluster();
119
120 const float absEta = fabs( pClus->eta() );
122
123
126 return false;
127 } else {
129 }
131 PassedCuts = PassedCuts + 1;
132
133 const float photon_eT = pClus->et();
134 mon_ET = photon_eT;
135
136
137 float avg_mu = 0;
138 SG::ReadDecorHandle<xAOD::EventInfo,float> eventInfoDecor(
m_avgMuKey);
139 if(eventInfoDecor.isPresent()) {
140 avg_mu = eventInfoDecor(0);
142 }
143 mon_mu = avg_mu;
144
145 float ptCone[
s_nCones] = {999, 999, 999};
146 float etCone[
s_nCones] = {999, 999, 999};
147 float topoEtCone[
s_nCones] = {999, 999, 999};
148
149 float relEtCone[
s_nCones] = {999, 999, 999};
150 float relTopoEtCone[
s_nCones] = {999, 999, 999};
151
156
160 }
161
162
164 relEtCone[
idx] = etCone[
idx] / photon_eT;
165 relTopoEtCone[
idx] = topoEtCone[
idx] / photon_eT;
166
169 }
170
171
172 mon_etcone20 = etCone[0];
173 mon_topoetcone20 = topoEtCone[0];
174 mon_relEtCone20 = relEtCone[0];
175 mon_relTopoEtCone20 = relTopoEtCone[0];
176
177 mon_etcone30 = etCone[1];
178 mon_topoetcone30 = topoEtCone[1];
179 mon_relEtCone30 = relEtCone[1];
180 mon_relTopoEtCone30 = relTopoEtCone[1];
181
182 mon_etcone40 = etCone[2];
183 mon_topoetcone40 = topoEtCone[2];
184 mon_relEtCone40 = relEtCone[2];
185 mon_relTopoEtCone40 = relTopoEtCone[2];
186
188 return true;
189 }
190
196 );
197
198 return pass_relTopoEtCone;
199
200}
bool absEta(const xAOD::TauJet &tau, float &out)
1.17.0