Fix: fix reduced mass

ariedel-cern · ariedel-cern · commit 8ea5acac2dd5 · 2025-10-31T14:53:49.000+01:00
diff --git a/PWGCF/Femto/Core/pairHistManager.h b/PWGCF/Femto/Core/pairHistManager.h
@@ -280,7 +280,7 @@ class PairHistManager
     mPdgMass1 = o2::analysis::femto::utils::getMass(PdgParticle1);
     mPdgMass2 = o2::analysis::femto::utils::getMass(PdgParticle2);
     mPdgAverageMass = (mPdgMass1 + mPdgMass2) / 2.f;
-    mPdgReducedMass = (mPdgMass1 * mPdgMass2) / (mPdgMass1 + mPdgMass2);
+    mPdgReducedMass = 2 * (mPdgMass1 * mPdgMass2) / (mPdgMass1 + mPdgMass2);
   }
   void setCharge(int chargeAbsParticle1, int chargeAbsParticle2)
   {

Original file line number	Diff line number	Diff line change
`@@ -280,7 +280,7 @@ class PairHistManager`
`280`	`280`	`mPdgMass1 = o2::analysis::femto::utils::getMass(PdgParticle1);`
`281`	`281`	`mPdgMass2 = o2::analysis::femto::utils::getMass(PdgParticle2);`
`282`	`282`	`mPdgAverageMass = (mPdgMass1 + mPdgMass2) / 2.f;`
`283`		`- mPdgReducedMass = (mPdgMass1 * mPdgMass2) / (mPdgMass1 + mPdgMass2);`
	`283`	`+ mPdgReducedMass = 2 * (mPdgMass1 * mPdgMass2) / (mPdgMass1 + mPdgMass2);`
`284`	`284`	`}`
`285`	`285`	`void setCharge(int chargeAbsParticle1, int chargeAbsParticle2)`
`286`	`286`	`{`