@@ -109,19 +109,19 @@ struct Config {
109109};
110110constexpr int NConfigs = 13 ;
111111constexpr int INominal = 0 ;
112- constexpr std::array<Config, NConfigs> Configs = {{{" NOM" , 1 , 3 , kCondOwn , 1 , 3 , kV1 },
113- {" GAP02" , 1 , 3 , kCondAll , 4 , 3 , kV1 },
114- {" GAP03" , 2 , 3 , kCondAll , 4 , 3 , kV1 },
115- {" GAP04" , 3 , 3 , kCondAll , 4 , 3 , kV1 },
116- {" GAP05" , 4 , 3 , kCondAll , 4 , 3 , kV1 },
117- {" PT10" , 1 , 1 , kCondUnwindowed , 1 , 1 , kV1 },
118- {" PT15" , 1 , 2 , kCondUnwindowed , 1 , 1 , kV1 },
119- {" PT20" , 1 , 3 , kCondUnwindowed , 1 , 1 , kV1 },
120- {" PT30" , 1 , 4 , kCondUnwindowed , 1 , 1 , kV1 },
121- {" PT50" , 1 , 5 , kCondUnwindowed , 1 , 1 , kV1 },
122- {" EVV0" , 1 , 3 , kCondOwn , 1 , 3 , kV0 },
123- {" EVV3" , 1 , 3 , kCondOwn , 1 , 3 , kV3 },
124- {" EVV6" , 1 , 3 , kCondOwn , 1 , 3 , kV6 }}};
112+ constexpr std::array<Config, NConfigs> Configs = {{{. name = " NOM" , . etaFirst = 1 , . ptEnd = 3 , . cond = kCondOwn , . validEtaFirst = 1 , . validPtEnd = 3 , . evVariant = kV1 },
113+ {. name = " GAP02" , . etaFirst = 1 , . ptEnd = 3 , . cond = kCondAll , . validEtaFirst = 4 , . validPtEnd = 3 , . evVariant = kV1 },
114+ {. name = " GAP03" , . etaFirst = 2 , . ptEnd = 3 , . cond = kCondAll , . validEtaFirst = 4 , . validPtEnd = 3 , . evVariant = kV1 },
115+ {. name = " GAP04" , . etaFirst = 3 , . ptEnd = 3 , . cond = kCondAll , . validEtaFirst = 4 , . validPtEnd = 3 , . evVariant = kV1 },
116+ {. name = " GAP05" , . etaFirst = 4 , . ptEnd = 3 , . cond = kCondAll , . validEtaFirst = 4 , . validPtEnd = 3 , . evVariant = kV1 },
117+ {. name = " PT10" , . etaFirst = 1 , . ptEnd = 1 , . cond = kCondUnwindowed , . validEtaFirst = 1 , . validPtEnd = 1 , . evVariant = kV1 },
118+ {. name = " PT15" , . etaFirst = 1 , . ptEnd = 2 , . cond = kCondUnwindowed , . validEtaFirst = 1 , . validPtEnd = 1 , . evVariant = kV1 },
119+ {. name = " PT20" , . etaFirst = 1 , . ptEnd = 3 , . cond = kCondUnwindowed , . validEtaFirst = 1 , . validPtEnd = 1 , . evVariant = kV1 },
120+ {. name = " PT30" , . etaFirst = 1 , . ptEnd = 4 , . cond = kCondUnwindowed , . validEtaFirst = 1 , . validPtEnd = 1 , . evVariant = kV1 },
121+ {. name = " PT50" , . etaFirst = 1 , . ptEnd = 5 , . cond = kCondUnwindowed , . validEtaFirst = 1 , . validPtEnd = 1 , . evVariant = kV1 },
122+ {. name = " EVV0" , . etaFirst = 1 , . ptEnd = 3 , . cond = kCondOwn , . validEtaFirst = 1 , . validPtEnd = 3 , . evVariant = kV0 },
123+ {. name = " EVV3" , . etaFirst = 1 , . ptEnd = 3 , . cond = kCondOwn , . validEtaFirst = 1 , . validPtEnd = 3 , . evVariant = kV3 },
124+ {. name = " EVV6" , . etaFirst = 1 , . ptEnd = 3 , . cond = kCondOwn , . validEtaFirst = 1 , . validPtEnd = 3 , . evVariant = kV6 }}};
125125
126126constexpr int NCountBins = 64 ; // larger counts in the overflow bin
127127constexpr int NSubCountBins = 32 ; // per sub-event
@@ -221,8 +221,8 @@ void monomials(double xF, double xB, double vF, double vB, int nMono, std::array
221221// / Moment sums of one reconstructed half or of the generated events
222222struct Sums {
223223 std::array<std::shared_ptr<TH3 >, NConfigs> config{}; // (N_fwd, stratum count, monomial)
224- std::shared_ptr<THn> strata3{}; // (N_fwd, n_F, n_B, first-order monomial)
225- std::shared_ptr<TH3 > occupancy{}; // (configuration, N_fwd, stratum count)
224+ std::shared_ptr<THn> strata3; // (N_fwd, n_F, n_B, first-order monomial)
225+ std::shared_ptr<TH3 > occupancy; // (configuration, N_fwd, stratum count)
226226};
227227
228228// / Nominal-configuration observables of an event
@@ -274,9 +274,18 @@ struct EventShapeCoex {
274274 const AxisSpec monoAxis{NMono, -0.5 , NMono - 0.5 , " monomial" };
275275 const AxisSpec mono1Axis{NMono1, -0.5 , NMono1 - 0.5 , " monomial" };
276276 const AxisSpec configAxis{NConfigs, -0.5 , NConfigs - 0.5 , " configuration" };
277- const std::string label = " [half " + half + " ] moment sums " ;
277+ std::string label = " [half " ;
278+ label += half;
279+ label += " ] moment sums " ;
278280 for (int c = 0 ; c < nConfigs; ++c) {
279- sums.config [c] = registry.add <TH3 >(dir + " /" + Configs[c].name + " /" + half, (label + Configs[c].name ).c_str (), HistType::kTH3D , {nFwdAxis, countAxis, monoAxis});
281+ std::string name = dir;
282+ name += " /" ;
283+ name += Configs[c].name ;
284+ name += " /" ;
285+ name += half;
286+ std::string title = label;
287+ title += Configs[c].name ;
288+ sums.config [c] = registry.add <TH3 >(name, title.c_str (), HistType::kTH3D , {nFwdAxis, countAxis, monoAxis});
280289 sums.config [c]->SetBit (TH1 ::kIsNotW ); // no Sumw2
281290 }
282291 sums.strata3 = registry.add <THn>(dir + " /NOM3/" + half, (label + " NOM in (n_F, n_B) strata" ).c_str (), HistType::kTHnD , {nFwdAxis, subCountAxis, subCountAxis, mono1Axis});
@@ -362,7 +371,7 @@ struct EventShapeCoex {
362371 }
363372 sums.occupancy ->Fill (c, nFwd, nCond);
364373 if (c == INominal) {
365- nominal = {true , nCond, xF, xB, vF, vB};
374+ nominal = {. valid = true , . nMid = nCond, . xF = xF, . xB = xB, . vF = vF, . vB = vB};
366375 }
367376 }
368377
0 commit comments