Skip to content

Commit e473593

Browse files
dmallick2dmallick2
authored
[PWGDQ] Add GJ axis definition and fill histograms (#18271)
1 parent be2cc33 commit e473593

4 files changed

Lines changed: 390 additions & 17 deletions

File tree

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 32 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1306,14 +1306,28 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
13061306
std::array<int, 5> bins = {50, 30, 10, 20, 10};
13071307
std::array<double, 5> xmin = {2., 0., 0, -1., 0.0};
13081308
std::array<double, 5> xmax = {4., 3., 100, 1., 6.28};
1309-
hm->AddHistogram(histClass, "MC_Dielectron_Mass_Pt_Cent_cosThetaHE", "", 5, varspTHE.data(), bins.data(), xmin.data(), xmax.data(), 0, -1, kFALSE);
1309+
hm->AddHistogram(histClass, "MC_Dielectron_Mass_Pt_Cent_cosThetaHE", "", 5, varspTHE.data(), bins.data(), xmin.data(), xmax.data(), nullptr, -1, kFALSE);
1310+
}
1311+
if (groupStr.CompareTo("polarization-dielectron-pbpb-midy-gjb1-gen") == 0) {
1312+
std::array<int, 5> varspTGJB1 = {VarManager::kMCMass, VarManager::kMCPt, VarManager::kMCEventCentrFT0C, VarManager::kMCCosThetaGJB1, VarManager::kMCPhiGJB1};
1313+
std::array<int, 5> bins = {50, 30, 10, 20, 10};
1314+
std::array<double, 5> xmin = {2., 0., 0, -1., 0.0};
1315+
std::array<double, 5> xmax = {4., 3., 100, 1., 6.28};
1316+
hm->AddHistogram(histClass, "MC_Dielectron_Mass_Pt_Cent_cosThetaGJB1", "", 5, varspTGJB1.data(), bins.data(), xmin.data(), xmax.data(), nullptr, -1, kFALSE);
1317+
}
1318+
if (groupStr.CompareTo("polarization-dielectron-pbpb-midy-gjb2-gen") == 0) {
1319+
std::array<int, 5> varspTGJB2 = {VarManager::kMCMass, VarManager::kMCPt, VarManager::kMCEventCentrFT0C, VarManager::kMCCosThetaGJB2, VarManager::kMCPhiGJB2};
1320+
std::array<int, 5> bins = {50, 30, 10, 20, 10};
1321+
std::array<double, 5> xmin = {2., 0., 0, -1., 0.0};
1322+
std::array<double, 5> xmax = {4., 3., 100, 1., 6.28};
1323+
hm->AddHistogram(histClass, "MC_Dielectron_Mass_Pt_Cent_cosThetaGJB2", "", 5, varspTGJB2.data(), bins.data(), xmin.data(), xmax.data(), nullptr, -1, kFALSE);
13101324
}
13111325
if (groupStr.CompareTo("polarization-dielectron-pbpb-midy-cs-gen") == 0) {
13121326
std::array<int, 5> varspTCS = {VarManager::kMCMass, VarManager::kMCPt, VarManager::kMCEventCentrFT0C, VarManager::kMCCosThetaCS, VarManager::kMCPhiCS};
13131327
std::array<int, 5> bins = {50, 30, 10, 20, 10};
13141328
std::array<double, 5> xmin = {2., 0., 0, -1., 0.0};
13151329
std::array<double, 5> xmax = {4., 3., 100, 1., 6.28};
1316-
hm->AddHistogram(histClass, "MC_Dielectron_Mass_Pt_Cent_cosThetaCS", "", 5, varspTCS.data(), bins.data(), xmin.data(), xmax.data(), 0, -1, kFALSE);
1330+
hm->AddHistogram(histClass, "MC_Dielectron_Mass_Pt_Cent_cosThetaCS", "", 5, varspTCS.data(), bins.data(), xmin.data(), xmax.data(), nullptr, -1, kFALSE);
13171331
}
13181332
if (groupStr.CompareTo("pair") == 0) {
13191333
if (subGroupStr.Contains("cepf")) {
@@ -1422,6 +1436,22 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
14221436
hm->AddHistogram(histClass, "Dielectron_Mass_Pt_Cent_cosThetaHE", "", 5, varsHEpbpb.data(), binspT.data(), xminpT.data(), xmaxpT.data(), nullptr, -1, kFALSE);
14231437
}
14241438

1439+
if (subGroupStr.Contains("dielectron-polarization-gjb1-pbpb")) {
1440+
std::array<int, 5> varsGJB1pbpb = {VarManager::kMass, VarManager::kPt, VarManager::kCentFT0C, VarManager::kCosThetaGJB1, VarManager::kPhiGJB1};
1441+
std::array<int, 5> binspT = {100, 30, 10, 10, 10};
1442+
std::array<double, 5> xminpT = {2., 0., 0, -1., 0.0};
1443+
std::array<double, 5> xmaxpT = {4.5, 3., 100, 1., 6.28};
1444+
hm->AddHistogram(histClass, "Dielectron_Mass_Pt_Cent_cosThetaGJB1", "", 5, varsGJB1pbpb.data(), binspT.data(), xminpT.data(), xmaxpT.data(), nullptr, -1, kFALSE);
1445+
}
1446+
1447+
if (subGroupStr.Contains("dielectron-polarization-gjb2-pbpb")) {
1448+
std::array<int, 5> varsGJB2pbpb = {VarManager::kMass, VarManager::kPt, VarManager::kCentFT0C, VarManager::kCosThetaGJB2, VarManager::kPhiGJB2};
1449+
std::array<int, 5> binspT = {100, 30, 10, 10, 10};
1450+
std::array<double, 5> xminpT = {2., 0., 0, -1., 0.0};
1451+
std::array<double, 5> xmaxpT = {4.5, 3., 100, 1., 6.28};
1452+
hm->AddHistogram(histClass, "Dielectron_Mass_Pt_Cent_cosThetaGJB2", "", 5, varsGJB2pbpb.data(), binspT.data(), xminpT.data(), xmaxpT.data(), nullptr, -1, kFALSE);
1453+
}
1454+
14251455
if (subGroupStr.Contains("dielectron-polarization-cs-pbpb")) {
14261456
std::array<int, 5> varsCSpbpb = {VarManager::kMass, VarManager::kPt, VarManager::kCentFT0C, VarManager::kCosThetaCS, VarManager::kPhiCS};
14271457
std::array<int, 5> binspT = {100, 30, 10, 10, 10};

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 36 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1273,6 +1273,18 @@ void VarManager::SetDefaultVarNames()
12731273
fgVariableUnits[kMCPhiCS] = "rad";
12741274
fgVariableNames[kMCPhiTildeCS] = "MC #tilde{#varphi}_{CS}";
12751275
fgVariableUnits[kMCPhiTildeCS] = "rad";
1276+
fgVariableNames[kMCCosThetaGJB1] = "MC cos#it{#theta}_{GJ1}";
1277+
fgVariableUnits[kMCCosThetaGJB1] = "";
1278+
fgVariableNames[kMCPhiGJB1] = "MC #varphi_{GJ1}";
1279+
fgVariableUnits[kMCPhiGJB1] = "rad";
1280+
fgVariableNames[kMCPhiTildeGJB1] = "MC #tilde{#varphi}_{GJ1}";
1281+
fgVariableUnits[kMCPhiTildeGJB1] = "rad";
1282+
fgVariableNames[kMCCosThetaGJB2] = "MC cos#it{#theta}_{GJ2}";
1283+
fgVariableUnits[kMCCosThetaGJB2] = "";
1284+
fgVariableNames[kMCPhiGJB2] = "MC #varphi_{GJ2}";
1285+
fgVariableUnits[kMCPhiGJB2] = "rad";
1286+
fgVariableNames[kMCPhiTildeGJB2] = "MC #tilde{#varphi}_{GJ2}";
1287+
fgVariableUnits[kMCPhiTildeGJB2] = "rad";
12761288
fgVariableNames[kMCCosThetaPP] = "MC cos(#theta_{PP})";
12771289
fgVariableUnits[kMCCosThetaPP] = "";
12781290
fgVariableNames[kMCPhiPP] = "MC #varphi_{PP}";
@@ -1715,6 +1727,18 @@ void VarManager::SetDefaultVarNames()
17151727
fgVariableUnits[kPhiCS] = "rad.";
17161728
fgVariableNames[kPhiTildeCS] = "#tilde{#varphi}_{CS}";
17171729
fgVariableUnits[kPhiTildeCS] = "rad.";
1730+
fgVariableNames[kCosThetaGJB1] = "cos#it{#theta}_{GJ1}";
1731+
fgVariableUnits[kCosThetaGJB1] = "";
1732+
fgVariableNames[kPhiGJB1] = "#varphi_{GJ1}";
1733+
fgVariableUnits[kPhiGJB1] = "rad.";
1734+
fgVariableNames[kPhiTildeGJB1] = "#tilde{#varphi}_{GJ1}";
1735+
fgVariableUnits[kPhiTildeGJB1] = "rad.";
1736+
fgVariableNames[kCosThetaGJB2] = "cos#it{#theta}_{GJ2}";
1737+
fgVariableUnits[kCosThetaGJB2] = "";
1738+
fgVariableNames[kPhiGJB2] = "#varphi_{GJ2}";
1739+
fgVariableUnits[kPhiGJB2] = "rad.";
1740+
fgVariableNames[kPhiTildeGJB2] = "#tilde{#varphi}_{GJ2}";
1741+
fgVariableUnits[kPhiTildeGJB2] = "rad.";
17181742
fgVariableNames[kCosThetaPP] = "cos#it{#theta}_{PP}";
17191743
fgVariableUnits[kCosThetaPP] = "";
17201744
fgVariableNames[kPhiPP] = "#varphi_{PP}";
@@ -2559,6 +2583,12 @@ void VarManager::SetDefaultVarNames()
25592583
fgVarNamesMap["kMCCosThetaCS"] = kMCCosThetaCS;
25602584
fgVarNamesMap["kMCPhiCS"] = kMCPhiCS;
25612585
fgVarNamesMap["kMCPhiTildeCS"] = kMCPhiTildeCS;
2586+
fgVarNamesMap["kMCCosThetaGJB1"] = kMCCosThetaGJB1;
2587+
fgVarNamesMap["kMCPhiGJB1"] = kMCPhiGJB1;
2588+
fgVarNamesMap["kMCPhiTildeGJB1"] = kMCPhiTildeGJB1;
2589+
fgVarNamesMap["kMCCosThetaGJB2"] = kMCCosThetaGJB2;
2590+
fgVarNamesMap["kMCPhiGJB2"] = kMCPhiGJB2;
2591+
fgVarNamesMap["kMCPhiTildeGJB2"] = kMCPhiTildeGJB2;
25622592
fgVarNamesMap["kMCCosThetaPP"] = kMCCosThetaPP;
25632593
fgVarNamesMap["kMCPhiPP"] = kMCPhiPP;
25642594
fgVarNamesMap["kMCPhiTildePP"] = kMCPhiTildePP;
@@ -2605,6 +2635,12 @@ void VarManager::SetDefaultVarNames()
26052635
fgVarNamesMap["kCosThetaCS"] = kCosThetaCS;
26062636
fgVarNamesMap["kPhiCS"] = kPhiCS;
26072637
fgVarNamesMap["kPhiTildeCS"] = kPhiTildeCS;
2638+
fgVarNamesMap["kCosThetaGJB1"] = kCosThetaGJB1;
2639+
fgVarNamesMap["kPhiGJB1"] = kPhiGJB1;
2640+
fgVarNamesMap["kPhiTildeGJB1"] = kPhiTildeGJB1;
2641+
fgVarNamesMap["kCosThetaGJB2"] = kCosThetaGJB2;
2642+
fgVarNamesMap["kPhiGJB2"] = kPhiGJB2;
2643+
fgVarNamesMap["kPhiTildeGJB2"] = kPhiTildeGJB2;
26082644
fgVarNamesMap["kCosThetaPP"] = kCosThetaPP;
26092645
fgVarNamesMap["kPhiPP"] = kPhiPP;
26102646
fgVarNamesMap["kPhiTildePP"] = kPhiTildePP;

0 commit comments

Comments
 (0)