Skip to content

Commit 2522c1c

Browse files
axti98axti98
committed
Updated ALICE 3 concerned cuts for D0
1 parent e7809d1 commit 2522c1c

3 files changed

Lines changed: 129 additions & 11 deletions

File tree

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 112 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -3707,6 +3707,11 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
37073707
return cut;
37083708
}
37093709

3710+
if (nameStr == "pairD0TopologyTight") {
3711+
cut->AddCut(GetAnalysisCut("pairD0TopologyTight"));
3712+
return cut;
3713+
}
3714+
37103715
if (nameStr == "pairD0HighPt1") {
37113716
cut->AddCut(GetAnalysisCut("pairLxyzProjected3sigma"));
37123717
cut->AddCut(GetAnalysisCut("pairPtLow5"));
@@ -3934,6 +3939,18 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
39343939
return cut;
39353940
}
39363941

3942+
if (nameStr == "alice3TOFEtaD0") {
3943+
cut->AddCut(GetAnalysisCut("alice3TOFEta"));
3944+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3945+
return cut;
3946+
}
3947+
3948+
if (nameStr == "alice3RICHEtaD0") {
3949+
cut->AddCut(GetAnalysisCut("alice3RICHEta"));
3950+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3951+
return cut;
3952+
}
3953+
39373954
if (nameStr == "alice3StandardTrackTOFAcceptanceD0") {
39383955
cut->AddCut(GetAnalysisCut("alice3KineTOFAcceptanceD0"));
39393956
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
@@ -4028,6 +4045,32 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
40284045
return cut;
40294046
}
40304047

4048+
if (nameStr == "alice3EitherTOFPIDPi") {
4049+
auto* pidCut = new AnalysisCompositeCut(
4050+
"alice3EitherTOFPIDPi",
4051+
"iTOF OR oTOF pion PID",
4052+
kFALSE);
4053+
4054+
pidCut->AddCut(GetAnalysisCut("alice3iTOFPIDPi"));
4055+
pidCut->AddCut(GetAnalysisCut("alice3oTOFPIDPi"));
4056+
4057+
delete cut;
4058+
return pidCut;
4059+
}
4060+
4061+
if (nameStr == "alice3EitherTOFPIDKa") {
4062+
auto* pidCut = new AnalysisCompositeCut(
4063+
"alice3EitherTOFPIDKa",
4064+
"iTOF OR oTOF kaon PID",
4065+
kFALSE);
4066+
4067+
pidCut->AddCut(GetAnalysisCut("alice3iTOFPIDKa"));
4068+
pidCut->AddCut(GetAnalysisCut("alice3oTOFPIDKa"));
4069+
4070+
delete cut;
4071+
return pidCut;
4072+
}
4073+
40314074
if (nameStr == "alice3RICHPIDEl") {
40324075
cut->AddCut(GetAnalysisCut("alice3RICHPIDEl"));
40334076
return cut;
@@ -4137,9 +4180,13 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
41374180
return cut;
41384181
}
41394182

4140-
if (nameStr == "alice3TOFPIDKa") {
4141-
cut->AddCut(GetAnalysisCut("alice3iTOFPIDKa"));
4142-
cut->AddCut(GetAnalysisCut("alice3oTOFPIDKa"));
4183+
if (nameStr == "alice3TOFPIDKaRej") {
4184+
cut->AddCut(GetAnalysisCut("alice3D0KaonTOFPID"));
4185+
return cut;
4186+
}
4187+
4188+
if (nameStr == "alice3TOFPIDPiRej") {
4189+
cut->AddCut(GetAnalysisCut("alice3D0PionTOFPID"));
41434190
return cut;
41444191
}
41454192

@@ -7127,6 +7174,26 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
71277174
return cut;
71287175
}
71297176

