Skip to content

Commit 91a5ab4

Browse files
henrylay97henrylay97
committed
Update high level CRT reco to have t0 & t1
1 parent df81a8a commit 91a5ab4

4 files changed

Lines changed: 66 additions & 36 deletions

File tree

‎sbnobj/SBND/CRT/CRTSpacePoint.cxx‎

Lines changed: 19 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -11,27 +11,34 @@ namespace sbnd {
1111
: fPos ({0., 0., 0.})
1212
, fPosErr ({0., 0., 0.})
1313
, fPE (0.)
14-
, fTime (0.)
15-
, fTimeErr (0.)
14+
, fT0 (0.)
15+
, fT0Err (0.)
16+
, fT1 (0.)
17+
, fT1Err (0.)
1618
, fComplete (false)
1719
{}
1820

1921
CRTSpacePoint::CRTSpacePoint(double _x, double _ex, double _y, double _ey, double _z, double _ez,
20-
double _pe, double _time, double _etime, bool _complete)
22+
double _pe, double _t0, double _et0, double _t1, double _et1, bool _complete)
2123
: fPos ({_x, _y, _z})
2224
, fPosErr ({_ex, _ey, _ez})
2325
, fPE (_pe)
24-
, fTime (_time)
25-
, fTimeErr (_etime)
26+
, fT0 (_t0)
27+
, fT0Err (_et0)
28+
, fT1 (_t1)
29+
, fT1Err (_et1)
2630
, fComplete (_complete)
2731
{}
2832

29-
CRTSpacePoint::CRTSpacePoint(geo::Point_t _pos, geo::Point_t _err, double _pe, double _time, double _etime, bool _complete)
33+
CRTSpacePoint::CRTSpacePoint(geo::Point_t _pos, geo::Point_t _err, double _pe, double _t0, double _et0,
34+
double _t1, double _et1, bool _complete)
3035
: fPos (_pos)
3136
, fPosErr (_err)
3237
, fPE (_pe)
33-
, fTime (_time)
34-
, fTimeErr (_etime)
38+
, fT0 (_t0)
39+
, fT0Err (_et0)
40+
, fT1 (_t1)
41+
, fT1Err (_et1)
3542
, fComplete (_complete)
3643
{}
3744

@@ -46,8 +53,10 @@ namespace sbnd {
4653
geo::Point_t CRTSpacePoint::Pos() const { return fPos; }
4754
geo::Point_t CRTSpacePoint::Err() const { return fPosErr; }
4855
double CRTSpacePoint::PE() const { return fPE; }
49-
double CRTSpacePoint::Time() const { return fTime; }
50-
double CRTSpacePoint::TimeErr() const { return fTimeErr; }
56+
double CRTSpacePoint::T0() const { return fT0; }
57+
double CRTSpacePoint::T0Err() const { return fT0Err; }
58+
double CRTSpacePoint::T1() const { return fT1; }
59+
double CRTSpacePoint::T1Err() const { return fT1Err; }
5160
bool CRTSpacePoint::Complete() const { return fComplete; }
5261
}
5362
}

‎sbnobj/SBND/CRT/CRTSpacePoint.hh‎

Lines changed: 11 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -19,18 +19,21 @@ namespace sbnd::crt {
1919
geo::Point_t fPos; // position [cm]
2020
geo::Point_t fPosErr; // positional error [cm]
2121
double fPE; // total PE
22-
double fTime; // time [ns]
23-
double fTimeErr; // time error [ns]
22+
double fT0; // time according to T0 clock [ns]
23+
double fT0Err; // error on time according to T0 clock [ns]
24+
double fT1; // time according to T1 clock [ns]
25+
double fT1Err; // error on time according to T1 clock [ns]
2426
bool fComplete; // whether or not the cluster was 3D and contained overlaps
2527

2628
public:
2729

2830
CRTSpacePoint();
2931

3032
CRTSpacePoint(double _x, double _ex, double _y, double _ey, double _z, double _ez, double _pe,
31-
double _time, double _etime, bool _complete);
33+
double _t0, double _et0, double _t1, double _et1, bool _complete);
3234

33-
CRTSpacePoint(geo::Point_t _pos, geo::Point_t _err, double _pe, double _time, double _etime, bool _complete);
35+
CRTSpacePoint(geo::Point_t _pos, geo::Point_t _err, double _pe, double _t0, double _et0,
36+
double _t1, double _et1, bool _complete);
3437

3538
virtual ~CRTSpacePoint();
3639

