Skip to content

Commit a3baf29

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent 46ba359 commit a3baf29

5 files changed

Lines changed: 99 additions & 104 deletions

File tree

‎PWGCF/DataModel/FemtoDerived.h‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -51,7 +51,7 @@ DECLARE_SOA_COLUMN(BitMaskTrackThree, bitmaskTrackThree, BitMaskType); //! Bit f
5151

5252
DECLARE_SOA_COLUMN(Downsample, downsample, bool); //! Flag for downsampling
5353

54-
DECLARE_SOA_COLUMN(QnBin, qnBin, int); //! Bin of qn-vector of the event
54+
DECLARE_SOA_COLUMN(QnBin, qnBin, int); //! Bin of qn-vector of the event
5555

5656
} // namespace femtodreamcollision
5757

‎PWGCF/FemtoDream/Core/femtoDreamCollisionSelection.h‎

Lines changed: 72 additions & 77 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

@@ -181,13 +182,8 @@ class FemtoDreamCollisionSelection
181182
bool noHighMultCollInPrevRof = true;
182183
bool noCollInTimeRangeStandard = true;
183184
// bool isVertexITSTPC = true;
184-
if ((noSameBunchPileup && !col.selection_bit(aod::evsel::kNoSameBunchPileup))
185-
|| (isGoodZvtxFT0vsPV && !col.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))
186-
|| (isGoodITSLayersAll && !col.selection_bit(aod::evsel::kIsGoodITSLayersAll))
187-
|| (noCollInRofStandard && !col.selection_bit(aod::evsel::kNoCollInRofStandard))
188-
|| (noHighMultCollInPrevRof && !col.selection_bit(aod::evsel::kNoHighMultCollInPrevRof))
189-
|| (noCollInTimeRangeStandard && !col.selection_bit(aod::evsel::kNoCollInTimeRangeStandard))
190-
// || (isVertexITSTPC && !col.selection_bit(aod::evsel::kIsVertexITSTPC))
185+
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))
186+
// || (isVertexITSTPC && !col.selection_bit(aod::evsel::kIsVertexITSTPC))
191187
) {
192188
return false;
193189
}
@@ -226,17 +222,17 @@ class FemtoDreamCollisionSelection
226222
mReQ2thisEvt = new TH2D("ReQ2thisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
227223
mImQ2thisEvt = new TH2D("ImQ2thisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
228224
mMQthisEvt = new TH2D("MQthisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
229-
mMQWeightthisEvt = new TH2D("MQWeightthisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
225+
mMQWeightthisEvt = new TH2D("MQWeightthisEvt", "", binPt, 0., 5., binEta, -0.8, 0.8);
230226
mHistogramQn = registry;
231227
mHistogramQn->add("Event/centFT0CBefore", "; cent", kTH1F, {{10, 0, 100}});
232228
mHistogramQn->add("Event/centFT0CAfter", "; cent", kTH1F, {{10, 0, 100}});
233229
mHistogramQn->add("Event/centVsqn", "; cent; qn", kTH2F, {{10, 0, 100}, {100, 0, 1000}});
234230
mHistogramQn->add("Event/centVsqnVsSpher", "; cent; qn; Sphericity", kTH3F, {{10, 0, 100}, {100, 0, 1000}, {100, 0, 1}});
235231
mHistogramQn->add("Event/qnBin", "; qnBin; entries", kTH1F, {{20, 0, 20}});
236-
232+
237233
mHistogramQn->add<TProfile>("Event/profileC22", "; cent; c22", kTProfile, {{10, 0, 100}});
238234
mHistogramQn->add<TProfile>("Event/profileC24", "; cent; c24", kTProfile, {{10, 0, 100}});
239-
if (doQnSeparation){
235+
if (doQnSeparation) {
240236
for (int iqn(0); iqn < mumQnBins; ++iqn) {
241237
mHistogramQn->add<TProfile>(("Qn/profileC22_" + std::to_string(iqn)).c_str(), "; cent; c22", kTProfile, {{10, 0, 100}});
242238
mHistogramQn->add(("Qn/mult_" + std::to_string(iqn)).c_str(), "; cent; c22", kTH1F, {{100, 0, 1000}});
@@ -331,14 +327,14 @@ class FemtoDreamCollisionSelection
331327
/// \todo to be implemented!
332328
/// \param col Collision
333329
/// \return the 1-d qn-vector separator to 2-d
334-
std::vector<std::vector<float>> getQnBinSeparator2D(std::vector<float> flat)
330+
std::vector<std::vector<float>> getQnBinSeparator2D(std::vector<float> flat)
335331
{
336332
constexpr size_t nBins = 11;
337333

338334
if (flat.empty() || flat.size() % nBins != 0) {
339-
LOGP(error, "ConfQnBinSeparator size = {} is not divisible by {}",
335+
LOGP(error, "ConfQnBinSeparator size = {} is not divisible by {}",
340336
flat.size(), nBins);
341-
return {{-999, -999}};
337+
return {{-999, -999}};
342338
}
343339

344340
size_t nCent = flat.size() / nBins;
@@ -367,14 +363,14 @@ class FemtoDreamCollisionSelection
367363
return -999; // safe fallback
368364
}
369365

370-
mHistogramQn->fill(HIST("Event/centFT0CBefore"), col.centFT0C());
366+
mHistogramQn->fill(HIST("Event/centFT0CBefore"), col.centFT0C());
371367
int qnBin = -999;
372368
float qn = computeqnVec(col);
373369
int mycentBin = static_cast<int>(col.centFT0C() / centBinWidth);
374370
if (mycentBin >= static_cast<int>(centMax / centBinWidth))
375371
return qnBin;
376372

377-
if (mycentBin > static_cast<int>(twoDSeparator.size()) -1)
373+
if (mycentBin > static_cast<int>(twoDSeparator.size()) - 1)
378374
return qnBin;
379375

380376
for (int iqn(0); iqn < static_cast<int>(twoDSeparator[mycentBin].size()) - 1; ++iqn) {
@@ -385,12 +381,12 @@ class FemtoDreamCollisionSelection
385381
continue;
386382
}
387383
}
388-
384+
389385
mHistogramQn->fill(HIST("Event/centFT0CAfter"), col.centFT0C());
390386
mHistogramQn->fill(HIST("Event/centVsqn"), col.centFT0C(), qn);
391387
mHistogramQn->fill(HIST("Event/centVsqnVsSpher"), col.centFT0C(), qn, fSpher);
392388
mHistogramQn->fill(HIST("Event/qnBin"), qnBin);
393-
if (qnBin >= 0 && qnBin < 10){
389+
if (qnBin >= 0 && qnBin < 10) {
394390
switch (qnBin) {
395391
case 0:
396392
mHistogramQn->fill(HIST("Qn/mult_") + HIST("0"), fMult);
@@ -399,8 +395,8 @@ class FemtoDreamCollisionSelection
399395
mHistogramQn->fill(HIST("Qn/mult_") + HIST("1"), fMult);
400396
break;
401397
case 2:
402-
mHistogramQn->fill(HIST("Qn/mult_") + HIST("2"), fMult);
403-
break;
398+
mHistogramQn->fill(HIST("Qn/mult_") + HIST("2"), fMult);
399+
break;
404400
case 3:
405401
mHistogramQn->fill(HIST("Qn/mult_") + HIST("3"), fMult);
406402
break;
@@ -421,9 +417,9 @@ class FemtoDreamCollisionSelection
421417
break;
422418
case 9:
423419
mHistogramQn->fill(HIST("Qn/mult_") + HIST("9"), fMult);
424-
break;
420+
break;
425421
default:
426-
return qnBin; // invalid qn bin
422+
return qnBin; // invalid qn bin
427423
}
428424
}
429425
return qnBin;
@@ -439,34 +435,34 @@ class FemtoDreamCollisionSelection
439435
/// \param tracks All tracks
440436
/// \return value of the sphericity of the event
441437
template <typename T1, typename T2>
442-
void fillCumulants(T1 const& col, T2 const& tracks, float fHarmonic=2.f)
443-
{
438+
void fillCumulants(T1 const& col, T2 const& tracks, float fHarmonic = 2.f)
439+
{
444440
int numOfTracks = col.numContrib();
445441
if (numOfTracks < 3)
446-
return ;
442+
return;
447443

448444
mReQthisEvt->Reset();
449445
mImQthisEvt->Reset();
450-
mReQ2thisEvt->Reset();
451-
mImQ2thisEvt->Reset();
446+
mReQ2thisEvt->Reset();
447+
mImQ2thisEvt->Reset();
452448
mMQthisEvt->Reset();
453-
mMQWeightthisEvt->Reset();
449+
mMQWeightthisEvt->Reset();
454450

455451
for (auto const& track : tracks) {
456-
double weight=1; // Will implement NUA&NUE correction
452+
double weight = 1; // Will implement NUA&NUE correction
457453
double phi = track.phi();
458454
double pt = track.pt();
459455
double eta = track.eta();
460-
double cosnphi = weight * TMath::Cos(fHarmonic*phi);
461-
double sinnphi = weight * TMath::Sin(fHarmonic*phi);
462-
double cos2nphi = weight * TMath::Cos(2*fHarmonic*phi);
463-
double sin2nphi = weight * TMath::Sin(2*fHarmonic*phi);
456+
double cosnphi = weight * TMath::Cos(fHarmonic * phi);
457+
double sinnphi = weight * TMath::Sin(fHarmonic * phi);
458+
double cos2nphi = weight * TMath::Cos(2 * fHarmonic * phi);
459+
double sin2nphi = weight * TMath::Sin(2 * fHarmonic * phi);
464460
mReQthisEvt->Fill(pt, eta, cosnphi);
465461
mImQthisEvt->Fill(pt, eta, sinnphi);
466-
mReQ2thisEvt->Fill(pt,eta,cos2nphi);
467-
mImQ2thisEvt->Fill(pt, eta, sin2nphi);
468-
mMQthisEvt ->Fill(pt, eta);
469-
mMQWeightthisEvt ->Fill(pt, eta, weight);
462+
mReQ2thisEvt->Fill(pt, eta, cos2nphi);
463+
mImQ2thisEvt->Fill(pt, eta, sin2nphi);
464+
mMQthisEvt->Fill(pt, eta);
465+
mMQWeightthisEvt->Fill(pt, eta, weight);
470466
}
471467

472468
return;
@@ -484,73 +480,72 @@ class FemtoDreamCollisionSelection
484480
template <typename T>
485481
void doCumulants(T const& col, bool doQnSeparation = false, int qnBin = -999, int mumQnBinNum = 10, float fEtaGap = 0.3f, int binPt = 100, int binEta = 32)
486482
{
487-
if (mMQthisEvt->Integral(1, binPt, 1, binEta) < 2)
483+
if (mMQthisEvt->Integral(1, binPt, 1, binEta) < 2)
488484
return;
489-
490-
double allReQ = mReQthisEvt ->Integral(1, binPt, 1, binEta);
491-
double allImQ = mImQthisEvt ->Integral(1, binPt, 1, binEta);
485+
486+
double allReQ = mReQthisEvt->Integral(1, binPt, 1, binEta);
487+
double allImQ = mImQthisEvt->Integral(1, binPt, 1, binEta);
492488
TComplex Q(allReQ, allImQ);
493489
TComplex QStar = TComplex::Conjugate(Q);
494-
495-
double posEtaRe = mReQthisEvt->Integral(1, binPt, mReQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta);
496-
double posEtaIm = mImQthisEvt->Integral(1, binPt, mImQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta);
497-
if (mMQthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta) < 2)
490+
491+
double posEtaRe = mReQthisEvt->Integral(1, binPt, mReQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta);
492+
double posEtaIm = mImQthisEvt->Integral(1, binPt, mImQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta);
493+
if (mMQthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta) < 2)
498494
return;
499-
float posEtaMQ = mMQWeightthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap+1e-6), binEta);
495+
float posEtaMQ = mMQWeightthisEvt->Integral(1, binPt, mMQthisEvt->GetYaxis()->FindBin(fEtaGap + 1e-6), binEta);
500496
TComplex posEtaQ = TComplex(posEtaRe, posEtaIm);
501497
TComplex posEtaQStar = TComplex::Conjugate(posEtaQ);
502-
503-
double negEtaRe = mReQthisEvt->Integral(1, binPt, 1, mReQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6));
504-
double negEtaIm = mImQthisEvt->Integral(1, binPt, 1, mImQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6));
505-
if (mMQthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6)) < 2)
498+
499+
double negEtaRe = mReQthisEvt->Integral(1, binPt, 1, mReQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6));
500+
double negEtaIm = mImQthisEvt->Integral(1, binPt, 1, mImQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6));
501+
if (mMQthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6)) < 2)
506502
return;
507-
float negEtaMQ = mMQWeightthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1*fEtaGap-1e-6));
503+
float negEtaMQ = mMQWeightthisEvt->Integral(1, binPt, 1, mMQthisEvt->GetYaxis()->FindBin(-1 * fEtaGap - 1e-6));
508504
TComplex negEtaQ = TComplex(negEtaRe, negEtaIm);
509505
TComplex negEtaQStar = TComplex::Conjugate(negEtaQ);
510506

