Skip to content

Commit 89cdbf9

Browse files
ypwanggypwangg
committed
Coherent Jpsi A2 definition using event plane at midrapidity, same event only
1 parent 21e9287 commit 89cdbf9

6 files changed

Lines changed: 551 additions & 7 deletions

File tree

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 48 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2973,6 +2973,17 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
29732973
return cut;
29742974
}
29752975

2976+
if (!nameStr.compare("jpsi_debug_TPCTOF3_rejBadTOF")) {
2977+
cut->AddCut(GetAnalysisCut("jpsiStandardKine5"));
2978+
cut->AddCut(GetAnalysisCut("electronStandardQualityTPCOnly3"));
2979+
cut->AddCut(GetAnalysisCut("SPDfirst"));
2980+
cut->AddCut(GetAnalysisCut("dcaCut1_ionut"));
2981+
cut->AddCut(GetAnalysisCut("pidJpsi_TPCpion0"));
2982+
cut->AddCut(GetAnalysisCut("pidJpsi_beta"));
2983+
cut->AddCut(GetAnalysisCut("pidJpsi_noTOF_prot"));
2984+
return cut;
2985+
}
2986+
29762987
// -------------------------------------------------------------------------------------------------
29772988
// lmee pair cuts
29782989

@@ -4113,6 +4124,15 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
41134124
return cut;
41144125
}
41154126

4127+
if (!nameStr.compare("eventStandardSel8NoPileup")) {
4128+
cut->AddCut(VarManager::kVtxZ, -10.0, 10.0);
4129+
cut->AddCut(VarManager::kIsSel8, 0.5, 1.5);
4130+
cut->AddCut(VarManager::kIsNoSameBunch, 0.5, 1.5);
4131+
cut->AddCut(VarManager::kIsGoodZvtxFT0vsPV, 0.5, 1.5);
4132+
cut->AddCut(VarManager::kNoCollInTimeRangeStandard, 0.5, 1.5);
4133+
return cut;
4134+
}
4135+
41164136
if (!nameStr.compare("eventStandardSel8PbPbQualityCent90")) {
41174137
cut->AddCut(VarManager::kVtxZ, -10.0, 10.0);
41184138
cut->AddCut(VarManager::kIsSel8, 0.5, 1.5);
@@ -4569,6 +4589,12 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
45694589
return cut;
45704590
}
45714591

4592+
if (!nameStr.compare("jpsiStandardKine5")) {
4593+
cut->AddCut(VarManager::kP, 1.0, 1000.0);
4594+
cut->AddCut(VarManager::kEta, -0.9, 0.9);
4595+
return cut;
4596+
}
4597+
45724598
if (!nameStr.compare("jpsiKineSkimmed")) {
45734599
cut->AddCut(VarManager::kPt, 0.7, 1000.0);
45744600
cut->AddCut(VarManager::kEta, -0.9, 0.9);
@@ -5052,6 +5078,13 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
50525078
return cut;
50535079
}
50545080

5081+
if (!nameStr.compare("electronStandardQualityTPCOnly3")) {
5082+
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0);
5083+
cut->AddCut(VarManager::kTPCncls, 120, 161.);
5084+
return cut;
5085+
}
5086+
5087+
50555088
if (!nameStr.compare("NoelectronStandardQualityTPCOnly")) {
50565089
cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0, true, VarManager::kTPCncls, 70, 161.);
50575090
return cut;
@@ -5270,6 +5303,21 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
52705303
return cut;
52715304
}
52725305

5306+
if (!nameStr.compare("pidJpsi_TPCpion0")) {
5307+
cut->AddCut(VarManager::kTPCnSigmaPi, 4.0, 1000.0);
5308+
return cut;
5309+
}
5310+
5311+
if (!nameStr.compare("pidJpsi_noTOF_prot")) {
5312+
cut->AddCut(VarManager::kTPCnSigmaPr, 3.5, 1000.0, false, VarManager::kHasTOF, -0.5, 0.5);
5313+
return cut;
5314+
}
5315+
5316+
if (!nameStr.compare("pidJpsi_beta")) {
5317+
cut->AddCut(VarManager::kTOFbeta, 0.98, 1.02, false, VarManager::kHasTOF, 0.5, 1.5);
5318+
return cut;
5319+
}
5320+
52735321
// Magnus cuts ----------------------------------------------------------
52745322

52755323
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: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -76,6 +76,8 @@ int VarManager::fgCalibrationType = 0; // 0 - no calibration, 1 -
7676
bool VarManager::fgUseInterpolatedCalibration = true; // use interpolated calibration histograms (default: true)
7777
int VarManager::fgEfficiencyType = 0; // type of efficiency to be applied, default is no efficiency
7878
TObject* VarManager::fgEfficiencyHist = nullptr; // histogram for efficiency
79+
TH3F* VarManager::fgObjQvec = nullptr;
80+
bool VarManager::fgApplyQVectorCorrection = false;
7981

8082
//__________________________________________________________________
8183
VarManager::VarManager() : TObject()
@@ -2484,6 +2486,12 @@ void VarManager::SetDefaultVarNames()
24842486
fgVarNamesMap["kDCATrackVtxProd"] = kDCATrackVtxProd;
24852487
fgVarNamesMap["kV2SP"] = kV2SP;
24862488
fgVarNamesMap["kV2EP"] = kV2EP;
2489+
fgVarNamesMap["kA2EP_PP_TPC"] = kA2EP_PP_TPC;
2490+
fgVarNamesMap["kA2EP_PP_FT0A"] = kA2EP_PP_FT0A;
2491+
fgVarNamesMap["kA2EP_PP_FT0C"] = kA2EP_PP_FT0C;
2492+
fgVarNamesMap["kA2EP_RP_TPC"] = kA2EP_RP_TPC;
2493+
fgVarNamesMap["kA2EP_RP_FT0A"] = kA2EP_RP_FT0A;
2494+
fgVarNamesMap["kA2EP_RP_FT0C"] = kA2EP_RP_FT0C;
24872495
fgVarNamesMap["kWV2SP"] = kWV2SP;
24882496
fgVarNamesMap["kWV2EP"] = kWV2EP;
24892497
fgVarNamesMap["kU2Q2"] = kU2Q2;

0 commit comments

Comments
 (0)