Skip to content

Commit 4993604

Browse files
prottayCMTProttay Das
prottayCMT
and
Prottay Das
authored
[PWGLF] added process function to check effect of spatial calibration on v1 (#18015)
Co-authored-by: Prottay Das <prottay@alipap1.cern.ch>
1 parent 7e385d4 commit 4993604

1 file changed

Lines changed: 139 additions & 73 deletions

File tree

‎PWGLF/Tasks/Strangeness/lambdapolsp.cxx‎

Lines changed: 139 additions & 73 deletions
Original file line numberDiff line numberDiff line change
@@ -14,6 +14,7 @@
1414
#include "PWGLF/DataModel/LFStrangenessPIDTables.h"
1515
#include "PWGLF/DataModel/LFStrangenessTables.h"
1616
#include "PWGLF/DataModel/SPCalibrationTables.h"
17+
#include "PWGLF/DataModel/ZDC2StageCalibrationTables.h"
1718

1819
#include "Common/CCDB/EventSelectionParams.h"
1920
#include "Common/CCDB/RCTSelectionFlags.h"
@@ -322,8 +323,6 @@ struct lambdapolsp {
322323
AxisSpec thnAxisres{binGrp.resNbins, binGrp.lbinres, binGrp.hbinres, "Reso"};
323324
AxisSpec thnAxisInvMass{binGrp.IMNbins, binGrp.lbinIM, binGrp.hbinIM, "#it{M} (GeV/#it{c}^{2})"};
324325
AxisSpec spAxis = {binGrp.spNbins, binGrp.lbinsp, binGrp.hbinsp, "Sp"};
325-
// AxisSpec qxZDCAxis = {binGrp.QxyNbins, binGrp.lbinQxy, binGrp.hbinQxy, "Qx"};
326-
// AxisSpec centAxis = {CentNbins, lbinCent, hbinCent, "V0M (%)"};
327326

328327
std::vector<AxisSpec> runaxes = {thnAxisInvMass, axisGrp.configthnAxispT, axisGrp.configthnAxisPol, axisGrp.configcentAxis};
329328
if (needetaaxis)
@@ -342,10 +341,6 @@ struct lambdapolsp {
342341

343342
runaxesSyst.push_back(systIDAxis);
344343

345-
// if (needetaaxis)
346-
// runaxes.insert(runaxes.end(), {axisGrp.configbinAxis});
347-
// std::vector<AxisSpec> runaxes2 = {thnAxisInvMass, axisGrp.configthnAxispT, axisGrp.configcentAxis};
348-
349344
if (checkwithpub) {
350345
if (useprofile == 2) {
351346
histos.add("hpuxQxpvscentpteta", "hpuxQxpvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
@@ -356,49 +351,8 @@ struct lambdapolsp {
356351
histos.add("hpuxyQxypvscentpteta", "hpuxyQxypvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
357352
histos.add("hpoddv1vscentpteta", "hpoddv1vscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
358353
histos.add("hpevenv1vscentpteta", "hpevenv1vscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
359-
/*histos.add("hpv21", "hpv21", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
360-
histos.add("hpv22", "hpv22", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
361-
histos.add("hpv23", "hpv23", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
362-
histos.add("hpx2Tx1Ax1Cvscentpteta", "hpx2Tx1Ax1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
363-
histos.add("hpx2Ty1Ay1Cvscentpteta", "hpx2Ty1Ay1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
364-
histos.add("hpy2Tx1Ay1Cvscentpteta", "hpy2Tx1Ay1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
365-
histos.add("hpy2Ty1Ax1Cvscentpteta", "hpy2Ty1Ax1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
366-
histos.add("hpx1Ax1Cvscentpteta", "hpx1Ax1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
367-
histos.add("hpy1Ay1Cvscentpteta", "hpy1Ay1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
368-
histos.add("hpx1Avscentpteta", "hpx1Avscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
369-
histos.add("hpx1Cvscentpteta", "hpx1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
370-
histos.add("hpy1Avscentpteta", "hpy1Avscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
371-
histos.add("hpy1Cvscentpteta", "hpy1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
372-
373-
histos.add("hpx2Tx1Avscentpteta", "hpx2Tx1Avscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
374-
histos.add("hpx2Tx1Cvscentpteta", "hpx2Tx1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
375-
histos.add("hpx2Ty1Avscentpteta", "hpx2Ty1Avscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
376-
histos.add("hpx2Ty1Cvscentpteta", "hpx2Ty1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
377-
histos.add("hpy2Tx1Avscentpteta", "hpy2Tx1Avscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
378-
histos.add("hpy2Ty1Cvscentpteta", "hpy2Ty1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
379-
histos.add("hpy2Ty1Avscentpteta", "hpy2Ty1Avscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
380-
histos.add("hpy2Tx1Cvscentpteta", "hpy2Tx1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
381-
histos.add("hpx1Ay1Cvscentpteta", "hpx1Ay1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
382-
histos.add("hpy1Ax1Cvscentpteta", "hpy1Ax1Cvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
383-
histos.add("hpx2Tvscentpteta", "hpx2Tvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
384-
histos.add("hpy2Tvscentpteta", "hpy2Tvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
385-
*/
386354
histos.add("hpuxvscentpteta", "hpuxvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
387355
histos.add("hpuyvscentpteta", "hpuyvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
388-
/*
389-
histos.add("hpuxvscentptetaneg", "hpuxvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
390-
histos.add("hpuyvscentptetaneg", "hpuyvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
391-
392-
histos.add("hpuxQxpvscentptetaneg", "hpuxQxpvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
393-
histos.add("hpuyQypvscentptetaneg", "hpuyQypvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
394-
histos.add("hpuxQxtvscentptetaneg", "hpuxQxtvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
395-
histos.add("hpuyQytvscentptetaneg", "hpuyQytvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
396-
histos.add("hpuxyQxytvscentptetaneg", "hpuxyQxytvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
397-
histos.add("hpuxyQxypvscentptetaneg", "hpuxyQxypvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
398-
histos.add("hpoddv1vscentptetaneg", "hpoddv1vscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
399-
histos.add("hpevenv1vscentptetaneg", "hpevenv1vscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, configthnAxispT, configetaAxis, spAxis}, true);
400-
*/
401-
402356
histos.add("hpQxtQxpvscent", "hpQxtQxpvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true);
403357
histos.add("hpQytQypvscent", "hpQytQypvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true);
404358
histos.add("hpQxytpvscent", "hpQxytpvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true);
@@ -446,17 +400,6 @@ struct lambdapolsp {
446400

447401
histos.add("hpuxvscentpteta", "hpuxvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
448402
histos.add("hpuyvscentpteta", "hpuyvscentpteta", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
449-
/*histos.add("hpuxvscentptetaneg", "hpuxvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
450-
histos.add("hpuyvscentptetaneg", "hpuyvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
451-
452-
histos.add("hpuxQxpvscentptetaneg", "hpuxQxpvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
453-
histos.add("hpuyQypvscentptetaneg", "hpuyQypvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
454-
histos.add("hpuxQxtvscentptetaneg", "hpuxQxtvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
455-
histos.add("hpuyQytvscentptetaneg", "hpuyQytvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
456-
histos.add("hpuxyQxytvscentptetaneg", "hpuxyQxytvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
457-
histos.add("hpuxyQxypvscentptetaneg", "hpuxyQxypvscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
458-
histos.add("hpoddv1vscentptetaneg", "hpoddv1vscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);
459-
histos.add("hpevenv1vscentptetaneg", "hpevenv1vscentptetaneg", HistType::kTHnSparseF, {axisGrp.configcentAxis, axisGrp.configthnAxispT, axisGrp.configetaAxis, spAxis}, true);*/
460403

461404
histos.add("hpQxtQxpvscent", "hpQxtQxpvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true);
462405
histos.add("hpQytQypvscent", "hpQytQypvscent", HistType::kTHnSparseF, {axisGrp.configcentAxis, spAxis}, true);
@@ -472,14 +415,6 @@ struct lambdapolsp {
472415
}
473416

474417
histos.add("hCentrality", "Centrality distribution", kTH1F, {{axisGrp.configcentAxis}});
475-
// histos.add("hpsiApsiC", "hpsiApsiC", kTHnSparseF, {psiACAxis, psiACAxis});
476-
// histos.add("hpsiApsiC", "hpsiApsiC", kTH2F, {psiACAxis, psiACAxis});
477-
// histos.add("hphiminuspsiA", "hphiminuspisA", kTH1F, {{50, 0, 6.28}}, true);
478-
// histos.add("hphiminuspsiC", "hphiminuspisC", kTH1F, {{50, 0, 6.28}}, true);
479-
// histos.add("hCentrality0", "Centrality distribution0", kTH1F, {{centAxis}});
480-
// histos.add("hCentrality1", "Centrality distribution1", kTH1F, {{centAxis}});
481-
// histos.add("hCentrality2", "Centrality distribution2", kTH1F, {{centAxis}});
482-
// histos.add("hCentrality3", "Centrality distribution3", kTH1F, {{centAxis}});
483418

484419
if (!checkwithpub) {
485420
// histos.add("hVtxZ", "Vertex distribution in Z;Z (cm)", kTH1F, {{20, -10.0, 10.0}});
@@ -1388,20 +1323,19 @@ struct lambdapolsp {
13881323
}
13891324

13901325
ROOT::Math::PxPyPzMVector Lambda, AntiLambda, Lambdadummy, AntiLambdadummy, Proton, Pion, AntiProton, AntiPion, fourVecDauCM, K0sdummy, K0s;
1391-
// double phiangle = 0.0;
1392-
// double angleLambda=0.0;
1393-
// double angleAntiLambda=0.0;
1394-
double massLambda = o2::constants::physics::MassLambda;
1395-
double massK0s = o2::constants::physics::MassK0Short;
1396-
double massPr = o2::constants::physics::MassProton;
1397-
double massPi = o2::constants::physics::MassPionCharged;
1326+
1327+
static constexpr double massLambda = o2::constants::physics::MassLambda;
1328+
static constexpr double massK0s = o2::constants::physics::MassK0Short;
1329+
static constexpr double massPr = o2::constants::physics::MassProton;
1330+
static constexpr double massPi = o2::constants::physics::MassPionCharged;
13981331

13991332
Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
14001333
Filter centralityFilter = (nabs(aod::cent::centFT0C) < cfgCutCentralityMax && nabs(aod::cent::centFT0C) > cfgCutCentralityMin);
14011334
Filter acceptanceFilter = (nabs(aod::track::eta) < cfgCutEta && nabs(aod::track::pt) > cfgCutPT);
14021335
Filter dcaCutFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz);
14031336

14041337
using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::SPCalibrationTables, aod::Mults>>;
1338+
using EventCandidatescalibstage = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::ZDC2StageCalibs, aod::Mults>>;
14051339
using EventCandidatesMC = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs>>;
14061340
using AllTrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTPCFullPr, aod::pidTPCFullKa>>;
14071341
using ResoV0s = aod::V0Datas;
@@ -1414,6 +1348,138 @@ struct lambdapolsp {
14141348

14151349
using BCsRun3 = soa::Join<aod::BCsWithTimestamps, aod::Run3MatchedToBCSparse>;
14161350

1351+
void processDatav1(EventCandidatescalibstage::iterator const& collision, AllTrackCandidates const& tracks, BCsRun3 const&)
1352+
{
1353+
1354+
if (!collision.sel8()) {
1355+
return;
1356+
}
1357+
double centrality = -999.;
1358+
if (centestim == 0)
1359+
centrality = collision.centFT0C();
1360+
1361+
if (!collision.triggerSP()) {
1362+
return;
1363+
}
1364+
1365+
if (evselGrp.additionalEvSel && (!collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
1366+
return;
1367+
}
1368+
if (evselGrp.additionalEvSel2 && (collision.trackOccupancyInTimeRange() > evselGrp.cfgMaxOccupancy || collision.trackOccupancyInTimeRange() < evselGrp.cfgMinOccupancy)) {
1369+
return;
1370+
}
1371+
if (evselGrp.additionalEvSel4 && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
1372+
return;
1373+
}
1374+
1375+
if (rctCut.requireRCTFlagChecker && !rctChecker(collision)) {
1376+
return;
1377+
}
1378+
1379+
auto qxZDCA = collision.qxZDCA();
1380+
auto qxZDCC = collision.qxZDCC();
1381+
auto qyZDCA = collision.qyZDCA();
1382+
auto qyZDCC = collision.qyZDCC();
1383+
1384+
histos.fill(HIST("hCentrality"), centrality);
1385+
1386+
auto QxtQxp = qxZDCA * qxZDCC;
1387+
auto QytQyp = qyZDCA * qyZDCC;
1388+
auto Qxytp = QxtQxp + QytQyp;
1389+
auto QxpQyt = qxZDCA * qyZDCC;
1390+
auto QxtQyp = qxZDCC * qyZDCA;
1391+
1392+
histos.fill(HIST("hpQxtQxpvscent"), centrality, QxtQxp);
1393+
histos.fill(HIST("hpQytQypvscent"), centrality, QytQyp);
1394+
histos.fill(HIST("hpQxytpvscent"), centrality, Qxytp);
1395+
histos.fill(HIST("hpQxpQytvscent"), centrality, QxpQyt);
1396+
histos.fill(HIST("hpQxtQypvscent"), centrality, QxtQyp);
1397+
1398+
int wNUA = 1;
1399+
1400+
for (const auto& track : tracks) {
1401+
if (!selectionTrack(track)) {
1402+
continue;
1403+
}
1404+
1405+
float sign = track.sign();
1406+
if (sign == 0.0) // removing neutral particles
1407+
continue;
1408+
1409+
auto ux = TMath::Cos(GetPhiInRange(track.phi()));
1410+
auto uy = TMath::Sin(GetPhiInRange(track.phi()));
1411+
1412+
auto uxQxp = ux * qxZDCA;
1413+
auto uyQyp = uy * qyZDCA;
1414+
auto uxyQxyp = uxQxp + uyQyp;
1415+
auto uxQxt = ux * qxZDCC;
1416+
auto uyQyt = uy * qyZDCC;
1417+
auto uxyQxyt = uxQxt + uyQyt;
1418+
auto oddv1 = ux * (qxZDCA - qxZDCC) + uy * (qyZDCA - qyZDCC);
1419+
auto evenv1 = ux * (qxZDCA + qxZDCC) + uy * (qyZDCA + qyZDCC);
1420+
1421+
if (globalpt) {
1422+
histos.fill(HIST("hpuxQxpvscentpteta"), centrality, track.pt(), track.eta(), uxQxp, wNUA);
1423+
histos.fill(HIST("hpuyQypvscentpteta"), centrality, track.pt(), track.eta(), uyQyp, wNUA);
1424+
histos.fill(HIST("hpuxQxtvscentpteta"), centrality, track.pt(), track.eta(), uxQxt, wNUA);
1425+
histos.fill(HIST("hpuyQytvscentpteta"), centrality, track.pt(), track.eta(), uyQyt, wNUA);
1426+
1427+
histos.fill(HIST("hpuxvscentpteta"), centrality, track.pt(), track.eta(), ux, wNUA);
1428+
histos.fill(HIST("hpuyvscentpteta"), centrality, track.pt(), track.eta(), uy, wNUA);
1429+
1430+
histos.fill(HIST("hpuxyQxytvscentpteta"), centrality, track.pt(), track.eta(), uxyQxyt, wNUA);
1431+
histos.fill(HIST("hpuxyQxypvscentpteta"), centrality, track.pt(), track.eta(), uxyQxyp, wNUA);
1432+
histos.fill(HIST("hpoddv1vscentpteta"), centrality, track.pt(), track.eta(), oddv1, wNUA);
1433+
histos.fill(HIST("hpevenv1vscentpteta"), centrality, track.pt(), track.eta(), evenv1, wNUA);
1434+
1435+
histos.fill(HIST("hpQxtQxpvscentptetaTrack"),
1436+
centrality, track.pt(), track.eta(), QxtQxp, wNUA);
1437+
1438+
histos.fill(HIST("hpQytQypvscentptetaTrack"),
1439+
centrality, track.pt(), track.eta(), QytQyp, wNUA);
1440+
1441+
histos.fill(HIST("hpQxytpvscentptetaTrack"),
1442+
centrality, track.pt(), track.eta(), Qxytp, wNUA);
1443+
1444+
histos.fill(HIST("hpQxpQytvscentptetaTrack"),
1445+
centrality, track.pt(), track.eta(), QxpQyt, wNUA);
1446+
1447+
histos.fill(HIST("hpQxtQypvscentptetaTrack"),
1448+
centrality, track.pt(), track.eta(), QxtQyp, wNUA);
1449+
} else {
1450+
histos.fill(HIST("hpuxQxpvscentpteta"), centrality, track.tpcInnerParam(), track.eta(), uxQxp, wNUA);
1451+
histos.fill(HIST("hpuyQypvscentpteta"), centrality, track.tpcInnerParam(), track.eta(), uyQyp, wNUA);
1452+
histos.fill(HIST("hpuxQxtvscentpteta"), centrality, track.tpcInnerParam(), track.eta(), uxQxt, wNUA);
1453+
histos.fill(HIST("hpuyQytvscentpteta"), centrality, track.tpcInnerParam(), track.eta(), uyQyt, wNUA);
1454+
1455+
histos.fill(HIST("hpuxvscentpteta"), centrality, track.pt(), track.eta(), ux, wNUA);
1456+
histos.fill(HIST("hpuyvscentpteta"), centrality, track.pt(), track.eta(), uy, wNUA);
1457+
1458+
histos.fill(HIST("hpuxyQxytvscentpteta"), centrality, track.tpcInnerParam(), track.eta(), uxyQxyt, wNUA);
1459+
histos.fill(HIST("hpuxyQxypvscentpteta"), centrality, track.tpcInnerParam(), track.eta(), uxyQxyp, wNUA);
1460+
histos.fill(HIST("hpoddv1vscentpteta"), centrality, track.pt(), track.eta(), oddv1, wNUA);
1461+
histos.fill(HIST("hpevenv1vscentpteta"), centrality, track.pt(), track.eta(), evenv1, wNUA);
1462+
1463+
histos.fill(HIST("hpQxtQxpvscentptetaTrack"),
1464+
centrality, track.pt(), track.eta(), QxtQxp, wNUA);
1465+
1466+
histos.fill(HIST("hpQytQypvscentptetaTrack"),
1467+
centrality, track.pt(), track.eta(), QytQyp, wNUA);
1468+
1469+
histos.fill(HIST("hpQxytpvscentptetaTrack"),
1470+
centrality, track.pt(), track.eta(), Qxytp, wNUA);
1471+
1472+
histos.fill(HIST("hpQxpQytvscentptetaTrack"),
1473+
centrality, track.pt(), track.eta(), QxpQyt, wNUA);
1474+
1475+
histos.fill(HIST("hpQxtQypvscentptetaTrack"),
1476+
centrality, track.pt(), track.eta(), QxtQyp, wNUA);
1477+
}
1478+
}
1479+
}
1480+
1481+
PROCESS_SWITCH(lambdapolsp, processDatav1, "Process datav1", false);
1482+
14171483
void processData(EventCandidates::iterator const& collision, AllTrackCandidates const& tracks, ResoV0s const& V0s, BCsRun3 const&)
14181484
{
14191485

0 commit comments

Comments
 (0)