1212#ifndef PWGLF_UTILS_NUCLEIUTILS_H_
1313#define PWGLF_UTILS_NUCLEIUTILS_H_
1414
15- #include " Framework/HistogramRegistry.h"
16- #include " Framework/HistogramSpec.h"
17-
1815#include " Common/DataModel/Centrality.h"
1916#include " Common/DataModel/EventSelection.h"
2017#include " Common/DataModel/PIDResponse.h"
2118#include " Common/DataModel/PIDResponseITS.h"
2219#include " Common/TableProducer/PID/pidTOFBase.h"
2320
2421#include " DataFormatsTPC/BetheBlochAleph.h"
22+ #include " Framework/HistogramRegistry.h"
23+ #include " Framework/HistogramSpec.h"
2524
25+ #include < algorithm>
2626#include < string>
2727#include < vector>
28- #include < algorithm>
2928
3029using namespace o2 ;
3130using namespace o2 ::framework;
@@ -435,15 +434,17 @@ class PidManager
435434 template <typename Ttrack>
436435 float getBetaTOF (const Ttrack& track)
437436 {
438- if (!track.hasTOF ()) return -999 .f ;
437+ if (!track.hasTOF ())
438+ return -999 .f ;
439439 float beta = o2::pid::tof::Beta::GetBeta (track);
440440 return std::min (1 .f - 1 .e -6f , std::max (1 .e -4f , beta)); // / sometimes beta > 1 or < 0, to be checked
441441 }
442442
443443 template <typename Ttrack>
444444 float getMassTOF (const Ttrack& track)
445445 {
446- if (!track.hasTOF ()) return -999 .f ;
446+ if (!track.hasTOF ())
447+ return -999 .f ;
447448 const float charge{1 .f + static_cast <float >(mSpecies == Species::kHe || mSpecies == Species::kAl )};
448449 const float beta = getBetaTOF (track);
449450 return track.tpcInnerParam () * charge * std::sqrt (1 .f / (beta * beta) - 1 .f );
@@ -452,15 +453,22 @@ class PidManager
452453 template <typename Ttrack>
453454 float getNSigmaTOF (const Ttrack& track)
454455 {
455- if (!track.hasTOF ()) return -999 .f ;
456+ if (!track.hasTOF ())
457+ return -999 .f ;
456458
457459 switch (mSpecies ) {
458- case Species::kPr : return track.tofNSigmaPr ();
459- case Species::kDe : return track.tofNSigmaDe ();
460- case Species::kTr : return track.tofNSigmaTr ();
461- case Species::kHe : return track.tofNSigmaHe ();
462- case Species::kAl : return track.tofNSigmaAl ();
463- default : return -999 .f ;
460+ case Species::kPr :
461+ return track.tofNSigmaPr ();
462+ case Species::kDe :
463+ return track.tofNSigmaDe ();
464+ case Species::kTr :
465+ return track.tofNSigmaTr ();
466+ case Species::kHe :
467+ return track.tofNSigmaHe ();
468+ case Species::kAl :
469+ return track.tofNSigmaAl ();
470+ default :
471+ return -999 .f ;
464472 }
465473 }
466474
@@ -497,7 +505,8 @@ class PidManager
497505 // TPC
498506 float getExpectedTPCsignal (const float p)
499507 {
500- if (!mUseTpcCentralCalibration ) return -999 .f ;
508+ if (!mUseTpcCentralCalibration )
509+ return -999 .f ;
501510
502511 float pScaled = p * mMomScaling [0 ] + mMomScaling [1 ];
503512 float betaGamma = pScaled / masses[mSpecies ];
@@ -526,12 +535,18 @@ class PidManager
526535 float getNSigmaTPCcentral (const Ttrack& track)
527536 {
528537 switch (mSpecies ) {
529- case Species::kPr : return track.tpcNSigmaPr ();
530- case Species::kDe : return track.tpcNSigmaDe ();
531- case Species::kTr : return track.tpcNSigmaTr ();
532- case Species::kHe : return track.tpcNSigmaHe ();
533- case Species::kAl : return track.tpcNSigmaAl ();
534- default : return -999 .f ;
538+ case Species::kPr :
539+ return track.tpcNSigmaPr ();
540+ case Species::kDe :
541+ return track.tpcNSigmaDe ();
542+ case Species::kTr :
543+ return track.tpcNSigmaTr ();
544+ case Species::kHe :
545+ return track.tpcNSigmaHe ();
546+ case Species::kAl :
547+ return track.tpcNSigmaAl ();
548+ default :
549+ return -999 .f ;
535550 }
536551 }
537552
0 commit comments