|
41 | 41 | #include <Framework/OutputObjHeader.h> |
42 | 42 | #include <Framework/runDataProcessing.h> |
43 | 43 | #include <ReconstructionDataFormats/DCA.h> |
| 44 | +#include <ReconstructionDataFormats/PID.h> |
44 | 45 |
|
45 | 46 | #include <Math/Vector4D.h> // IWYU pragma: keep (do not replace with Math/Vector4Dfwd.h) |
46 | 47 | #include <Math/Vector4Dfwd.h> |
@@ -96,7 +97,7 @@ struct skimmerPrimaryElectronFromDalitzEE { |
96 | 97 | Preslice<aod::TrackAssoc> trackIndicesPerCollision = aod::track_association::collisionId; |
97 | 98 | Produces<aod::EMPrimaryElectronsFromDalitz> emprimaryelectrons; |
98 | 99 | Produces<aod::EMPrimaryElectronsDeDxMC> emprimaryelectronsDeDxMC; |
99 | | - Service<o2::pid::tof::TOFResponse> mTOFResponse; |
| 100 | + Service<o2::pid::tof::TOFResponse> mTOFResponse{}; |
100 | 101 |
|
101 | 102 | Produces<aod::EMTOFNSigmas> emtofs; |
102 | 103 |
|
@@ -290,7 +291,7 @@ struct skimmerPrimaryElectronFromDalitzEE { |
290 | 291 | if (track.hasTOF() && track.has_collision()) { // TTCA may use orphan tracks. |
291 | 292 | auto bcTrack = track.template collision_as<TCollisions>().template bc_as<TBCs>(); |
292 | 293 | float tofNSigmaEl = mTOFResponse->nSigma<o2::track::PID::Electron>(track.tofSignalInAnotherBC(bcTrack.globalBC(), bcCollision.globalBC()), track.tofExpMom(), track.length(), track.p(), track.eta(), mapCollisionTime[collision.globalIndex()], mapCollisionTimeError[collision.globalIndex()]); |
293 | | - float beta = track.length() / (track.tofSignalInAnotherBC(bcTrack.globalBC(), bcCollision.globalBC()) - mapCollisionTime[collision.globalIndex()]) / (TMath::C() * 1e+2 * 1e-12); |
| 294 | + float beta = track.length() / (track.tofSignalInAnotherBC(bcTrack.globalBC(), bcCollision.globalBC()) - mapCollisionTime[collision.globalIndex()]) / (o2::constants::physics::LightSpeedCm2S * 1e-12); |
294 | 295 | mapTOFNsigmaReassociated[std::make_pair(collision.globalIndex(), track.globalIndex())] = tofNSigmaEl; |
295 | 296 | mapTOFBetaReassociated[std::make_pair(collision.globalIndex(), track.globalIndex())] = beta; |
296 | 297 | emtofs(collision.globalIndex(), track.globalIndex(), mapTOFBetaReassociated[std::make_pair(collision.globalIndex(), track.globalIndex())], mapTOFNsigmaReassociated[std::make_pair(collision.globalIndex(), track.globalIndex())]); |
|
0 commit comments