@@ -43,8 +46,10 @@ namespace sbnd::crt {
4346
geo::Point_t Pos() const;
4447
geo::Point_t Err() const;
4548
double PE() const;
46-
double Time() const;
47-
double TimeErr() const;
49+
double T0() const;
50+
double T0Err() const;
51+
double T1() const;
52+
double T1Err() const;
4853
bool Complete() const;
4954

5055
CRTSpacePoint& operator= (CRTSpacePoint const&) = default;

‎sbnobj/SBND/CRT/CRTTrack.cxx‎

Lines changed: 22 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -9,28 +9,36 @@ namespace sbnd {
99

1010
CRTTrack::CRTTrack()
1111
: fPoints ({})
12-
, fTime (0.)
13-
, fTimeErr (0.)
12+
, fT0 (0.)
13+
, fT0Err (0.)
14+
, fT1 (0.)
15+
, fT1Err (0.)
1416
, fPE (0.)
1517
, fToF (0.)
1618
, fTaggers ({})
1719
{}
1820

19-
CRTTrack::CRTTrack(const geo::Point_t &_start, const geo::Point_t &_end, const double &_time, const double &_etime,
20-
const double &_pe, const double &_tof, const std::set<CRTTagger> &_taggers)
21+
CRTTrack::CRTTrack(const geo::Point_t &_start, const geo::Point_t &_end, const double &_t0, const double &_et0,
22+
const double &_t1, const double &_et1, const double &_pe, const double &_tof,
23+
const std::set<CRTTagger> &_taggers)
2124
: fPoints ({_start, _end})
22-
, fTime (_time)
23-
, fTimeErr (_etime)
25+
, fT0 (_t0)
26+
, fT0Err (_et0)
27+
, fT1 (_t1)
28+
, fT1Err (_et1)
2429
, fPE (_pe)
2530
, fToF (_tof)
2631
, fTaggers (_taggers)
2732
{}
2833

29-
CRTTrack::CRTTrack(const std::vector<geo::Point_t> &_points, const double &_time, const double &_etime,
30-
const double &_pe, const double &_tof, const std::set<CRTTagger> &_taggers)
34+
CRTTrack::CRTTrack(const std::vector<geo::Point_t> &_points, const double &_t0, const double &_et0,
35+
const double &_t1, const double &_et1, const double &_pe, const double &_tof,
36+
const std::set<CRTTagger> &_taggers)
3137
: fPoints (_points)
32-
, fTime (_time)
33-
, fTimeErr (_etime)
38+
, fT0 (_t0)
39+
, fT0Err (_et0)
40+
, fT1 (_t1)
41+
, fT1Err (_et1)
3442
, fPE (_pe)
3543
, fToF (_tof)
3644
, fTaggers (_taggers)
@@ -39,8 +47,10 @@ namespace sbnd {
3947
CRTTrack::~CRTTrack() {}
4048

4149
std::vector<geo::Point_t> CRTTrack::Points() const { return fPoints; }
42-
double CRTTrack::Time() const {return fTime; }
43-
double CRTTrack::TimeErr() const { return fTimeErr; }
50+
double CRTTrack::T0() const {return fT0; }
51+
double CRTTrack::T0Err() const { return fT0Err; }
52+
double CRTTrack::T1() const {return fT1; }
53+
double CRTTrack::T1Err() const { return fT1Err; }
4454
double CRTTrack::PE() const { return fPE; }
4555
double CRTTrack::ToF() const { return fToF; }
4656
std::set<CRTTagger> CRTTrack::Taggers() const { return fTaggers; }

‎sbnobj/SBND/CRT/CRTTrack.hh‎

Lines changed: 14 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -21,8 +21,10 @@ namespace sbnd::crt {
2121
class CRTTrack {
2222

2323
std::vector<geo::Point_t> fPoints; // the fitted track points at each tagger [cm]
24-
double fTime; // average time [ns]
25-
double fTimeErr; // average time error [ns]
24+
double fT0; // average time according to T0 clock [ns]
25+
double fT0Err; // average time error according to T0 clock [ns]
26+
double fT1; // average time according to T1 clock [ns]
27+
double fT1Err; // average time error according to T1 clock [ns]
2628
double fPE; // total PE
2729
double fToF; // time from first space point to last [ns]
2830
std::set<CRTTagger> fTaggers; // which taggers were used to create the track
@@ -31,17 +33,21 @@ namespace sbnd::crt {
3133

3234
CRTTrack();
3335

34-
CRTTrack(const geo::Point_t &_start, const geo::Point_t &_end, const double &_time, const double &_etime,
35-
const double &_pe, const double &_tof, const std::set<CRTTagger> &_taggers);
36+
CRTTrack(const geo::Point_t &_start, const geo::Point_t &_end, const double &_t0, const double &_et0,
37+
const double &_t1, const double &_et1, const double &_pe, const double &_tof,
38+
const std::set<CRTTagger> &_taggers);
3639

37-
CRTTrack(const std::vector<geo::Point_t> &_points, const double &_time, const double &_etime,
38-
const double &_pe, const double &_tof, const std::set<CRTTagger> &_taggers);
40+
CRTTrack(const std::vector<geo::Point_t> &_points, const double &_t0, const double &_et0,
41+
const double &_t1, const double &_et1, const double &_pe, const double &_tof,
42+
const std::set<CRTTagger> &_taggers);
3943

4044
virtual ~CRTTrack();
4145

4246
std::vector<geo::Point_t> Points() const;
43-
double Time() const;
44-
double TimeErr() const;
47+
double T0() const;
48+
double T0Err() const;
49+
double T1() const;
50+
double T1Err() const;
4551
double PE() const;
4652
double ToF() const;
4753
std::set<CRTTagger> Taggers() const;

0 commit comments

Comments
 (0)