Skip to content

Commit 9a1b6f0

Browse files
ilikmetailikmeta
committed
[PWGCF] Add eta gap on c3{4}
1 parent 92b0296 commit 9a1b6f0

File tree

1 file changed

+33
-2
lines changed

1 file changed

+33
-2
lines changed

‎PWGCF/Flow/Tasks/flowGfwTask.cxx‎

Lines changed: 33 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -172,6 +172,7 @@ struct FlowGfwTask {
172172
kc32,
173173
kc32etagap,
174174
kc34,
175+
kc34etagap,
175176
kc22Nch,
176177
kc24Nch,
177178
kc26Nch,
@@ -180,6 +181,7 @@ struct FlowGfwTask {
180181
kc32Nch,
181182
kc32Nchetagap,
182183
kc34Nch,
184+
kc34Nchetagap,
183185
kc22Nch05,
184186
kc24Nch05,
185187
kc26Nch05,
@@ -188,9 +190,13 @@ struct FlowGfwTask {
188190
kc32Nch05,
189191
kc32Nch05etagap,
190192
kc34Nch05,
193+
kc34Nch05etagap,
194+
kc22ft0c,
191195
kc22etagapft0c,
196+
kc32ft0c,
192197
kc32etagapft0c,
193198
kc34ft0c,
199+
kc34etagapft0c,
194200

195201
// Count the total number of enum
196202
kCount_ExtraProfile
@@ -387,6 +393,7 @@ struct FlowGfwTask {
387393
registry.add("c32", ";Centrality (%) ; C_{3}{2} ", {HistType::kTProfile, {axisCentrality}});
388394
registry.add("c32etagap", ";Centrality (%) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisCentrality}});
389395
registry.add("c34", ";Centrality (%) ; C_{3}{4} ", {HistType::kTProfile, {axisCentrality}});
396+
registry.add("c34etagap", ";Centrality (%) ; C_{3}{4} (|#eta| < 0.8) ", {HistType::kTProfile, {axisCentrality}});
390397

391398
registry.add("c22Nch", ";N_{ch}(|#eta| < 0.8) ; C_{2}{2} ", {HistType::kTProfile, {axisNch}});
392399
registry.add("c24Nch", ";N_{ch}(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}});
@@ -396,6 +403,7 @@ struct FlowGfwTask {
396403
registry.add("c32Nch", ";N_{ch}(|#eta| < 0.8) ; C_{3}{2} ", {HistType::kTProfile, {axisNch}});
397404
registry.add("c32Nchetagap", ";N_ch(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisNch}});
398405
registry.add("c34Nch", ";N_{ch}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}});
406+
registry.add("c34Nchetagap", ";N_{ch}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}});
399407

400408
registry.add("c22Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{2}{2} ", {HistType::kTProfile, {axisNch}});
401409
registry.add("c24Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}});
@@ -405,9 +413,13 @@ struct FlowGfwTask {
405413
registry.add("c32Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{2} ", {HistType::kTProfile, {axisNch}});
406414
registry.add("c32Nch05etagap", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisNch}});
407415
registry.add("c34Nch05", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}});
416+
registry.add("c34Nch05etagap", ";N_{ch 0-5%}(|#eta| < 0.8) ; C_{3}{4} ", {HistType::kTProfile, {axisNch}});
408417

418+
registry.add("c22ft0c", ";FT0C Amplitude ; C_{2}{2} ", {HistType::kTProfile, {axisFT0CAmp}});
409419
registry.add("c22etagapft0c", ";FT0C Amplitude ; C_{2}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}});
420+
registry.add("c32ft0c", ";FT0C Amplitude ; C_{2}{2} ", {HistType::kTProfile, {axisFT0CAmp}});
410421
registry.add("c32etagapft0c", ";FT0C Amplitude ; C_{3}{2} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}});
422+
registry.add("c34ft0c", ";FT0C Amplitude ; C_{3}{4} ", {HistType::kTProfile, {axisFT0CAmp}});
411423
registry.add("c34etagapft0c", ";FT0C Amplitude ; C_{3}{4} (|#eta| < 0.8) ", {HistType::kTProfile, {axisFT0CAmp}});
412424
} // End doprocessData
413425

@@ -471,6 +483,7 @@ struct FlowGfwTask {
471483
bootstrapArray[i][kc32] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32", i), ";Centrality (%) ; C_{3}{2}", {HistType::kTProfile, {axisCentrality}});
472484
bootstrapArray[i][kc32etagap] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32etagap", i), ";Centrality (%) ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisCentrality}});
473485
bootstrapArray[i][kc34] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34", i), ";Centrality (%) ; C_{3}{4}", {HistType::kTProfile, {axisCentrality}});
486+
bootstrapArray[i][kc34etagap] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34etagap", i), ";Centrality (%) ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisCentrality}});
474487

475488
bootstrapArray[i][kc22Nch] = registry.add<TProfile>(Form("BootstrapContainer_%d/c22Nch", i), ";N_ch(|#eta| < 0.8) ; C_{2}{2}", {HistType::kTProfile, {axisNch}});
476489
bootstrapArray[i][kc24Nch] = registry.add<TProfile>(Form("BootstrapContainer_%d/c24Nch", i), ";N_ch(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}});
@@ -480,6 +493,7 @@ struct FlowGfwTask {
480493
bootstrapArray[i][kc32Nch] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32Nch", i), ";N_ch(|#eta| < 0.8) ; C_{3}{2}", {HistType::kTProfile, {axisNch}});
481494
bootstrapArray[i][kc32Nchetagap] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32Nchetagap", i), ";N_ch(|#eta| < 0.8) ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}});
482495
bootstrapArray[i][kc34Nch] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34Nch", i), ";N_ch(|#eta| < 0.8) ; C_{3}{4}", {HistType::kTProfile, {axisNch}});
496+
bootstrapArray[i][kc34Nchetagap] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34Nchetagap", i), ";N_ch(|#eta| < 0.8) ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}});
483497

484498
bootstrapArray[i][kc22Nch05] = registry.add<TProfile>(Form("BootstrapContainer_%d/c22Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{2}{2}", {HistType::kTProfile, {axisNch}});
485499
bootstrapArray[i][kc24Nch05] = registry.add<TProfile>(Form("BootstrapContainer_%d/c24Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{2}{4}", {HistType::kTProfile, {axisNch}});
@@ -489,10 +503,14 @@ struct FlowGfwTask {
489503
bootstrapArray[i][kc32Nch05] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{2}", {HistType::kTProfile, {axisNch}});
490504
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}});
491505
bootstrapArray[i][kc34Nch05] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34Nch05", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{4}", {HistType::kTProfile, {axisNch}});
506+
bootstrapArray[i][kc34Nch05etagap] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34Nch05etagap", i), ";N_ch05(|#eta| < 0.8) ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisNch}});
492507

508+
bootstrapArray[i][kc22ft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c22ftoc", i), ";FT0C Amplitude ; C_{2}{2}", {HistType::kTProfile, {axisFT0CAmp}});
493509
bootstrapArray[i][kc22etagapft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c22etagapftoc", i), ";FT0C Amplitude ; C_{2}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}});
510+
bootstrapArray[i][kc32ft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32ftoc", i), ";FT0C Amplitude ; C_{3}{2}", {HistType::kTProfile, {axisFT0CAmp}});
494511
bootstrapArray[i][kc32etagapft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c32etagapftoc", i), ";FT0C Amplitude ; C_{3}{2} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}});
495-
bootstrapArray[i][kc34ft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34etagapftoc", i), ";FT0C Amplitude ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}});
512+
bootstrapArray[i][kc34ft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34ftoc", i), ";FT0C Amplitude ; C_{3}{4}", {HistType::kTProfile, {axisFT0CAmp}});
513+
bootstrapArray[i][kc34etagapft0c] = registry.add<TProfile>(Form("BootstrapContainer_%d/c34ftocetagap", i), ";FT0C Amplitude ; C_{3}{4} (|#eta| < 0.8)", {HistType::kTProfile, {axisFT0CAmp}});
496514
}
497515

498516
o2::framework::AxisSpec axis = axisPt;
@@ -524,6 +542,7 @@ struct FlowGfwTask {
524542
corrconfigs.push_back(fGFW->GetCorrelatorConfig("full {3 -3}", "ChFull32", kFALSE));
525543
corrconfigs.push_back(fGFW->GetCorrelatorConfig("refN10 {3} refP10 {-3}", "Ch10Gap32", kFALSE));
526544
corrconfigs.push_back(fGFW->GetCorrelatorConfig("full {3 3 -3 -3}", "ChFull34", kFALSE));
545+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("refN10 {3 3} refP10 {-3 -3}", "Ch10Gap34", kFALSE));
527546
fGFW->CreateRegions(); // finalize the initialization
528547

529548
if (cfgUseAdditionalEventCut) {
@@ -1078,6 +1097,7 @@ struct FlowGfwTask {
10781097
fillProfile(corrconfigs.at(5), HIST("c32"), centrality);
10791098
fillProfile(corrconfigs.at(6), HIST("c32etagap"), centrality);
10801099
fillProfile(corrconfigs.at(7), HIST("c34"), centrality);
1100+
fillProfile(corrconfigs.at(8), HIST("c34etagap"), centrality);
10811101

10821102
fillProfile(corrconfigs.at(0), HIST("c22Nch"), nch);
10831103
fillProfile(corrconfigs.at(1), HIST("c24Nch"), nch);
@@ -1087,6 +1107,7 @@ struct FlowGfwTask {
10871107
fillProfile(corrconfigs.at(5), HIST("c32Nch"), nch);
10881108
fillProfile(corrconfigs.at(6), HIST("c32Nchetagap"), nch);
10891109
fillProfile(corrconfigs.at(7), HIST("c34Nch"), nch);
1110+
fillProfile(corrconfigs.at(8), HIST("c34Nchetagap"), nch);
10901111

10911112
// 0-5% centrality Nch
10921113
float zero = 0, five = 5;
@@ -1099,12 +1120,16 @@ struct FlowGfwTask {
10991120
fillProfile(corrconfigs.at(5), HIST("c32Nch05"), nch);
11001121
fillProfile(corrconfigs.at(6), HIST("c32Nch05etagap"), nch);
11011122
fillProfile(corrconfigs.at(7), HIST("c34Nch05"), nch);
1123+
fillProfile(corrconfigs.at(8), HIST("c34Nch05etagap"), nch);
11021124
}
11031125

11041126
// C22, C32 and C34 vs FT0C amplitude
1127+
fillProfile(corrconfigs.at(0), HIST("c22ft0c"), ft0cAmp);
11051128
fillProfile(corrconfigs.at(4), HIST("c22etagapft0c"), ft0cAmp);
1129+
fillProfile(corrconfigs.at(5), HIST("c32ft0c"), ft0cAmp);
11061130
fillProfile(corrconfigs.at(6), HIST("c32etagapft0c"), ft0cAmp);
1107-
fillProfile(corrconfigs.at(7), HIST("c34etagapft0c"), ft0cAmp);
1131+
fillProfile(corrconfigs.at(7), HIST("c34ft0c"), ft0cAmp);
1132+
fillProfile(corrconfigs.at(8), HIST("c34etagapft0c"), ft0cAmp);
11081133

11091134
// Filling Bootstrap Samples
11101135
int sampleIndex = static_cast<int>(cfgNbootstrap * lRandom);
@@ -1116,6 +1141,7 @@ struct FlowGfwTask {
11161141
fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32], centrality);
11171142
fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32etagap], centrality);
11181143
fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34], centrality);
1144+
fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34etagap], centrality);
11191145

11201146
fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22Nch], nch);
11211147
fillProfile(corrconfigs.at(1), bootstrapArray[sampleIndex][kc24Nch], nch);
@@ -1125,6 +1151,7 @@ struct FlowGfwTask {
11251151
fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32Nch], nch);
11261152
fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32Nchetagap], nch);
11271153
fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34Nch], nch);
1154+
fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34Nchetagap], nch);
11281155

