Skip to content
Snippets Groups Projects
Unverified Commit fc1b695f authored by Eric Pooser's avatar Eric Pooser Committed by GitHub
Browse files

Cut files restructure (#391)

* Remove ROOT file

* Utilize #include directives in SHMS

* Update SHMS replay scripts accordingly

* Remove trailing white space

* Do the same for the HMS

* Replecate changes for coincidence replays
parent 8c990e0d
No related branches found
No related tags found
No related merge requests found
Showing
with 181 additions and 377 deletions
Block: RawDecode
Pedestal_event g.evtyp == 99
SHMS_event g.evtyp == 1
HMS_event g.evtyp == 2
BOTH_SINGLES_event g.evtyp == 3
COIN_ONLY_event g.evtyp == 4
SHMS_COIN_event g.evtyp == 5
HMS_COIN_event g.evtyp == 6
COIN_ALL_event g.evtyp == 7
ALL_events SHMS_event || HMS_event || BOTH_SINGLES_event || COIN_ONLY_event || COIN_ALL_event
RawDecode_master ALL_events
Block: Decode
pTRIG1_ROC2 T.coin.pTRIG1_ROC2_tdcTimeRaw > 0
......@@ -57,36 +43,4 @@ pcut_good_S1_time T.coin.p1T_tdcTime > 241 && T.coin.p1T_tdcTime < 243
pcut_good_S2X_time T.coin.p2X_tdcTime > 180 && T.coin.p2X_tdcTime < 220
pcut_good_S1_S2X_time pcut_good_S1_time && pcut_good_S2X_time
ALL_COIN_NO_EDTM T.coin.hEDTM_tdcTime==0.0&&T.coin.pEDTM_tdcTime==0.0&&ALL_events
Decode_master ALL_COIN_NO_EDTM
Block: CoarseTracking
CoarseTracking_master ALL_COIN_NO_EDTM
Block: CoarseReconstruct
CoarseReconstruct_master ALL_COIN_NO_EDTM
Block: Tracking
Tracking_master ALL_COIN_NO_EDTM
Block: Reconstruct
Reconstruct_master ALL_COIN_NO_EDTM
hcleantrack H.gtr.index > -1
pcleantrack P.gtr.index > -1
ccleantrack hcleantrack && pcleantrack
coincut ccleantrack && goodctime
tight_coin ccleantrack && tightctime
ran_coincut ccleantrack && ranctime
e_cut_cal P.cal.etottracknorm > 0.6 && P.cal.etottracknorm < 1.6
e_cut_cer P.hgcer.npeSum > 0.5
select_e_cut e_cut_cal && e_cut_cer
betacut_proton H.gtr.beta > 0.6 && H.gtr.beta < 0.8
coincut_ep ccleantrack && select_e_cut && betacut_proton
#include "DEF-files/HMS/PRODUCTION/hstackana_production_cuts.def"
#include "DEF-files/SHMS/PRODUCTION/pstackana_production_cuts.def"
ALL_COIN_NO_EDTM T.coin.hEDTM_tdcTime==0.0&&T.coin.pEDTM_tdcTime==0.0&&ALL_events
\ No newline at end of file
Block: RawDecode
#include "DEF-files/COIN/PRODUCTION/CUTS/coin_rawdecode_cuts.def"
RawDecode_master ALL_events
Block: Decode
#include "DEF-files/COIN/PRODUCTION/CUTS/coin_decode_cuts.def"
Decode_master ALL_COIN_NO_EDTM
Block: CoarseTracking
#include "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_coarsetracking_cuts.def"
#include "DEF-files/HMS/PRODUCTION/CUTS/hstackana_coarsetracking_cuts.def"
CoarseTracking_master ALL_COIN_NO_EDTM
Block: CoarseReconstruct
CoarseReconstruct_master ALL_COIN_NO_EDTM
Block: Tracking
Tracking_master ALL_COIN_NO_EDTM
Block: Reconstruct
#include "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_reconstruct_cuts.def"
#include "DEF-files/HMS/PRODUCTION/CUTS/hstackana_reconstruct_cuts.def"
#include "DEF-files/COIN/PRODUCTION/CUTS/coin_reconstruct_cuts.def"
Reconstruct_master ALL_COIN_NO_EDTM
\ No newline at end of file
Block: RawDecode
Pedestal_event g.evtyp == 99
SHMS_event g.evtyp == 1
HMS_event g.evtyp == 2
BOTH_SINGLES_event g.evtyp == 3
COIN_ONLY_event g.evtyp == 4
SHMS_COIN_event g.evtyp == 5
HMS_COIN_event g.evtyp == 6
COIN_ALL_event g.evtyp == 7
ALL_HMS_events HMS_event || COIN_ONLY_event || HMS_COIN_event || COIN_ALL_event
ALL_SHMS_events SHMS_event || COIN_ONLY_event || SHMS_COIN_event || COIN_ALL_event
ALL_events SHMS_event || HMS_event || BOTH_SINGLES_event || COIN_ONLY_event || COIN_ALL_event
\ No newline at end of file
Block: Reconstruct
hcleantrack H.gtr.index > -1
pcleantrack P.gtr.index > -1
ccleantrack hcleantrack && pcleantrack
coincut ccleantrack && goodctime
tight_coin ccleantrack && tightctime
ran_coincut ccleantrack && ranctime
he_cut_cal H.cal.etottracknorm > 0.6 && H.cal.etottracknorm < 1.6
he_cut_cer H.cer.npeSum > 0.5
hselect_e_cut he_cut_cal && he_cut_cer
hbetacut_proton P.gtr.beta > 0.6 && P.gtr.beta < 0.8
hcoincut_ep ccleantrack && hselect_e_cut && hbetacut_proton
pe_cut_cal P.cal.etottracknorm > 0.6 && P.cal.etottracknorm < 1.6
pe_cut_cer P.hgcer.npeSum > 0.5
pselect_e_cut pe_cut_cal && pe_cut_cer
pbetacut_proton H.gtr.beta > 0.6 && H.gtr.beta < 0.8
pcoincut_ep ccleantrack && pselect_e_cut && pbetacut_proton
\ No newline at end of file
......@@ -6,9 +6,9 @@
formula mmiss P.kin.secondary.emiss_nuc*P.kin.secondary.emiss_nuc-P.kin.secondary.pmiss*P.kin.secondary.pmiss
TH1F pkin_mmiss 'SHMS M^{2}_{miss} (ccleantrack); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 ccleantrack
TH1F pkin_mmiss_coincut 'SHMS M^{2}_{miss} (coincut); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 coincut
TH1F pkin_mmiss_coincut_ep 'SHMS M^{2}_{miss} (coincut_ep); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 coincut_ep
TH1F pkin_mmiss_tightcut 'SHMS M^{2}_{miss} (tight_coin); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 tight_coin&&select_e_cut&&betacut_proton
TH1F pkin_mmiss_ran_coincut 'SHMS M^{2}_{miss} (ran_coincut); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 ran_coincut&&select_e_cut&&betacut_proton
TH1F pkin_mmiss_coincut_ep 'SHMS M^{2}_{miss} (coincut_ep); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 hcoincut_ep
TH1F pkin_mmiss_tightcut 'SHMS M^{2}_{miss} (tight_coin); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 tight_coin&&hselect_e_cut&&hbetacut_proton
TH1F pkin_mmiss_ran_coincut 'SHMS M^{2}_{miss} (ran_coincut); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 ran_coincut&&hselect_e_cut&&hbetacut_proton
############################
......@@ -16,17 +16,17 @@ TH1F pkin_mmiss_ran_coincut 'SHMS M^{2}_{miss} (ran_coincut); M^{2}_{miss} (GeV
############################
TH1F pkin_pmiss 'SHMS p_{miss} Magnitude (ccleantrack); p_{miss} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss 200 -2 2 ccleantrack
TH1F pkin_pmiss_coincut 'SHMS p_{miss} Magnitude (coincut); p_{miss} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss 200 -2 2 coincut
TH1F pkin_pmiss_pid 'SHMS p_{miss} Magnitude (coincut_ep); p_{miss} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss 200 -2 2 coincut_ep
TH1F pkin_pmiss_x 'SHMS p_{miss,x} (coincut_ep); p_{miss,x} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss_x 200 -2 2 coincut_ep
TH1F pkin_pmiss_y 'SHMS p_{miss,y} (coincut_ep); p_{miss,y} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss_y 200 -2 2 coincut_ep
TH1F pkin_pmiss_z 'SHMS p_{miss,z} (coincut_ep); p_{miss,z} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss_z 200 -2 2 coincut_ep
TH1F pkin_pmiss_pid 'SHMS p_{miss} Magnitude (coincut_ep); p_{miss} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss 200 -2 2 hcoincut_ep
TH1F pkin_pmiss_x 'SHMS p_{miss,x} (coincut_ep); p_{miss,x} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss_x 200 -2 2 hcoincut_ep
TH1F pkin_pmiss_y 'SHMS p_{miss,y} (coincut_ep); p_{miss,y} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss_y 200 -2 2 hcoincut_ep
TH1F pkin_pmiss_z 'SHMS p_{miss,z} (coincut_ep); p_{miss,z} (GeV); Count / 0.02 GeV;' P.kin.secondary.pmiss_z 200 -2 2 hcoincut_ep
############################
# MISSING ENERGY
############################
TH1F pkin_Emiss 'SHMS E_{miss} (ccleantrack); E_{miss} (GeV); Count / 0.01 GeV;' P.kin.secondary.emiss_nuc 400 -1 3 ccleantrack
TH1F pkin_Emiss_coincut 'SHMS E_{miss} (coincut); E_{miss} (GeV); Count / 0.01 GeV;' P.kin.secondary.emiss_nuc 400 -1 3 coincut
TH1F pkin_Emiss_pid 'SHMS E_{miss} (coincut_ep); E_{miss} (GeV); Count / 0.01 GeV;' P.kin.secondary.emiss_nuc 400 -1 3 coincut_ep
TH1F pkin_Emiss_pid 'SHMS E_{miss} (coincut_ep); E_{miss} (GeV); Count / 0.01 GeV;' P.kin.secondary.emiss_nuc 400 -1 3 hcoincut_ep
############################
......
......@@ -6,9 +6,9 @@
formula mmiss H.kin.secondary.emiss_nuc*H.kin.secondary.emiss_nuc-H.kin.secondary.pmiss*H.kin.secondary.pmiss
TH1F hkin_mmiss 'HMS M^{2}_{miss} (ccleantrack); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 ccleantrack
TH1F hkin_mmiss_coincut 'HMS M^{2}_{miss} (coincut); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 coincut
TH1F hkin_mmiss_coincut_ep 'HMS M^{2}_{miss} (coincut_ep); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 coincut_ep
TH1F hkin_mmiss_tightcut 'HMS M^{2}_{miss} (tight_coin); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 tight_coin&&select_e_cut&&betacut_proton
TH1F hkin_mmiss_ran_coincut 'HMS M^{2}_{miss} (ran_coincut); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 ran_coincut&&select_e_cut&&betacut_proton
TH1F hkin_mmiss_coincut_ep 'HMS M^{2}_{miss} (coincut_ep); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 pcoincut_ep
TH1F hkin_mmiss_tightcut 'HMS M^{2}_{miss} (tight_coin); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 tight_coin&&pselect_e_cut&&pbetacut_proton
TH1F hkin_mmiss_ran_coincut 'HMS M^{2}_{miss} (ran_coincut); M^{2}_{miss} (GeV); Count / 0.01 GeV' mmiss 400 -2 2 ran_coincut&&pselect_e_cut&&pbetacut_proton
############################
......@@ -16,17 +16,17 @@ TH1F hkin_mmiss_ran_coincut 'HMS M^{2}_{miss} (ran_coincut); M^{2}_{miss} (GeV)
############################
TH1F hkin_pmiss 'HMS p_{miss} Magnitude (ccleantrack); p_{miss} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss 200 -2 2 ccleantrack
TH1F hkin_pmiss_coincut 'HMS p_{miss} Magnitude (coincut); p_{miss} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss 200 -2 2 coincut
TH1F hkin_pmiss_pid 'HMS p_{miss} Magnitude (coincut_ep); p_{miss} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss 200 -2 2 coincut_ep
TH1F hkin_pmiss_x 'HMS p_{miss,x} (coincut_ep); p_{miss,x} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss_x 200 -2 2 coincut_ep
TH1F hkin_pmiss_y 'HMS p_{miss,y} (coincut_ep); p_{miss,y} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss_y 200 -2 2 coincut_ep
TH1F hkin_pmiss_z 'HMS p_{miss,z} (coincut_ep); p_{miss,z} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss_z 200 -2 2 coincut_ep
TH1F hkin_pmiss_pid 'HMS p_{miss} Magnitude (coincut_ep); p_{miss} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss 200 -2 2 pcoincut_ep
TH1F hkin_pmiss_x 'HMS p_{miss,x} (coincut_ep); p_{miss,x} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss_x 200 -2 2 pcoincut_ep
TH1F hkin_pmiss_y 'HMS p_{miss,y} (coincut_ep); p_{miss,y} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss_y 200 -2 2 pcoincut_ep
TH1F hkin_pmiss_z 'HMS p_{miss,z} (coincut_ep); p_{miss,z} (GeV); Count / 0.02 GeV;' H.kin.secondary.pmiss_z 200 -2 2 pcoincut_ep
############################
# MISSING ENERGY
############################
TH1F hkin_Emiss 'HMS E_{miss} (ccleantrack); E_{miss} (GeV); Count / 0.01 GeV;' H.kin.secondary.emiss_nuc 400 -1 3 ccleantrack
TH1F hkin_Emiss_coincut 'HMS E_{miss} (coincut); E_{miss} (GeV); Count / 0.01 GeV;' H.kin.secondary.emiss_nuc 400 -1 3 coincut
TH1F hkin_Emiss_pid 'HMS E_{miss} (coincut_ep); E_{miss} (GeV); Count / 0.01 GeV;' H.kin.secondary.emiss_nuc 400 -1 3 coincut_ep
TH1F hkin_Emiss_pid 'HMS E_{miss} (coincut_ep); E_{miss} (GeV); Count / 0.01 GeV;' H.kin.secondary.emiss_nuc 400 -1 3 pcoincut_ep
############################
......
......@@ -12,10 +12,10 @@ TH2F pgtr_beta_vs_htrig_cointime 'SHMS Beta vs pTRIG1_ROC1-pTRIG4_ROC1; cointime
############################
TH1F ptrig_cointime_ccleantrack 'SHMS cointime (ccleantrack); SHMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC2_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC2_tdcTimeRaw*0.1 1600 -800 800 ccleantrack
TH1F ptrig_cointime_coincut 'SHMS cointime (coincut); SHMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC2_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC2_tdcTimeRaw*0.1 1600 -800 800 coincut
TH1F ptrig_cointime_coincut_ep 'SHMS cointime (coincut_ep); SHMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC2_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC2_tdcTimeRaw*0.1 1600 -800 800 coincut_ep
TH1F ptrig_cointime_coincut_ep 'SHMS cointime (coincut_ep); SHMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC2_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC2_tdcTimeRaw*0.1 1600 -800 800 pcoincut_ep
TH1F htrig_cointime_ccleantrack 'HMS cointime (ccleantrack); HMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC1_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC1_tdcTimeRaw*0.1 1600 -800 800 ccleantrack
TH1F htrig_cointime_coincut 'HMS cointime (coincut); HMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC1_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC1_tdcTimeRaw*0.1 1600 -800 800 coincut
TH1F htrig_cointime_coincut_ep 'HMS cointime (coincut_ep); HMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC1_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC1_tdcTimeRaw*0.1 1600 -800 800 coincut_ep
TH1F htrig_cointime_coincut_ep 'HMS cointime (coincut_ep); HMS cointime (ns); Count / 1 ns;' T.coin.pTRIG1_ROC1_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC1_tdcTimeRaw*0.1 1600 -800 800 hcoincut_ep
############################
......
Block: RawDecode
Pedestal_event g.evtyp == 99
SHMS_event g.evtyp == 1
HMS_event g.evtyp == 2
BOTH_SINGLES_event g.evtyp == 3
COIN_ONLY_event g.evtyp == 4
SHMS_COIN_event g.evtyp == 5
HMS_COIN_event g.evtyp == 6
COIN_ALL_event g.evtyp == 7
ALL_events SHMS_event || HMS_event || BOTH_SINGLES_event || COIN_ONLY_event || COIN_ALL_event
RawDecode_master ALL_events
Block: Decode
pTRIG1_ROC2 T.coin.pTRIG1_ROC2_tdcTimeRaw > 0
pTRIG2_ROC2 T.coin.pTRIG2_ROC2_tdcTimeRaw > 0
pTRIG3_ROC2 T.coin.pTRIG3_ROC2_tdcTimeRaw > 0
pTRIG4_ROC2 T.coin.pTRIG4_ROC2_tdcTimeRaw > 0
pTRIG5_ROC2 T.coin.pTRIG5_ROC2_tdcTimeRaw > 0
pTRIG6_ROC2 T.coin.pTRIG6_ROC2_tdcTimeRaw > 0
pTRIG1_ROC1 T.coin.pTRIG1_ROC1_tdcTimeRaw > 0
pTRIG2_ROC1 T.coin.pTRIG2_ROC1_tdcTimeRaw > 0
pTRIG3_ROC1 T.coin.pTRIG3_ROC1_tdcTimeRaw > 0
pTRIG4_ROC1 T.coin.pTRIG4_ROC1_tdcTimeRaw > 0
pTRIG5_ROC1 T.coin.pTRIG5_ROC1_tdcTimeRaw > 0
pTRIG6_ROC1 T.coin.pTRIG6_ROC1_tdcTimeRaw > 0
hTRIG1_ROC2 T.coin.hTRIG1_ROC2_tdcTimeRaw > 0
hTRIG2_ROC2 T.coin.hTRIG2_ROC2_tdcTimeRaw > 0
hTRIG3_ROC2 T.coin.hTRIG3_ROC2_tdcTimeRaw > 0
hTRIG4_ROC2 T.coin.hTRIG4_ROC2_tdcTimeRaw > 0
hTRIG5_ROC2 T.coin.hTRIG5_ROC2_tdcTimeRaw > 0
hTRIG6_ROC2 T.coin.hTRIG6_ROC2_tdcTimeRaw > 0
hTRIG1_ROC1 T.coin.hTRIG1_ROC1_tdcTimeRaw > 0
hTRIG2_ROC1 T.coin.hTRIG2_ROC1_tdcTimeRaw > 0
hTRIG3_ROC1 T.coin.hTRIG3_ROC1_tdcTimeRaw > 0
hTRIG4_ROC1 T.coin.hTRIG4_ROC1_tdcTimeRaw > 0
hTRIG5_ROC1 T.coin.hTRIG5_ROC1_tdcTimeRaw > 0
hTRIG6_ROC1 T.coin.hTRIG6_ROC1_tdcTimeRaw > 0
shms_trigs SHMS_event
hms_trigs HMS_event
SHMS_edtm_accepted T.coin.pEDTM_tdcTime != 0.0
SHMS_physics_accepted T.coin.pEDTM_tdcTime == 0.0
HMS_edtm_accepted T.coin.hEDTM_tdcTime != 0.0
HMS_physics_accepted T.coin.hEDTM_tdcTime == 0.0
goodctime ALL_events #LOOSE CUT AROUND COINTIME PEAK (3 beam buckets)
tightctime ALL_events #TIGHT CUT AROUND COINTIME PEAK (1 beam buckets)
ranctime ALL_events #CUT AWAY FROM COINTIME PEAK (4 beam buckets)
pcut_good_S1_time T.coin.p1T_tdcTime > 241 && T.coin.p1T_tdcTime < 243
pcut_good_S2X_time T.coin.p2X_tdcTime > 180 && T.coin.p2X_tdcTime < 220
pcut_good_S1_S2X_time pcut_good_S1_time && pcut_good_S2X_time
ALL_COIN_NO_EDTM T.coin.hEDTM_tdcTime==0.0&&T.coin.pEDTM_tdcTime==0.0&&ALL_events
Decode_master ALL_COIN_NO_EDTM
Block: CoarseTracking
#include "DEF-files/SHMS/PRODUCTION/pstackana_coarsetracking_cuts.def"
#include "DEF-files/HMS/PRODUCTION/hstackana_coarsetracking_cuts.def"
CoarseTracking_master ALL_COIN_NO_EDTM
Block: CoarseReconstruct
CoarseReconstruct_master ALL_COIN_NO_EDTM
Block: Tracking
Tracking_master ALL_COIN_NO_EDTM
Block: Reconstruct
#include "DEF-files/SHMS/PRODUCTION/pstackana_reconstruction_cuts.def"
#include "DEF-files/HMS/PRODUCTION/hstackana_reconstruction_cuts.def"
Reconstruct_master ALL_COIN_NO_EDTM
hcleantrack H.gtr.index > -1
pcleantrack P.gtr.index > -1
ccleantrack hcleantrack && pcleantrack
coincut ccleantrack && goodctime
tight_coin ccleantrack && tightctime
ran_coincut ccleantrack && ranctime
e_cut_cal H.cal.etottracknorm > 0.6 && H.cal.etottracknorm < 1.6
e_cut_cer H.cer.npeSum > 0.5
select_e_cut e_cut_cal && e_cut_cer
betacut_proton P.gtr.beta > 0.6 && P.gtr.beta < 0.8
coincut_ep ccleantrack && select_e_cut && betacut_proton
Block: Decode
hcut_TRIG1 T.hms.hTRIG1_tdcTimeRaw > 0
hcut_TRIG2 T.hms.hTRIG2_tdcTimeRaw > 0
hcut_TRIG3 T.hms.hTRIG3_tdcTimeRaw > 0
HMS_Pretrig hcut_TRIG1 || hcut_TRIG2 || hcut_TRIG3
hms_trigs HMS_event
hcut_edtm_accepted T.hms.hEDTM_tdcTime != 0.0
hcut_physics_accepted T.hms.hEDTM_tdcTime == 0.0
ALL_HMS_NO_EDTM ALL_HMS_events && hcut_physics_accepted
\ No newline at end of file
# Report file for HMS stack
Block: RawDecode
#include "DEF-files/HMS/PRODUCTION/CUTS/hstackana_rawdecode_cuts.def"
RawDecode_master 1
Block: Decode
#include "DEF-files/HMS/PRODUCTION/CUTS/hstackana_decode_cuts.def"
Decode_master ALL_HMS_NO_EDTM
Block: CoarseTracking
#include "DEF-files/HMS/PRODUCTION/CUTS/hstackana_coarsetracking_cuts.def"
CoarseTracking_master ALL_HMS_events
Block: CoarseReconstruct
CoarseReconstruct_master ALL_HMS_NO_EDTM
Block: Tracking
Tracking_master ALL_HMS_NO_EDTM
Block: Reconstruct
#include "DEF-files/HMS/PRODUCTION/CUTS/hstackana_reconstruct_cuts.def"
Reconstruct_master ALL_HMS_NO_EDTM
Block: Physics
Physics_master ALL_HMS_NO_EDTM
Block: RawDecode
Pedestal_event g.evtyp == 99
SHMS_event g.evtyp == 1
HMS_event g.evtyp == 2
BOTH_SINGLES_event g.evtyp == 3
COIN_ONLY_event g.evtyp == 4
SHMS_COIN_event g.evtyp == 5
HMS_COIN_event g.evtyp == 6
COIN_ALL_event g.evtyp == 7
ALL_events SHMS_event || HMS_event || BOTH_SINGLES_event || COIN_ONLY_event || COIN_ALL_event
ALL_HMS_events HMS_event || COIN_ONLY_event || HMS_COIN_event || COIN_ALL_event
\ No newline at end of file
# Report file for HMS stack
Block: RawDecode
Pedestal_event g.evtyp == 99
SHMS_event g.evtyp == 1
HMS_event g.evtyp == 2
BOTH_SINGLES_event g.evtyp == 3
COIN_ONLY_event g.evtyp == 4
SHMS_COIN_event g.evtyp == 5
HMS_COIN_event g.evtyp == 6
COIN_ALL_event g.evtyp == 7
ALL_events SHMS_event || HMS_event || BOTH_SINGLES_event || COIN_ONLY_event || COIN_ALL_event
ALL_HMS_events HMS_event || COIN_ONLY_event || HMS_COIN_event || COIN_ALL_event
RawDecode_master 1
Block: Decode
hcut_TRIG1 T.hms.hTRIG1_tdcTimeRaw > 0
hcut_TRIG2 T.hms.hTRIG2_tdcTimeRaw > 0
hcut_TRIG3 T.hms.hTRIG3_tdcTimeRaw > 0
HMS_Pretrig hcut_TRIG1 || hcut_TRIG2 || hcut_TRIG3
hms_trigs HMS_event
hcut_edtm_accepted T.hms.hEDTM_tdcTime != 0.0
hcut_physics_accepted T.hms.hEDTM_tdcTime == 0.0
ALL_HMS_NO_EDTM ALL_HMS_events && hcut_physics_accepted
Decode_master ALL_HMS_NO_EDTM
Block: CoarseTracking
CoarseTracking_master ALL_HMS_events
add cuts to stricter times and apply them to drift distances (HMS 6 GeV Drift Chambers)
hcut_time_cut1 H.dc.1x1.time > 0 && H.dc.1x1.time < 250
hcut_time_cut2 H.dc.1v2.time > 0 && H.dc.1v2.time < 250
hcut_time_cut3 H.dc.1u1.time > 0 && H.dc.1u1.time < 250
hcut_time_cut4 H.dc.1v1.time > 0 && H.dc.1v1.time < 250
hcut_time_cut5 H.dc.1u2.time > 0 && H.dc.1u2.time < 250
hcut_time_cut6 H.dc.1x2.time > 0 && H.dc.1x2.time < 250
hcut_time_cut7 H.dc.2x1.time > 0 && H.dc.2x1.time < 250
hcut_time_cut8 H.dc.2v2.time > 0 && H.dc.1v2.time < 250
hcut_time_cut9 H.dc.2u1.time > 0 && H.dc.2u1.time < 250
hcut_time_cut10 H.dc.2v1.time > 0 && H.dc.2v1.time < 250
hcut_time_cut11 H.dc.2u2.time > 0 && H.dc.2u2.time < 250
hcut_time_cut12 H.dc.2x2.time > 0 && H.dc.2x2.time < 250
hmsDC1Planes_large (H.dc.1x1.nhit+H.dc.1v2.nhit+H.dc.1u1.nhit+H.dc.1v1.nhit+H.dc.1x2.nhit+H.dc.1u2.nhit) >20
hmsDC2Planes_large (H.dc.2x1.nhit+H.dc.2v2.nhit+H.dc.2u1.nhit+H.dc.2v1.nhit+H.dc.2x2.nhit+H.dc.2u2.nhit) >20
hmsDCany_large hmsDC1Planes_large || hmsDC2Planes_large
hmsDC1_1hit_x1 H.dc.1x1.nhit == 1
hmsDC1_1hit_v2 H.dc.1v2.nhit == 1
hmsDC1_1hit_u1 H.dc.1u1.nhit == 1
hmsDC1_1hit_v1 H.dc.1v1.nhit == 1
hmsDC1_1hit_u2 H.dc.1u2.nhit == 1
hmsDC1_1hit_x2 H.dc.1x2.nhit == 1
hmsDC2_1hit_x1 H.dc.2x1.nhit == 1
hmsDC2_1hit_v2 H.dc.2v2.nhit == 1
hmsDC2_1hit_u1 H.dc.2u1.nhit == 1
hmsDC2_1hit_v1 H.dc.2v1.nhit == 1
hmsDC2_1hit_u2 H.dc.2u2.nhit == 1
hmsDC2_1hit_x2 H.dc.2x2.nhit == 1
h1hit1 H.dc.1x1.nhit >= 1
h1hit2 H.dc.1v2.nhit >= 1
h1hit3 H.dc.1u1.nhit >= 1
h1hit4 H.dc.1v1.nhit >= 1
h1hit5 H.dc.1u2.nhit >= 1
h1hit6 H.dc.1x2.nhit >= 1
h2hit1 H.dc.2x1.nhit >= 1
h2hit2 H.dc.2v2.nhit >= 1
h2hit3 H.dc.2u1.nhit >= 1
h2hit4 H.dc.2v1.nhit >= 1
h2hit5 H.dc.2u2.nhit >= 1
h2hit6 H.dc.2x2.nhit >= 1
hmsDC1Planes6hits (H.dc.1x1.nhit+H.dc.1v2.nhit+H.dc.1u1.nhit+H.dc.1v1.nhit+H.dc.1x2.nhit+H.dc.1u2.nhit) ==6
hmsDC2Planes6hits (H.dc.2x1.nhit+H.dc.2v2.nhit+H.dc.2u1.nhit+H.dc.2v1.nhit+H.dc.2x2.nhit+H.dc.2u2.nhit )==6
hmsDC1_5hits_x1 hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_6hits_x1 h1hit1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_5hits_v2 hmsDC1_1hit_x1&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_6hits_v2 h1hit2&&hmsDC1_1hit_x1&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_5hits_u1 hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_6hits_u1 h1hit3&&hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_5hits_v1 hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_6hits_v1 h1hit4&&hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_x2&&hmsDC1_1hit_u2
hmsDC1_5hits_x2 hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_u2
hmsDC1_6hits_x2 h1hit6&&hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_u2
hmsDC1_5hits_u2 hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2
hmsDC1_6hits_u2 h1hit5&&hmsDC1_1hit_x1&&hmsDC1_1hit_v2&&hmsDC1_1hit_u1&&hmsDC1_1hit_v1&&hmsDC1_1hit_x2
hmsDC2_5hits_x1 hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_6hits_x1 h2hit1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_5hits_v2 hmsDC2_1hit_x1&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_6hits_v2 h2hit2&&hmsDC2_1hit_x1&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_5hits_u1 hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_6hits_u1 h2hit3&&hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_5hits_v1 hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_6hits_v1 h2hit4&&hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_x2&&hmsDC2_1hit_u2
hmsDC2_5hits_x2 hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_u2
hmsDC2_6hits_x2 h2hit6&&hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_u2
hmsDC2_5hits_u2 hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2
hmsDC2_6hits_u2 h2hit5&&hmsDC2_1hit_x1&&hmsDC2_1hit_v2&&hmsDC2_1hit_u1&&hmsDC2_1hit_v1&&hmsDC2_1hit_x2
hms1HitsLt H.dc.Ch1.nhit <= H.dc.Ch1.maxhits && g.evtyp==1
hms2HitsLt H.dc.Ch2.nhit <= H.dc.Ch2.maxhits && g.evtyp==1
hmsHitsLt H.dc.Ch1.nhit <= H.dc.Ch1.maxhits && H.dc.Ch2.nhit <= H.dc.Ch2.maxhits && g.evtyp==1
hmsDC1PlanesGT (h1hit1 + h1hit2 + h1hit3 + h1hit4 + h1hit5 + h1hit6 )>=5
hmsDC2PlanesGT (h2hit1 + h2hit2 + h2hit3 + h2hit4 + h2hit5 + h2hit6 )>=5
hmsPlanesGT hmsDC1PlanesGT && hmsDC2PlanesGT
hmsHitsPlanes (H.dc.Ch1.nhit <= H.dc.Ch1.maxhits) && (H.dc.Ch2.nhit <= H.dc.Ch2.maxhits) && hmsPlanesGT
hSpacePoints H.dc.Ch1.spacepoints >= 1 && H.dc.Ch2.spacepoints >=1
hSpacePointsStub H.dc.stubtest==1 && H.dc.Ch1.spacepoints >=1 && H.dc.Ch2.spacepoints >=1
hFoundTrack H.dc.ntrack>0
hStubLT H.dc.stubtest==1
f1HSpacePoints hms1HitsLt && hmsDC1PlanesGT && H.dc.Ch1.spacepoints==0 && g.evtyp==1
f2HSpacePoints hms2HitsLt && hmsDC2PlanesGT && H.dc.Ch2.spacepoints==0 && g.evtyp==1
hTest1 hmsHitsPlanes && (!hSpacePoints)
hTest2 hSpacePoints && (!hStubLT)
Block: CoarseReconstruct
CoarseReconstruct_master ALL_HMS_NO_EDTM
Block: Tracking
Tracking_master ALL_HMS_NO_EDTM
Block: Reconstruct
Reconstruct_master ALL_HMS_NO_EDTM
hcut_cer_pmt1_elec H.cer.npe[0]>.5
hcut_cer_pmt2_elec H.cer.npe[0]>.5
hcut_cer_elec H.cer.npeSum>.5
hcut_cer_pi H.cer.npeSum<.5
HMSScinGood H.hod.goodscinhit == 1
HMSGoodBetanotrk H.hod.betanotrack > 0.8 && H.hod.betanotrack < 1.3
HMSScinShould HMSScinGood && HMSGoodBetanotrk && !hmsDCany_large
HMSScinShoulde HMSScinShould && H.cal.etotnorm > 0.6&& H.cal.etotnorm < 2.0&& H.cer.npeSum > 0.5
HMSScinShouldh HMSScinGood && H.cal.etotnorm <0.6&& H.cal.etotnorm>0.0&& H.cer.npeSum < 0.5
HMSScinDid HMSScinShould && H.dc.ntrack > 0
HMSScinDide HMSScinShoulde && H.dc.ntrack > 0
HMSScinDidh HMSScinShouldh && H.dc.ntrack > 0
hcut_goodHDC1x1 H.dc.1x1.nhit > 0 && H.dc.1x1.nhit < 3
hcut_goodHDC1v2 H.dc.1v2.nhit > 0 && H.dc.1v2.nhit < 3
hcut_goodHDC1u1 H.dc.1u1.nhit > 0 && H.dc.1u1.nhit < 3
hcut_goodHDC1v1 H.dc.1v1.nhit > 0 && H.dc.1v1.nhit < 3
hcut_goodHDC1u2 H.dc.1u2.nhit > 0 && H.dc.1u2.nhit < 3
hcut_goodHDC1x2 H.dc.1x2.nhit > 0 && H.dc.1x2.nhit < 3
hcut_goodHDC2x1 H.dc.2x1.nhit > 0 && H.dc.2x1.nhit < 3
hcut_goodHDC2v2 H.dc.2v2.nhit > 0 && H.dc.2v2.nhit < 3
hcut_goodHDC2u1 H.dc.2u1.nhit > 0 && H.dc.2u1.nhit < 3
hcut_goodHDC2v1 H.dc.2v1.nhit > 0 && H.dc.2v1.nhit < 3
hcut_goodHDC2u2 H.dc.2u2.nhit > 0 && H.dc.2u2.nhit < 3
hcut_goodHDC2x2 H.dc.2x2.nhit > 0 && H.dc.2x2.nhit < 3
hcut_goodHDC1 hcut_goodHDC1x1 && hcut_goodHDC1v2 && hcut_goodHDC1u1 && hcut_goodHDC1v1 && hcut_goodHDC1u2 && hcut_goodHDC1x2
hcut_goodHDC2 hcut_goodHDC2x1 && hcut_goodHDC2v2 && hcut_goodHDC2u1 && hcut_goodHDC2v1 && hcut_goodHDC2u2 && hcut_goodHDC2x2
hcut_bothGood hcut_goodHDC1 && hcut_goodHDC2
hcut_realhdc1x1 hcut_goodHDC1x1 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc1v2 hcut_goodHDC1v2 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc1u1 hcut_goodHDC1u1 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc1v1 hcut_goodHDC1v1 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc1u2 hcut_goodHDC1u2 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc1x2 hcut_goodHDC1x2 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc2x1 hcut_goodHDC2x1 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc2v2 hcut_goodHDC2v2 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc2u1 hcut_goodHDC2u1 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc2v1 hcut_goodHDC2v1 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc2u2 hcut_goodHDC2u2 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_realhdc2x2 hcut_goodHDC2x2 && ((H.dc.Ch1.spacepoints+H.dc.Ch2.spacepoints)>0)
hcut_FoundTrack H.dc.ntrack > 0
hcut_Found1Track H.dc.ntrack == 1
hcut_Found2Track H.dc.ntrack == 2
hcut_Found3Track H.dc.ntrack == 3
hcut_Found4Track H.dc.ntrack == 4
hcut_CleanTrack H.gtr.index > -1
hcut_Clean1Track H.gtr.index == 0
hcut_Clean2Track H.gtr.index == 1
hcut_Clean3Track H.gtr.index == 2
hcut_Clean4Track H.gtr.index == 3
hcut_anys1x H.hod.1x.nhits > 0
hcut_anys1y H.hod.1y.nhits > 0
hcut_anys2x H.hod.2x.nhits > 0
hcut_anys2y H.hod.2y.nhits > 0
hcut_goods1x H.hod.1x.nhits > 0 && H.hod.1x.nhits < 3
hcut_goods1y H.hod.1y.nhits > 0 && H.hod.1y.nhits < 3
hcut_goods1 hcut_goods1x && hcut_goods1y
hcut_goods2x H.hod.2x.nhits > 0 && H.hod.2x.nhits < 3
hcut_goods2y H.hod.2y.nhits > 0 && H.hod.2y.nhits < 3
hcut_goods2 hcut_goods2x && hcut_goods2y
hcut_goods1s2 hcut_goods1 && hcut_goods2
HMS_hcer_track_matched_region_1 H.cer.numTracksMatched[0]>0
HMS_hcer_track_fired_region_1 H.cer.numTracksFired[0]>0
HMS_hcer_track_matched_region_2 H.cer.numTracksMatched[1]>0
HMS_hcer_track_fired_region_2 H.cer.numTracksFired[1]>0
HMS_hcer_track_matched_tot H.cer.totNumTracksMatched>0
HMS_hcer_track_fired_tot H.cer.totNumTracksFired>0
Block: Physics
Physics_master ALL_HMS_NO_EDTM
Block: Decode
pcut_TRIG1 T.shms.pTRIG1_tdcTimeRaw > 0
pcut_TRIG2 T.shms.pTRIG2_tdcTimeRaw > 0
pcut_TRIG3 T.shms.pTRIG3_tdcTimeRaw > 0
SHMS_Pretrig pcut_TRIG1 || pcut_TRIG2 || pcut_TRIG3
shms_trigs SHMS_event
pcut_edtm_accepted T.shms.pEDTM_tdcTime != 0.0
pcut_physics_accepted T.shms.pEDTM_tdcTime == 0.0
pcut_good_S1_time T.shms.p1T_tdcTime > 241 && T.shms.p1T_tdcTime < 243
pcut_good_S2X_time T.shms.p2X_tdcTime > 180 && T.shms.p2X_tdcTime < 220
pcut_good_S1_S2X_time pcut_good_S1_time && pcut_good_S2X_time
ALL_SHMS_NO_EDTM ALL_events && pcut_physics_accepted
\ No newline at end of file
# Report file for SHMS stack
Block: RawDecode
#include "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_rawdecode_cuts.def"
RawDecode_master 1
Block: Decode
#include "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_decode_cuts.def"
Decode_master ALL_SHMS_NO_EDTM
Block: CoarseTracking
#include "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_coarsetracking_cuts.def"
CoarseTracking_master ALL_SHMS_NO_EDTM
Block: CoarseReconstruct
CoarseReconstruct_master ALL_SHMS_NO_EDTM
Block: Tracking
Tracking_master ALL_SHMS_NO_EDTM
Block: Reconstruct
#include "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_reconstruct_cuts.def"
Reconstruct_master ALL_SHMS_NO_EDTM
Block: Physics
Physics_master ALL_SHMS_NO_EDTM
Block: RawDecode
Pedestal_event g.evtyp == 99
SHMS_event g.evtyp == 1
HMS_event g.evtyp == 2
BOTH_SINGLES_event g.evtyp == 3
COIN_ONLY_event g.evtyp == 4
SHMS_COIN_event g.evtyp == 5
HMS_COIN_event g.evtyp == 6
COIN_ALL_event g.evtyp == 7
ALL_events SHMS_event || HMS_event || BOTH_SINGLES_event || COIN_ONLY_event || COIN_ALL_event
ALL_SHMS_events SHMS_event || COIN_ONLY_event || SHMS_COIN_event || COIN_ALL_event
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment