Skip to content

Commit b075c38

Browse files
shreyasiacharyashreyasiacharya
authored
[PWGDQ] Add Pole-mass to pseudoproper decay time and TOF-PID to cut library (#9881)
1 parent 56f2e4d commit b075c38

File tree

6 files changed

+27
-6
lines changed

6 files changed

+27
-6
lines changed

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 14 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -153,7 +153,13 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
153153
cut->AddCut(GetAnalysisCut("electronPIDnsigmaMedium"));
154154
return cut;
155155
}
156-
156+
if (!nameStr.compare("electronSelection1_ionut_withTOFPID")) {
157+
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
158+
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
159+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
160+
cut->AddCut(GetAnalysisCut("electronPIDnsigmaMedium_withLargeTOFPID"));
161+
return cut;
162+
}
157163
if (!nameStr.compare("electronSelection1_idstoreh")) { // same as electronSelection1_ionut, but with kIsSPDAny -> kIsITSibAny
158164
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
159165
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug4"));
@@ -5406,7 +5412,13 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
54065412
cut->AddCut(VarManager::kTPCnSigmaPi, 2.7, 3000.0);
54075413
return cut;
54085414
}
5409-
5415+
if (!nameStr.compare("electronPIDnsigmaMedium_withLargeTOFPID")) {
5416+
cut->AddCut(VarManager::kTPCnSigmaEl, -3.0, 3.0);
5417+
cut->AddCut(VarManager::kTPCnSigmaPr, 2.7, 3000.0);
5418+
cut->AddCut(VarManager::kTPCnSigmaPi, 2.7, 3000.0);
5419+
cut->AddCut(VarManager::kTOFnSigmaEl, -5.0, 5.0);
5420+
return cut;
5421+
}
54105422
if (!nameStr.compare("electronPIDnsigmaSkewed")) {
54115423
cut->AddCut(VarManager::kTPCnSigmaEl, -2.0, 3.0);
54125424
cut->AddCut(VarManager::kTPCnSigmaPr, 3.5, 3000.0);

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -912,6 +912,7 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
912912
}
913913
hm->AddHistogram(histClass, "Mass_PtFine", "", false, 75, massBins, VarManager::kMass, 69, ptBins, VarManager::kPt);
914914
hm->AddHistogram(histClass, "Eta_Pt", "", false, 40, -2.0, 2.0, VarManager::kEta, 40, 0.0, 20.0, VarManager::kPt);
915+
hm->AddHistogram(histClass, "Y_Pt", "", false, 40, -2.0, 2.0, VarManager::kRap, 40, 0.0, 20.0, VarManager::kPt);
915916
hm->AddHistogram(histClass, "Mass_VtxZ", "", true, 30, -15.0, 15.0, VarManager::kVtxZ, 500, 0.0, 5.0, VarManager::kMass);
916917
if (subGroupStr.Contains("pbpb")) {
917918
hm->AddHistogram(histClass, "Mass_CentFT0C", "", false, 125, 0.0, 5.0, VarManager::kMass, 20, 0.0, 100.0, VarManager::kCentFT0C);

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -639,6 +639,8 @@ void VarManager::SetDefaultVarNames()
639639
fgVariableUnits[kVertexingTauzProjected] = "ns";
640640
fgVariableNames[kVertexingTauxyProjected] = "Pair pseudo-proper Tauxy";
641641
fgVariableUnits[kVertexingTauxyProjected] = "ns";
642+
fgVariableNames[kVertexingTauxyProjectedPoleJPsiMass] = "Pair pseudo-proper Tauxy (with pole JPsi mass)";
643+
fgVariableUnits[kVertexingTauxyProjectedPoleJPsiMass] = "ns";
642644
fgVariableNames[kVertexingTauxyzProjected] = "Pair pseudo-proper Tauxyz";
643645
fgVariableUnits[kVertexingTauxyzProjected] = "ns";
644646
fgVariableNames[kCosPointingAngle] = "cos(#theta_{pointing})";
@@ -1521,6 +1523,7 @@ void VarManager::SetDefaultVarNames()
15211523
fgVarNamesMap["kVertexingLxyzProjected"] = kVertexingLxyzProjected;
15221524
fgVarNamesMap["kVertexingTauzProjected"] = kVertexingTauzProjected;
15231525
fgVarNamesMap["kVertexingTauxyProjected"] = kVertexingTauxyProjected;
1526+
fgVarNamesMap["kVertexingTauxyProjectedPoleJPsiMass"] = kVertexingTauxyProjectedPoleJPsiMass;
15241527
fgVarNamesMap["kVertexingTauxyProjectedNs"] = kVertexingTauxyProjectedNs;
15251528
fgVarNamesMap["kVertexingTauxyzProjected"] = kVertexingTauxyzProjected;
15261529
fgVarNamesMap["kVertexingTauz"] = kVertexingTauz;

‎PWGDQ/Core/VarManager.h‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -635,6 +635,7 @@ class VarManager : public TObject
635635
kVertexingLxyzProjected,
636636
kVertexingTauzProjected,
637637
kVertexingTauxyProjected,
638+
kVertexingTauxyProjectedPoleJPsiMass,
638639
kVertexingTauxyProjectedNs,
639640
kVertexingTauxyzProjected,
640641
kVertexingTauz,
@@ -3406,6 +3407,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
34063407
values[kVertexingLxyzProjected] = ((secondaryVertex[0] - collision.posX()) * v12.Px()) + ((secondaryVertex[1] - collision.posY()) * v12.Py()) + ((secondaryVertex[2] - collision.posZ()) * v12.Pz());
34073408
values[kVertexingLxyzProjected] = values[kVertexingLxyzProjected] / TMath::Sqrt((v12.Px() * v12.Px()) + (v12.Py() * v12.Py()) + (v12.Pz() * v12.Pz()));
34083409
values[kVertexingTauxyProjected] = values[kVertexingLxyProjected] * v12.M() / (v12.Pt());
3410+
values[kVertexingTauxyProjectedPoleJPsiMass] = values[kVertexingLxyProjected] * o2::constants::physics::MassJPsi / (v12.Pt());
34093411
values[kVertexingTauxyProjectedNs] = values[kVertexingTauxyProjected] / o2::constants::physics::LightSpeedCm2NS;
34103412
values[kVertexingTauzProjected] = values[kVertexingLzProjected] * v12.M() / TMath::Abs(v12.Pz());
34113413
values[kVertexingTauxyzProjected] = values[kVertexingLxyzProjected] * v12.M() / (v12.P());
@@ -3494,6 +3496,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
34943496
values[kVertexingLxyzProjected] = (dxPair2PV * KFGeoTwoProng.GetPx()) + (dyPair2PV * KFGeoTwoProng.GetPy()) + (dzPair2PV * KFGeoTwoProng.GetPz());
34953497
values[kVertexingLxyzProjected] = values[kVertexingLxyzProjected] / TMath::Sqrt((KFGeoTwoProng.GetPx() * KFGeoTwoProng.GetPx()) + (KFGeoTwoProng.GetPy() * KFGeoTwoProng.GetPy()) + (KFGeoTwoProng.GetPz() * KFGeoTwoProng.GetPz()));
34963498
values[kVertexingTauxyProjected] = values[kVertexingLxyProjected] * KFGeoTwoProng.GetMass() / (KFGeoTwoProng.GetPt());
3499+
values[kVertexingTauxyProjectedPoleJPsiMass] = values[kVertexingLxyProjected] * o2::constants::physics::MassJPsi / (KFGeoTwoProng.GetPt());
34973500
values[kVertexingTauxyProjectedNs] = values[kVertexingTauxyProjected] / o2::constants::physics::LightSpeedCm2NS;
34983501
values[kVertexingTauzProjected] = values[kVertexingLzProjected] * KFGeoTwoProng.GetMass() / TMath::Abs(KFGeoTwoProng.GetPz());
34993502
}

