Skip to content

Commit 4fb2f92

Browse files
hanseoparkhanseopark
committed
update multiplicity of b-jet
1 parent 3c60d82 commit 4fb2f92

File tree

1 file changed

+99
-39
lines changed

1 file changed

+99
-39
lines changed

‎PWGJE/Tasks/jetTaggerHFQA.cxx‎

Lines changed: 99 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -77,6 +77,8 @@ struct JetTaggerHFQA {
7777
Configurable<bool> checkMcCollisionIsMatched{"checkMcCollisionIsMatched", false, "0: count whole MCcollisions, 1: select MCcollisions which only have their correspond collisions"};
7878
Configurable<int> trackOccupancyInTimeRangeMax{"trackOccupancyInTimeRangeMax", 999999, "maximum occupancy of tracks in neighbouring collisions in a given time range; only applied to reconstructed collisions (data and mcd jets), not mc collisions (mcp jets)"};
7979
Configurable<int> trackOccupancyInTimeRangeMin{"trackOccupancyInTimeRangeMin", -999999, "minimum occupancy of tracks in neighbouring collisions in a given time range; only applied to reconstructed collisions (data and mcd jets), not mc collisions (mcp jets)"};
80+
Configurable<float> meanFT0A{"meanFT0A", -1., "Mean value of FT0A signal"};
81+
Configurable<float> meanFT0C{"meanFT0C", -1., "Mean value of FT0C signal"};
8082

8183
Configurable<std::string> eventSelections{"eventSelections", "sel8", "choose event selection"};
8284
Configurable<std::string> trackSelections{"trackSelections", "globalTracks", "set track selections"};
@@ -106,6 +108,7 @@ struct JetTaggerHFQA {
106108
ConfigurableAxis binMass{"binMass", {50, 0, 10.f}, ""};
107109
ConfigurableAxis binSigmaLxy{"binSigmaLxy", {100, 0., 0.1}, ""};
108110
ConfigurableAxis binSigmaLxyz{"binSigmaLxyz", {100, 0., 0.1}, ""};
111+
ConfigurableAxis binMultScaledFT0M{"binMultScaledFT0M", {VARIABLE_WIDTH, 0, 0.2, 0.3, 0.4, 0.6, 0.8, 1., 1.4, 1.8, 2.4, 3.6, 5., 20.}, "Percentiles of scaled FT0M: 100-90%, 90-80%, 80-70%, 70-60%, 60-50%, 50-40%, 40-30%, 30-20%, 20-10%, 10-1%, 1-0.1%"};
109112

110113
int numberOfJetFlavourSpecies = 6;
111114
std::vector<int> eventSelectionBits;
@@ -147,7 +150,7 @@ struct JetTaggerHFQA {
147150
AxisSpec axisSigmaLxy = {binSigmaLxy, "#sigma_{L_{XY}} [cm]"};
148151
AxisSpec axisSigmaLxyz = {binSigmaLxyz, "#sigma_{L_{XYZ}} [cm]"};
149152
AxisSpec axisFracSecPt = {100, 0, 1, "#frac{#Sigma#it{p}_{T}^{secondary track}}{#it{p}_{T, jet}}"};
150-
AxisSpec axisMult = {3500, 0, 3500, "Multiplicity"};
153+
AxisSpec axisMultScaledFT0M = {binMultScaledFT0M, "Multiplicity classes"};
151154

152155
if (doprocessTracksDca) {
153156
if (fillIPxy) {
@@ -389,17 +392,23 @@ struct JetTaggerHFQA {
389392
registry.add("h2_taggedjet_pt_3prong_mass_xyz_N1", "", {HistType::kTH2F, {{axisJetPt}, {axisMass}}});
390393
}
391394
if (doprocessSV3ProngDataMult) {
392-
registry.add("h2_jet_pt_mult", "", {HistType::kTH2F, {{axisJetPt}, {axisMult}}});
393-
registry.add("h2_jet_eta_mult", "", {HistType::kTH2F, {{axisEta}, {axisMult}}});
394-
registry.add("h2_jet_phi_mult", "", {HistType::kTH2F, {{axisPhi}, {axisMult}}});
395+
registry.add("h_event_mult", "", {HistType::kTH1F, {{axisMultScaledFT0M}}});
396+
registry.add("h2_jet_pt_mult", "", {HistType::kTH2F, {{axisJetPt}, {axisMultScaledFT0M}}});
397+
registry.add("h2_jet_eta_mult", "", {HistType::kTH2F, {{axisEta}, {axisMultScaledFT0M}}});
398+
registry.add("h2_jet_phi_mult", "", {HistType::kTH2F, {{axisPhi}, {axisMultScaledFT0M}}});
395399
if (fillGeneralSVQA) {
396-
registry.add("h2_3prong_nprongs_mult", "", {HistType::kTH2F, {{axisNprongs}, {axisMult}}});
397-
registry.add("hn_jet_3prong_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxy}, {axisSigmaLxy}, {axisSxy}, {axisLxyz}, {axisSigmaLxyz}, {axisSxyz}, {axisMult}}});
400+
registry.add("h2_3prong_nprongs_mult", "", {HistType::kTH2F, {{axisNprongs}, {axisMultScaledFT0M}}});
401+
registry.add("hn_jet_3prong_Sxy_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxy}, {axisSigmaLxy}, {axisSxy}, {axisMultScaledFT0M}}});
402+
if (fillSVxyz) {
403+
registry.add("hn_jet_3prong_Sxyz_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxyz}, {axisSigmaLxyz}, {axisSxyz}, {axisMultScaledFT0M}}});
404+
}
405+
}
406+
registry.add("hn_jet_3prong_Sxy_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMultScaledFT0M}}});
407+
registry.add("hn_taggedjet_3prong_Sxy_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMultScaledFT0M}}});
408+
if (fillSVxyz) {
409+
registry.add("hn_jet_3prong_Sxyz_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}}});
410+
registry.add("hn_taggedjet_3prong_Sxyz_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}}});
398411
}
399-
registry.add("hn_jet_3prong_Sxy_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMult}}});
400-
registry.add("hn_jet_3prong_Sxyz_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMult}}});
401-
registry.add("hn_taggedjet_3prong_Sxy_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMult}}});
402-
registry.add("hn_taggedjet_3prong_Sxyz_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMult}}});
403412
}
404413
if (doprocessSV2ProngMCD || doprocessSV2ProngMCDWeighted || doprocessSV2ProngMCPMCDMatched || doprocessSV2ProngMCPMCDMatchedWeighted) {
405414
if (!(doprocessIPsMCD || doprocessIPsMCDWeighted || doprocessIPsMCPMCDMatched || doprocessIPsMCPMCDMatchedWeighted) && !(doprocessJPMCD || doprocessJPMCDWeighted || doprocessJPMCPMCDMatched || doprocessJPMCPMCDMatchedWeighted) && !(doprocessSV3ProngMCD || doprocessSV3ProngMCDWeighted || doprocessSV3ProngMCPMCDMatched || doprocessSV3ProngMCPMCDMatchedWeighted)) {
@@ -450,17 +459,23 @@ struct JetTaggerHFQA {
450459
registry.add("h3_taggedjet_pt_3prong_mass_xyz_N1_flavour", "", {HistType::kTH3F, {{axisJetPt}, {axisMass}, {axisJetFlavour}}});
451460
}
452461
if (doprocessSV3ProngMCDMult || doprocessSV3ProngMCDMultWeighted || doprocessSV3ProngMCPMCDMatchedMult || doprocessSV3ProngMCPMCDMatchedMultWeighted) {
453-
registry.add("h3_jet_pt_mult_flavour", "", {HistType::kTH3F, {{axisJetPt}, {axisMult}, {axisJetFlavour}}});
454-
registry.add("h3_jet_eta_mult_flavour", "", {HistType::kTH3F, {{axisEta}, {axisMult}, {axisJetFlavour}}});
455-
registry.add("h3_jet_phi_mult_flavour", "", {HistType::kTH3F, {{axisPhi}, {axisMult}, {axisJetFlavour}}});
462+
registry.add("h_event_mult", "", {HistType::kTH1F, {{axisMultScaledFT0M}}});
463+
registry.add("h3_jet_pt_mult_flavour", "", {HistType::kTH3F, {{axisJetPt}, {axisMultScaledFT0M}, {axisJetFlavour}}});
464+
registry.add("h3_jet_eta_mult_flavour", "", {HistType::kTH3F, {{axisEta}, {axisMultScaledFT0M}, {axisJetFlavour}}});
465+
registry.add("h3_jet_phi_mult_flavour", "", {HistType::kTH3F, {{axisPhi}, {axisMultScaledFT0M}, {axisJetFlavour}}});
456466
if (fillGeneralSVQA) {
457-
registry.add("h3_3prong_nprongs_mult_flavour", "", {HistType::kTH3F, {{axisNprongs}, {axisMult}, {axisJetFlavour}}});
458-
registry.add("hn_jet_3prong_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxy}, {axisSigmaLxy}, {axisSxy}, {axisLxyz}, {axisSigmaLxyz}, {axisSxyz}, {axisMult}, {axisJetFlavour}}});
467+
registry.add("h3_3prong_nprongs_mult_flavour", "", {HistType::kTH3F, {{axisNprongs}, {axisMultScaledFT0M}, {axisJetFlavour}}});
468+
registry.add("hn_jet_3prong_Sxy_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxy}, {axisSigmaLxy}, {axisSxy}, {axisMultScaledFT0M}, {axisJetFlavour}}});
469+
if (fillSVxyz) {
470+
registry.add("hn_jet_3prong_Sxyz_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxyz}, {axisSigmaLxyz}, {axisSxyz}, {axisMultScaledFT0M}, {axisJetFlavour}}});
471+
}
472+
}
473+
registry.add("hn_jet_3prong_Sxy_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
474+
registry.add("hn_taggedjet_3prong_Sxy_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
475+
if (fillSVxyz) {
476+
registry.add("hn_jet_3prong_Sxyz_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
477+
registry.add("hn_taggedjet_3prong_Sxyz_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
459478
}
460-
registry.add("hn_jet_3prong_Sxy_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMult}, {axisJetFlavour}}});
461-
registry.add("hn_jet_3prong_Sxyz_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMult}, {axisJetFlavour}}});
462-
registry.add("hn_taggedjet_3prong_Sxy_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisSxyz}, {axisMass}, {axisMult}, {axisJetFlavour}}});
463-
registry.add("hn_taggedjet_3prong_Sxyz_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisSxyz}, {axisMass}, {axisMult}, {axisJetFlavour}}});
464479
}
465480
}
466481

@@ -524,6 +539,21 @@ struct JetTaggerHFQA {
524539
return true;
525540
}
526541

542+
float getScaledFT0A(const float multFT0A)
543+
{
544+
return multFT0A / meanFT0A;
545+
}
546+
547+
float getScaledFT0C(const float multFT0C)
548+
{
549+
return multFT0C / meanFT0C;
550+
}
551+
552+
float getScaledFT0M(const float multFT0A, const float multFT0C)
553+
{
554+
return 0.5 * (getScaledFT0A(multFT0A) + getScaledFT0C(multFT0C));
555+
}
556+
527557
template <typename U, typename T, typename V, typename W, typename X>
528558
void fillValidationFlavourDefMCD(T const& mcdjet, V const& tracks, W const& particles, X const& particlesPerColl, float eventWeight = 1.0)
529559
{
@@ -1021,24 +1051,29 @@ struct JetTaggerHFQA {
10211051
{
10221052
if (jet.template secondaryVertices_as<V>().size() < 1)
10231053
return;
1024-
auto mult = 0.5 * (collision.multFT0A() + collision.multFT0C());
1025-
registry.fill(HIST("h2_jet_pt_mult"), jet.pt(), mult);
1026-
registry.fill(HIST("h2_jet_eta_mult"), jet.eta(), mult);
1027-
registry.fill(HIST("h2_jet_phi_mult"), jet.phi(), mult);
1054+
float multFT0A = collision.multFT0A();
1055+
float multFT0C = collision.multFT0C();
1056+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1057+
registry.fill(HIST("h2_jet_pt_mult"), jet.pt(), scaledFT0M);
1058+
registry.fill(HIST("h2_jet_eta_mult"), jet.eta(), scaledFT0M);
1059+
registry.fill(HIST("h2_jet_phi_mult"), jet.phi(), scaledFT0M);
10281060
if (fillGeneralSVQA) {
1029-
registry.fill(HIST("h2_3prong_nprongs_mult"), jet.template secondaryVertices_as<V>().size(), mult);
1061+
registry.fill(HIST("h2_3prong_nprongs_mult"), jet.template secondaryVertices_as<V>().size(), scaledFT0M);
10301062
for (const auto& prong : jet.template secondaryVertices_as<V>()) {
1031-
registry.fill(HIST("hn_jet_3prong_mult"), jet.pt(), prong.decayLengthXY(), prong.errorDecayLengthXY(), prong.decayLengthXY() / prong.errorDecayLengthXY(), prong.decayLength(), prong.errorDecayLength(), prong.decayLength() / prong.errorDecayLengthXY(), mult);
1063+
registry.fill(HIST("hn_jet_3prong_Sxy_mult"), jet.pt(), prong.decayLengthXY(), prong.errorDecayLengthXY(), prong.decayLengthXY() / prong.errorDecayLengthXY(), scaledFT0M);
1064+
if (fillSVxyz) {
1065+
registry.fill(HIST("hn_jet_3prong_Sxyz_mult"), jet.pt(), prong.decayLength(), prong.errorDecayLength(), prong.decayLength() / prong.errorDecayLength(), scaledFT0M);
1066+
}
10321067
}
10331068
}
10341069
bool checkSv = false;
10351070
auto bjetCand = jettaggingutilities::jetFromProngMaxDecayLength<V>(jet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(2), prongCuts->at(4), prongCuts->at(5), false, &checkSv);
10361071
if (checkSv && jettaggingutilities::svAcceptance(bjetCand, svDispersionMax)) {
10371072
auto maxSxy = bjetCand.decayLengthXY() / bjetCand.errorDecayLengthXY();
10381073
auto massSV = bjetCand.m();
1039-
registry.fill(HIST("hn_jet_3prong_Sxy_N1_mult"), jet.pt(), maxSxy, massSV, mult);
1074+
registry.fill(HIST("hn_jet_3prong_Sxy_N1_mult"), jet.pt(), maxSxy, massSV, scaledFT0M);
10401075
if (jet.isTagged(BJetTaggingMethod::SV)) {
1041-
registry.fill(HIST("hn_taggedjet_3prong_Sxy_N1_mult"), jet.pt(), maxSxy, massSV, mult);
1076+
registry.fill(HIST("hn_taggedjet_3prong_Sxy_N1_mult"), jet.pt(), maxSxy, massSV, scaledFT0M);
10421077
}
10431078
}
10441079
if (fillSVxyz) {
@@ -1047,9 +1082,9 @@ struct JetTaggerHFQA {
10471082
if (checkSv && jettaggingutilities::svAcceptance(bjetCandXYZ, svDispersionMax)) {
10481083
auto maxSxyz = bjetCandXYZ.decayLength() / bjetCandXYZ.errorDecayLength();
10491084
auto massSV = bjetCandXYZ.m();
1050-
registry.fill(HIST("hn_jet_3prong_Sxyz_N1_mult"), jet.pt(), maxSxyz, massSV, mult);
1085+
registry.fill(HIST("hn_jet_3prong_Sxyz_N1_mult"), jet.pt(), maxSxyz, massSV, scaledFT0M);
10511086
if (jet.isTagged(BJetTaggingMethod::SV3D)) {
1052-
registry.fill(HIST("hn_taggedjet_3prong_Sxyz_N1_mult"), jet.pt(), maxSxyz, massSV, mult);
1087+
registry.fill(HIST("hn_taggedjet_3prong_Sxyz_N1_mult"), jet.pt(), maxSxyz, massSV, scaledFT0M);
10531088
}
10541089
}
10551090
}
@@ -1175,36 +1210,41 @@ struct JetTaggerHFQA {
11751210
return;
11761211
}
11771212
auto origin = mcdjet.origin();
1178-
auto mult = 0.5 * (collision.multFT0A() + collision.multFT0C());
1213+
float multFT0A = collision.multFT0A();
1214+
float multFT0C = collision.multFT0C();
1215+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
11791216
if (mcdjet.template secondaryVertices_as<V>().size() < 1)
11801217
return;
1181-
registry.fill(HIST("h3_jet_pt_mult_flavour"), mcdjet.pt(), mult, origin, eventWeight);
1182-
registry.fill(HIST("h3_jet_eta_mult_flavour"), mcdjet.eta(), mult, origin, eventWeight);
1183-
registry.fill(HIST("h3_jet_phi_mult_flavour"), mcdjet.phi(), mult, origin, eventWeight);
1218+
registry.fill(HIST("h3_jet_pt_mult_flavour"), mcdjet.pt(), scaledFT0M, origin, eventWeight);
1219+
registry.fill(HIST("h3_jet_eta_mult_flavour"), mcdjet.eta(), scaledFT0M, origin, eventWeight);
1220+
registry.fill(HIST("h3_jet_phi_mult_flavour"), mcdjet.phi(), scaledFT0M, origin, eventWeight);
11841221
if (fillGeneralSVQA) {
1185-
registry.fill(HIST("h3_3prong_nprongs_mult_flavour"), mcdjet.template secondaryVertices_as<V>().size(), mult, origin, eventWeight);
1222+
registry.fill(HIST("h3_3prong_nprongs_mult_flavour"), mcdjet.template secondaryVertices_as<V>().size(), scaledFT0M, origin, eventWeight);
11861223
for (const auto& prong : mcdjet.template secondaryVertices_as<V>()) {
1187-
registry.fill(HIST("hn_jet_3prong_mult_flavour"), mcdjet.pt(), prong.decayLengthXY(), prong.errorDecayLengthXY(), prong.decayLengthXY() / prong.errorDecayLengthXY(), prong.decayLength(), prong.errorDecayLengthXY(), prong.decayLength() / prong.errorDecayLengthXY(), origin, eventWeight);
1224+
registry.fill(HIST("hn_jet_3prong_Sxy_mult_flavour"), mcdjet.pt(), prong.decayLengthXY(), prong.errorDecayLengthXY(), prong.decayLengthXY() / prong.errorDecayLengthXY(), scaledFT0M, origin, eventWeight);
1225+
if (fillSVxyz) {
1226+
registry.fill(HIST("hn_jet_3prong_Sxyz_mult_flavour"), mcdjet.pt(), prong.decayLength(), prong.errorDecayLength(), prong.decayLength() / prong.errorDecayLength(), scaledFT0M, origin, eventWeight);
1227+
}
11881228
}
11891229
}
11901230
bool checkSv = false;
11911231
auto bjetCand = jettaggingutilities::jetFromProngMaxDecayLength<V>(mcdjet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(2), prongCuts->at(4), prongCuts->at(5), false, &checkSv);
11921232
if (checkSv && jettaggingutilities::svAcceptance(bjetCand, svDispersionMax)) {
11931233
auto maxSxy = bjetCand.decayLengthXY() / bjetCand.errorDecayLengthXY();
11941234
auto massSV = bjetCand.m();
1195-
registry.fill(HIST("hn_jet_3prong_Sxy_N1_mult_flavour"), mcdjet.pt(), maxSxy, massSV, mult, origin, eventWeight);
1235+
registry.fill(HIST("hn_jet_3prong_Sxy_N1_mult_flavour"), mcdjet.pt(), maxSxy, massSV, scaledFT0M, origin, eventWeight);
11961236
if (mcdjet.isTagged(BJetTaggingMethod::SV)) {
1197-
registry.fill(HIST("hn_taggedjet_3prong_Sxy_N1_mult_flavour"), mcdjet.pt(), maxSxy, massSV, mult, origin, eventWeight);
1237+
registry.fill(HIST("hn_taggedjet_3prong_Sxy_N1_mult_flavour"), mcdjet.pt(), maxSxy, massSV, scaledFT0M, origin, eventWeight);
11981238
}
11991239
}
12001240
if (fillSVxyz) {
12011241
auto bjetCandXYZ = jettaggingutilities::jetFromProngMaxDecayLength<V>(mcdjet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(3), prongCuts->at(4), prongCuts->at(5), true, &checkSv);
12021242
if (checkSv && jettaggingutilities::svAcceptance(bjetCandXYZ, svDispersionMax)) {
12031243
auto maxSxyz = bjetCandXYZ.decayLength() / bjetCandXYZ.errorDecayLength();
12041244
auto massSV = bjetCandXYZ.m();
1205-
registry.fill(HIST("hn_jet_3prong_Sxyz_N1_mult_flavour"), mcdjet.pt(), maxSxyz, massSV, mult, origin, eventWeight);
1245+
registry.fill(HIST("hn_jet_3prong_Sxyz_N1_mult_flavour"), mcdjet.pt(), maxSxyz, massSV, scaledFT0M, origin, eventWeight);
12061246
if (mcdjet.isTagged(BJetTaggingMethod::SV3D)) {
1207-
registry.fill(HIST("hn_taggedjet_3prong_Sxyz_N1_mult_flavour"), mcdjet.pt(), maxSxyz, massSV, mult, origin, eventWeight);
1247+
registry.fill(HIST("hn_taggedjet_3prong_Sxyz_N1_mult_flavour"), mcdjet.pt(), maxSxyz, massSV, scaledFT0M, origin, eventWeight);
12081248
}
12091249
}
12101250
}
@@ -1620,6 +1660,10 @@ struct JetTaggerHFQA {
16201660
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
16211661
return;
16221662
}
1663+
float multFT0A = collision.multFT0A();
1664+
float multFT0C = collision.multFT0C();
1665+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1666+
registry.fill(HIST("h_event_mult"), scaledFT0M);
16231667
for (auto const& jet : jets) {
16241668
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
16251669
continue;
@@ -1789,6 +1833,10 @@ struct JetTaggerHFQA {
17891833
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
17901834
return;
17911835
}
1836+
float multFT0A = collision.multFT0A();
1837+
float multFT0C = collision.multFT0C();
1838+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1839+
registry.fill(HIST("h_event_mult"), scaledFT0M);
17921840
for (auto const& mcdjet : mcdjets) {
17931841
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
17941842
continue;
@@ -1806,6 +1854,10 @@ struct JetTaggerHFQA {
18061854
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
18071855
return;
18081856
}
1857+
float multFT0A = collision.multFT0A();
1858+
float multFT0C = collision.multFT0C();
1859+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1860+
registry.fill(HIST("h_event_mult"), scaledFT0M);
18091861
for (auto const& mcdjet : mcdjets) {
18101862
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
18111863
continue;
@@ -1823,6 +1875,10 @@ struct JetTaggerHFQA {
18231875
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
18241876
return;
18251877
}
1878+
float multFT0A = collision.multFT0A();
1879+
float multFT0C = collision.multFT0C();
1880+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1881+
registry.fill(HIST("h_event_mult"), scaledFT0M);
18261882
for (auto const& mcdjet : mcdjets) {
18271883
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
18281884
continue;
@@ -1844,6 +1900,10 @@ struct JetTaggerHFQA {
18441900
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
18451901
return;
18461902
}
1903+
float multFT0A = collision.multFT0A();
1904+
float multFT0C = collision.multFT0C();
1905+
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1906+
registry.fill(HIST("h_event_mult"), scaledFT0M);
18471907
for (auto const& mcdjet : mcdjets) {
18481908
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
18491909
continue;

0 commit comments

Comments
 (0)