1) some info is written in tuda/Documents/READMEs/readme.MIDASsort
MIDAS Main Menu -> Sorting -> MIDASsort
enter sort.f file in Sortfile window and compile with Compile Sortfile
- example file sort3.f with the variables.dat file that contains gain, offset, calibration and other variables, is attached below.
- the base code which just unpacks data without sorting it, sort.f is also attached
- make sure you have the correct path for the variables.dat file in sort3.f
- Note that this variables.dat file was downloaded on a Mac OSX and which adds extra invisible characters at the end of each line and therefore won't be read in properly in MIDASsort.
for online sorting go to Input Media -> Online
for offline sorting go to Input Media -> Disc, the Disc Info. Browse file, Add it and Save and Quit
hit Start button to run sort file. the output is printed in the terminal running MIDAS-session
in Spectrum Directory Browser or Spectrum Viewer go to Resource -> SortSas the Redisplay. here you'll see all spectra created in sort.f
you can add a spectrum title files like the example given below in the Spectrum Directory Browser via Actions -> Apply title and then Options -> Details
- you can add titles for the Resource -> hist spectra with the another title file
if the spectra don't refresh properly Select all and Actions -> Delete in Spectrum Directory Browser and then Start the sorting program
2) you can view online spectra in Spectrum Directory Browser -> Resource -> hist
remember to Select all and Actions -> zero and deselect if you want to view only the current run, otherwise you are seeing the spectra of incremented runs
3) to get gain, offset and calibration coefficients go to tuda/Documents/READMEs/readme.calibration
you need to have pulsar runs with amplitude varied from 90,000 to 10,000 in step size of 10,000 and the peaks should span the entire spectrum if possible and all contain ~1000 counts
you need 3-alpha calibration source runs
- gain and offsets are calculated for ADCs and offsets only for TDCs since clock is the same and there are offsets due to cable lengths, hence delay times
- Tom did this for us since we couldn't figure out the first step in the readme file using the peakfind program. The correct syntax for future knowledge is:
peakfind -spectrum <spectrum file name> -limit <low limit> <high limit> -fwhm <fwhm> -accept <accept>
4) to view scalars during run go to VME Module Control -> CEN V560 scalar and wait 10 s for refresh rate
5) meaning of c and w variables in sort3.f
c(0) and c(1) are calibration constants
c(8) and c(9) are used in equation pi (particle identification) which uses the total energy of S2 delE and E
c(8)=b is based on the semi-empirical range formula of charged particles in the detector r=(a*E)^b. You can find the values in the Knoll book of F. Goulding paper.
c(9) stretches the spectrum and can be a value between 0.5 and 5
w(0) and w(1) are windows for the (front strip - back strip) S2 delE values which should be a Gaussian peak centered around 2048, an arbitrary value so that the center is not around 0.
w(2) and w(3) are the same but for S2 E
w(17) and w(18) are windows used to select a particle type from the equation pi
w(19) and w(20) ignore. was used for a second delE and E detector in a previous experiment
6) summary of spectra created and used in sort3.f:
0-351 1D ADC histograms
3100-3115 1D scalar histograms
3000 1D channels w/ hits
3003 2D ADC ch vs ADC histo ch
3006 1D multiplicity + (3000 + number of good delE hits) + (3100 + number of good E hits)
3050 2D S2 delE front vs back
3060 1D S2 delE (front - back)
3051 2D S2 delE good front vs back
3061 1D S2 delE good (front - back)
3052 2D S2 E front vs back
3062 1D S2 E (front - back)
3053 2D S2 E good front vs back
3063 1D S2 E good (front - back)
4000 2D S2 energy vs annulus strip + (48 + energy vs annulus strip)
4010 2D S2 delE vs E
4011 1D pi
4012 2D pi vs tot energy
4013 1D total energy of selected particle using gates w(17) and w(18)
7) unpacking raw data in sort.f
ADC channels are put into channels 0-999
ADCs have 32 channels each and the first module is in the right slot of the VME crate.
TDC channels are put into channels 1000-1999 and ADC channel 1 corresponds to TDC channel 1001
TDC have 128 channels but only 112 channels are used. The first 16 are empty. So channels 16-127 from TDC#1 are put into channel 1000-1111 and from TDC#2 into channel 1112-1223, etc |