‎PWGDQ/Tasks/dqEfficiency_withAssoc.cxx‎

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -87,6 +87,7 @@ DECLARE_SOA_COLUMN(McFlag, mcFlag, int8_t);
8787
DECLARE_SOA_COLUMN(Massee, massee, float);
8888
DECLARE_SOA_COLUMN(Ptee, ptee, float);
8989
DECLARE_SOA_COLUMN(Lxyee, lxyee, float);
90+
DECLARE_SOA_COLUMN(LxyeePoleMass, lxyeepolemass, float);
9091
DECLARE_SOA_COLUMN(Lzee, lzee, float);
9192
DECLARE_SOA_COLUMN(AmbiguousInBunchPairs, AmbiguousJpsiPairsInBunch, bool);
9293
DECLARE_SOA_COLUMN(AmbiguousOutOfBunchPairs, AmbiguousJpsiPairsOutOfBunch, bool);
@@ -100,7 +101,7 @@ DECLARE_SOA_TABLE(MuonTrackCuts, "AOD", "DQANAMUONCUTS", dqanalysisflags::IsMuon
100101
DECLARE_SOA_TABLE(MuonAmbiguities, "AOD", "DQMUONAMB", dqanalysisflags::MuonAmbiguityInBunch, dqanalysisflags::MuonAmbiguityOutOfBunch); //! joinable to ReducedMuonTracks
101102
DECLARE_SOA_TABLE(Prefilter, "AOD", "DQPREFILTER", dqanalysisflags::IsBarrelSelectedPrefilter); //! joinable to ReducedTracksAssoc
102103
DECLARE_SOA_TABLE(BmesonCandidates, "AOD", "DQBMESONS", dqanalysisflags::massBcandidate, dqanalysisflags::deltaMassBcandidate, dqanalysisflags::pTBcandidate, dqanalysisflags::LxyBcandidate, dqanalysisflags::LxyzBcandidate, dqanalysisflags::LzBcandidate, dqanalysisflags::TauxyBcandidate, dqanalysisflags::TauzBcandidate, dqanalysisflags::DCAxyzBetweenProngs, dqanalysisflags::CosPBcandidate, dqanalysisflags::Chi2Bcandidate, dqanalysisflags::McFlag);
103-
DECLARE_SOA_TABLE(JPsieeCandidates, "AOD", "DQPSEUDOPROPER", dqanalysisflags::Massee, dqanalysisflags::Ptee, dqanalysisflags::Lxyee, dqanalysisflags::Lzee, dqanalysisflags::AmbiguousInBunchPairs, dqanalysisflags::AmbiguousOutOfBunchPairs, dqanalysisflags::Corrassoc);
104+
DECLARE_SOA_TABLE(JPsieeCandidates, "AOD", "DQPSEUDOPROPER", dqanalysisflags::Massee, dqanalysisflags::Ptee, dqanalysisflags::Lxyee, dqanalysisflags::LxyeePoleMass, dqanalysisflags::Lzee, dqanalysisflags::AmbiguousInBunchPairs, dqanalysisflags::AmbiguousOutOfBunchPairs, dqanalysisflags::Corrassoc);
104105
} // namespace o2::aod
105106

