2222#include " PWGEM/PhotonMeson/Utils/TrackSelection.h"
2323
2424#include < CCDB/CcdbApi.h>
25+ #include < CommonConstants/MathConstants.h>
2526#include < Framework/ASoA.h>
2627#include < Framework/Array2D.h>
2728
4142#include < utility>
4243#include < vector>
4344
44- using namespace o2 ::pwgem::photonmeson;
45-
4645namespace o2 ::analysis
4746{
4847
@@ -263,7 +262,7 @@ class V0PhotonCut : public TNamed
263262 auto pos = v0.template posTrack_as <TLeg>();
264263 auto ele = v0.template negTrack_as <TLeg>();
265264
266- for (auto & track : {pos, ele}) {
265+ for (const auto & track : {pos, ele}) {
267266 if (!IsSelectedTrack (track, V0PhotonCuts::kTrackPtRange )) {
268267 return false ;
269268 }
@@ -279,24 +278,24 @@ class V0PhotonCut : public TNamed
279278 if (!track.hasITS () && !track.hasTPC ()) { // track has to be ITSonly or TPConly or ITS-TPC
280279 return false ;
281280 }
282- if (mDisableITSonly && isITSonlyTrack (track)) {
281+ if (mDisableITSonly && o2::pwgem::photonmeson:: isITSonlyTrack (track)) {
283282 return false ;
284283 }
285- if (mDisableTPConly && isTPConlyTrack (track)) {
284+ if (mDisableTPConly && o2::pwgem::photonmeson:: isTPConlyTrack (track)) {
286285 return false ;
287286 }
288287
289288 if (mRejectITSib ) {
290289 auto hits_ib = std::count_if (its_ib_Requirement.second .begin (), its_ib_Requirement.second .end (), [&](auto && requiredLayer) { return track.itsClusterMap () & (1 << requiredLayer); });
291290 auto hits_ob = std::count_if (its_ob_Requirement.second .begin (), its_ob_Requirement.second .end (), [&](auto && requiredLayer) { return track.itsClusterMap () & (1 << requiredLayer); });
292291 bool its_ob_only = (hits_ib <= its_ib_Requirement.first ) && (hits_ob >= its_ob_Requirement.first );
293- if (isITSonlyTrack (track) && !its_ob_only) { // ITSonly tracks should not have any ITSib hits.
292+ if (o2::pwgem::photonmeson:: isITSonlyTrack (track) && !its_ob_only) { // ITSonly tracks should not have any ITSib hits.
294293 return false ;
295294 }
296295
297296 auto hits_ob_itstpc = std::count_if (its_ob_Requirement_ITSTPC.second .begin (), its_ob_Requirement_ITSTPC.second .end (), [&](auto && requiredLayer) { return track.itsClusterMap () & (1 << requiredLayer); });
298297 bool its_ob_only_itstpc = (hits_ib <= its_ib_Requirement.first ) && (hits_ob_itstpc >= its_ob_Requirement_ITSTPC.first );
299- if (isITSTPCTrack (track) && !its_ob_only_itstpc) { // ITSTPC tracks should not have any ITSib hits.
298+ if (o2::pwgem::photonmeson:: isITSTPCTrack (track) && !its_ob_only_itstpc) { // ITSTPC tracks should not have any ITSib hits.
300299 return false ;
301300 }
302301 }
@@ -416,7 +415,7 @@ class V0PhotonCut : public TNamed
416415 return v0.mGamma () < mMaxQt * 2 .f ;
417416
418417 case V0PhotonCuts::kAP :
419- return pow (v0.alpha () / mMaxAlpha , 2 ) + pow (v0.qtarm () / mMaxQt , 2 ) < 1.0 ;
418+ return std:: pow (v0.alpha () / mMaxAlpha , 2 ) + std:: pow (v0.qtarm () / mMaxQt , 2 ) < 1.0 ;
420419
421420 case V0PhotonCuts::kPsiPair :
422421 return true ;
@@ -454,24 +453,24 @@ class V0PhotonCut : public TNamed
454453 float y = v0.vy (); // cm, measured secondary vertex of gamma->ee
455454 float z = v0.vz (); // cm, measured secondary vertex of gamma->ee
456455
457- float rxy = sqrt (x * x + y * y);
456+ float rxy = std:: sqrt (x * x + y * y);
458457 if (rxy < 7.0 || 14.0 < rxy) {
459458 return false ;
460459 }
461460
462- // r = 0.192 * z + 8.88 (cm) expected wire position in RZ plane.TMath::Tan(10.86 * TMath::DegToRad() ) = 0.192
461+ // r = 0.192 * z + 8.88 (cm) expected wire position in RZ plane.TMath::Tan(10.86 * o2::constants::math::Deg2Rad ) = 0.192
463462 if (rxy > 0.192 * z + 14.0 ) { // upper limit
464463 return false ;
465464 }
466465
467- float dxy = std::fabs (1.0 * y - x * std::tan (-8.52 * TMath::DegToRad ())) / sqrt (pow (1.0 , 2 ) + pow (std::tan (-8.52 * TMath::DegToRad () ), 2 ));
466+ float dxy = std::fabs (1.0 * y - x * std::tan (-8 .52f * o2::constants::math::Deg2Rad)) / std:: sqrt (std:: pow (1.0 , 2 ) + std:: pow (std::tan (-8 .52f * o2::constants::math::Deg2Rad ), 2 ));
468467 return !(dxy > margin_xy);
469468 }
470469 case V0PhotonCuts::kOnWwireOB : {
471- const float margin_xy = 1.0 ; // cm
472- const float rxy_exp = 30.8 ; // cm
473- const float x_exp = rxy_exp * std::cos (-1.3 * TMath::DegToRad () ); // cm, expected position x of W wire
474- const float y_exp = rxy_exp * std::sin (-1.3 * TMath::DegToRad () ); // cm, expected position y of W wire
470+ const float margin_xy = 1.0 ; // cm
471+ const float rxy_exp = 30.8 ; // cm
472+ const float x_exp = rxy_exp * std::cos (-1 .3f * o2::constants::math::Deg2Rad ); // cm, expected position x of W wire
473+ const float y_exp = rxy_exp * std::sin (-1 .3f * o2::constants::math::Deg2Rad ); // cm, expected position y of W wire
475474 // const float z_min = -47.0; // cm
476475 // const float z_max = +47.0; // cm
477476 float x = v0.vx (); // cm, measured secondary vertex of gamma->ee
@@ -536,26 +535,26 @@ class V0PhotonCut : public TNamed
536535 return mMinChi2PerClusterITS < track.itsChi2NCl () && track.itsChi2NCl () < mMaxChi2PerClusterITS ;
537536
538537 case V0PhotonCuts::kITSClusterSize : {
539- if (!isITSonlyTrack (track)) {
538+ if (!o2::pwgem::photonmeson:: isITSonlyTrack (track)) {
540539 return true ;
541540 }
542541 return mMinMeanClusterSizeITS < track.meanClusterSizeITSob () * std::cos (std::atan (track.tgl ())) && track.meanClusterSizeITSob () * std::cos (std::atan (track.tgl ())) < mMaxMeanClusterSizeITS ;
543542 }
544543
545544 case V0PhotonCuts::kRequireITSTPC :
546- return isITSTPCTrack (track);
545+ return o2::pwgem::photonmeson:: isITSTPCTrack (track);
547546
548547 case V0PhotonCuts::kRequireITSonly :
549- return isITSonlyTrack (track);
548+ return o2::pwgem::photonmeson:: isITSonlyTrack (track);
550549
551550 case V0PhotonCuts::kRequireTPConly :
552- return isTPConlyTrack (track);
551+ return o2::pwgem::photonmeson:: isTPConlyTrack (track);
553552
554553 case V0PhotonCuts::kRequireTPCTRD :
555- return isTPCTRDTrack (track);
554+ return o2::pwgem::photonmeson:: isTPCTRDTrack (track);
556555
557556 case V0PhotonCuts::kRequireTPCTOF :
558- return isTPCTOFTrack (track);
557+ return o2::pwgem::photonmeson:: isTPCTOFTrack (track);
559558
560559 default :
561560 return false ;
0 commit comments