Skip to content

Commit 91793e0

Browse files
ilikmetailikmeta
committed
Add plots with FT0C aplitude for UCC
1 parent 686c3db commit 91793e0

File tree

1 file changed

+36
-12
lines changed

1 file changed

+36
-12
lines changed

‎PWGCF/Flow/Tasks/flowGfwTask.cxx‎

Lines changed: 36 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -75,11 +75,11 @@ struct FlowGfwTask {
7575
O2_DEFINE_CONFIGURABLE(cfgTrackSel, bool, false, "ITS and TPC cluster selection")
7676
O2_DEFINE_CONFIGURABLE(cfgMinCentFT0C, float, 0.0f, "Minimum FT0C Centrality")
7777
O2_DEFINE_CONFIGURABLE(cfgMaxCentFT0C, float, 100.0f, "Maximum FT0C Centrality")
78-
O2_DEFINE_CONFIGURABLE(cfgcentEstFt0c, bool, false, "Centrality estimator based on FT0C signal")
79-
O2_DEFINE_CONFIGURABLE(cfgcentEstFt0a, bool, false, "Centrality estimator based on FT0A signal")
80-
O2_DEFINE_CONFIGURABLE(cfgcentEstFt0m, bool, false, " A centrality estimator based on FT0A+FT0C signals.")
81-
O2_DEFINE_CONFIGURABLE(cfgcentEstFv0a, bool, false, "Centrality estimator based on FV0A signal")
82-
O2_DEFINE_CONFIGURABLE(cfgcentEstFt0cVariant1, bool, false, "A variant of FT0C")
78+
O2_DEFINE_CONFIGURABLE(cfgCentEstFt0c, bool, false, "Centrality estimator based on FT0C signal")
79+
O2_DEFINE_CONFIGURABLE(cfgCentEstFt0a, bool, false, "Centrality estimator based on FT0A signal")
80+
O2_DEFINE_CONFIGURABLE(cfgCentEstFt0m, bool, false, " A centrality estimator based on FT0A+FT0C signals.")
81+
O2_DEFINE_CONFIGURABLE(cfgCentEstFv0a, bool, false, "Centrality estimator based on FV0A signal")
82+
O2_DEFINE_CONFIGURABLE(cfgCentEstFt0cVariant1, bool, false, "A variant of FT0C")
8383
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalEventCut, bool, false, "Use additional event cut on mult correlations")
8484
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalTrackCut, bool, false, "Use additional track cut on phi")
8585
O2_DEFINE_CONFIGURABLE(cfgTrackSelRun3ITSMatch, bool, false, "Track selection for ITS matches")
@@ -94,7 +94,7 @@ struct FlowGfwTask {
9494
O2_DEFINE_CONFIGURABLE(cfgCutOccupancyLow, int, 0, "Low cut on TPC occupancy")
9595
O2_DEFINE_CONFIGURABLE(cfgCutDCAz, float, 2, "Custom DCA Z cut")
9696
O2_DEFINE_CONFIGURABLE(cfgCutDCAxy, float, 0.2f, "Custom DCA XY cut")
97-
O2_DEFINE_CONFIGURABLE(cfgDCAzPt, bool, false, "switch for DCAz pt dependent")
97+
O2_DEFINE_CONFIGURABLE(cfgDCAzPt, bool, true, "switch for DCAz pt dependent")
9898
O2_DEFINE_CONFIGURABLE(cfgNoTimeFrameBorder, bool, false, "kNoTimeFrameBorder");
9999
O2_DEFINE_CONFIGURABLE(cfgNoITSROFrameBorder, bool, false, "kNoITSROFrameBorder");
100100
O2_DEFINE_CONFIGURABLE(cfgNoSameBunchPileup, bool, false, "kNoSameBunchPileup");
@@ -187,6 +187,8 @@ struct FlowGfwTask {
187187
kc32Nch05,
188188
kc32Nch05etagap,
189189
kc34Nch05,
190+
kc22etagapft0c,
191+
kc32etagapft0c,
190192

191193
// Count the total number of enum
192194
kCount_ExtraProfile
@@ -400,6 +402,9 @@ struct FlowGfwTask {
400402
registry.add("c32Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{2} ", {HistType::kTProfile, {axisNch}});
401403
registry.add("c32Nch05etagap", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisNch}});
402404
registry.add("c34Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}});
405+
406+
registry.add("c22etagapft0c", ";FT0C Amplitude ; C_{2}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}});
407+
registry.add("c32etagapft0c", ";FT0C Amplitude ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}});
403408
} // End doprocessData
404409

405410
const AxisSpec axisZpos{48, -12., 12., "Vtx_{z} (cm)"};
@@ -480,6 +485,9 @@ struct FlowGfwTask {
480485
bootstrapArray[i][kc32Nch05] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{2}", {HistType::kTProfile, {axisNch}});
481486
bootstrapArray[i][kc32Nch05etagap] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32Nch05etagap", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}});
482487
bootstrapArray[i][kc34Nch05] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{4}", {HistType::kTProfile, {axisNch}});
488+
489+
bootstrapArray[i][kc22etagapft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c22etagapftoc", i), ";FT0C Amplitude ; C_{2}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}});
490+
bootstrapArray[i][kc32etagapft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32etagapftoc", i), ";FT0C Amplitude ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}});
483491
}
484492

