Skip to content
Snippets Groups Projects
Commit e6f32ebc authored by Jure Bericic's avatar Jure Bericic Committed by GitHub
Browse files

Merge pull request #85 from JureBericic/fixes

Minor fixes.
parents 4b6a3ce9 19adf222
No related branches found
No related tags found
No related merge requests found
...@@ -30,21 +30,21 @@ Running code ...@@ -30,21 +30,21 @@ Running code
* Edit the parameters in tofcal.param into either PARAM/HMS/HODO/hhodo.param or PARAM/SHMS/HODO/phodo.param * Edit the parameters in tofcal.param into either PARAM/HMS/HODO/hhodo.param or PARAM/SHMS/HODO/phodo.param
Brief decription of code Brief decription of code
------------------------ ------------------------
* The code fits tcorr = time - offset - path * velocity - adccor / sqrt(ADC) * The code fits tcorr = time - offset - path * velocity - adccor / sqrt(ADC)
* Uses paddle 10 positive PMT on plane 1 as the reference PMT ( offset fixed at 0) * Uses paddle 10 positive PMT on plane 1 as the reference PMT ( offset fixed at 0)
* To get data file name. Reads-in the input file: tofcal.inp * To get data file name. Reads-in the input file: tofcal.inp
* Opens data file * Opens data file
** "0" line separates events * "0" line separates events
** Should have pair of neg and pos PMTs for each paddle * Should have pair of neg and pos PMTs for each paddle
** Format is : iside ipl ipad tdcval_uncorr pathl zcor tdcval_corr adcval * Format is : iside ipl ipad tdcval_uncorr pathl zcor tdcval_corr adcval
*** iside =1,2 for pos and neg PMT * iside = 1,2 for pos and neg PMT
*** ipl = plane number * ipl = plane number
*** tdcval_uncorr = raw time in ns. The program uses this in optimization. * tdcval_uncorr = raw time in ns. The program uses this in optimization.
*** pathl = distance for focal plane to the paddle using the track * pathl = distance for focal plane to the paddle using the track
*** zcor = time (ns) for the particle to go from focal plane to paddle. Correction subtracted from tdcval_uncrr * zcor = time (ns) for the particle to go from focal plane to paddle. Correction subtracted from tdcval_uncrr
*** tdcval_corr = the corrected time in ns. Not used by program. Just gives the initial chi2. * tdcval_corr = the corrected time in ns. Not used by program. Just gives the initial chi2.
*** adcval = ADC value used in fit * adcval = ADC value used in fit
* Loops through data determines the number of hits in each PMT. Will only include PMT if the number of hits > 100. * Loops through data determines the number of hits in each PMT. Will only include PMT if the number of hits > 100.
* Again loops through the data filling array with time difference between all pairs PMTs, pathl and adcval. * Again loops through the data filling array with time difference between all pairs PMTs, pathl and adcval.
* Use the CERNLIB routine deqn to invert the matrix and determine the fit parameters. * Use the CERNLIB routine deqn to invert the matrix and determine the fit parameters.
\ No newline at end of file
...@@ -35,18 +35,18 @@ pdc_wire_counting = 0, 0, 0, 0, 0, 0 ...@@ -35,18 +35,18 @@ pdc_wire_counting = 0, 0, 0, 0, 0, 0
; z position in cm of each plane measured from focal plane. ; z position in cm of each plane measured from focal plane.
; From technical drawings. ; From technical drawings.
pdc_zpos = pdc_1_zpos - 0.6875*cminch pdc_zpos = pdc_1_zpos - 0.68701*cminch
pdc_1_zpos - 0.4375*cminch pdc_1_zpos - 0.43701*cminch
pdc_1_zpos - 0.1875*cminch pdc_1_zpos - 0.18701*cminch
pdc_1_zpos + 0.1875*cminch pdc_1_zpos + 0.18701*cminch
pdc_1_zpos + 0.4375*cminch pdc_1_zpos + 0.43701*cminch
pdc_1_zpos + 0.6875*cminch pdc_1_zpos + 0.68701*cminch
pdc_2_zpos - 0.6875*cminch pdc_2_zpos - 0.68701*cminch
pdc_2_zpos - 0.4375*cminch pdc_2_zpos - 0.43701*cminch
pdc_2_zpos - 0.1875*cminch pdc_2_zpos - 0.18701*cminch
pdc_2_zpos + 0.1875*cminch pdc_2_zpos + 0.18701*cminch
pdc_2_zpos + 0.4375*cminch pdc_2_zpos + 0.43701*cminch
pdc_2_zpos + 0.6875*cminch pdc_2_zpos + 0.68701*cminch
; Alpha angle of wires for each plane in radians. ; Alpha angle of wires for each plane in radians.
; From technical drawings. ; From technical drawings.
......
...@@ -14,7 +14,7 @@ pdc_tdc_max_win = 55000, 55000, 55000, 55000, 55000, 55000 ; TODO - wide limits ...@@ -14,7 +14,7 @@ pdc_tdc_max_win = 55000, 55000, 55000, 55000, 55000, 55000 ; TODO - wide limits
pdc_tdc_time_per_channel = -0.10 pdc_tdc_time_per_channel = -0.10
; Zero time correction for each plane in ns that is added to TDC time. ; Zero time correction for each plane in ns that is added to TDC time.
pdc_plane_time_zero = 1260.0, 1260.0, 1260.0, 1260.0, 1260.0, 1260.0 ; TODO - SOS values for now pdc_plane_time_zero = 1260.0, 1260.0, 1260.0, 1260.0, 1260.0, 1260.0
1260.0, 1260.0, 1260.0, 1260.0, 1260.0, 1260.0 1260.0, 1260.0, 1260.0, 1260.0, 1260.0, 1260.0
; For wire velocity corrections. ; For wire velocity corrections.
......
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