60 unsigned int nEntries(6);
62 ROOT::RDataFrame rdf(nEntries);
64 std::vector<double> pdgIds(nEntries,
static_cast<double>(
m_testHypo.getPDGCode()));
65 std::vector<double> pMinEdges = {1.0, 1.0, 1.0, 1.5, 1.5, 1.5};
66 std::vector<double> pMaxEdges = {1.5, 1.5, 1.5, 3.0, 3.0, 3.0};
67 std::vector<double> thetaMinEdges = {0.22, 0.56, 2.23, 0.22, 0.56, 2.23};
68 std::vector<double> thetaMaxEdges = {0.56, 2.23, 2.71, 0.56, 2.23, 2.71};
69 std::vector<double> pBinIdxs = {1.0, 1.0, 1.0, 2.0, 2.0, 2.0};
70 std::vector<double> thetaBinIdxs = {1.0, 2.0, 3.0, 1.0, 2.0, 3.0};
71 std::vector<std::string> regLabels = {
"FWD",
"Barrel",
"BWD",
"FWD",
"Barrel",
"BWD"};
72 std::vector<double> weights_SVD = {0., 0., -0.003, 0.031, 0.017, 0.077};
73 std::vector<double> weights_CDC = {-0.061, -0.084, -0.162, -0.073, -0.073, -0.05, -0.121};
74 std::vector<double> weights_TOP = {0., 0., 0., 0., 0., 0.};
75 std::vector<double> weights_ARICH = {0., 0., 0., -0.004, 0., 0.};
76 std::vector<double> weights_ECL = {-0.062, -0.118, -0.136, -0.22, -0.255, -0.377};
77 std::vector<double> weights_KLM = {0., 0., 0.006, -0.003, 0.003, 0.002};
79 unsigned int iEntry(0);
80 rdf.Define(
"pdgId", [&]() {
auto x = pdgIds[iEntry];
return x; })
81 .Define(
"p_min", [&]() {
auto x = pMinEdges[iEntry];
return x; })
82 .Define(
"p_max", [&]() {
auto x = pMaxEdges[iEntry];
return x; })
83 .Define(
"theta_min", [&]() {
auto x = thetaMinEdges[iEntry];
return x; })
84 .Define(
"theta_max", [&]() {
auto x = thetaMaxEdges[iEntry];
return x; })
85 .Define(
"p_bin_idx", [&]() {
auto x = pBinIdxs[iEntry];
return x; })
86 .Define(
"theta_bin_idx", [&]() {
auto x = thetaBinIdxs[iEntry];
return x; })
87 .Define(
"reg_label", [&]() {
auto x = regLabels[iEntry];
return x; })
88 .Define(
"ablat_s_SVD", [&]() {
auto x = weights_SVD[iEntry];
return x; })
89 .Define(
"ablat_s_CDC", [&]() {
auto x = weights_CDC[iEntry];
return x; })
90 .Define(
"ablat_s_TOP", [&]() {
auto x = weights_TOP[iEntry];
return x; })
91 .Define(
"ablat_s_ARICH", [&]() {
auto x = weights_ARICH[iEntry];
return x; })
92 .Define(
"ablat_s_ECL", [&]() {
auto x = weights_ECL[iEntry];
return x; })
93 .Define(
"ablat_s_KLM", [&]() {
auto x = weights_KLM[iEntry]; ++iEntry;
return x; })
105 unsigned int nEntries(6);
107 ROOT::RDataFrame rdf(nEntries);
109 std::vector<double> pdgIds(nEntries,
static_cast<double>(
m_testHypo.getPDGCode()));
110 std::vector<double> pMinEdges = {1.0, 1.0, 1.0, 1.5, 1.5, 1.5};
111 std::vector<double> pMaxEdges = {1.8, 1.8, 1.8, 3.0, 3.0, 3.0};
112 std::vector<double> thetaMinEdges = {0.22, 0.56, 2.23, 0.22, 0.56, 2.23};
113 std::vector<double> thetaMaxEdges = {0.56, 2.23, 2.71, 0.56, 2.23, 2.71};
114 std::vector<double> pBinIdxs = {1.0, 1.0, 1.0, 2.0, 2.0, 2.0};
115 std::vector<double> thetaBinIdxs = {1.0, 2.0, 3.0, 1.0, 2.0, 3.0};
118 ROOT::RDF::RSnapshotOptions opt;
119 opt.fMode =
"UPDATE";
121 unsigned int iEntry(0);
122 rdf.Define(
"pdgId", [&]() {
auto x = pdgIds[iEntry];
return x; })
123 .Define(
"p_min", [&]() {
auto x = pMinEdges[iEntry];
return x; })
124 .Define(
"p_max", [&]() {
auto x = pMaxEdges[iEntry];
return x; })
125 .Define(
"theta_min", [&]() {
auto x = thetaMinEdges[iEntry];
return x; })
126 .Define(
"theta_max", [&]() {
auto x = thetaMaxEdges[iEntry];
return x; })
127 .Define(
"p_bin_idx", [&]() {
auto x = pBinIdxs[iEntry];
return x; })
128 .Define(
"theta_bin_idx", [&]() {
auto x = thetaBinIdxs[iEntry]; ++iEntry;
return x; })