Skip to content

Commit 4713a7e

Browse files
zjxiongOvOzjxiongOvO
authored
[PWGDQ] Updated the definition of polarization variables (#11375)
1 parent eb026ea commit 4713a7e

File tree

3 files changed

+320
-102
lines changed

3 files changed

+320
-102
lines changed

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 28 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -972,14 +972,34 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
972972
hm->AddHistogram(histClass, "VtxZ_VtxNcontribReal", "VtxZ vs VtxNcontribReal", false, 240, -12.0, 12.0, VarManager::kVtxZ, 200, 0, 200.0, VarManager::kVtxNcontribReal);
973973
}
974974
if (subGroupStr.Contains("polarization")) {
975-
hm->AddHistogram(histClass, "cosThetaHE", "", false, 100, -1., 1., VarManager::kCosThetaHE);
976-
hm->AddHistogram(histClass, "cosThetaCS", "", false, 100, -1., 1., VarManager::kCosThetaCS);
977-
hm->AddHistogram(histClass, "PhiHE", "", false, 100, -o2::constants::math::PI, o2::constants::math::PI, VarManager::kPhiHE);
978-
hm->AddHistogram(histClass, "PhiCS", "", false, 100, -o2::constants::math::PI, o2::constants::math::PI, VarManager::kPhiCS);
979-
hm->AddHistogram(histClass, "Mass_Pt_cosThetaHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 250, 0.0, 25.0, VarManager::kPt, 40, -1., 1., VarManager::kCosThetaHE);
980-
hm->AddHistogram(histClass, "Mass_Pt_cosThetaCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 250, 0.0, 25.0, VarManager::kPt, 40, -1., 1., VarManager::kCosThetaCS);
981-
hm->AddHistogram(histClass, "Mass_Pt_PhiHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 250, 0.0, 25.0, VarManager::kPt, 40, -o2::constants::math::PI, o2::constants::math::PI, VarManager::kPhiHE);
982-
hm->AddHistogram(histClass, "Mass_Pt_PhiCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 250, 0.0, 25.0, VarManager::kPt, 40, -o2::constants::math::PI, o2::constants::math::PI, VarManager::kPhiCS);
975+
if (subGroupStr.Contains("helicity")) {
976+
hm->AddHistogram(histClass, "cosThetaHE", "", false, 100, -1., 1., VarManager::kCosThetaHE);
977+
hm->AddHistogram(histClass, "phiHE", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiHE);
978+
hm->AddHistogram(histClass, "phitildeHE", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiTildeHE);
979+
hm->AddHistogram(histClass, "Mass_Pt_CosThetaHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaHE);
980+
hm->AddHistogram(histClass, "Mass_Pt_PhiHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiHE);
981+
hm->AddHistogram(histClass, "Mass_Pt_PhiTildeHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiTildeHE);
982+
}
983+
if (subGroupStr.Contains("collins-soper")) {
984+
hm->AddHistogram(histClass, "cosThetaCS", "", false, 100, -1., 1., VarManager::kCosThetaCS);
985+
hm->AddHistogram(histClass, "phiCS", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiCS);
986+
hm->AddHistogram(histClass, "phitildeCS", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiTildeCS);
987+
hm->AddHistogram(histClass, "Mass_Pt_CosThetaCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaCS);
988+
hm->AddHistogram(histClass, "Mass_Pt_PhiCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiCS);
989+
hm->AddHistogram(histClass, "Mass_Pt_PhiTildeCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiTildeCS);
990+
}
991+
if (subGroupStr.Contains("production")) {
992+
hm->AddHistogram(histClass, "cosThetaPP", "", false, 100, -1., 1., VarManager::kCosThetaPP);
993+
hm->AddHistogram(histClass, "phiPP", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiPP);
994+
hm->AddHistogram(histClass, "phitildePP", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiTildePP);
995+
hm->AddHistogram(histClass, "Mass_Pt_CosThetaPP", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaPP);
996+
hm->AddHistogram(histClass, "Mass_Pt_PhiPP", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiPP);
997+
hm->AddHistogram(histClass, "Mass_Pt_PhiTildePP", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiTildePP);
998+
}
999+
if (subGroupStr.Contains("random")) {
1000+
hm->AddHistogram(histClass, "cosThetaRM", "", false, 100, -1., 1., VarManager::kCosThetaRM);
1001+
hm->AddHistogram(histClass, "Mass_Pt_CosThetaRM", "", false, 200, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaRM);
1002+
}
9831003
}
9841004
if (subGroupStr.Contains("upsilon")) {
9851005
hm->AddHistogram(histClass, "MassUpsilon_Pt", "", false, 500, 7.0, 12.0, VarManager::kMass, 400, 0.0, 40.0, VarManager::kPt);

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 38 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -542,18 +542,26 @@ void VarManager::SetDefaultVarNames()
542542
fgVariableUnits[kMCVy] = "cm"; // TODO: check the unit
543543
fgVariableNames[kMCVz] = "MC vz";
544544
fgVariableUnits[kMCVz] = "cm"; // TODO: check the unit
545-
fgVariableNames[kMCCosThetaHE] = "cos#it{#theta}^{MC}_{HE}";
545+
fgVariableNames[kMCCosThetaHE] = "MC cos(#theta_{HE})";
546546
fgVariableUnits[kMCCosThetaHE] = "";
547-
fgVariableNames[kMCPhiHE] = "#varphi^{MC}_{HE}";
548-
fgVariableUnits[kMCPhiHE] = "rad.";
549-
fgVariableNames[kMCCosThetaCS] = "cos#it{#theta}^{MC}_{CS}";
547+
fgVariableNames[kMCPhiHE] = "MC #varphi_{HE}";
548+
fgVariableUnits[kMCPhiHE] = "rad";
549+
fgVariableNames[kMCPhiTildeHE] = "MC #tilde{#varphi}_{HE}";
550+
fgVariableUnits[kMCPhiTildeHE] = "rad";
551+
fgVariableNames[kMCCosThetaCS] = "MC cos(#theta_{CS})";
550552
fgVariableUnits[kMCCosThetaCS] = "";
551-
fgVariableNames[kMCPhiCS] = "#varphi^{MC}_{CS}";
552-
fgVariableUnits[kMCPhiCS] = "rad.";
553-
fgVariableNames[kMCCosThetaPP] = "cos#it{#theta}^{MC}_{PP}";
553+
fgVariableNames[kMCPhiCS] = "MC #varphi_{CS}";
554+
fgVariableUnits[kMCPhiCS] = "rad";
555+
fgVariableNames[kMCPhiTildeCS] = "MC #tilde{#varphi}_{CS}";
556+
fgVariableUnits[kMCPhiTildeCS] = "rad";
557+
fgVariableNames[kMCCosThetaPP] = "MC cos(#theta_{PP})";
554558
fgVariableUnits[kMCCosThetaPP] = "";
555-
fgVariableNames[kMCPhiPP] = "#varphi^{MC}_{PP}";
556-
fgVariableUnits[kMCPhiPP] = "rad.";
559+
fgVariableNames[kMCPhiPP] = "MC #varphi_{PP}";
560+
fgVariableUnits[kMCPhiPP] = "rad";
561+
fgVariableNames[kMCPhiTildePP] = "MC #tilde{#varphi}_{PP}";
562+
fgVariableUnits[kMCPhiTildePP] = "rad";
563+
fgVariableNames[kMCCosThetaRM] = "MC cos(#theta_{RM})";
564+
fgVariableUnits[kMCCosThetaRM] = "";
557565
fgVariableNames[kCandidateId] = "";
558566
fgVariableUnits[kCandidateId] = "";
559567
fgVariableNames[kPairType] = "Pair type";
@@ -934,18 +942,26 @@ void VarManager::SetDefaultVarNames()
934942
fgVariableUnits[kDeltaPhi] = "rad.";
935943
fgVariableNames[kDeltaPhiSym] = "#Delta#phi";
936944
fgVariableUnits[kDeltaPhiSym] = "rad.";
937-
fgVariableNames[kCosThetaHE] = "cos#it{#theta}_{HE}";
945+
fgVariableNames[kCosThetaHE] = "cos#it{#theta}";
938946
fgVariableUnits[kCosThetaHE] = "";
939947
fgVariableNames[kPhiHE] = "#varphi_{HE}";
940948
fgVariableUnits[kPhiHE] = "rad.";
949+
fgVariableNames[kPhiTildeHE] = "#tilde{#varphi}_{HE}";
950+
fgVariableUnits[kPhiTildeHE] = "rad.";
941951
fgVariableNames[kCosThetaCS] = "cos#it{#theta}_{CS}";
942952
fgVariableUnits[kCosThetaCS] = "";
943953
fgVariableNames[kPhiCS] = "#varphi_{CS}";
944954
fgVariableUnits[kPhiCS] = "rad.";
955+
fgVariableNames[kPhiTildeCS] = "#tilde{#varphi}_{CS}";
956+
fgVariableUnits[kPhiTildeCS] = "rad.";
945957
fgVariableNames[kCosThetaPP] = "cos#it{#theta}_{PP}";
946958
fgVariableUnits[kCosThetaPP] = "";
947959
fgVariableNames[kPhiPP] = "#varphi_{PP}";
948960
fgVariableUnits[kPhiPP] = "rad.";
961+
fgVariableNames[kPhiTildePP] = "#tilde{#varphi}_{PP}";
962+
fgVariableUnits[kPhiTildePP] = "rad.";
963+
fgVariableNames[kCosThetaRM] = "cos#it{#theta}_{RM}";
964+
fgVariableUnits[kCosThetaRM] = "";
949965
fgVariableNames[kCosPhiVP] = "cos#it{#varphi}_{VP}";
950966
fgVariableUnits[kCosPhiVP] = "";
951967
fgVariableNames[kPhiVP] = "#varphi_{VP} - #Psi_{2}";
@@ -1462,11 +1478,15 @@ void VarManager::SetDefaultVarNames()
14621478
fgVarNamesMap["kMCEta"] = kMCEta;
14631479
fgVarNamesMap["kMCY"] = kMCY;
14641480
fgVarNamesMap["kMCCosThetaHE"] = kMCCosThetaHE;
1465-
fgVarNamesMap["kMCCosThetaCS"] = kMCCosThetaCS;
1466-
fgVarNamesMap["kMCCosThetaPP"] = kMCCosThetaPP;
14671481
fgVarNamesMap["kMCPhiHE"] = kMCPhiHE;
1482+
fgVarNamesMap["kMCPhiTildeHE"] = kMCPhiTildeHE;
1483+
fgVarNamesMap["kMCCosThetaCS"] = kMCCosThetaCS;
14681484
fgVarNamesMap["kMCPhiCS"] = kMCPhiCS;
1485+
fgVarNamesMap["kMCPhiTildeCS"] = kMCPhiTildeCS;
1486+
fgVarNamesMap["kMCCosThetaPP"] = kMCCosThetaPP;
14691487
fgVarNamesMap["kMCPhiPP"] = kMCPhiPP;
1488+
fgVarNamesMap["kMCPhiTildePP"] = kMCPhiTildePP;
1489+
fgVarNamesMap["kMCCosThetaRM"] = kMCCosThetaRM;
14701490
fgVarNamesMap["kMCParticleGeneratorId"] = kMCParticleGeneratorId;
14711491
fgVarNamesMap["kNMCParticleVariables"] = kNMCParticleVariables;
14721492
fgVarNamesMap["kMCMotherPdgCode"] = kMCMotherPdgCode;
@@ -1496,11 +1516,15 @@ void VarManager::SetDefaultVarNames()
14961516
fgVarNamesMap["kVertexingProcCode"] = kVertexingProcCode;
14971517
fgVarNamesMap["kVertexingChi2PCA"] = kVertexingChi2PCA;
14981518
fgVarNamesMap["kCosThetaHE"] = kCosThetaHE;
1499-
fgVarNamesMap["kCosThetaCS"] = kCosThetaCS;
1500-
fgVarNamesMap["kCosThetaPP"] = kCosThetaPP;
15011519
fgVarNamesMap["kPhiHE"] = kPhiHE;
1520+
fgVarNamesMap["kPhiTildeHE"] = kPhiTildeHE;
1521+
fgVarNamesMap["kCosThetaCS"] = kCosThetaCS;
15021522
fgVarNamesMap["kPhiCS"] = kPhiCS;
1523+
fgVarNamesMap["kPhiTildeCS"] = kPhiTildeCS;
1524+
fgVarNamesMap["kCosThetaPP"] = kCosThetaPP;
15031525
fgVarNamesMap["kPhiPP"] = kPhiPP;
1526+
fgVarNamesMap["kPhiTildePP"] = kPhiTildePP;
1527+
fgVarNamesMap["kCosThetaRM"] = kCosThetaRM;
15041528
fgVarNamesMap["kCosPhiVP"] = kCosPhiVP;
15051529
fgVarNamesMap["kPhiVP"] = kPhiVP;
15061530
fgVarNamesMap["kDeltaPhiPair2"] = kDeltaPhiPair2;

0 commit comments

Comments
 (0)