11291156
if (centrality >= zero && centrality <= five) {
11301157
fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22Nch05], nch);
@@ -1135,16 +1162,20 @@ struct FlowGfwTask {
11351162
fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32Nch05], nch);
11361163
fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32Nch05etagap], nch);
11371164
fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34Nch05], nch);
1165+
fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34Nch05etagap], nch);
11381166

11391167
registry.fill(HIST("Nch05"), nch);
11401168
}
11411169

11421170
registry.fill(HIST("Nch"), nch);
11431171

11441172
// Filling Bootstrap Samples for FT0C Amplitudes
1173+
fillProfile(corrconfigs.at(0), bootstrapArray[sampleIndex][kc22ft0c], ft0cAmp);
11451174
fillProfile(corrconfigs.at(4), bootstrapArray[sampleIndex][kc22etagapft0c], ft0cAmp);
1175+
fillProfile(corrconfigs.at(5), bootstrapArray[sampleIndex][kc32ft0c], ft0cAmp);
11461176
fillProfile(corrconfigs.at(6), bootstrapArray[sampleIndex][kc32etagapft0c], ft0cAmp);
11471177
fillProfile(corrconfigs.at(7), bootstrapArray[sampleIndex][kc34ft0c], ft0cAmp);
1178+
fillProfile(corrconfigs.at(8), bootstrapArray[sampleIndex][kc34etagapft0c], ft0cAmp);
11481179

11491180
// Filling Flow Container
11501181
for (uint l_ind = 0; l_ind < corrconfigs.size(); l_ind++) {

0 commit comments

Comments
 (0)