Skip to content

Commit 558feec

Browse files
axti98axti98
authored
[PWGDQ] Update ALICE 3 heavy-flavor cuts (#18603)
1 parent 586e49c commit 558feec

1 file changed

Lines changed: 128 additions & 32 deletions

File tree

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 128 additions & 32 deletions
Original file line numberDiff line numberDiff line change
@@ -3707,11 +3707,6 @@ 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-
37153710
if (nameStr == "pairD0HighPt1") {
37163711
cut->AddCut(GetAnalysisCut("pairLxyzProjected3sigma"));
37173712
cut->AddCut(GetAnalysisCut("pairPtLow5"));
@@ -3939,15 +3934,15 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
39393934
return cut;
39403935
}
39413936

3942-
if (nameStr == "alice3TOFEtaD0") {
3937+
if (nameStr == "alice3TOFEta") {
39433938
cut->AddCut(GetAnalysisCut("alice3TOFEta"));
3944-
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3939+
cut->AddCut(GetAnalysisCut("alice3TrackQualityHF"));
39453940
return cut;
39463941
}
39473942

3948-
if (nameStr == "alice3RICHEtaD0") {
3943+
if (nameStr == "alice3RICHEta") {
39493944
cut->AddCut(GetAnalysisCut("alice3RICHEta"));
3950-
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3945+
cut->AddCut(GetAnalysisCut("alice3TrackQualityHF"));
39513946
return cut;
39523947
}
39533948

@@ -4119,6 +4114,30 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
41194114
return cut;
41204115
}
41214116

4117+
if (nameStr == "alice3KaonTOFRICH") {
4118+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
4119+
cut->AddCut(GetAnalysisCut("alice3KaonTOFRICH"));
4120+
return cut;
4121+
}
4122+
4123+
if (nameStr == "alice3PionTOFRICH") {
4124+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
4125+
cut->AddCut(GetAnalysisCut("alice3PionTOFRICH"));
4126+
return cut;
4127+
}
4128+
4129+
if (nameStr == "alice3ProtonTOFRICH") {
4130+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
4131+
cut->AddCut(GetAnalysisCut("alice3ProtonTOFRICH"));
4132+
return cut;
4133+
}
4134+
4135+
if (nameStr == "alice3ElectronTOFRICH") {
4136+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
4137+
cut->AddCut(GetAnalysisCut("alice3ElectronTOFRICH"));
4138+
return cut;
4139+
}
4140+
41224141
if (nameStr == "alice3DielectronPID") {
41234142
cut->AddCut(GetAnalysisCut("alice3JpsiKine"));
41244143
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
@@ -7176,31 +7195,27 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
71767195
}
71777196

71787197
if (nameStr == "pairD0Topology") {
7179-
cut->AddCut(VarManager::kCosPointingAngle, 0.9, 1.0);
7180-
cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.03);
7181-
cut->AddCut(VarManager::kD0ImpParProduct, -1000.0, 0.0);
7182-
cut->AddCut(VarManager::kDCAxy1, -0.1, 0.1);
7183-
cut->AddCut(VarManager::kDCAxy2, -0.1, 0.1);
7184-
return cut;
7185-
}
7186-
7187-
if (nameStr == "pairD0TopologyTight") {
71887198
cut->AddCut(VarManager::kCosPointingAngle, 0.8, 1.0, false, VarManager::kPt, 0.0, 1.0);
71897199
cut->AddCut(VarManager::kCosPointingAngle, 0.9, 1.0, false, VarManager::kPt, 1.0, 4.0);
71907200
cut->AddCut(VarManager::kCosPointingAngle, 0.85, 1.0, false, VarManager::kPt, 4.0, 50.0);
7191-
cut->AddCut(VarManager::kD0DcaDaughters, 0.035, 1000.0, false, VarManager::kPt, 0.0, 1.0);
7192-
cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.03, false, VarManager::kPt, 1.0, 50.0);
7193-
cut->AddCut(VarManager::kD0ImpParProduct, -1000.0, 0.0);
7194-
cut->AddCut(VarManager::kDCAxy1, -0.1, 0.1);
7195-
cut->AddCut(VarManager::kDCAxy2, -0.1, 0.1);
7196-
cut->AddCut(VarManager::kPt1, 0.5, 100.0, false, VarManager::kPt, 0.0, 1.0);
7197-
cut->AddCut(VarManager::kPt1, 0.3, 100.0, false, VarManager::kPt, 1.0, 1.5);
7198-
cut->AddCut(VarManager::kPt1, 0.4, 100.0, false, VarManager::kPt, 1.5, 2.0);
7199-
cut->AddCut(VarManager::kPt1, 0.7, 100.0, false, VarManager::kPt, 2.0, 50.0);
7200-
cut->AddCut(VarManager::kPt2, 0.5, 100.0, false, VarManager::kPt, 0.0, 1.0);
7201-
cut->AddCut(VarManager::kPt2, 0.3, 100.0, false, VarManager::kPt, 1.0, 1.5);
7202-
cut->AddCut(VarManager::kPt2, 0.4, 100.0, false, VarManager::kPt, 1.5, 2.0);
7203-
cut->AddCut(VarManager::kPt2, 0.7, 100.0, false, VarManager::kPt, 2.0, 50.0);
7201+
cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.01, false, VarManager::kPt, 0.0, 50.0);
7202+
cut->AddCut(VarManager::kVertexingLxy, 0.01, 10.0);
7203+
cut->AddCut(VarManager::kRap, -2.0, 2.0);
7204+
return cut;
7205+
}
7206+
7207+
if (nameStr == "pairLambdaCTopology") {
7208+
cut->AddCut(VarManager::kPt, 0.0, 50.0);
7209+
cut->AddCut(VarManager::kRap, -2.0, 2.0);
7210+
cut->AddCut(VarManager::kVertexingProcCode, 0.5, 1.e6);
7211+
cut->AddCut(VarManager::kCosPointingAngle, 0.80, 1.0, false, VarManager::kPt, 0.0, 1.0);
7212+
cut->AddCut(VarManager::kCosPointingAngle, 0.90, 1.0, false, VarManager::kPt, 1.0, 4.0);
7213+
cut->AddCut(VarManager::kCosPointingAngle, 0.95, 1.0, false, VarManager::kPt, 4.0, 50.0);
7214+
cut->AddCut(VarManager::kVertexingLxy, 0.0, 10.0);
7215+
cut->AddCut(VarManager::kVertexingLxy, 0.002, 10.0, false, VarManager::kPt, 0.0, 1.0);
7216+
cut->AddCut(VarManager::kVertexingLxy, 0.003, 10.0, false, VarManager::kPt, 1.0, 4.0);
7217+
cut->AddCut(VarManager::kVertexingLxy, 0.005, 10.0, false, VarManager::kPt, 4.0, 50.0);
7218+
72047219
return cut;
72057220
}
72067221

