Skip to content
Snippets Groups Projects

Introduction of chi2_time and ndf_time to AT converter

Merged Frederic Julian Linz requested to merge f.kornas/cbmroot:chi2 into master
2 files
+ 10
2
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -61,11 +61,13 @@ void CbmStsTracksConverter::Init()
@@ -61,11 +61,13 @@ void CbmStsTracksConverter::Init()
InitInput();
InitInput();
AnalysisTree::BranchConfig vtx_tracks_config(out_branch_, AnalysisTree::DetType::kTrack);
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.AddField<float>("vtx_chi2", "chi2 to to the primary vertex");
vtx_tracks_config.AddFields<float>({"dcax", "dcay", "dcaz"},
vtx_tracks_config.AddFields<float>({"dcax", "dcay", "dcaz"},
"not actuall Distance of Closest Approach, but extrapolated to z=z_vtx");
"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>("q", "charge");
vtx_tracks_config.AddField<int>("nhits", "number of hits (total MVD+STS)");
vtx_tracks_config.AddField<int>("nhits", "number of hits (total MVD+STS)");
vtx_tracks_config.AddField<int>("nhits_mvd", "number of hits in MVD");
vtx_tracks_config.AddField<int>("nhits_mvd", "number of hits in MVD");
@@ -74,7 +76,9 @@ void CbmStsTracksConverter::Init()
@@ -74,7 +76,9 @@ void CbmStsTracksConverter::Init()
iq_ = vtx_tracks_config.GetFieldId("q");
iq_ = vtx_tracks_config.GetFieldId("q");
indf_ = vtx_tracks_config.GetFieldId("ndf");
indf_ = vtx_tracks_config.GetFieldId("ndf");
 
indf_time_ = vtx_tracks_config.GetFieldId("ndf_time");
ichi2_ = vtx_tracks_config.GetFieldId("chi2");
ichi2_ = vtx_tracks_config.GetFieldId("chi2");
 
ichi2_time_ = vtx_tracks_config.GetFieldId("chi2_time");
inhits_ = vtx_tracks_config.GetFieldId("nhits");
inhits_ = vtx_tracks_config.GetFieldId("nhits");
inhits_mvd_ = vtx_tracks_config.GetFieldId("nhits_mvd");
inhits_mvd_ = vtx_tracks_config.GetFieldId("nhits_mvd");
idcax_ = vtx_tracks_config.GetFieldId("dcax");
idcax_ = vtx_tracks_config.GetFieldId("dcax");
@@ -162,7 +166,9 @@ void CbmStsTracksConverter::ReadVertexTracks(CbmEvent* event)
@@ -162,7 +166,9 @@ void CbmStsTracksConverter::ReadVertexTracks(CbmEvent* event)
track.SetMomentum3(momRec);
track.SetMomentum3(momRec);
track.SetField(int(q), iq_);
track.SetField(int(q), iq_);
track.SetField(int(sts_track->GetNDF()), indf_);
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->GetChiSq()), ichi2_);
 
track.SetField(float(sts_track->GetChiSqTime()), ichi2_time_);
track.SetField(int(sts_track->GetNofHits()), inhits_);
track.SetField(int(sts_track->GetNofHits()), inhits_);
track.SetField(float(trackParamFirst->GetX() - cbm_prim_vertex_->GetX()), idcax_);
track.SetField(float(trackParamFirst->GetX() - cbm_prim_vertex_->GetX()), idcax_);
track.SetField(float(trackParamFirst->GetY() - cbm_prim_vertex_->GetY()), idcax_ + 1);
track.SetField(float(trackParamFirst->GetY() - cbm_prim_vertex_->GetY()), idcax_ + 1);
Loading