source: JEM-EUSO/equalization_gain/trunk/README @ 197

*************************************************************
*
* README file for gain equalisation
* SDC le 25/09/2013
*
**********************************************************

Repository
===========

to see the projects for JEM-EUSO just dou :

==> svn ls https://svn.lal.in2p3.fr/projects/JEM-EUSO

==> svn ls https://svn.lal.in2p3.fr/projects/JEM-EUSO/equalization_gain
branches/
tags/
trunk/


The Gain Equalisation files are inside the analysis directory.

The analysis directory contains :
- this Top level README
- the MySetup.sh file


1) configure the MySetUp file to define the proper environnement variables
   - srcDIR : repository which contains the analysis code and the scripts files 
   - dataDIR : contains the data

2) the analysis directory further contains 3 subdirectories
   - the src directory that contains the sources
   - the data directory that contains the data to analyses
   - the work directory (or any other name) that should be created for the analysis



SRC directory
===============

Makefile : to check the coding is correct from C++ compiler point of view

PlotScurvesAllMain.cc : The main for the compilation used from now on for the analysis
PlotScurvesAllMain.h

PlotScurvesAll-diffFit.cc : The only file that does really the gain equalisation inside the PlotScurve function
PlotScurvesAll-diffFitConst.h : define the constant used by PlotScurve function in the file PlotScurvesAll-diffFit.cc
PlotScurvesAll-diffFit.h


Please do

make 
-----
to check that PlotScurvesAll-diffFit.cc syntax is correct for the C++ compiler

makeDataFile.sh : Extract the data from the EC-ASIC board read by the Test-EC-ASIC into one file of data per ASIC channel
                  Create directories, one per ASIC.
                  Link *.h and *.cc files locally for analysis 
makeGainTable.sh : extract the table for the ASIC obtained by the output result of the execution of PlotScurve() inside PlotScurvesAll-diffFit.cc 

PlotScurvesAll-diffFit.C ! Old file uncompiled by Hiroko (23/09/2013)




DATA directory
================
/Users/dagoret-campagnesylvie/MacOsX/EUSO-BALLOON/GainEqualisation/analysis/data>ls -l 
total 16
20130923-APC contains S-Curves
20130924-LAL contains S-Curves
gain-default-ASIC_C.txt : Default ASIC Gain
gain-default-ASIC_D.txt

Procedures
==========

From Hiroko : Monday 23 09 2013

First of all, things are not fully automatick yet as we are still
looking for the best way to carry it out.


----------------------------------------------------------------------------
<STEP 1> Organize your directories (I prefer to use bash.sh instead of csh,tcsh):
% mkdir YOUR_DIR/analysis
% mkdir YOUR_DIR/analysis/src
% mkdir YOUR_DIR/analysis/data
% export YOUR_DIR/analysis/src srcDIR
% export YOUR_DIR/analysis/data dataDIR

make your data directory for the day. e.g.,:
% mkdir $dataDIR/yYYmMMdDD-APC
% cd yYYmMMdDD-APC
----------------------------------------------------------------------------

----------------------------------------------------------------------------
<STEP 2> Modify data file obtained by LabView to fit and make
data directories for the analysis:
There are many ways to do this, and you can do as you wish.
Ask Pierre to make the data file which includes DAC values
and data for ASIC_C and ASIC_D at the beginning as he knows
what I mean.

If you want, the data obtained by LabView can be used as itself
with modifying the parameter on my makeDataFile.sh script.
I leave it for right now because we haven't fixed the data format
yet, at least by the last time I was there.


% cp -p data_s_curves_runXX.data .
% sh $srcDIR/makeDataFile.sh data_s_curves_runXX

This sh automatically makes two directories data_s_curves_runXX_C
and data_s_curves_runXX_D for the data of ASIC_C and ASIC_D respectively.

From next steps, it totally depends on which data you are taking.
Assume you want to analyse the data of ASIC_C.
----------------------------------------------------------------------------

