Skip to content

Commit ac21ceb

Browse files
Elros60Elros60
authored
[PWGDQ] Adding profile histogram of correlations for cumulant evaluation (#7003)
* [PWGDQ] Adding profile histogram of correlations for cumulant evaluation * Moving correlation profiles to new subgroup
1 parent 87e1c6c commit ac21ceb

3 files changed

Lines changed: 45 additions & 13 deletions

File tree

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1009,6 +1009,21 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
10091009
double maxBins[6] = {5.0, 30.0, 4.0, 90.0, 10.0, 2.0};
10101010
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_V2", "", 6, varV2, bins, minBins, maxBins, 0, -1, kTRUE);
10111011
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_V3", "", 6, varV3, bins, minBins, maxBins, 0, -1, kTRUE);
1012+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2REF, VarManager::kM11REFoverMp);
1013+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4REF, VarManager::kM1111REFoverMp);
1014+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2POIMp, VarManager::kM01POIoverMp);
1015+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4POIMp, VarManager::kM0111POIoverMp);
1016+
}
1017+
if (subGroupStr.Contains("cumulant")) {
1018+
int var[4] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C};
1019+
int bins[4] = {250, 60, 6, 18};
1020+
double minBins[4] = {0.0, 0.0, 2.5, 0.0};
1021+
double maxBins[4] = {5.0, 30.0, 4.0, 90.0};
1022+
hm->AddHistogram(histClass, "Mass_Pt_Rapidity_CentFT0C", "", 4, var, bins, minBins, maxBins, 0, -1, kTRUE);
1023+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2REF, VarManager::kM11REFoverMp);
1024+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4REF, VarManager::kM1111REFoverMp);
1025+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2POIMp, VarManager::kM01POIoverMp);
1026+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4POIMp, VarManager::kM0111POIoverMp);
10121027
}
10131028
if (subGroupStr.Contains("res-flow-dimuon")) {
10141029
int varV2[6] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kR2SP_AB, VarManager::kR2EP_AB};

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 16 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -739,26 +739,32 @@ void VarManager::SetDefaultVarNames()
739739
fgVariableUnits[kS31A] = "";
740740
fgVariableNames[kM11REF] = "M_{11}^{REF} ";
741741
fgVariableUnits[kM11REF] = "";
742-
fgVariableNames[kM01POI] = "M_{01}^{POI} ";
742+
fgVariableNames[kM01POI] = "M^{'}_{01}^{POI} ";
743743
fgVariableUnits[kM01POI] = "";
744744
fgVariableNames[kM1111REF] = "M_{1111}^{REF} ";
745745
fgVariableUnits[kM1111REF] = "";
746-
fgVariableNames[kM0111POI] = "M_{0111}^{POI} ";
746+
fgVariableNames[kM0111POI] = "M^{'}_{0111}^{POI} ";
747747
fgVariableUnits[kM0111POI] = "";
748748
fgVariableNames[kCORR2REF] = "<2> ";
749749
fgVariableUnits[kCORR2REF] = "";
750-
fgVariableNames[kCORR2POI] = " <2'> ";
750+
fgVariableNames[kCORR2POI] = "<2'> ";
751751
fgVariableUnits[kCORR2POI] = "";
752-
fgVariableNames[kCORR4REF] = " <4>";
752+
fgVariableNames[kCORR4REF] = "<4> ";
753753
fgVariableUnits[kCORR4REF] = "";
754754
fgVariableNames[kCORR4POI] = "<4'> ";
755755
fgVariableUnits[kCORR4POI] = "";
756-
fgVariableNames[kC4REF] = "c_{2}(4)";
757-
fgVariableUnits[kC4REF] = "";
758-
fgVariableNames[kC4POI] = "d_{2}(4)";
759-
fgVariableUnits[kC4POI] = "";
760-
fgVariableNames[kV4] = "v_{2}(4)";
761-
fgVariableUnits[kV4] = "";
756+
fgVariableNames[kM11REFoverMp] = "M_{11}^{REF}/M_{p} ";
757+
fgVariableUnits[kM11REFoverMp] = "";
758+
fgVariableNames[kM01POIoverMp] = "M^{'}_{01}^{POI}/M_{p} ";
759+
fgVariableUnits[kM01POIoverMp] = "";
760+
fgVariableNames[kM1111REFoverMp] = "M_{1111}^{REF}/M_{p} ";
761+
fgVariableUnits[kM1111REFoverMp] = "";
762+
fgVariableNames[kM0111POIoverMp] = "M^{'}_{0111}^{POI}/M_{p} ";
763+
fgVariableUnits[kM0111POIoverMp] = "";
764+
fgVariableNames[kCORR2POIMp] = "<2'> M_{p} ";
765+
fgVariableUnits[kCORR2POIMp] = "";
766+
fgVariableNames[kCORR4POIMp] = "<4'> M_{p} ";
767+
fgVariableUnits[kCORR4POIMp] = "";
762768
fgVariableNames[kCos2DeltaPhi] = "cos 2(#varphi-#Psi_{2}^{A}) ";
763769
fgVariableUnits[kCos2DeltaPhi] = "";
764770
fgVariableNames[kCos3DeltaPhi] = "cos 3(#varphi-#Psi_{3}^{A}) ";

