Skip to content
Snippets Groups Projects
Commit 71cea08a authored by Viktor's avatar Viktor
Browse files

add some titles, remove muons from TRD probabilities

parent 17b99c4e
No related branches found
No related tags found
1 merge request!272Update AnalysisTree
...@@ -19,7 +19,7 @@ ClassImp(CbmSimEventHeaderConverter) ...@@ -19,7 +19,7 @@ ClassImp(CbmSimEventHeaderConverter)
AnalysisTree::BranchConfig SimEventHeaderBranch( AnalysisTree::BranchConfig SimEventHeaderBranch(
"SimEventHeader", AnalysisTree::DetType::kEventHeader); "SimEventHeader", AnalysisTree::DetType::kEventHeader);
SimEventHeaderBranch.AddField<float>("psi_RP", "reaction plane orientation"); SimEventHeaderBranch.AddField<float>("psi_RP", "reaction plane orientation");
SimEventHeaderBranch.AddField<float>("b", "impact parameter"); SimEventHeaderBranch.AddField<float>("b", "fm, impact parameter");
SimEventHeaderBranch.AddField<float>("time", "ns, event time"); SimEventHeaderBranch.AddField<float>("time", "ns, event time");
SimEventHeaderBranch.AddField<int>("run_id", "run identifier"); SimEventHeaderBranch.AddField<int>("run_id", "run identifier");
SimEventHeaderBranch.AddField<int>("event_id", "event identifier"); SimEventHeaderBranch.AddField<int>("event_id", "event identifier");
......
...@@ -32,13 +32,13 @@ ClassImp(CbmTofHitsConverter) ...@@ -32,13 +32,13 @@ ClassImp(CbmTofHitsConverter)
AnalysisTree::BranchConfig tof_branch(out_branch_, AnalysisTree::BranchConfig tof_branch(out_branch_,
AnalysisTree::DetType::kHit); AnalysisTree::DetType::kHit);
tof_branch.AddField<float>("mass2"); tof_branch.AddField<float>("mass2", "Mass squared");
tof_branch.AddField<float>("l"); tof_branch.AddField<float>("l", "Track lenght");
tof_branch.AddField<float>("t"); tof_branch.AddField<float>("t", "ps(?), Measured time ");
tof_branch.AddField<float>("qp_tof"); tof_branch.AddField<float>("qp_tof", "charge * momentum extrapoleted to TOF");
tof_branch.AddFields<float>({"dx", "dy", "dz"}); tof_branch.AddFields<float>({"dx", "dy", "dz"}, "Distance between TOF hit and extrapolated global track, cm");
tof_branch.AddField<int>("mc_pdg"); tof_branch.AddField<int>("mc_pdg", "MC-true PDG code of particle with highest contribution to TOF hit");
tof_branch.AddField<bool>("is_correct_match"); tof_branch.AddField<bool>("is_correct_match", "is the matched track corresponds to MC-true track from TOF hit");
auto* man = AnalysisTree::TaskManager::GetInstance(); auto* man = AnalysisTree::TaskManager::GetInstance();
man->AddBranch(out_branch_, tof_hits_, tof_branch); man->AddBranch(out_branch_, tof_hits_, tof_branch);
......
...@@ -32,12 +32,12 @@ ClassImp(CbmTrdTracksConverter) ...@@ -32,12 +32,12 @@ ClassImp(CbmTrdTracksConverter)
AnalysisTree::BranchConfig trd_branch(out_branch_, AnalysisTree::BranchConfig trd_branch(out_branch_,
AnalysisTree::DetType::kTrack); AnalysisTree::DetType::kTrack);
trd_branch.AddFields<float>( trd_branch.AddFields<float>(
{"energy_loss_0", "energy_loss_1", "energy_loss_2", "energy_loss_3"}); {"energy_loss_0", "energy_loss_1", "energy_loss_2", "energy_loss_3"}, "Energy loss per TRD station");
trd_branch.AddFields<float>( trd_branch.AddFields<float>(
{"pid_like_e", "pid_like_pi", "pid_like_k", "pid_like_p", "pid_like_mu"}); {"pid_like_e", "pid_like_pi", "pid_like_k", "pid_like_p"}, "Probability to be a given particle specie");
trd_branch.AddField<float>("chi2_ov_ndf"); trd_branch.AddField<float>("chi2_ov_ndf", "chi2 divided by NDF of the track fit");
trd_branch.AddFields<float>({"pT_out", "p_out"}); trd_branch.AddFields<float>({"pT_out", "p_out"}, "Momentum at last point (?)");
trd_branch.AddField<int>("n_hits"); trd_branch.AddField<int>("n_hits", "Number of hits");
auto* man = AnalysisTree::TaskManager::GetInstance(); auto* man = AnalysisTree::TaskManager::GetInstance();
man->AddBranch(out_branch_, trd_tracks_, trd_branch); man->AddBranch(out_branch_, trd_tracks_, trd_branch);
...@@ -94,7 +94,7 @@ void CbmTrdTracksConverter::FillTrdTracks() { ...@@ -94,7 +94,7 @@ void CbmTrdTracksConverter::FillTrdTracks() {
track.SetField(float(trd_track->GetPidLikePI()), i_pid_like + 1); track.SetField(float(trd_track->GetPidLikePI()), i_pid_like + 1);
track.SetField(float(trd_track->GetPidLikeKA()), i_pid_like + 2); track.SetField(float(trd_track->GetPidLikeKA()), i_pid_like + 2);
track.SetField(float(trd_track->GetPidLikePR()), i_pid_like + 3); track.SetField(float(trd_track->GetPidLikePR()), i_pid_like + 3);
track.SetField(float(trd_track->GetPidLikeMU()), i_pid_like + 4); // track.SetField(float(trd_track->GetPidLikeMU()), i_pid_like + 4);
track.SetField(float(trd_track->GetNDF() > 0. track.SetField(float(trd_track->GetNDF() > 0.
? trd_track->GetChiSq() / trd_track->GetNDF() ? trd_track->GetChiSq() / trd_track->GetNDF()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment