diff --git a/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx b/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx index 59a2be962c55923976f578c7e403cd2aae64a899..c0761ca3ec1ad0751d8328aa9fad6fa9c6ef7f1a 100644 --- a/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx +++ b/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx @@ -18,6 +18,8 @@ #include <TH1F.h> #include <numeric> +using namespace Acts::UnitLiterals; + NtupleDumperAlg::NtupleDumperAlg(const std::string &name, ISvcLocator *pSvcLocator) : AthReentrantAlgorithm(name, pSvcLocator), @@ -1027,7 +1029,7 @@ StatusCode NtupleDumperAlg::execute(const EventContext &ctx) const params_up[Acts::eBoundLoc1] = position_up.x(); params_up[Acts::eBoundPhi] = momentum_up.phi(); params_up[Acts::eBoundTheta] = momentum_up.theta(); - params_up[Acts::eBoundQOverP] = upstreamParameters->charge() / momentum_up.mag(); + params_up[Acts::eBoundQOverP] = upstreamParameters->charge() / (momentum_up.mag() * 1_MeV); params_up[Acts::eBoundTime] = 0; auto startSurface_up = Acts::Surface::makeShared<Acts::PlaneSurface>(Acts::Vector3(0, 0, position_up.z()), Acts::Vector3(0, 0, 1)); Acts::BoundTrackParameters startParameters_up(std::move(startSurface_up), params_up, upstreamParameters->charge()); @@ -1040,7 +1042,7 @@ StatusCode NtupleDumperAlg::execute(const EventContext &ctx) const params_down[Acts::eBoundLoc1] = position_down.x(); params_down[Acts::eBoundPhi] = momentum_down.phi(); params_down[Acts::eBoundTheta] = momentum_down.theta(); - params_down[Acts::eBoundQOverP] = downstreamParameters->charge() / momentum_down.mag(); + params_down[Acts::eBoundQOverP] = downstreamParameters->charge() / (momentum_down.mag() * 1_MeV); params_down[Acts::eBoundTime] = 0; auto startSurface_down = Acts::Surface::makeShared<Acts::PlaneSurface>(Acts::Vector3(0, 0, position_down.z()), Acts::Vector3(0, 0, 1)); Acts::BoundTrackParameters startParameters_down(std::move(startSurface_down), params_down, downstreamParameters->charge());