Skip to content

Commit 1cd69c5

Browse files
fuchuncuifuchuncui
authored
Add files via upload
1 parent 267de0c commit 1cd69c5

File tree

1 file changed

+86
-84
lines changed

1 file changed

+86
-84
lines changed

‎PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx‎

Lines changed: 86 additions & 84 deletions
Original file line numberDiff line numberDiff line change
@@ -9,8 +9,10 @@
99
// granted to it by virtue of its status as an Intergovernmental Organization
1010
// or submit itself to any jurisdiction.
1111

12-
/// In case of questions please write to:
12+
/// \file flowGFWOmegaXi.cxx
1313
/// \author Fuchun Cui([email protected])
14+
/// \since Sep/13/2024
15+
/// \brief This task is to caculate V0s and cascades flow by GenericFramework
1416

1517
#include <CCDB/BasicCCDBManager.h>
1618
#include <vector>
@@ -50,16 +52,16 @@ using namespace o2::framework;
5052
using namespace o2::framework::expressions;
5153
namespace
5254
{
53-
std::shared_ptr<TProfile> REFc22[10];
54-
std::shared_ptr<TProfile> REFc24[10];
55-
std::shared_ptr<TProfile3D> K0sc22[10];
56-
std::shared_ptr<TProfile3D> K0sc24[10];
57-
std::shared_ptr<TProfile3D> Lambdac22[10];
58-
std::shared_ptr<TProfile3D> Lambdac24[10];
59-
std::shared_ptr<TProfile3D> Xic22[10];
60-
std::shared_ptr<TProfile3D> Xic24[10];
61-
std::shared_ptr<TProfile3D> Omegac22[10];
62-
std::shared_ptr<TProfile3D> Omegac24[10];
55+
std::shared_ptr<TProfile> refc22[10];
56+
std::shared_ptr<TProfile> refc24[10];
57+
std::shared_ptr<TProfile3D> k0sc22[10];
58+
std::shared_ptr<TProfile3D> k0sc24[10];
59+
std::shared_ptr<TProfile3D> lambdac22[10];
60+
std::shared_ptr<TProfile3D> lambdac24[10];
61+
std::shared_ptr<TProfile3D> xic22[10];
62+
std::shared_ptr<TProfile3D> xic24[10];
63+
std::shared_ptr<TProfile3D> omegac22[10];
64+
std::shared_ptr<TProfile3D> omegac24[10];
6365
} // namespace
6466