@@ -7467,12 +7482,13 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
74677482
}
74687483

74697484
if (nameStr == "alice3TOFEta") {
7470-
cut->AddCut(VarManager::kPt, 0.7, 1000.0);
7485+
cut->AddCut(VarManager::kPt, 0.3, 1000.0);
74717486
cut->AddCut(VarManager::kEta, -2.0, 2.0); // Total tracker acceptance in v3b geomety
74727487
return cut;
74737488
}
74747489

74757490
if (nameStr == "alice3RICHEta") {
7491+
cut->AddCut(VarManager::kPt, 0.3, 1000.0);
74767492
cut->AddCut(VarManager::kEta, -0.8, 0.8); // Total tracker acceptance in v3b geomety
74777493
return cut;
74787494
}
@@ -7519,6 +7535,14 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
75197535
return cut;
75207536
}
75217537

7538+
if (nameStr == "alice3TrackQualityHF") {
7539+
cut->AddCut(VarManager::kIsReconstructed, 0.5, 1.5);
7540+
cut->AddCut(VarManager::kNSiliconHits, 5.0, 12.0);
7541+
cut->AddCut(VarManager::kTrackDCAxy, -0.05, 0.05);
7542+
cut->AddCut(VarManager::kTrackDCAz, -0.05, 0.05);
7543+
return cut;
7544+
}
7545+
75227546
if (nameStr == "alice3TrackQualityTightDCA") {
75237547
cut->AddCut(VarManager::kIsReconstructed, 0.5, 1.5);
75247548
cut->AddCut(VarManager::kNSiliconHits, 6.0, 12.0);
@@ -7643,6 +7667,78 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
76437667
return cut;
76447668
}
76457669

