Skip to content

Commit 86d46c5

Browse files
henrylay97henrylay97
committed
Ensure ToF modules compile - and that their behaviour is unaffected for MC
1 parent 1ffeacb commit 86d46c5

2 files changed

Lines changed: 104 additions & 104 deletions

File tree

‎sbndcode/ToFStudies/ToFAnalyzer_module.cc‎

Lines changed: 55 additions & 55 deletions
Original file line numberDiff line numberDiff line change
@@ -672,7 +672,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
672672
std::sort(trksps.begin(),trksps.end(),
673673
[](const art::Ptr<CRTSpacePoint>& a, const art::Ptr<CRTSpacePoint>& b)->bool
674674
{
675-
return a->Time() < b->Time();
675+
return a->Ts1() < b->Ts1();
676676
});
677677
tracksps.push_back(trksps);
678678
} // Crt space points coming from are ordered on ascending order by looking into ts1_ns variable
@@ -717,7 +717,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
717717

718718
art::Ptr<sbnd::crt::CRTCluster> cluster = spToCluster.at(crt.key());
719719

720-
if(!(crt->Time() >= fBeamLow && crt->Time()<= fBeamUp)) continue;
720+
if(!(crt->Ts1() >= fBeamLow && crt->Ts1()<= fBeamUp)) continue;
721721
if(crt->PE() < fCRTSpacePointThresh) continue;
722722

723723
bool frm_trk=false;
@@ -760,7 +760,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
760760
auto const opDetPos = cryo0.OpDet(cryo0.GetClosestOpDet(point)).GetCenter();
761761
double dprop=(opDetPos - point).R();
762762
double tprop=pos.T() + dprop*LAR_PROP_DELAY;
763-
fTrue_TOF.push_back(crt->Time()-tprop);
763+
fTrue_TOF.push_back(crt->Ts1()-tprop);
764764
fTrue_TOF_hit.push_back(true);
765765
fTrue_TOF_pdg.push_back(particle->PdgCode());
766766
fTrue_TOF_part_ID.push_back(particle->TrackId());
@@ -786,7 +786,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
786786
auto const opDetPos = cryo0.OpDet(cryo0.GetClosestOpDet(point)).GetCenter();
787787
double dprop=(opDetPos - point).R();
788788
double tprop=pos.T() + dprop*LAR_PROP_DELAY;
789-
fTrue_TOF.push_back(crt->Time()-tprop);
789+
fTrue_TOF.push_back(crt->Ts1()-tprop);
790790
fTrue_TOF_hit.push_back(true);
791791
fTrue_TOF_pdg.push_back(particle->PdgCode());
792792
fTrue_TOF_part_ID.push_back(particle->TrackId());
@@ -827,7 +827,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
827827
if(hit->PE()<fHitPeThresh) continue;
828828
double thit = hit->PeakTime()*1e3-fOpDelay;
829829