7177+
if (nameStr == "pairD0TopologyTight") {
7178+
cut->AddCut(VarManager::kCosPointingAngle, 0.8, 1.0, false, VarManager::kPt, 0.0, 1.0);
7179+
cut->AddCut(VarManager::kCosPointingAngle, 0.9, 1.0, false, VarManager::kPt, 1.0, 4.0);
7180+
cut->AddCut(VarManager::kCosPointingAngle, 0.85, 1.0, false, VarManager::kPt, 4.0, 50.0);
7181+
cut->AddCut(VarManager::kD0DcaDaughters, 0.035, 1000.0, false, VarManager::kPt, 0.0, 1.0);
7182+
cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.03, false, VarManager::kPt, 1.0, 50.0);
7183+
cut->AddCut(VarManager::kD0ImpParProduct, -1000.0, 0.0);
7184+
cut->AddCut(VarManager::kDCAxy1, -0.1, 0.1);
7185+
cut->AddCut(VarManager::kDCAxy2, -0.1, 0.1);
7186+
cut->AddCut(VarManager::kPt1, 0.5, 100.0, false, VarManager::kPt, 0.0, 1.0);
7187+
cut->AddCut(VarManager::kPt1, 0.3, 100.0, false, VarManager::kPt, 1.0, 1.5);
7188+
cut->AddCut(VarManager::kPt1, 0.4, 100.0, false, VarManager::kPt, 1.5, 2.0);
7189+
cut->AddCut(VarManager::kPt1, 0.7, 100.0, false, VarManager::kPt, 2.0, 50.0);
7190+
cut->AddCut(VarManager::kPt2, 0.5, 100.0, false, VarManager::kPt, 0.0, 1.0);
7191+
cut->AddCut(VarManager::kPt2, 0.3, 100.0, false, VarManager::kPt, 1.0, 1.5);
7192+
cut->AddCut(VarManager::kPt2, 0.4, 100.0, false, VarManager::kPt, 1.5, 2.0);
7193+
cut->AddCut(VarManager::kPt2, 0.7, 100.0, false, VarManager::kPt, 2.0, 50.0);
7194+
return cut;
7195+
}
7196+
71307197
if (nameStr == "pairJpsi") {
71317198
cut->AddCut(VarManager::kMass, 2.8, 3.3);
71327199
return cut;
@@ -7389,6 +7456,16 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
73897456
return cut;
73907457
}
73917458

7459+
if (nameStr == "alice3TOFEta") {
7460+
cut->AddCut(VarManager::kEta, -2.0, 2.0); // Total tracker acceptance in v3b geomety
7461+
return cut;
7462+
}
7463+
7464+
if (nameStr == "alice3RICHEta") {
7465+
cut->AddCut(VarManager::kEta, -0.8, 0.8); // Total tracker acceptance in v3b geomety
7466+
return cut;
7467+
}
7468+
73927469
if (nameStr == "alice3KineTOFAcceptanceD0") {
73937470
cut->AddCut(VarManager::kPt, 0.7, 1000.0);
73947471
cut->AddCut(VarManager::kEta, -2.0, 2.0); // TOF acceptance in v3b geomety
@@ -7460,42 +7537,66 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
74607537
}
74617538

74627539
if (nameStr == "alice3iTOFPIDEl") {
7463-
cut->AddCut(VarManager::kInnerTOFnSigmaEl, -3.0, 3.0);
7540+
cut->AddCut(VarManager::kInnerTOFnSigmaEl, -3.0, 3.0, false, VarManager::kP, 0.0, 1.5);
74647541
return cut;
74657542
}
74667543

74677544
if (nameStr == "alice3iTOFPIDPi") {
7468-
cut->AddCut(VarManager::kInnerTOFnSigmaPi, -3.0, 3.0);
7545+
cut->AddCut(VarManager::kInnerTOFnSigmaPi, -3.0, 3.0, false, VarManager::kP, 0.0, 1.5);
74697546
return cut;
74707547
}
74717548

74727549
if (nameStr == "alice3iTOFPIDKa") {
7473-
cut->AddCut(VarManager::kInnerTOFnSigmaKa, -3.0, 3.0);
7550+
cut->AddCut(VarManager::kInnerTOFnSigmaKa, -3.0, 3.0, false, VarManager::kP, 0.0, 1.5);
74747551
return cut;
74757552
}
74767553

74777554
if (nameStr == "alice3iTOFPIDPr") {
7478-
cut->AddCut(VarManager::kInnerTOFnSigmaPr, -3.0, 3.0);
7555+
cut->AddCut(VarManager::kInnerTOFnSigmaPr, -3.0, 3.0, false, VarManager::kP, 0.0, 1.5);
74797556
return cut;
74807557
}
74817558

74827559
if (nameStr == "alice3oTOFPIDEl") {
7483-
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -3.0, 3.0);
7560+
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -3.0, 3.0, false, VarManager::kP, 0.0, 2.5);
74847561
return cut;
74857562
}
74867563

74877564
if (nameStr == "alice3oTOFPIDPi") {
7488-
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0);
7565+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0, false, VarManager::kP, 0.0, 2.5);
74897566
return cut;
74907567
}
74917568

74927569
if (nameStr == "alice3oTOFPIDKa") {
7493-
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0);
7570+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0, false, VarManager::kP, 0.0, 2.5);
74947571
return cut;
74957572
}
74967573

