Skip to content

Commit d2d6efc

Browse files
Shunsuke-KuritaShunsuke-Kurita
authored
Merge branch 'AliceO2Group:master' into master
2 parents 144bb96 + 760886f commit d2d6efc

17 files changed

Lines changed: 655 additions & 519 deletions

File tree

‎Common/TableProducer/qVectorsTable.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -522,14 +522,14 @@ struct qVectorsTable {
522522
coll.centFV0A()};
523523
cent = centAllEstim[cfgCentEsti];
524524
bool IsCalibrated = true;
525-
if (cent < 0. || cent > 80.) {
525+
if (cent < 0. || cent > 90.) {
526526
cent = 110.;
527527
IsCalibrated = false;
528528
}
529529
for (std::size_t id = 0; id < cfgnMods->size(); id++) {
530530
int ind = cfgnMods->at(id);
531531
CalQvec(ind, coll, tracks, qvecRe, qvecIm, qvecAmp, TrkBPosLabel, TrkBNegLabel, TrkBTotLabel);
532-
if (cent < 80) {
532+
if (cent < 90.) {
533533
for (auto i{0u}; i < kBTot + 1; i++) {
534534
helperEP.DoRecenter(qvecRe[(kBTot + 1) * 4 * id + i * 4 + 1], qvecIm[(kBTot + 1) * 4 * id + i * 4 + 1],
535535
objQvec.at(id)->GetBinContent(static_cast<int>(cent) + 1, 1, i + 1), objQvec.at(id)->GetBinContent(static_cast<int>(cent) + 1, 2, i + 1));

‎DPG/Tasks/AOTTrack/qaMatchEff.cxx‎

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -87,14 +87,15 @@ struct qaMatchEff {
8787
Configurable<bool> isPbPb{"isPbPb", false, "Boolean to tag if the data is PbPb collisions. If false, it is pp"};
8888
Configurable<bool> isEnableEventSelection{"isEnableEventSelection", true, "Boolean to switch the event selection on/off."};
8989
Configurable<bool> isCentralityRequired{"isCentralityRequired", false, "Boolean to switch the centrality selection on/off."};
90+
Configurable<bool> isRejectNearByEvent{"isRejectNearByEvent", false, "Boolean to switch the rejection of near by events on/off."};
9091
Configurable<bool> isEnableOccupancyCut{"isEnableOccupancyCut", false, "Boolean to switch the occupancy cut on/off."};
9192
struct : ConfigurableGroup {
9293
Configurable<float> centralityMinCut{"centralityMinCut", 0.0f, "Minimum centrality"};
9394
Configurable<float> centralityMaxCut{"centralityMaxCut", 100.0f, "Maximum centrality"};
9495
} centralityCuts;
9596
struct : ConfigurableGroup {
9697
Configurable<int> minTracksInTimeRange{"minTracksInTimeRange", 0, "Minimum number of tracks in the time range"};
97-
Configurable<int> maxTracksInTimeRange{"maxTracksInTimeRange", 1000, "Maximum number of tracks in the time range"};
98+
Configurable<int> maxTracksInTimeRange{"maxTracksInTimeRange", 999999, "Maximum number of tracks in the time range"};
9899
} occupancyCuts;
99100
//
100101
// Track selections
@@ -3363,6 +3364,13 @@ struct qaMatchEff {
33633364
return;
33643365
}
33653366
}
3367+
if (isRejectNearByEvent) {
3368+
if (!collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
3369+
if (doDebug)
3370+
LOGF(info, "Nearby event found, skipping...");
3371+
return;
3372+
}
3373+
}
33663374
if (isEnableOccupancyCut) {
33673375
if (occupancy < occupancyCuts.minTracksInTimeRange || occupancy > occupancyCuts.maxTracksInTimeRange) {
33683376
if (doDebug)

‎PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx‎

Lines changed: 39 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -312,10 +312,10 @@ struct femtoUniverseProducerTask {
312312

313313
void init(InitContext&)
314314
{
315-
if ((doprocessFullData || doprocessTrackPhiData || doprocessTrackData || doprocessTrackV0 || doprocessTrackD0mesonData || doprocessTrackCentRun2Data || doprocessTrackCentRun3Data || doprocessTrackV0CentRun3) == false && (doprocessFullMC || doprocessTrackMC || doprocessTrackMCTruth) == false) {
315+
if ((doprocessFullData || doprocessTrackPhiData || doprocessTrackData || doprocessTrackV0 || doprocessTrackD0mesonData || doprocessTrackCentRun2Data || doprocessTrackCentRun3Data || doprocessTrackV0CentRun3) == false && (doprocessFullMC || doprocessTrackMC || doprocessTrackMCTruth || doprocessTruthAndFullMC) == false) {
316316
LOGF(fatal, "Neither processFullData nor processFullMC enabled. Please choose one.");
317317
}
318-
if ((doprocessFullData || doprocessTrackPhiData || doprocessTrackData || doprocessTrackV0 || doprocessTrackD0mesonData || doprocessTrackCentRun2Data || doprocessTrackCentRun3Data || doprocessTrackV0CentRun3) == true && (doprocessFullMC || doprocessTrackMC || doprocessTrackMCTruth) == true) {
318+
if ((doprocessFullData || doprocessTrackPhiData || doprocessTrackData || doprocessTrackV0 || doprocessTrackD0mesonData || doprocessTrackCentRun2Data || doprocessTrackCentRun3Data || doprocessTrackV0CentRun3) == true && (doprocessFullMC || doprocessTrackMC || doprocessTrackMCTruth || doprocessTruthAndFullMC) == true) {
319319
LOGF(fatal,
320320
"Cannot enable process Data and process MC at the same time. "
321321
"Please choose one.");
@@ -999,7 +999,7 @@ struct femtoUniverseProducerTask {
999999
}
10001000
}
10011001

1002-
template <typename TrackType>
1002+
template <typename TrackType, bool transientLabels = false>
10031003
void fillParticles(TrackType const& tracks)
10041004
{
10051005
std::vector<int> childIDs = {0, 0}; // these IDs are necessary to keep track of the children
@@ -1058,6 +1058,12 @@ struct femtoUniverseProducerTask {
10581058
if (ConfIsDebug) {
10591059
fillDebugParticle<false, true>(particle);
10601060
}
1061+
1062+
// Workaround to keep the FDParticles and MC label tables
1063+
// aligned, so that they can be joined in the task.
1064+
if constexpr (transientLabels) {
1065+
outputPartsMCLabels(-1);
1066+
}
10611067
}
10621068
}
10631069

@@ -1211,6 +1217,36 @@ struct femtoUniverseProducerTask {
12111217
}
12121218
PROCESS_SWITCH(femtoUniverseProducerTask, processTrackMCTruth, "Provide MC data for MC truth track analysis", false);
12131219

1220+
Preslice<aod::McParticles> perMCCollision = aod::mcparticle::mcCollisionId;
1221+
PresliceUnsorted<soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::McCollisionLabels>> recoCollsPerMCColl = aod::mcparticle::mcCollisionId;
1222+
Preslice<soa::Join<aod::FemtoFullTracks, aod::McTrackLabels>> perCollisionTracks = aod::track::collisionId;
1223+
Preslice<soa::Join<o2::aod::V0Datas, aod::McV0Labels>> perCollisionV0s = aod::track::collisionId;
1224+
void processTruthAndFullMC(
1225+
aod::McCollisions const& mccols,
1226+
aod::McParticles const& mcParticles,
1227+
soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::McCollisionLabels> const& collisions,
1228+
soa::Join<aod::FemtoFullTracks, aod::McTrackLabels> const& tracks,
1229+
soa::Join<o2::aod::V0Datas, aod::McV0Labels> const& fullV0s,
1230+
aod::BCsWithTimestamps const&)
1231+
{
1232+
// truth
1233+
for (auto& mccol : mccols) {
1234+
auto groupedMCParticles = mcParticles.sliceBy(perMCCollision, mccol.globalIndex());
1235+
auto groupedCollisions = collisions.sliceBy(recoCollsPerMCColl, mccol.globalIndex());
1236+
fillMCTruthCollisions(groupedCollisions, groupedMCParticles); // fills the reco collisions for mc collision
1237+
fillParticles<decltype(groupedMCParticles), true>(groupedMCParticles); // fills mc particles
1238+
}
1239+
1240+
// recos
1241+
for (auto& col : collisions) {
1242+
auto groupedTracks = tracks.sliceBy(perCollisionTracks, col.globalIndex());
1243+
auto groupedV0s = fullV0s.sliceBy(perCollisionV0s, col.globalIndex());
1244+
getMagneticFieldTesla(col.bc_as<aod::BCsWithTimestamps>());
1245+
fillCollisionsAndTracksAndV0AndPhi<true>(col, groupedTracks, groupedV0s);
1246+
}
1247+
}
1248+
PROCESS_SWITCH(femtoUniverseProducerTask, processTruthAndFullMC, "Provide both MC truth and reco for tracks and V0s", false);
1249+
12141250
void processTrackCentRun2Data(aod::FemtoFullCollisionCentRun2 const& col,
12151251
aod::BCsWithTimestamps const&,
12161252
aod::FemtoFullTracks const& tracks)

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 14 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -392,6 +392,13 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
392392
return cut;
393393
}
394394

395+
if (!nameStr.compare("emu_electronCuts")) {
396+
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
397+
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug4"));
398+
cut->AddCut(GetAnalysisCut("electronPIDnsigmaSkewed"));
399+
return cut;
400+
}
401+
395402
if (!nameStr.compare("jpsiKineAndQuality")) {
396403
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
397404
cut->AddCut(GetAnalysisCut("electronStandardQuality"));
@@ -3763,6 +3770,13 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
37633770
return cut;
37643771
}
37653772

3773+
if (!nameStr.compare("electronStandardQualityForO2MCdebug4")) {
3774+
cut->AddCut(VarManager::kIsITSibAny, 0.5, 1.5);
3775+
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
3776+
cut->AddCut(VarManager::kTPCncls, 70, 161.);
3777+
return cut;
3778+
}
3779+
37663780
if (!nameStr.compare("electronStandardQualityITSOnly")) {
37673781
cut->AddCut(VarManager::kIsSPDany, 0.5, 1.5);
37683782
cut->AddCut(VarManager::kITSchi2, 0.0, 5.0);

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1009,6 +1009,21 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
10091009
double maxBins[6] = {5.0, 30.0, 4.0, 90.0, 10.0, 2.0};
10101010
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_V2", "", 6, varV2, bins, minBins, maxBins, 0, -1, kTRUE);
10111011
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_V3", "", 6, varV3, bins, minBins, maxBins, 0, -1, kTRUE);
1012+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2REF, VarManager::kM11REFoverMp);
1013+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4REF, VarManager::kM1111REFoverMp);
1014+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2POIMp, VarManager::kM01POIoverMp);
1015+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4POIMp, VarManager::kM0111POIoverMp);
1016+
}
1017+
if (subGroupStr.Contains("cumulant")) {
1018+
int var[4] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C};
1019+
int bins[4] = {250, 60, 6, 18};
1020+
double minBins[4] = {0.0, 0.0, 2.5, 0.0};
1021+
double maxBins[4] = {5.0, 30.0, 4.0, 90.0};
1022+
hm->AddHistogram(histClass, "Mass_Pt_Rapidity_CentFT0C", "", 4, var, bins, minBins, maxBins, 0, -1, kTRUE);
1023+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2REF, VarManager::kM11REFoverMp);
1024+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4REF", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4REF, VarManager::kM1111REFoverMp);
1025+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr2POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR2POIMp, VarManager::kM01POIoverMp);
1026+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_Corr4POI", "", true, 250, 0.0, 5.0, VarManager::kMass, 60, 0.0, 30.0, VarManager::kPt, 18, 0.0, 90.0, VarManager::kCentFT0C, "", "", "", VarManager::kCORR4POIMp, VarManager::kM0111POIoverMp);
10121027
}
10131028
if (subGroupStr.Contains("res-flow-dimuon")) {
10141029
int varV2[6] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kR2SP_AB, VarManager::kR2EP_AB};

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 16 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -739,26 +739,32 @@ void VarManager::SetDefaultVarNames()
739739
fgVariableUnits[kS31A] = "";
740740
fgVariableNames[kM11REF] = "M_{11}^{REF} ";
741741
fgVariableUnits[kM11REF] = "";
742-
fgVariableNames[kM01POI] = "M_{01}^{POI} ";
742+
fgVariableNames[kM01POI] = "M^{'}_{01}^{POI} ";
743743
fgVariableUnits[kM01POI] = "";
744744
fgVariableNames[kM1111REF] = "M_{1111}^{REF} ";
745745
fgVariableUnits[kM1111REF] = "";
746-
fgVariableNames[kM0111POI] = "M_{0111}^{POI} ";
746+
fgVariableNames[kM0111POI] = "M^{'}_{0111}^{POI} ";
747747
fgVariableUnits[kM0111POI] = "";
748748
fgVariableNames[kCORR2REF] = "<2> ";
749749
fgVariableUnits[kCORR2REF] = "";
750-
fgVariableNames[kCORR2POI] = " <2'> ";
750+
fgVariableNames[kCORR2POI] = "<2'> ";
751751
fgVariableUnits[kCORR2POI] = "";
752-
fgVariableNames[kCORR4REF] = " <4>";
752+
fgVariableNames[kCORR4REF] = "<4> ";
753753
fgVariableUnits[kCORR4REF] = "";
754754
fgVariableNames[kCORR4POI] = "<4'> ";
755755
fgVariableUnits[kCORR4POI] = "";
756-
fgVariableNames[kC4REF] = "c_{2}(4)";
757-
fgVariableUnits[kC4REF] = "";
758-
fgVariableNames[kC4POI] = "d_{2}(4)";
759-
fgVariableUnits[kC4POI] = "";
760-
fgVariableNames[kV4] = "v_{2}(4)";
761-
fgVariableUnits[kV4] = "";
756+
fgVariableNames[kM11REFoverMp] = "M_{11}^{REF}/M_{p} ";
757+
fgVariableUnits[kM11REFoverMp] = "";
758+
fgVariableNames[kM01POIoverMp] = "M^{'}_{01}^{POI}/M_{p} ";
759+
fgVariableUnits[kM01POIoverMp] = "";
760+
fgVariableNames[kM1111REFoverMp] = "M_{1111}^{REF}/M_{p} ";
761+
fgVariableUnits[kM1111REFoverMp] = "";
762+
fgVariableNames[kM0111POIoverMp] = "M^{'}_{0111}^{POI}/M_{p} ";
763+
fgVariableUnits[kM0111POIoverMp] = "";
764+
fgVariableNames[kCORR2POIMp] = "<2'> M_{p} ";
765+
fgVariableUnits[kCORR2POIMp] = "";
766+
fgVariableNames[kCORR4POIMp] = "<4'> M_{p} ";
767+
fgVariableUnits[kCORR4POIMp] = "";
762768
fgVariableNames[kCos2DeltaPhi] = "cos 2(#varphi-#Psi_{2}^{A}) ";
763769
fgVariableUnits[kCos2DeltaPhi] = "";
764770
fgVariableNames[kCos3DeltaPhi] = "cos 3(#varphi-#Psi_{3}^{A}) ";

‎PWGDQ/Core/VarManager.h‎

Lines changed: 14 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -601,9 +601,12 @@ class VarManager : public TObject
601601
kCORR2POI,
602602
kCORR4REF,
603603
kCORR4POI,
604-
kC4REF,
605-
kC4POI,
606-
kV4,
604+
kM11REFoverMp,
605+
kM01POIoverMp,
606+
kM1111REFoverMp,
607+
kM0111POIoverMp,
608+
kCORR2POIMp,
609+
kCORR4POIMp,
607610
kPsi2A,
608611
kPsi2APOS,
609612
kPsi2ANEG,
@@ -3914,6 +3917,14 @@ void VarManager::FillPairVn(T1 const& t1, T2 const& t2, float* values)
39143917
values[kM0111POI] = values[kMultDimuons] * (values[kS31A] - 3. * values[kS11A] * values[kS12A] + 2. * values[kS13A]);
39153918
values[kCORR2POI] = (P2 * conj(Q21)).real() / values[kM01POI];
39163919
values[kCORR4POI] = (P2 * Q21 * conj(Q21) * conj(Q21) - P2 * Q21 * conj(Q42) - 2. * values[kS12A] * P2 * conj(Q21) + 2. * P2 * conj(Q23)).real() / values[kM0111POI];
3920+
values[kM01POIoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM01POI]) || isinf(values[kM01POI])) ? values[kM01POI] / values[kMultDimuons] : 0;
3921+
values[kM0111POIoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM0111POI]) || isinf(values[kM0111POI])) ? values[kM0111POI] / values[kMultDimuons] : 0;
3922+
values[kM11REFoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM11REF]) || isinf(values[kM11REF])) ? values[kM11REF] / values[kMultDimuons] : 0;
3923+
values[kM1111REFoverMp] = values[kMultDimuons] > 0 && !(isnan(values[kM1111REF]) || isinf(values[kM1111REF])) ? values[kM1111REF] / values[kMultDimuons] : 0;
3924+
values[kCORR2POIMp] = isnan(values[kCORR2POI]) || isinf(values[kCORR2POI]) ? 0 : values[kCORR2POI] * values[kMultDimuons];
3925+
values[kCORR4POIMp] = isnan(values[kCORR4POI]) || isinf(values[kCORR4POI]) ? 0 : values[kCORR4POI] * values[kMultDimuons];
3926+
values[kCORR2REF] = isnan(values[kCORR2REF]) || isinf(values[kCORR2REF]) ? 0 : values[kCORR2REF];
3927+
values[kCORR4REF] = isnan(values[kCORR4REF]) || isinf(values[kCORR4REF]) ? 0 : values[kCORR4REF];
39173928
}
39183929
}
39193930

