Skip to content

Commit 568ec0c

Browse files
prchakraprchakra
authored
PWGCF / FemtoUniverse : Fixing bug in two-particle pair kinematics (#7628)
* Fixing bug in two-particle pair kinematics * Fixing bug in two-particle pair kinematics * Fixing bug in two-particle pair kinematics * Fixing bug in two-particle pair kinematics * Fixing bug in two-particle pair kinematics * Fixing bug in two-particle pair kinematics * Fixing error in clang-format * Fixing error in clang-format
1 parent 977f730 commit 568ec0c

6 files changed

Lines changed: 75 additions & 98 deletions

‎PWGCF/FemtoUniverse/Core/FemtoUniverse3DContainer.h‎

Lines changed: 2 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -74,7 +74,6 @@ class FemtoUniverse3DContainer
7474
/// \param mTAxis axis object for the mT axis
7575
/// \param use3dplots Flag to fill 3D plots
7676
/// \param isiden Identical or non-identical particle pair
77-
/// \param islcms LCMS or PRF
7877
template <typename T>
7978
void init_base(std::string folderName, std::string femtoObs1D, std::string femtoObsKout, std::string femtoObsKside, std::string femtoObsKlong, T femtoObsAxis1D, T femtoObsAxisOut, T femtoObsAxisSide, T femtoObsAxisLong, T multAxis, T kTAxis, T mTAxis, T multAxis3D, T mTAxis3D, bool use3dplots, bool isiden)
8079
{
@@ -214,13 +213,13 @@ class FemtoUniverse3DContainer
214213
/// \param isiden Choosing identical or non-identical pairs
215214
/// \param islcm Choosing LCMS or PRF
216215
template <bool isMC, typename T>
217-
void setPair(T const& part1, T const& part2, const int mult, bool use3dplots, bool isiden, bool islcms)
216+
void setPair(T const& part1, T const& part2, const int mult, bool use3dplots, bool isiden)
218217
{
219218
std::vector<double> f3d;
220219
const float kT = FemtoUniverseMath::getkT(part1, mMassOne, part2, mMassTwo);
221220
const float mT = FemtoUniverseMath::getmT(part1, mMassOne, part2, mMassTwo);
222221

223-
f3d = FemtoUniverseMath::getpairmom3d(part1, mMassOne, part2, mMassTwo, isiden, islcms);
222+
f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo, isiden);
224223

225224
const float femtoObs1D = f3d[0];
226225
const float femtoObsKout = f3d[1];

‎PWGCF/FemtoUniverse/Core/FemtoUniverseMath.h‎

Lines changed: 39 additions & 58 deletions
Original file line numberDiff line numberDiff line change
@@ -149,9 +149,8 @@ class FemtoUniverseMath
149149
/// \param part2 Particle 2
150150
/// \param mass2 Mass of particle 2
151151
/// \param isiden Identical or non-identical particle pair
152-
/// \param islcms LCMS or PRF
153152
template <typename T>
154-
static std::vector<double> getpairmom3d(const T& part1, const float mass1, const T& part2, const float mass2, bool isiden, bool islcms)
153+
static std::vector<double> newpairfunc(const T& part1, const float mass1, const T& part2, const float mass2, bool isiden)
155154
{
156155
const double E1 = sqrt(pow(part1.px(), 2) + pow(part1.py(), 2) + pow(part1.pz(), 2) + pow(mass1, 2));
157156
const double E2 = sqrt(pow(part2.px(), 2) + pow(part2.py(), 2) + pow(part2.pz(), 2) + pow(mass2, 2));
@@ -165,80 +164,62 @@ class FemtoUniverseMath
165164
const double tPx = trackSum.px();
166165
const double tPy = trackSum.py();
167166
const double tPz = trackSum.pz();
168-
const double tPE = trackSum.E();
167+
const double tE = trackSum.E();
169168

170-
const double tPt = trackSum.pt();
171-
const double tMt = trackSum.mt();
172-
const double tPinv = std::sqrt((tMt * tMt) - (tPt * tPt));
169+
double tPt = (tPx * tPx + tPy * tPy);
170+
double tMt = (tE * tE - tPz * tPz);
171+
double tM = sqrt(tMt - tPt);
172+
tMt = sqrt(tMt);
173+
tPt = sqrt(tPt);
173174

174-
float nullmass = 0.0;
175-
const double m1 = std::max(nullmass, mass1);
176-
const double m2 = std::max(nullmass, mass2);
177-
178-
const double tQinvL = std::pow((E1 - E2), 2) - std::pow((part1.px() - part2.px()), 2) -
179-
std::pow((part1.py() - part2.py()), 2) - std::pow((part1.pz() - part2.pz()), 2);
180-
181-
double tQ = (m1 - m2) / tPinv;
182-
tQ = ::sqrt(tQ * tQ - tQinvL);
183-
184-
const double fKStarCalc = tQ / 2.0;
185-
vect.push_back(fKStarCalc);
175+
double fDKOutLCMS, fDKSideLCMS, fDKLongLCMS;
176+
double fDKOut, fDKSide, fDKLong, fDE;
177+
double px1LCMS, py1LCMS, pz1LCMS;
178+
double px2LCMS, py2LCMS, pz2LCMS;
179+
double kstar;
186180

187181
// Boost to LCMS
188182

189-
const double beta = tPz / tPE;
190-
const double gamma = tPE / tMt;
183+
const double beta = tPz / tE;
184+
const double gamma = tE / tMt;
191185

192-
const double px1L = (part1.px() * tPx + part1.py() * tPy) / tPt;
193-
const double py1L = (-part1.px() * tPy + part1.py() * tPx) / tPt;
194-
const double pz1L = gamma * (part1.pz() - beta * E1);
195-
const double pE1L = gamma * (E1 - beta * part1.pz());
186+
fDKOut = (part1.px() * tPx + part1.py() * tPy) / tPt;
187+
fDKSide = (-part1.px() * tPy + part1.py() * tPx) / tPt;
188+
fDKLong = gamma * (part1.pz() - beta * E1);
189+
fDE = gamma * (E1 - beta * part1.pz());
196190

197-
const double px2L = (part2.px() * tPx + part2.py() * tPy) / tPt;
198-
const double py2L = (-part2.px() * tPy + part2.py() * tPx) / tPt;
199-
const double pz2L = gamma * (part2.pz() - beta * E2);
200-
const double pE2L = gamma * (E2 - beta * part2.pz());
191+
px1LCMS = fDKOut;
192+
py1LCMS = fDKSide;
193+
pz1LCMS = fDKLong;
194+
// pE1LCMS = fDE;
201195

202-
double fDKOutLCMS;
203-
double fDKSideLCMS;
204-
double fDKLongLCMS;
196+
px2LCMS = (part2.px() * tPx + part2.py() * tPy) / tPt;
197+
py2LCMS = (part2.py() * tPx - part2.px() * tPy) / tPt;
198+
pz2LCMS = gamma * (part2.pz() - beta * E2);
199+
// pE2LCMS = gamma * (E2 - beta * part2.pz());
205200

206-
double fDKOutPRF;
207-
double fDKSidePRF;
208-
double fDKLongPRF;
209-
210-
if (!isiden) {
211-
fDKOutLCMS = px1L;
212-
fDKSideLCMS = py1L;
213-
fDKLongLCMS = pz1L;
214-
} else {
215-
fDKOutLCMS = px1L - px2L;
216-
fDKSideLCMS = py1L - py2L;
217-
fDKLongLCMS = pz1L - pz2L;
218-
}
201+
fDKOutLCMS = px1LCMS - px2LCMS;
202+
fDKSideLCMS = py1LCMS - py2LCMS;
203+
fDKLongLCMS = pz1LCMS - pz2LCMS;
219204

220205
// Boost to PRF
206+
221207
const double betaOut = tPt / tMt;
222-
const double gammaOut = tMt / tPinv;
208+
const double gammaOut = tMt / tM;
223209

224-
if (!isiden) {
225-
fDKOutPRF = gammaOut * (fDKOutLCMS - betaOut * pE1L);
226-
fDKSidePRF = fDKSideLCMS;
227-
fDKLongPRF = fDKLongLCMS;
228-
} else {
229-
fDKOutPRF = gammaOut * (fDKOutLCMS - betaOut * (pE1L - pE2L));
230-
fDKSidePRF = fDKSideLCMS;
231-
fDKLongPRF = fDKLongLCMS;
232-
}
210+
fDKOut = gammaOut * (fDKOut - betaOut * fDE);
211+
kstar = sqrt(fDKOut * fDKOut + fDKSide * fDKSide + fDKLong * fDKLong);
233212

234-
if (islcms) {
213+
if (isiden) {
214+
vect.push_back(2.0 * (kstar));
235215
vect.push_back(fDKOutLCMS);
236216
vect.push_back(fDKSideLCMS);
237217
vect.push_back(fDKLongLCMS);
238218
} else {
239-
vect.push_back(fDKOutPRF);
240-
vect.push_back(fDKSidePRF);
241-
vect.push_back(fDKLongPRF);
219+
vect.push_back(kstar);
220+
vect.push_back(fDKOut);
221+
vect.push_back(fDKSide);
222+
vect.push_back(fDKLong);
242223
}
243224
return vect;
244225
}

‎PWGCF/FemtoUniverse/Core/FemtoUniversePairSHCentMultKt.h‎

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -184,7 +184,7 @@ class PairSHCentMultKt
184184
/// \param ktval kT value
185185
template <typename T>
186186
void fill_mult_NumDen(T const& part1, T const& part2, uint8_t ChosenEventType,
187-
int maxl, int multval, float ktval)
187+
int maxl, int multval, float ktval, bool isiden)
188188
{
189189
int multbinval;
190190
int absmultval = multval;
@@ -201,7 +201,7 @@ class PairSHCentMultKt
201201
return;
202202
}
203203
// std::cout<<"multbinval "<<multbinval<<std::endl;
204-
fill_kT_NumDen(part1, part2, ChosenEventType, maxl, multbinval, ktval);
204+
fill_kT_NumDen(part1, part2, ChosenEventType, maxl, multbinval, ktval, isiden);
205205
}
206206

207207
/// Templated function to access different kT directory and call AddEventPair
@@ -213,7 +213,7 @@ class PairSHCentMultKt
213213
/// \param ktval kT value
214214
template <typename T>
215215
void fill_kT_NumDen(T const& part1, T const& part2, uint8_t ChosenEventType,
216-
int maxl, int multval, float ktval)
216+
int maxl, int multval, float ktval, bool isiden)
217217
{
218218
int ktbinval = -1;
219219
if ((ktval >= KtBins[0]) && (ktval < KtBins[1])) {
@@ -227,7 +227,7 @@ class PairSHCentMultKt
227227
} else {
228228
return;
229229
}
230-
AddEventPair(part1, part2, ChosenEventType, maxl, multval, ktbinval);
230+
AddEventPair(part1, part2, ChosenEventType, maxl, multval, ktbinval, isiden);
231231
}
232232

