Skip to content

Commit ba76a2a

Browse files
ThePhDaneThePhDane
authored
[PWGDQ] Added a few PID cuts to the cuts library (AliceO2Group#9362)
1 parent 96b0592 commit ba76a2a

File tree

1 file changed

+180
-0
lines changed

1 file changed

+180
-0
lines changed

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 180 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1016,6 +1016,153 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
10161016
}
10171017
}
10181018

1019+
// Magnus composite cuts -----------------------------------------------------------------------------------------------------------------
1020+
1021+
AnalysisCompositeCut* magnus_PID111 = new AnalysisCompositeCut("magnus_PID111", "");
1022+
magnus_PID111->AddCut(GetAnalysisCut("pidJpsi_magnus_ele1"));
1023+
magnus_PID111->AddCut(GetAnalysisCut("pidJpsi_magnus_pion1"));
1024+
magnus_PID111->AddCut(GetAnalysisCut("pidJpsi_magnus_prot1"));
1025+
if (!nameStr.compare("MagnussOptimization111")) {
1026+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1027+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1028+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1029+
cut->AddCut(magnus_PID111);
1030+
return cut;
1031+
}
1032+
1033+
AnalysisCompositeCut* magnus_PID211 = new AnalysisCompositeCut("magnus_PID211", "");
1034+
magnus_PID211->AddCut(GetAnalysisCut("pidJpsi_magnus_ele2"));
1035+
magnus_PID211->AddCut(GetAnalysisCut("pidJpsi_magnus_pion1"));
1036+
magnus_PID211->AddCut(GetAnalysisCut("pidJpsi_magnus_prot1"));
1037+
if (!nameStr.compare("MagnussOptimization211")) {
1038+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1039+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1040+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1041+
cut->AddCut(magnus_PID211);
1042+
return cut;
1043+
}
1044+
1045+
AnalysisCompositeCut* magnus_PID311 = new AnalysisCompositeCut("magnus_PID311", "");
1046+
magnus_PID311->AddCut(GetAnalysisCut("pidJpsi_magnus_ele3"));
1047+
magnus_PID311->AddCut(GetAnalysisCut("pidJpsi_magnus_pion1"));
1048+
magnus_PID311->AddCut(GetAnalysisCut("pidJpsi_magnus_prot1"));
1049+
if (!nameStr.compare("MagnussOptimization311")) {
1050+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1051+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1052+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1053+
cut->AddCut(magnus_PID311);
1054+
return cut;
1055+
}
1056+
1057+
AnalysisCompositeCut* magnus_PID121 = new AnalysisCompositeCut("magnus_PID121", "");
1058+
magnus_PID121->AddCut(GetAnalysisCut("pidJpsi_magnus_ele1"));
1059+
magnus_PID121->AddCut(GetAnalysisCut("pidJpsi_magnus_pion2"));
1060+
magnus_PID121->AddCut(GetAnalysisCut("pidJpsi_magnus_prot1"));
1061+
if (!nameStr.compare("MagnussOptimization121")) {
1062+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1063+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1064+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1065+
cut->AddCut(magnus_PID121);
1066+
return cut;
1067+
}
1068+
1069+
AnalysisCompositeCut* magnus_PID112 = new AnalysisCompositeCut("magnus_PID112", "");
1070+
magnus_PID112->AddCut(GetAnalysisCut("pidJpsi_magnus_ele1"));
1071+
magnus_PID112->AddCut(GetAnalysisCut("pidJpsi_magnus_pion1"));
1072+
magnus_PID112->AddCut(GetAnalysisCut("pidJpsi_magnus_prot2"));
1073+
if (!nameStr.compare("MagnussOptimization112")) {
1074+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1075+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1076+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1077+
cut->AddCut(magnus_PID112);
1078+
return cut;
1079+
}
1080+
1081+
AnalysisCompositeCut* magnus_PID122 = new AnalysisCompositeCut("magnus_PID122", "");
1082+
magnus_PID122->AddCut(GetAnalysisCut("pidJpsi_magnus_ele1"));
1083+
magnus_PID122->AddCut(GetAnalysisCut("pidJpsi_magnus_pion2"));
1084+
magnus_PID122->AddCut(GetAnalysisCut("pidJpsi_magnus_prot2"));
1085+
if (!nameStr.compare("MagnussOptimization122")) {
1086+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1087+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1088+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1089+
cut->AddCut(magnus_PID122);
1090+
return cut;
1091+
}
1092+
1093+
AnalysisCompositeCut* magnus_PID222 = new AnalysisCompositeCut("magnus_PID222", "");
1094+
magnus_PID222->AddCut(GetAnalysisCut("pidJpsi_magnus_ele2"));
1095+
magnus_PID222->AddCut(GetAnalysisCut("pidJpsi_magnus_pion2"));
1096+
magnus_PID222->AddCut(GetAnalysisCut("pidJpsi_magnus_prot2"));
1097+
if (!nameStr.compare("MagnussOptimization222")) {
1098+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1099+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1100+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1101+
cut->AddCut(magnus_PID222);
1102+
return cut;
1103+
}
1104+
1105+
AnalysisCompositeCut* magnus_PID212 = new AnalysisCompositeCut("magnus_PID212", "");
1106+
magnus_PID212->AddCut(GetAnalysisCut("pidJpsi_magnus_ele2"));
1107+
magnus_PID212->AddCut(GetAnalysisCut("pidJpsi_magnus_pion1"));
1108+
magnus_PID212->AddCut(GetAnalysisCut("pidJpsi_magnus_prot2"));
1109+
if (!nameStr.compare("MagnussOptimization212")) {
1110+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1111+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1112+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1113+
cut->AddCut(magnus_PID212);
1114+
return cut;
1115+
}
1116+
1117+
AnalysisCompositeCut* magnus_PID221 = new AnalysisCompositeCut("magnus_PID221", "");
1118+
magnus_PID221->AddCut(GetAnalysisCut("pidJpsi_magnus_ele2"));
1119+
magnus_PID221->AddCut(GetAnalysisCut("pidJpsi_magnus_pion2"));
1120+
magnus_PID221->AddCut(GetAnalysisCut("pidJpsi_magnus_prot1"));
1121+
if (!nameStr.compare("MagnussOptimization221")) {
1122+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1123+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1124+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1125+
cut->AddCut(magnus_PID221);
1126+
return cut;
1127+
}
1128+
1129+
AnalysisCompositeCut* magnus_PID321 = new AnalysisCompositeCut("magnus_PID321", "");
1130+
magnus_PID321->AddCut(GetAnalysisCut("pidJpsi_magnus_ele3"));
1131+
magnus_PID321->AddCut(GetAnalysisCut("pidJpsi_magnus_pion2"));
1132+
magnus_PID321->AddCut(GetAnalysisCut("pidJpsi_magnus_prot1"));
1133+
if (!nameStr.compare("MagnussOptimization321")) {
1134+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1135+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1136+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1137+
cut->AddCut(magnus_PID321);
1138+
return cut;
1139+
}
1140+
1141+
AnalysisCompositeCut* magnus_PID312 = new AnalysisCompositeCut("magnus_PID312", "");
1142+
magnus_PID312->AddCut(GetAnalysisCut("pidJpsi_magnus_ele3"));
1143+
magnus_PID312->AddCut(GetAnalysisCut("pidJpsi_magnus_pion1"));
1144+
magnus_PID312->AddCut(GetAnalysisCut("pidJpsi_magnus_prot2"));
1145+
if (!nameStr.compare("MagnussOptimization312")) {
1146+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1147+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1148+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1149+
cut->AddCut(magnus_PID312);
1150+
return cut;
1151+
}
1152+
1153+
AnalysisCompositeCut* magnus_PID322 = new AnalysisCompositeCut("magnus_PID322", "");
1154+
magnus_PID322->AddCut(GetAnalysisCut("pidJpsi_magnus_ele1"));
1155+
magnus_PID322->AddCut(GetAnalysisCut("pidJpsi_magnus_pion2"));
1156+
magnus_PID322->AddCut(GetAnalysisCut("pidJpsi_magnus_prot2"));
1157+
if (!nameStr.compare("MagnussOptimization322")) {
1158+
cut->AddCut(GetAnalysisCut("kineJpsiEle_ionut"));
1159+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
1160+
cut->AddCut(GetAnalysisCut("trackQuality_ionut"));
1161+
cut->AddCut(magnus_PID322);
1162+
return cut;
1163+
}
1164+
//-------------------------------------------------------------------------------------------------------
1165+
10191166
//---------------------------------------------------------------
10201167
// Cuts for the selection of legs from dalitz decay
10211168
//
@@ -4431,6 +4578,39 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
44314578
return cut;
44324579
}
44334580

