Rretrieve basic parameters, mainly four-vectors.
114 {
115
117
119
124
125
126
127
134
137
138
141 DataVect qualityTile; qualityTile.reserve(
jets->size());
143 DataVect timeClusters; timeClusters.reserve(
jets->size());
152 DataVect fcorJetForCell; fcorJetForCell.reserve(
jets->size());
154 DataVect fracSamplingMax; fracSamplingMax.reserve(
jets->size());
156 DataVect OutOfTimeEfrac; OutOfTimeEfrac.reserve(
jets->size());
161
162 int id = 0;
163
165 for (; itr !=
jets->end(); ++itr) {
168 et.push_back(
DataType((*itr)->et()/CLHEP::GeV));
169 pt.push_back(
DataType((*itr)->pt()/CLHEP::GeV));
171
172
173
174
175 bTagName.push_back(
DataType(
"JetFitterCOMBNN" ));
176 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"JetFitterCOMBNN") ));
177 bTagName.push_back(
DataType(
"JetFitterTagNN" ));
178 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"JetFitterTagNN") ));
179 bTagName.push_back(
DataType(
"IP3D+SV1" ));
180 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight() ));
181 bTagName.push_back(
DataType(
"IP2D" ));
182 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"IP2D") ));
183 bTagName.push_back(
DataType(
"IP3D" ));
184 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"IP3D") ));
185 bTagName.push_back(
DataType(
"SV1" ));
186 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"SV1") ));
187 bTagName.push_back(
DataType(
"SV2" ));
188 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"SV2") ));
189 bTagName.push_back(
DataType(
"MV1" ));
190 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"MV1") ));
191 bTagName.push_back(
DataType(
"MV2" ));
192 bTagValue.push_back(
DataType( (*itr)->getFlavourTagWeight(
"MV2") ));
193
194
195 quality.push_back(
DataType((*itr)->getMoment(
"LArQuality")));
196
201
202
203 jvf.push_back(
DataType((*itr)->getMoment(
"JVF") ));
204
207 px.push_back(
DataType((*itr)->px()/CLHEP::GeV ) );
208 py.push_back(
DataType((*itr)->py()/CLHEP::GeV ) );
209 pz.push_back(
DataType((*itr)->pz()/CLHEP::GeV ) );
210 }
211
212
213 const auto n =
phi.size();
218 DataMap[
"id"] = std::move(idVec);
219
220 DataMap[
"bTagName multiple=\"9\""] = std::move(bTagName);
221 DataMap[
"bTagValue multiple=\"9\""] = std::move(bTagValue);
222
223
224 DataMap[
"quality"] = std::move(quality);
225 DataMap[
"isGood"] = std::move(isGood);
226 DataMap[
"isBad"] = std::move(isBad);
227 DataMap[
"isUgly"] = std::move(isUgly);
228 DataMap[
"emfrac"] = std::move(emfrac);
229
230 DataMap[
"jvf"] = std::move(jvf);
231
233 DataMap[
"qualityLAr"] = std::move(qualityLAr);
234 DataMap[
"qualityTile"] = std::move(qualityTile);
235 DataMap[
"time"] = std::move(time);
236 DataMap[
"timeClusters"] = std::move(timeClusters);
237 DataMap[
"n90cells"] = std::move(n90cells);
238 DataMap[
"n90const"] = std::move(n90const);
239 DataMap[
"hecf"] = std::move(hecf);
240 DataMap[
"tileGap3f"] = std::move(tileGap3f);
241 DataMap[
"fcorCell"] = std::move(fcorCell);
242 DataMap[
"fcorDotx"] = std::move(fcorDotx);
243 DataMap[
"fcorJet"] = std::move(fcorJet);
244 DataMap[
"fcorJetForCell"] = std::move(fcorJetForCell);
245 DataMap[
"nbadcells"] = std::move(nbadcells);
246 DataMap[
"fracSamplingMax"] = std::move(fracSamplingMax);
247 DataMap[
"sMax"] = std::move(sMax);
248 DataMap[
"OutOfTimeEfrac"] = std::move(OutOfTimeEfrac);
249 }
250
251
252
253
254 DataMap[
"mass"] = std::move(mass);
258 DataMap[
"energy"] = std::move(energy);
259
262 msg(MSG::DEBUG) <<
" retrieved with " <<
n <<
" entries"<<
endmsg;
263 }
264
265
267
268 }
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
float et(const xAOD::jFexSRJetRoI *j)
bool msgLvl(const MSG::Level lvl) const
DataModel_detail::const_iterator< DataVector > const_iterator
virtual std::string dataTypeName() const
Return the name of the data type.
time(flags, cells_name, *args, **kw)
std::map< std::string, DataVect > DataMap
std::vector< DataType > DataVect
Defines a map with a key and a vector of DataType objects e.g.