233233
/// Set the PDG codes of the two particles involved
@@ -262,14 +262,14 @@ class PairSHCentMultKt
262262
/// \param ktval kT value
263263
template <typename T>
264264
void AddEventPair(T const& part1, T const& part2, uint8_t ChosenEventType,
265-
int /*maxl*/, int multval, int ktval)
265+
int /*maxl*/, int multval, int ktval, bool isiden)
266266
{
267267
int fMultBin = multval;
268268
int fKtBin = ktval;
269269
std::vector<std::complex<double>> fYlmBuffer(fMaxJM);
270270
std::vector<double> f3d;
271-
f3d = FemtoUniverseMath::getpairmom3d(part1, mMassOne, part2, mMassTwo,
272-
true, true);
271+
f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo,
272+
isiden);
273273

274274
const float qout = f3d[1];
275275
const float qside = f3d[2];

‎PWGCF/FemtoUniverse/Core/FemtoUniverseSHContainer.h‎

Lines changed: 3 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -159,20 +159,17 @@ class FemtoUniverseSHContainer
159159
/// \param ChosenEventType same or mixed event
160160
/// \param maxl Maximum valie of L component of the spherical harmonics
161161
template <bool isMC, typename T>
162-
void AddEventPair(T const& part1, T const& part2, uint8_t ChosenEventType, int /*maxl*/)
162+
void AddEventPair(T const& part1, T const& part2, uint8_t ChosenEventType, int /*maxl*/, bool isiden)
163163
{
164-
// int fMaxL = 2;
165-
// int fMaxJM = (2+1)*(2+1);
166164
std::vector<std::complex<double>> fYlmBuffer(fMaxJM);
167165
std::vector<double> f3d;
168-
f3d = FemtoUniverseMath::getpairmom3d(part1, mMassOne, part2, mMassTwo, true, true);
166+
f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo, isiden);
169167

