Skip to content

Commit b58008d

Browse files
ypwanggalibuild
ypwangg
and
alibuild
authored
[PWGDQ] Code development for coherent Jpsi A2 analysis (#16698)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 0b36181 commit b58008d

6 files changed

Lines changed: 457 additions & 5 deletions

File tree

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 75 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1593,6 +1593,34 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
15931593
cut->AddCut(GetAnalysisCut("electronPID1shiftDown"));
15941594
return cut;
15951595
}
1596+
// -------------------------------------------------------------------------------------------------
1597+
//
1598+
// Q vector contributor cut
1599+
//
1600+
if (!nameStr.compare("selTPCCentral")) {
1601+
AnalysisCut* kineCut = new AnalysisCut("kineCut", "kine cut");
1602+
kineCut->AddCut(VarManager::kEta, -0.8, 0.8);
1603+
kineCut->AddCut(VarManager::kPt, 0.15, 5);
1604+
1605+
AnalysisCut* qualityCuts = new AnalysisCut("qualityCuts", "quality cuts");
1606+
qualityCuts->AddCut(VarManager::kTPCchi2, 0., 4.);
1607+
qualityCuts->AddCut(VarManager::kTPCnCRoverFindCls, 0.8, 1.);
1608+
qualityCuts->AddCut(VarManager::kIsITSibAny, 0.5, 1.5);
1609+
qualityCuts->AddCut(VarManager::kITSchi2, 0., 36.);
1610+
1611+
AnalysisCut* dcaCuts = new AnalysisCut("dcaCuts", "DCA cuts");
1612+
std::shared_ptr<TF1> f1dcaxyHigh = std::make_shared<TF1>("f1dcaxy", "[0]+[1]/pow(x,[2])", 0., 10.);
1613+
f1dcaxyHigh->SetParameters(0.0105, 0.035, 1.1);
1614+
std::shared_ptr<TF1> f1dcaxyLow = std::make_shared<TF1>("f1dcaxy_low", "[0]+[1]/pow(x,[2])", 0., 10.);
1615+
f1dcaxyLow->SetParameters(-0.0105, -0.035, 1.1);
1616+
dcaCuts->AddCut(VarManager::kTrackDCAxy, f1dcaxyLow, f1dcaxyHigh);
1617+
dcaCuts->AddCut(VarManager::kTrackDCAz, -2., 2.);
1618+
1619+
cut->AddCut(kineCut);
1620+
cut->AddCut(qualityCuts);
1621+
cut->AddCut(dcaCuts);
1622+
}
1623+
15961624
// -------------------------------------------------------------------------------------------------
15971625
//
15981626
// LMee cuts
@@ -3001,6 +3029,17 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
30013029
return cut;
30023030
}
30033031

3032+
if (!nameStr.compare("jpsi_debug_TPCTOF3_rejBadTOF")) {
3033+
cut->AddCut(GetAnalysisCut("jpsiStandardKine5"));
3034+
cut->AddCut(GetAnalysisCut("electronStandardQualityTPCOnly3"));
3035+
cut->AddCut(GetAnalysisCut("SPDfirst"));
3036+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
3037+
cut->AddCut(GetAnalysisCut("pidJpsi_TPCpion0"));
3038+
cut->AddCut(GetAnalysisCut("pidJpsi_beta"));
3039+
cut->AddCut(GetAnalysisCut("pidJpsi_noTOF_prot"));
3040+
return cut;
3041+
}
3042+
30043043
// -------------------------------------------------------------------------------------------------
30053044
// lmee pair cuts
30063045

@@ -4141,6 +4180,15 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
41414180
return cut;
41424181
}
41434182

4183+
if (!nameStr.compare("eventStandardSel8NoPileup")) {
4184+
cut->AddCut(VarManager::kVtxZ, -10.0, 10.0);
4185+
cut->AddCut(VarManager::kIsSel8, 0.5, 1.5);
4186+
cut->AddCut(VarManager::kIsNoSameBunch, 0.5, 1.5);
4187+
cut->AddCut(VarManager::kIsGoodZvtxFT0vsPV, 0.5, 1.5);
4188+
cut->AddCut(VarManager::kNoCollInTimeRangeStandard, 0.5, 1.5);
4189+
return cut;
4190+
}
4191+
41444192
if (!nameStr.compare("eventStandardSel8PbPbQualityCent90")) {
41454193
cut->AddCut(VarManager::kVtxZ, -10.0, 10.0);
41464194
cut->AddCut(VarManager::kIsSel8, 0.5, 1.5);
@@ -4597,6 +4645,12 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
45974645
return cut;
45984646
}
45994647

