Skip to content

Commit fb23a9a

Browse files
hanseoparkhanseopark
committed
delete multiplicity analysis
1 parent d40ba5e commit fb23a9a

File tree

1 file changed

+0
-268
lines changed

1 file changed

+0
-268
lines changed

‎PWGJE/Tasks/jetTaggerHFQA.cxx‎

Lines changed: 0 additions & 268 deletions
Original file line numberDiff line numberDiff line change
@@ -21,8 +21,6 @@
2121
#include "PWGJE/DataModel/JetReducedData.h"
2222
#include "PWGJE/DataModel/JetTagging.h"
2323

24-
#include "Common/DataModel/Multiplicity.h"
25-
2624
#include "Framework/ASoA.h"
2725
#include "Framework/AnalysisDataModel.h"
2826
#include "Framework/AnalysisTask.h"
@@ -109,7 +107,6 @@ struct JetTaggerHFQA {
109107
ConfigurableAxis binMass{"binMass", {50, 0, 10.f}, ""};
110108
ConfigurableAxis binSigmaLxy{"binSigmaLxy", {100, 0., 0.1}, ""};
111109
ConfigurableAxis binSigmaLxyz{"binSigmaLxyz", {100, 0., 0.1}, ""};
112-
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%"};
113110

114111
int numberOfJetFlavourSpecies = 6;
115112
std::vector<int> eventSelectionBits;
@@ -151,7 +148,6 @@ struct JetTaggerHFQA {
151148
AxisSpec axisSigmaLxy = {binSigmaLxy, "#sigma_{L_{XY}} [cm]"};
152149
AxisSpec axisSigmaLxyz = {binSigmaLxyz, "#sigma_{L_{XYZ}} [cm]"};
153150
AxisSpec axisFracSecPt = {100, 0, 1, "#frac{#Sigma#it{p}_{T}^{secondary track}}{#it{p}_{T, jet}}"};
154-
AxisSpec axisMultScaledFT0M = {binMultScaledFT0M, "Multiplicity classes"};
155151

156152
registry.add("h_collision_events", "data;mcd;mcp evnets", {HistType::kTH1F, {{4, 0.0, 4.0}}});
157153
if (doprocessTracksDca) {
@@ -397,25 +393,6 @@ struct JetTaggerHFQA {
397393
registry.add("h2_taggedjet_pt_3prong_mass_N1", "", {HistType::kTH2F, {{axisJetPt}, {axisMass}}});
398394
registry.add("h2_taggedjet_pt_3prong_mass_xyz_N1", "", {HistType::kTH2F, {{axisJetPt}, {axisMass}}});
399395
}
400-
if (doprocessSV3ProngDataMult) {
401-
registry.add("h_event_mult", "", {HistType::kTH1F, {{axisMultScaledFT0M}}});
402-
registry.add("h2_jet_pt_mult", "", {HistType::kTH2F, {{axisJetPt}, {axisMultScaledFT0M}}});
403-
registry.add("h2_jet_eta_mult", "", {HistType::kTH2F, {{axisEta}, {axisMultScaledFT0M}}});
404-
registry.add("h2_jet_phi_mult", "", {HistType::kTH2F, {{axisPhi}, {axisMultScaledFT0M}}});
405-
if (fillGeneralSVQA) {
406-
registry.add("h2_3prong_nprongs_mult", "", {HistType::kTH2F, {{axisNprongs}, {axisMultScaledFT0M}}});
407-
registry.add("hn_jet_3prong_Sxy_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxy}, {axisSigmaLxy}, {axisSxy}, {axisMultScaledFT0M}}});
408-
if (fillSVxyz) {
409-
registry.add("hn_jet_3prong_Sxyz_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxyz}, {axisSigmaLxyz}, {axisSxyz}, {axisMultScaledFT0M}}});
410-
}
411-
}
412-
registry.add("hn_jet_3prong_Sxy_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMultScaledFT0M}}});
413-
registry.add("hn_taggedjet_3prong_Sxy_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMultScaledFT0M}}});
414-
if (fillSVxyz) {
415-
registry.add("hn_jet_3prong_Sxyz_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}}});
416-
registry.add("hn_taggedjet_3prong_Sxyz_N1_mult", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}}});
417-
}
418-
}
419396
if (doprocessSV2ProngMCD || doprocessSV2ProngMCDWeighted || doprocessSV2ProngMCPMCDMatched || doprocessSV2ProngMCPMCDMatchedWeighted) {
420397
if (!(doprocessIPsMCD || doprocessIPsMCDWeighted || doprocessIPsMCPMCDMatched || doprocessIPsMCPMCDMatchedWeighted) && !(doprocessJPMCD || doprocessJPMCDWeighted || doprocessJPMCPMCDMatched || doprocessJPMCPMCDMatchedWeighted) && !(doprocessSV3ProngMCD || doprocessSV3ProngMCDWeighted || doprocessSV3ProngMCPMCDMatched || doprocessSV3ProngMCPMCDMatchedWeighted)) {
421398
registry.add("h2_jet_pt_flavour", "", {HistType::kTH2F, {{axisJetPt}, {axisJetFlavour}}});
@@ -464,25 +441,6 @@ struct JetTaggerHFQA {
464441
registry.add("h3_taggedjet_pt_3prong_mass_N1_flavour", "", {HistType::kTH3F, {{axisJetPt}, {axisMass}, {axisJetFlavour}}});
465442
registry.add("h3_taggedjet_pt_3prong_mass_xyz_N1_flavour", "", {HistType::kTH3F, {{axisJetPt}, {axisMass}, {axisJetFlavour}}});
466443
}
467-
if (doprocessSV3ProngMCDMult || doprocessSV3ProngMCDMultWeighted || doprocessSV3ProngMCPMCDMatchedMult || doprocessSV3ProngMCPMCDMatchedMultWeighted) {
468-
registry.add("h_event_mult", "", {HistType::kTH1F, {{axisMultScaledFT0M}}});
469-
registry.add("h3_jet_pt_mult_flavour", "", {HistType::kTH3F, {{axisJetPt}, {axisMultScaledFT0M}, {axisJetFlavour}}});
470-
registry.add("h3_jet_eta_mult_flavour", "", {HistType::kTH3F, {{axisEta}, {axisMultScaledFT0M}, {axisJetFlavour}}});
471-
registry.add("h3_jet_phi_mult_flavour", "", {HistType::kTH3F, {{axisPhi}, {axisMultScaledFT0M}, {axisJetFlavour}}});
472-
if (fillGeneralSVQA) {
473-
registry.add("h3_3prong_nprongs_mult_flavour", "", {HistType::kTH3F, {{axisNprongs}, {axisMultScaledFT0M}, {axisJetFlavour}}});
474-
registry.add("hn_jet_3prong_Sxy_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxy}, {axisSigmaLxy}, {axisSxy}, {axisMultScaledFT0M}, {axisJetFlavour}}});
475-
if (fillSVxyz) {
476-
registry.add("hn_jet_3prong_Sxyz_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisLxyz}, {axisSigmaLxyz}, {axisSxyz}, {axisMultScaledFT0M}, {axisJetFlavour}}});
477-
}
478-
}
479-
registry.add("hn_jet_3prong_Sxy_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
480-
registry.add("hn_taggedjet_3prong_Sxy_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
481-
if (fillSVxyz) {
482-
registry.add("hn_jet_3prong_Sxyz_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
483-
registry.add("hn_taggedjet_3prong_Sxyz_N1_mult_flavour", "", {HistType::kTHnSparseF, {{axisJetPt}, {axisSxy}, {axisSxyz}, {axisMass}, {axisMultScaledFT0M}, {axisJetFlavour}}});
484-
}
485-
}
486444
}
487445

488446
// Filter trackCuts = (aod::jtrack::pt >= trackCuts->at(0) && aod::jtrack::pt < trackCuts->at(1) && aod::jtrack::eta > trackCuts->at(2) && aod::jtrack::eta < trackCuts->at(3));
@@ -546,21 +504,6 @@ struct JetTaggerHFQA {
546504
return true;
547505
}
548506

549-
float getScaledFT0A(const float multFT0A)
550-
{
551-
return multFT0A / meanFT0A;
552-
}
553-
554-
float getScaledFT0C(const float multFT0C)
555-
{
556-
return multFT0C / meanFT0C;
557-
}
558-
559-
float getScaledFT0M(const float multFT0A, const float multFT0C)
560-
{
561-
return 0.5 * (getScaledFT0A(multFT0A) + getScaledFT0C(multFT0C));
562-
}
563-
564507
template <typename U, typename T, typename V, typename W, typename X>
565508
void fillValidationFlavourDefMCD(T const& mcdjet, V const& tracks, W const& particles, X const& particlesPerColl, float eventWeight = 1.0)
566509
{
@@ -1052,50 +995,6 @@ struct JetTaggerHFQA {
1052995
}
1053996
}
1054997

1055-
template <typename T, typename U, typename V>
1056-
void fillHistogramSV3ProngDataMult(T const& collision, U const& jet, V const& /*prongs*/)
1057-
{
1058-
if (jet.template secondaryVertices_as<V>().size() < 1)
1059-
return;
1060-
float multFT0A = collision.multFT0A();
1061-
float multFT0C = collision.multFT0C();
1062-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1063-
registry.fill(HIST("h2_jet_pt_mult"), jet.pt(), scaledFT0M);
1064-
registry.fill(HIST("h2_jet_eta_mult"), jet.eta(), scaledFT0M);
1065-
registry.fill(HIST("h2_jet_phi_mult"), jet.phi(), scaledFT0M);
1066-
if (fillGeneralSVQA) {
1067-
registry.fill(HIST("h2_3prong_nprongs_mult"), jet.template secondaryVertices_as<V>().size(), scaledFT0M);
1068-
for (const auto& prong : jet.template secondaryVertices_as<V>()) {
1069-
registry.fill(HIST("hn_jet_3prong_Sxy_mult"), jet.pt(), prong.decayLengthXY(), prong.errorDecayLengthXY(), prong.decayLengthXY() / prong.errorDecayLengthXY(), scaledFT0M);
1070-
if (fillSVxyz) {
1071-
registry.fill(HIST("hn_jet_3prong_Sxyz_mult"), jet.pt(), prong.decayLength(), prong.errorDecayLength(), prong.decayLength() / prong.errorDecayLength(), scaledFT0M);
1072-
}
1073-
}
1074-
}
1075-
bool checkSv = false;
1076-
auto bjetCand = jettaggingutilities::jetFromProngMaxDecayLength<V>(jet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(2), prongCuts->at(4), prongCuts->at(5), false, &checkSv);
1077-
if (checkSv && jettaggingutilities::svAcceptance(bjetCand, svDispersionMax)) {
1078-
auto maxSxy = bjetCand.decayLengthXY() / bjetCand.errorDecayLengthXY();
1079-
auto massSV = bjetCand.m();
1080-
registry.fill(HIST("hn_jet_3prong_Sxy_N1_mult"), jet.pt(), maxSxy, massSV, scaledFT0M);
1081-
if (jet.isTagged(BJetTaggingMethod::SV)) {
1082-
registry.fill(HIST("hn_taggedjet_3prong_Sxy_N1_mult"), jet.pt(), maxSxy, massSV, scaledFT0M);
1083-
}
1084-
}
1085-
if (fillSVxyz) {
1086-
checkSv = false;
1087-
auto bjetCandXYZ = jettaggingutilities::jetFromProngMaxDecayLength<V>(jet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(3), prongCuts->at(4), prongCuts->at(5), true, &checkSv);
1088-
if (checkSv && jettaggingutilities::svAcceptance(bjetCandXYZ, svDispersionMax)) {
1089-
auto maxSxyz = bjetCandXYZ.decayLength() / bjetCandXYZ.errorDecayLength();
1090-
auto massSV = bjetCandXYZ.m();
1091-
registry.fill(HIST("hn_jet_3prong_Sxyz_N1_mult"), jet.pt(), maxSxyz, massSV, scaledFT0M);
1092-
if (jet.isTagged(BJetTaggingMethod::SV3D)) {
1093-
registry.fill(HIST("hn_taggedjet_3prong_Sxyz_N1_mult"), jet.pt(), maxSxyz, massSV, scaledFT0M);
1094-
}
1095-
}
1096-
}
1097-
}
1098-
1099998
template <typename T, typename U>
1100999
void fillHistogramSV2ProngMCD(T const& mcdjet, U const& /*prongs*/, float eventWeight = 1.0)
11011000
{
@@ -1208,54 +1107,6 @@ struct JetTaggerHFQA {
12081107
}
12091108
}
12101109

1211-
template <typename T, typename U, typename V>
1212-
void fillHistogramSV3ProngMCDMult(T const& collision, U const& mcdjet, V const& /*prongs*/, float eventWeight = 1.0)
1213-
{
1214-
float pTHat = 10. / (std::pow(eventWeight, 1.0 / pTHatExponent));
1215-
if (mcdjet.pt() > pTHatMaxMCD * pTHat) {
1216-
return;
1217-
}
1218-
auto origin = mcdjet.origin();
1219-
float multFT0A = collision.multFT0A();
1220-
float multFT0C = collision.multFT0C();
1221-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1222-
if (mcdjet.template secondaryVertices_as<V>().size() < 1)
1223-
return;
1224-
registry.fill(HIST("h3_jet_pt_mult_flavour"), mcdjet.pt(), scaledFT0M, origin, eventWeight);
1225-
registry.fill(HIST("h3_jet_eta_mult_flavour"), mcdjet.eta(), scaledFT0M, origin, eventWeight);
1226-
registry.fill(HIST("h3_jet_phi_mult_flavour"), mcdjet.phi(), scaledFT0M, origin, eventWeight);
1227-
if (fillGeneralSVQA) {
1228-
registry.fill(HIST("h3_3prong_nprongs_mult_flavour"), mcdjet.template secondaryVertices_as<V>().size(), scaledFT0M, origin, eventWeight);
1229-
for (const auto& prong : mcdjet.template secondaryVertices_as<V>()) {
1230-
registry.fill(HIST("hn_jet_3prong_Sxy_mult_flavour"), mcdjet.pt(), prong.decayLengthXY(), prong.errorDecayLengthXY(), prong.decayLengthXY() / prong.errorDecayLengthXY(), scaledFT0M, origin, eventWeight);
1231-
if (fillSVxyz) {
1232-
registry.fill(HIST("hn_jet_3prong_Sxyz_mult_flavour"), mcdjet.pt(), prong.decayLength(), prong.errorDecayLength(), prong.decayLength() / prong.errorDecayLength(), scaledFT0M, origin, eventWeight);
1233-
}
1234-
}
1235-
}
1236-
bool checkSv = false;
1237-
auto bjetCand = jettaggingutilities::jetFromProngMaxDecayLength<V>(mcdjet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(2), prongCuts->at(4), prongCuts->at(5), false, &checkSv);
1238-
if (checkSv && jettaggingutilities::svAcceptance(bjetCand, svDispersionMax)) {
1239-
auto maxSxy = bjetCand.decayLengthXY() / bjetCand.errorDecayLengthXY();
1240-
auto massSV = bjetCand.m();
1241-
registry.fill(HIST("hn_jet_3prong_Sxy_N1_mult_flavour"), mcdjet.pt(), maxSxy, massSV, scaledFT0M, origin, eventWeight);
1242-
if (mcdjet.isTagged(BJetTaggingMethod::SV)) {
1243-
registry.fill(HIST("hn_taggedjet_3prong_Sxy_N1_mult_flavour"), mcdjet.pt(), maxSxy, massSV, scaledFT0M, origin, eventWeight);
1244-
}
1245-
}
1246-
if (fillSVxyz) {
1247-
auto bjetCandXYZ = jettaggingutilities::jetFromProngMaxDecayLength<V>(mcdjet, prongCuts->at(0), prongCuts->at(1), prongCuts->at(3), prongCuts->at(4), prongCuts->at(5), true, &checkSv);
1248-
if (checkSv && jettaggingutilities::svAcceptance(bjetCandXYZ, svDispersionMax)) {
1249-
auto maxSxyz = bjetCandXYZ.decayLength() / bjetCandXYZ.errorDecayLength();
1250-
auto massSV = bjetCandXYZ.m();
1251-
registry.fill(HIST("hn_jet_3prong_Sxyz_N1_mult_flavour"), mcdjet.pt(), maxSxyz, massSV, scaledFT0M, origin, eventWeight);
1252-
if (mcdjet.isTagged(BJetTaggingMethod::SV3D)) {
1253-
registry.fill(HIST("hn_taggedjet_3prong_Sxyz_N1_mult_flavour"), mcdjet.pt(), maxSxyz, massSV, scaledFT0M, origin, eventWeight);
1254-
}
1255-
}
1256-
}
1257-
}
1258-
12591110
void processDummy(aod::Collision const&, aod::Tracks const&)
12601111
{
12611112
}
@@ -1749,27 +1600,6 @@ struct JetTaggerHFQA {
17491600
}
17501601
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngData, "Fill 3prong imformation for data jets", false);
17511602

1752-
void processSV3ProngDataMult(soa::Filtered<aod::JetCollisions>::iterator const& collision, soa::Join<JetTableData, TagTableData, aod::DataSecondaryVertex3ProngIndices> const& jets, aod::DataSecondaryVertex3Prongs const& prongs)
1753-
{
1754-
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
1755-
return;
1756-
}
1757-
float multFT0A = collision.multFT0A();
1758-
float multFT0C = collision.multFT0C();
1759-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1760-
registry.fill(HIST("h_event_mult"), scaledFT0M);
1761-
for (auto const& jet : jets) {
1762-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
1763-
continue;
1764-
}
1765-
if (!isAcceptedJet<aod::JetTracks>(jet)) {
1766-
continue;
1767-
}
1768-
fillHistogramSV3ProngDataMult(collision, jet, prongs);
1769-
}
1770-
}
1771-
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngDataMult, "Fill 3prong imformation for data jets", false);
1772-
17731603
void processSV2ProngMCD(soa::Filtered<aod::JetCollisions>::iterator const& collision, soa::Join<JetTableMCD, TagTableMCD, aod::MCDSecondaryVertex2ProngIndices> const& mcdjets, aod::MCDSecondaryVertex2Prongs const& prongs)
17741604
{
17751605
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
@@ -1941,104 +1771,6 @@ struct JetTaggerHFQA {
19411771
}
19421772
}
19431773
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngMCPMCDMatchedWeighted, "Fill 3prong imformation for mcd jets", false);
1944-
1945-
void processSV3ProngMCDMult(soa::Filtered<aod::JetCollisions>::iterator const& collision, soa::Join<JetTableMCD, TagTableMCD, aod::MCDSecondaryVertex3ProngIndices> const& mcdjets, aod::MCDSecondaryVertex3Prongs const& prongs)
1946-
{
1947-
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
1948-
return;
1949-
}
1950-
float multFT0A = collision.multFT0A();
1951-
float multFT0C = collision.multFT0C();
1952-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1953-
registry.fill(HIST("h_event_mult"), scaledFT0M);
1954-
for (auto const& mcdjet : mcdjets) {
1955-
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
1956-
continue;
1957-
}
1958-
if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
1959-
continue;
1960-
}
1961-
fillHistogramSV3ProngMCDMult(collision, mcdjet, prongs);
1962-
}
1963-
}
1964-
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngMCDMult, "Fill 3prong imformation for mcd jets with multiplicity", false);
1965-
1966-
void processSV3ProngMCDMultWeighted(soa::Filtered<soa::Join<aod::JetCollisionsMCD, aod::JCollisionOutliers>>::iterator const& collision, soa::Join<JetTableMCD, TagTableMCD, aod::MCDSecondaryVertex3ProngIndices> const& mcdjets, aod::MCDSecondaryVertex3Prongs const& prongs)
1967-
{
1968-
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
1969-
return;
1970-
}
1971-
if (collision.isOutlier()) {
1972-
return;
1973-
}
1974-
float multFT0A = collision.multFT0A();
1975-
float multFT0C = collision.multFT0C();
1976-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1977-
registry.fill(HIST("h_event_mult"), scaledFT0M);
1978-
for (auto const& mcdjet : mcdjets) {
1979-
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
1980-
continue;
1981-
}
1982-
if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
1983-
continue;
1984-
}
1985-
fillHistogramSV3ProngMCDMult(collision, mcdjet, prongs, collision.weight());
1986-
}
1987-
}
1988-
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngMCDMultWeighted, "Fill 3prong imformation for mcd jets with multiplicity weighted", false);
1989-
1990-
void processSV3ProngMCPMCDMatchedMult(soa::Filtered<soa::Join<aod::JetCollisions, aod::JCollisionPIs, aod::JMcCollisionLbs>>::iterator const& collision, soa::Join<JetTableMCD, TagTableMCD, JetTableMCDMCP, aod::MCDSecondaryVertex3ProngIndices> const& mcdjets, soa::Join<JetTableMCP, JetTableMCPMCD> const& /*mcpjets*/, aod::MCDSecondaryVertex3Prongs const& prongs)
1991-
{
1992-
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
1993-
return;
1994-
}
1995-
float multFT0A = collision.multFT0A();
1996-
float multFT0C = collision.multFT0C();
1997-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
1998-
registry.fill(HIST("h_event_mult"), scaledFT0M);
1999-
for (auto const& mcdjet : mcdjets) {
2000-
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
2001-
continue;
2002-
}
2003-
if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
2004-
continue;
2005-
}
2006-
if (!mcdjet.has_matchedJetGeo()) {
2007-
continue;
2008-
}
2009-
if (!doprocessSV3ProngMCDMult)
2010-
fillHistogramSV3ProngMCDMult(collision, mcdjet, prongs);
2011-
}
2012-
}
2013-
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngMCPMCDMatchedMult, "Fill 3prong imformation for mcd jets matched with multiplicity", false);
2014-
2015-
void processSV3ProngMCPMCDMatchedMultWeighted(soa::Filtered<soa::Join<aod::JetCollisionsMCD, aod::JCollisionPIs, aod::JCollisionOutliers>>::iterator const& collision, soa::Join<JetTableMCD, TagTableMCD, JetTableMCDMCP, aod::MCDSecondaryVertex3ProngIndices> const& mcdjets, soa::Join<JetTableMCP, JetTableMCPMCD> const& /*mcpjets*/, aod::MCDSecondaryVertex3Prongs const& prongs)
2016-
{
2017-
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
2018-
return;
2019-
}
2020-
if (collision.isOutlier()) {
2021-
return;
2022-
}
2023-
float multFT0A = collision.multFT0A();
2024-
float multFT0C = collision.multFT0C();
2025-
float scaledFT0M = getScaledFT0M(multFT0A, multFT0C);
2026-
registry.fill(HIST("h_event_mult"), scaledFT0M);
2027-
for (auto const& mcdjet : mcdjets) {
2028-
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaCuts->at(0), jetEtaCuts->at(1), trackCuts->at(2), trackCuts->at(3))) {
2029-
continue;
2030-
}
2031-
if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
2032-
continue;
2033-
}
2034-
if (!mcdjet.has_matchedJetGeo()) {
2035-
continue;
2036-
}
2037-
if (!doprocessSV3ProngMCDMultWeighted)
2038-
fillHistogramSV3ProngMCDMult(collision, mcdjet, prongs, collision.weight());
2039-
}
2040-
}
2041-
PROCESS_SWITCH(JetTaggerHFQA, processSV3ProngMCPMCDMatchedMultWeighted, "Fill 3prong imformation for mcd jets matched with multiplicity weightd", false);
20421774
};
20431775

20441776
using JetTaggerQAChargedDataJets = soa::Join<aod::ChargedJets, aod::ChargedJetConstituents>;

0 commit comments

Comments
 (0)