170-
// const float qstar = f3d[0];
168+
const float kv = f3d[0];
171169
const float qout = f3d[1];
172170
const float qside = f3d[2];
173171
const float qlong = f3d[3];
174172

175-
double kv = sqrt(qout * qout + qside * qside + qlong * qlong);
176173
int nqbin = fbinctn->GetXaxis()->FindFixBin(kv) - 1;
177174

178175
FemtoUniverseSpherHarMath Ylm;

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrack3DMultKtExtended.cxx‎

Lines changed: 16 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -476,9 +476,9 @@ struct femtoUniversePairTaskTrackTrack3DMultKtExtended {
476476
float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2);
477477

478478
if (!cfgProcessMultBins) {
479-
sameEventCont.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
479+
sameEventCont.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
480480
} else {
481-
std::vector<double> k3d = FemtoUniverseMath::getpairmom3d(p1, mass1, p2, mass2, ConfIsIden, ConfIsLCMS);
481+
std::vector<double> k3d = FemtoUniverseMath::newpairfunc(p1, mass1, p2, mass2, ConfIsIden);
482482
sameEventMultCont.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
483483
}
484484
}
@@ -515,16 +515,16 @@ struct femtoUniversePairTaskTrackTrack3DMultKtExtended {
515515

516516
if (rand > 0.5) {
517517
if (!cfgProcessMultBins) {
518-
sameEventContPP.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
518+
sameEventContPP.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
519519
} else {
520-
k3d = FemtoUniverseMath::getpairmom3d(p1, mass1, p2, mass2, ConfIsIden, ConfIsLCMS);
520+
k3d = FemtoUniverseMath::newpairfunc(p1, mass1, p2, mass2, ConfIsIden);
521521
sameEventMultContPP.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
522522
}
523523
} else {
524524
if (!cfgProcessMultBins) {
525-
sameEventContPP.setPair<isMC>(p2, p1, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
525+
sameEventContPP.setPair<isMC>(p2, p1, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
526526
} else {
527-
k3d = FemtoUniverseMath::getpairmom3d(p2, mass2, p1, mass1, ConfIsIden, ConfIsLCMS);
527+
k3d = FemtoUniverseMath::newpairfunc(p2, mass2, p1, mass1, ConfIsIden);
528528
sameEventMultContPP.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
529529
}
530530
}
@@ -538,16 +538,16 @@ struct femtoUniversePairTaskTrackTrack3DMultKtExtended {
538538

539539
if (rand > 0.5) {
540540
if (!cfgProcessMultBins) {
541-
sameEventContMM.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
541+
sameEventContMM.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
542542
} else {
543-
k3d = FemtoUniverseMath::getpairmom3d(p1, mass1, p2, mass2, ConfIsIden, ConfIsLCMS);
543+
k3d = FemtoUniverseMath::newpairfunc(p1, mass1, p2, mass2, ConfIsIden);
544544
sameEventMultContMM.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
545545
}
546546
} else {
547547
if (!cfgProcessMultBins) {
548-
sameEventContMM.setPair<isMC>(p2, p1, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
548+
sameEventContMM.setPair<isMC>(p2, p1, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
549549
} else {
550-
k3d = FemtoUniverseMath::getpairmom3d(p2, mass2, p1, mass1, ConfIsIden, ConfIsLCMS);
550+
k3d = FemtoUniverseMath::newpairfunc(p2, mass2, p1, mass1, ConfIsIden);
551551
sameEventMultContMM.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
552552
}
553553
}
@@ -650,19 +650,19 @@ struct femtoUniversePairTaskTrackTrack3DMultKtExtended {
650650
case 1: {
651651
float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2);
652652
if (!cfgProcessMultBins) {
653-
mixedEventCont.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
653+
mixedEventCont.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
654654
} else {
655-
std::vector<double> k3d = FemtoUniverseMath::getpairmom3d(p1, mass1, p2, mass2, ConfIsIden, ConfIsLCMS);
655+
std::vector<double> k3d = FemtoUniverseMath::newpairfunc(p1, mass1, p2, mass2, ConfIsIden);
656656
mixedEventMultCont.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
657657
}
658658
break;
659659
}
660660
case 2: {
661661
float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass1);
662662
if (!cfgProcessMultBins) {
663-
mixedEventContPP.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
663+
mixedEventContPP.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
664664
} else {
665-
std::vector<double> k3d = FemtoUniverseMath::getpairmom3d(p1, mass1, p2, mass2, ConfIsIden, ConfIsLCMS);
665+
std::vector<double> k3d = FemtoUniverseMath::newpairfunc(p1, mass1, p2, mass2, ConfIsIden);
666666
mixedEventMultContPP.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
667667
}
668668
break;
@@ -671,9 +671,9 @@ struct femtoUniversePairTaskTrackTrack3DMultKtExtended {
671671
case 3: {
672672
float kT = FemtoUniverseMath::getkT(p1, mass2, p2, mass2);
673673
if (!cfgProcessMultBins) {
674-
mixedEventContMM.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden, ConfIsLCMS);
674+
mixedEventContMM.setPair<isMC>(p1, p2, multCol, twotracksconfigs.ConfUse3D, ConfIsIden);
675675
} else {
676-
std::vector<double> k3d = FemtoUniverseMath::getpairmom3d(p1, mass1, p2, mass2, ConfIsIden, ConfIsLCMS);
676+
std::vector<double> k3d = FemtoUniverseMath::newpairfunc(p1, mass1, p2, mass2, ConfIsIden);
677677
mixedEventMultContMM.fill_3D<float>(k3d[1], k3d[2], k3d[3], multCol, kT);
678678
}
679679
break;

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx‎

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -475,7 +475,7 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
475475
continue;
476476
}
477477
float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2);
478-
sameEventMultCont.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::same, 2, multCol, kT);
478+
sameEventMultCont.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::same, 2, multCol, kT, ConfIsIden);
479479
}
480480
} else {
481481
for (auto& [p1, p2] : combinations(CombinationsStrictlyUpperIndexPolicy(groupPartsOne, groupPartsOne))) {
@@ -507,19 +507,19 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
507507
case 2: {
508508
rand = randgen->Rndm();
509509
if (rand > 0.5) {
510-
sameEventMultContPP.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::same, 2, multCol, kT);
510+
sameEventMultContPP.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::same, 2, multCol, kT, ConfIsIden);
511511
} else if (rand <= 0.5) {
512-
sameEventMultContPP.fill_mult_NumDen(p2, p1, femtoUniverseSHContainer::EventType::same, 2, multCol, kT);
512+
sameEventMultContPP.fill_mult_NumDen(p2, p1, femtoUniverseSHContainer::EventType::same, 2, multCol, kT, ConfIsIden);
513513
}
514514
break;
515515
}
516516

517517
case 3: {
518518
rand = randgen->Rndm();
519519
if (rand > 0.5) {
520-
sameEventMultContMM.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::same, 2, multCol, kT);
520+
sameEventMultContMM.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::same, 2, multCol, kT, ConfIsIden);
521521
} else if (rand <= 0.5) {
522-
sameEventMultContMM.fill_mult_NumDen(p2, p1, femtoUniverseSHContainer::EventType::same, 2, multCol, kT);
522+
sameEventMultContMM.fill_mult_NumDen(p2, p1, femtoUniverseSHContainer::EventType::same, 2, multCol, kT, ConfIsIden);
523523
}
524524
break;
525525
}
@@ -606,15 +606,15 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
606606
float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2);
607607
switch (ContType) {
608608
case 1: {
609-
mixedEventMultCont.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::mixed, 2, multCol, kT);
609+
mixedEventMultCont.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::mixed, 2, multCol, kT, ConfIsIden);
610610
break;
611611
}
612612
case 2: {
613-
mixedEventMultContPP.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::mixed, 2, multCol, kT);
613+
mixedEventMultContPP.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::mixed, 2, multCol, kT, ConfIsIden);
614614
break;
615615
}
616616
case 3: {
617-
mixedEventMultContMM.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::mixed, 2, multCol, kT);
617+
mixedEventMultContMM.fill_mult_NumDen(p1, p2, femtoUniverseSHContainer::EventType::mixed, 2, multCol, kT, ConfIsIden);
618618
break;
619619
}
620620
default:

0 commit comments

Comments
 (0)