7670+
if (nameStr == "alice3KaonTOFRICH") {
7671+
cut->AddCut(VarManager::kEta, -2.0, 2.0);
7672+
cut->AddCut(VarManager::kP, 0.3, 1000.0);
7673+
7674+
cut->AddCut(VarManager::kTOFEventTimeErr, 0.0, 30.0, false, VarManager::kP, 0.0, 2.0);
7675+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0, false, VarManager::kP, 0.0f, 2.0);
7676+
7677+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -2.0, 2.0, true, VarManager::kP, 0.0, 1.9);
7678+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -2.0, 2.0, true, VarManager::kP, 0.0, 3.2);
7679+
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 2.0, true, VarManager::kP, 0.0, 2.0);
7680+
7681+
cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 2.0, 1000.0);
7682+
cut->AddCut(VarManager::kRICHnSigmaKa, 0.5, 1.5, false, VarManager::kP, 2.0, 1000.0);
7683+
cut->AddCut(VarManager::kRICHnSigmaKa, -3.0, 3.0, false, VarManager::kP, 2.0, 1000.0);
7684+
7685+
return cut;
7686+
}
7687+
7688+
if (nameStr == "alice3PionTOFRICH") {
7689+
cut->AddCut(VarManager::kEta, -2.0, 2.0);
7690+
cut->AddCut(VarManager::kP, 0.3, 1000.0);
7691+
7692+
cut->AddCut(VarManager::kTOFEventTimeErr, 0.0, 30.0, false, VarManager::kP, 0.0, 2.0);
7693+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0, false, VarManager::kP, 0.0, 2.0);
7694+
7695+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -2.0, 2.0, true, VarManager::kP, 0.0, 1.9);
7696+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -2.0, 2.0, true, VarManager::kP, 0.0, 3.8);
7697+
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 2.0, true, VarManager::kP, 0.0, 0.55);
7698+
7699+
cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 2.0, 1000.0);
7700+
cut->AddCut(VarManager::kRICHnSigmaPi, 0.5, 1.5, false, VarManager::kP, 2.0, 1000.0);
7701+
cut->AddCut(VarManager::kRICHnSigmaPi, -3.0, 3.0, false, VarManager::kP, 2.0, 1000.0);
7702+
7703+
return cut;
7704+
}
7705+
7706+
if (nameStr == "alice3ProtonTOFRICH") {
7707+
cut->AddCut(VarManager::kEta, -2.0, 2.0);
7708+
cut->AddCut(VarManager::kP, 0.3f, 1000.0);
7709+
7710+
cut->AddCut(VarManager::kTOFEventTimeErr, 0.0f, 30.0f, false, VarManager::kP, 0.0f, 4.5);
7711+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0f, 3.0f, false, VarManager::kP, 0.0f, 4.5);
7712+
7713+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -2.0, 2.0, true, VarManager::kP, 0.0, 3.2);
7714+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -2.0, 2.0, true, VarManager::kP, 0.0, 3.8);
7715+
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 2.0, true, VarManager::kP, 0.0, 3.8);
7716+
7717+
cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 4.5, 1000.0);
7718+
cut->AddCut(VarManager::kHasRICHSigPr, 0.5f, 1.5f, false, VarManager::kP, 4.5, 1000.0);
7719+
cut->AddCut(VarManager::kRICHnSigmaPr, -3.0f, 3.0f, false, VarManager::kP, 4.5, 1000.0);
7720+
7721+
return cut;
7722+
}
7723+
7724+
if (nameStr == "alice3ElectronTOFRICH") {
7725+
cut->AddCut(VarManager::kEta, -2.0, 2.0);
7726+
cut->AddCut(VarManager::kP, 0.3f, 1000.0);
7727+
7728+
cut->AddCut(VarManager::kTOFEventTimeErr, 0.0f, 30.0f, false, VarManager::kP, 0.0f, 0.6);
7729+
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -3.0f, 3.0f, false, VarManager::kP, 0.0f, 0.6);
7730+
7731+
cut->AddCut(VarManager::kOuterTOFnSigmaPi, -2.0f, 2.0f, true, VarManager::kP, 0.0f, 0.6);
7732+
cut->AddCut(VarManager::kOuterTOFnSigmaKa, -2.0f, 2.0f, true, VarManager::kP, 0.0f, 0.6);
7733+
cut->AddCut(VarManager::kOuterTOFnSigmaPr, -2.0f, 2.0f, true, VarManager::kP, 0.0f, 0.6);
7734+
7735+
cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 0.6, 1000.0);
7736+
cut->AddCut(VarManager::kHasRICHSigEl, 0.5f, 1.5f, false, VarManager::kP, 0.6, 1000.0);
7737+
cut->AddCut(VarManager::kRICHnSigmaEl, -3.0f, 3.0f, false, VarManager::kP, 0.6, 1000.0);
7738+
7739+
return cut;
7740+
}
7741+
76467742
if (nameStr == "alice3CharmoniumPID") {
76477743
cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 3.0, false, VarManager::kP, 0.0, 1.2);
76487744
cut->AddCut(VarManager::kRICHnSigmaEl, -2.0, 3.0, false, VarManager::kHasRICHSigEl, 0.5, 1.5, false, VarManager::kP, 0.7, 1000.0);

0 commit comments

Comments
 (0)