106107
// Declarations of various short names
@@ -1839,7 +1840,7 @@ struct AnalysisSameEventPairing {
18391840
fHistMan->FillHistClass(histNames[icut][0].Data(), VarManager::fgValues); // reconstructed, unmatched
18401841
for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { // loop over MC signals
18411842
if (mcDecision & (static_cast<uint32_t>(1) << isig)) {
1842-
PromptNonPromptSepTable(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt], VarManager::fgValues[VarManager::kVertexingTauxyProjected], VarManager::fgValues[VarManager::kVertexingTauzProjected], isAmbiInBunch, isAmbiOutOfBunch, isCorrect_pair);
1843+
PromptNonPromptSepTable(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt], VarManager::fgValues[VarManager::kVertexingTauxyProjected], VarManager::fgValues[VarManager::kVertexingTauxyProjectedPoleJPsiMass], VarManager::fgValues[VarManager::kVertexingTauzProjected], isAmbiInBunch, isAmbiOutOfBunch, isCorrect_pair);
18431844
fHistMan->FillHistClass(histNamesMC[icut * fRecMCSignals.size() + isig][0].Data(), VarManager::fgValues); // matched signal
18441845
if (fConfigQA) {
18451846
if (isCorrectAssoc_leg1 && isCorrectAssoc_leg2) { // correct track-collision association

‎PWGDQ/Tasks/tableReader_withAssoc.cxx‎

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -93,6 +93,7 @@ DECLARE_SOA_COLUMN(Chi2Bcandidate, chi2Bcandidate, float);
9393
DECLARE_SOA_COLUMN(Massee, massJPsi2ee, float);
9494
DECLARE_SOA_COLUMN(Ptee, ptJPsi2ee, float);
9595
DECLARE_SOA_COLUMN(Lxyee, lxyJPsi2ee, float);
96+
DECLARE_SOA_COLUMN(LxyeePoleMass, lxyJPsi2eePoleMass, float);
9697
DECLARE_SOA_COLUMN(Lzee, lzJPsi2ee, float);
9798
DECLARE_SOA_COLUMN(AmbiguousInBunchPairs, AmbiguousJpsiPairsInBunch, bool);
9899
DECLARE_SOA_COLUMN(AmbiguousOutOfBunchPairs, AmbiguousJpsiPairsOutOfBunch, bool);
@@ -106,7 +107,7 @@ DECLARE_SOA_TABLE(MuonTrackCuts, "AOD", "DQANAMUONCUTSA", dqanalysisflags::IsMuo
106107
DECLARE_SOA_TABLE(MuonAmbiguities, "AOD", "DQMUONAMBA", dqanalysisflags::MuonAmbiguityInBunch, dqanalysisflags::MuonAmbiguityOutOfBunch); //! joinable to ReducedMuonTracks
107108
DECLARE_SOA_TABLE(Prefilter, "AOD", "DQPREFILTERA", dqanalysisflags::IsBarrelSelectedPrefilter); //! joinable to ReducedTracksAssoc
108109
DECLARE_SOA_TABLE(BmesonCandidates, "AOD", "DQBMESONSA", dqanalysisflags::massBcandidate, dqanalysisflags::deltamassBcandidate, dqanalysisflags::pTBcandidate, dqanalysisflags::LxyBcandidate, dqanalysisflags::LxyzBcandidate, dqanalysisflags::LzBcandidate, dqanalysisflags::TauxyBcandidate, dqanalysisflags::TauzBcandidate, dqanalysisflags::CosPBcandidate, dqanalysisflags::Chi2Bcandidate);
109-
DECLARE_SOA_TABLE(JPsieeCandidates, "AOD", "DQPSEUDOPROPER", dqanalysisflags::Massee, dqanalysisflags::Ptee, dqanalysisflags::Lxyee, dqanalysisflags::Lzee, dqanalysisflags::AmbiguousInBunchPairs, dqanalysisflags::AmbiguousOutOfBunchPairs);
110+
DECLARE_SOA_TABLE(JPsieeCandidates, "AOD", "DQPSEUDOPROPER", dqanalysisflags::Massee, dqanalysisflags::Ptee, dqanalysisflags::Lxyee, dqanalysisflags::LxyeePoleMass, dqanalysisflags::Lzee, dqanalysisflags::AmbiguousInBunchPairs, dqanalysisflags::AmbiguousOutOfBunchPairs);
110111
} // namespace o2::aod
111112

112113
// Declarations of various short names
@@ -1703,7 +1704,7 @@ struct AnalysisSameEventPairing {
17031704
isUnambiguous = !(isAmbiInBunch || isAmbiOutOfBunch);
17041705
if (sign1 * sign2 < 0) {
17051706
fHistMan->FillHistClass(histNames[icut][0].Data(), VarManager::fgValues);
1706-
PromptNonPromptSepTable(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt], VarManager::fgValues[VarManager::kVertexingTauxyProjected], VarManager::fgValues[VarManager::kVertexingTauzProjected], isAmbiInBunch, isAmbiOutOfBunch);
1707+
PromptNonPromptSepTable(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt], VarManager::fgValues[VarManager::kVertexingTauxyProjected], VarManager::fgValues[VarManager::kVertexingTauxyProjectedPoleJPsiMass], VarManager::fgValues[VarManager::kVertexingTauzProjected], isAmbiInBunch, isAmbiOutOfBunch);
17071708
if (isAmbiInBunch) {
17081709
fHistMan->FillHistClass(histNames[icut][3 + histIdxOffset].Data(), VarManager::fgValues);
17091710
}

0 commit comments

Comments
 (0)