Skip to content

PWGDQ : Adjustment of Cuts in the CutsLibrary.cxx #4128

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
iarsene merged 1 commit into from
Dec 11, 2023
Merged

PWGDQ : Adjustment of Cuts in the CutsLibrary.cxx #4128

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
179 changes: 74 additions & 105 deletions PWGDQ/Core/CutsLibrary.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -28,27 +28,53 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
// These are the Cuts used in the CEFP Task //
// to select tracks in the event selection //
// ///////////////////////////////////////////////
if (!nameStr.compare("Electron")) {
if (!nameStr.compare("Electron2022")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_debug5"));
return cut;
}
if (!nameStr.compare("ElectronTight")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine3"));
if (!nameStr.compare("Electron2023")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_debug1"));
cut->AddCut(GetAnalysisCut("pidCut_lowP_Corr"));

AnalysisCompositeCut* pidCut_highP = new AnalysisCompositeCut("pidCut_highP", "pidCut_highP", kFALSE);
pidCut_highP->AddCut(GetAnalysisCut("EleInclusion_highP_Corr"));
pidCut_highP->AddCut(GetAnalysisCut("PionExclusion_highP_Corr"));
cut->AddCut(pidCut_highP);
return cut;
}
if (!nameStr.compare("Electron2023_Tight")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("pidCut_lowP_Corr"));

AnalysisCompositeCut* pidCut_highP = new AnalysisCompositeCut("pidCut_highP", "pidCut_highP", kFALSE);
pidCut_highP->AddCut(GetAnalysisCut("EleInclusion_highP2_Corr"));
pidCut_highP->AddCut(GetAnalysisCut("PionExclusion_highP_Corr"));
cut->AddCut(pidCut_highP);
return cut;
}
if (!nameStr.compare("MuonLow")) {
if (!nameStr.compare("MuonLow2022")) {
cut->AddCut(GetAnalysisCut("muonLowPt2"));
cut->AddCut(GetAnalysisCut("muonQualityCuts"));
return cut;
}
if (!nameStr.compare("MuonLow2023")) {
cut->AddCut(GetAnalysisCut("muonLowPt2"));
cut->AddCut(GetAnalysisCut("muonQualityCuts"));
cut->AddCut(GetAnalysisCut("MCHMID"));
return cut;
}
if (!nameStr.compare("MuonHigh")) {
if (!nameStr.compare("MuonHigh2022")) {
cut->AddCut(GetAnalysisCut("muonHighPt2"));
cut->AddCut(GetAnalysisCut("muonQualityCuts"));
return cut;
}
if (!nameStr.compare("MuonHigh2023")) {
cut->AddCut(GetAnalysisCut("muonHighPt6"));
cut->AddCut(GetAnalysisCut("muonQualityCuts"));
cut->AddCut(GetAnalysisCut("MCHMID"));
return cut;
}
Expand Down Expand Up @@ -214,95 +240,53 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
cut->AddCut(GetAnalysisCut("electronPIDnsigmaVeryLoose")); // with 3 sigma El TOF
return cut;
}
//
// New cuts to test for jpsi2ee electron selection (potentially not useful)
//
// With pT > 1.0 GeV/c as primary cut
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_1")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_1")); // Loose pion rejection at high p with no correction maps

return cut;
}
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_1_Corr")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_1_Corr")); // Loose pion rejection at high p with correction maps

return cut;
}
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_2")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_2")); // No pion rejection at high p and asymmetrical e selection with no correction maps

return cut;
}
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_2_Corr")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_2_Corr")); // No pion rejection at high p and asymmetrical e selection with correction maps

return cut;
}

