Skip to content
Snippets Groups Projects
Commit a568d3a6 authored by John Matter's avatar John Matter Committed by Eric Pooser
Browse files

Update coincidence report template and replay scripts (#365) 

* Update coincidence report template and replay scripts

* Fix BCM label typo
parent 06b19cf1
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 1279 additions and 31 deletions
DEF-files/COIN/PRODUCTION/GTR/coin_gtr_histos.def
+ 3
− 2
View file @ a568d3a6
############################
# TARGET
############################
TH2F hgtrz_vs_pgtrz 'HMS z_{tar} vs SHMS z_{tar}; SHMS z_{tar} / 0.1 cm; HMS z_{tar} / 0.1 cm;' P.gtr.z H.gtr.z 300 -15 15 300 -15 15 coincut
TH2F hgtry_vs_pgtry 'HMS y_{tar} vs SHMS y_{tar}; SHMS y_{tar} / 0.1 cm; HMS y_{tar} / 0.1 cm;' P.gtr.y H.gtr.y 300 -15 15 300 -15 15 coincut
TH2F hgtry_vs_pgtry 'HMS y_{tar} vs SHMS y_{tar}; SHMS y_{tar} / 0.1 cm; HMS y_{tar} / 0.1 cm;' P.gtr.y H.gtr.y 200 -15 15 200 -15 15 coincut
TH2F hgtryp_vs_pgtryp 'HMS y'_{tar} vs SHMS y'_{tar}; SHMS y'_{tar} / 0.1 cm; HMS y'_{tar} / 0.1 cm;' P.gtr.ph H.gtr.ph 200 -0.1 0.1 200 -0.1 0.1 coincut
TH2F hgtrxp_vs_pgtrxp 'HMS x'_{tar} vs SHMS x'_{tar}; SHMS x'_{tar} / 0.1 cm; HMS x'_{tar} / 0.1 cm;' P.gtr.th H.gtr.th 200 -0.1 0.1 200 -0.1 0.1 coincut
DEF-files/COIN/PRODUCTION/TRIG/coin_trig_histos.def
+ 1
− 1
View file @ a568d3a6
......@@ -32,6 +32,6 @@ TH1F c_pathcorr 'Cointime pathcorr; Cointime pathcorr (ns); Count / 1 ns;' T.
TH1F h_htrig_tdc 'pTRIG4_ROC1' T.coin.hTRIG1_ROC1_tdcTimeRaw*0.1 1600 -800 800
TH1F h_ptrig_tdc 'pTRIG4_ROC1' T.coin.pTRIG1_ROC1_tdcTimeRaw*0.1 1600 -800 800
TH1F p_ptrig_tdc 'pTRIG1_ROC2' T.coin.pTRIG1_ROC2_tdcTimeRaw*0.1 1600 -800 800
TH1F p_htrig_tdc 'pTRIG1_ROC2' T.coin.hTRIG1_ROC2_tdcTimeRaw*0.1 1600 -800 800&&shms_trigs
TH1F p_htrig_tdc 'pTRIG1_ROC2' T.coin.hTRIG1_ROC2_tdcTimeRaw*0.1 1600 -800 800
TH1F h_ctrig_tdc 'pTRIG1_ROC1-pTRIG4_ROC1' T.coin.pTRIG1_ROC1_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC1_tdcTimeRaw 1600 -800 800
TH1F p_ctrig_tdc 'pTRIG1_ROC2-pTRIG4_ROC2' T.coin.pTRIG1_ROC2_tdcTimeRaw*0.1-T.coin.pTRIG4_ROC2_tdcTimeRaw 1600 -800 800
DEF-files/COIN/PRODUCTION/coin_hElec_pProt_production_cuts.def → DEF-files/COIN/PRODUCTION/coin_production_hElec_pProt_cuts.def
+ 89
− 0
View file @ a568d3a6
#include "DEF-files/HMS/PRODUCTION/hstackana_production_cuts.def"
#include "DEF-files/SHMS/PRODUCTION/pstackana_production_cuts.def"
Block: RawDecode
Pedestal_event g.evtyp == 99
......@@ -17,10 +14,49 @@ 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
Decode_master ALL_events
Block: CoarseTracking
......@@ -43,8 +79,11 @@ coincut ccleantrack && goodctime
tight_coin ccleantrack && tightctime
ran_coincut ccleantrack && ranctime
e_cut_cal H.cal.etottracknorm > 0.6 && H.cal.etottracknorm < 2.0
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
#include "DEF-files/HMS/PRODUCTION/hstackana_production_cuts.def"
#include "DEF-files/SHMS/PRODUCTION/pstackana_production_cuts.def"
DEF-files/COIN/PRODUCTION/coin_production_pElec_hProt_cuts.def
+ 44
− 5
View file @ a568d3a6
#include "DEF-files/HMS/PRODUCTION/hstackana_production_cuts.def"
#include "DEF-files/SHMS/PRODUCTION/pstackana_production_cuts.def"
Block: RawDecode
Pedestal_event g.evtyp == 99
......@@ -17,10 +14,49 @@ 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
Decode_master ALL_events
Block: CoarseTracking
......@@ -43,8 +79,11 @@ coincut ccleantrack && goodctime
tight_coin ccleantrack && tightctime
ran_coincut ccleantrack && ranctime
e_cut_cal P.cal.etottracknorm > 0.6 && P.cal.etottracknorm < 2.0
e_cut_cer P.hgcer.npeSum > 5.0 && P.ngcer.npeSum > 15.0
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"
SCRIPTS/COIN/PRODUCTION/replay_production_coin_hElec_pProt.C
+ 5
− 5
View file @ a568d3a6
......@@ -146,9 +146,9 @@ void replay_production_coin_hElec_pProt (Int_t RunNumber = 0, Int_t MaxEvent = 0
THcScalerEvtHandler *hscaler = new THcScalerEvtHandler("H", "Hall C scaler event type 4");
hscaler->AddEvtType(2);
hscaler->AddEvtType(4);
pscaler->AddEvtType(5);
pscaler->AddEvtType(6);
pscaler->AddEvtType(7);
hscaler->AddEvtType(5);
hscaler->AddEvtType(6);
hscaler->AddEvtType(7);
hscaler->AddEvtType(129);
hscaler->SetDelayedType(129);
hscaler->SetUseFirstEvent(kTRUE);
......@@ -238,9 +238,9 @@ void replay_production_coin_hElec_pProt (Int_t RunNumber = 0, Int_t MaxEvent = 0
// Define output ROOT file
analyzer->SetOutFile(ROOTFileName.Data());
// Define DEF-file+
analyzer->SetOdefFile("DEF-files/COIN/PRODUCTION/coin_hElec_pProt_production.def");
analyzer->SetOdefFile("DEF-files/COIN/PRODUCTION/coin_production_hElec_pProt.def");
// Define cuts file
analyzer->SetCutFile("DEF-files/COIN/PRODUCTION/coin_hElec_pProt_production_cuts.def"); // optional
analyzer->SetCutFile("DEF-files/COIN/PRODUCTION/coin_production_hElec_pProt_cuts.def"); // optional
// File to record accounting information for cuts
analyzer->SetSummaryFile(Form("REPORT_OUTPUT/COIN/PRODUCTION/summary_production_%d_%d.report", RunNumber, MaxEvent)); // optional
// Start the actual analysis.
......
SCRIPTS/COIN/PRODUCTION/replay_production_coin_pElec_hProt.C
+ 3
− 3
View file @ a568d3a6
......@@ -146,9 +146,9 @@ void replay_production_coin_pElec_hProt (Int_t RunNumber = 0, Int_t MaxEvent = 0
THcScalerEvtHandler *hscaler = new THcScalerEvtHandler("H", "Hall C scaler event type 4");
hscaler->AddEvtType(2);
hscaler->AddEvtType(4);
pscaler->AddEvtType(5);
pscaler->AddEvtType(6);
pscaler->AddEvtType(7);
hscaler->AddEvtType(5);
hscaler->AddEvtType(6);
hscaler->AddEvtType(7);
hscaler->AddEvtType(129);
hscaler->SetDelayedType(129);
hscaler->SetUseFirstEvent(kTRUE);
......
TEMPLATES/COIN/PRODUCTION/coin_production.template
+ 1134
− 15
View file @ a568d3a6
This diff is collapsed.
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment