Skip to content

Commit 6d1f085

Browse files
YazhenLinYazhenLin
authored
Add files via upload
1 parent 8af0b91 commit 6d1f085

File tree

2 files changed

+57
-78
lines changed

2 files changed

+57
-78
lines changed

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 42 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -906,42 +906,42 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
906906
}
907907

908908
if (!groupStr.CompareTo("mctruth_triple")) {
909-
hm->AddHistogram(histClass, "Eta_Pt", "", false, 100, -2.0, 2.0, VarManager::kPairEta, 200, 0.0, 30.0, VarManager::kPairPt);
910-
hm->AddHistogram(histClass, "Eta_Pt_lepton1", "", false, 100, -2.0, 2.0, VarManager::kEta1, 200, 0.0, 30.0, VarManager::kPt1);
911-
hm->AddHistogram(histClass, "Eta_Pt_lepton2", "", false, 100, -2.0, 2.0, VarManager::kEta2, 200, 0.0, 30.0, VarManager::kPt2);
912-
hm->AddHistogram(histClass, "Eta_Pt_Assoc", "", false, 100, -2.0, 2.0, VarManager::kEta, 200, 0.0, 30.0, VarManager::kPt);
909+
hm->AddHistogram(histClass, "Eta_Pt", "", false, 100, -2.0, 2.0, VarManager::kPairEta, 200, 0.0, 20.0, VarManager::kPairPt);
910+
hm->AddHistogram(histClass, "Eta_Pt_lepton1", "", false, 100, -2.0, 2.0, VarManager::kEta1, 200, 0.0, 20.0, VarManager::kPt1);
911+
hm->AddHistogram(histClass, "Eta_Pt_lepton2", "", false, 100, -2.0, 2.0, VarManager::kEta2, 200, 0.0, 20.0, VarManager::kPt2);
912+
hm->AddHistogram(histClass, "Eta_Pt_Photon", "", false, 100, -2.0, 2.0, VarManager::kEta, 200, 0.0, 20.0, VarManager::kPt);
913913
hm->AddHistogram(histClass, "Phi_Eta", "#phi vs #eta distribution", false, 200, -5.0, 5.0, VarManager::kPairEta, 200, -2. * o2::constants::math::PI, 2. * o2::constants::math::PI, VarManager::kPairPhi);
914914
hm->AddHistogram(histClass, "Mass_Dilepton", "", false, 4500, 0.0, 4.5, VarManager::kPairMassDau);
915-
hm->AddHistogram(histClass, "Mass_Assoc", "", false, 500, 0.0, 1.0, VarManager::kMassDau);
916-
hm->AddHistogram(histClass, "Mass_Dilepton_Mass_Assoc", "", false, 500, 0.0, 5.0, VarManager::kPairMassDau, 500, 0.0, 5.0, VarManager::kMassDau);
917-
hm->AddHistogram(histClass, "Pt_Dilepton", "", false, 2000, 0.0, 30.0, VarManager::kPairPtDau);
918-
hm->AddHistogram(histClass, "Pt_Assoc", "", false, 500, 0.0, 10.0, VarManager::kPt);
919-
hm->AddHistogram(histClass, "Mass_DileptonAssoc", "", false, 4500, 0.0, 30.0, VarManager::kPairMass);
920-
hm->AddHistogram(histClass, "Pt_DileptonAssoc", "", false, 2000, 0.0, 30.0, VarManager::kPairPt);
921-
hm->AddHistogram(histClass, "Mass_Pt_DileptonAssoc", "", false, 500, 0.0, 30.0, VarManager::kPairMass, 200, 0.0, 30.0, VarManager::kPairPt);
922-
hm->AddHistogram(histClass, "Pt_Rap_DileptonAssoc", "", false, 5120, 0.0, 30.0, VarManager::kPairPt, 1000, -5.0, 5.0, VarManager::kRap);
923-
hm->AddHistogram(histClass, "DeltaMass", "", false, 1500, 0.0, 5.5, VarManager::kDeltaMass);
924-
hm->AddHistogram(histClass, "DeltaMass_ptdileptonassoc", "", false, 1000, 0.0, 1.0, VarManager::kDeltaMass, 3000, 0.0, 30.0, VarManager::kPairPt);
915+
hm->AddHistogram(histClass, "Mass_Photon", "", false, 500, 0.0, 0.1, VarManager::kMassDau);
916+
hm->AddHistogram(histClass, "Mass_Dilepton_Mass_Photon", "", false, 500, 0.0, 5.0, VarManager::kPairMassDau, 500, 0.0, 5.0, VarManager::kMassDau);
917+
hm->AddHistogram(histClass, "Pt_Dilepton", "", false, 2000, 0.0, 20.0, VarManager::kPairPtDau);
918+
hm->AddHistogram(histClass, "Pt_Photon", "", false, 500, 0.0, 5.0, VarManager::kPt);
919+
hm->AddHistogram(histClass, "Mass_DileptonPhoton", "", false, 4500, 0.0, 4.5, VarManager::kPairMass);
920+
hm->AddHistogram(histClass, "Pt_DileptonPhoton", "", false, 2000, 0.0, 20.0, VarManager::kPairPt);
921+
hm->AddHistogram(histClass, "Mass_Pt_DileptonPhoton", "", false, 500, 0.0, 5.0, VarManager::kPairMass, 200, 0.0, 20.0, VarManager::kPairPt);
922+
hm->AddHistogram(histClass, "DeltaMass", "", false, 1500, 0.0, 1.5, VarManager::kDeltaMass);
923+
hm->AddHistogram(histClass, "DeltaMass_ptdileptonphoton", "", false, 1000, 0.0, 1.0, VarManager::kDeltaMass, 3000, 0.0, 30.0, VarManager::kPairPt);
925924
hm->AddHistogram(histClass, "DeltaMass_Jpsi", "", false, 1500, 3, 4.5, (VarManager::kDeltaMass_jpsi));
926-
hm->AddHistogram(histClass, "Rapidity", "", false, 400, -5.0, 5.0, VarManager::kRap);
925+
hm->AddHistogram(histClass, "Rapidity", "", false, 400, -4.0, 4.0, VarManager::kRap);
927926
}
928927
if (!groupStr.CompareTo("mctruth_pair")) {
929-
hm->AddHistogram(histClass, "Mass_Pt", "", false, 500, 0.0, 15.0, VarManager::kMCMass, 40, 0.0, 20.0, VarManager::kMCPt);
930-
hm->AddHistogram(histClass, "Pt", "", false, 200, 0.0, 20.0, VarManager::kMCPt);
928+
hm->AddHistogram(histClass, "Mass_Pt", "", false, 500, 0.0, 15.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt);
929+
hm->AddHistogram(histClass, "Pt", "", false, 200, 0.0, 20.0, VarManager::kPt);
931930
hm->AddHistogram(histClass, "Pt_Dilepton", "", false, 200, 0.0, 20.0, VarManager::kPairPtDau);
932-
hm->AddHistogram(histClass, "Mass", "", false, 500, 0.0, 15.0, VarManager::kMCMass);
933-
hm->AddHistogram(histClass, "Rapidity", "", false, 100, -5.0, 5.0, VarManager::kMCY);
934-
hm->AddHistogram(histClass, "Eta_Pt", "", false, 40, -2.0, 2.0, VarManager::kMCEta, 200, 0.0, 20.0, VarManager::kMCPt);
935-
hm->AddHistogram(histClass, "Phi_Eta", "#phi vs #eta distribution", false, 200, -5.0, 5.0, VarManager::kMCEta, 200, -2. * o2::constants::math::PI, 2. * o2::constants::math::PI, VarManager::kMCPhi);
936-
if (subGroupStr.Contains("polarization")) {
937-
int varspTHE[4] = {VarManager::kMCPt, VarManager::kMCCosThetaHE, VarManager::kMCPhiHE, VarManager::kMCPhiTildeHE};
938-
int varspTCS[4] = {VarManager::kMCPt, VarManager::kMCCosThetaCS, VarManager::kMCPhiCS, VarManager::kMCPhiTildeCS};
939-
int bins[4] = {20, 20, 20, 20};
940-
double xmin[4] = {0., -1., 0., 0.};
941-
double xmax[4] = {20., 1., 2. * o2::constants::math::PI, 2. * o2::constants::math::PI};
942-
hm->AddHistogram(histClass, "Pt_cosThetaHE_phiHE_phiTildeHE", "", 4, varspTHE, bins, xmin, xmax, 0, -1, kFALSE);
943-
hm->AddHistogram(histClass, "Pt_cosThetaCS_phiCS_phiTildeCS", "", 4, varspTCS, bins, xmin, xmax, 0, -1, kFALSE);
944-
}
931+
hm->AddHistogram(histClass, "Eta_Pt_lepton1", "", false, 100, -2.0, 2.0, VarManager::kEta1, 200, 0.0, 20.0, VarManager::kPt1);
932+
hm->AddHistogram(histClass, "Eta_Pt_lepton2", "", false, 100, -2.0, 2.0, VarManager::kEta2, 200, 0.0, 20.0, VarManager::kPt2);
933+
hm->AddHistogram(histClass, "Mass", "", false, 500, 0.0, 15.0, VarManager::kMass);
934+
hm->AddHistogram(histClass, "Eta_Pt", "", false, 40, -2.0, 2.0, VarManager::kEta, 200, 0.0, 20.0, VarManager::kPt);
935+
hm->AddHistogram(histClass, "Phi_Eta", "#phi vs #eta distribution", false, 200, -5.0, 5.0, VarManager::kEta, 200, -2. * o2::constants::math::PI, 2. * o2::constants::math::PI, VarManager::kPhi);
936+
int varspTHE[3] = {VarManager::kMCPt, VarManager::kMCCosThetaHE, VarManager::kMCPhiHE};
937+
int varspTCS[3] = {VarManager::kMCPt, VarManager::kMCCosThetaCS, VarManager::kMCPhiCS};
938+
int varspTPP[3] = {VarManager::kMCPt, VarManager::kMCCosThetaPP, VarManager::kMCPhiPP};
939+
int binspT[3] = {40, 20, 20};
940+
double xminpT[3] = {0., -1., -3.14};
941+
double xmaxpT[3] = {20., 1., +3.14};
942+
hm->AddHistogram(histClass, "Pt_cosThetaHE_phiHE", "", 3, varspTHE, binspT, xminpT, xmaxpT, 0, -1, kFALSE);
943+
hm->AddHistogram(histClass, "Pt_cosThetaCS_phiCS", "", 3, varspTCS, binspT, xminpT, xmaxpT, 0, -1, kFALSE);
944+
hm->AddHistogram(histClass, "Pt_cosThetaPP_phiPP", "", 3, varspTPP, binspT, xminpT, xmaxpT, 0, -1, kFALSE);
945945
}
946946
if (!groupStr.CompareTo("mctruth_quad")) {
947947
hm->AddHistogram(histClass, "hMass_defaultDileptonMass", "", false, 1000, 3.0, 5.0, VarManager::kQuadDefaultDileptonMass);
@@ -1052,8 +1052,10 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
10521052
hm->AddHistogram(histClass, "Eta_Pt", "", false, 40, -2.0, 2.0, VarManager::kEta, 40, 0.0, 20.0, VarManager::kPt);
10531053
hm->AddHistogram(histClass, "Y_Pt", "", false, 40, -2.0, 2.0, VarManager::kRap, 40, 0.0, 20.0, VarManager::kPt);
10541054
hm->AddHistogram(histClass, "Mass_VtxZ", "", true, 30, -15.0, 15.0, VarManager::kVtxZ, 500, 0.0, 5.0, VarManager::kMass);
1055+
if (subGroupStr.Contains("energy-correlator")) {
10551056
hm->AddHistogram(histClass, "Mass_Y_Pt", "", false, 125, 0.0, 5.0, VarManager::kMass, 40, -2.0, 2.0, VarManager::kRap, 40, 0.0, 20.0, VarManager::kPt);
10561057
hm->AddHistogram(histClass, "Y", ";y", false, 40, -2.0, 2.0, VarManager::kRap);
1058+
}
10571059
if (subGroupStr.Contains("pbpb")) {
10581060
hm->AddHistogram(histClass, "Mass_CentFT0C", "", false, 125, 0.0, 5.0, VarManager::kMass, 20, 0.0, 100.0, VarManager::kCentFT0C);
10591061
hm->AddHistogram(histClass, "Pt_CentFT0C", "", false, 100, 0.0, 10.0, VarManager::kPt, 20, 0.0, 100.0, VarManager::kCentFT0C);
@@ -1082,21 +1084,31 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
10821084
}
10831085
if (subGroupStr.Contains("polarization")) {
10841086
if (subGroupStr.Contains("helicity")) {
1087+
hm->AddHistogram(histClass, "cosThetaHE", "", false, 100, -1., 1., VarManager::kCosThetaHE);
1088+
hm->AddHistogram(histClass, "phiHE", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiHE);
1089+
hm->AddHistogram(histClass, "phitildeHE", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiTildeHE);
10851090
hm->AddHistogram(histClass, "Mass_Pt_CosThetaHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaHE);
10861091
hm->AddHistogram(histClass, "Mass_Pt_PhiHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiHE);
10871092
hm->AddHistogram(histClass, "Mass_Pt_PhiTildeHE", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiTildeHE);
10881093
}
10891094
if (subGroupStr.Contains("collins-soper")) {
1095+
hm->AddHistogram(histClass, "cosThetaCS", "", false, 100, -1., 1., VarManager::kCosThetaCS);
1096+
hm->AddHistogram(histClass, "phiCS", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiCS);
1097+
hm->AddHistogram(histClass, "phitildeCS", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiTildeCS);
10901098
hm->AddHistogram(histClass, "Mass_Pt_CosThetaCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaCS);
10911099
hm->AddHistogram(histClass, "Mass_Pt_PhiCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiCS);
10921100
hm->AddHistogram(histClass, "Mass_Pt_PhiTildeCS", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiTildeCS);
10931101
}
10941102
if (subGroupStr.Contains("production")) {
1103+
hm->AddHistogram(histClass, "cosThetaPP", "", false, 100, -1., 1., VarManager::kCosThetaPP);
1104+
hm->AddHistogram(histClass, "phiPP", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiPP);
1105+
hm->AddHistogram(histClass, "phitildePP", "", false, 100, 0, 2 * o2::constants::math::PI, VarManager::kPhiTildePP);
10951106
hm->AddHistogram(histClass, "Mass_Pt_CosThetaPP", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaPP);
10961107
hm->AddHistogram(histClass, "Mass_Pt_PhiPP", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiPP);
10971108
hm->AddHistogram(histClass, "Mass_Pt_PhiTildePP", "", false, 100, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, 0., 2 * o2::constants::math::PI, VarManager::kPhiTildePP);
10981109
}
10991110
if (subGroupStr.Contains("random")) {
1111+
hm->AddHistogram(histClass, "cosThetaRM", "", false, 100, -1., 1., VarManager::kCosThetaRM);
11001112
hm->AddHistogram(histClass, "Mass_Pt_CosThetaRM", "", false, 200, 1.0, 5.0, VarManager::kMass, 40, 0.0, 20.0, VarManager::kPt, 20, -1., 1., VarManager::kCosThetaRM);
11011113
}
11021114
}

‎PWGDQ/Core/VarManager.h‎

Lines changed: 15 additions & 48 deletions
Original file line numberDiff line numberDiff line change
@@ -156,13 +156,9 @@ class VarManager : public TObject
156156
kDecayToKPi, // e.g. D0 -> K+ pi- or cc.
157157
kTripleCandidateToKPiPi, // e.g. D+ -> K- pi+ pi+
158158
kTripleCandidateToPKPi, // e.g. Lambda_c -> p K- pi+
159-
kNMaxCandidateTypes
160-
};
161-
162-
enum HadronMassCandidateType {
163-
// The mass of the associated hadron
164-
kRealHadronMass = 0, // using the real hadron mass
165-
kTreatAsPion // treat the hadron as pion
159+
kNMaxCandidateTypes,
160+
kJpsiHadronMass, // using the real hadron mass
161+
kJpsiPionMass // treat the hadron as pion
166162
};
167163

168164
enum BarrelTrackFilteringBits {
@@ -1156,8 +1152,8 @@ class VarManager : public TObject
11561152
static void FillTrackCollisionMatCorr(T const& track, C const& collision, M const& materialCorr, P const& propagator, float* values = nullptr);
11571153
template <typename U, typename T>
11581154
static void FillTrackMC(const U& mcStack, T const& track, float* values = nullptr);
1159-
template <int massType, typename T, typename T1>
1160-
static void FillEnergyCorrelators(T const& track, T1 const& t1, float* values = nullptr);
1155+
template <int pairType, typename T, typename T1>
1156+
static void FillEnergyCorrelatorsMC(T const& track, T1 const& t1, float* values = nullptr);
11611157
template <uint32_t fillMap, typename T1, typename T2, typename C>
11621158
static void FillPairPropagateMuon(T1 const& muon1, T2 const& muon2, const C& collision, float* values = nullptr);
11631159
template <uint32_t fillMap, typename T1, typename T2, typename C>
@@ -1176,8 +1172,8 @@ class VarManager : public TObject
11761172
static void FillPairME(T1 const& t1, T2 const& t2, float* values = nullptr);
11771173
template <int pairType, typename T1, typename T2>
11781174
static void FillPairMC(T1 const& t1, T2 const& t2, float* values = nullptr);
1179-
template <int candidateType, typename T1, typename T2, typename T3>
1180-
static void FillTripleMC(T1 const& t1, T2 const& t2, T3 const& t3, float* values = nullptr);
1175+
template <typename T1, typename T2, typename T3>
1176+
static void FillTripleMC(T1 const& t1, T2 const& t2, T3 const& t3, float* values = nullptr, PairCandidateType pairType = kTripleCandidateToEEPhoton);
11811177
template <int candidateType, typename T1, typename T2>
11821178
static void FillQuadMC(T1 const& t1, T2 const& t2, T2 const& t3, float* values = nullptr);
11831179
template <int pairType, uint32_t collFillMap, uint32_t fillMap, typename C, typename T>
@@ -2825,16 +2821,16 @@ void VarManager::FillTrackMC(const U& mcStack, T const& track, float* values)
28252821
FillTrackDerived(values);
28262822
}
28272823

2828-
template <int massType, typename T, typename T1>
2829-
void VarManager::FillEnergyCorrelators(T const& track, T1 const& t1, float* values)
2824+
template <int pairType, typename T, typename T1>
2825+
void VarManager::FillEnergyCorrelatorsMC(T const& track, T1 const& t1, float* values)
28302826
{
28312827
// energy correlators
28322828
float MassHadron;
2833-
if constexpr (massType == kRealHadronMass) {
2829+
if constexpr (pairType == kJpsiHadronMass) {
28342830
MassHadron = TMath::Sqrt(t1.e() * t1.e() - t1.p() * t1.p());
28352831
;
28362832
}
2837-
if constexpr (massType == kTreatAsPion) {
2833+
if constexpr (pairType == kJpsiPionMass) {
28382834
MassHadron = o2::constants::physics::MassPionCharged;
28392835
}
28402836
ROOT::Math::PtEtaPhiMVector v1(track.pt(), track.eta(), track.phi(), o2::constants::physics::MassJPsi);
@@ -3797,14 +3793,14 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
37973793
}
37983794
}
37993795

3800-
template <int candidateType, typename T1, typename T2, typename T3>
3801-
void VarManager::FillTripleMC(T1 const& t1, T2 const& t2, T3 const& t3, float* values)
3796+
template <typename T1, typename T2, typename T3>
3797+
void VarManager::FillTripleMC(T1 const& t1, T2 const& t2, T3 const& t3, float* values, PairCandidateType pairType)
38023798
{
38033799
if (!values) {
38043800
values = fgValues;
38053801
}
38063802

3807-
if constexpr (candidateType == kTripleCandidateToEEPhoton) {
3803+
if (pairType == kTripleCandidateToEEPhoton) {
38083804
float m1 = o2::constants::physics::MassElectron;
38093805
float m2 = o2::constants::physics::MassElectron;
38103806
float m3 = o2::constants::physics::MassPhoton;
@@ -3833,7 +3829,7 @@ void VarManager::FillTripleMC(T1 const& t1, T2 const& t2, T3 const& t3, float* v
38333829
values[kPt2] = t2.pt();
38343830
}
38353831

3836-
if constexpr (candidateType == kTripleCandidateToKPiPi) {
3832+
if (pairType == kTripleCandidateToKPiPi) {
38373833
float m1 = o2::constants::physics::MassKaonCharged;
38383834
float m2 = o2::constants::physics::MassPionCharged;
38393835

@@ -3851,35 +3847,6 @@ void VarManager::FillTripleMC(T1 const& t1, T2 const& t2, T3 const& t3, float* v
38513847
values[kS13] = (v1 + v3).M2();
38523848
values[kS23] = (v2 + v3).M2();
38533849
}
3854-
if constexpr (candidateType == kBtoJpsiEEK) {
3855-
float m1 = o2::constants::physics::MassElectron;
3856-
float m2 = o2::constants::physics::MassElectron;
3857-
float m3 = o2::constants::physics::MassKaonCharged;
3858-
float m4 = o2::constants::physics::MassJPsi;
3859-
ROOT::Math::PtEtaPhiMVector v1(t1.pt(), t1.eta(), t1.phi(), m1);
3860-
ROOT::Math::PtEtaPhiMVector v2(t2.pt(), t2.eta(), t2.phi(), m2);
3861-
ROOT::Math::PtEtaPhiMVector v12 = v1 + v2;
3862-
ROOT::Math::PtEtaPhiMVector v3(t3.pt(), t3.eta(), t3.phi(), m3);
3863-
ROOT::Math::PtEtaPhiMVector v123 = v12 + v3;
3864-
values[kPairMass] = v123.M();
3865-
values[kPairPt] = v123.Pt();
3866-
values[kPairEta] = v123.Eta();
3867-
values[kPhi] = v123.Phi();
3868-
values[kMCY] = -v123.Rapidity();
3869-
values[kPairMassDau] = v12.M();
3870-
values[kPairPtDau] = v12.Pt();
3871-
values[kRap] = -v123.Rapidity();
3872-
values[kMassDau] = m3;
3873-
values[VarManager::kDeltaMass] = v123.M() - v12.M();
3874-
values[VarManager::kDeltaMass_jpsi] = v123.M() - v12.M() + m4;
3875-
values[kPt] = t3.pt();
3876-
values[kEta] = t3.eta();
3877-
values[kEta1] = t1.eta();
3878-
values[kEta2] = t2.eta();
3879-
values[kDeltaEta] = v12.Eta();
3880-
values[kPt1] = t1.pt();
3881-
values[kPt2] = t2.pt();
3882-
}
38833850
}
38843851

38853852
template <int pairType, uint32_t collFillMap, uint32_t fillMap, typename C, typename T>

0 commit comments

Comments
 (0)