Skip to content

Commit 509b9e4

Browse files
authored
Modify chi2Bins and add chi2 histograms
Updated chi2Bins range and added new histograms for chi2 over NCls for TPC and ITS.
1 parent 9605a27 commit 509b9e4

File tree

1 file changed

+11
-3
lines changed

1 file changed

+11
-3
lines changed

PWGLF/Tasks/QC/systematicsMapping.cxx

Lines changed: 11 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -43,7 +43,7 @@ struct SystematicsMapping {
4343
ConfigurableAxis tpcCrossedRowsBins{"tpcCrossedRowsBins", {5, 70, 100, 120, 135, 150}, "Min TPC clusters for tracks"};
4444
ConfigurableAxis itsClustersBins{"itsClustersBins", {5, 0, 6}, "Min ITS clusters for tracks"};
4545
ConfigurableAxis dcaBins{"dcaBins", {100, 0.f, 5.f}, "Binning for DCA (cm)"};
46-
ConfigurableAxis chi2Bins{"chi2Bins", {100, 0.f, 10.f}, "Binning for chi2"};
46+
ConfigurableAxis chi2Bins{"chi2Bins", {100, 0.f, 100.f}, "Binning for chi2"};
4747
// Selection configurables
4848
Configurable<float> selectionPosZ{"selectionPosZ", 10.f, "Max |z| of the primary vertex"};
4949

@@ -68,7 +68,8 @@ struct SystematicsMapping {
6868
registry.add("K/hITSClusters", "", HistType::kTH1F, {{10, 0, 10}});
6969
registry.add("K/hDCAxy", "", HistType::kTH1F, {dcaBins});
7070
registry.add("K/hDCAz", "", HistType::kTH1F, {dcaBins});
71-
registry.add("K/hChi2", "", HistType::kTH1F, {chi2Bins});
71+
registry.add("K/hChi2OverNCLsTPC", "", HistType::kTH1F, {chi2Bins});
72+
registry.add("K/hChi2OverNCLsITS", "", HistType::kTH1F, {chi2Bins});
7273
registry.addClone("K/", "K0s/");
7374

7475
// Add the signal histograms
@@ -103,7 +104,8 @@ struct SystematicsMapping {
103104
registry.fill(HIST("K/hITSClusters"), track.itsNCls());
104105
registry.fill(HIST("K/hDCAxy"), track.dcaXY());
105106
registry.fill(HIST("K/hDCAz"), track.dcaZ());
106-
registry.fill(HIST("K/hChi2"), track.tpcChi2NCl());
107+
registry.fill(HIST("K/hChi2OverNCLsTPC"), track.tpcChi2NCl());
108+
registry.fill(HIST("K/hChi2OverNCLsITS"), track.itsChi2NCl());
107109
if (track.sign() > 0)
108110
registry.fill(HIST("K/SignalPositive"), track.pt(), track.eta(), track.phi(), track.tpcNSigmaKa(), track.tpcNClsCrossedRows(), track.itsNCls(), track.dcaXY(), track.dcaZ(), track.tpcChi2NCl());
109111
else
@@ -116,6 +118,12 @@ struct SystematicsMapping {
116118
continue;
117119
const auto& posTrack = v0.posTrack_as<TrackType>();
118120
const auto& negTrack = v0.negTrack_as<TrackType>();
121+
registry.fill(HIST("K0s/hTPCCrossedRows"), std::min(posTrack.tpcNClsCrossedRows(), negTrack.tpcNClsCrossedRows()));
122+
registry.fill(HIST("K0s/hITSClusters"), std::min(posTrack.itsNCls(), negTrack.itsNCls()));
123+
registry.fill(HIST("K0s/hDCAxy"), std::min(posTrack.dcaXY(), negTrack.dcaXY()));
124+
registry.fill(HIST("K0s/hDCAz"), std::min(posTrack.dcaZ(), negTrack.dcaZ()));
125+
registry.fill(HIST("K0s/hChi2OverNCLsTPC"), std::min(posTrack.tpcChi2NCl(), negTrack.tpcChi2NCl()));
126+
registry.fill(HIST("K0s/hChi2OverNCLsITS"), std::min(posTrack.itsChi2NCl(), negTrack.itsChi2NCl()));
119127
registry.fill(HIST("K0s/Signal"), v0.pt(), v0.eta(), v0.phi(), v0.mK0Short() - constants::physics::MassK0Short,
120128
std::min(posTrack.tpcNClsCrossedRows(), negTrack.tpcNClsCrossedRows()),
121129
std::min(posTrack.itsNCls(), negTrack.itsNCls()),

0 commit comments

Comments
 (0)