‎PWGDQ/Core/VarManager.h‎

Lines changed: 14 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -601,9 +601,12 @@ class VarManager : public TObject
601601
kCORR2POI,
602602
kCORR4REF,
603603
kCORR4POI,
604-
kC4REF,
605-
kC4POI,
606-
kV4,
604+
kM11REFoverMp,
605+
kM01POIoverMp,
606+
kM1111REFoverMp,
607+
kM0111POIoverMp,
608+
kCORR2POIMp,
609+
kCORR4POIMp,
607610
kPsi2A,
608611
kPsi2APOS,
609612
kPsi2ANEG,
@@ -3914,6 +3917,14 @@ void VarManager::FillPairVn(T1 const& t1, T2 const& t2, float* values)
39143917
values[kM0111POI] = values[kMultDimuons] * (values[kS31A] - 3. * values[kS11A] * values[kS12A] + 2. * values[kS13A]);
39153918
values[kCORR2POI] = (P2 * conj(Q21)).real() / values[kM01POI];
39163919
values[kCORR4POI] = (P2 * Q21 * conj(Q21) * conj(Q21) - P2 * Q21 * conj(Q42) - 2. * values[kS12A] * P2 * conj(Q21) + 2. * P2 * conj(Q23)).real() / values[kM0111POI];
3920+
values[kM01POIoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM01POI]) || isinf(values[kM01POI])) ? values[kM01POI] / values[kMultDimuons] : 0;
3921+
values[kM0111POIoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM0111POI]) || isinf(values[kM0111POI])) ? values[kM0111POI] / values[kMultDimuons] : 0;
3922+
values[kM11REFoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM11REF]) || isinf(values[kM11REF])) ? values[kM11REF] / values[kMultDimuons] : 0;
3923+
values[kM1111REFoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM1111REF]) || isinf(values[kM1111REF])) ? values[kM1111REF] / values[kMultDimuons] : 0;
3924+
values[kCORR2POIMp] = isnan(values[kCORR2POI]) || isinf(values[kCORR2POI]) ? 0 : values[kCORR2POI] * values[kMultDimuons];
3925+
values[kCORR4POIMp] = isnan(values[kCORR4POI]) || isinf(values[kCORR4POI]) ? 0 : values[kCORR4POI] * values[kMultDimuons];
3926+
values[kCORR2REF] = isnan(values[kCORR2REF]) || isinf(values[kCORR2REF]) ? 0 : values[kCORR2REF];
3927+
values[kCORR4REF] = isnan(values[kCORR4REF]) || isinf(values[kCORR4REF]) ? 0 : values[kCORR4REF];
39173928
}
39183929
}
39193930

0 commit comments

Comments
 (0)