@@ -49,7 +49,7 @@ using namespace o2::framework::expressions;
4949struct FlowSP {
5050 // QA Plots
5151 O2_DEFINE_CONFIGURABLE (cfgFillQAHistos, bool , true , " Fill histograms for event and track QA" );
52- // Flags to make and fill histograms
52+ // Flags to make and fill histograms
5353 O2_DEFINE_CONFIGURABLE (cfgFillMixedHarmonics, bool , true , " Flag to make and fill histos for mixed harmonics" );
5454 O2_DEFINE_CONFIGURABLE (cfgFillEventPlane, bool , true , " Flag to make and fill histos with Event Plane" );
5555 O2_DEFINE_CONFIGURABLE (cfgFillXandYterms, bool , true , " Flag to make and fill histos for with separate x and y terms for SPM" );
@@ -246,7 +246,7 @@ struct FlowSP {
246246 }
247247
248248 if (doprocessData) {
249-
249+
250250 // track properties per centrality and per eta, pt bin
251251 registry.add <TProfile>(" incl/vnC_eta" , " " , kTProfile , {axisEtaVn});
252252 registry.add <TProfile>(" incl/vnA_eta" , " " , kTProfile , {axisEtaVn});
@@ -261,7 +261,7 @@ struct FlowSP {
261261 registry.add <TProfile>(" incl/vnC_cent_plusEta" , " " , kTProfile , {axisCent});
262262 registry.add <TProfile>(" incl/vnA_cent_plusEta" , " " , kTProfile , {axisCent});
263263
264- if (cfgFillXandYterms){
264+ if (cfgFillXandYterms) {
265265 registry.add <TProfile>(" incl/vnAx_eta" , " " , kTProfile , {axisEtaVn});
266266 registry.add <TProfile>(" incl/vnAy_eta" , " " , kTProfile , {axisEtaVn});
267267 registry.add <TProfile>(" incl/vnCx_eta" , " " , kTProfile , {axisEtaVn});
@@ -276,7 +276,7 @@ struct FlowSP {
276276 registry.add <TProfile>(" incl/vnAy_pt_odd" , " " , kTProfile , {axisPt});
277277 }
278278
279- if (cfgFillMixedHarmonics){
279+ if (cfgFillMixedHarmonics) {
280280 registry.add <TProfile>(" incl/vnAxCxUx_pt_MH" , " " , kTProfile , {axisPt});
281281 registry.add <TProfile>(" incl/vnAxCyUx_pt_MH" , " " , kTProfile , {axisPt});
282282 registry.add <TProfile>(" incl/vnAxCyUy_pt_MH" , " " , kTProfile , {axisPt});
@@ -293,7 +293,7 @@ struct FlowSP {
293293 registry.add <TProfile>(" incl/vnAyCxUy_eta_MH" , " " , kTProfile , {axisEtaVn});
294294 }
295295
296- if (cfgFillEventPlane){
296+ if (cfgFillEventPlane) {
297297 registry.add <TProfile>(" incl/vnA_cent_EP" , " " , kTProfile , {axisCent});
298298 registry.add <TProfile>(" incl/vnC_cent_EP" , " " , kTProfile , {axisCent});
299299 registry.add <TProfile>(" incl/vnFull_cent_EP" , " " , kTProfile , {axisCent});
@@ -339,7 +339,7 @@ struct FlowSP {
339339 registry.add (" QA/after/PsiC_vs_Vz" , " " , {HistType::kTH2D , {axisPhiPlane, axisVz}});
340340 registry.add (" QA/after/PsiFull_vs_Vz" , " " , {HistType::kTH2D , {axisPhiPlane, axisVz}});
341341
342- registry.addClone (" QA/after/" , " QA/before/" );
342+ registry.addClone (" QA/after/" , " QA/before/" );
343343 // track QA for pos, neg, incl
344344 registry.add <TH1>(" incl/QA/hPt" , " " , kTH1D , {axisPt});
345345 registry.add <TH1>(" incl/QA/hPhi" , " " , kTH1D , {axisPhi});
@@ -350,11 +350,10 @@ struct FlowSP {
350350 registry.add <TH2>(" incl/QA/hDCAz_pt" , " " , kTH2D , {axisPt, axisDCAz});
351351 registry.add (" incl/QA/hSharedClusters_pt" , " " , {HistType::kTH2D , {axisPt, axisShCl}});
352352 registry.add (" incl/QA/hCrossedRows_pt" , " " , {HistType::kTH2D , {axisPt, axisCl}});
353-
354353 }
355354 }
356-
357- if (cfgFillChargeDependence){
355+
356+ if (cfgFillChargeDependence) {
358357 registry.addClone (" incl/" , " pos/" );
359358 registry.addClone (" incl/" , " neg/" );
360359 }
@@ -421,7 +420,6 @@ struct FlowSP {
421420 fPhiCutLow = new TF1 (" fPhiCutLow" , " 0.06/x+pi/18.0-0.06" , 0 , 100 );
422421 fPhiCutHigh = new TF1 (" fPhiCutHigh" , " 0.1/x+pi/18.0+0.06" , 0 , 100 );
423422 }
424-
425423 }
426424
427425 int getMagneticField (uint64_t timestamp)
@@ -680,7 +678,7 @@ struct FlowSP {
680678 registry.fill (HIST (Charge[ct]) + HIST (" vnA_pt" ), track.pt (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)), wacc * weff);
681679 registry.fill (HIST (Charge[ct]) + HIST (" vnC_pt" ), track.pt (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)), wacc * weff);
682680
683- if (cfgFillMixedHarmonics){
681+ if (cfgFillMixedHarmonics) {
684682 registry.fill (HIST (Charge[ct]) + HIST (" vnAxCxUx_eta_MH" ), track.eta (), (uxMH * qxA * qxC) / corrQQx, wacc * weff);
685683 registry.fill (HIST (Charge[ct]) + HIST (" vnAxCyUx_eta_MH" ), track.eta (), (uxMH * qyA * qyC) / corrQQy, wacc * weff);
686684 registry.fill (HIST (Charge[ct]) + HIST (" vnAxCyUy_eta_MH" ), track.eta (), (uyMH * qxA * qyC) / corrQQx, wacc * weff);
@@ -697,7 +695,7 @@ struct FlowSP {
697695 registry.fill (HIST (Charge[ct]) + HIST (" vnAyCxUy_cent_MH" ), centrality, (uyMH * qyA * qxC) / corrQQy, wacc * weff);
698696 }
699697
700- if (cfgFillXandYterms){
698+ if (cfgFillXandYterms) {
701699 registry.fill (HIST (Charge[ct]) + HIST (" vnAx_eta" ), track.eta (), (ux * qxA) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
702700 registry.fill (HIST (Charge[ct]) + HIST (" vnAy_eta" ), track.eta (), (uy * qyA) / std::sqrt (std::fabs (corrQQy)), wacc * weff);
703701 registry.fill (HIST (Charge[ct]) + HIST (" vnCx_eta" ), track.eta (), (ux * qxC) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
@@ -707,10 +705,9 @@ struct FlowSP {
707705 registry.fill (HIST (Charge[ct]) + HIST (" vnAy_pt" ), track.pt (), (uy * qyA) / std::sqrt (std::fabs (corrQQy)), wacc * weff);
708706 registry.fill (HIST (Charge[ct]) + HIST (" vnCx_pt" ), track.pt (), (ux * qxC) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
709707 registry.fill (HIST (Charge[ct]) + HIST (" vnCy_pt" ), track.pt (), (uy * qyC) / std::sqrt (std::fabs (corrQQy)), wacc * weff);
710-
711708 }
712709
713- if (cfgFillEventPlane){
710+ if (cfgFillEventPlane) {
714711 registry.fill (HIST (Charge[ct]) + HIST (" vnA_eta_EP" ), track.eta (), vnA, wacc * weff);
715712 registry.fill (HIST (Charge[ct]) + HIST (" vnC_eta_EP" ), track.eta (), vnC, wacc * weff);
716713 registry.fill (HIST (Charge[ct]) + HIST (" vnFull_eta_EP" ), track.eta (), vnFull, wacc * weff);
@@ -733,7 +730,7 @@ struct FlowSP {
733730 registry.fill (HIST (Charge[ct]) + HIST (" vnA_pt_odd" ), track.pt (), -1.0 * (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)), wacc * weff);
734731 registry.fill (HIST (Charge[ct]) + HIST (" vnC_pt_odd" ), track.pt (), -1.0 * (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)), wacc * weff);
735732
736- if (cfgFillXandYterms){
733+ if (cfgFillXandYterms) {
737734 registry.fill (HIST (Charge[ct]) + HIST (" vnAx_pt_odd" ), track.pt (), -1.0 * (ux * qxA) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
738735 registry.fill (HIST (Charge[ct]) + HIST (" vnAy_pt_odd" ), track.pt (), -1.0 * (uy * qyA) / std::sqrt (std::fabs (corrQQy)), wacc * weff);
739736 registry.fill (HIST (Charge[ct]) + HIST (" vnCx_pt_odd" ), track.pt (), -1.0 * (ux * qxC) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
@@ -746,7 +743,7 @@ struct FlowSP {
746743 registry.fill (HIST (Charge[ct]) + HIST (" vnA_pt_odd" ), track.pt (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)), wacc * weff);
747744 registry.fill (HIST (Charge[ct]) + HIST (" vnC_pt_odd" ), track.pt (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)), wacc * weff);
748745
749- if (cfgFillXandYterms){
746+ if (cfgFillXandYterms) {
750747 registry.fill (HIST (Charge[ct]) + HIST (" vnAx_pt_odd" ), track.pt (), (ux * qxA) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
751748 registry.fill (HIST (Charge[ct]) + HIST (" vnAy_pt_odd" ), track.pt (), (uy * qyA) / std::sqrt (std::fabs (corrQQy)), wacc * weff);
752749 registry.fill (HIST (Charge[ct]) + HIST (" vnCx_pt_odd" ), track.pt (), (ux * qxC) / std::sqrt (std::fabs (corrQQx)), wacc * weff);
@@ -815,7 +812,7 @@ struct FlowSP {
815812 // https://twiki.cern.ch/twiki/pub/ALICE/DirectedFlowAnalysisNote/vn_ZDC_ALICE_INT_NOTE_version02.pdf
816813 double psiFull = 1.0 * std::atan2 (qyA + qyC, qxA + qxC);
817814
818- if (cfgFillQAHistos){
815+ if (cfgFillQAHistos) {
819816 fillEventQA<kAfter >(collision, tracks);
820817
821818 registry.fill (HIST (" QA/hSPplaneA" ), psiA, 1 );
@@ -924,10 +921,10 @@ struct FlowSP {
924921
925922 fillHistograms<kInclusive >(track, wacc, weff, ux, uy, uxMH, uyMH, qxA, qyA, qxC, qyC, corrQQx, corrQQy, corrQQ, vnA, vnC, vnFull, centrality);
926923
927- if (cfgFillQAHistos)
924+ if (cfgFillQAHistos)
928925 fillTrackQA<kInclusive >(track, vtxz, wacc, weff);
929926
930- if (cfgFillChargeDependence){
927+ if (cfgFillChargeDependence) {
931928 if (pos) {
932929 fillHistograms<kPositive >(track, waccP, weffP, ux, uy, uxMH, uyMH, qxA, qyA, qxC, qyC, corrQQx, corrQQy, corrQQ, vnA, vnC, vnFull, centrality);
933930 fillTrackQA<kPositive >(track, vtxz, wacc, weff);
0 commit comments