@@ -460,28 +460,27 @@ struct FlowSP {
460460 fMultPVCutLow = new TF1 (" fMultPVCutLow" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
461461
462462 // Variables from fitting distribution with mean+1sigma
463- double fitParamLowPV1 = 2942.55 ;
463+ double fitParamLowPV1 = 2942.55 ;
464464 double fitParamLowPV2 = -103.111 ;
465465 double fitParamLowPV3 = 1.4397 ;
466466 double fitParamLowPV4 = -0.00974862 ;
467467 double fitParamLowPV5 = 2.71433e-05 ;
468468
469- if (cfgnSigmaMultCuts == 2 ){
469+ if (cfgnSigmaMultCuts == 2 ) {
470470 fitParamLowPV1 = 2665.68 ;
471471 fitParamLowPV2 = -93.3784 ;
472472 fitParamLowPV3 = 1.27137 ;
473473 fitParamLowPV4 = -0.00818936 ;
474474 fitParamLowPV5 = 2.115e-05 ;
475- }
476- else if (cfgnSigmaMultCuts == 3 ){
475+ } else if (cfgnSigmaMultCuts == 3 ) {
477476 fitParamLowPV1 = 2389.99 ;
478477 fitParamLowPV2 = -83.8483 ;
479478 fitParamLowPV3 = 1.11062 ;
480479 fitParamLowPV4 = -0.00672263 ;
481480 fitParamLowPV5 = 1.54725e-05 ;
482481 }
483482
484- fMultPVCutLow ->SetParameters ( fitParamLowPV1, fitParamLowPV2, fitParamLowPV3, fitParamLowPV4, fitParamLowPV5);
483+ fMultPVCutLow ->SetParameters (fitParamLowPV1, fitParamLowPV2, fitParamLowPV3, fitParamLowPV4, fitParamLowPV5);
485484
486485 fMultPVCutHigh = new TF1 (" fMultPVCutHigh" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
487486
@@ -492,14 +491,13 @@ struct FlowSP {
492491 double fitParamHighPV4 = -0.0145343 ;
493492 double fitParamHighPV5 = 4.80688e-05 ;
494493
495- if (cfgnSigmaMultCuts == 2 ){
494+ if (cfgnSigmaMultCuts == 2 ) {
496495 fitParamHighPV1 = 3787.93 ;
497496 fitParamHighPV2 = -135.184 ;
498497 fitParamHighPV3 = 2.07683 ;
499498 fitParamHighPV4 = -0.0165997 ;
500499 fitParamHighPV5 = 5.68725e-05 ;
501- }
502- else if (cfgnSigmaMultCuts == 3 ){
500+ } else if (cfgnSigmaMultCuts == 3 ) {
503501 fitParamHighPV1 = 4067.4 ;
504502 fitParamHighPV2 = -145.485 ;
505503 fitParamHighPV3 = 2.27273 ;
@@ -516,14 +514,13 @@ struct FlowSP {
516514 double fitParamLow4 = -0.00235284 ;
517515 double fitParamLow5 = 3.01132e-06 ;
518516
519- if (cfgnSigmaMultCuts == 2 ){
517+ if (cfgnSigmaMultCuts == 2 ) {
520518 fitParamLow1 = 1307.92 ;
521519 fitParamLow2 = -39.9168 ;
522520 fitParamLow3 = 0.412675 ;
523521 fitParamLow4 = -0.00148081 ;
524522 fitParamLow5 = 1.10868e-07 ;
525- }
526- else if (cfgnSigmaMultCuts == 3 ){
523+ } else if (cfgnSigmaMultCuts == 3 ) {
527524 fitParamLow1 = 1048.48 ;
528525 fitParamLow2 = -31.4568 ;
529526 fitParamLow3 = 0.287794 ;
@@ -540,14 +537,13 @@ struct FlowSP {
540537 double fitParamHigh4 = -0.00496563 ;
541538 double fitParamHigh5 = 1.34314e-05 ;
542539
543- if (cfgnSigmaMultCuts == 2 ){
540+ if (cfgnSigmaMultCuts == 2 ) {
544541 fitParamHigh1 = 2350.39 ;
545542 fitParamHigh2 = -74.6939 ;
546543 fitParamHigh3 = 0.953287 ;
547544 fitParamHigh4 = -0.006162 ;
548545 fitParamHigh5 = 1.80808e-05 ;
549- }
550- else if (cfgnSigmaMultCuts == 3 ){
546+ } else if (cfgnSigmaMultCuts == 3 ) {
551547 fitParamHigh1 = 2610.98 ;
552548 fitParamHigh2 = -83.3983 ;
553549 fitParamHigh3 = 1.0893 ;
@@ -557,14 +553,13 @@ struct FlowSP {
557553
558554 fMultCutHigh ->SetParameters (fitParamHigh1, fitParamHigh2, fitParamHigh3, fitParamHigh4, fitParamHigh5);
559555
560- if (cfgManualEventParameters) {
561- fMultPVCutLow ->SetParameters ((cfgMultPv.value )[0 ], (cfgMultPv.value )[1 ], (cfgMultPv.value )[2 ], (cfgMultPv.value )[3 ], (cfgMultPv.value )[4 ]);
562- fMultPVCutHigh ->SetParameters ((cfgMultPv.value )[5 ], (cfgMultPv.value )[6 ], (cfgMultPv.value )[7 ], (cfgMultPv.value )[8 ], (cfgMultPv.value )[9 ]);
563- fMultCutLow ->SetParameters ((cfgMult.value )[0 ], (cfgMult.value )[1 ], (cfgMult.value )[2 ], (cfgMult.value )[3 ], (cfgMult.value )[4 ]);
564- fMultCutHigh ->SetParameters ((cfgMult.value )[5 ], (cfgMult.value )[6 ], (cfgMult.value )[7 ], (cfgMult.value )[8 ], (cfgMult.value )[9 ]);
556+ if (cfgManualEventParameters) {
557+ fMultPVCutLow ->SetParameters ((cfgMultPv.value )[0 ], (cfgMultPv.value )[1 ], (cfgMultPv.value )[2 ], (cfgMultPv.value )[3 ], (cfgMultPv.value )[4 ]);
558+ fMultPVCutHigh ->SetParameters ((cfgMultPv.value )[5 ], (cfgMultPv.value )[6 ], (cfgMultPv.value )[7 ], (cfgMultPv.value )[8 ], (cfgMultPv.value )[9 ]);
559+ fMultCutLow ->SetParameters ((cfgMult.value )[0 ], (cfgMult.value )[1 ], (cfgMult.value )[2 ], (cfgMult.value )[3 ], (cfgMult.value )[4 ]);
560+ fMultCutHigh ->SetParameters ((cfgMult.value )[5 ], (cfgMult.value )[6 ], (cfgMult.value )[7 ], (cfgMult.value )[8 ], (cfgMult.value )[9 ]);
561+ }
565562 }
566-
567- }
568563
569564 if (cfgUseAdditionalTrackCut) {
570565 fPhiCutLow = new TF1 (" fPhiCutLow" , " 0.06/x+pi/18.0-0.06" , 0 , 100 );
@@ -950,23 +945,23 @@ struct FlowSP {
950945 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" neg/hPt_hadron" ), track.pt ());
951946 }
952947
953- if (pdgCode == kPiPlus || pdgCode == kPiMinus ){
948+ if (pdgCode == kPiPlus || pdgCode == kPiMinus ) {
954949 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" incl/hPt_pion" ), track.pt ());
955- if (pdgCode == kPiPlus ) {
950+ if (pdgCode == kPiPlus ) {
956951 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" pos/hPt_pion" ), track.pt ());
957952 } else {
958953 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" neg/hPt_pion" ), track.pt ());
959954 }
960- } else if (pdgCode == kKPlus || pdgCode == kKMinus ) {
955+ } else if (pdgCode == kKPlus || pdgCode == kKMinus ) {
961956 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" incl/hPt_kaon" ), track.pt ());
962- if (pdgCode == kKPlus ) {
957+ if (pdgCode == kKPlus ) {
963958 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" pos/hPt_kaon" ), track.pt ());
964959 } else {
965960 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" neg/hPt_kaon" ), track.pt ());
966961 }
967- } else if (pdgCode == kProton || pdgCode == kProtonBar ) {
962+ } else if (pdgCode == kProton || pdgCode == kProtonBar ) {
968963 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" incl/hPt_proton" ), track.pt ());
969- if (pdgCode == kProton ) {
964+ if (pdgCode == kProton ) {
970965 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" pos/hPt_proton" ), track.pt ());
971966 } else {
972967 registry.fill (HIST (" trackMC" ) + HIST (Mode[md]) + HIST (Time[ft]) + HIST (" neg/hPt_proton" ), track.pt ());
0 commit comments