4648+
if (!nameStr.compare("jpsiStandardKine5")) {
4649+
cut->AddCut(VarManager::kP, 1.0, 1000.0);
4650+
cut->AddCut(VarManager::kEta, -0.9, 0.9);
4651+
return cut;
4652+
}
4653+
46004654
if (!nameStr.compare("jpsiKineSkimmed")) {
46014655
cut->AddCut(VarManager::kPt, 0.7, 1000.0);
46024656
cut->AddCut(VarManager::kEta, -0.9, 0.9);
@@ -5080,6 +5134,12 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
50805134
return cut;
50815135
}
50825136

5137+
if (!nameStr.compare("electronStandardQualityTPCOnly3")) {
5138+
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
5139+
cut->AddCut(VarManager::kTPCncls, 120, 161.);
5140+
return cut;
5141+
}
5142+
50835143
if (!nameStr.compare("NoelectronStandardQualityTPCOnly")) {
50845144
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0, true, VarManager::kTPCncls, 70, 161.);
50855145
return cut;
@@ -5298,6 +5358,21 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
52985358
return cut;
52995359
}
53005360

5361+
if (!nameStr.compare("pidJpsi_TPCpion0")) {
5362+
cut->AddCut(VarManager::kTPCnSigmaPi, 4.0, 1000.0);
5363+
return cut;
5364+
}
5365+
5366+
if (!nameStr.compare("pidJpsi_noTOF_prot")) {
5367+
cut->AddCut(VarManager::kTPCnSigmaPr, 3.5, 1000.0, false, VarManager::kHasTOF, -0.5, 0.5);
5368+
return cut;
5369+
}
5370+
5371+
if (!nameStr.compare("pidJpsi_beta")) {
5372+
cut->AddCut(VarManager::kTOFbeta, 0.98, 1.02, false, VarManager::kHasTOF, 0.5, 1.5);
5373+
return cut;
5374+
}
5375+
53015376
// Magnus cuts ----------------------------------------------------------
53025377

