Skip to content

Commit ba15dd3

Browse files
SwatiSaha-1997SwatiSaha-1997
authored
[PWGCF] added more event selection cuts and choice for centrality estimation (#12276)
1 parent ea47759 commit ba15dd3

File tree

2 files changed

+39
-11
lines changed

2 files changed

+39
-11
lines changed

‎PWGCF/EbyEFluctuations/Tasks/MeanptFluctuations.cxx‎

Lines changed: 23 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -78,6 +78,9 @@ struct MeanptFluctuations_QA_QnTable {
7878
ConfigurableAxis nchAxis{"nchAxis", {5000, 0.5, 5000.5}, ""};
7979
Configurable<bool> cfgEvSelkNoSameBunchPileup{"cfgEvSelkNoSameBunchPileup", true, "Pileup removal"};
8080
Configurable<bool> cfgUseGoodITSLayerAllCut{"cfgUseGoodITSLayerAllCut", true, "Remove time interval with dead ITS zone"};
81+
Configurable<bool> cfgEvSelkNoITSROFrameBorder{"cfgEvSelkNoITSROFrameBorder", true, "ITSROFrame border event selection cut"};
82+
Configurable<bool> cfgEvSelkNoTimeFrameBorder{"cfgEvSelkNoTimeFrameBorder", true, "TimeFrame border event selection cut"};
83+
Configurable<int> cfgCentralityEstimator{"cfgCentralityEstimator", 1, "Centrlaity estimatore choice: 1-->FT0C, 2-->FT0A; 3-->FT0M, 4-->FV0A"};
8184

8285
O2_DEFINE_CONFIGURABLE(cfgUse22sEventCut, bool, true, "Use 22s event cut on mult correlations")
8386

@@ -93,7 +96,7 @@ struct MeanptFluctuations_QA_QnTable {
9396
HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject};
9497

9598
// filtering collisions and tracks***********
96-
using aodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::Mults>>;
99+
using aodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::CentFT0Ms, aod::CentFT0As, aod::CentFV0As, aod::Mults>>;
97100
// using aodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels>>;
98101
using aodTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA>>;
99102

@@ -192,18 +195,35 @@ struct MeanptFluctuations_QA_QnTable {
192195
if (cfgEvSelkNoSameBunchPileup && !(coll.selection_bit(o2::aod::evsel::kNoSameBunchPileup))) {
193196
return;
194197
}
198+
if (cfgEvSelkNoITSROFrameBorder && !(coll.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))) {
199+
return;
200+
}
201+
if (cfgEvSelkNoTimeFrameBorder && !(coll.selection_bit(o2::aod::evsel::kNoTimeFrameBorder))) {
202+
return;
203+
}
195204

196205
const auto CentralityFT0C = coll.centFT0C();
197206
if (cfgUse22sEventCut && !eventSelected(coll, inputTracks.size(), CentralityFT0C))
198207
return;
199208

200209
histos.fill(HIST("hZvtx_after_sel"), coll.posZ());
201-
histos.fill(HIST("hCentrality"), coll.centFT0C());
210+
211+
double cent = 0.0;
212+
if (cfgCentralityEstimator == 1)
213+
cent = coll.centFT0C();
214+
else if (cfgCentralityEstimator == 2)
215+
cent = coll.centFT0A();
216+
else if (cfgCentralityEstimator == 3)
217+
cent = coll.centFT0M();
218+
else if (cfgCentralityEstimator == 4)
219+
cent = coll.centFV0A();
220+
221+
histos.fill(HIST("hCentrality"), cent);
222+
202223
histos.fill(HIST("Hist2D_globalTracks_PVTracks"), coll.multNTracksPV(), inputTracks.size());
203224
histos.fill(HIST("Hist2D_cent_nch"), inputTracks.size(), CentralityFT0C);
204225

205226
// variables
206-
double cent = coll.centFT0C();
207227
double pT_sum = 0.0;
208228
double N = 0.0;
209229

‎PWGCF/EbyEFluctuations/Tasks/v0ptHadPiKaProt.cxx‎

Lines changed: 16 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -87,9 +87,11 @@ struct V0ptHadPiKaProt {
8787
Configurable<float> cfgCutEtaLeft{"cfgCutEtaLeft", 0.8f, "Left end of eta gap"};
8888
Configurable<float> cfgCutEtaRight{"cfgCutEtaRight", 0.8f, "Right end of eta gap"};
8989
Configurable<int> cfgNSubsample{"cfgNSubsample", 10, "Number of subsamples"};
90-
Configurable<int> cfgCentralityChoice{"cfgCentralityChoice", 1, "Which centrality estimator? 0-->FT0M, 1-->FT0C"};
90+
Configurable<int> cfgCentralityChoice{"cfgCentralityChoice", 1, "Which centrality estimator? 1-->FT0C, 2-->FT0A, 3-->FT0M, 4-->FV0A"};
9191
Configurable<bool> cfgEvSelkNoSameBunchPileup{"cfgEvSelkNoSameBunchPileup", true, "Pileup removal"};
9292
Configurable<bool> cfgUseGoodITSLayerAllCut{"cfgUseGoodITSLayerAllCut", true, "Remove time interval with dead ITS zone"};
93+
Configurable<bool> cfgEvSelkNoITSROFrameBorder{"cfgEvSelkNoITSROFrameBorder", true, "ITSROFrame border event selection cut"};
94+
Configurable<bool> cfgEvSelkNoTimeFrameBorder{"cfgEvSelkNoTimeFrameBorder", true, "TimeFrame border event selection cut"};
9395

9496
// Connect to ccdb
9597
Service<ccdb::BasicCCDBManager> ccdb;
@@ -342,13 +344,23 @@ struct V0ptHadPiKaProt {
342344
if (cfgEvSelkNoSameBunchPileup && !(coll.selection_bit(o2::aod::evsel::kNoSameBunchPileup))) {
343345
return;
344346
}
347+
if (cfgEvSelkNoITSROFrameBorder && !(coll.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))) {
348+
return;
349+
}
350+
if (cfgEvSelkNoTimeFrameBorder && !(coll.selection_bit(o2::aod::evsel::kNoTimeFrameBorder))) {
351+
return;
352+
}
345353

346354
// Centrality
347355
double cent = 0.0;
348-
if (cfgCentralityChoice == 0)
349-
cent = coll.centFT0M();
350-
else if (cfgCentralityChoice == 1)
356+
if (cfgCentralityChoice == 1)
357+
cent = coll.centFT0C();
358+
else if (cfgCentralityChoice == 2)
359+
cent = coll.centFT0A();
360+
else if (cfgCentralityChoice == 3)
351361
cent = coll.centFT0M();
362+
else if (cfgCentralityChoice == 4)
363+
cent = coll.centFV0A();
352364

353365
histos.fill(HIST("hZvtx_after_sel"), coll.posZ());
354366
histos.fill(HIST("hCentrality"), cent);
@@ -483,7 +495,6 @@ struct V0ptHadPiKaProt {
483495
subSample[sampleIndex][0]->Fill(cent, fPtProfileHad->GetBinCenter(i + 1), (fPtProfileHad->GetBinContent(i + 1) / nSumEtaLeftHad));
484496
subSample[sampleIndex][1]->Fill(cent, fPtProfileHad->GetBinCenter(i + 1), ((fPtProfileHad->GetBinContent(i + 1) / nSumEtaLeftHad) * (pTsumEtaRightHad / nSumEtaRightHad)));
485497
subSample[sampleIndex][2]->Fill(cent, 0.5, ((pTsumEtaLeftHad / nSumEtaLeftHad) * (pTsumEtaRightHad / nSumEtaRightHad)));
486-
;
487498
subSample[sampleIndex][3]->Fill(cent, 0.5, (pTsumEtaLeftHad / nSumEtaLeftHad));
488499
subSample[sampleIndex][4]->Fill(cent, 0.5, (pTsumEtaRightHad / nSumEtaRightHad));
489500
}
@@ -500,7 +511,6 @@ struct V0ptHadPiKaProt {
500511
subSample[sampleIndex][5]->Fill(cent, fPtProfilePi->GetBinCenter(i + 1), (fPtProfilePi->GetBinContent(i + 1) / nSumEtaLeftPi));
501512
subSample[sampleIndex][6]->Fill(cent, fPtProfilePi->GetBinCenter(i + 1), ((fPtProfilePi->GetBinContent(i + 1) / nSumEtaLeftPi) * (pTsumEtaRightHad / nSumEtaRightHad)));
502513
subSample[sampleIndex][7]->Fill(cent, 0.5, ((pTsumEtaLeftHad / nSumEtaLeftHad) * (pTsumEtaRightHad / nSumEtaRightHad)));
503-
;
504514
subSample[sampleIndex][8]->Fill(cent, 0.5, (pTsumEtaLeftHad / nSumEtaLeftHad));
505515
subSample[sampleIndex][9]->Fill(cent, 0.5, (pTsumEtaRightHad / nSumEtaRightHad));
506516
}
@@ -517,7 +527,6 @@ struct V0ptHadPiKaProt {
517527
subSample[sampleIndex][10]->Fill(cent, fPtProfileKa->GetBinCenter(i + 1), (fPtProfileKa->GetBinContent(i + 1) / nSumEtaLeftKa));
518528
subSample[sampleIndex][11]->Fill(cent, fPtProfileKa->GetBinCenter(i + 1), ((fPtProfileKa->GetBinContent(i + 1) / nSumEtaLeftKa) * (pTsumEtaRightHad / nSumEtaRightHad)));
519529
subSample[sampleIndex][12]->Fill(cent, 0.5, ((pTsumEtaLeftHad / nSumEtaLeftHad) * (pTsumEtaRightHad / nSumEtaRightHad)));
520-
;
521530
subSample[sampleIndex][13]->Fill(cent, 0.5, (pTsumEtaLeftHad / nSumEtaLeftHad));
522531
subSample[sampleIndex][14]->Fill(cent, 0.5, (pTsumEtaRightHad / nSumEtaRightHad));
523532
}
@@ -533,7 +542,6 @@ struct V0ptHadPiKaProt {
533542
subSample[sampleIndex][15]->Fill(cent, fPtProfileProt->GetBinCenter(i + 1), (fPtProfileProt->GetBinContent(i + 1) / nSumEtaLeftProt));
534543
subSample[sampleIndex][16]->Fill(cent, fPtProfileProt->GetBinCenter(i + 1), ((fPtProfileProt->GetBinContent(i + 1) / nSumEtaLeftProt) * (pTsumEtaRightHad / nSumEtaRightHad)));
535544
subSample[sampleIndex][17]->Fill(cent, 0.5, ((pTsumEtaLeftHad / nSumEtaLeftHad) * (pTsumEtaRightHad / nSumEtaRightHad)));
536-
;
537545
subSample[sampleIndex][18]->Fill(cent, 0.5, (pTsumEtaLeftHad / nSumEtaLeftHad));
538546
subSample[sampleIndex][19]->Fill(cent, 0.5, (pTsumEtaRightHad / nSumEtaRightHad));
539547
}

0 commit comments

Comments
 (0)