// With p > 1.5 GeV/c as primary cut
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_3")) {
if (!nameStr.compare("jpsiO2MCdebugCuts_Pdependent_Corr")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_1")); // Loose pion rejection at high p with no correction maps
cut->AddCut(GetAnalysisCut("pidCut_lowP_Corr"));

AnalysisCompositeCut* pidCut_highP = new AnalysisCompositeCut("pidCut_highP", "pidCut_highP", kFALSE);
pidCut_highP->AddCut(GetAnalysisCut("EleInclusion_highP_Corr"));
pidCut_highP->AddCut(GetAnalysisCut("PionExclusion_highP_Corr"));
cut->AddCut(pidCut_highP);
return cut;
}
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_3_Corr")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_1_Corr")); // Loose pion rejection at high p with correction maps

return cut;
}
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_4")) {
if (!nameStr.compare("jpsiO2MCdebugCuts_Pdependent")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_2")); // No pion rejection at high p and asymmetrical e selection with no correction maps
cut->AddCut(GetAnalysisCut("pidCut_lowP"));

AnalysisCompositeCut* pidCut_highP = new AnalysisCompositeCut("pidCut_highP", "pidCut_highP", kFALSE);
pidCut_highP->AddCut(GetAnalysisCut("EleInclusion_highP"));
pidCut_highP->AddCut(GetAnalysisCut("PionExclusion_highP"));
cut->AddCut(pidCut_highP);
return cut;
}
if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_4_Corr")) {
if (!nameStr.compare("jpsiO2MCdebugCuts_Pdependent2_Corr")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("jpsi_TPCPID_TEST_2_Corr")); // No pion rejection at high p and asymmetrical e selection with correction maps

return cut;
}

if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_6_Corr")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("pidCut_lowP"));
cut->AddCut(GetAnalysisCut("pidCut_lowP_Corr"));

AnalysisCompositeCut* pidCut_highP = new AnalysisCompositeCut("pidCut_highP", "pidCut_highP", kFALSE);
pidCut_highP->AddCut(GetAnalysisCut("MySpecial_EleInclusion"));
pidCut_highP->AddCut(GetAnalysisCut("MySpecial_PionExclusion"));
pidCut_highP->AddCut(GetAnalysisCut("EleInclusion_highP2_Corr"));
pidCut_highP->AddCut(GetAnalysisCut("PionExclusion_highP_Corr"));
cut->AddCut(pidCut_highP);
return cut;
}

if (!nameStr.compare("jpsiO2MCdebugCuts_TEST_5_Corr")) {
if (!nameStr.compare("jpsiO2MCdebugCuts_Pdependent2")) {
cut->AddCut(GetAnalysisCut("jpsiStandardKine4"));
cut->AddCut(GetAnalysisCut("electronStandardQualityForO2MCdebug"));
cut->AddCut(GetAnalysisCut("pidCut_lowP"));

AnalysisCompositeCut* pidCut_highP = new AnalysisCompositeCut("pidCut_highP", "pidCut_highP", kFALSE);
pidCut_highP->AddCut(GetAnalysisCut("MySpecial_EleInclusion"));
pidCut_highP->AddCut(GetAnalysisCut("MySpecial_PionExclusion"));
pidCut_highP->AddCut(GetAnalysisCut("EleInclusion_highP2"));
pidCut_highP->AddCut(GetAnalysisCut("PionExclusion_highP"));
cut->AddCut(pidCut_highP);
return cut;
}
//
// end new test cuts
//

if (!nameStr.compare("JpsiPWGSkimmedCuts1")) { // please do not remove or modify, this is used for the common Skimmed tree production, (Xiaozhi Bai)
cut->AddCut(GetAnalysisCut("jpsiKineSkimmed"));
Expand Down Expand Up @@ -2921,60 +2905,45 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
cut->AddCut(VarManager::kTPCnSigmaPr, 3.5, 999);
return cut;
}
//
// New test for jpsi2ee PID selection
//
if (!nameStr.compare("jpsi_TPCPID_TEST_1_Corr")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -3.0, 3.0);
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 2.0, 999, false, VarManager::kP, 4.0, 999);
cut->AddCut(VarManager::kTPCnSigmaPr_Corr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPr_Corr, 2.0, 999, false, VarManager::kP, 4.0, 999);

