Skip to content
Snippets Groups Projects
Commit c0422ac6 authored by fjlinz's avatar fjlinz
Browse files

Introduction of chi2_time and ndf_time to AT converter

parent 72a5cac2
No related branches found
No related tags found
1 merge request!1469Introduction of chi2_time and ndf_time to AT converter
Pipeline #25252 passed
......@@ -61,11 +61,13 @@ void CbmStsTracksConverter::Init()
InitInput();
AnalysisTree::BranchConfig vtx_tracks_config(out_branch_, AnalysisTree::DetType::kTrack);
vtx_tracks_config.AddField<float>("chi2", "chi2 of the track fit");
vtx_tracks_config.AddField<float>("chi2", "spatial chi2 of the track fit");
vtx_tracks_config.AddField<float>("chi2_time", "time chi2 of the track fit");
vtx_tracks_config.AddField<float>("vtx_chi2", "chi2 to to the primary vertex");
vtx_tracks_config.AddFields<float>({"dcax", "dcay", "dcaz"},
"not actuall Distance of Closest Approach, but extrapolated to z=z_vtx");
vtx_tracks_config.AddField<int>("ndf", "number degrees of freedom");
vtx_tracks_config.AddField<int>("ndf", "spatial number degrees of freedom");
vtx_tracks_config.AddField<int>("ndf_time", "time number degrees of freedom");
vtx_tracks_config.AddField<int>("q", "charge");
vtx_tracks_config.AddField<int>("nhits", "number of hits (total MVD+STS)");
vtx_tracks_config.AddField<int>("nhits_mvd", "number of hits in MVD");
......@@ -74,7 +76,9 @@ void CbmStsTracksConverter::Init()
iq_ = vtx_tracks_config.GetFieldId("q");
indf_ = vtx_tracks_config.GetFieldId("ndf");
indf_time_ = vtx_tracks_config.GetFieldId("ndf_time");
ichi2_ = vtx_tracks_config.GetFieldId("chi2");
ichi2_time_ = vtx_tracks_config.GetFieldId("chi2_time");
inhits_ = vtx_tracks_config.GetFieldId("nhits");
inhits_mvd_ = vtx_tracks_config.GetFieldId("nhits_mvd");
idcax_ = vtx_tracks_config.GetFieldId("dcax");
......@@ -162,7 +166,9 @@ void CbmStsTracksConverter::ReadVertexTracks(CbmEvent* event)
track.SetMomentum3(momRec);
track.SetField(int(q), iq_);
track.SetField(int(sts_track->GetNDF()), indf_);
track.SetField(int(sts_track->GetNdfTime()), indf_time_);
track.SetField(float(sts_track->GetChiSq()), ichi2_);
track.SetField(float(sts_track->GetChiSqTime()), ichi2_time_);
track.SetField(int(sts_track->GetNofHits()), inhits_);
track.SetField(float(trackParamFirst->GetX() - cbm_prim_vertex_->GetX()), idcax_);
track.SetField(float(trackParamFirst->GetY() - cbm_prim_vertex_->GetY()), idcax_ + 1);
......
......@@ -63,7 +63,9 @@ private:
int iq_ {AnalysisTree::UndefValueInt};
int indf_ {AnalysisTree::UndefValueInt};
int indf_time_ {AnalysisTree::UndefValueInt};
int ichi2_ {AnalysisTree::UndefValueInt};
int ichi2_time_ {AnalysisTree::UndefValueInt};
int inhits_ {AnalysisTree::UndefValueInt};
int inhits_mvd_ {AnalysisTree::UndefValueInt};
int idcax_ {AnalysisTree::UndefValueInt};
......
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