Hi all-
It looks like this ADC sparsification macro is nearly done. I say
nearly as it still may require fine tuning but that is not apparent yet.
I have most recently incorporated the LADS into this macro. To review, it
works like this:
-the code calculates pedestals by requiring tdc <0 or tdc > 4095
-it then fits the pedestal with a gaussian and determines sigma
-the threshold is then determined as the pedestal mean plus
sigma*MULTIPLIER where the multiplier is some number which is easily
adjustable in the first few lines of the code for each type of detector.
i.e., from the code,
///////////////////////////////////////////////////////
//MULTIPLIERS TO CALCULATE THRESHOLDS
//i.e. THRESHHOLD = PEDESTAL + MULTPLIER*PEDESTAL_SIGMA
//BELOW ARE THE MULTIPLIERS, ADJUST AS REQ.
int TOF_MULT = 3;
int NC_MULT = 3;
int CC_MULT = 3;
int L20_MULT = 3;
int L15_MULT = 3;
///////////////////////////////////////////////////////
-empty channels get a default threshold of 1000
The code scrolls a list of detectors and thresholds including
spare channels so that it can be compared to the current electronics map.
Four files are written to the current directory with the run number(s) as
part of the filename. E.g. for run 4100, the code will produce:
LADCsparsify:4100
RADCsparsify:4100
L_LADS_ADCsparsify:4100
R_LADS_ADCsparsify:4100
where the top two files are for the SC, CC, and NC ADCs
and the bottom two files are for the LADS ADCs
these files can then be renamed and copied into the appropriate place so
that they can be used in taking ADC data and hopefully reduce our disk
space usage!
I put the 1100 line macro (it's a long macro! s/b compiled code?) into:
/home/daq/blast/pro2003/analysis/macros/Raw/ADCsparsify.C
Pete
PS Eventually the BATS should become part of this code.
----------------------------------------------
Pete Karpius
Graduate Research Assistant
Nuclear Physics Group
University of New Hampshire
phone: (603)862-1220
FAX: (603)862-2998
email: karpiusp@einstein.unh.edu
http://pubpages.unh.edu/~pkarpius/homepage.htm
----------------------------------------------
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:30 EST