@@ -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 ();
@@ -1304,9 +1302,9 @@ struct FlowSP {
13041302
13051303 fillEventQA<kAfter >(collision, tracks);
13061304
1307- TProfile2D* hRelEtaPt = new TProfile2D (" hRelEtaPt" , " hRelEtaPt" , 8 ,-.8 ,.8 , 3 ,0 , 3 );
1308- TProfile2D* ptV1A = new TProfile2D (" ptV1A" , " ptV1A" , 8 ,-.8 ,.8 , 3 ,0 , 3 );
1309- TProfile2D* ptV1C = new TProfile2D (" ptV1C" , " ptV1C" , 8 ,-.8 ,.8 , 3 ,0 , 3 );
1305+ TProfile2D* hRelEtaPt = new TProfile2D (" hRelEtaPt" , " hRelEtaPt" , 8 , -.8 , .8 , 3 , 0 , 3 );
1306+ TProfile2D* ptV1A = new TProfile2D (" ptV1A" , " ptV1A" , 8 , -.8 , .8 , 3 , 0 , 3 );
1307+ TProfile2D* ptV1C = new TProfile2D (" ptV1C" , " ptV1C" , 8 , -.8 , .8 , 3 , 0 , 3 );
13101308
13111309 for (const auto & track : tracks) {
13121310
@@ -1382,7 +1380,6 @@ struct FlowSP {
13821380 ptV1C->Fill (track.eta (), kInclusive , track.pt () * ((spm.uy * spm.qyC + spm.ux * spm.qxC ) / std::sqrt (std::fabs (spm.corrQQ ))), weight);
13831381 ptV1C->Fill (track.eta (), spm.charge , track.pt () * ((spm.uy * spm.qyC + spm.ux * spm.qxC ) / std::sqrt (std::fabs (spm.corrQQ ))), weight_charged);
13841382
1385-
13861383 fillAllQA<kAfter , kUnidentified >(track);
13871384 if (cfgFillPIDQA) {
13881385 switch (trackPID) {
@@ -1481,10 +1478,9 @@ struct FlowSP {
14811478
14821479 } // end of track loop
14831480
1484-
1485- // Now we want to fill the final relPt histogram
1481+ // Now we want to fill the final relPt histogram
14861482 // Loop over all eta and fill bins
1487- if (cfgFillMeanPT) {
1483+ if (cfgFillMeanPT) {
14881484 for (int i = 0 ; i < hRelEtaPt->GetNbinsX (); i++) {
14891485 double eta = hRelEtaPt->GetXaxis ()->GetBinCenter (i);
14901486
@@ -1493,32 +1489,34 @@ struct FlowSP {
14931489 double drelPt = hRelEtaPt->GetBinContent (bin);
14941490 double dptV1A = ptV1A->GetBinContent (bin);
14951491 double dptV1C = ptV1C->GetBinContent (bin);
1496- if (drelPt) registry.fill (HIST (" incl/meanPT/meanRelPtA" ), eta, spm.centrality , dptV1A / drelPt,1 );
1497- if (drelPt) registry.fill (HIST (" incl/meanPT/meanRelPtC" ), eta, spm.centrality , dptV1C / drelPt,1 );
1492+ if (drelPt)
1493+ registry.fill (HIST (" incl/meanPT/meanRelPtA" ), eta, spm.centrality , dptV1A / drelPt, 1 );
1494+ if (drelPt)
1495+ registry.fill (HIST (" incl/meanPT/meanRelPtC" ), eta, spm.centrality , dptV1C / drelPt, 1 );
14981496
14991497 bin = hRelEtaPt->FindBin (eta, kPositive );
15001498 double drelPt_pos = hRelEtaPt->GetBinContent (bin);
15011499 double dptV1A_pos = ptV1A->GetBinContent (bin);
15021500 double dptV1C_pos = ptV1C->GetBinContent (bin);
1503- if (drelPt_pos) registry.fill (HIST (" pos/meanPT/meanRelPtA" ), eta, spm.centrality , dptV1A_pos / drelPt_pos,1 );
1504- if (drelPt_pos) registry.fill (HIST (" pos/meanPT/meanRelPtC" ), eta, spm.centrality , dptV1C_pos / drelPt_pos,1 );
1501+ if (drelPt_pos)
1502+ registry.fill (HIST (" pos/meanPT/meanRelPtA" ), eta, spm.centrality , dptV1A_pos / drelPt_pos, 1 );
1503+ if (drelPt_pos)
1504+ registry.fill (HIST (" pos/meanPT/meanRelPtC" ), eta, spm.centrality , dptV1C_pos / drelPt_pos, 1 );
15051505
15061506 bin = hRelEtaPt->FindBin (eta, kNegative );
15071507 double drelPt_neg = hRelEtaPt->GetBinContent (bin);
15081508 double dptV1A_neg = ptV1A->GetBinContent (bin);
15091509 double dptV1C_neg = ptV1C->GetBinContent (bin);
1510- if (drelPt_neg) registry.fill (HIST (" neg/meanPT/meanRelPtA" ), eta, spm.centrality , dptV1A_neg / drelPt_neg,1 );
1511- if (drelPt_neg) registry.fill (HIST (" neg/meanPT/meanRelPtC" ), eta, spm.centrality , dptV1C_neg / drelPt_neg,1 );
1510+ if (drelPt_neg)
1511+ registry.fill (HIST (" neg/meanPT/meanRelPtA" ), eta, spm.centrality , dptV1A_neg / drelPt_neg, 1 );
1512+ if (drelPt_neg)
1513+ registry.fill (HIST (" neg/meanPT/meanRelPtC" ), eta, spm.centrality , dptV1C_neg / drelPt_neg, 1 );
15121514 }
15131515 }
15141516
1515-
1516-
1517- delete hRelEtaPt;
1518- delete ptV1A;
1519- delete ptV1C;
1520-
1521-
1517+ delete hRelEtaPt;
1518+ delete ptV1A;
1519+ delete ptV1C;
15221520 }
15231521
15241522 PROCESS_SWITCH (FlowSP, processData, " Process analysis for non-derived data" , true );
0 commit comments