@@ -49,7 +49,7 @@ struct centralityStudy {
4949 int mRunNumber ;
5050 uint64_t startOfRunTimestamp;
5151
52- // vertex Z equalization
52+ // vertex Z equalization
5353 TList* hCalibObjects;
5454 TProfile* hVtxZFV0A;
5555 TProfile* hVtxZFT0A;
@@ -304,8 +304,8 @@ struct centralityStudy {
304304 o2::parameters::GRPECSObject* grpo = ccdb->getForRun <o2::parameters::GRPECSObject>(pathGRPECSObject, mRunNumber );
305305 startOfRunTimestamp = grpo->getTimeStart ();
306306
307- if (applyVertexZEqualization.value ){
308- // acquire vertex-Z equalization histograms if requested
307+ if (applyVertexZEqualization.value ) {
308+ // acquire vertex-Z equalization histograms if requested
309309 LOGF (info, " Acquiring vertex-Z profiles for run %i" , mRunNumber );
310310 hCalibObjects = ccdb->getForRun <TList>(pathVertexZ, mRunNumber );
311311
@@ -349,13 +349,13 @@ struct centralityStudy {
349349 histPointers.insert ({histPath + " hNMFTTracks" , histos.add ((histPath + " hNMFTTracks" ).c_str (), " hNMFTTracks" , {kTH1D , {{axisMultUltraFineMFTTracks}}})});
350350 histPointers.insert ({histPath + " hNPVContributors" , histos.add ((histPath + " hNPVContributors" ).c_str (), " hNPVContributors" , {kTH1D , {{axisMultUltraFinePVContributors}}})});
351351
352- if (applyVertexZEqualization){
352+ if (applyVertexZEqualization) {
353353 histPointers.insert ({histPath + " hFT0C_Collisions_Unequalized" , histos.add ((histPath + " hFT0C_Collisions_Unequalized" ).c_str (), " hFT0C_Collisions_Unequalized" , {kTH1D , {{axisMultUltraFineFT0C}}})});
354354 histPointers.insert ({histPath + " hFT0M_Collisions_Unequalized" , histos.add ((histPath + " hFT0M_Collisions_Unequalized" ).c_str (), " hFT0M_Collisions_Unequalized" , {kTH1D , {{axisMultUltraFineFT0M}}})});
355355 histPointers.insert ({histPath + " hFV0A_Collisions_Unequalized" , histos.add ((histPath + " hFV0A_Collisions_Unequalized" ).c_str (), " hFV0A_Collisions_Unequalized" , {kTH1D , {{axisMultUltraFineFV0A}}})});
356356 histPointers.insert ({histPath + " hNGlobalTracks_Unequalized" , histos.add ((histPath + " hNGlobalTracks_Unequalized" ).c_str (), " hNGlobalTracks_Unequalized" , {kTH1D , {{axisMultUltraFineGlobalTracks}}})});
357357 histPointers.insert ({histPath + " hNMFTTracks_Unequalized" , histos.add ((histPath + " hNMFTTracks_Unequalized" ).c_str (), " hNMFTTracks_Unequalized" , {kTH1D , {{axisMultUltraFineMFTTracks}}})});
358- histPointers.insert ({histPath + " hNPVContributors_Unequalized" , histos.add ((histPath + " hNPVContributors_Unequalized" ).c_str (), " hNPVContributors_Unequalized" , {kTH1D , {{axisMultUltraFinePVContributors}}})});
358+ histPointers.insert ({histPath + " hNPVContributors_Unequalized" , histos.add ((histPath + " hNPVContributors_Unequalized" ).c_str (), " hNPVContributors_Unequalized" , {kTH1D , {{axisMultUltraFinePVContributors}}})});
359359 }
360360
361361 histPointers.insert ({histPath + " hFT0CvsPVz_Collisions_All" , histos.add ((histPath + " hFT0CvsPVz_Collisions_All" ).c_str (), " hFT0CvsPVz_Collisions_All" , {kTProfile , {{axisPVz}}})});
@@ -421,38 +421,38 @@ struct centralityStudy {
421421 histos.fill (HIST (" hCollisionSelection" ), 1 );
422422 getHist (TH1, histPath + " hCollisionSelection" )->Fill (1 );
423423
424- // calculate vertex-Z-equalized quantities if desired
425- float multFV0A = collision.multFV0A ();
426- float multFT0A = collision.multFT0A ();
424+ // calculate vertex-Z-equalized quantities if desired
425+ float multFV0A = collision.multFV0A ();
426+ float multFT0A = collision.multFT0A ();
427427 float multFT0C = collision.multFT0C ();
428- float multNTracksGlobal = collision.multNTracksGlobal ();
429- float mftNtracks = collision.mftNtracks ();
430- float multNTracksPV = collision.multNTracksPV ();
431- if (applyVertexZEqualization){
428+ float multNTracksGlobal = collision.multNTracksGlobal ();
429+ float mftNtracks = collision.mftNtracks ();
430+ float multNTracksPV = collision.multNTracksPV ();
431+ if (applyVertexZEqualization) {
432432 float epsilon = 1e-2 ; // average value after which this collision will be disregarded
433433 multFV0A = -1 .0f ;
434434 multFT0A = -1 .0f ;
435435 multFT0C = -1 .0f ;
436436 multNTracksGlobal = -1 .0f ;
437- mftNtracks = -1 .0f ;
438- multNTracksPV = -1 .0f ;
439-
440- if (hVtxZFV0A->Interpolate (collision.multPVz ()) > epsilon){
437+ mftNtracks = -1 .0f ;
438+ multNTracksPV = -1 .0f ;
439+
440+ if (hVtxZFV0A->Interpolate (collision.multPVz ()) > epsilon) {
441441 multFV0A = hVtxZFV0A->Interpolate (0.0 ) * collision.multFV0A () / hVtxZFV0A->Interpolate (collision.multPVz ());
442442 }
443- if (hVtxZFT0A->Interpolate (collision.multPVz ()) > epsilon){
443+ if (hVtxZFT0A->Interpolate (collision.multPVz ()) > epsilon) {
444444 multFT0A = hVtxZFT0A->Interpolate (0.0 ) * collision.multFT0A () / hVtxZFT0A->Interpolate (collision.multPVz ());
445445 }
446- if (hVtxZFT0C->Interpolate (collision.multPVz ()) > epsilon){
446+ if (hVtxZFT0C->Interpolate (collision.multPVz ()) > epsilon) {
447447 multFT0C = hVtxZFT0C->Interpolate (0.0 ) * collision.multFT0C () / hVtxZFT0C->Interpolate (collision.multPVz ());
448448 }
449- if (hVtxZNGlobals->Interpolate (collision.multPVz ()) > epsilon){
449+ if (hVtxZNGlobals->Interpolate (collision.multPVz ()) > epsilon) {
450450 multNTracksGlobal = hVtxZNGlobals->Interpolate (0.0 ) * collision.multNTracksGlobal () / hVtxZNGlobals->Interpolate (collision.multPVz ());
451451 }
452- if (hVtxZMFT->Interpolate (collision.multPVz ()) > epsilon){
452+ if (hVtxZMFT->Interpolate (collision.multPVz ()) > epsilon) {
453453 mftNtracks = hVtxZMFT->Interpolate (0.0 ) * collision.mftNtracks () / hVtxZMFT->Interpolate (collision.multPVz ());
454454 }
455- if (hVtxZNTracks->Interpolate (collision.multPVz ()) > epsilon){
455+ if (hVtxZNTracks->Interpolate (collision.multPVz ()) > epsilon) {
456456 multNTracksPV = hVtxZNTracks->Interpolate (0.0 ) * collision.multNTracksPV () / hVtxZNTracks->Interpolate (collision.multPVz ());
457457 }
458458 }
@@ -470,9 +470,9 @@ struct centralityStudy {
470470 bool passRejectCollInTimeRangeNarrow = !(rejectCollInTimeRangeNarrow && !collision.selection_bit (o2::aod::evsel::kNoCollInTimeRangeNarrow ));
471471 // _______________________________________________________
472472 // sidestep vertex-Z rejection for vertex-Z profile histograms
473- if (passRejectITSROFBorder && passRejectTFBorder && passRequireIsVertexITSTPC && passRequireIsGoodZvtxFT0VsPV &&
474- passRequireIsVertexTOFmatched && passRequireIsVertexTRDmatched && passRejectSameBunchPileup && passRejectITSinROFpileupStandard && passRejectITSinROFpileupStrict &&
475- passSelectUPCcollisions && passRejectCollInTimeRangeNarrow){
473+ if (passRejectITSROFBorder && passRejectTFBorder && passRequireIsVertexITSTPC && passRequireIsGoodZvtxFT0VsPV &&
474+ passRequireIsVertexTOFmatched && passRequireIsVertexTRDmatched && passRejectSameBunchPileup && passRejectITSinROFpileupStandard && passRejectITSinROFpileupStrict &&
475+ passSelectUPCcollisions && passRejectCollInTimeRangeNarrow) {
476476 getHist (TProfile, histPath + " hFT0CvsPVz_Collisions_All" )->Fill (collision.multPVz (), multFT0C * scaleSignalFT0C);
477477 getHist (TProfile, histPath + " hFT0CvsPVz_Collisions" )->Fill (collision.multPVz (), multFT0C * scaleSignalFT0C);
478478 getHist (TProfile, histPath + " hFT0AvsPVz_Collisions" )->Fill (collision.multPVz (), multFT0A * scaleSignalFT0C);
@@ -606,7 +606,7 @@ struct centralityStudy {
606606 getHist (TH1, histPath + " hNGlobalTracks" )->Fill (multNTracksGlobal);
607607 getHist (TH1, histPath + " hNMFTTracks" )->Fill (mftNtracks);
608608
609- if (applyVertexZEqualization.value ) {
609+ if (applyVertexZEqualization.value ) {
610610 // save unequalized for cross-checks
611611 getHist (TH1, histPath + " hNPVContributors_Unequalized" )->Fill (collision.multNTracksPV ());
612612 getHist (TH1, histPath + " hFT0C_Collisions_Unequalized" )->Fill (collision.multFT0C () * scaleSignalFT0C);
0 commit comments