Difference between revisions of "FY2022 PID analysis"

Latest revision as of 15:25, 2 November 2023

general software for FY2022 beamtest

real data decoder https://github.com/JixieZhang/beamtest_hallc_decoder
real data decoder with maroc https://github.com/zhaozhiwen/beamtest_hallc_decoder
real data analysis ??? (Ye and Michael can add more)
sim and analysis https://github.com/JeffersonLab/solid_gemc/tree/master/prototype/ec_pro/HallC2022

code and document for 2022 beamtest PID analysis

real data root file at  (Ye and Michael can add more about the root file used)
pi file:
electron file:

real data csv file (converted from root file) at ???

sim root file at /cache/halla/solid/sim/solid_gemc/ec_pro_JLAB_VERSION_2.5/pass1/18deg_reduce/
Pencil Beam: root file?
Full Sim:  root file?
Full sim, extra pi-:  root file?
Background events:  root file?
(Ye can add more about the root file used)

sim csv file (converted from root file) at /work/halla/solid/dwupton/Beamtest/
Pencil Beam: root file?, Pencil_Beam/Sim_CSV/Sim_Pencil_AllEvents_TID1.csv
Full Sim:  root file?, Full_ML/Sim_CSV/FullSim_WithTracking.csv
Full sim, extra pi-:  root file?, Full_ML/Sim_CSV/pim_more_All_new.csv
Background events:  root file?, Full_ML/Sim_CSV/BeamOnTarget.csv

Darren's ML PID code:
(Zhiwen will add more details)

Spencer's traditional PID code:
Pencil Beam code: TS3_Pencil_Beam_PID.ipynb
Full Sim Code: Sim_Full_PID.ipynb

Spencer Opatrny working dir and report
/w/halla-scshelf2102/solid/sopatrny/Spencer_Summer_2023
https://solid.jlab.org/cgi-bin/public/ShowDocument?docid=496

Darren Upton's working dir and report
/work/halla/solid/dwupton/Beamtest
https://solid.jlab.org/cgi-bin/public/ShowDocument?docid=495

=== how to merge background in sim ============

here is an example what Zhiwen think how to merge signal and background in sim
// 1.assume beam current 15ua and 50ns time window. we have 15e-6/1.6e-19*50e-9=4.7e6 beam e- within 50ns time window
// 2.assume the background "beamontarget" sim file is from 1e10 e- with 10x actual target density and it has 1e5 number of events with hits in at least one detector.
// 3.beamontarget file has 1e10/4.7e6 = 2128 number of time window.
// 4. for each signal particle in sim, we should add 1e5/10/2127=5 background event from beamontarget file randomly for every detector channel, then output in csv
// 5.for GEM, we need to use 300ns time window. all other detector use 50ns. so GEM needs to add 6 time more background events for each signal event