if (!nameStr.compare("pidCut_lowP_Corr")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -3.0, 3.0, false, VarManager::kP, 0.0, 5.0);
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 3.0, 999, false, VarManager::kP, 0.0, 5.0);
cut->AddCut(VarManager::kTPCnSigmaPr_Corr, 2.5, 999, false, VarManager::kP, 0.0, 5.0);
return cut;
}
if (!nameStr.compare("jpsi_TPCPID_TEST_1")) {
cut->AddCut(VarManager::kTPCnSigmaEl, -3.0, 3.0);
cut->AddCut(VarManager::kTPCnSigmaPi, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPi, 2.0, 999, false, VarManager::kP, 4.0, 999);
cut->AddCut(VarManager::kTPCnSigmaPr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPr, 2.0, 999, false, VarManager::kP, 4.0, 999);

if (!nameStr.compare("EleInclusion_highP_Corr")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -1.0, 4.0, false, VarManager::kP, 5.0, 999.0);
return cut;
}
if (!nameStr.compare("jpsi_TPCPID_TEST_2_Corr")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -3.0, 3.0, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -1.0, 4.0, false, VarManager::kP, 4.0, 999);
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 2.0, 999, false, VarManager::kP, 4.0, 999);
cut->AddCut(VarManager::kTPCnSigmaPr_Corr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPr_Corr, 2.0, 999, false, VarManager::kP, 4.0, 999);
if (!nameStr.compare("EleInclusion_highP2_Corr")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -0.5, 4.0, false, VarManager::kP, 5.0, 999.0);
return cut;
}
if (!nameStr.compare("jpsi_TPCPID_TEST_2")) {
cut->AddCut(VarManager::kTPCnSigmaEl, -3.0, 3.0, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaEl, -1.0, 4.0, false, VarManager::kP, 4.0, 999.0);
cut->AddCut(VarManager::kTPCnSigmaPi, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPi, 0, 999, false, VarManager::kP, 4.0, 999);
cut->AddCut(VarManager::kTPCnSigmaPr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPr, 0, 999, false, VarManager::kP, 4.0, 999);
if (!nameStr.compare("PionExclusion_highP_Corr")) {
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 2.0, 999, false, VarManager::kP, 5.0, 999.0);
return cut;
}

if (!nameStr.compare("pidCut_lowP")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -3.0, 3.0, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaPr_Corr, 3.0, 999, false, VarManager::kP, 0.0, 4.0);
cut->AddCut(VarManager::kTPCnSigmaEl, -3.0, 3.0, false, VarManager::kP, 0.0, 5.0);
cut->AddCut(VarManager::kTPCnSigmaPi, 3.0, 999, false, VarManager::kP, 0.0, 5.0);
cut->AddCut(VarManager::kTPCnSigmaPr, 2.5, 999, false, VarManager::kP, 0.0, 5.0);
return cut;
}

if (!nameStr.compare("MySpecial_EleInclusion")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -1.0, 4.0, false, VarManager::kP, 4.0, 999.0);
if (!nameStr.compare("EleInclusion_highP")) {
cut->AddCut(VarManager::kTPCnSigmaEl, -1.0, 4.0, false, VarManager::kP, 5.0, 999.0);
return cut;
}
if (!nameStr.compare("EleInclusion_highP2")) {
cut->AddCut(VarManager::kTPCnSigmaEl, -0.5, 4.0, false, VarManager::kP, 5.0, 999.0);
return cut;
}
if (!nameStr.compare("MySpecial_PionExclusion")) {
cut->AddCut(VarManager::kTPCnSigmaPi_Corr, 2.0, 999, false, VarManager::kP, 4.0, 999.0);
if (!nameStr.compare("PionExclusion_highP")) {
cut->AddCut(VarManager::kTPCnSigmaPi, 2.0, 999, false, VarManager::kP, 5.0, 999.0);
return cut;
}
// end of new tests cuts

if (!nameStr.compare("lmee_TPCPID_debug1")) {
cut->AddCut(VarManager::kTPCnSigmaEl_Corr, -5.0, 5.0);
Expand Down