----------------------------------------------------------------------------
<STEP 3> Run ROOT, search the suitable parameters for each data
% cd data_s_curves_runXX_C
% ln -s ../FORMER_GAIN_DATA.txt gain-org.txt

Usually, the following command is not necessary as long as the data
includes pedestal, otherwise:
% ln -s ../$pedestal_data_directory_D/*.dat .

Check the parameters in header file or directory input the parameter
in main function as you can see in the source code.
This totally depends on the data. You have to modify it to suitable
parameters.

% root -l

Set the parameters on the top of the program depending on your analysis.
• General policy: parameters starting with "k" at the beginning is
fixed gloval parameter and usually set in the header file, never change
through the program.
• kDrawCh: channel just you want to check on display. doens't effect
the analysis. 
• fname, fname_ped (kFname, kFnamePed): "YOUR DATA FILE" without ".data".
See examples in the header file.
• dac_half_spe_fix: set DAC value only when you want to adjust the gain
to particular pulse height (DAC value) instead of either average of the
pulse height of reference pixel. Usually just leave it -1.
• fit_min, fit_max (kFitMin, kFitMax): leave them as default,
it is not currently used for the analysis.
• fit_min_scurve, fit_max_scurve (kFitMinS, kFitMaxS): define the range
of s-curve fitting. search proper values by yourself with running root
program a few times.
• kBoard: don't forget to change 0(C) or 1(D) depending on which ASIC
you're analysing
• kChRef (kRefCh): a pixel whose gain remains 64 thorough the gain
adjustment.
• checkAve: for the first run, set this parameter 1, and find a pixel
which remains 64. In the case of parameter 1, all the gains are adjusted
to the average pulse height.
Once you choose the reference pixel, set the parameter 0, as well as
set kRefCh the pixel number. (better to keep in the header file for the
record of each PMTs.)

Once you set all the parameters above, then
[].L PlotScurvesAll-diffFit.C
[]PlotScurve()

or do in another way according Sylvie's 
from the shell (not inside root), execute directly the binary

$(srcDIR)/PlotScurvesAllMain 


Repeat until you define all the parameters above.
----------------------------------------------------------------------------


----------------------------------------------------------------------------
<STEP 4> Make the gain table
Once you decide the parameters and modify the header file, then run
% root -l > tmp.data
[].L PlotScurvesAll-diffFit.C
[]PlotScurve()
[].q

If you have the data of ASIC_D, go to the directory of ASIC_D
and do the same.

% cd ../
Check the paramters in $srcDIR/makeGainTable.sh and run
You should add/remove #s depending on which data you're analysing.
% sh makeGainTable.sh
Or type by your hands the commands in the script file.

If everything has been done properly, you should get a gain
table with a proper format.
I would recommend that when you make tmp.data file, to include
the value of threshold used, and the number of iteration, e.g.,
tmp-th180-it0.data
And then do the iteration of the analysis with the data with
the obtained gain table until you get reasonable uniformity.

NOTE that, sometimes PMT gain or efficiency is too low to obtain
the inflection point, in this case the program outputs the -1 for
flag_fit. If it remains -1 after some trials, you should adjust
by hand. Ask Philippe for detailes.
----------------------------------------------------------------------------


----------------------------------------------------------------------------
<STEP 5> Adjust gain of EC_unit
Once you succeed to obtain all the gain tables:
1) check the pulse height (dac_sfit[kChRef]-dac_ped[kChRef])
2) choose 1 of 4 PMTs in EC_unit
3) take the ratio of pulse height between reference pixels
4) multiply the gain table by the ratio above
5) combine all the gain tables
----------------------------------------------------------------------------



Good luck!


Logbook
=======


up to run95, the definition of iteration:

T1, T2: board1, board2
no description: data with ASIC gain 64 everywhere,
it0: data after the 1st adjustment,
it1: data after the 2nd adjustment, ...
it2: data after the 3rd adjustment, ...

should be changed afterward:
T1, T3: board1, board2
it0: data with ASIC gain 64 everywhere, 
it1: data after the 1st adjustment,
it2: data after the 2nd adjustment, ...





run60: 2 ch (44,56), 0.562
run62: 2ch it0 (by hand)
run63: 64ch it0 (by file, sfit)
run64: 64ch it1 strange noise default gain
run65: 384 ch (128ch HVon default gain)
run66: 64 ch ASIC_D default gain
run67: sfit it0 th180 100 events
run68: sfit it0 th180 modified 10000 events ch62:103->92, ch60:109->95
run69: sfit it1 th240 100 events
run70: sfit it1 th240 10000 events
run71: 64 ch ASIC_D default gain


run72: ASIC_D gain 64 others 0 100 events, step 10, 
run73: with adjust gain for ASIC_D th240 it1 100 events, step 10
74: light off, HV ON ped step 2 ASC_C
75: light off, HV ON ped step 1 ASC_C -> same as 76
? forgot to save
76: 10000 events                   

78: T2 (board2) ASIC_C default 10000 events
79: T2 (board2) ASIC_D default 10000 events

80: T2 ASIC_C it0 100 events
81: T2 ASIC_C it0 10000 events
82: T2 ASIC_D it0 10000 events
83: T1 ASIC_C it0 10000 events
84: T2 ASIC_C it1 100 events
85: T2 ASIC_C it1 10000 events
86: T2 ASIC_D it1 10000 events
87: T2 ASIC_C it1 10000 events

89: T1 ASIC_C it1 100 events
90: T1 ASIC_C it1 10000 events
91: T1 ASIC_C it2 10000 events
92: T1 ASIC_C it2 10000 events, lower luminosity (50%)
93: T1 ASIC_C it3 (with using 91) 10000 events, lower luminosity (50%) 

94: T1 ASIC_C it3 (with using 91) 10000 events, lower luminosity (50%),
modified the gain of ch39,ch43,ch47 by hand
aiming threshold at 150
95: T1 ASIC_C it3 (with using 91) 10000 events, lower luminosity (50%),
modified the gain of ch39,ch43,ch47 by hand
with HV of CW corresponding to 900V


tmp tables for each PMTs:
Board  ASIC    dir incl. table_fin, chRef, data\w table_fin, ped, peak, PH, ratio
T1     ASIC_C: run91, chRef 31, run 94&95, 72.5 232.5, 160, 1.27
T1     ASIC_D: run68, chRef 31, run70,     74.5 277.5, 203, 1
T2(T3) ASIC_C: run87, chRef 37, run85&87,  72.5 287.5, 215, 0.94
T2(T3) ASIC_D: run86, chRef 32, run86,     67.5 252.5, 180, 1.13


History
=======

------------------------------------------------------------------------------
Up to 23/09/2013 : code developped by Hiroko
------------------------------------------------------------------------------
24/09/2013 : Sylvie introduces PlotScurvesAll-diffFit.cc and makes it compilable with the makefile. 
------------------------------------------------------------------------------
24/09/2013 : From Pierre :

So you have a pptx file which is the summary of the data taken yesterday.
The you have two types of data file:
- s_curves_runXXX.data (where XXX is the run number) : this is the raw data with the first column being the DAC value et then you have 384 columns corresponding to the 64 channels of ASIC A, B, C, D, E and F
- data_s_curves_runXXX.data (where XXX is the run number) : this is the data from the two used ASICs (C and D). 1st column is still the the DAC values, then 128 columns corresponding to the 64 channels of ASICs C and D

I don't remember if Hiroko uses first or second type of file for her analysis.
cheers
pierre 
------------------------------------------------------------------------------
25/09/2013 : Improve a little bit PlotScurvesAll-diffFit.cc with welcome comments inside the code.
             add the PlotScurvesAll-diffFitConst.h
             Put Canvas name and title to understand better what is seen on screen
------------------------------------------------------------------------------