Skip to content
Snippets Groups Projects
Commit e13f4027 authored by Mark Jones's avatar Mark Jones
Browse files

Merge conflicts in tshms.param with raster update

parents 95f3fc76 66fa962e
No related branches found
No related tags found
No related merge requests found
; SHMS default calibration parameter files
#include "PARAM/SHMS/NGCER/CALIB/pngcer_calib.param"
#include "PARAM/SHMS/DC/CALIB/pdc_calib.param"
#include "PARAM/SHMS/DC/CALIB/pdc_tzero_per_wire.param"
#include "PARAM/SHMS/HODO/CALIB/phodo_calib.param"
#include "PARAM/SHMS/HGCER/CALIB/phgcer_calib.param"
#include "PARAM/SHMS/AERO/CALIB/paero_calib.param"
#include "PARAM/SHMS/CAL/CALIB/pcal_calib.param"
0-99999 0-99999
g_ctp_parm_filename ="DBASE/SHMS/GEN/general.param" g_ctp_parm_filename ="DBASE/SHMS/GEN/general.param"
g_ctp_kinematics_filename ="DBASE/SHMS/STD/standard.kinematics" g_ctp_kinematics_filename ="DBASE/SHMS/STD/standard.kinematics"
g_ctp_calib_filename ="DBASE/SHMS/GEN/default_det_calib_files.param"
464-464 464-464
g_ctp_calib_filename ="DBASE/SHMS/GEN/KPP_Spring_2017/calib_464.param" g_ctp_calib_filename ="DBASE/SHMS/GEN/KPP_Spring_2017/calib_464.param"
......
...@@ -14,6 +14,7 @@ block P.aero.* ...@@ -14,6 +14,7 @@ block P.aero.*
block P.cal.* block P.cal.*
block P.gtr.* block P.gtr.*
block P.kin.* block P.kin.*
block P.rb.*
#********************** #**********************
...@@ -396,7 +397,6 @@ TH1F phodo_beta 'SHMS Hodo Beta; Beta; Number of Entries / 0.02' P.hod.beta 100 ...@@ -396,7 +397,6 @@ TH1F phodo_beta 'SHMS Hodo Beta; Beta; Number of Entries / 0.02' P.hod.beta 100
TH1F phodo_fpHitsTime 'SHMS Hodo Focal Plane Time; Focal Plane Time (ns); Number of Entries / 1 ns' P.hod.fpHitsTime 90 0 90 TH1F phodo_fpHitsTime 'SHMS Hodo Focal Plane Time; Focal Plane Time (ns); Number of Entries / 1 ns' P.hod.fpHitsTime 90 0 90
TH1F phodo_starttime 'SHMS Hodo Start Time; Start Time (ns); Number of Entries / 1ns' P.hod.starttime 90 0 90 TH1F phodo_starttime 'SHMS Hodo Start Time; Start Time (ns); Number of Entries / 1ns' P.hod.starttime 90 0 90
TH1F phodo_good_startime 'SHMS Hodo Good Start Time (ns); Good Start Time (ns); Number of Entries / 1ns' P.hod.goodstarttime 90 0 90 TH1F phodo_good_startime 'SHMS Hodo Good Start Time (ns); Good Start Time (ns); Number of Entries / 1ns' P.hod.goodstarttime 90 0 90
TH2F phodo_xhod_vs_yhod 'SHMS Hodo Projected X vs. Y; Y_hodo(cm) / 1mm; X_hodo(cm) / 1mm' P.hod.y_hodo P.hod.x_hodo 1000 -50 50 1000 -50 50
#********************** #**********************
# Heavy Gas Cherenkov * # Heavy Gas Cherenkov *
...@@ -668,18 +668,18 @@ TH1F ptrig_pdc_ref10 'SHMS DC Reference Time Slot 15 Channel 47; Raw TDC Time (T ...@@ -668,18 +668,18 @@ TH1F ptrig_pdc_ref10 'SHMS DC Reference Time Slot 15 Channel 47; Raw TDC Time (T
# Assume FADC range is set to 1V and the integration is 25 channels # Assume FADC range is set to 1V and the integration is 25 channels
# Offset in the FADC is about 0.122V or 500 channels. # Offset in the FADC is about 0.122V or 500 channels.
formula FRXA_V T.shms.pFRXA_adcPulseIntRaw*(1./4096.)/25. formula FRXA_V P.rb.raster.frxaRawAdc*(1./4096.)/25.
formula FRXB_V T.shms.pFRXB_adcPulseIntRaw*(1./4096.)/25. formula FRXB_V P.rb.raster.frxbRawAdc*(1./4096.)/25.
formula FRYA_V T.shms.pFRYA_adcPulseIntRaw*(1./4096.)/25. formula FRYA_V P.rb.raster.fryaRawAdc*(1./4096.)/25.
formula FRYB_V T.shms.pFRYB_adcPulseIntRaw*(1./4096.)/25. formula FRYB_V P.rb.raster.frybRawAdc*(1./4096.)/25.
TH1F pFRXA 'SHMS FRXA; FRXA RAW ADC (chan); Counts' T.shms.pFRXA_adcPulseIntRaw 1000 0 100000 TH1F pFRXA 'SHMS FRXA; FRXA RAW ADC (chan); Counts' P.rb.raster.frxaRawAdc 1000 0 100000
TH1F pFRXA_V 'SHMS FRXA; FRXA RAW ADC (Volts); Counts' FRXA_V 1000 0 1.0 TH1F pFRXA_V 'SHMS FRXA; FRXA RAW ADC (Volts); Counts' FRXA_V 1000 0 1.0
TH1F pFRXB 'SHMS FRXB; FRXB RAW ADC (chan); Counts' T.shms.pFRXB_adcPulseIntRaw 1000 0 100000 TH1F pFRXB 'SHMS FRXB; FRXB RAW ADC (chan); Counts' P.rb.raster.frxbRawAdc 1000 0 100000
TH1F pFRXB_V 'SHMS FRXB; FRXB RAW ADC (Volts); Counts' FRXB_V 1000 0 1.0 TH1F pFRXB_V 'SHMS FRXB; FRXB RAW ADC (Volts); Counts' FRXB_V 1000 0 1.0
TH1F pFRYA 'SHMS FRYA; FRYA RAW ADC (chan); Counts' T.shms.pFRYA_adcPulseIntRaw 1000 0 100000 TH1F pFRYA 'SHMS FRYA; FRYA RAW ADC (chan); Counts' P.rb.raster.fryaRawAdc 1000 0 100000
TH1F pFRYA_V 'SHMS FRYA; FRYA RAW ADC (Volts); Counts' FRYA_V 1000 0 1.0 TH1F pFRYA_V 'SHMS FRYA; FRYA RAW ADC (Volts); Counts' FRYA_V 1000 0 1.0
TH1F pFRYB 'SHMS FRYB; FRYB RAW ADC (chan); Counts' T.shms.pFRYB_adcPulseIntRaw 1000 0 100000 TH1F pFRYB 'SHMS FRYB; FRYB RAW ADC (chan); Counts' P.rb.raster.frybRawAdc 1000 0 100000
TH1F pFRYB_V 'SHMS FRYB; FRYB RAW ADC (Volts); Counts' FRYB_V 1000 0 1.0 TH1F pFRYB_V 'SHMS FRYB; FRYB RAW ADC (Volts); Counts' FRYB_V 1000 0 1.0
TH2F pFRA_XvsY 'HMSFRA_XvsY; FRXA RAW ADC (Volts); FRYA RAW ADC (Volts)' FRXA_V FRYA_V 1000 0 1.0 1000 0 1.0 TH2F pFRA_XvsY 'HMSFRA_XvsY; FRXA RAW ADC (Volts); FRYA RAW ADC (Volts)' FRXA_V FRYA_V 1000 0 1.0 1000 0 1.0
TH2F pFRB_XvsY 'HMSFRB_XvsY; FRXB RAW ADC (Volts); FRYB RAW ADC (Volts)' FRXB_V FRYB_V 1000 0 1.0 1000 0 1.0 TH2F pFRB_XvsY 'HMSFRB_XvsY; FRXB RAW ADC (Volts); FRYB RAW ADC (Volts)' FRXB_V FRYB_V 1000 0 1.0 1000 0 1.0
...@@ -711,4 +711,4 @@ iunser ...@@ -711,4 +711,4 @@ iunser
itov3out itov3out
itov4out itov4out
end epics end epics
\ No newline at end of file
!HMS: !HMS:
! HRASTER_ID=16 ADC ! HRASTER_ID=18 ADC
! channel is physics channel in ADC
! Plane = 1 is not used
! sig = 0 means data in Positive ADC
! bar will be a counter
! Plane = 1 >> ADC Only detector=18
detector=16
ROC=1 ROC=1
......
!HMS:
! THMS_ID=10 ADC,TDC
! Plane = 1 >> ADC Only
! Plane = 2 >> TDC Only
detector=10
ROC=1
!F250 module
SLOT=18
! chan plane bar sig
12, 1, 1, 0! FR-XA
13, 1, 2, 0! FR-XB
14, 1, 3, 0! FR-YA
15, 1, 4, 0! FR-YB
!SHMS: !SHMS:
! PRASTER_ID=27 ADC ! PRASTER_ID=28 ADC
! channel is physics channel in ADC
! Plane = 1 is not used
! sig = 0 means data in Positive ADC
! bar will be a counter
! Plane = 1 >> ADC Only detector=28
detector=27
ROC=2 ROC=2
...@@ -15,3 +17,4 @@ SLOT=14 ...@@ -15,3 +17,4 @@ SLOT=14
14, 1, 3, 0! FR-YA 14, 1, 3, 0! FR-YA
15, 1, 4, 0! FR-YB 15, 1, 4, 0! FR-YB
!SHMS:
! TSHMS_ID=20 ADC,TDC
! Plane = 1 >> ADC Only
! Plane = 2 >> TDC Only
detector=20
ROC=2
!F250 module
SLOT=14
! chan plane bar sig
12, 1, 1, 0! FR-XA
13, 1, 2, 0! FR-XB
14, 1, 3, 0! FR-YA
15, 1, 4, 0! FR-YB
...@@ -5,7 +5,17 @@ ...@@ -5,7 +5,17 @@
! PHGCER_ID=24 :: ADC ! PHGCER_ID=24 :: ADC
! PAERO_ID=25 :: ,ADC+,ADC- ! PAERO_ID=25 :: ,ADC+,ADC-
! PCAL_ID=26 :: ADC ! PCAL_ID=26 :: ADC
! PRASTER_ID=28 :: ADC
DETECTOR=28
ROC=2
! chan,plane,counter,sig=0 means Pos ADC
SLOT=14
12, 1, 1, 0 ! FR-XA
13, 1, 2, 0 ! FR-XB
14, 1, 3, 0 ! FR-YA
15, 1, 4, 0 ! FR-YB
DETECTOR=20 DETECTOR=20
...@@ -15,10 +25,6 @@ SLOT=14 ...@@ -15,10 +25,6 @@ SLOT=14
0, 1, 1, 0 ! pAERSUM 0, 1, 1, 0 ! pAERSUM
1, 1, 2, 0 ! pHGCERSUM 1, 1, 2, 0 ! pHGCERSUM
2, 1, 3, 0 ! pNGCERSUM 2, 1, 3, 0 ! pNGCERSUM
12, 1, 4, 0 ! FR-XA
13, 1, 5, 0 ! FR-XB
14, 1, 6, 0 ! FR-YA
15, 1, 7, 0 ! FR-YB
SLOT=20 SLOT=20
15, 2, 1, 1 ! pT1 15, 2, 1, 1 ! pT1
......
...@@ -19,8 +19,8 @@ gbpm_kappa = 2.00 ,1.95,1.85 ; modified mkj 4/9/03 ...@@ -19,8 +19,8 @@ gbpm_kappa = 2.00 ,1.95,1.85 ; modified mkj 4/9/03
; gbpm_x_off = -0.011,-0.004,0. ; in cm: survey according to Paul's note ; gbpm_x_off = -0.011,-0.004,0. ; in cm: survey according to Paul's note
; gbpm_y_off = +0.052,+0.056,0. ; gbpm_y_off = +0.052,+0.056,0.
; From Dahlberg survey DT_C853 Apr 7, 2003 ; From Dahlberg survey DT_C853 Apr 7, 2003
gbpm_x_off = 0.017+.412,0.061-.353,0.037 gbpm_x_off = 0.017+.412-.094 , 0.061-.353-.055 ,0.037-0.09
gbpm_y_off = 0.037+.381,-0.04-.211,0.028 gbpm_y_off = 0.037+.381+.003 ,-0.040-.211-.033 ,0.028+.219
; average beam positions: only used if 'guse_bpm_in_recon=0' ; average beam positions: only used if 'guse_bpm_in_recon=0'
; gbeam_xoff = +0.18 ; gbeam_xoff = +0.18
...@@ -31,7 +31,7 @@ gbpm_y_off = 0.037+.381,-0.04-.211,0.028 ...@@ -31,7 +31,7 @@ gbpm_y_off = 0.037+.381,-0.04-.211,0.028
; spectrometers would like to see this positions (from optics studies) ; spectrometers would like to see this positions (from optics studies)
; we treat SOS and HMS the same, however, we may change this if necessary ; we treat SOS and HMS the same, however, we may change this if necessary
; gspec_xoff = +0.18 ; gspec_xoff = +0.18
gspec_xoff = +0.00 gspec_xoff = +0.00
gspec_xpoff = +0.00 gspec_xpoff = +0.00
gspec_yoff = +0.00 gspec_yoff = +0.00
...@@ -55,23 +55,31 @@ gbpm_y_off = 0.037+.381,-0.04-.211,0.028 ...@@ -55,23 +55,31 @@ gbpm_y_off = 0.037+.381,-0.04-.211,0.028
gusefr = 1 ; if 1 correct for FRY in reconstruction gusefr = 1 ; if 1 correct for FRY in reconstruction
guse_frdefault = 1 ; if 1 do no phase correction (default) guse_frdefault = 1 ; if 1 do no phase correction (default)
;if 0 apply phase correction gfr_cal_mom = 6.4 ; = beam momnetum during calibration run
; gfr_cal_mom = 3.245 ; = beam momnetum during calibration run
; gfrx_adcpercm = 11328. ; = FR channels per cm deflection on target ; gfrx_adcpercm = 11328. ; = FR channels per cm deflection on target
; gfry_adcpercm = 10708. ; from run 9981, August 25, 1996, book X-22 ; gfry_adcpercm = 10708. ; from run 9981, August 25, 1996, book X-22
gfr_cal_mom = 2.038 ; = beam momnetum during calibration run ;gfr_cal_mom = 2.038 ; = beam momnetum during calibration run
gfrx_adcpercm = 4364.7 ; = FR channels per cm deflection on target ;gfrx_adcpercm = 4364.7 ; = FR channels per cm deflection on target
gfry_adcpercm = 5471.9 ; from harp scan (H00A), June 19, 2004 ;gfry_adcpercm = 5471.9 ; from harp scan (H00A), June 19, 2004
gfrxa_adcpercm = 97666.7;
gfrxb_adcpercm = 97333.3;
gfrya_adcpercm = 113667;
gfryb_adcpercm = 114333;
gfrx_dphase = 0.0 ; phase shift (never measured for x so far) gfrxa_adc_zero_offset = 67800;
gFrx_synccut = 0. ; should be zero if applied on centered distr. gfrxb_adc_zero_offset = 69050;
gfrya_adc_zero_offset = 67700;
gfryb_adc_zero_offset = 67000;
gfrx_adcmax = 1000 ; ADC amplitude in channels. gfrx_adcmax = 1000 ; ADC amplitude in channels.
gfrx_maxsize = 0.1 ; fast raster amplitude in centimeter. gfrx_maxsize = 0.1 ; fast raster amplitude in centimeter.
gfry_dphase = 5.8 ; phase shift
gFry_synccut = 0. ; should be zero if applied on centered distr.
gfry_adcmax = 1000 ; ADC amplitude in channels. gfry_adcmax = 1000 ; ADC amplitude in channels.
gfry_maxsize = 0.1 ; fast raster amplitude in centimeter. gfry_maxsize = 0.1 ; fast raster amplitude in centimeter.
...@@ -80,8 +88,8 @@ gbpm_y_off = 0.037+.381,-0.04-.211,0.028 ...@@ -80,8 +88,8 @@ gbpm_y_off = 0.037+.381,-0.04-.211,0.028
; was determined to be 5.8 degree. ; was determined to be 5.8 degree.
; positions of FR magnets relative to target ; positions of FR magnets relative to target
gfrx_dist = 2111 ; cm gfrx_dist =1375 ; cm
gfry_dist = 2071 ; cm gfry_dist = 1337 ; cm
......
t_shms_numAdc = 7 t_shms_numAdc = 3
t_shms_numTdc = 23 t_shms_numTdc = 23
t_shms_tdcoffset = -300. t_shms_tdcoffset = -300.
t_shms_tdcchanperns = 0.1 t_shms_tdcchanperns = 0.1
; bar num: 1 2 3 4 5 6 7 ; bar num: 1 2 3
t_shms_adcNames = "pAERSUM pHGCSUM pNGCSUM pFRXA pFRYA pFRXB pFRYB" t_shms_adcNames = "pAERSUM pHGCSUM pNGCSUM"
; bar num: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 ; bar num: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
t_shms_tdcNames = "pT1 pT2 p1X p1Y p2X p2Y p1T p2T pT3 pAERSUM pHGCSUM pNGCSUM pDCREF1 pDCREF2 pDCREF3 pDCREF4 pDCREF5 pDCREF6 pDCREF7 pDCREF8 pDCREF9 pDCREF10 pHODO_EDTM" t_shms_tdcNames = "pT1 pT2 p1X p1Y p2X p2Y p1T p2T pT3 pAERSUM pHGCSUM pNGCSUM pDCREF1 pDCREF2 pDCREF3 pDCREF4 pDCREF5 pDCREF6 pDCREF7 pDCREF8 pDCREF9 pDCREF10 pHODO_EDTM"
\ No newline at end of file
...@@ -39,6 +39,8 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) { ...@@ -39,6 +39,8 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
// Add trigger apparatus // Add trigger apparatus
THaApparatus* TRG = new THcTrigApp("T", "TRG"); THaApparatus* TRG = new THcTrigApp("T", "TRG");
gHaApps->Add(TRG); gHaApps->Add(TRG);
// Add trigger detector to trigger apparatus // Add trigger detector to trigger apparatus
THcTrigDet* shms = new THcTrigDet("shms", "SHMS Trigger Information"); THcTrigDet* shms = new THcTrigDet("shms", "SHMS Trigger Information");
TRG->AddDetector(shms); TRG->AddDetector(shms);
...@@ -58,6 +60,7 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) { ...@@ -58,6 +60,7 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
// Add Heavy Gas Cherenkov to SHMS apparatus // Add Heavy Gas Cherenkov to SHMS apparatus
THcCherenkov* hgcer = new THcCherenkov("hgcer", "Heavy Gas Cherenkov"); THcCherenkov* hgcer = new THcCherenkov("hgcer", "Heavy Gas Cherenkov");
SHMS->AddDetector(hgcer); SHMS->AddDetector(hgcer);
// Add Aerogel Cherenkov to SHMS apparatus // Add Aerogel Cherenkov to SHMS apparatus
THcAerogel* aero = new THcAerogel("aero", "Aerogel"); THcAerogel* aero = new THcAerogel("aero", "Aerogel");
SHMS->AddDetector(aero); SHMS->AddDetector(aero);
...@@ -68,13 +71,15 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) { ...@@ -68,13 +71,15 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
// Include golden track information // Include golden track information
THaGoldenTrack* gtr = new THaGoldenTrack("P.gtr", "SHMS Golden Track", "P"); THaGoldenTrack* gtr = new THaGoldenTrack("P.gtr", "SHMS Golden Track", "P");
gHaPhysics->Add(gtr); gHaPhysics->Add(gtr);
// Add Ideal Beam Apparatus // Add Rastered Beam Apparatus
THaApparatus* beam = new THaIdealBeam("IB", "Ideal Beamline"); THaApparatus* beam = new THcRasteredBeam("P.rb", "Rastered Beamline");
gHaApps->Add(beam); gHaApps->Add(beam);
// Add Physics Module to calculate primary (scattered beam - usually electrons) kinematics // Add Physics Module to calculate primary (scattered beam - usually electrons) kinematics
THcPrimaryKine* kin = new THcPrimaryKine("P.kin", "SHMS Single Arm Kinematics", "P", "IB"); THcPrimaryKine* kin = new THcPrimaryKine("P.kin", "SHMS Single Arm Kinematics", "P", "P.rb");
gHaPhysics->Add(kin); gHaPhysics->Add(kin);
// Add event handler for prestart event 125. // Add event handler for prestart event 125.
THcConfigEvtHandler* ev125 = new THcConfigEvtHandler("HC", "Config Event type 125"); THcConfigEvtHandler* ev125 = new THcConfigEvtHandler("HC", "Config Event type 125");
gHaEvtHandlers->Add(ev125); gHaEvtHandlers->Add(ev125);
......
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