Skip to content

Commit 11646b4

Browse files
glegrasglegras
committed
define TF1 only when needed
1 parent 8926f02 commit 11646b4

1 file changed

Lines changed: 67 additions & 63 deletions

File tree

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 67 additions & 63 deletions
Original file line numberDiff line numberDiff line change
@@ -5095,40 +5095,42 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
50955095
return cut;
50965096
}
50975097

5098-
TF1* fDCAxyresLow = new TF1("fDCAxyresLow", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5099-
TF1* fDCAzresLow = new TF1("fDCAzresLow", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5100-
TF1* fDCAxyresUp = new TF1("fDCAxyresUp", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5101-
TF1* fDCAzresUp = new TF1("fDCAzresUp", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5102-
5103-
if (!nameStr.compare("electronPrimary_dca3sigma")) {
5104-
// DCAxy and DCAz 3 sigma cut. DCA resolution vs pt extracted from fits of Users/m/mfaggin/test/inputsTrackTuner/pp2024/pass1_minBias/vsPhi (used for the track tuner)
5105-
// we add in addition a term for the misalignment of the mean of the distribution, which seems to be at most 20 mum for DCAxy and 10 mum for DCAz
5106-
fDCAxyresLow->SetParameters(-3 * 8.7e-4 - 20e-4, -3 * 25.4e-4, 0.79); // res is 8.7 + 25.4/pt^0.79 mum
5107-
fDCAzresLow->SetParameters(-3 * 9.4e-4 - 10e-4, -3 * 26.5e-4, 0.79); // res is 9.4 + 26.5/pt^0.79 mum
5108-
fDCAxyresUp->SetParameters(3 * 8.7e-4 + 20e-4, 3 * 25.4e-4, 0.79); // res is 8.7 + 25.4/pt^0.79 mum
5109-
fDCAzresUp->SetParameters(3 * 9.4e-4 + 10e-4, 3 * 26.5e-4, 0.79); // res is 9.4 + 26.5/pt^0.79 mum
5110-
cut->AddCut(VarManager::kTrackDCAxy, fDCAxyresLow, fDCAxyresUp, false, VarManager::kPt, 0.2, 1000.);
5111-
cut->AddCut(VarManager::kTrackDCAz, fDCAzresLow, fDCAzresUp, false, VarManager::kPt, 0.2, 1000.);
5112-
cut->AddCut(VarManager::kIsITSibFirst, 0.5, 1.5);
5113-
cut->AddCut(VarManager::kITSchi2, 0.0, 5.0);
5114-
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
5115-
cut->AddCut(VarManager::kTPCncls, 100, 161.);
5116-
return cut;
5117-
}
5098+
if (!nameStr.compare("electronPrimary_dca3sigma") || !nameStr.compare("electronPrimary_dca7sigma")) {
5099+
TF1* fDCAxyresLow = new TF1("fDCAxyresLow", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5100+
TF1* fDCAzresLow = new TF1("fDCAzresLow", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5101+
TF1* fDCAxyresUp = new TF1("fDCAxyresUp", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5102+
TF1* fDCAzresUp = new TF1("fDCAzresUp", "[0] + [1] * pow(x, -[2])", 0.1, 1000.);
5103+
5104+
if (!nameStr.compare("electronPrimary_dca3sigma")) {
5105+
// DCAxy and DCAz 3 sigma cut. DCA resolution vs pt extracted from fits of Users/m/mfaggin/test/inputsTrackTuner/pp2024/pass1_minBias/vsPhi (used for the track tuner)
5106+
// we add in addition a term for the misalignment of the mean of the distribution, which seems to be at most 20 mum for DCAxy and 10 mum for DCAz
5107+
fDCAxyresLow->SetParameters(-3 * 8.7e-4 - 20e-4, -3 * 25.4e-4, 0.79); // res is 8.7 + 25.4/pt^0.79 mum
5108+
fDCAzresLow->SetParameters(-3 * 9.4e-4 - 10e-4, -3 * 26.5e-4, 0.79); // res is 9.4 + 26.5/pt^0.79 mum
5109+
fDCAxyresUp->SetParameters(3 * 8.7e-4 + 20e-4, 3 * 25.4e-4, 0.79); // res is 8.7 + 25.4/pt^0.79 mum
5110+
fDCAzresUp->SetParameters(3 * 9.4e-4 + 10e-4, 3 * 26.5e-4, 0.79); // res is 9.4 + 26.5/pt^0.79 mum
5111+
cut->AddCut(VarManager::kTrackDCAxy, fDCAxyresLow, fDCAxyresUp, false, VarManager::kPt, 0.2, 1000.);
5112+
cut->AddCut(VarManager::kTrackDCAz, fDCAzresLow, fDCAzresUp, false, VarManager::kPt, 0.2, 1000.);
5113+
cut->AddCut(VarManager::kIsITSibFirst, 0.5, 1.5);
5114+
cut->AddCut(VarManager::kITSchi2, 0.0, 5.0);
5115+
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
5116+
cut->AddCut(VarManager::kTPCncls, 100, 161.);
5117+
return cut;
5118+
}
51185119

5119-
if (!nameStr.compare("electronPrimary_dca7sigma")) {
5120-
// DCAxy and DCAz 7 sigma cut
5121-
fDCAxyresLow->SetParameters(-7 * 8.7e-4 - 20e-4, -7 * 25.4e-4, 0.79);
5122-
fDCAzresLow->SetParameters(-7 * 9.4e-4 - 10e-4, -7 * 26.5e-4, 0.79);
5123-
fDCAxyresUp->SetParameters(7 * 8.7e-4 + 20e-4, 7 * 25.4e-4, 0.79);
5124-
fDCAzresUp->SetParameters(7 * 9.4e-4 + 10e-4, 7 * 26.5e-4, 0.79);
5125-
cut->AddCut(VarManager::kTrackDCAxy, fDCAxyresLow, fDCAxyresUp, false, VarManager::kPt, 0.2, 1000.);
5126-
cut->AddCut(VarManager::kTrackDCAz, fDCAzresLow, fDCAzresUp, false, VarManager::kPt, 0.2, 1000.);
5127-
cut->AddCut(VarManager::kIsITSibFirst, 0.5, 1.5);
5128-
cut->AddCut(VarManager::kITSchi2, 0.0, 5.0);
5129-
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
5130-
cut->AddCut(VarManager::kTPCncls, 100, 161.);
5131-
return cut;
5120+
if (!nameStr.compare("electronPrimary_dca7sigma")) {
5121+
// DCAxy and DCAz 7 sigma cut
5122+
fDCAxyresLow->SetParameters(-7 * 8.7e-4 - 20e-4, -7 * 25.4e-4, 0.79);
5123+
fDCAzresLow->SetParameters(-7 * 9.4e-4 - 10e-4, -7 * 26.5e-4, 0.79);
5124+
fDCAxyresUp->SetParameters(7 * 8.7e-4 + 20e-4, 7 * 25.4e-4, 0.79);
5125+
fDCAzresUp->SetParameters(7 * 9.4e-4 + 10e-4, 7 * 26.5e-4, 0.79);
5126+
cut->AddCut(VarManager::kTrackDCAxy, fDCAxyresLow, fDCAxyresUp, false, VarManager::kPt, 0.2, 1000.);
5127+
cut->AddCut(VarManager::kTrackDCAz, fDCAzresLow, fDCAzresUp, false, VarManager::kPt, 0.2, 1000.);
5128+
cut->AddCut(VarManager::kIsITSibFirst, 0.5, 1.5);
5129+
cut->AddCut(VarManager::kITSchi2, 0.0, 5.0);
5130+
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
5131+
cut->AddCut(VarManager::kTPCncls, 100, 161.);
5132+
return cut;
5133+
}
51325134
}
51335135

51345136
if (!nameStr.compare("dcaCut1_ionut")) {
@@ -5487,40 +5489,42 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
54875489
return cut;
54885490
}
54895491

5490-
TF1* cutLow1 = new TF1("cutLow1", "pol1", 0., 10.);
5491-
if (!nameStr.compare("electronPID1")) {
5492-
cutLow1->SetParameters(130., -40.0);
5493-
cut->AddCut(VarManager::kTPCsignal, 70., 100.);
5494-
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0, false, VarManager::kPin, 0.5, 3.0);
5495-
return cut;
5496-
}
5492+
if (!nameStr.compare("electronPID1") || !nameStr.compare("electronPID1shiftUp") || !nameStr.compare("electronPID1shiftDown") || !nameStr.compare("electronPID2") || !nameStr.compare("electronPID3")) {
5493+
TF1* cutLow1 = new TF1("cutLow1", "pol1", 0., 10.);
5494+
if (!nameStr.compare("electronPID1")) {
5495+
cutLow1->SetParameters(130., -40.0);
5496+
cut->AddCut(VarManager::kTPCsignal, 70., 100.);
5497+
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0, false, VarManager::kPin, 0.5, 3.0);
5498+
return cut;
5499+
}
54975500

5498-
if (!nameStr.compare("electronPID1shiftUp")) {
5499-
cut->AddCut(VarManager::kTPCsignal, 70. - 0.85, 100. - 0.85);
5500-
cutLow1->SetParameters(130. - 0.85, -40.0);
5501-
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0 - 0.85, false, VarManager::kPin, 0.5, 3.0);
5502-
return cut;
5503-
}
5501+
if (!nameStr.compare("electronPID1shiftUp")) {
5502+
cut->AddCut(VarManager::kTPCsignal, 70. - 0.85, 100. - 0.85);
5503+
cutLow1->SetParameters(130. - 0.85, -40.0);
5504+
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0 - 0.85, false, VarManager::kPin, 0.5, 3.0);
5505+
return cut;
5506+
}
55045507

5505-
if (!nameStr.compare("electronPID1shiftDown")) {
5506-
cut->AddCut(VarManager::kTPCsignal, 70.0 + 0.85, 100.0 + 0.85);
5507-
cutLow1->SetParameters(130. + 0.85, -40.0);
5508-
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0 + 0.85, false, VarManager::kPin, 0.5, 3.0);
5509-
return cut;
5510-
}
5508+
if (!nameStr.compare("electronPID1shiftDown")) {
5509+
cut->AddCut(VarManager::kTPCsignal, 70.0 + 0.85, 100.0 + 0.85);
5510+
cutLow1->SetParameters(130. + 0.85, -40.0);
5511+
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0 + 0.85, false, VarManager::kPin, 0.5, 3.0);
5512+
return cut;
5513+
}
55115514

5512-
if (!nameStr.compare("electronPID2")) {
5513-
cutLow1->SetParameters(130., -40.0);
5514-
cut->AddCut(VarManager::kTPCsignal, 73., 100.);
5515-
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0, false, VarManager::kPin, 0.5, 3.0);
5516-
return cut;
5517-
}
5515+
if (!nameStr.compare("electronPID2")) {
5516+
cutLow1->SetParameters(130., -40.0);
5517+
cut->AddCut(VarManager::kTPCsignal, 73., 100.);
5518+
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0, false, VarManager::kPin, 0.5, 3.0);
5519+
return cut;
5520+
}
55185521

5519-
if (!nameStr.compare("electronPID3")) {
5520-
cutLow1->SetParameters(130., -40.0);
5521-
cut->AddCut(VarManager::kTPCsignal, 60., 110.);
5522-
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0, false, VarManager::kPin, 0.5, 3.0);
5523-
return cut;
5522+
if (!nameStr.compare("electronPID3")) {
5523+
cutLow1->SetParameters(130., -40.0);
5524+
cut->AddCut(VarManager::kTPCsignal, 60., 110.);
5525+
cut->AddCut(VarManager::kTPCsignal, cutLow1, 100.0, false, VarManager::kPin, 0.5, 3.0);
5526+
return cut;
5527+
}
55245528
}
55255529

55265530
if (!nameStr.compare("electronPIDnsigma")) {

0 commit comments

Comments
 (0)