53035378
if (!nameStr.compare("pidJpsi_magnus_ele1")) {

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 39 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1389,6 +1389,45 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
13891389
hm->AddHistogram(histClass, "CosThetaStarMC", "", false, 100, -1.0, 1.0, VarManager::kMCCosThetaStar);
13901390
}
13911391
}
1392+
if (subGroupStr.Contains("flow-jpsi-ep")) {
1393+
int bins_A2[5] = {50, 20, 20, 9, 200};
1394+
double minBins_A2[5] = {2.0, 0.0, -1., 0.0, -20.0};
1395+
double maxBins_A2[5] = {4.0, 2.0, 1.0, 90.0, 20.0};
1396+
int bins_DeltaPhi[5] = {50, 20, 20, 9, 10};
1397+
double minBins_DeltaPhi[5] = {2.0, 0.0, -1., 0.0, 0};
1398+
double maxBins_DeltaPhi[5] = {4.0, 2.0, 1.0, 90.0, 3.14};
1399+
TString labels[5] = {"kMass", "kPt", "kRapidity", "kCentFT0C", "kA2EP"};
1400+
if (subGroupStr.Contains("tpc")) {
1401+
int varA2_TPC_PP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kA2EP_PP_TPC};
1402+
int varA2_TPC_RP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kA2EP_RP_TPC};
1403+
int varDeltaPhi_TPC_PP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kDeltaPhiPP_TPC};
1404+
int varDeltaPhi_TPC_RP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kDeltaPhiRP_TPC};
1405+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_A2PP_TPC", "", 5, varA2_TPC_PP, bins_A2, minBins_A2, maxBins_A2, 0, -1, kTRUE);
1406+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_A2RP_TPC", "", 5, varA2_TPC_RP, bins_A2, minBins_A2, maxBins_A2, 0, -1, kTRUE);
1407+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_DeltaPhiPP_TPC", "", 5, varDeltaPhi_TPC_PP, bins_DeltaPhi, minBins_DeltaPhi, maxBins_DeltaPhi, 0, -1, kTRUE);
1408+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_DeltaPhiRP_TPC", "", 5, varDeltaPhi_TPC_RP, bins_DeltaPhi, minBins_DeltaPhi, maxBins_DeltaPhi, 0, -1, kTRUE);
1409+
}
1410+
if (subGroupStr.Contains("ft0c")) {
1411+
int varA2_FT0C_PP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kA2EP_PP_FT0C};
1412+
int varA2_FT0C_RP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kA2EP_RP_FT0C};
1413+
int varDeltaPhi_FT0C_PP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kDeltaPhiPP_FT0C};
1414+
int varDeltaPhi_FT0C_RP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kDeltaPhiRP_FT0C};
1415+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_A2PP_FT0C", "", 5, varA2_FT0C_PP, bins_A2, minBins_A2, maxBins_A2, 0, -1, kTRUE);
1416+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_A2RP_FT0C", "", 5, varA2_FT0C_RP, bins_A2, minBins_A2, maxBins_A2, 0, -1, kTRUE);
1417+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_DeltaPhiPP_FT0C", "", 5, varDeltaPhi_FT0C_PP, bins_DeltaPhi, minBins_DeltaPhi, maxBins_DeltaPhi, 0, -1, kTRUE);
1418+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_DeltaPhiRP_FT0C", "", 5, varDeltaPhi_FT0C_RP, bins_DeltaPhi, minBins_DeltaPhi, maxBins_DeltaPhi, 0, -1, kTRUE);
1419+
}
1420+
if (subGroupStr.Contains("ft0a")) {
1421+
int varA2_FT0A_PP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kA2EP_PP_FT0A};
1422+
int varA2_FT0A_RP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kA2EP_RP_FT0A};
1423+
int varDeltaPhi_FT0A_PP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kDeltaPhiPP_FT0A};
1424+
int varDeltaPhi_FT0A_RP[5] = {VarManager::kMass, VarManager::kPt, VarManager::kRap, VarManager::kCentFT0C, VarManager::kDeltaPhiRP_FT0A};
1425+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_A2PP_FT0A", "", 5, varA2_FT0A_PP, bins_A2, minBins_A2, maxBins_A2, 0, -1, kTRUE);
1426+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_A2RP_FT0A", "", 5, varA2_FT0A_RP, bins_A2, minBins_A2, maxBins_A2, 0, -1, kTRUE);
1427+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_DeltaPhiPP_FT0A", "", 5, varDeltaPhi_FT0A_PP, bins_DeltaPhi, minBins_DeltaPhi, maxBins_DeltaPhi, 0, -1, kTRUE);
1428+
hm->AddHistogram(histClass, "Mass_Pt_centrFT0C_DeltaPhiRP_FT0A", "", 5, varDeltaPhi_FT0A_RP, bins_DeltaPhi, minBins_DeltaPhi, maxBins_DeltaPhi, 0, -1, kTRUE);
1429+
}
1430+
}
13921431
if (subGroupStr.Contains("upsilon")) {
13931432
hm->AddHistogram(histClass, "MassUpsilon_Pt", "", false, 500, 7.0, 12.0, VarManager::kMass, 400, 0.0, 40.0, VarManager::kPt);
13941433
}

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2485,6 +2485,12 @@ void VarManager::SetDefaultVarNames()
24852485
fgVarNamesMap["kDCATrackVtxProd"] = kDCATrackVtxProd;
24862486
fgVarNamesMap["kV2SP"] = kV2SP;
24872487
fgVarNamesMap["kV2EP"] = kV2EP;
2488+
fgVarNamesMap["kA2EP_PP_TPC"] = kA2EP_PP_TPC;
2489+
fgVarNamesMap["kA2EP_PP_FT0A"] = kA2EP_PP_FT0A;
2490+
fgVarNamesMap["kA2EP_PP_FT0C"] = kA2EP_PP_FT0C;
2491+
fgVarNamesMap["kA2EP_RP_TPC"] = kA2EP_RP_TPC;
2492+
fgVarNamesMap["kA2EP_RP_FT0A"] = kA2EP_RP_FT0A;
2493+
fgVarNamesMap["kA2EP_RP_FT0C"] = kA2EP_RP_FT0C;
24882494
fgVarNamesMap["kWV2SP"] = kWV2SP;
24892495
fgVarNamesMap["kWV2EP"] = kWV2EP;
24902496
fgVarNamesMap["kU2Q2"] = kU2Q2;

0 commit comments

Comments
 (0)