Skip to content

Commit d88275c

Browse files
axti98axti98
committed
Added additional TOF only cut, added alice 3 track histograms.
1 parent 85604b4 commit d88275c

File tree

3 files changed

+120
-6
lines changed

3 files changed

+120
-6
lines changed

‎PWGDQ/Core/CutsLibrary.cxx‎

Lines changed: 75 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -3529,6 +3529,11 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
35293529
return cut;
35303530
}
35313531

3532+
if (!nameStr.compare("pairD0")) {
3533+
cut->AddCut(GetAnalysisCut("pairD0"));
3534+
return cut;
3535+
}
3536+
35323537
if (!nameStr.compare("pairD0HighPt1")) {
35333538
cut->AddCut(GetAnalysisCut("pairLxyzProjected3sigma"));
35343539
cut->AddCut(GetAnalysisCut("pairPtLow5"));
@@ -3881,7 +3886,49 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
38813886
}
38823887

38833888
if (!nameStr.compare("alice3DielectronPID")) {
3884-
cut->AddCut(GetAnalysisCut("alice3StandardKine"));
3889+
cut->AddCut(GetAnalysisCut("alice3JpsiKine"));
3890+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3891+
cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
3892+
cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
3893+
cut->AddCut(GetAnalysisCut("alice3RICHPIDEl"));
3894+
return cut;
3895+
}
3896+
3897+
if (!nameStr.compare("alice3DielectronPIDTOFOnly")) {
3898+
cut->AddCut(GetAnalysisCut("alice3JpsiKineTOFAcceptance"));
3899+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3900+
cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
3901+
cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
3902+
return cut;
3903+
}
3904+
3905+
if (!nameStr.compare("alice3DielectronPIDRICHOnly")) {
3906+
cut->AddCut(GetAnalysisCut("alice3JpsiKineTOFAcceptance"));
3907+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3908+
cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
3909+
cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
3910+
return cut;
3911+
}
3912+
3913+
if (!nameStr.compare("alice3DielectronPIDTOFOnly")) {
3914+
cut->AddCut(GetAnalysisCut("alice3JpsiKine"));
3915+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3916+
cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
3917+
cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
3918+
return cut;
3919+
}
3920+
3921+
if (!nameStr.compare("alice3DielectronPIDTOFAcceptance")) {
3922+
cut->AddCut(GetAnalysisCut("alice3JpsiKineTOFAcceptance"));
3923+
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
3924+
cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
3925+
cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
3926+
cut->AddCut(GetAnalysisCut("alice3RICHPIDEl"));
3927+
return cut;
3928+
}
3929+
3930+
if (!nameStr.compare("alice3DielectronPIDRICHAcceptance")) {
3931+
cut->AddCut(GetAnalysisCut("alice3JpsiKineRICHAcceptance"));
38853932
cut->AddCut(GetAnalysisCut("alice3TrackQuality"));
38863933
cut->AddCut(GetAnalysisCut("alice3iTOFPIDEl"));
38873934
cut->AddCut(GetAnalysisCut("alice3oTOFPIDEl"));
@@ -6732,6 +6779,11 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
67326779
return cut;
67336780
}
67346781

6782+
if (!nameStr.compare("pairD0")) {
6783+
cut->AddCut(VarManager::kMass, 1.814, 1.914);
6784+
return cut;
6785+
}
6786+
67356787
if (!nameStr.compare("pairJpsi")) {
67366788
cut->AddCut(VarManager::kMass, 2.8, 3.3);
67376789
return cut;
@@ -6975,13 +7027,35 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
69757027
return cut;
69767028
}
69777029

