Skip to content

Commit 19f4831

Browse files
RafRaf11alibuild
andauthored
Update; NucleiHistTask (AliceO2Group#2432)
* Changed ptBinning For ability to analyse lower pT ranges * Feature; Changes to the derived table and small bugfixes. * Bugfix; Fixed bug with filling 3D centrality histograms. * Please consider the following formatting changes * UpdateReduced table size; * Update NucleiHistTask.cxx * Update NucleiHistTask.cxx Feature; match variables to offline version. * Please consider the following formatting changes * Feature; Removed Table + added new species * Please consider the following formatting changes * Fix; fixed getters for histogram filling * Bugfix; * Feature; implemented PVcontributor cuts * Please consider the following formatting changes * NucleiHistTask; Update added TRD track selection cuts * Please consider the following formatting changes * Update; Added histogram for centrality distribution, added global isPvcontributor cut for fillCentHistorgrams. --------- Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent ede3d7d commit 19f4831

1 file changed

Lines changed: 9 additions & 1 deletion

File tree

PWGLF/Tasks/NucleiHistTask.cxx

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -64,6 +64,7 @@ struct NucleiHistTask {
6464
AxisSpec ptAxis = {ptBinning, "#it{p}_{T} (GeV/#it{c})"};
6565
AxisSpec centAxis = {centBinning, "V0M (%)"};
6666
AxisSpec centralityAxis = {100, 0.0, 100.0, "VT0C (%)"};
67+
AxisSpec centralityAxis_extended = {105, 0.0, 105.0, "VT0C (%)"};
6768

6869
// QA histograms
6970
spectra.add("histRecVtxZData", "collision z position", HistType::kTH1F, {{200, -20., +20., "z position (cm)"}});
@@ -78,6 +79,7 @@ struct NucleiHistTask {
7879
spectra.add("histNClusterITS", "Number of Clusters in ITS vs Pt", HistType::kTH2F, {ptAxis, {10, 0.0, 10.0, "nCluster"}});
7980
spectra.add("histChi2TPC", "chi^2 TPC vs Pt", HistType::kTH2F, {ptAxis, {100, 0.0, 5.0, "chi^2"}});
8081
spectra.add("histChi2ITS", "chi^2 ITS vs Pt", HistType::kTH2F, {ptAxis, {500, 0.0, 50.0, "chi^2"}});
82+
spectra.add("histCentrality", "Centrality", HistType::kTH1F, {centralityAxis_extended});
8183

8284
// histograms for Proton
8385
proton_erg.add("histKeepEventData", "skimming histogram (p)", HistType::kTH1F, {{2, -0.5, +1.5, "true: keep event, false: reject event"}});
@@ -949,6 +951,8 @@ struct NucleiHistTask {
949951
void fillCentHistorgrams(const CollisionType& event, const TracksType& tracks)
950952
{
951953

954+
spectra.fill(HIST("histCentrality"), event.centFT0C());
955+
952956
for (auto track : tracks) { // start loop over tracks
953957

954958
float TPCnumberClsFound = track.tpcNClsFound();
@@ -958,7 +962,11 @@ struct NucleiHistTask {
958962
float Chi2perClusterITS = track.itsChi2NCl();
959963

960964
// track cuts
961-
if (TPCnumberClsFound < minTPCnClsFound || TPC_nCls_Crossed_Rows < minNCrossedRowsTPC || RatioCrossedRowsOverFindableTPC < minRatioCrossedRowsTPC || RatioCrossedRowsOverFindableTPC > maxRatioCrossedRowsTPC || Chi2perClusterTPC > maxChi2TPC || Chi2perClusterITS > maxChi2ITS || !(track.passedTPCRefit()) || !(track.passedITSRefit()) || (track.itsNCls()) < minReqClusterITS || !(track.isPVContributor())) {
965+
if (enable_PVcontributor_global && !(track.isPVContributor())) {
966+
continue;
967+
}
968+
969+
if (TPCnumberClsFound < minTPCnClsFound || TPC_nCls_Crossed_Rows < minNCrossedRowsTPC || RatioCrossedRowsOverFindableTPC < minRatioCrossedRowsTPC || RatioCrossedRowsOverFindableTPC > maxRatioCrossedRowsTPC || Chi2perClusterTPC > maxChi2TPC || Chi2perClusterITS > maxChi2ITS || !(track.passedTPCRefit()) || !(track.passedITSRefit()) || (track.itsNCls()) < minReqClusterITS) {
962970
continue;
963971
}
964972

0 commit comments

Comments
 (0)