6567
#define O2_DEFINE_CONFIGURABLE(NAME, TYPE, DEFAULT, HELP) Configurable<TYPE> NAME{#NAME, DEFAULT, HELP};
@@ -114,17 +116,17 @@ struct FlowGFWOmegaXi {
114116
O2_DEFINE_CONFIGURABLE(cfgAcceptancePath, std::vector<std::string>, (std::vector<std::string>{"Users/f/fcui/NUA/NUAREFPartical", "Users/f/fcui/NUA/NUAK0s", "Users/f/fcui/NUA/NUALambda", "Users/f/fcui/NUA/NUAXi", "Users/f/fcui/NUA/NUAOmega"}), "CCDB path to acceptance object")
115117
O2_DEFINE_CONFIGURABLE(cfgEfficiencyPath, std::vector<std::string>, (std::vector<std::string>{"PathtoRef"}), "CCDB path to efficiency object")
116118

117-
ConfigurableAxis cfgaxisVertex{"axisVertex", {20, -10, 10}, "vertex axis for histograms"};
118-
ConfigurableAxis cfgaxisPhi{"axisPhi", {60, 0.0, constants::math::TwoPI}, "phi axis for histograms"};
119-
ConfigurableAxis cfgaxisEta{"axisEta", {40, -1., 1.}, "eta axis for histograms"};
120-
ConfigurableAxis cfgaxisPt{"axisPtREF", {VARIABLE_WIDTH, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.20, 2.40, 2.60, 2.80, 3.00, 3.50, 4.00, 4.50, 5.00, 5.50, 6.00, 10.0}, "pt (GeV)"};
121-
ConfigurableAxis cfgaxisPtXi{"axisPtXi", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
122-
ConfigurableAxis cfgaxisPtOmega{"axisPtOmega", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
123-
ConfigurableAxis cfgaxisPtV0{"axisPtV0", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
124-
ConfigurableAxis cfgaxisOmegaminusMassforflow{"axismassOmegaFlow", {16, 1.63f, 1.71f}, "Inv. Mass (GeV)"};
125-
ConfigurableAxis cfgaxisXiminusMassforflow{"axismassXiFlow", {14, 1.3f, 1.37f}, "Inv. Mass (GeV)"};
126-
ConfigurableAxis cfgaxisK0sMassforflow{"axismassK0sFlow", {40, 0.4f, 0.6f}, "Inv. Mass (GeV)"};
127-
ConfigurableAxis cfgaxisLambdaMassforflow{"axismassLambdaFlow", {32, 1.08f, 1.16f}, "Inv. Mass (GeV)"};
119+
ConfigurableAxis cfgaxisVertex{"cfgaxisVertex", {20, -10, 10}, "vertex axis for histograms"};
120+
ConfigurableAxis cfgaxisPhi{"cfgaxisPhi", {60, 0.0, constants::math::TwoPI}, "phi axis for histograms"};
121+
ConfigurableAxis cfgaxisEta{"cfgaxisEta", {40, -1., 1.}, "eta axis for histograms"};
122+
ConfigurableAxis cfgaxisPt{"cfgaxisPt", {VARIABLE_WIDTH, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.20, 2.40, 2.60, 2.80, 3.00, 3.50, 4.00, 4.50, 5.00, 5.50, 6.00, 10.0}, "pt (GeV)"};
123+
ConfigurableAxis cfgaxisPtXi{"cfgaxisPtXi", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
124+
ConfigurableAxis cfgaxisPtOmega{"cfgaxisPtOmega", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
125+
ConfigurableAxis cfgaxisPtV0{"cfgaxisPtV0", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
126+
ConfigurableAxis cfgaxisOmegaminusMassforflow{"cfgaxisOmegaminusMassforflow", {16, 1.63f, 1.71f}, "Inv. Mass (GeV)"};
127+
ConfigurableAxis cfgaxisXiminusMassforflow{"cfgaxisXiminusMassforflow", {14, 1.3f, 1.37f}, "Inv. Mass (GeV)"};
128+
ConfigurableAxis cfgaxisK0sMassforflow{"cfgaxisK0sMassforflow", {40, 0.4f, 0.6f}, "Inv. Mass (GeV)"};
129+
ConfigurableAxis cfgaxisLambdaMassforflow{"cfgaxisLambdaMassforflow", {32, 1.08f, 1.16f}, "Inv. Mass (GeV)"};
128130

129131
AxisSpec axisMultiplicity{{0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90}, "Centrality (%)"};
130132
AxisSpec axisOmegaminusMass = {80, 1.63f, 1.71f, "Inv. Mass (GeV)"};
@@ -165,10 +167,10 @@ struct FlowGFWOmegaXi {
165167
TF1* fT0AV0AMean = nullptr;
166168
TF1* fT0AV0ASigma = nullptr;
167169

168-
using TracksPID = soa::Join<aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr, aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr>;
169-
using aodTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, TracksPID>>; // tracks filter
170+
using tracksPID = soa::Join<aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr, aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr>;
171+
using aodTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, tracksPID>>; // tracks filter
170172
using aodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::Mults>>; // collisions filter
171-
using DaughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, TracksPID>;
173+
using daughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, tracksPID>;
172174

173175
// Declare the pt, mult and phi Axis;
174176
int nPtBins = 0;
@@ -255,16 +257,16 @@ struct FlowGFWOmegaXi {
255257
registry.add("Lambdac24dpt", ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisLambdaMassforflow, axisMultiplicity}});
256258
// for Jackknife
257259
for (int i = 1; i <= 10; i++) {
258-
REFc22[i - 1] = registry.add<TProfile>(Form("Jackknife/REF/c22_%d", i), ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
259-
REFc24[i - 1] = registry.add<TProfile>(Form("Jackknife/REF/c24_%d", i), ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
260-
Xic22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Xi/Xic22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisXiminusMassforflow, axisMultiplicity}});
261-
Omegac22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Omega/Omegac22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisOmegaminusMassforflow, axisMultiplicity}});
262-
K0sc22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/K0s/K0sc22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisK0sMassforflow, axisMultiplicity}});
263-
Lambdac22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Lambda/Lambdac22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisLambdaMassforflow, axisMultiplicity}});
264-
Xic24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Xi/Xic24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisXiminusMassforflow, axisMultiplicity}});
265-
Omegac24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Omega/Omegac24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisOmegaminusMassforflow, axisMultiplicity}});
266-
K0sc24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/K0s/K0sc24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisK0sMassforflow, axisMultiplicity}});
267-
Lambdac24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Lambda/Lambdac24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisLambdaMassforflow, axisMultiplicity}});
260+
refc22[i - 1] = registry.add<TProfile>(Form("Jackknife/REF/c22_%d", i), ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
261+
refc24[i - 1] = registry.add<TProfile>(Form("Jackknife/REF/c24_%d", i), ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
262+
xic22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Xi/Xic22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisXiminusMassforflow, axisMultiplicity}});
263+
omegac22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Omega/Omegac22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisOmegaminusMassforflow, axisMultiplicity}});
264+
k0sc22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/K0s/K0sc22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisK0sMassforflow, axisMultiplicity}});
265+
lambdac22[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Lambda/Lambdac22dpt_%d", i), ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisLambdaMassforflow, axisMultiplicity}});
266+
xic24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Xi/Xic24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisXiminusMassforflow, axisMultiplicity}});
267+
omegac24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Omega/Omegac24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtXi, cfgaxisOmegaminusMassforflow, axisMultiplicity}});
268+
k0sc24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/K0s/K0sc24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisK0sMassforflow, axisMultiplicity}});
269+
lambdac24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Lambda/Lambdac24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisLambdaMassforflow, axisMultiplicity}});
268270
}
269271
// InvMass(GeV) of casc and v0
270272
registry.add("InvMassXi_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
@@ -389,7 +391,7 @@ struct FlowGFWOmegaXi {
389391

390392
// input HIST("name")
391393
template <char... chars>
392-
void FillProfile(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const double& cent)
394+
void fillProfile(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const double& cent)
393395
{
394396
double dnx, val;
395397
dnx = fGFW->Calculate(corrconf, 0, kTRUE).real();
@@ -405,7 +407,7 @@ struct FlowGFWOmegaXi {
405407
}
406408

407409
// input shared_ptr<TProfile>
408-
void FillProfile(const GFW::CorrConfig& corrconf, std::shared_ptr<TProfile> TProfile, const double& cent)
410+
void fillProfile(const GFW::CorrConfig& corrconf, std::shared_ptr<TProfile> TProfile, const double& cent)
409411
{
410412
double dnx, val;
411413
dnx = fGFW->Calculate(corrconf, 0, kTRUE).real();
@@ -421,7 +423,7 @@ struct FlowGFWOmegaXi {
421423
}
422424

423425
template <char... chars>
424-
void FillProfilepT(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const int& ptbin, const double& cent)
426+
void fillProfilepT(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const int& ptbin, const double& cent)
425427
{
426428
float dnx = 0;
427429
float val = 0;
@@ -437,7 +439,7 @@ struct FlowGFWOmegaXi {
437439

438440
// input HIST("name")
439441
template <char... chars>
440-
void FillProfilepTMass(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const int& ptbin, const int& PDGCode, const float& cent)
442+
void fillProfilepTMass(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const int& ptbin, const int& PDGCode, const float& cent)
441443
{
442444
int nMassBins = 0;
443445
int nptbins = 0;
@@ -482,7 +484,7 @@ struct FlowGFWOmegaXi {
482484
}
483485

484486
// input shared_ptr<TProfile3D>
485-
void FillProfilepTMass(const GFW::CorrConfig& corrconf, std::shared_ptr<TProfile3D> TProfile3D, const int& ptbin, const int& PDGCode, const float& cent)
487+
void fillProfilepTMass(const GFW::CorrConfig& corrconf, std::shared_ptr<TProfile3D> TProfile3D, const int& ptbin, const int& PDGCode, const float& cent)
486488
{
487489
int nMassBins = 0;
488490
int nptbins = 0;
@@ -634,7 +636,7 @@ struct FlowGFWOmegaXi {
634636
return true;
635637
}
636638

637-
void process(aodCollisions::iterator const& collision, aod::BCsWithTimestamps const&, aodTracks const& tracks, aod::CascDataExt const& Cascades, aod::V0Datas const& V0s, DaughterTracks&)
639+
void process(aodCollisions::iterator const& collision, aod::BCsWithTimestamps const&, aodTracks const& tracks, aod::CascDataExt const& Cascades, aod::V0Datas const& V0s, daughterTracks const&)
638640
{
639641
int Ntot = tracks.size();
640642
int CandNum_all[4] = {0, 0, 0, 0};
@@ -663,7 +665,7 @@ struct FlowGFWOmegaXi {
663665
float weff = 1;
664666
float wacc = 1;
665667
// fill GFW ref flow
666-
for (auto& track : tracks) {
668+
for (const auto& track : tracks) {
667669
if (!setCurrentParticleWeights(weff, wacc, track, vtxz, 0))
668670
continue;
669671
registry.fill(HIST("hPhi"), track.phi());
@@ -682,9 +684,9 @@ struct FlowGFWOmegaXi {
682684
fWeightsREF->Fill(track.phi(), track.eta(), vtxz, track.pt(), cent, 0);
683685
}
684686
// fill GFW of V0 flow
685-
for (auto& v0 : V0s) {
686-
auto v0posdau = v0.posTrack_as<DaughterTracks>();
687-
auto v0negdau = v0.negTrack_as<DaughterTracks>();
687+
for (const auto& v0 : V0s) {
688+
auto v0posdau = v0.posTrack_as<daughterTracks>();
689+
auto v0negdau = v0.negTrack_as<daughterTracks>();
688690
// check tpc
689691
bool isK0s = false;
690692
bool isLambda = false;
@@ -777,10 +779,10 @@ struct FlowGFWOmegaXi {
777779
}
778780
}
779781
// fill GFW of casc flow
780-
for (auto& casc : Cascades) {
781-
auto bachelor = casc.bachelor_as<DaughterTracks>();
782-
auto posdau = casc.posTrack_as<DaughterTracks>();
783-
auto negdau = casc.negTrack_as<DaughterTracks>();
782+
for (const auto& casc : Cascades) {
783+
auto bachelor = casc.bachelor_as<daughterTracks>();
784+
auto posdau = casc.posTrack_as<daughterTracks>();
785+
auto negdau = casc.negTrack_as<daughterTracks>();
784786
// check TPC
785787
if (cfgcheckDauTPC && (!posdau.hasTPC() || !negdau.hasTPC() || !bachelor.hasTPC())) {
786788
continue;
@@ -879,56 +881,56 @@ struct FlowGFWOmegaXi {
879881
}
880882
}
881883
// Filling cumulant with ROOT TProfile and loop for all ptBins
882-
FillProfile(corrconfigs.at(15), HIST("c22"), cent);
883-
FillProfile(corrconfigs.at(16), HIST("c24"), cent);
884-
FillProfile(corrconfigs.at(17), HIST("K0sc22"), cent);
885-
FillProfile(corrconfigs.at(18), HIST("K0sc22"), cent);
886-
FillProfile(corrconfigs.at(19), HIST("Lambdac22"), cent);
887-
FillProfile(corrconfigs.at(20), HIST("Lambdac22"), cent);
884+
fillProfile(corrconfigs.at(15), HIST("c22"), cent);
885+
fillProfile(corrconfigs.at(16), HIST("c24"), cent);
886+
fillProfile(corrconfigs.at(17), HIST("K0sc22"), cent);
887+
fillProfile(corrconfigs.at(18), HIST("K0sc22"), cent);
888+
fillProfile(corrconfigs.at(19), HIST("Lambdac22"), cent);
889+
fillProfile(corrconfigs.at(20), HIST("Lambdac22"), cent);
888890
for (int i = 1; i <= nPtBins; i++) {
889-
FillProfilepT(corrconfigs.at(0), HIST("c22dpt"), i, cent);
890-
FillProfilepT(corrconfigs.at(1), HIST("c22dpt"), i, cent);
891-
FillProfilepT(corrconfigs.at(2), HIST("c24dpt"), i, cent);
891+
fillProfilepT(corrconfigs.at(0), HIST("c22dpt"), i, cent);
892+
fillProfilepT(corrconfigs.at(1), HIST("c22dpt"), i, cent);
893+
fillProfilepT(corrconfigs.at(2), HIST("c24dpt"), i, cent);
892894
}
893895
for (int i = 1; i <= nV0PtBins; i++) {
894-
FillProfilepTMass(corrconfigs.at(9), HIST("K0sc22dpt"), i, kK0Short, cent);
895-
FillProfilepTMass(corrconfigs.at(10), HIST("K0sc22dpt"), i, kK0Short, cent);
896-
FillProfilepTMass(corrconfigs.at(11), HIST("K0sc24dpt"), i, kK0Short, cent);
897-
FillProfilepTMass(corrconfigs.at(12), HIST("Lambdac22dpt"), i, kLambda0, cent);
898-
FillProfilepTMass(corrconfigs.at(13), HIST("Lambdac22dpt"), i, kLambda0, cent);
899-
FillProfilepTMass(corrconfigs.at(14), HIST("Lambdac24dpt"), i, kLambda0, cent);
896+
fillProfilepTMass(corrconfigs.at(9), HIST("K0sc22dpt"), i, kK0Short, cent);
897+
fillProfilepTMass(corrconfigs.at(10), HIST("K0sc22dpt"), i, kK0Short, cent);
898+
fillProfilepTMass(corrconfigs.at(11), HIST("K0sc24dpt"), i, kK0Short, cent);
899+
fillProfilepTMass(corrconfigs.at(12), HIST("Lambdac22dpt"), i, kLambda0, cent);
900+
fillProfilepTMass(corrconfigs.at(13), HIST("Lambdac22dpt"), i, kLambda0, cent);
901+
fillProfilepTMass(corrconfigs.at(14), HIST("Lambdac24dpt"), i, kLambda0, cent);
900902
}
901903
for (int i = 1; i <= nXiPtBins; i++) {
902-
FillProfilepTMass(corrconfigs.at(3), HIST("Xic22dpt"), i, kXiMinus, cent);
903-
FillProfilepTMass(corrconfigs.at(4), HIST("Xic22dpt"), i, kXiMinus, cent);
904-
FillProfilepTMass(corrconfigs.at(5), HIST("Xic24dpt"), i, kXiMinus, cent);
905-
FillProfilepTMass(corrconfigs.at(6), HIST("Omegac22dpt"), i, kOmegaMinus, cent);
906-
FillProfilepTMass(corrconfigs.at(7), HIST("Omegac22dpt"), i, kOmegaMinus, cent);
907-
FillProfilepTMass(corrconfigs.at(8), HIST("Omegac24dpt"), i, kOmegaMinus, cent);
904+
fillProfilepTMass(corrconfigs.at(3), HIST("Xic22dpt"), i, kXiMinus, cent);
905+
fillProfilepTMass(corrconfigs.at(4), HIST("Xic22dpt"), i, kXiMinus, cent);
906+
fillProfilepTMass(corrconfigs.at(5), HIST("Xic24dpt"), i, kXiMinus, cent);
907+
fillProfilepTMass(corrconfigs.at(6), HIST("Omegac22dpt"), i, kOmegaMinus, cent);
908+
fillProfilepTMass(corrconfigs.at(7), HIST("Omegac22dpt"), i, kOmegaMinus, cent);
909+
fillProfilepTMass(corrconfigs.at(8), HIST("Omegac24dpt"), i, kOmegaMinus, cent);
908910
}
909911
// Fill subevents flow
910912
TRandom3* fRdm = new TRandom3(0);
911913
double Eventrdm = 10 * fRdm->Rndm();
912914
for (int j = 1; j <= 10; j++) {
913915
if (Eventrdm > (j - 1) && Eventrdm < j)
914916
continue;
915-
FillProfile(corrconfigs.at(15), REFc22[j - 1], cent);
916-
FillProfile(corrconfigs.at(16), REFc24[j - 1], cent);
917+
fillProfile(corrconfigs.at(15), refc22[j - 1], cent);
918+
fillProfile(corrconfigs.at(16), refc24[j - 1], cent);
917919
for (int i = 1; i <= nV0PtBins; i++) {
918-
FillProfilepTMass(corrconfigs.at(9), K0sc22[j - 1], i, kK0Short, cent);
919-
FillProfilepTMass(corrconfigs.at(10), K0sc22[j - 1], i, kK0Short, cent);
920-
FillProfilepTMass(corrconfigs.at(11), K0sc24[j - 1], i, kK0Short, cent);
921-
FillProfilepTMass(corrconfigs.at(12), Lambdac22[j - 1], i, kLambda0, cent);
922-
FillProfilepTMass(corrconfigs.at(13), Lambdac22[j - 1], i, kLambda0, cent);
923-
FillProfilepTMass(corrconfigs.at(14), Lambdac24[j - 1], i, kLambda0, cent);
920+
fillProfilepTMass(corrconfigs.at(9), k0sc22[j - 1], i, kK0Short, cent);
921+
fillProfilepTMass(corrconfigs.at(10), k0sc22[j - 1], i, kK0Short, cent);
922+
fillProfilepTMass(corrconfigs.at(11), k0sc24[j - 1], i, kK0Short, cent);
923+
fillProfilepTMass(corrconfigs.at(12), lambdac22[j - 1], i, kLambda0, cent);
924+
fillProfilepTMass(corrconfigs.at(13), lambdac22[j - 1], i, kLambda0, cent);
925+
fillProfilepTMass(corrconfigs.at(14), lambdac24[j - 1], i, kLambda0, cent);
924926
}
925927
for (int i = 1; i <= nXiPtBins; i++) {
926-
FillProfilepTMass(corrconfigs.at(3), Xic22[j - 1], i, kXiMinus, cent);
927-
FillProfilepTMass(corrconfigs.at(4), Xic22[j - 1], i, kXiMinus, cent);
928-
FillProfilepTMass(corrconfigs.at(5), Xic24[j - 1], i, kXiMinus, cent);
929-
FillProfilepTMass(corrconfigs.at(6), Omegac22[j - 1], i, kOmegaMinus, cent);
930-
FillProfilepTMass(corrconfigs.at(7), Omegac22[j - 1], i, kOmegaMinus, cent);
931-
FillProfilepTMass(corrconfigs.at(8), Omegac24[j - 1], i, kOmegaMinus, cent);
928+
fillProfilepTMass(corrconfigs.at(3), xic22[j - 1], i, kXiMinus, cent);
929+
fillProfilepTMass(corrconfigs.at(4), xic22[j - 1], i, kXiMinus, cent);
930+
fillProfilepTMass(corrconfigs.at(5), xic24[j - 1], i, kXiMinus, cent);
931+
fillProfilepTMass(corrconfigs.at(6), omegac22[j - 1], i, kOmegaMinus, cent);
932+
fillProfilepTMass(corrconfigs.at(7), omegac22[j - 1], i, kOmegaMinus, cent);
933+
fillProfilepTMass(corrconfigs.at(8), omegac24[j - 1], i, kOmegaMinus, cent);
932934
}
933935
}
934936
}

0 commit comments

Comments
 (0)