Skip to content

Commit 80ef429

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent 0ff30ed commit 80ef429

4 files changed

Lines changed: 116 additions & 121 deletions

File tree

‎PWGCF/FemtoDream/Core/femtoDreamCollisionSelection.h‎

Lines changed: 79 additions & 84 deletions
Original file line numberDiff line numberDiff line change
@@ -16,16 +16,17 @@
1616
#ifndef PWGCF_FEMTODREAM_CORE_FEMTODREAMCOLLISIONSELECTION_H_
1717
#define PWGCF_FEMTODREAM_CORE_FEMTODREAMCOLLISIONSELECTION_H_
1818

19-
#include <string>
20-
#include <iostream>
21-
#include <vector>
2219
#include "Common/CCDB/TriggerAliases.h"
20+
#include "Common/Core/EventPlaneHelper.h"
2321
#include "Common/DataModel/EventSelection.h"
22+
#include "Common/DataModel/Qvectors.h"
23+
2424
#include "Framework/HistogramRegistry.h"
2525
#include "Framework/Logger.h"
2626

27-
#include "Common/Core/EventPlaneHelper.h"
28-
#include "Common/DataModel/Qvectors.h"
27+
#include <iostream>
28+
#include <string>
29+
#include <vector>
2930

3031
using namespace o2::framework;
3132

