Skip to content

Add TrackTime #123

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
njacazio merged 1 commit into from
Sep 19, 2021
Merged

Add TrackTime #123

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
6 changes: 6 additions & 0 deletions examples/aod/createO2tables.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -419,14 +419,20 @@ int createO2tables(const char* inputFile = "delphes.root",
}

aod_track.fLength = track->L * 0.1; // [cm]
aod_track.fTOFChi2 = 1.f; // Negative if TOF is not available
aod_track.fTOFSignal = track->TOuter * 1.e12; // [ps]
aod_track.fTrackTime = track->TOuter * 1.e9; // [ns]
aod_track.fTrackTimeRes = 200 * 1.e9; // [ns]
aod_track.fTOFExpMom = track->P * 0.029979246;
// if primary push to TOF tracks
if (fabs(aod_track.fY) < 3. * aod_track.fSigmaY && fabs(aod_track.fZ) < 3. * aod_track.fSigmaZ)
tof_tracks.push_back(track);
} else {
aod_track.fTOFChi2 = -1.f;
aod_track.fLength = -999.f;
aod_track.fTOFSignal = -999.f;
aod_track.fTrackTime = -999.f;
aod_track.fTrackTimeRes = 2000 * 1.e9;
aod_track.fTOFExpMom = -999.f;
}

Expand Down
16 changes: 10 additions & 6 deletions examples/aod/createO2tables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -282,11 +282,13 @@ struct {
Float_t fTOFChi2 = -999.f; /// chi2 TOF match (?)

// PID
Float_t fTPCSignal = -999.f; /// dE/dX TPC
Float_t fTRDSignal = -999.f; /// dE/dX TRD
Float_t fTOFSignal = -999.f; /// TOFsignal
Float_t fLength = -999.f; /// Int.Lenght @ TOF
Float_t fTOFExpMom = -999.f; /// TOF Expected momentum based on the expected time of pions
Float_t fTPCSignal = -999.f; /// dE/dX TPC
Float_t fTRDSignal = -999.f; /// dE/dX TRD
Float_t fTOFSignal = -999.f; /// TOFsignal
Float_t fTrackTime = -999.f; /// fTrackTime
Float_t fTrackTimeRes = -999.f; /// fTrackTimeRes
Float_t fLength = -999.f; /// Int.Lenght @ TOF
Float_t fTOFExpMom = -999.f; /// TOF Expected momentum based on the expected time of pions

// Track extrapolation to EMCAL surface
Float_t fTrackEtaEMCAL = -999.f; /// Track eta at the EMCAL surface
Expand Down Expand Up @@ -349,6 +351,8 @@ void MakeTreeO2trackExtra()
tTracksExtra->Branch("fTPCSignal", &aod_track.fTPCSignal, "fTPCSignal/F");
tTracksExtra->Branch("fTRDSignal", &aod_track.fTRDSignal, "fTRDSignal/F");
tTracksExtra->Branch("fTOFSignal", &aod_track.fTOFSignal, "fTOFSignal/F");
tTracksExtra->Branch("fTrackTime", &aod_track.fTrackTime, "fTrackTime/F");
tTracksExtra->Branch("fTrackTimeRes", &aod_track.fTrackTimeRes, "fTrackTimeRes/F");
tTracksExtra->Branch("fLength", &aod_track.fLength, "fLength/F");
tTracksExtra->Branch("fTOFExpMom", &aod_track.fTOFExpMom, "fTOFExpMom/F");
tTracksExtra->Branch("fTrackEtaEMCAL", &aod_track.fTrackEtaEMCAL, "fTrackEtaEMCAL/F");
Expand All @@ -373,7 +377,7 @@ struct {
Float_t fRICHNsigmaPi = -999.f; /// Nsigma for Pi
Float_t fRICHNsigmaKa = -999.f; /// Nsigma for Ka
Float_t fRICHNsigmaPr = -999.f; /// Nsigma for Pr
} rich, frich; //! structure to keep RICH info
} rich, frich; //! structure to keep RICH info

void MakeTreeO2rich()
{
Expand Down