7030+
if (!nameStr.compare("alice3JpsiKine")) {
7031+
cut->AddCut(VarManager::kPt, 1.0, 1000.0);
7032+
cut->AddCut(VarManager::kEta, -2.5, 2.5); // Total tracker acceptance in v3b geomety
7033+
return cut;
7034+
}
7035+
7036+
if (!nameStr.compare("alice3JpsiKineTOFAcceptance")) {
7037+
cut->AddCut(VarManager::kPt, 1.0, 1000.0);
7038+
cut->AddCut(VarManager::kEta, -2.0, 2.0); // TOF acceptance in v3b geomety
7039+
return cut;
7040+
}
7041+
7042+
if (!nameStr.compare("alice3JpsiKineRICHAcceptance")) {
7043+
cut->AddCut(VarManager::kPt, 1.0, 1000.0);
7044+
cut->AddCut(VarManager::kEta, -0.8, 0.8); // RICH acceptance in v3b geomety
7045+
return cut;
7046+
}
7047+
69787048
if (!nameStr.compare("alice3TrackQuality")) {
7049+
cut->AddCut(VarManager::kIsReconstructed, 0.5, 1.5);
7050+
cut->AddCut(VarManager::kNSiliconHits, 6.0, 12.0);
69797051
cut->AddCut(VarManager::kTrackDCAxy, -3.0, 3.0);
69807052
cut->AddCut(VarManager::kTrackDCAz, -3.0, 3.0);
69817053
return cut;
69827054
}
69837055

