From 9a1b6f07cb4ebde2244422fc2912163c8a82c7b8 Mon Sep 17 00:00:00 2001 From: ilikmeta Date: Wed, 5 Nov 2025 13:32:45 -0600 Subject: [PATCH] [PWGCF] Add eta gap on c3{4} --- PWGCF/Flow/Tasks/flowGfwTask.cxx | 35 ++++++++++++++++++++++++++++++-- 1 file changed, 33 insertions(+), 2 deletions(-) diff --git a/PWGCF/Flow/Tasks/flowGfwTask.cxx b/PWGCF/Flow/Tasks/flowGfwTask.cxx index 868d53197e4..6e5525dc8ea 100644 --- a/PWGCF/Flow/Tasks/flowGfwTask.cxx +++ b/PWGCF/Flow/Tasks/flowGfwTask.cxx @@ -172,6 +172,7 @@ struct FlowGfwTask { kc32, kc32etagap, kc34, + kc34etagap, kc22Nch, kc24Nch, kc26Nch, @@ -180,6 +181,7 @@ struct FlowGfwTask { kc32Nch, kc32Nchetagap, kc34Nch, + kc34Nchetagap, kc22Nch05, kc24Nch05, kc26Nch05, @@ -188,9 +190,13 @@ struct FlowGfwTask { kc32Nch05, kc32Nch05etagap, kc34Nch05, + kc34Nch05etagap, + kc22ft0c, kc22etagapft0c, + kc32ft0c, kc32etagapft0c, kc34ft0c, + kc34etagapft0c, // Count the total number of enum kCount_ExtraProfile @@ -387,6 +393,7 @@ struct FlowGfwTask { registry.add("c32", ";Centrality (%) ; C_{3}{2} ", {HistType::kTProfile, {axisCentrality}}); registry.add("c32etagap", ";Centrality (%) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisCentrality}}); registry.add("c34", ";Centrality (%) ; C_{3}{4} ", {HistType::kTProfile, {axisCentrality}}); + registry.add("c34etagap", ";Centrality (%) ; C_{3}{4} (|#eta| < 0.8) ", {HistType::kTProfile, {axisCentrality}}); registry.add("c22Nch", ";N_{ch}(|#eta| < 0.8) ; C_{2}{2} ", {HistType::kTProfile, {axisNch}}); registry.add("c24Nch", ";N_{ch}(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}}); @@ -396,6 +403,7 @@ struct FlowGfwTask { registry.add("c32Nch", ";N_{ch}(|#eta| < 0.8) ; C_{3}{2} ", {HistType::kTProfile, {axisNch}}); registry.add("c32Nchetagap", ";N_ch(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisNch}}); registry.add("c34Nch", ";N_{ch}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}}); + registry.add("c34Nchetagap", ";N_{ch}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}}); registry.add("c22Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{2}{2} ", {HistType::kTProfile, {axisNch}}); registry.add("c24Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}}); @@ -405,9 +413,13 @@ struct FlowGfwTask { registry.add("c32Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{2} ", {HistType::kTProfile, {axisNch}}); registry.add("c32Nch05etagap", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisNch}}); registry.add("c34Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}}); + registry.add("c34Nch05etagap", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}}); + registry.add("c22ft0c", ";FT0C Amplitude ; C_{2}{2} ", {HistType::kTProfile, {axisFT0CAmp}}); registry.add("c22etagapft0c", ";FT0C Amplitude ; C_{2}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}}); + registry.add("c32ft0c", ";FT0C Amplitude ; C_{2}{2} ", {HistType::kTProfile, {axisFT0CAmp}}); registry.add("c32etagapft0c", ";FT0C Amplitude ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}}); + registry.add("c34ft0c", ";FT0C Amplitude ; C_{3}{4} ", {HistType::kTProfile, {axisFT0CAmp}}); registry.add("c34etagapft0c", ";FT0C Amplitude ; C_{3}{4} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}}); } // End doprocessData @@ -471,6 +483,7 @@ struct FlowGfwTask { bootstrapArray[i][kc32] = registry.add(Form("BootstrapContainer_%d/c32", i), ";Centrality (%) ; C_{3}{2}", {HistType::kTProfile, {axisCentrality}}); bootstrapArray[i][kc32etagap] = registry.add(Form("BootstrapContainer_%d/c32etagap", i), ";Centrality (%) ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisCentrality}}); bootstrapArray[i][kc34] = registry.add(Form("BootstrapContainer_%d/c34", i), ";Centrality (%) ; C_{3}{4}", {HistType::kTProfile, {axisCentrality}}); + bootstrapArray[i][kc34etagap] = registry.add(Form("BootstrapContainer_%d/c34etagap", i), ";Centrality (%) ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisCentrality}}); bootstrapArray[i][kc22Nch] = registry.add(Form("BootstrapContainer_%d/c22Nch", i), ";N_ch(|#eta| < 0.8) ; C_{2}{2}", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc24Nch] = registry.add(Form("BootstrapContainer_%d/c24Nch", i), ";N_ch(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}}); @@ -480,6 +493,7 @@ struct FlowGfwTask { bootstrapArray[i][kc32Nch] = registry.add(Form("BootstrapContainer_%d/c32Nch", i), ";N_ch(|#eta| < 0.8) ; C_{3}{2}", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc32Nchetagap] = registry.add(Form("BootstrapContainer_%d/c32Nchetagap", i), ";N_ch(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc34Nch] = registry.add(Form("BootstrapContainer_%d/c34Nch", i), ";N_ch(|#eta| < 0.8) ; C_{3}{4}", {HistType::kTProfile, {axisNch}}); + bootstrapArray[i][kc34Nchetagap] = registry.add(Form("BootstrapContainer_%d/c34Nchetagap", i), ";N_ch(|#eta| < 0.8) ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc22Nch05] = registry.add(Form("BootstrapContainer_%d/c22Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{2}{2}", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc24Nch05] = registry.add(Form("BootstrapContainer_%d/c24Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}}); @@ -489,10 +503,14 @@ struct FlowGfwTask { bootstrapArray[i][kc32Nch05] = registry.add(Form("BootstrapContainer_%d/c32Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{2}", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc32Nch05etagap] = registry.add(Form("BootstrapContainer_%d/c32Nch05etagap", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}}); bootstrapArray[i][kc34Nch05] = registry.add(Form("BootstrapContainer_%d/c34Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{4}", {HistType::kTProfile, {axisNch}}); + bootstrapArray[i][kc34Nch05etagap] = registry.add(Form("BootstrapContainer_%d/c34Nch05etagap", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}}); + bootstrapArray[i][kc22ft0c] = registry.add(Form("BootstrapContainer_%d/c22ftoc", i), ";FT0C Amplitude ; C_{2}{2}", {HistType::kTProfile, {axisFT0CAmp}}); bootstrapArray[i][kc22etagapft0c] = registry.add(Form("BootstrapContainer_%d/c22etagapftoc", i), ";FT0C Amplitude ; C_{2}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}}); + bootstrapArray[i][kc32ft0c] = registry.add(Form("BootstrapContainer_%d/c32ftoc", i), ";FT0C Amplitude ; C_{3}{2}", {HistType::kTProfile, {axisFT0CAmp}}); bootstrapArray[i][kc32etagapft0c] = registry.add(Form("BootstrapContainer_%d/c32etagapftoc", i), ";FT0C Amplitude ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}}); - bootstrapArray[i][kc34ft0c] = registry.add(Form("BootstrapContainer_%d/c34etagapftoc", i), ";FT0C Amplitude ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}}); + bootstrapArray[i][kc34ft0c] = registry.add(Form("BootstrapContainer_%d/c34ftoc", i), ";FT0C Amplitude ; C_{3}{4}", {HistType::kTProfile, {axisFT0CAmp}}); + bootstrapArray[i][kc34etagapft0c] = registry.add(Form("BootstrapContainer_%d/c34ftocetagap", i), ";FT0C Amplitude ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}}); } o2::framework::AxisSpec axis = axisPt; @@ -524,6 +542,7 @@ struct FlowGfwTask { corrconfigs.push_back(fGFW->GetCorrelatorConfig("full {3 -3}", "ChFull32", kFALSE)); corrconfigs.push_back(fGFW->GetCorrelatorConfig("refN10 {3} refP10 {-3}", "Ch10Gap32", kFALSE)); corrconfigs.push_back(fGFW->GetCorrelatorConfig("full {3 3 -3 -3}", "ChFull34", kFALSE)); + corrconfigs.push_back(fGFW->GetCorrelatorConfig("refN10 {3 3} refP10 {-3 -3}", "Ch10Gap34", kFALSE)); fGFW->CreateRegions(); // finalize the initialization if (cfgUseAdditionalEventCut) { @@ -1078,6 +1097,7 @@ struct FlowGfwTask { fillProfile(corrconfigs.at(5), HIST("c32"), centrality); fillProfile(corrconfigs.at(6), HIST("c32etagap"), centrality); fillProfile(corrconfigs.at(7), HIST("c34"), centrality); + fillProfile(corrconfigs.at(8), HIST("c34etagap"), centrality); fillProfile(corrconfigs.at(0), HIST("c22Nch"), nch); fillProfile(corrconfigs.at(1), HIST("c24Nch"), nch); @@ -1087,6 +1107,7 @@ struct FlowGfwTask { fillProfile(corrconfigs.at(5), HIST("c32Nch"), nch); fillProfile(corrconfigs.at(6), HIST("c32Nchetagap"), nch); fillProfile(corrconfigs.at(7), HIST("c34Nch"), nch); + fillProfile(corrconfigs.at(8), HIST("c34Nchetagap"), nch); // 0-5% centrality Nch float zero = 0, five = 5; @@ -1099,12 +1120,16 @@ struct FlowGfwTask { fillProfile(corrconfigs.at(5), HIST("c32Nch05"), nch); fillProfile(corrconfigs.at(6), HIST("c32Nch05etagap"), nch); fillProfile(corrconfigs.at(7), HIST("c34Nch05"), nch); + fillProfile(corrconfigs.at(8), HIST("c34Nch05etagap"), nch); } // C22, C32 and C34 vs FT0C amplitude + fillProfile(corrconfigs.at(0), HIST("c22ft0c"), ft0cAmp); fillProfile(corrconfigs.at(4), HIST("c22etagapft0c"), ft0cAmp); + fillProfile(corrconfigs.at(5), HIST("c32ft0c"), ft0cAmp); fillProfile(corrconfigs.at(6), HIST("c32etagapft0c"), ft0cAmp); - fillProfile(corrconfigs.at(7), HIST("c34etagapft0c"), ft0cAmp); + fillProfile(corrconfigs.at(7), HIST("c34ft0c"), ft0cAmp); + fillProfile(corrconfigs.at(8), HIST("c34etagapft0c"), ft0cAmp); // Filling Bootstrap Samples int sampleIndex = static_cast(cfgNbootstrap * lRandom); @@ -1116,6 +1141,7 @@ struct FlowGfwTask { fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32], centrality); fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32etagap], centrality); fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34], centrality); + fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34etagap], centrality); fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22Nch], nch); fillProfile(corrconfigs.at(1), bootstrapArray[sampleIndex][kc24Nch], nch); @@ -1125,6 +1151,7 @@ struct FlowGfwTask { fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32Nch], nch); fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32Nchetagap], nch); fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34Nch], nch); + fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34Nchetagap], nch); if (centrality >= zero && centrality <= five) { fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22Nch05], nch); @@ -1135,6 +1162,7 @@ struct FlowGfwTask { fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32Nch05], nch); fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32Nch05etagap], nch); fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34Nch05], nch); + fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34Nch05etagap], nch); registry.fill(HIST("Nch05"), nch); } @@ -1142,9 +1170,12 @@ struct FlowGfwTask { registry.fill(HIST("Nch"), nch); // Filling Bootstrap Samples for FT0C Amplitudes + fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22ft0c], ft0cAmp); fillProfile(corrconfigs.at(4), bootstrapArray[sampleIndex][kc22etagapft0c], ft0cAmp); + fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32ft0c], ft0cAmp); fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32etagapft0c], ft0cAmp); fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34ft0c], ft0cAmp); + fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34etagapft0c], ft0cAmp); // Filling Flow Container for (uint l_ind = 0; l_ind < corrconfigs.size(); l_ind++) {