4581+
// Magnus cuts ----------------------------------------------------------
4582+
4583+
if (!nameStr.compare("pidJpsi_magnus_ele1")) {
4584+
cut->AddCut(VarManager::kTPCnSigmaEl, -3.0, 4.0);
4585+
return cut;
4586+
}
4587+
if (!nameStr.compare("pidJpsi_magnus_ele2")) {
4588+
cut->AddCut(VarManager::kTPCnSigmaEl, -2.0, 4.0);
4589+
return cut;
4590+
}
4591+
if (!nameStr.compare("pidJpsi_magnus_ele3")) {
4592+
cut->AddCut(VarManager::kTPCnSigmaEl, -1.0, 4.0);
4593+
return cut;
4594+
}
4595+
if (!nameStr.compare("pidJpsi_magnus_prot1")) {
4596+
cut->AddCut(VarManager::kTPCnSigmaPr, 3.0, 1000.0);
4597+
return cut;
4598+
}
4599+
if (!nameStr.compare("pidJpsi_magnus_prot2")) {
4600+
cut->AddCut(VarManager::kTPCnSigmaPr, 3.5, 1000.0);
4601+
return cut;
4602+
}
4603+
if (!nameStr.compare("pidJpsi_magnus_pion1")) {
4604+
cut->AddCut(VarManager::kTPCnSigmaPi, 3.0, 1000.0);
4605+
return cut;
4606+
}
4607+
if (!nameStr.compare("pidJpsi_magnus_pion2")) {
4608+
cut->AddCut(VarManager::kTPCnSigmaPi, 3.5, 1000.0);
4609+
return cut;
4610+
}
4611+
4612+
// ----------------------------------------------------------------------------------
4613+
44344614
if (!nameStr.compare("standardPrimaryTrackDCAz")) {
44354615
cut->AddCut(VarManager::kTrackDCAxy, -3.0, 3.0);
44364616
cut->AddCut(VarManager::kTrackDCAz, -1.0, 1.0);

0 commit comments

Comments
 (0)