830-
if(abs(crt->Time()-thit)<fCoinWindow && hit->PE()>pehit_max){
830+
if(abs(crt->Ts1()-thit)<fCoinWindow && hit->PE()>pehit_max){
831831
pehit_max = hit->PE();
832832
ophit_index = hit.key();
833833
found_tof = true;
@@ -841,12 +841,12 @@ void ToFAnalyzer::analyze(art::Event const& evt)
841841
}
842842

843843
else{
844-
fLhit_tof_vec.push_back(crt->Time() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
844+
fLhit_tof_vec.push_back(crt->Ts1() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
845845
fLhit_frmtrk_vec.push_back(false);
846846
fLhit_frmhit_vec.push_back(true);
847847
fLhit_crtspkey_vec.push_back(crt.key());
848848
fLhit_crttrkkey_vec.push_back(-9999);
849-
fLhit_crttime_t1_vec.push_back(crt->Time());
849+
fLhit_crttime_t1_vec.push_back(crt->Ts1());
850850
// fLhit_crttime_t0_vec.push_back(crt->ts0_ns);
851851
fLhit_crtpos_X_vec.push_back(crt->X());
852852
fLhit_crtpos_Y_vec.push_back(crt->Y());
@@ -883,8 +883,8 @@ void ToFAnalyzer::analyze(art::Event const& evt)
883883
if(hit->PE()<fHitPeThresh) continue;
884884
double thit = hit->PeakTime()*1e3-fOpDelay;
885885

886-
if(abs(crt->Time()-thit)<fCoinWindow && abs(crt->Time()-thit)<ophit_minTOF){
887-
ophit_minTOF = abs(crt->Time()-thit);
886+
if(abs(crt->Ts1()-thit)<fCoinWindow && abs(crt->Ts1()-thit)<ophit_minTOF){
887+
ophit_minTOF = abs(crt->Ts1()-thit);
888888
ophit_index = hit.key();
889889
found_tof = true;
890890
}
@@ -897,12 +897,12 @@ void ToFAnalyzer::analyze(art::Event const& evt)
897897
}
898898

899899
else{
900-
fChit_tof_vec.push_back(crt->Time() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
900+
fChit_tof_vec.push_back(crt->Ts1() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
901901
fChit_frmtrk_vec.push_back(false);
902902
fChit_frmhit_vec.push_back(true);
903903
fChit_crtspkey_vec.push_back(crt.key());
904904
fChit_crttrkkey_vec.push_back(-9999);
905-
fChit_crttime_t1_vec.push_back(crt->Time());
905+
fChit_crttime_t1_vec.push_back(crt->Ts1());
906906
// fChit_crttime_t0_vec.push_back(crt->ts0_ns);
907907
fChit_crtpos_X_vec.push_back(crt->X());
908908
fChit_crtpos_Y_vec.push_back(crt->Y());
@@ -943,7 +943,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
943943
auto const& flash = flashList.second[iflash];
944944
if(flash->TotalPE()<fFlashPeThresh) continue;
945945
double tflash = flash->AbsTime()*1e3-fOpDelay;
946-
if(abs(crt->Time()-tflash)<fCoinWindow && flash->TotalPE()>peflash_max){
946+
if(abs(crt->Ts1()-tflash)<fCoinWindow && flash->TotalPE()>peflash_max){
947947
peflash_max=flash->TotalPE();
948948
opflash_index = flash.key();
949949
found_tof = true;
@@ -959,12 +959,12 @@ void ToFAnalyzer::analyze(art::Event const& evt)
959959
Lflsh_tof_op_tpc[index].push_back(flash_tpc);
960960
}
961961
else{
962-
fLflsh_tof_vec.push_back(crt->Time() - (opFlashLists[flash_tpc][opflash_index]->Time()*1e3-fOpDelay));
962+
fLflsh_tof_vec.push_back(crt->Ts1() - (opFlashLists[flash_tpc][opflash_index]->Time()*1e3-fOpDelay));
963963
fLflsh_frmtrk_vec.push_back(false);
964964
fLflsh_frmhit_vec.push_back(true);
965965
fLflsh_crtspkey_vec.push_back(crt.key());
966966
fLflsh_crttrkkey_vec.push_back(-9999);
967-
fLflsh_crttime_t1_vec.push_back(crt->Time());
967+
fLflsh_crttime_t1_vec.push_back(crt->Ts1());
968968
// fLflsh_crttime_t0_vec.push_back(crt->ts0_ns);
969969
fLflsh_crtpos_X_vec.push_back(crt->X());
970970
fLflsh_crtpos_Y_vec.push_back(crt->Y());
@@ -1004,8 +1004,8 @@ void ToFAnalyzer::analyze(art::Event const& evt)
10041004
auto const& flash = flashList.second[iflash];
10051005
if(flash->TotalPE()<fFlashPeThresh) continue;
10061006
double tflash = flash->Time()*1e3-fOpDelay;
1007-
if(abs(crt->Time()-tflash)<fCoinWindow && abs(crt->Time()-tflash)<flash_minTOF){
1008-
flash_minTOF= abs(crt->Time()-tflash);
1007+
if(abs(crt->Ts1()-tflash)<fCoinWindow && abs(crt->Ts1()-tflash)<flash_minTOF){
1008+
flash_minTOF= abs(crt->Ts1()-tflash);
10091009
opflash_index = flash.key();
10101010
found_tof = true;
10111011
flash_tpc = flashList.first;
@@ -1020,12 +1020,12 @@ void ToFAnalyzer::analyze(art::Event const& evt)
10201020
Cflsh_tof_op_tpc[index].push_back(flash_tpc);
10211021
}
10221022
else{
1023-
fCflsh_tof_vec.push_back(crt->Time() - (opFlashLists[flash_tpc][opflash_index]->Time()*1e3-fOpDelay));
1023+
fCflsh_tof_vec.push_back(crt->Ts1() - (opFlashLists[flash_tpc][opflash_index]->Time()*1e3-fOpDelay));
10241024
fCflsh_frmtrk_vec.push_back(false);
10251025
fCflsh_frmhit_vec.push_back(true);
10261026
fCflsh_crtspkey_vec.push_back(crt.key());
10271027
fCflsh_crttrkkey_vec.push_back(-9999);
1028-
fCflsh_crttime_t1_vec.push_back(crt->Time());
1028+
fCflsh_crttime_t1_vec.push_back(crt->Ts1());
10291029
// fCflsh_crttime_t0_vec.push_back(crt->ts0_ns);
10301030
fCflsh_crtpos_X_vec.push_back(crt->X());
10311031
fCflsh_crtpos_Y_vec.push_back(crt->Y());
@@ -1072,7 +1072,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
10721072
auto const& flash = flashList.second[iflash];
10731073
if(flash->TotalPE()<fFlashPeThresh) continue;
10741074
double tflash = flash->AbsTime()*1e3-fOpDelay;
1075-
if(abs(crt->Time()-tflash)<fCoinWindow && flash->TotalPE()>peflash_max){
1075+
if(abs(crt->Ts1()-tflash)<fCoinWindow && flash->TotalPE()>peflash_max){
10761076
peflash_max=flash->TotalPE();
10771077
opflash_index = flash.key();
10781078
found_tof = true;
@@ -1099,12 +1099,12 @@ void ToFAnalyzer::analyze(art::Event const& evt)
10991099
}
11001100

11011101
else{
1102-
fLflshhit_tof_vec.push_back(crt->Time() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
1102+
fLflshhit_tof_vec.push_back(crt->Ts1() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
11031103
fLflshhit_frmtrk_vec.push_back(false);
11041104
fLflshhit_frmhit_vec.push_back(true);
11051105
fLflshhit_crtspkey_vec.push_back(crt.key());
11061106
fLflshhit_crttrkkey_vec.push_back(-9999);
1107-
fLflshhit_crttime_t1_vec.push_back(crt->Time());
1107+
fLflshhit_crttime_t1_vec.push_back(crt->Ts1());
11081108
// fLflshhit_crttime_t0_vec.push_back(crt->ts0_ns);
11091109
fLflshhit_crtpos_X_vec.push_back(crt->X());
11101110
fLflshhit_crtpos_Y_vec.push_back(crt->Y());
@@ -1154,8 +1154,8 @@ void ToFAnalyzer::analyze(art::Event const& evt)
11541154
auto const& flash = flashList.second[iflash];
11551155
if(flash->TotalPE()<fFlashPeThresh) continue;
11561156
double tflash = flash->AbsTime()*1e3-fOpDelay;
1157-
if(abs(crt->Time()-tflash)<fCoinWindow && abs(crt->Time()-tflash)<flash_minTOF){
1158-
flash_minTOF= abs(crt->Time()-tflash);
1157+
if(abs(crt->Ts1()-tflash)<fCoinWindow && abs(crt->Ts1()-tflash)<flash_minTOF){
1158+
flash_minTOF= abs(crt->Ts1()-tflash);
11591159
opflash_index = flash.key();
11601160
found_tof = true;
11611161
flash_tpc = flashList.first;
@@ -1181,12 +1181,12 @@ void ToFAnalyzer::analyze(art::Event const& evt)
11811181
}
11821182

11831183
else{
1184-
fCflshhit_tof_vec.push_back(crt->Time() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
1184+
fCflshhit_tof_vec.push_back(crt->Ts1() - (opHitList[ophit_index]->PeakTime()*1e3-fOpDelay));
11851185
fCflshhit_frmtrk_vec.push_back(false);
11861186
fCflshhit_frmhit_vec.push_back(true);
11871187
fCflshhit_crtspkey_vec.push_back(crt.key());
11881188
fCflshhit_crttrkkey_vec.push_back(-9999);
1189-
fCflshhit_crttime_t1_vec.push_back(crt->Time());
1189+
fCflshhit_crttime_t1_vec.push_back(crt->Ts1());
11901190
// fCflshhit_crttime_t0_vec.push_back(crt->ts0_ns);
11911191
fCflshhit_crtpos_X_vec.push_back(crt->X());
11921192
fCflshhit_crtpos_Y_vec.push_back(crt->Y());
@@ -1220,19 +1220,19 @@ void ToFAnalyzer::analyze(art::Event const& evt)
12201220
int all_index = 0;
12211221
int min_index = 0;
12221222
for (auto const& hit: ele.second){
1223-
if(hit->Time() < min_time){
1224-
min_time = hit->Time();
1223+
if(hit->Ts1() < min_time){
1224+
min_time = hit->Ts1();
12251225
min_index = all_index;
12261226
}
12271227
all_index++;
12281228
}
12291229

1230-
fLhit_tof_vec.push_back(tracksps[ele.first].front()->Time() - (Lhit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
1230+
fLhit_tof_vec.push_back(tracksps[ele.first].front()->Ts1() - (Lhit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
12311231
fLhit_frmtrk_vec.push_back(true);
12321232
fLhit_frmhit_vec.push_back(false);
12331233
fLhit_crtspkey_vec.push_back(tracksps[ele.first].front().key());
12341234
fLhit_crttrkkey_vec.push_back(ele.first);
1235-
fLhit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Time());
1235+
fLhit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Ts1());
12361236
// fLhit_crttime_t0_vec.push_back(tracksps[ele.first].front()->ts0_ns);
12371237
fLhit_crtpos_X_vec.push_back(tracksps[ele.first].front()->X());
12381238
fLhit_crtpos_Y_vec.push_back(tracksps[ele.first].front()->Y());
@@ -1262,19 +1262,19 @@ void ToFAnalyzer::analyze(art::Event const& evt)
12621262
int all_index = 0;
12631263
int min_index = 0;
12641264
for (auto const& hit: ele.second){
1265-
if(hit->Time() < min_time){
1266-
min_time = hit->Time();
1265+
if(hit->Ts1() < min_time){
1266+
min_time = hit->Ts1();
12671267
min_index = all_index;
12681268
}
12691269
all_index++;
12701270
}
12711271

1272-
fChit_tof_vec.push_back(tracksps[ele.first].front()->Time() - (Chit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
1272+
fChit_tof_vec.push_back(tracksps[ele.first].front()->Ts1() - (Chit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
12731273
fChit_frmtrk_vec.push_back(true);
12741274
fChit_frmhit_vec.push_back(false);
12751275
fChit_crtspkey_vec.push_back(tracksps[ele.first].front().key());
12761276
fChit_crttrkkey_vec.push_back(ele.first);
1277-
fChit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Time());
1277+
fChit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Ts1());
12781278
// fChit_crttime_t0_vec.push_back(tracksps[ele.first].front()->ts0_ns);
12791279
fChit_crtpos_X_vec.push_back(tracksps[ele.first].front()->X());
12801280
fChit_crtpos_Y_vec.push_back(tracksps[ele.first].front()->Y());
@@ -1304,19 +1304,19 @@ void ToFAnalyzer::analyze(art::Event const& evt)
13041304
int all_index = 0;
13051305
int min_index = 0;
13061306
for (auto const& hit: ele.second){
1307-
if(hit->Time() < min_time){
1308-
min_time = hit->Time();
1307+
if(hit->Ts1() < min_time){
1308+
min_time = hit->Ts1();
13091309
min_index = all_index;
13101310
}
13111311
all_index++;
13121312
}
13131313

1314-
fLflsh_tof_vec.push_back(tracksps[ele.first].front()->Time() - (Lflsh_tof_op_flashes[ele.first][min_index]->Time()*1e3-fOpDelay));
1314+
fLflsh_tof_vec.push_back(tracksps[ele.first].front()->Ts1() - (Lflsh_tof_op_flashes[ele.first][min_index]->Time()*1e3-fOpDelay));
13151315
fLflsh_frmtrk_vec.push_back(true);
13161316
fLflsh_frmhit_vec.push_back(false);
13171317
fLflsh_crtspkey_vec.push_back(tracksps[ele.first].front().key());
13181318
fLflsh_crttrkkey_vec.push_back(ele.first);
1319-
fLflsh_crttime_t1_vec.push_back(tracksps[ele.first].front()->Time());
1319+
fLflsh_crttime_t1_vec.push_back(tracksps[ele.first].front()->Ts1());
13201320
// fLflsh_crttime_t0_vec.push_back(tracksps[ele.first].front()->ts0_ns);
13211321
fLflsh_crtpos_X_vec.push_back(tracksps[ele.first].front()->X());
13221322
fLflsh_crtpos_Y_vec.push_back(tracksps[ele.first].front()->Y());
@@ -1345,19 +1345,19 @@ void ToFAnalyzer::analyze(art::Event const& evt)
13451345
int all_index = 0;
13461346
int min_index = 0;
13471347
for (auto const& hit: ele.second){
1348-
if(hit->Time() < min_time){
1349-
min_time = hit->Time();
1348+
if(hit->Ts1() < min_time){
1349+
min_time = hit->Ts1();
13501350
min_index = all_index;
13511351
}
13521352
all_index++;
13531353
}
13541354

1355-
fCflsh_tof_vec.push_back(tracksps[ele.first].front()->Time() - (Cflsh_tof_op_flashes[ele.first][min_index]->Time()*1e3-fOpDelay));
1355+
fCflsh_tof_vec.push_back(tracksps[ele.first].front()->Ts1() - (Cflsh_tof_op_flashes[ele.first][min_index]->Time()*1e3-fOpDelay));
13561356
fCflsh_frmtrk_vec.push_back(true);
13571357
fCflsh_frmhit_vec.push_back(false);
13581358
fCflsh_crtspkey_vec.push_back(tracksps[ele.first].front().key());
13591359
fCflsh_crttrkkey_vec.push_back(ele.first);
1360-
fCflsh_crttime_t1_vec.push_back(tracksps[ele.first].front()->Time());
1360+
fCflsh_crttime_t1_vec.push_back(tracksps[ele.first].front()->Ts1());
13611361
// fCflsh_crttime_t0_vec.push_back(tracksps[ele.first].front()->ts0_ns);
13621362
fCflsh_crtpos_X_vec.push_back(tracksps[ele.first].front()->X());
13631363
fCflsh_crtpos_Y_vec.push_back(tracksps[ele.first].front()->Y());
@@ -1386,19 +1386,19 @@ void ToFAnalyzer::analyze(art::Event const& evt)
13861386
int all_index = 0;
13871387
int min_index = 0;
13881388
for (auto const& hit: ele.second){
1389-
if(hit->Time() < min_time){
1390-
min_time = hit->Time();
1389+
if(hit->Ts1() < min_time){
1390+
min_time = hit->Ts1();
13911391
min_index = all_index;
13921392
}
13931393
all_index++;
13941394
}
13951395

1396-
fLflshhit_tof_vec.push_back(tracksps[ele.first].front()->Time() - (Lflshhit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
1396+
fLflshhit_tof_vec.push_back(tracksps[ele.first].front()->Ts1() - (Lflshhit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
13971397
fLflshhit_frmtrk_vec.push_back(true);
13981398
fLflshhit_frmhit_vec.push_back(false);
13991399
fLflshhit_crtspkey_vec.push_back(tracksps[ele.first].front().key());
14001400
fLflshhit_crttrkkey_vec.push_back(ele.first);
1401-
fLflshhit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Time());
1401+
fLflshhit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Ts1());
14021402
// fLflshhit_crttime_t0_vec.push_back(tracksps[ele.first].front()->ts0_ns);
14031403
fLflshhit_crtpos_X_vec.push_back(tracksps[ele.first].front()->X());
14041404
fLflshhit_crtpos_Y_vec.push_back(tracksps[ele.first].front()->Y());
@@ -1430,19 +1430,19 @@ void ToFAnalyzer::analyze(art::Event const& evt)
14301430
int all_index = 0;
14311431
int min_index = 0;
14321432
for (auto const& hit: ele.second){
1433-
if(hit->Time() < min_time){
1434-
min_time = hit->Time();
1433+
if(hit->Ts1() < min_time){
1434+
min_time = hit->Ts1();
14351435
min_index = all_index;
14361436
}
14371437
all_index++;
14381438
}
14391439

1440-
fCflshhit_tof_vec.push_back(tracksps[ele.first].front()->Time() - (Cflshhit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
1440+
fCflshhit_tof_vec.push_back(tracksps[ele.first].front()->Ts1() - (Cflshhit_tof_op_hits[ele.first][min_index]->PeakTime()*1e3-fOpDelay));
14411441
fCflshhit_frmtrk_vec.push_back(true);
14421442
fCflshhit_frmhit_vec.push_back(false);
14431443
fCflshhit_crtspkey_vec.push_back(tracksps[ele.first].front().key());
14441444
fCflshhit_crttrkkey_vec.push_back(ele.first);
1445-
fCflshhit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Time());
1445+
fCflshhit_crttime_t1_vec.push_back(tracksps[ele.first].front()->Ts1());
14461446
// fCflshhit_crttime_t0_vec.push_back(tracksps[ele.first].front()->ts0_ns);
14471447
fCflshhit_crtpos_X_vec.push_back(tracksps[ele.first].front()->X());
14481448
fCflshhit_crtpos_Y_vec.push_back(tracksps[ele.first].front()->Y());
@@ -1475,8 +1475,8 @@ void ToFAnalyzer::analyze(art::Event const& evt)
14751475
int all_index = 0;
14761476
int min_index = 0;
14771477
for (auto const& hit: ele.second){
1478-
if(hit->Time() < min_time){
1479-
min_time = hit->Time();
1478+
if(hit->Ts1() < min_time){
1479+
min_time = hit->Ts1();
14801480
min_index = all_index;
14811481
}
14821482
all_index++;
@@ -1490,7 +1490,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
14901490
auto const opDetPos = cryo0.OpDet(cryo0.GetClosestOpDet(point)).GetCenter();
14911491
double dprop=(opDetPos - point).R();
14921492
double tprop=pos.T() + dprop*LAR_PROP_DELAY;
1493-
fTrue_TOF.push_back(tracksps[ele.first].front()->Time()-tprop);
1493+
fTrue_TOF.push_back(tracksps[ele.first].front()->Ts1()-tprop);
14941494
fTrue_TOF_hit.push_back(false);
14951495
fTrue_TOF_pdg.push_back(True_tof_sim_particles[ele.first][min_index]->PdgCode());
14961496
fTrue_TOF_part_ID.push_back(True_tof_sim_particles[ele.first][min_index]->TrackId());
@@ -1516,7 +1516,7 @@ void ToFAnalyzer::analyze(art::Event const& evt)
15161516
auto const opDetPos = cryo0.OpDet(cryo0.GetClosestOpDet(point)).GetCenter();
15171517
double dprop=(opDetPos - point).R();
15181518
double tprop=pos.T() + dprop*LAR_PROP_DELAY;
1519-
fTrue_TOF.push_back(tracksps[ele.first].front()->Time()-tprop);
1519+
fTrue_TOF.push_back(tracksps[ele.first].front()->Ts1()-tprop);
15201520
fTrue_TOF_hit.push_back(false);
15211521
fTrue_TOF_pdg.push_back(True_tof_sim_particles[ele.first][min_index]->PdgCode());
15221522
fTrue_TOF_part_ID.push_back(True_tof_sim_particles[ele.first][min_index]->TrackId());
@@ -1544,11 +1544,11 @@ void ToFAnalyzer::analyze(art::Event const& evt)
15441544

15451545
//==========================================================================================
15461546
bool ToFAnalyzer::SpacePointCompare(const art::Ptr<CRTSpacePoint>& sp1, const art::Ptr<CRTSpacePoint>& sp2) {
1547-
if(sp1->Time() != sp2->Time()) return false;
1547+
if(sp1->Ts1() != sp2->Ts1()) return false;
15481548
if(sp1->Pos() != sp2->Pos()) return false;
15491549
if(sp1->Err() != sp2->Err()) return false;
15501550
if(sp1->PE() != sp2->PE()) return false;
1551-
if(sp1->TimeErr() != sp2->TimeErr()) return false;
1551+
if(sp1->Ts1Err() != sp2->Ts1Err()) return false;
15521552
if(sp1->Complete() != sp2->Complete()) return false;
15531553

15541554
return true;

0 commit comments

Comments
 (0)