diff --git a/PWGDQ/Core/CutsLibrary.cxx b/PWGDQ/Core/CutsLibrary.cxx index 467b4888b0e..8c617a96e59 100644 --- a/PWGDQ/Core/CutsLibrary.cxx +++ b/PWGDQ/Core/CutsLibrary.cxx @@ -3707,11 +3707,6 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName) return cut; } - if (nameStr == "pairD0TopologyTight") { - cut->AddCut(GetAnalysisCut("pairD0TopologyTight")); - return cut; - } - if (nameStr == "pairD0HighPt1") { cut->AddCut(GetAnalysisCut("pairLxyzProjected3sigma")); cut->AddCut(GetAnalysisCut("pairPtLow5")); @@ -3939,15 +3934,15 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName) return cut; } - if (nameStr == "alice3TOFEtaD0") { + if (nameStr == "alice3TOFEta") { cut->AddCut(GetAnalysisCut("alice3TOFEta")); - cut->AddCut(GetAnalysisCut("alice3TrackQuality")); + cut->AddCut(GetAnalysisCut("alice3TrackQualityHF")); return cut; } - if (nameStr == "alice3RICHEtaD0") { + if (nameStr == "alice3RICHEta") { cut->AddCut(GetAnalysisCut("alice3RICHEta")); - cut->AddCut(GetAnalysisCut("alice3TrackQuality")); + cut->AddCut(GetAnalysisCut("alice3TrackQualityHF")); return cut; } @@ -4119,6 +4114,30 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName) return cut; } + if (nameStr == "alice3KaonTOFRICH") { + cut->AddCut(GetAnalysisCut("alice3TrackQuality")); + cut->AddCut(GetAnalysisCut("alice3KaonTOFRICH")); + return cut; + } + + if (nameStr == "alice3PionTOFRICH") { + cut->AddCut(GetAnalysisCut("alice3TrackQuality")); + cut->AddCut(GetAnalysisCut("alice3PionTOFRICH")); + return cut; + } + + if (nameStr == "alice3ProtonTOFRICH") { + cut->AddCut(GetAnalysisCut("alice3TrackQuality")); + cut->AddCut(GetAnalysisCut("alice3ProtonTOFRICH")); + return cut; + } + + if (nameStr == "alice3ElectronTOFRICH") { + cut->AddCut(GetAnalysisCut("alice3TrackQuality")); + cut->AddCut(GetAnalysisCut("alice3ElectronTOFRICH")); + return cut; + } + if (nameStr == "alice3DielectronPID") { cut->AddCut(GetAnalysisCut("alice3JpsiKine")); cut->AddCut(GetAnalysisCut("alice3TrackQuality")); @@ -7176,31 +7195,27 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName) } if (nameStr == "pairD0Topology") { - cut->AddCut(VarManager::kCosPointingAngle, 0.9, 1.0); - cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.03); - cut->AddCut(VarManager::kD0ImpParProduct, -1000.0, 0.0); - cut->AddCut(VarManager::kDCAxy1, -0.1, 0.1); - cut->AddCut(VarManager::kDCAxy2, -0.1, 0.1); - return cut; - } - - if (nameStr == "pairD0TopologyTight") { cut->AddCut(VarManager::kCosPointingAngle, 0.8, 1.0, false, VarManager::kPt, 0.0, 1.0); cut->AddCut(VarManager::kCosPointingAngle, 0.9, 1.0, false, VarManager::kPt, 1.0, 4.0); cut->AddCut(VarManager::kCosPointingAngle, 0.85, 1.0, false, VarManager::kPt, 4.0, 50.0); - cut->AddCut(VarManager::kD0DcaDaughters, 0.035, 1000.0, false, VarManager::kPt, 0.0, 1.0); - cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.03, false, VarManager::kPt, 1.0, 50.0); - cut->AddCut(VarManager::kD0ImpParProduct, -1000.0, 0.0); - cut->AddCut(VarManager::kDCAxy1, -0.1, 0.1); - cut->AddCut(VarManager::kDCAxy2, -0.1, 0.1); - cut->AddCut(VarManager::kPt1, 0.5, 100.0, false, VarManager::kPt, 0.0, 1.0); - cut->AddCut(VarManager::kPt1, 0.3, 100.0, false, VarManager::kPt, 1.0, 1.5); - cut->AddCut(VarManager::kPt1, 0.4, 100.0, false, VarManager::kPt, 1.5, 2.0); - cut->AddCut(VarManager::kPt1, 0.7, 100.0, false, VarManager::kPt, 2.0, 50.0); - cut->AddCut(VarManager::kPt2, 0.5, 100.0, false, VarManager::kPt, 0.0, 1.0); - cut->AddCut(VarManager::kPt2, 0.3, 100.0, false, VarManager::kPt, 1.0, 1.5); - cut->AddCut(VarManager::kPt2, 0.4, 100.0, false, VarManager::kPt, 1.5, 2.0); - cut->AddCut(VarManager::kPt2, 0.7, 100.0, false, VarManager::kPt, 2.0, 50.0); + cut->AddCut(VarManager::kD0DcaDaughters, 0.0, 0.01, false, VarManager::kPt, 0.0, 50.0); + cut->AddCut(VarManager::kVertexingLxy, 0.01, 10.0); + cut->AddCut(VarManager::kRap, -2.0, 2.0); + return cut; + } + + if (nameStr == "pairLambdaCTopology") { + cut->AddCut(VarManager::kPt, 0.0, 50.0); + cut->AddCut(VarManager::kRap, -2.0, 2.0); + cut->AddCut(VarManager::kVertexingProcCode, 0.5, 1.e6); + cut->AddCut(VarManager::kCosPointingAngle, 0.80, 1.0, false, VarManager::kPt, 0.0, 1.0); + cut->AddCut(VarManager::kCosPointingAngle, 0.90, 1.0, false, VarManager::kPt, 1.0, 4.0); + cut->AddCut(VarManager::kCosPointingAngle, 0.95, 1.0, false, VarManager::kPt, 4.0, 50.0); + cut->AddCut(VarManager::kVertexingLxy, 0.0, 10.0); + cut->AddCut(VarManager::kVertexingLxy, 0.002, 10.0, false, VarManager::kPt, 0.0, 1.0); + cut->AddCut(VarManager::kVertexingLxy, 0.003, 10.0, false, VarManager::kPt, 1.0, 4.0); + cut->AddCut(VarManager::kVertexingLxy, 0.005, 10.0, false, VarManager::kPt, 4.0, 50.0); + return cut; } @@ -7467,12 +7482,13 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName) } if (nameStr == "alice3TOFEta") { - cut->AddCut(VarManager::kPt, 0.7, 1000.0); + cut->AddCut(VarManager::kPt, 0.3, 1000.0); cut->AddCut(VarManager::kEta, -2.0, 2.0); // Total tracker acceptance in v3b geomety return cut; } if (nameStr == "alice3RICHEta") { + cut->AddCut(VarManager::kPt, 0.3, 1000.0); cut->AddCut(VarManager::kEta, -0.8, 0.8); // Total tracker acceptance in v3b geomety return cut; } @@ -7519,6 +7535,14 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName) return cut; } + if (nameStr == "alice3TrackQualityHF") { + cut->AddCut(VarManager::kIsReconstructed, 0.5, 1.5); + cut->AddCut(VarManager::kNSiliconHits, 5.0, 12.0); + cut->AddCut(VarManager::kTrackDCAxy, -0.05, 0.05); + cut->AddCut(VarManager::kTrackDCAz, -0.05, 0.05); + return cut; + } + if (nameStr == "alice3TrackQualityTightDCA") { cut->AddCut(VarManager::kIsReconstructed, 0.5, 1.5); cut->AddCut(VarManager::kNSiliconHits, 6.0, 12.0); @@ -7643,6 +7667,78 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName) return cut; } + if (nameStr == "alice3KaonTOFRICH") { + cut->AddCut(VarManager::kEta, -2.0, 2.0); + cut->AddCut(VarManager::kP, 0.3, 1000.0); + + cut->AddCut(VarManager::kTOFEventTimeErr, 0.0, 30.0, false, VarManager::kP, 0.0, 2.0); + cut->AddCut(VarManager::kOuterTOFnSigmaKa, -3.0, 3.0, false, VarManager::kP, 0.0f, 2.0); + + cut->AddCut(VarManager::kOuterTOFnSigmaPi, -2.0, 2.0, true, VarManager::kP, 0.0, 1.9); + cut->AddCut(VarManager::kOuterTOFnSigmaPr, -2.0, 2.0, true, VarManager::kP, 0.0, 3.2); + cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 2.0, true, VarManager::kP, 0.0, 2.0); + + cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 2.0, 1000.0); + cut->AddCut(VarManager::kRICHnSigmaKa, 0.5, 1.5, false, VarManager::kP, 2.0, 1000.0); + cut->AddCut(VarManager::kRICHnSigmaKa, -3.0, 3.0, false, VarManager::kP, 2.0, 1000.0); + + return cut; + } + + if (nameStr == "alice3PionTOFRICH") { + cut->AddCut(VarManager::kEta, -2.0, 2.0); + cut->AddCut(VarManager::kP, 0.3, 1000.0); + + cut->AddCut(VarManager::kTOFEventTimeErr, 0.0, 30.0, false, VarManager::kP, 0.0, 2.0); + cut->AddCut(VarManager::kOuterTOFnSigmaPi, -3.0, 3.0, false, VarManager::kP, 0.0, 2.0); + + cut->AddCut(VarManager::kOuterTOFnSigmaKa, -2.0, 2.0, true, VarManager::kP, 0.0, 1.9); + cut->AddCut(VarManager::kOuterTOFnSigmaPr, -2.0, 2.0, true, VarManager::kP, 0.0, 3.8); + cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 2.0, true, VarManager::kP, 0.0, 0.55); + + cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 2.0, 1000.0); + cut->AddCut(VarManager::kRICHnSigmaPi, 0.5, 1.5, false, VarManager::kP, 2.0, 1000.0); + cut->AddCut(VarManager::kRICHnSigmaPi, -3.0, 3.0, false, VarManager::kP, 2.0, 1000.0); + + return cut; + } + + if (nameStr == "alice3ProtonTOFRICH") { + cut->AddCut(VarManager::kEta, -2.0, 2.0); + cut->AddCut(VarManager::kP, 0.3f, 1000.0); + + cut->AddCut(VarManager::kTOFEventTimeErr, 0.0f, 30.0f, false, VarManager::kP, 0.0f, 4.5); + cut->AddCut(VarManager::kOuterTOFnSigmaPr, -3.0f, 3.0f, false, VarManager::kP, 0.0f, 4.5); + + cut->AddCut(VarManager::kOuterTOFnSigmaKa, -2.0, 2.0, true, VarManager::kP, 0.0, 3.2); + cut->AddCut(VarManager::kOuterTOFnSigmaPi, -2.0, 2.0, true, VarManager::kP, 0.0, 3.8); + cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 2.0, true, VarManager::kP, 0.0, 3.8); + + cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 4.5, 1000.0); + cut->AddCut(VarManager::kHasRICHSigPr, 0.5f, 1.5f, false, VarManager::kP, 4.5, 1000.0); + cut->AddCut(VarManager::kRICHnSigmaPr, -3.0f, 3.0f, false, VarManager::kP, 4.5, 1000.0); + + return cut; + } + + if (nameStr == "alice3ElectronTOFRICH") { + cut->AddCut(VarManager::kEta, -2.0, 2.0); + cut->AddCut(VarManager::kP, 0.3f, 1000.0); + + cut->AddCut(VarManager::kTOFEventTimeErr, 0.0f, 30.0f, false, VarManager::kP, 0.0f, 0.6); + cut->AddCut(VarManager::kOuterTOFnSigmaEl, -3.0f, 3.0f, false, VarManager::kP, 0.0f, 0.6); + + cut->AddCut(VarManager::kOuterTOFnSigmaPi, -2.0f, 2.0f, true, VarManager::kP, 0.0f, 0.6); + cut->AddCut(VarManager::kOuterTOFnSigmaKa, -2.0f, 2.0f, true, VarManager::kP, 0.0f, 0.6); + cut->AddCut(VarManager::kOuterTOFnSigmaPr, -2.0f, 2.0f, true, VarManager::kP, 0.0f, 0.6); + + cut->AddCut(VarManager::kEta, -0.8, 0.8, false, VarManager::kP, 0.6, 1000.0); + cut->AddCut(VarManager::kHasRICHSigEl, 0.5f, 1.5f, false, VarManager::kP, 0.6, 1000.0); + cut->AddCut(VarManager::kRICHnSigmaEl, -3.0f, 3.0f, false, VarManager::kP, 0.6, 1000.0); + + return cut; + } + if (nameStr == "alice3CharmoniumPID") { cut->AddCut(VarManager::kOuterTOFnSigmaEl, -2.0, 3.0, false, VarManager::kP, 0.0, 1.2); cut->AddCut(VarManager::kRICHnSigmaEl, -2.0, 3.0, false, VarManager::kHasRICHSigEl, 0.5, 1.5, false, VarManager::kP, 0.7, 1000.0);