@@ -165,24 +166,19 @@ class FemtoDreamCollisionSelection
165166
/// \param col Collision
166167
/// \return whether or not the collisions fulfills the specified selections
167168
template <typename C>
168-
bool isPileUpCollisionPbPb(C const& col,
169-
bool noSameBunchPileup, bool isGoodZvtxFT0vsPV,
170-
bool isGoodITSLayersAll, bool noCollInRofStandard,
171-
bool noHighMultCollInPrevRof, bool noCollInTimeRangeStandard,
172-
bool /*isVertexITSTPC*/,
173-
int tpcOccupancyMin, int tpcOccupancyMax)
169+
bool isPileUpCollisionPbPb(C const& col,
170+
bool noSameBunchPileup, bool isGoodZvtxFT0vsPV,
171+
bool isGoodITSLayersAll, bool noCollInRofStandard,
172+
bool noHighMultCollInPrevRof, bool noCollInTimeRangeStandard,
173+
bool /*isVertexITSTPC*/,
174+
int tpcOccupancyMin, int tpcOccupancyMax)
174175
{
175176
const auto occupancy = col.trackOccupancyInTimeRange();
176177
if ((occupancy < tpcOccupancyMin || occupancy > tpcOccupancyMax)) {
177178
return false;
178179
}
179-
if ((noSameBunchPileup && !col.selection_bit(aod::evsel::kNoSameBunchPileup))
180-
|| (isGoodZvtxFT0vsPV && !col.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))
181-
|| (isGoodITSLayersAll && !col.selection_bit(aod::evsel::kIsGoodITSLayersAll))
182-
|| (noCollInRofStandard && !col.selection_bit(aod::evsel::kNoCollInRofStandard))
183-
|| (noHighMultCollInPrevRof && !col.selection_bit(aod::evsel::kNoHighMultCollInPrevRof))
184-
|| (noCollInTimeRangeStandard && !col.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))
185-
// || (isVertexITSTPC && !col.selection_bit(aod::evsel::kIsVertexITSTPC))
180+
if ((noSameBunchPileup && !col.selection_bit(aod::evsel::kNoSameBunchPileup)) || (isGoodZvtxFT0vsPV && !col.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)) || (isGoodITSLayersAll && !col.selection_bit(aod::evsel::kIsGoodITSLayersAll)) || (noCollInRofStandard && !col.selection_bit(aod::evsel::kNoCollInRofStandard)) || (noHighMultCollInPrevRof && !col.selection_bit(aod::evsel::kNoHighMultCollInPrevRof)) || (noCollInTimeRangeStandard && !col.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))
181+
// || (isVertexITSTPC && !col.selection_bit(aod::evsel::kIsVertexITSTPC))
186182
) {
187183
return false;
188184
}
@@ -221,17 +217,17 @@ class FemtoDreamCollisionSelection
221217
mReQ2thisEvt = new TH2D("ReQ2thisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
222218
mImQ2thisEvt = new TH2D("ImQ2thisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
223219
mMQthisEvt = new TH2D("MQthisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
224-
mMQWeightthisEvt = new TH2D("MQWeightthisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
220+
mMQWeightthisEvt = new TH2D("MQWeightthisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
225221
mHistogramQn = registry;
226222
mHistogramQn->add("Event/centFT0CBefore", "; cent", kTH1F, {{10, 0, 100}});
227223
mHistogramQn->add("Event/centFT0CAfter", "; cent", kTH1F, {{10, 0, 100}});
228224
mHistogramQn->add("Event/centVsqn", "; cent; qn", kTH2F, {{10, 0, 100}, {100, 0, 1000}});
229225
mHistogramQn->add("Event/centVsqnVsSpher", "; cent; qn; Sphericity", kTH3F, {{10, 0, 100}, {100, 0, 1000}, {100, 0, 1}});
230226
mHistogramQn->add("Event/qnBin", "; qnBin; entries", kTH1F, {{20, 0, 20}});
231-
227+
232228
mHistogramQn->add<TProfile>("Event/profileC22", "; cent; c22", kTProfile, {{10, 0, 100}});
233229
mHistogramQn->add<TProfile>("Event/profileC24", "; cent; c24", kTProfile, {{10, 0, 100}});
234-
if (doQnSeparation){
230+
if (doQnSeparation) {
235231
for (int iqn(0); iqn < mumQnBins; ++iqn) {
236232
mHistogramQn->add<TProfile>(("Qn/profileC22_" + std::to_string(iqn)).c_str(), "; cent; c22", kTProfile, {{10, 0, 100}});
237233
mHistogramQn->add(("Qn/mult_" + std::to_string(iqn)).c_str(), "; cent; c22", kTH1F, {{100, 0, 1000}});
@@ -325,14 +321,14 @@ class FemtoDreamCollisionSelection
325321

326322
/// \todo to be implemented!
327323
/// \return the 1-d qn-vector separator to 2-d
328-
std::vector<std::vector<float>> getQnBinSeparator2D(std::vector<float> flat)
324+
std::vector<std::vector<float>> getQnBinSeparator2D(std::vector<float> flat)
329325
{
330326
constexpr size_t nBins = 11;
331327

332328
if (flat.empty() || flat.size() % nBins != 0) {
333-
LOGP(error, "ConfQnBinSeparator size = {} is not divisible by {}",
329+
LOGP(error, "ConfQnBinSeparator size = {} is not divisible by {}",
334330
flat.size(), nBins);
335-
return {{-999, -999}};
331+
return {{-999, -999}};
336332
}
337333

338334
size_t nCent = flat.size() / nBins;
@@ -361,14 +357,14 @@ class FemtoDreamCollisionSelection
361357
return -999; // safe fallback
362358
}
363359

364-
mHistogramQn->fill(HIST("Event/centFT0CBefore"), col.centFT0C());
360+
mHistogramQn->fill(HIST("Event/centFT0CBefore"), col.centFT0C());
365361
int qnBin = -999;
366362
float qn = computeqnVec(col);
367363
int mycentBin = static_cast<int>(col.centFT0C() / centBinWidth);
368364
if (mycentBin >= static_cast<int>(centMax / centBinWidth))
369365
return qnBin;
370366

371-
if (mycentBin > static_cast<int>(twoDSeparator.size()) -1)
367+
if (mycentBin > static_cast<int>(twoDSeparator.size()) - 1)
372368
return qnBin;
373369

374370
for (int iqn(0); iqn < static_cast<int>(twoDSeparator[mycentBin].size()) - 1; ++iqn) {
@@ -379,12 +375,12 @@ class FemtoDreamCollisionSelection
379375
continue;
380376
}
381377
}
382-
378+
383379
mHistogramQn->fill(HIST("Event/centFT0CAfter"), col.centFT0C());
384380
mHistogramQn->fill(HIST("Event/centVsqn"), col.centFT0C(), qn);
385381
mHistogramQn->fill(HIST("Event/centVsqnVsSpher"), col.centFT0C(), qn, fSpher);
386382
mHistogramQn->fill(HIST("Event/qnBin"), qnBin);
387-
if (qnBin >= 0 && qnBin < 10){
383+
if (qnBin >= 0 && qnBin < 10) {
388384
switch (qnBin) {
389385
case 0:
390386
mHistogramQn->fill(HIST("Qn/mult_") + HIST("0"), fMult);
@@ -393,8 +389,8 @@ class FemtoDreamCollisionSelection
393389
mHistogramQn->fill(HIST("Qn/mult_") + HIST("1"), fMult);
394390
break;
395391
case 2:
396-
mHistogramQn->fill(HIST("Qn/mult_") + HIST("2"), fMult);
397-
break;
392+
mHistogramQn->fill(HIST("Qn/mult_") + HIST("2"), fMult);
393+
break;
398394
case 3:
399395
mHistogramQn->fill(HIST("Qn/mult_") + HIST("3"), fMult);
400396
break;
@@ -415,9 +411,9 @@ class FemtoDreamCollisionSelection
415411
break;
416412
case 9:
417413
mHistogramQn->fill(HIST("Qn/mult_") + HIST("9"), fMult);
418-
break;
414+
break;
419415
default:
420-
return qnBin; // invalid qn bin
416+
return qnBin; // invalid qn bin
421417
}
422418
}
423419
return qnBin;
@@ -430,34 +426,34 @@ class FemtoDreamCollisionSelection
430426
/// \tparam T2 type of the tracks
431427
/// \param tracks All tracks
432428
template <typename T1, typename T2>
433-
void fillCumulants(T1 const& col, T2 const& tracks, float fHarmonic=2.f)
434-
{
429+
void fillCumulants(T1 const& col, T2 const& tracks, float fHarmonic = 2.f)
430+
{
435431
int numOfTracks = col.numContrib();
436432
if (numOfTracks < 3)
437-
return ;
433+
return;
438434

439435
mReQthisEvt->Reset();
440436
mImQthisEvt->Reset();
441-
mReQ2thisEvt->Reset();
442-
mImQ2thisEvt->Reset();
437+
mReQ2thisEvt->Reset();
438+
mImQ2thisEvt->Reset();
443439
mMQthisEvt->Reset();
444-
mMQWeightthisEvt->Reset();
440+
mMQWeightthisEvt->Reset();
445441

446442
for (auto const& track : tracks) {
447-
double weight=1; // Will implement NUA&NUE correction
443+
double weight = 1; // Will implement NUA&NUE correction
448444
double phi = track.phi();
449445
double pt = track.pt();
450446
double eta = track.eta();
451-
double cosnphi = weight * TMath::Cos(fHarmonic*phi);
452-
double sinnphi = weight * TMath::Sin(fHarmonic*phi);
453-
double cos2nphi = weight * TMath::Cos(2*fHarmonic*phi);
454-
double sin2nphi = weight * TMath::Sin(2*fHarmonic*phi);
447+
double cosnphi = weight * TMath::Cos(fHarmonic * phi);
448+
double sinnphi = weight * TMath::Sin(fHarmonic * phi);
449+
double cos2nphi = weight * TMath::Cos(2 * fHarmonic * phi);
450+
double sin2nphi = weight * TMath::Sin(2 * fHarmonic * phi);
455451
mReQthisEvt->Fill(pt, eta, cosnphi);
456452
mImQthisEvt->Fill(pt, eta, sinnphi);
457-
mReQ2thisEvt->Fill(pt,eta,cos2nphi);
458-
mImQ2thisEvt->Fill(pt, eta, sin2nphi);
459-
mMQthisEvt ->Fill(pt, eta);
460-
mMQWeightthisEvt ->Fill(pt, eta, weight);
453+
mReQ2thisEvt->Fill(pt, eta, cos2nphi);
454+
mImQ2thisEvt->Fill(pt, eta, sin2nphi);
455+
mMQthisEvt->Fill(pt, eta);
456+
mMQWeightthisEvt->Fill(pt, eta, weight);
461457
}
462458
return;
463459
}
@@ -471,73 +467,72 @@ class FemtoDreamCollisionSelection
471467
template <typename T>
472468
void doCumulants(T const& col, bool doQnSeparation = false, int qnBin = -999, int mumQnBinNum = 10, float fEtaGap = 0.3f, int binPt = 100, int binEta = 32)
473469
{
474-
if (mMQthisEvt->Integral(1, binPt, 1, binEta) < 2)
470+
if (mMQthisEvt->Integral(1, binPt, 1, binEta) < 2)
475471
return;
476-
477-
double allReQ = mReQthisEvt ->Integral(1, binPt, 1, binEta);
478-
double allImQ = mImQthisEvt ->Integral(1, binPt, 1, binEta);
472+
473+
double allReQ = mReQthisEvt->Integral(1, binPt, 1, binEta);
474+
double allImQ = mImQthisEvt->Integral(1, binPt, 1, binEta);
479475
TComplex Q(allReQ, allImQ);
480476
TComplex QStar = TComplex::Conjugate(Q);
481-
482-
double posEtaRe = mReQthisEvt->Integral(1, binPt, mReQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta);
483-
double posEtaIm = mImQthisEvt->Integral(1, binPt, mImQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta);
484-
if (mMQthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta) < 2)
477+
478+
double posEtaRe = mReQthisEvt->Integral(1, binPt, mReQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta);
479+
double posEtaIm = mImQthisEvt->Integral(1, binPt, mImQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta);
480+
if (mMQthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta) < 2)
485481
return;
486-
float posEtaMQ = mMQWeightthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta);
482+
float posEtaMQ = mMQWeightthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta);
487483
TComplex posEtaQ = TComplex(posEtaRe, posEtaIm);
488484
TComplex posEtaQStar = TComplex::Conjugate(posEtaQ);
489-
490-
double negEtaRe = mReQthisEvt->Integral(1, binPt, 1, mReQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6));
491-
double negEtaIm = mImQthisEvt->Integral(1, binPt, 1, mImQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6));
492-
if (mMQthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6)) < 2)
485+
486+
double negEtaRe = mReQthisEvt->Integral(1, binPt, 1, mReQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6));
487+
double negEtaIm = mImQthisEvt->Integral(1, binPt, 1, mImQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6));
488+
if (mMQthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6)) < 2)
493489
return;
494-
float negEtaMQ = mMQWeightthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6));
490+
float negEtaMQ = mMQWeightthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6));
495491
TComplex negEtaQ = TComplex(negEtaRe, negEtaIm);
496492
TComplex negEtaQStar = TComplex::Conjugate(negEtaQ);
497493

