Sign of sigma shouldn't flip with bz

mpuccio · shahor02 · commit 2bcf8d5d3529 · 2022-09-20T15:32:22.000+02:00
Spotted by @wang-yuanzhe
diff --git a/Detectors/Vertexing/src/SVertexHypothesis.cxx b/Detectors/Vertexing/src/SVertexHypothesis.cxx
@@ -25,7 +25,8 @@ void SVertexHypothesis::set(PID v0, PID ppos, PID pneg, float sig, float nSig, f
   mPars[SigmaM] = sig;
   mPars[NSigmaM] = nSig;
   mPars[MarginM] = margin;
-  mPars[CPt] = std::abs(bz) > 1e-3 ? cpt * 5.0066791 / bz : 0.; // assume that pT dependent sigma is linear with B
+  float absBz{std::abs(bz)};
+  mPars[CPt] = absBz > 1e-3 ? cpt * 5.0066791 / absBz : 0.; // assume that pT dependent sigma is linear with B
 }
 
 void SVertexHypothesis::set(PID v0, PID ppos, PID pneg, const float pars[NPIDParams], float bz)

Original file line number	Diff line number	Diff line change
`@@ -25,7 +25,8 @@ void SVertexHypothesis::set(PID v0, PID ppos, PID pneg, float sig, float nSig, f`
`25`	`25`	`mPars[SigmaM] = sig;`
`26`	`26`	`mPars[NSigmaM] = nSig;`
`27`	`27`	`mPars[MarginM] = margin;`
`28`		`- mPars[CPt] = std::abs(bz) > 1e-3 ? cpt * 5.0066791 / bz : 0.; // assume that pT dependent sigma is linear with B`
	`28`	`+ float absBz{std::abs(bz)};`
	`29`	`+ mPars[CPt] = absBz > 1e-3 ? cpt * 5.0066791 / absBz : 0.; // assume that pT dependent sigma is linear with B`
`29`	`30`	`}`
`30`	`31`
`31`	`32`	`void SVertexHypothesis::set(PID v0, PID ppos, PID pneg, const float pars[NPIDParams], float bz)`