@@ -347,7 +347,7 @@ struct FlowSP {
347347 {
348348
349349 if (track.tpcNClsFound () < cfgNcls)
350- return false ;
350+ return false ;
351351
352352 double phimodn = track.phi ();
353353 if (field < 0 ) // for negative polarity field
@@ -509,10 +509,10 @@ struct FlowSP {
509509 registry.fill (HIST (" incl/vnA_eta" ), track.eta (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
510510 registry.fill (HIST (" incl/vnC_eta" ), track.eta (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
511511
512- registry.fill (HIST (" incl/vnAxCxUx_eta_MH" ), track.eta (), (ux * qxA * qxC)/ corrQQx);
513- registry.fill (HIST (" incl/vnAxCyUx_eta_MH" ), track.eta (), (ux * qyA * qyC)/ corrQQy);
514- registry.fill (HIST (" incl/vnAxCyUy_eta_MH" ), track.eta (), (uy * qxA * qyC)/ corrQQx);
515- registry.fill (HIST (" incl/vnAyCxUy_eta_MH" ), track.eta (), (uy * qyA * qxC)/ corrQQy);
512+ registry.fill (HIST (" incl/vnAxCxUx_eta_MH" ), track.eta (), (ux * qxA * qxC) / corrQQx);
513+ registry.fill (HIST (" incl/vnAxCyUx_eta_MH" ), track.eta (), (ux * qyA * qyC) / corrQQy);
514+ registry.fill (HIST (" incl/vnAxCyUy_eta_MH" ), track.eta (), (uy * qxA * qyC) / corrQQx);
515+ registry.fill (HIST (" incl/vnAyCxUy_eta_MH" ), track.eta (), (uy * qyA * qxC) / corrQQy);
516516
517517 registry.fill (HIST (" incl/vnAx_pt" ), track.pt (), (ux * qxA) / std::sqrt (std::fabs (corrQQx)));
518518 registry.fill (HIST (" incl/vnAy_pt" ), track.pt (), (uy * qyA) / std::sqrt (std::fabs (corrQQy)));
@@ -521,12 +521,11 @@ struct FlowSP {
521521 registry.fill (HIST (" incl/vnA_pt" ), track.pt (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
522522 registry.fill (HIST (" incl/vnC_pt" ), track.pt (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
523523
524- registry.fill (HIST (" incl/vnAxCxUx_pt_MH" ), track.pt (), (ux * qxA * qxC)/ corrQQx);
525- registry.fill (HIST (" incl/vnAxCyUx_pt_MH" ), track.pt (), (ux * qyA * qyC)/ corrQQy);
526- registry.fill (HIST (" incl/vnAxCyUy_pt_MH" ), track.pt (), (uy * qxA * qyC)/ corrQQx);
527- registry.fill (HIST (" incl/vnAyCxUy_pt_MH" ), track.pt (), (uy * qyA * qxC)/ corrQQy);
524+ registry.fill (HIST (" incl/vnAxCxUx_pt_MH" ), track.pt (), (ux * qxA * qxC) / corrQQx);
525+ registry.fill (HIST (" incl/vnAxCyUx_pt_MH" ), track.pt (), (ux * qyA * qyC) / corrQQy);
526+ registry.fill (HIST (" incl/vnAxCyUy_pt_MH" ), track.pt (), (uy * qxA * qyC) / corrQQx);
527+ registry.fill (HIST (" incl/vnAyCxUy_pt_MH" ), track.pt (), (uy * qyA * qxC) / corrQQy);
528528
529-
530529 registry.fill (HIST (" incl/vnA_eta_EP" ), track.eta (), vnA);
531530 registry.fill (HIST (" incl/vnC_eta_EP" ), track.eta (), vnC);
532531 registry.fill (HIST (" incl/vnFull_eta_EP" ), track.eta (), vnFull);
@@ -535,9 +534,9 @@ struct FlowSP {
535534 registry.fill (HIST (" incl/vnC_pt_EP" ), track.pt (), vnC);
536535 registry.fill (HIST (" incl/vnFull_pt_EP" ), track.pt (), vnFull);
537536
538- // For integrated v1 take only tracks from eta>0.
537+ // For integrated v1 take only tracks from eta>0.
539538 // Following https://arxiv.org/pdf/1306.4145
540- if (cfgHarm== 1 && track.eta ()< 0 )
539+ if (cfgHarm == 1 && track.eta () < 0 )
541540 continue ;
542541
543542 registry.fill (HIST (" incl/vnAx_cent" ), centrality, (ux * qxA) / std::sqrt (std::fabs (corrQQx)));
@@ -547,10 +546,10 @@ struct FlowSP {
547546 registry.fill (HIST (" incl/vnA_cent" ), centrality, (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
548547 registry.fill (HIST (" incl/vnC_cent" ), centrality, (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
549548
550- registry.fill (HIST (" incl/vnAxCxUx_cent_MH" ), centrality, (ux * qxA * qxC)/ corrQQx);
551- registry.fill (HIST (" incl/vnAxCyUx_cent_MH" ), centrality, (ux * qyA * qyC)/ corrQQy);
552- registry.fill (HIST (" incl/vnAxCyUy_cent_MH" ), centrality, (uy * qxA * qyC)/ corrQQx);
553- registry.fill (HIST (" incl/vnAyCxUy_cent_MH" ), centrality, (uy * qyA * qxC)/ corrQQy);
549+ registry.fill (HIST (" incl/vnAxCxUx_cent_MH" ), centrality, (ux * qxA * qxC) / corrQQx);
550+ registry.fill (HIST (" incl/vnAxCyUx_cent_MH" ), centrality, (ux * qyA * qyC) / corrQQy);
551+ registry.fill (HIST (" incl/vnAxCyUy_cent_MH" ), centrality, (uy * qxA * qyC) / corrQQx);
552+ registry.fill (HIST (" incl/vnAyCxUy_cent_MH" ), centrality, (uy * qyA * qxC) / corrQQy);
554553
555554 registry.fill (HIST (" incl/vnA_cent_EP" ), centrality, vnA);
556555 registry.fill (HIST (" incl/vnC_cent_EP" ), centrality, vnC);
@@ -564,10 +563,10 @@ struct FlowSP {
564563 registry.fill (HIST (" pos/vnA_eta" ), track.eta (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
565564 registry.fill (HIST (" pos/vnC_eta" ), track.eta (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
566565
567- registry.fill (HIST (" pos/vnAxCxUx_eta_MH" ), track.eta (), (ux * qxA * qxC)/ corrQQx);
568- registry.fill (HIST (" pos/vnAxCyUx_eta_MH" ), track.eta (), (ux * qyA * qyC)/ corrQQy);
569- registry.fill (HIST (" pos/vnAxCyUy_eta_MH" ), track.eta (), (uy * qxA * qyC)/ corrQQx);
570- registry.fill (HIST (" pos/vnAyCxUy_eta_MH" ), track.eta (), (uy * qyA * qxC)/ corrQQy);
566+ registry.fill (HIST (" pos/vnAxCxUx_eta_MH" ), track.eta (), (ux * qxA * qxC) / corrQQx);
567+ registry.fill (HIST (" pos/vnAxCyUx_eta_MH" ), track.eta (), (ux * qyA * qyC) / corrQQy);
568+ registry.fill (HIST (" pos/vnAxCyUy_eta_MH" ), track.eta (), (uy * qxA * qyC) / corrQQx);
569+ registry.fill (HIST (" pos/vnAyCxUy_eta_MH" ), track.eta (), (uy * qyA * qxC) / corrQQy);
571570
572571 registry.fill (HIST (" pos/vnAx_pt" ), track.pt (), (ux * qxA) / std::sqrt (std::fabs (corrQQx)));
573572 registry.fill (HIST (" pos/vnAy_pt" ), track.pt (), (uy * qyA) / std::sqrt (std::fabs (corrQQy)));
@@ -576,12 +575,11 @@ struct FlowSP {
576575 registry.fill (HIST (" pos/vnA_pt" ), track.pt (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
577576 registry.fill (HIST (" pos/vnC_pt" ), track.pt (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
578577
579- registry.fill (HIST (" pos/vnAxCxUx_pt_MH" ), track.pt (), (ux * qxA * qxC)/ corrQQx);
580- registry.fill (HIST (" pos/vnAxCyUx_pt_MH" ), track.pt (), (ux * qyA * qyC)/ corrQQy);
581- registry.fill (HIST (" pos/vnAxCyUy_pt_MH" ), track.pt (), (uy * qxA * qyC)/ corrQQx);
582- registry.fill (HIST (" pos/vnAyCxUy_pt_MH" ), track.pt (), (uy * qyA * qxC)/ corrQQy);
578+ registry.fill (HIST (" pos/vnAxCxUx_pt_MH" ), track.pt (), (ux * qxA * qxC) / corrQQx);
579+ registry.fill (HIST (" pos/vnAxCyUx_pt_MH" ), track.pt (), (ux * qyA * qyC) / corrQQy);
580+ registry.fill (HIST (" pos/vnAxCyUy_pt_MH" ), track.pt (), (uy * qxA * qyC) / corrQQx);
581+ registry.fill (HIST (" pos/vnAyCxUy_pt_MH" ), track.pt (), (uy * qyA * qxC) / corrQQy);
583582
584-
585583 registry.fill (HIST (" pos/vnA_eta_EP" ), track.eta (), vnA);
586584 registry.fill (HIST (" pos/vnC_eta_EP" ), track.eta (), vnC);
587585 registry.fill (HIST (" pos/vnFull_eta_EP" ), track.eta (), vnFull);
@@ -590,9 +588,9 @@ struct FlowSP {
590588 registry.fill (HIST (" pos/vnC_pt_EP" ), track.pt (), vnC);
591589 registry.fill (HIST (" pos/vnFull_pt_EP" ), track.pt (), vnFull);
592590
593- // For integrated v1 take only tracks from eta>0.
591+ // For integrated v1 take only tracks from eta>0.
594592 // Following https://arxiv.org/pdf/1306.4145
595- if (cfgHarm== 1 && track.eta ()< 0 )
593+ if (cfgHarm == 1 && track.eta () < 0 )
596594 continue ;
597595
598596 registry.fill (HIST (" pos/vnAx_cent" ), centrality, (ux * qxA) / std::sqrt (std::fabs (corrQQx)));
@@ -602,10 +600,10 @@ struct FlowSP {
602600 registry.fill (HIST (" pos/vnA_cent" ), centrality, (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
603601 registry.fill (HIST (" pos/vnC_cent" ), centrality, (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
604602
605- registry.fill (HIST (" pos/vnAxCxUx_cent_MH" ), centrality, (ux * qxA * qxC)/ corrQQx);
606- registry.fill (HIST (" pos/vnAxCyUx_cent_MH" ), centrality, (ux * qyA * qyC)/ corrQQy);
607- registry.fill (HIST (" pos/vnAxCyUy_cent_MH" ), centrality, (uy * qxA * qyC)/ corrQQx);
608- registry.fill (HIST (" pos/vnAyCxUy_cent_MH" ), centrality, (uy * qyA * qxC)/ corrQQy);
603+ registry.fill (HIST (" pos/vnAxCxUx_cent_MH" ), centrality, (ux * qxA * qxC) / corrQQx);
604+ registry.fill (HIST (" pos/vnAxCyUx_cent_MH" ), centrality, (ux * qyA * qyC) / corrQQy);
605+ registry.fill (HIST (" pos/vnAxCyUy_cent_MH" ), centrality, (uy * qxA * qyC) / corrQQx);
606+ registry.fill (HIST (" pos/vnAyCxUy_cent_MH" ), centrality, (uy * qyA * qxC) / corrQQy);
609607
610608 registry.fill (HIST (" pos/vnA_cent_EP" ), centrality, vnA);
611609 registry.fill (HIST (" pos/vnC_cent_EP" ), centrality, vnC);
@@ -618,10 +616,10 @@ struct FlowSP {
618616 registry.fill (HIST (" neg/vnA_eta" ), track.eta (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
619617 registry.fill (HIST (" neg/vnC_eta" ), track.eta (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
620618
621- registry.fill (HIST (" neg/vnAxCxUx_eta_MH" ), track.eta (), (ux * qxA * qxC)/ corrQQx);
622- registry.fill (HIST (" neg/vnAxCyUx_eta_MH" ), track.eta (), (ux * qyA * qyC)/ corrQQy);
623- registry.fill (HIST (" neg/vnAxCyUy_eta_MH" ), track.eta (), (uy * qxA * qyC)/ corrQQx);
624- registry.fill (HIST (" neg/vnAyCxUy_eta_MH" ), track.eta (), (uy * qyA * qxC)/ corrQQy);
619+ registry.fill (HIST (" neg/vnAxCxUx_eta_MH" ), track.eta (), (ux * qxA * qxC) / corrQQx);
620+ registry.fill (HIST (" neg/vnAxCyUx_eta_MH" ), track.eta (), (ux * qyA * qyC) / corrQQy);
621+ registry.fill (HIST (" neg/vnAxCyUy_eta_MH" ), track.eta (), (uy * qxA * qyC) / corrQQx);
622+ registry.fill (HIST (" neg/vnAyCxUy_eta_MH" ), track.eta (), (uy * qyA * qxC) / corrQQy);
625623
626624 registry.fill (HIST (" neg/vnAx_pt" ), track.pt (), (ux * qxA) / std::sqrt (std::fabs (corrQQx)));
627625 registry.fill (HIST (" neg/vnAy_pt" ), track.pt (), (uy * qyA) / std::sqrt (std::fabs (corrQQy)));
@@ -630,12 +628,11 @@ struct FlowSP {
630628 registry.fill (HIST (" neg/vnA_pt" ), track.pt (), (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
631629 registry.fill (HIST (" neg/vnC_pt" ), track.pt (), (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
632630
633- registry.fill (HIST (" neg/vnAxCxUx_pt_MH" ), track.pt (), (ux * qxA * qxC)/ corrQQx);
634- registry.fill (HIST (" neg/vnAxCyUx_pt_MH" ), track.pt (), (ux * qyA * qyC)/ corrQQy);
635- registry.fill (HIST (" neg/vnAxCyUy_pt_MH" ), track.pt (), (uy * qxA * qyC)/ corrQQx);
636- registry.fill (HIST (" neg/vnAyCxUy_pt_MH" ), track.pt (), (uy * qyA * qxC)/ corrQQy);
631+ registry.fill (HIST (" neg/vnAxCxUx_pt_MH" ), track.pt (), (ux * qxA * qxC) / corrQQx);
632+ registry.fill (HIST (" neg/vnAxCyUx_pt_MH" ), track.pt (), (ux * qyA * qyC) / corrQQy);
633+ registry.fill (HIST (" neg/vnAxCyUy_pt_MH" ), track.pt (), (uy * qxA * qyC) / corrQQx);
634+ registry.fill (HIST (" neg/vnAyCxUy_pt_MH" ), track.pt (), (uy * qyA * qxC) / corrQQy);
637635
638-
639636 registry.fill (HIST (" neg/vnA_eta_EP" ), track.eta (), vnA);
640637 registry.fill (HIST (" neg/vnC_eta_EP" ), track.eta (), vnC);
641638 registry.fill (HIST (" neg/vnFull_eta_EP" ), track.eta (), vnFull);
@@ -644,9 +641,9 @@ struct FlowSP {
644641 registry.fill (HIST (" neg/vnC_pt_EP" ), track.pt (), vnC);
645642 registry.fill (HIST (" neg/vnFull_pt_EP" ), track.pt (), vnFull);
646643
647- // For integrated v1 take only tracks from eta>0.
644+ // For integrated v1 take only tracks from eta>0.
648645 // Following https://arxiv.org/pdf/1306.4145
649- if (cfgHarm== 1 && track.eta ()< 0 )
646+ if (cfgHarm == 1 && track.eta () < 0 )
650647 continue ;
651648
652649 registry.fill (HIST (" neg/vnAx_cent" ), centrality, (ux * qxA) / std::sqrt (std::fabs (corrQQx)));
@@ -656,10 +653,10 @@ struct FlowSP {
656653 registry.fill (HIST (" neg/vnA_cent" ), centrality, (uy * qyA + ux * qxA) / std::sqrt (std::fabs (corrQQ)));
657654 registry.fill (HIST (" neg/vnC_cent" ), centrality, (uy * qyC + ux * qxC) / std::sqrt (std::fabs (corrQQ)));
658655
659- registry.fill (HIST (" neg/vnAxCxUx_cent_MH" ), centrality, (ux * qxA * qxC)/ corrQQx);
660- registry.fill (HIST (" neg/vnAxCyUx_cent_MH" ), centrality, (ux * qyA * qyC)/ corrQQy);
661- registry.fill (HIST (" neg/vnAxCyUy_cent_MH" ), centrality, (uy * qxA * qyC)/ corrQQx);
662- registry.fill (HIST (" neg/vnAyCxUy_cent_MH" ), centrality, (uy * qyA * qxC)/ corrQQy);
656+ registry.fill (HIST (" neg/vnAxCxUx_cent_MH" ), centrality, (ux * qxA * qxC) / corrQQx);
657+ registry.fill (HIST (" neg/vnAxCyUx_cent_MH" ), centrality, (ux * qyA * qyC) / corrQQy);
658+ registry.fill (HIST (" neg/vnAxCyUy_cent_MH" ), centrality, (uy * qxA * qyC) / corrQQx);
659+ registry.fill (HIST (" neg/vnAyCxUy_cent_MH" ), centrality, (uy * qyA * qxC) / corrQQy);
663660
664661 registry.fill (HIST (" neg/vnA_cent_EP" ), centrality, vnA);
665662 registry.fill (HIST (" neg/vnC_cent_EP" ), centrality, vnC);
0 commit comments