Skip to content

Commit 1e874c7

Browse files
authored
Update histogram configurations for TPC and ITS
1 parent 088b554 commit 1e874c7

File tree

1 file changed

+11
-11
lines changed

1 file changed

+11
-11
lines changed

PWGLF/Tasks/QC/systematicsMapping.cxx

Lines changed: 11 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -40,8 +40,8 @@ struct SystematicsMapping {
4040
ConfigurableAxis invariantMassBins{"invariantMassBins", {100, -0.1f, 0.1f}, "Binning for the invariant mass (GeV/c^2)"};
4141
ConfigurableAxis nsigmaBins{"nsigmaBins", {100, -10.f, 10.f}, "Binning for nSigma"};
4242
// Selection bins
43-
ConfigurableAxis tpcClusterBins{"tpcClusterBins", {5, 70, 100, 120, 135, 150}, "Min TPC clusters for tracks"};
44-
ConfigurableAxis itsClustersBins{"itsClustersBins", {5, 0, 6}, "Min ITS clusters for tracks"};
43+
ConfigurableAxis tpcCrossedRowsBins{"tpcCrossedRowsBins", {5, 70, 100, 120, 135, 150}, "TPC Crossed rows bins for tracks"};
44+
ConfigurableAxis itsClustersBins{"itsClustersBins", {5, 0, 6}, "ITS Crossed rows bins for tracks"};
4545
// Selection configurables
4646
Configurable<float> selectionPosZ{"selectionPosZ", 10.f, "Max |z| of the primary vertex"};
4747

@@ -61,15 +61,15 @@ struct SystematicsMapping {
6161

6262
if (doprocessData) {
6363

64-
// First we define the histograms on which we are cutting (tpc clusters, its clusters, ..)
65-
registry.add("K/hTPCClusters", "", HistType::kTH1F, {{100, 0, 200}});
64+
// First we define the histograms on which we are cutting (tpc crossed rows, its clusters, ..)
65+
registry.add("K/hTPCCrossedRows", "", HistType::kTH1F, {{100, 0, 200}});
6666
registry.add("K/hITSClusters", "", HistType::kTH1F, {{10, 0, 10}});
6767
registry.addClone("K/", "K0s/");
6868

6969
// Add the signal histograms
70-
registry.add("K/SignalPositive", "", HistType::kTHnSparseF, {ptBins, etaBins, phiBins, nsigmaBins, tpcClusterBins, itsClustersBins});
71-
registry.add("K/SignalNegative", "", HistType::kTHnSparseF, {ptBins, etaBins, phiBins, nsigmaBins, tpcClusterBins, itsClustersBins});
72-
registry.add("K0s/Signal", "", HistType::kTHnSparseF, {ptBins, etaBins, phiBins, invariantMassBins, tpcClusterBins, itsClustersBins});
70+
registry.add("K/SignalPositive", "", HistType::kTHnSparseF, {ptBins, etaBins, phiBins, nsigmaBins, tpcCrossedRowsBins, itsClustersBins});
71+
registry.add("K/SignalNegative", "", HistType::kTHnSparseF, {ptBins, etaBins, phiBins, nsigmaBins, tpcCrossedRowsBins, itsClustersBins});
72+
registry.add("K0s/Signal", "", HistType::kTHnSparseF, {ptBins, etaBins, phiBins, invariantMassBins, tpcCrossedRowsBins, itsClustersBins});
7373
}
7474

7575
if (doprocessMc) {
@@ -94,12 +94,12 @@ struct SystematicsMapping {
9494
for (const auto& track : tracks) {
9595
if (track.collisionId() != collision.globalIndex())
9696
continue;
97-
registry.fill(HIST("hTPCClusters"), track.tpcNClsFound());
97+
registry.fill(HIST("hTPCCrossedRows"), track.tpcCrossedRows());
9898
registry.fill(HIST("hITSClusters"), track.itsNCls());
9999
if (track.sign() > 0)
100-
registry.fill(HIST("K/SignalPositive"), track.pt(), track.eta(), track.phi(), track.tpcNSigmaKa(), track.tpcNClsFound(), track.itsNCls());
100+
registry.fill(HIST("K/SignalPositive"), track.pt(), track.eta(), track.phi(), track.tpcNSigmaKa(), track.tpcCrossedRows(), track.itsNCls());
101101
else
102-
registry.fill(HIST("K/SignalNegative"), track.pt(), track.eta(), track.phi(), track.tpcNSigmaKa(), track.tpcNClsFound(), track.itsNCls());
102+
registry.fill(HIST("K/SignalNegative"), track.pt(), track.eta(), track.phi(), track.tpcNSigmaKa(), track.tpcCrossedRows(), track.itsNCls());
103103
}
104104

105105
// K0s loop
@@ -108,7 +108,7 @@ struct SystematicsMapping {
108108
continue;
109109
const auto& posTrack = v0.posTrack_as<TrackType>();
110110
const auto& negTrack = v0.negTrack_as<TrackType>();
111-
registry.fill(HIST("K0s/Signal"), v0.pt(), v0.eta(), v0.phi(), v0.mK0Short() - constants::physics::MassK0Short, std::min(posTrack.tpcNClsFound(), negTrack.tpcNClsFound()), std::min(posTrack.itsNCls(), negTrack.itsNCls()));
111+
registry.fill(HIST("K0s/Signal"), v0.pt(), v0.eta(), v0.phi(), v0.mK0Short() - constants::physics::MassK0Short, std::min(posTrack.tpcCrossedRows(), negTrack.tpcCrossedRows()), std::min(posTrack.itsNCls(), negTrack.itsNCls()));
112112
}
113113
}
114114
}

0 commit comments

Comments
 (0)