34 const ROOT::Math::XYZVector& resultMomentum = result->getMomentum();
35 var<
named(
"pt")>() = resultMomentum.Rho();
36 var<
named(
"theta")>() = resultMomentum.Theta();
38 const std::vector<const SpacePoint*>& spacePoints = result->getHits();
40 var<
named(
"number_of_hits")>() = spacePoints.size();
42 std::vector<bool> layerUsed;
43 layerUsed.resize(7,
false);
45 for (
const SpacePoint* spacePoint : spacePoints) {
46 layerUsed[spacePoint->getVxdID().getLayerNumber()] =
true;
52 var<
named(
"number_of_holes")>() = std::count(layerUsed.begin(), layerUsed.end(),
true);
54 var<
named(
"has_missing_layer_1")>() = layerUsed[1] == 0;
55 var<
named(
"has_missing_layer_2")>() = layerUsed[2] == 0;
56 var<
named(
"has_missing_layer_3")>() = layerUsed[3] == 0;
57 var<
named(
"has_missing_layer_4")>() = layerUsed[4] == 0;
58 var<
named(
"has_missing_layer_5")>() = layerUsed[5] == 0;
59 var<
named(
"has_missing_layer_6")>() = layerUsed[6] == 0;
65 var<
named(
"number_of_overlap_hits")>() = std::count(layerUsed.begin(), layerUsed.end(), 2);
67 if (spacePoints.empty()) {
69 var<
named(
"first_hit_layer")>() = -1;
71 var<
named(
"last_hit_layer")>() = spacePoints.back()->getVxdID().getLayerNumber();
72 var<
named(
"first_hit_layer")>() = spacePoints.front()->getVxdID().getLayerNumber();
75 genfit::MeasuredStateOnPlane mSoP = result->getMSoP();
76 const genfit::MeasuredStateOnPlane& firstCDCHit = result->getSeedMSoP();
77 m_advancer.extrapolateToPlane(mSoP, firstCDCHit.getPlane());
79 const auto& distance = mSoP.getPos() - firstCDCHit.getPos();
80 var<
named(
"distance_to_seed_track")>() = distance.Mag();
81 var<
named(
"distance_to_seed_track_xy")>() = distance.Pt();
83 const RecoTrack* seedTrack = result->getSeed();
85 var<
named(
"chi2")>() = result->getChi2();
86 var<
named(
"chi2_vxd_max")>() = result->getMaximalChi2();
87 var<
named(
"chi2_vxd_min")>() = result->getMinimalChi2();
89 var<
named(
"weight_sum")>() = result->getWeightSum();