Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions PWGEM/Dilepton/Utils/EMTrack.h
Original file line number Diff line number Diff line change
Expand Up @@ -53,7 +53,7 @@
float cZY() const { return fCZY; }
float cZZ() const { return fCZZ; }

float p() const { return 1.f / std::fabs(fSigned1Pt) * std::cosh(fEta); }
float p() const { return 1.f / std::fabs(fSigned1Pt) * std::cosh(fEta); }
float px() const { return 1.f / std::fabs(fSigned1Pt) * std::cos(fPhi); }
float py() const { return 1.f / std::fabs(fSigned1Pt) * std::sin(fPhi); }
float pz() const { return 1.f / std::fabs(fSigned1Pt) * std::sinh(fEta); }
Expand Down Expand Up @@ -163,7 +163,7 @@

~EMPair() {}
float mass() const { return fMass; }
float rapidity() const { return std::log((std::sqrt(std::pow(fMass, 2) + std::pow(1.f/ std::fabs(fSigned1Pt) * std::cosh(fEta), 2)) + 1.f/ std::fabs(fSigned1Pt) * std::sinh(fEta)) / std::sqrt(std::pow(fMass, 2) + std::pow(1.f/ std::fabs(fSigned1Pt), 2))); }
float rapidity() const { return std::log((std::sqrt(std::pow(fMass, 2) + std::pow(1.f / std::fabs(fSigned1Pt) * std::cosh(fEta), 2)) + 1.f / std::fabs(fSigned1Pt) * std::sinh(fEta)) / std::sqrt(std::pow(fMass, 2) + std::pow(1.f / std::fabs(fSigned1Pt), 2))); }

void setPairDCA(float dca) { fPairDCA = dca; }
float getPairDCA() const { return fPairDCA; }
Expand All @@ -189,7 +189,7 @@
float eta = std::atanh(pz / std::sqrt(std::pow(px, 2) + std::pow(py, 2) + std::pow(pz, 2)));
float phi = std::atan2(py, px);
if (phi < 0.f) {
phi += 2.f * M_PI;

Check failure on line 192 in PWGEM/Dilepton/Utils/EMTrack.h

View workflow job for this annotation

GitHub Actions / O2 linter

[two-pi-add-subtract]

Use RecoDecay::constrainAngle to restrict angle to a given range.
}

fVPos.SetPt(pt);
Expand All @@ -203,7 +203,7 @@
float eta = std::atanh(pz / std::sqrt(std::pow(px, 2) + std::pow(py, 2) + std::pow(pz, 2)));
float phi = std::atan2(py, px);
if (phi < 0.f) {
phi += 2.f * M_PI;

Check failure on line 206 in PWGEM/Dilepton/Utils/EMTrack.h

View workflow job for this annotation

GitHub Actions / O2 linter

[two-pi-add-subtract]

Use RecoDecay::constrainAngle to restrict angle to a given range.
}

fVNeg.SetPt(pt);
Expand Down
2 changes: 1 addition & 1 deletion PWGEM/PhotonMeson/Utils/EMPhoton.h
Original file line number Diff line number Diff line change
Expand Up @@ -42,7 +42,7 @@ class EMPhoton
float eta() const { return fEta; }
float phi() const { return fPhi; }

float p() const { return fPt * std::cosh(fEta); }
float p() const { return fPt * std::cosh(fEta); }
float px() const { return fPt * std::cos(fPhi); }
float py() const { return fPt * std::sin(fPhi); }
float pz() const { return fPt * std::sinh(fEta); }
Expand Down
Loading