Skip to content

Adding CRT and PDS info in calib ntuple #422

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
fjnicolas merged 16 commits into from
Mar 22, 2024
Merged

Adding CRT and PDS info in calib ntuple #422

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
16 commits
Select commit Hold shift + click to select a range
b6c0926
A commit to add CRT and PDS info in calib ntuple using updated reco2_…
sungbinoh Feb 29, 2024
234018d
Merge branch 'develop' into feature/sungbinoh_crt_pds_for_calib
fjnicolas Mar 1, 2024
c4c238b
Re-organize the fcl files for calibration skimmers
linyan-w Mar 11, 2024
e1e23ad
Re-organize the calibration skimmer fcls
linyan-w Mar 11, 2024
41f0e90
Bug fix for the last commit...
linyan-w Mar 12, 2024
c3af2ef
Merge branch 'feature/sungbinoh_crt_pds_for_calib' of https://github.…
linyan-w Mar 12, 2024
1ffea3b
Add comments to explain the CRT variables
linyan-w Mar 12, 2024
c2876f2
Remove some comments to tidy up
linyan-w Mar 12, 2024
5429cc2
Merge branch 'develop' into feature/sungbinoh_crt_pds_for_calib
linyan-w Mar 12, 2024
e2c6a16
Move PDS (calibration) to Calibration
linyan-w Mar 12, 2024
cb68545
Re-organize CRT fcl
linyan-w Mar 12, 2024
2bf4e75
Re-organize CRT fcl
linyan-w Mar 12, 2024
c03f524
Move crt config out of reco_sbnd
linyan-w Mar 12, 2024
0249bcd
Clock etc
linyan-w Mar 12, 2024
da5e56d
Update sbndcode/JobConfigurations/standard/reco/reco2_sce.fcl
linyan-w Mar 12, 2024
990167f
Update sbndcode/JobConfigurations/standard/reco/reco2_sce.fcl
linyan-w Mar 12, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions sbndcode/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ add_subdirectory(OpDetSim)
add_subdirectory(OpDetReco)
add_subdirectory(OpDetAnalyzer)
add_subdirectory(OpT0Finder)
#add_subdirectory(PDS)
#add_subdirectory(CosmicId)
add_subdirectory(ShowerAna)
add_subdirectory(Commissioning)
Expand Down
7 changes: 7 additions & 0 deletions sbndcode/CRT/CRTAna/CRTAnalysis_module.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -96,6 +96,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
int _subrun;
int _event;

//mc truth
std::vector<int16_t> _mc_trackid;
std::vector<int16_t> _mc_pdg;
std::vector<int16_t> _mc_status;
Expand All @@ -118,6 +119,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
std::vector<double> _mc_endpz;
std::vector<double> _mc_ende;

//G4 detector id
std::vector<int16_t> _ide_trackid;
std::vector<float> _ide_e;
std::vector<float> _ide_entryx;
Expand All @@ -129,6 +131,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
std::vector<float> _ide_exitz;
std::vector<float> _ide_exitt;

//front end mother board
std::vector<uint16_t> _feb_mac5;
std::vector<uint16_t> _feb_flags;
std::vector<uint32_t> _feb_ts0;
Expand All @@ -137,6 +140,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
std::vector<std::vector<uint16_t>> _feb_adc;
std::vector<uint32_t> _feb_coinc;

//strip hit to select the strip which has ADC above threshold
std::vector<uint32_t> _sh_channel;
std::vector<uint32_t> _sh_ts0;
std::vector<uint32_t> _sh_ts1;
Expand All @@ -154,6 +158,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
std::vector<double> _sh_truth_energy;
std::vector<double> _sh_truth_time;

//cluster from x-y coincidence for CRTSpacePoint , this is what we normally call a CRT hit
std::vector<uint32_t> _cl_ts0;
std::vector<uint32_t> _cl_ts1;
std::vector<uint32_t> _cl_unixs;
Expand Down Expand Up @@ -188,6 +193,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
std::vector<double> _cl_sp_etime;
std::vector<bool> _cl_sp_complete;

//backtrack truth information from reco level
std::vector<int> _td_tag_trackid;
std::vector<int> _td_tag_pdg;
std::vector<int16_t> _td_tag_tagger;
Expand All @@ -206,6 +212,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
std::vector<bool> _td_reco_status;
std::vector<bool> _td_reco_triple;

//track level information
std::vector<double> _tr_start_x;
std::vector<double> _tr_start_y;
std::vector<double> _tr_start_z;
Expand Down
Loading