Skip to content

Commit 02ba036

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent 6f71560 commit 02ba036

File tree

1 file changed

+26
-30
lines changed

1 file changed

+26
-30
lines changed

‎PWGCF/Flow/Tasks/flowSP.cxx‎

Lines changed: 26 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -150,7 +150,7 @@ struct FlowSP {
150150
Configurable<std::vector<double>> cfgEvSelsMult{"cfgEvSelsMult", std::vector<double>{1301.56, -41.4615, 0.478224, -0.00239449, 4.46966e-06, 2967.6, -102.927, 1.47488, -0.0106534, 3.28622e-05}, "Multiplicity cuts (Global) first 5 parameters cutLOW last 5 cutHIGH (Default is +-2sigma pass5) "};
151151

152152
Filter collisionFilter = nabs(aod::collision::posZ) < cfgEvSelsVtxZ;
153-
Filter trackFilter = nabs(aod::track::eta) < cfgTrackSelsEta && aod::track::pt > cfgTrackSelsPtmin&& aod::track::pt < cfgTrackSelsPtmax && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t) true)) && nabs(aod::track::dcaXY) < cfgTrackSelsDCAxy && nabs(aod::track::dcaZ) < cfgTrackSelsDCAz;
153+
Filter trackFilter = nabs(aod::track::eta) < cfgTrackSelsEta && aod::track::pt > cfgTrackSelsPtmin&& aod::track::pt < cfgTrackSelsPtmax && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t) true)) && nabs(aod::track::dcaXY) < cfgTrackSelsDCAxy&& nabs(aod::track::dcaZ) < cfgTrackSelsDCAz;
154154
Filter trackFilterMC = nabs(aod::mcparticle::eta) < cfgTrackSelsEta && aod::mcparticle::pt > cfgTrackSelsPtmin&& aod::mcparticle::pt < cfgTrackSelsPtmax;
155155
using GeneralCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Cs, aod::CentFT0CVariant1s, aod::CentFT0Ms, aod::CentFV0As, aod::CentNGlobals>;
156156
using UnfilteredTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelection, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr, aod::pidTOFbeta, aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>;
@@ -523,7 +523,7 @@ struct FlowSP {
523523
registry.add<TProfile3D>("incl/vnC", "", kTProfile3D, {axisPt, axisEtaVn, axisCentrality});
524524
registry.add<TProfile3D>("incl/vnA", "", kTProfile3D, {axisPt, axisEtaVn, axisCentrality});
525525
}
526-
if(cfgFillMeanPT) {
526+
if (cfgFillMeanPT) {
527527
registry.add<TProfile2D>("incl/meanPT/meanRelPtA", "", kTProfile2D, {axisEtaVn, axisCentrality});
528528
registry.add<TProfile2D>("incl/meanPT/meanRelPtC", "", kTProfile2D, {axisEtaVn, axisCentrality});
529529
}
@@ -791,9 +791,9 @@ struct FlowSP {
791791
int sizeEff = cfg.mEfficiency.size();
792792
if (sizeEff > pID) {
793793
if (cfgUseNUE2D) {
794-
int binx;
795-
int biny;
796-
if(cfgUseNUE2Deta) {
794+
int binx;
795+
int biny;
796+
if (cfgUseNUE2Deta) {
797797
biny = cfg.mEfficiency2D[pID]->GetYaxis()->FindBin(pt);
798798
binx = cfg.mEfficiency2D[pID]->GetXaxis()->FindBin(eta);
799799
} else {
@@ -1220,8 +1220,6 @@ struct FlowSP {
12201220
return;
12211221
histos.fill(HIST("hEventCount"), evSel_isSelectedZDC);
12221222

1223-
1224-
12251223
spm.qxA = collision.qxA();
12261224
spm.qyA = collision.qyA();
12271225
spm.qxC = collision.qxC();
@@ -1308,15 +1306,13 @@ struct FlowSP {
13081306
float trackPx = 0;
13091307
int trackSize = 0;
13101308

1311-
TProfile2D* hRelEtaPt = new TProfile2D("hRelEtaPt", "hRelEtaPt", 8,-.8,.8, 3,0,3);
1312-
TProfile2D* ptV1A = new TProfile2D("ptV1A", "ptV1A", 8,-.8,.8, 3,0,3);
1313-
TProfile2D* ptV1C = new TProfile2D("ptV1C", "ptV1C", 8,-.8,.8, 3,0,3);
1309+
TProfile2D* hRelEtaPt = new TProfile2D("hRelEtaPt", "hRelEtaPt", 8, -.8, .8, 3, 0, 3);
1310+
TProfile2D* ptV1A = new TProfile2D("ptV1A", "ptV1A", 8, -.8, .8, 3, 0, 3);
1311+
TProfile2D* ptV1C = new TProfile2D("ptV1C", "ptV1C", 8, -.8, .8, 3, 0, 3);
13141312

1315-
double sumPt = 0;
1313+
double sumPt = 0;
13161314
int ptCounter = 0;
13171315

1318-
1319-
13201316
for (const auto& track : tracks) {
13211317

13221318
ParticleType trackPID = (cfgFillPID || cfgFillPIDQA) ? getTrackPID(track) : kUnidentified;
@@ -1391,7 +1387,6 @@ struct FlowSP {
13911387
ptV1C->Fill(track.eta(), kInclusive, track.pt() * ((spm.uy * spm.qyC + spm.ux * spm.qxC) / std::sqrt(std::fabs(spm.corrQQ))), weight);
13921388
ptV1C->Fill(track.eta(), spm.charge, track.pt() * ((spm.uy * spm.qyC + spm.ux * spm.qxC) / std::sqrt(std::fabs(spm.corrQQ))), weight_charged);
13931389

1394-
13951390
fillAllQA<kAfter, kUnidentified>(track);
13961391
if (cfgFillPIDQA) {
13971392
switch (trackPID) {
@@ -1490,10 +1485,9 @@ struct FlowSP {
14901485

14911486
} // end of track loop
14921487

1493-
1494-
// Now we want to fill the final relPt histogram
1488+
// Now we want to fill the final relPt histogram
14951489
// Loop over all eta and fill bins
1496-
if(cfgFillMeanPT) {
1490+
if (cfgFillMeanPT) {
14971491
for (int i = 0; i < hRelEtaPt->GetNbinsX(); i++) {
14981492
double eta = hRelEtaPt->GetXaxis()->GetBinCenter(i);
14991493
int nEntries = hRelEtaPt->GetEntries();
@@ -1503,32 +1497,34 @@ struct FlowSP {
15031497
double drelPt = hRelEtaPt->GetBinContent(bin);
15041498
double dptV1A = ptV1A->GetBinContent(bin);
15051499
double dptV1C = ptV1C->GetBinContent(bin);
1506-
if(drelPt) registry.fill(HIST("incl/meanPT/meanRelPtA"), eta, spm.centrality, dptV1A / drelPt,1);
1507-
if(drelPt) registry.fill(HIST("incl/meanPT/meanRelPtC"), eta, spm.centrality, dptV1C / drelPt,1);
1500+
if (drelPt)
1501+
registry.fill(HIST("incl/meanPT/meanRelPtA"), eta, spm.centrality, dptV1A / drelPt, 1);
1502+
if (drelPt)
1503+
registry.fill(HIST("incl/meanPT/meanRelPtC"), eta, spm.centrality, dptV1C / drelPt, 1);
15081504

15091505
bin = hRelEtaPt->FindBin(eta, kPositive);
15101506
double drelPt_pos = hRelEtaPt->GetBinContent(bin);
15111507
double dptV1A_pos = ptV1A->GetBinContent(bin);
15121508
double dptV1C_pos = ptV1C->GetBinContent(bin);
1513-
if(drelPt_pos) registry.fill(HIST("pos/meanPT/meanRelPtA"), eta, spm.centrality, dptV1A_pos / drelPt_pos,1);
1514-
if(drelPt_pos) registry.fill(HIST("pos/meanPT/meanRelPtC"), eta, spm.centrality, dptV1C_pos / drelPt_pos,1);
1509+
if (drelPt_pos)
1510+
registry.fill(HIST("pos/meanPT/meanRelPtA"), eta, spm.centrality, dptV1A_pos / drelPt_pos, 1);
1511+
if (drelPt_pos)
1512+
registry.fill(HIST("pos/meanPT/meanRelPtC"), eta, spm.centrality, dptV1C_pos / drelPt_pos, 1);
15151513

15161514
bin = hRelEtaPt->FindBin(eta, kNegative);
15171515
double drelPt_neg = hRelEtaPt->GetBinContent(bin);
15181516
double dptV1A_neg = ptV1A->GetBinContent(bin);
15191517
double dptV1C_neg = ptV1C->GetBinContent(bin);
1520-
if(drelPt_neg) registry.fill(HIST("neg/meanPT/meanRelPtA"), eta, spm.centrality, dptV1A_neg / drelPt_neg,1);
1521-
if(drelPt_neg) registry.fill(HIST("neg/meanPT/meanRelPtC"), eta, spm.centrality, dptV1C_neg / drelPt_neg,1);
1518+
if (drelPt_neg)
1519+
registry.fill(HIST("neg/meanPT/meanRelPtA"), eta, spm.centrality, dptV1A_neg / drelPt_neg, 1);
1520+
if (drelPt_neg)
1521+
registry.fill(HIST("neg/meanPT/meanRelPtC"), eta, spm.centrality, dptV1C_neg / drelPt_neg, 1);
15221522
}
15231523
}
15241524

1525-
1526-
1527-
delete hRelEtaPt;
1528-
delete ptV1A;
1529-
delete ptV1C;
1530-
1531-
1525+
delete hRelEtaPt;
1526+
delete ptV1A;
1527+
delete ptV1C;
15321528
}
15331529

15341530
PROCESS_SWITCH(FlowSP, processData, "Process analysis for non-derived data", true);

0 commit comments

Comments
 (0)