485493
o2::framework::AxisSpec axis = axisPt;
@@ -828,27 +836,27 @@ struct FlowGfwTask {
828836

829837
// Choose centrality estimator -- Only one can be true
830838
auto centrality = -1;
831-
if (cfgcentEstFt0c) {
839+
if (cfgCentEstFt0c) {
832840
centrality = collision.centFT0C();
833841
registry.fill(HIST("hCentEstimators"), kCentFT0C);
834842
registry.fill(HIST("hCentFT0C"), centrality);
835843
}
836-
if (cfgcentEstFt0a) {
844+
if (cfgCentEstFt0a) {
837845
centrality = collision.centFT0A();
838846
registry.fill(HIST("hCentEstimators"), kCentFT0A);
839847
registry.fill(HIST("hCentFT0A"), centrality);
840848
}
841-
if (cfgcentEstFt0m) {
849+
if (cfgCentEstFt0m) {
842850
centrality = collision.centFT0M();
843851
registry.fill(HIST("hCentEstimators"), kCentFT0M);
844852
registry.fill(HIST("hCentFT0M"), centrality);
845853
}
846-
if (cfgcentEstFv0a) {
854+
if (cfgCentEstFv0a) {
847855
centrality = collision.centFV0A();
848856
registry.fill(HIST("hCentEstimators"), kCentFV0A);
849857
registry.fill(HIST("hCentFV0A"), centrality);
850858
}
851-
if (cfgcentEstFt0cVariant1) {
859+
if (cfgCentEstFt0cVariant1) {
852860
centrality = collision.centFT0CVariant1();
853861
registry.fill(HIST("hCentEstimators"), kCentFT0CVariant1);
854862
registry.fill(HIST("hCentFT0CVariant1"), centrality);
@@ -938,6 +946,14 @@ struct FlowGfwTask {
938946
}
939947
} // End of ZDC
940948

949+
// Use for c22 vs ft0 amplitude
950+
double ft0cAmp = 0;
951+
if (foundBC.has_ft0()) {
952+
for (const auto& amplitude : foundBC.ft0().amplitudeC()) {
953+
ft0cAmp += amplitude;
954+
}
955+
}
956+
941957
float vtxz = collision.posZ();
942958
float lRandom = fRndm->Rndm();
943959
registry.fill(HIST("hVtxZ"), vtxz);
@@ -1046,7 +1062,7 @@ struct FlowGfwTask {
10461062
registry.fill(HIST("Events_per_Centrality_Bin"), centrality);
10471063
registry.fill(HIST("Tracks_per_Centrality_Bin"), centrality, nch);
10481064

1049-
// Filling c22 with ROOT TProfile
1065+
// Filling Cumulants with ROOT TProfile
10501066
fillProfile(corrconfigs.at(0), HIST("c22"), centrality);
10511067
fillProfile(corrconfigs.at(1), HIST("c24"), centrality);
10521068
fillProfile(corrconfigs.at(2), HIST("c26"), centrality);
@@ -1077,6 +1093,10 @@ struct FlowGfwTask {
10771093
fillProfile(corrconfigs.at(7), HIST("c34Nch05"), nch);
10781094
}
10791095

1096+
// C22 and C32 vs FT0C amplitude
1097+
fillProfile(corrconfigs.at(4), HIST("c22etagapft0c"), ft0cAmp);
1098+
fillProfile(corrconfigs.at(6), HIST("c32etagapft0c"), ft0cAmp);
1099+
10801100
// Filling Bootstrap Samples
10811101
int sampleIndex = static_cast<int>(cfgNbootstrap * lRandom);
10821102
fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22], centrality);
@@ -1112,6 +1132,10 @@ struct FlowGfwTask {
11121132

11131133
registry.fill(HIST("Nch"), nch);
11141134

1135+
// Filling Bootstrap Samples for FT0C Amplitudes
1136+
fillProfile(corrconfigs.at(4), bootstrapArray[sampleIndex][kc22etagapft0c], ft0cAmp);
1137+
fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32etagapft0c], ft0cAmp);
1138+
11151139
// Filling Flow Container
11161140
for (uint l_ind = 0; l_ind < corrconfigs.size(); l_ind++) {
11171141
fillFC(corrconfigs.at(l_ind), centrality, lRandom);

0 commit comments

Comments
 (0)