74977574
if (nameStr == "alice3oTOFPIDPr") {
7498-
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0, 3.0);
7575+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0, 3.0, false, VarManager::kP, 0.0, 2.5);
7576+
return cut;
7577+
}
7578+
7579+
if (nameStr == "alice3D0KaonTOFPID") {
7580+
cut->AddCut(VarManager::kInnerTOFnSigmaKa, -3.0, 3.0, false, VarManager::kP, 0.0, 1.5);
7581+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0, false, VarManager::kP, 0.0, 2.5);
7582+
7583+
cut->AddCut(VarManager::kInnerTOFnSigmaPi, -3.0, 3.0, true, VarManager::kP, 0.0, 1.5);
7584+
cut->AddCut(VarManager::kInnerTOFnSigmaPr, -3.0, 3.0, true, VarManager::kP, 0.0, 1.5);
7585+
7586+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0, true, VarManager::kP, 0.0, 2.5);
7587+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0, 3.0, true, VarManager::kP, 0.0, 2.5);
7588+
return cut;
7589+
}
7590+
7591+
if (nameStr == "alice3D0PionTOFPID") {
7592+
cut->AddCut(VarManager::kInnerTOFnSigmaPi, -3.0, 3.0, false, VarManager::kP, 0.0, 1.5);
7593+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0, false, VarManager::kP, 0.0, 2.5);
7594+
7595+
cut->AddCut(VarManager::kInnerTOFnSigmaKa, -3.0, 3.0, true, VarManager::kP, 0.0, 1.5);
7596+
cut->AddCut(VarManager::kInnerTOFnSigmaPr, -3.0, 3.0, true, VarManager::kP, 0.0, 1.5);
7597+
7598+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0, true, VarManager::kP, 0.0, 2.5);
7599+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0, 3.0, true, VarManager::kP, 0.0, 2.5);
74997600
return cut;
75007601
}
75017602

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1403,6 +1403,13 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
14031403
}
14041404
hm->AddHistogram(histClass, "Mass_PtLong", "", false, 500, alice3MassBins.data(), VarManager::kMass, 20, pTBins.data(), VarManager::kPt);
14051405
}
1406+
if (subGroupStr.Contains("alice3pair")) {
1407+
hm->AddHistogram(histClass, "CosinePointingAngle", "cos(#theta_{p})", false, 2000.0, -1.0, 1.0, VarManager::kCosPointingAngle);
1408+
hm->AddHistogram(histClass, "DCADaughters", "DCA_{dau}", false, 2000.0, -10.0, 10.0, VarManager::kD0DcaDaughters);
1409+
hm->AddHistogram(histClass, "ImpParProduct", "d_{k} #times d_{#pi}", false, 2000.0, -100.0, 100.0, VarManager::kD0ImpParProduct);
1410+
hm->AddHistogram(histClass, "DCAxy1", "DCA_{xy, 1}", false, 2000.0, -10.0, 10.0, VarManager::kDCAxy1);
1411+
hm->AddHistogram(histClass, "DCAxy2", "DCA_{xy, 2}", false, 2000.0, -10.0, 10.0, VarManager::kDCAxy2);
1412+
}
14061413
if (subGroupStr.Contains("alice3pbpb")) {
14071414
hm->AddHistogram(histClass, "Mass_Cent", "", false, 125, 0.0, 5.0, VarManager::kMass, 20, 0.0, 100.0, VarManager::kCent);
14081415
hm->AddHistogram(histClass, "Pt_CentFT0C", "", false, 100, 0.0, 10.0, VarManager::kPt, 20, 0.0, 100.0, VarManager::kCent);

‎PWGDQ/Core/MCSignalLibrary.cxx‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1754,6 +1754,16 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
17541754
signal = new MCSignal(name, "Kaons from D+/- decays", {prong}, {-1});
17551755
return signal;
17561756
}
1757+
if (nameStr == "KFromD0") {
1758+
MCProng prong(2, {321, Pdg::kD0}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false});
1759+
signal = new MCSignal(name, "Kaons from D0 decays", {prong}, {-1});
1760+
return signal;
1761+
}
1762+
if (nameStr == "PiFromD0") {
1763+
MCProng prong(2, {211, Pdg::kD0}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false});
1764+
signal = new MCSignal(name, "Pions from D0 decays", {prong}, {-1});
1765+
return signal;
1766+
}
17571767
if (nameStr == "LambdaC") {
17581768
MCProng prong(1, {Pdg::kLambdaCPlus}, {true}, {false}, {0}, {0}, {false});
17591769
signal = new MCSignal(name, "Lambda_c", {prong}, {-1});

0 commit comments

Comments
 (0)