498-
mHistogramQn->get<TProfile>(HIST("Event/profileC22"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
499-
if (doQnSeparation && qnBin >= 0 && qnBin < mumQnBinNum){
494+
mHistogramQn->get<TProfile>(HIST("Event/profileC22"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
495+
if (doQnSeparation && qnBin >= 0 && qnBin < mumQnBinNum) {
500496
switch (qnBin) {
501497
case 0:
502-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("0"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
498+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("0"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
503499
break;
504500
case 1:
505-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("1"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
501+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("1"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
506502
break;
507503
case 2:
508-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("2"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
504+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("2"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
509505
break;
510506
case 3:
511-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("3"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
507+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("3"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
512508
break;
513509
case 4:
514-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("4"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
510+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("4"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
515511
break;
516512
case 5:
517-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("5"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
513+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("5"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
518514
break;
519515
case 6:
520-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("6"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
516+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("6"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
521517
break;
522518
case 7:
523-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("7"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
519+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("7"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
524520
break;
525521
case 8:
526-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("8"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
522+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("8"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
527523
break;
528524
case 9:
529-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("9"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
530-
break;
525+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("9"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
526+
break;
531527
default:
532-
return; // invalid qn bin
528+
return; // invalid qn bin
533529
}
534530
}
535531
return;
536532
}
537533

538-
539534
private:
540-
HistogramRegistry* mHistogramRegistry = nullptr; ///< For QA output
535+
HistogramRegistry* mHistogramRegistry = nullptr; ///< For QA output
541536
bool mCutsSet = false; ///< Protection against running without cuts
542537
bool mCheckTrigger = false; ///< Check for trigger
543538
bool mCheckOffline = false; ///< Check for offline criteria (might change)
@@ -547,13 +542,13 @@ class FemtoDreamCollisionSelection
547542
float mZvtxMax = 999.f; ///< Maximal deviation from nominal z-vertex (cm)
548543
float mMinSphericity = 0.f;
549544
float mSphericityPtmin = 0.f;
550-
HistogramRegistry* mHistogramQn = nullptr; ///< For flow cumulant output
551-
TH2D* mReQthisEvt = nullptr; ///< For flow cumulant in an event
552-
TH2D* mImQthisEvt = nullptr; ///< For flow cumulant in an event
553-
TH2D* mReQ2thisEvt = nullptr; ///< For flow cumulant in an event
554-
TH2D* mImQ2thisEvt = nullptr; ///< For flow cumulant in an event
555-
TH2D* mMQthisEvt = nullptr; ///< For flow cumulant in an event
556-
TH2D* mMQWeightthisEvt = nullptr; ///< For flow cumulant in an event
545+
HistogramRegistry* mHistogramQn = nullptr; ///< For flow cumulant output
546+
TH2D* mReQthisEvt = nullptr; ///< For flow cumulant in an event
547+
TH2D* mImQthisEvt = nullptr; ///< For flow cumulant in an event
548+
TH2D* mReQ2thisEvt = nullptr; ///< For flow cumulant in an event
549+
TH2D* mImQ2thisEvt = nullptr; ///< For flow cumulant in an event
550+
TH2D* mMQthisEvt = nullptr; ///< For flow cumulant in an event
551+
TH2D* mMQWeightthisEvt = nullptr; ///< For flow cumulant in an event
557552
};
558553
} // namespace o2::analysis::femtoDream
559554

0 commit comments

Comments
 (0)