Tuda logbooks General E1104 ISACII E1195 E1103 E1147 S1287 S1203 S1233 S1284 S1447 23Na(a,p) S1540 S1381  S1847 Catania
  S1284, Page 1 of 1  ELOG logo
IDdown Date Author Type Category Subject
  3   Fri Nov 29 18:20:08 2013 Naomi, Titti, Angelo  delE vs E spectrum

 for run R54,55,56 spectrum number S4010 in sort3.f

  2   Fri Nov 29 16:00:09 2013 Naomi GeneralMIDASUK sort

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

  1   Tue Nov 26 14:16:54 2013 Naomi, Manuela, Diana, Mariangela Generalalpha particle run screen shots
ELOG V2.9.2-2455