# File: partialfitter.mac # Brief: This uses the partialFitter. # # Date: 2022-03-22 # Contact: W. Parker, # Revisions: # # This macro runs the partialFitter on simulated partial fill data. # It's based on the scintFitter rat test macro # Load a partial fill phase geometry /rat/db/set DETECTOR geo_file "geo/snoplus_partial.geo" /rat/db/set GEO[inner_av] material_top "labppo_0p5_scintillator" # Measured 0.53 g/L PPO /rat/db/set GEO[inner_av] split_z -5950 /rat/physics_list/OmitMuonicProcesses true /rat/physics_list/OmitHadronicProcesses true # Choose an Nhit threshold #/rat/db/set DAQ_RUN_LEVEL n100_hi_thresh 20.0 #/rat/db/set DAQ_RUN_LEVEL n100_med_thresh 18.0 #/rat/db/set DAQ_RUN_LEVEL n100_lo_thresh 16.0 #/rat/db/set DAQ_RUN_LEVEL n20_thresh 20.0 #/rat/db/set DAQ_RUN_LEVEL n20_lb_thresh 20.0 #/rat/db/set DAQ_RUN_LEVEL eshi_thresh 160.0 #/rat/db/set DAQ_RUN_LEVEL eslo_thresh 120.0 /rat/db/set MC thin_factor 2.5 /run/initialize /rat/proc frontend /rat/proc trigger /rat/proc eventbuilder /rat/proc count /rat/procset update 10 /rat/proc calibratePMT /rat/proc partialFitter /rat/proc outroot /rat/procset file "partialfitter.root" # Generate 3 classes of events: # monoenergetic at (0,0,1)m # monoenergetic distributed in inner_av_top # energy range at (0,0,1)m /generator/add combo gun:point /generator/vtx/set e- 0 0 0 1.0 /generator/pos/set 0 0 1000 /generator/add combo gun:fill /generator/vtx/set e- 0 0 0 1.0 /generator/pos/set 0 0 1000 /generator/add combo gun2:point /generator/vtx/set e- 0 0 0 0 1.0 3.5 /generator/pos/set 0 0 1000 /rat/run/start 1000