69847056
if (!nameStr.compare("alice3TrackQualityTightDCA")) {
7057+
cut->AddCut(VarManager::kIsReconstructed, 0.5, 1.5);
7058+
cut->AddCut(VarManager::kNSiliconHits, 6.0, 12.0);
69857059
cut->AddCut(VarManager::kTrackDCAxy, -1.0, 1.0);
69867060
cut->AddCut(VarManager::kTrackDCAz, -1.0, 1.0);
69877061
return cut;

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 42 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -1041,19 +1041,21 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
10411041
hm->AddHistogram(histClass, "Rapidity", "", false, 400, -4.0, 4.0, VarManager::kRap);
10421042
}
10431043
if (subGroupStr.Contains("alice3tracks")) {
1044-
hm->AddHistogram(histClass, "Pt_vs_PtMC", "pT vs MC pT", false, 200, 0.0, 20.0, VarManager::kPt, 200, 0.0, 20.0, VarManager::kMCPt);
1044+
hm->AddHistogram(histClass, "IsReconstructed", "IsReconstructed", false, 2, 0.5, 1.5, VarManager::kIsReconstructed);
1045+
hm->AddHistogram(histClass, "nSiliconHits", "Num silicon hits", false, 12, 0.0, 12.0, VarManager::kNSiliconHits);
1046+
hm->AddHistogram(histClass, "Pt_vs_PtMC", "pT vs MC pT", false, 200, 0.0, 100.0, VarManager::kPt, 200, 0.0, 20.0, VarManager::kMCPt);
10451047
hm->AddHistogram(histClass, "Eta_vs_EtaMC", "#eta vs MC #eta", false, 150, -4.0, 4.0, VarManager::kEta, 150, -4.0, 4.0, VarManager::kMCEta);
10461048
hm->AddHistogram(histClass, "Phi_vs_PhiMC", "#varphi vs MC #varphi", false, 50, 0.0, 2. * o2::constants::math::PI, VarManager::kPhi, 50, 0.0, 2. * o2::constants::math::PI, VarManager::kMCPhi);
10471049
hm->AddHistogram(histClass, "Eta_DCAxy", "#eta vs DCA_{xy}", false, 80, -4.0, 4.0, VarManager::kEta, 400, -2.0, 2.0, VarManager::kTrackDCAxy);
10481050
hm->AddHistogram(histClass, "Eta_DCAz", "#eta vs DCA_{z}", false, 80, -4.0, 4.0, VarManager::kEta, 800, -4.0, 4.0, VarManager::kTrackDCAz);
10491051
}
10501052
if (subGroupStr.Contains("resolutions")) {
1051-
hm->AddHistogram(histClass, "PtMC_vs_DeltaPt", "MC pT vs DeltaPt", false, 200, 0.0, 20.0, VarManager::kMCPt, 2000, -1.0, 1.0, VarManager::kDeltaPt);
1052-
hm->AddHistogram(histClass, "Pt_vs_DeltaPt", "pT vs pTRes", false, 200, 0.0, 20.0, VarManager::kPt, 2000, -1.0, 1.0, VarManager::kDeltaPt);
1053+
hm->AddHistogram(histClass, "PtMC_vs_DeltaPt", "MC pT vs DeltaPt", false, 200, 0.0, 100.0, VarManager::kMCPt, 2000, -1.0, 1.0, VarManager::kDeltaPt);
1054+
hm->AddHistogram(histClass, "Pt_vs_DeltaPt", "pT vs pTRes", false, 200, 0.0, 100.0, VarManager::kPt, 2000, -1.0, 1.0, VarManager::kDeltaPt);
10531055
hm->AddHistogram(histClass, "EtaMC_vs_DeltaPt", "MC Eta vs DeltaPt", false, 200, -6.0, 6.0, VarManager::kMCEta, 2000, -1.0, 1.0, VarManager::kDeltaPt);
10541056
hm->AddHistogram(histClass, "Eta_vs_DeltaPt", "MC Eta vs DeltaPt", false, 200, -6.0, 6.0, VarManager::kEta, 2000, -1.0, 1.0, VarManager::kDeltaPt);
1055-
hm->AddHistogram(histClass, "PtMC_vs_Res", "MC pT vs DeltaPt", false, 200, 0.0, 20.0, VarManager::kMCPt, 2000, -1.0, 1.0, VarManager::kPtResolution);
1056-
hm->AddHistogram(histClass, "Pt_vs_Res", "pT vs pTRes", false, 200, 0.0, 20.0, VarManager::kPt, 2000, -1.0, 1.0, VarManager::kPtResolution);
1057+
hm->AddHistogram(histClass, "PtMC_vs_Res", "MC pT vs DeltaPt", false, 200, 0.0, 100.0, VarManager::kMCPt, 2000, -1.0, 1.0, VarManager::kPtResolution);
1058+
hm->AddHistogram(histClass, "Pt_vs_Res", "pT vs pTRes", false, 200, 0.0, 100.0, VarManager::kPt, 2000, -1.0, 1.0, VarManager::kPtResolution);
10571059
hm->AddHistogram(histClass, "EtaMC_vs_Res", "MC Eta vs pTRes", false, 200, -6.0, 6.0, VarManager::kMCEta, 2000, -1.0, 1.0, VarManager::kPtResolution);
10581060
hm->AddHistogram(histClass, "Eta_vs_Res", "MC Eta vs pTRes", false, 200, -6.0, 6.0, VarManager::kEta, 2000, -1.0, 1.0, VarManager::kPtResolution);
10591061
hm->AddHistogram(histClass, "PtRes", "pT Resolution", false, 2000, -1.0, 1.0, VarManager::kPtResolution);
@@ -1207,6 +1209,29 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
12071209
hm->AddHistogram(histClass, "Mass_MultFDDA", "Mass vs MultFDDA", false, 200, 2.0, 5.0, VarManager::kMass, 1000, 0, 25000.0, VarManager::kMultFDDA);
12081210
hm->AddHistogram(histClass, "Mass_MultFDDC", "Mass vs MultFDDC", false, 200, 2.0, 5.0, VarManager::kMass, 1000, 0, 25000.0, VarManager::kMultFDDC);
12091211
}
1212+
if (subGroupStr.Contains("mass_mult_pv")) {
1213+
double multAnalysisBins[4] = {0.0, 4.0, 8.0, 1000.0};
1214+
1215+
double multAnalsisMassBins[201] = {0.0};
1216+
for (int i = 0; i <= 200; i++) {
1217+
multAnalsisMassBins[i] = 2.0 + i * 0.015;
1218+
}
1219+
1220+
double multAnalysisPVBins[151] = {0.0};
1221+
for (int i = 0; i <= 150; i++) {
1222+
multAnalysisPVBins[i] = i * 1.0;
1223+
}
1224+
1225+
double multAnalysisFV0ABins[1001] = {0.0};
1226+
for (int i = 0; i <= 1000; i++) {
1227+
multAnalysisFV0ABins[i] = i * 1.0;
1228+
}
1229+
1230+
hm->AddHistogram(histClass, "Mass_VtxNcontribReal_Pt", "Mass vs VtxNcontribReal vs Pt", false, 200, multAnalsisMassBins, VarManager::kMass, 150, multAnalysisPVBins, VarManager::kVtxNcontribReal, 3, multAnalysisBins, VarManager::kPt);
1231+
hm->AddHistogram(histClass, "Mass_MultFV0A_Pt", "Mass vs MultFV0A vs Pt", false, 200, multAnalsisMassBins, VarManager::kMass, 1000, multAnalysisFV0ABins, VarManager::kMultFV0A, 3, multAnalysisBins, VarManager::kPt);
1232+
hm->AddHistogram(histClass, "Mass_MultFT0A_Pt", "Mass vs MultFT0A vs Pt", false, 200, multAnalsisMassBins, VarManager::kMass, 1000, multAnalysisFV0ABins, VarManager::kMultFT0A, 3, multAnalysisBins, VarManager::kPt);
1233+
hm->AddHistogram(histClass, "Mass_MultFT0C_Pt", "Mass vs MultFT0C vs Pt", false, 200, multAnalsisMassBins, VarManager::kMass, 1000, multAnalysisFV0ABins, VarManager::kMultFT0C, 3, multAnalysisBins, VarManager::kPt);
1234+
}
12101235
if (subGroupStr.Contains("barrel")) {
12111236
hm->AddHistogram(histClass, "Mass", "", false, 500, 0.0, 5.0, VarManager::kMass);
12121237
hm->AddHistogram(histClass, "Mass_HighRange", "", false, 375, 0.0, 15.0, VarManager::kMass);
@@ -1241,6 +1266,18 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
12411266
hm->AddHistogram(histClass, "Mass_VtxNcontribReal_VtxZ", "Mass vs VtxNcontribReal vs VtxZ", false, 200, 2.0, 5.0, VarManager::kMass, 150, 0, 150.0, VarManager::kVtxNcontribReal, 20, -10.0, 10.0, VarManager::kVtxZ);
12421267
hm->AddHistogram(histClass, "VtxZ_VtxNcontribReal", "VtxZ vs VtxNcontribReal", false, 240, -12.0, 12.0, VarManager::kVtxZ, 200, 0, 200.0, VarManager::kVtxNcontribReal);
12431268
}
1269+
if (subGroupStr.Contains("alice3barrel")) {
1270+
double pTBins[21] = {
1271+
0.1, 0.5, 1.0, 1.25, 1.5, 2.0, 2.5, 3.0,
1272+
4.0, 5.0, 6.5, 8.0, 10.0, 12.0, 15.0, 20.0,
1273+
30.0, 40.0, 60.0, 80.0, 100.0};
1274+
1275+
double massBins[501];
1276+
for (int i = 0; i < 501; i++) {
1277+
massBins[i] = 0.01 * (i);
1278+
}
1279+
hm->AddHistogram(histClass, "Mass_PtLong", "", false, 500, massBins, VarManager::kMass, 20, pTBins, VarManager::kPt);
1280+
}
12441281
if (subGroupStr.Contains("dielectron-polarization-he-pbpb")) {
12451282
int varsHEpbpb[5] = {VarManager::kMass, VarManager::kPt, VarManager::kCentFT0C, VarManager::kCosThetaHE, VarManager::kPhiHE};
12461283
int binspT[5] = {100, 30, 10, 10, 10};

‎PWGDQ/Core/VarManager.h‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -70,6 +70,7 @@
7070
#include <cstdint>
7171
#include <iostream>
7272
#include <map>
73+
#include <tuple>
7374
#include <utility>
7475
#include <vector>
7576

@@ -6618,6 +6619,8 @@ void VarManager::FillTrackAlice3(T const& track, float* values)
66186619
values[kITSClusterMap] = track.itsClusterMap();
66196620
}
66206621

6622+
values[kIsReconstructed] = track.isReconstructed();
6623+
values[kNSiliconHits] = track.nSiliconHits();
66216624
values[kITSchi2] = track.itsChi2NCl();
66226625
}
66236626

0 commit comments

Comments
 (0)