‎PWGEM/Dilepton/Core/SingleTrackQC.h‎

Lines changed: 13 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -101,8 +101,8 @@ struct SingleTrackQC {
101101
Configurable<float> cfg_max_chi2its{"cfg_max_chi2its", 5.0, "max chi2/NclsITS"};
102102
Configurable<float> cfg_max_dcaxy{"cfg_max_dcaxy", 1.0, "max dca XY for single track in cm"};
103103
Configurable<float> cfg_max_dcaz{"cfg_max_dcaz", 1.0, "max dca Z for single track in cm"};
104-
Configurable<bool> cfg_require_itsib_any{"cfg_require_itsib_any", true, "flag to require ITS ib any hits"};
105-
Configurable<bool> cfg_require_itsib_1st{"cfg_require_itsib_1st", false, "flag to require ITS ib 1st hit"};
104+
Configurable<bool> cfg_require_itsib_any{"cfg_require_itsib_any", false, "flag to require ITS ib any hits"};
105+
Configurable<bool> cfg_require_itsib_1st{"cfg_require_itsib_1st", true, "flag to require ITS ib 1st hit"};
106106

107107
Configurable<int> cfg_pid_scheme{"cfg_pid_scheme", static_cast<int>(DielectronCut::PIDSchemes::kTPChadrejORTOFreq), "pid scheme [kTOFreq : 0, kTPChadrej : 1, kTPChadrejORTOFreq : 2, kTPConly : 3]"};
108108
Configurable<float> cfg_min_TPCNsigmaEl{"cfg_min_TPCNsigmaEl", -2.0, "min. TPC n sigma for electron inclusion"};
@@ -176,17 +176,17 @@ struct SingleTrackQC {
176176
const AxisSpec axis_phi{36, 0.0, 2 * M_PI, "#varphi_{e} (rad.)"};
177177
const AxisSpec axis_pin{{0.00, 0.05, 0.06, 0.07, 0.08, 0.09, 0.10, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90, 1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.00}, "p_{TPC inner wall} (GeV/c)"};
178178

179-
const AxisSpec axis_tpc_nsigma_el{100, -5.f, +5.f, "n #sigma_{e}^{TPC}"};
180-
const AxisSpec axis_tpc_nsigma_mu{100, -5.f, +5.f, "n #sigma_{#mu}^{TPC}"};
181-
const AxisSpec axis_tpc_nsigma_pi{100, -5.f, +5.f, "n #sigma_{#pi}^{TPC}"};
182-
const AxisSpec axis_tpc_nsigma_ka{100, -5.f, +5.f, "n #sigma_{K}^{TPC}"};
183-
const AxisSpec axis_tpc_nsigma_pr{100, -5.f, +5.f, "n #sigma_{p}^{TPC}"};
184-
185-
const AxisSpec axis_tof_nsigma_el{100, -5.f, +5.f, "n #sigma_{e}^{TOF}"};
186-
const AxisSpec axis_tof_nsigma_mu{100, -5.f, +5.f, "n #sigma_{#mu}^{TOF}"};
187-
const AxisSpec axis_tof_nsigma_pi{100, -5.f, +5.f, "n #sigma_{#pi}^{TOF}"};
188-
const AxisSpec axis_tof_nsigma_ka{100, -5.f, +5.f, "n #sigma_{K}^{TOF}"};
189-
const AxisSpec axis_tof_nsigma_pr{100, -5.f, +5.f, "n #sigma_{p}^{TOF}"};
179+
const AxisSpec axis_tpc_nsigma_el{80, -4.f, +4.f, "n #sigma_{e}^{TPC}"};
180+
const AxisSpec axis_tpc_nsigma_mu{80, -4.f, +4.f, "n #sigma_{#mu}^{TPC}"};
181+
const AxisSpec axis_tpc_nsigma_pi{80, -4.f, +4.f, "n #sigma_{#pi}^{TPC}"};
182+
const AxisSpec axis_tpc_nsigma_ka{80, -4.f, +4.f, "n #sigma_{K}^{TPC}"};
183+
const AxisSpec axis_tpc_nsigma_pr{80, -4.f, +4.f, "n #sigma_{p}^{TPC}"};
184+
185+
const AxisSpec axis_tof_nsigma_el{80, -4.f, +4.f, "n #sigma_{e}^{TOF}"};
186+
const AxisSpec axis_tof_nsigma_mu{80, -4.f, +4.f, "n #sigma_{#mu}^{TOF}"};
187+
const AxisSpec axis_tof_nsigma_pi{80, -4.f, +4.f, "n #sigma_{#pi}^{TOF}"};
188+
const AxisSpec axis_tof_nsigma_ka{80, -4.f, +4.f, "n #sigma_{K}^{TOF}"};
189+
const AxisSpec axis_tof_nsigma_pr{80, -4.f, +4.f, "n #sigma_{p}^{TOF}"};
190190

191191
// track info
192192
fRegistry.add("Track/positive/hs", "rec. single electron", kTHnSparseF, {axis_pt, axis_eta, axis_phi}, true);

0 commit comments

Comments
 (0)