511-
mHistogramQn->get<TProfile>(HIST("Event/profileC22"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
512-
if (doQnSeparation && qnBin >= 0 && qnBin < mumQnBinNum){
507+
mHistogramQn->get<TProfile>(HIST("Event/profileC22"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
508+
if (doQnSeparation && qnBin >= 0 && qnBin < mumQnBinNum) {
513509
switch (qnBin) {
514510
case 0:
515-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("0"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
511+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("0"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
516512
break;
517513
case 1:
518-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("1"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
514+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("1"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
519515
break;
520516
case 2:
521-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("2"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
517+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("2"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
522518
break;
523519
case 3:
524-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("3"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
520+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("3"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
525521
break;
526522
case 4:
527-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("4"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
523+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("4"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
528524
break;
529525
case 5:
530-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("5"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
526+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("5"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
531527
break;
532528
case 6:
533-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("6"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
529+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("6"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
534530
break;
535531
case 7:
536-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("7"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
532+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("7"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
537533
break;
538534
case 8:
539-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("8"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
535+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("8"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
540536
break;
541537
case 9:
542-
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("9"))->Fill(col.centFT0C(), (negEtaQ*posEtaQStar).Re()/(negEtaMQ*posEtaMQ), (negEtaMQ*posEtaMQ));
543-
break;
538+
mHistogramQn->get<TProfile>(HIST("Qn/profileC22_") + HIST("9"))->Fill(col.centFT0C(), (negEtaQ * posEtaQStar).Re() / (negEtaMQ * posEtaMQ), (negEtaMQ * posEtaMQ));
539+
break;
544540
default:
545-
return; // invalid qn bin
541+
return; // invalid qn bin
546542
}
547543
}
548544
return;
549545
}
550546

551-
552547
private:
553-
HistogramRegistry* mHistogramRegistry = nullptr; ///< For QA output
548+
HistogramRegistry* mHistogramRegistry = nullptr; ///< For QA output
554549
bool mCutsSet = false; ///< Protection against running without cuts
555550
bool mCheckTrigger = false; ///< Check for trigger
556551
bool mCheckOffline = false; ///< Check for offline criteria (might change)
@@ -560,13 +555,13 @@ class FemtoDreamCollisionSelection
560555
float mZvtxMax = 999.f; ///< Maximal deviation from nominal z-vertex (cm)
561556
float mMinSphericity = 0.f;
562557
float mSphericityPtmin = 0.f;
563-
HistogramRegistry* mHistogramQn = nullptr; ///< For flow cumulant output
564-
TH2D* mReQthisEvt = nullptr; ///< For flow cumulant in an event
565-
TH2D* mImQthisEvt = nullptr; ///< For flow cumulant in an event
566-
TH2D* mReQ2thisEvt = nullptr; ///< For flow cumulant in an event
567-
TH2D* mImQ2thisEvt = nullptr; ///< For flow cumulant in an event
558+
HistogramRegistry* mHistogramQn = nullptr; ///< For flow cumulant output
559+
TH2D* mReQthisEvt = nullptr; ///< For flow cumulant in an event
560+
TH2D* mImQthisEvt = nullptr; ///< For flow cumulant in an event
561+
TH2D* mReQ2thisEvt = nullptr; ///< For flow cumulant in an event
562+
TH2D* mImQ2thisEvt = nullptr; ///< For flow cumulant in an event
568563
TH2D* mMQthisEvt = nullptr;
569-
TH2D* mMQWeightthisEvt = nullptr; ///< For flow cumulant in an event
564+
TH2D* mMQWeightthisEvt = nullptr; ///< For flow cumulant in an event
570565
};
571566
} // namespace o2::analysis::femtoDream
572567

‎PWGCF/FemtoDream/Core/femtoDreamContainer.h‎

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -183,15 +183,15 @@ class FemtoDreamContainer
183183
/// Initialize the histograms for pairs in divided qn bins
184184
template <typename T>
185185
void init_base_qn(std::string folderName, std::string femtoObs,
186-
T& femtoObsAxis, T& mTAxi4D, T& multPercentileAxis4D, T& qnAxis4D)
186+
T& femtoObsAxis, T& mTAxi4D, T& multPercentileAxis4D, T& qnAxis4D)
187187
{
188188
mHistogramRegistry->add((folderName + "/relPairkstarmTMultMultPercentileQn").c_str(), ("; " + femtoObs + "; #it{m}_{T} (GeV/#it{c}); Centrality; qn").c_str(), kTHnSparseF, {femtoObsAxis, mTAxi4D, multPercentileAxis4D, qnAxis4D});
189189
}
190190

191191
template <typename T>
192192
void init_qn(HistogramRegistry* registry,
193-
T& kstarBins4D, T& mTBins4D, T& multPercentileBins4D,
194-
bool isMC, float highkstarCut, ConfigurableAxis qnBins4D = {"qnBins4D", {10,0,10}, "qn binning"})
193+
T& kstarBins4D, T& mTBins4D, T& multPercentileBins4D,
194+
bool isMC, float highkstarCut, ConfigurableAxis qnBins4D = {"qnBins4D", {10, 0, 10}, "qn binning"})
195195
{
196196
mHistogramRegistry = registry;
197197
std::string femtoObs;
@@ -207,13 +207,13 @@ class FemtoDreamContainer
207207

208208
std::string folderName = static_cast<std::string>(mFolderSuffix[mEventType]) + static_cast<std::string>(o2::aod::femtodreamMCparticle::MCTypeName[o2::aod::femtodreamMCparticle::MCType::kRecon]) + static_cast<std::string>("_qn");
209209

210-
init_base_qn(folderName, femtoObs,
211-
kstarAxis4D, mTAxis4D, multPercentileAxis4D, qnAxis4D);
210+
init_base_qn(folderName, femtoObs,
211+
kstarAxis4D, mTAxis4D, multPercentileAxis4D, qnAxis4D);
212212

213213
if (isMC) {
214214
folderName = static_cast<std::string>(mFolderSuffix[mEventType]) + static_cast<std::string>(o2::aod::femtodreamMCparticle::MCTypeName[o2::aod::femtodreamMCparticle::MCType::kTruth]) + static_cast<std::string>("_qn");
215215
init_base_qn(folderName, femtoObs,
216-
kstarAxis4D, mTAxis4D, multPercentileAxis4D, qnAxis4D);
216+
kstarAxis4D, mTAxis4D, multPercentileAxis4D, qnAxis4D);
217217
}
218218
}
219219

@@ -359,7 +359,7 @@ class FemtoDreamContainer
359359
void setPair_qn_base(const float femtoObs, const float mT, const float multPercentile, const int myQnBin, const int numQnBins = 10)
360360
{
361361
if (myQnBin >= 0 && myQnBin < numQnBins) {
362-
mHistogramRegistry->fill(HIST(mFolderSuffix[mEventType]) + HIST(o2::aod::femtodreamMCparticle::MCTypeName[mc]) + HIST("_qn") + HIST("/relPairkstarmTMultMultPercentileQn"), femtoObs, mT, multPercentile, myQnBin);
362+
mHistogramRegistry->fill(HIST(mFolderSuffix[mEventType]) + HIST(o2::aod::femtodreamMCparticle::MCTypeName[mc]) + HIST("_qn") + HIST("/relPairkstarmTMultMultPercentileQn"), femtoObs, mT, multPercentile, myQnBin);
363363
} else {
364364
return;
365365
}

0 commit comments

Comments
 (0)