Difference between revisions of "Gams reconstruction"
Line 41: | Line 41: | ||
'''adcgam_bk.inc''' Reconstruction common looks like: | '''adcgam_bk.inc''' Reconstruction common looks like: | ||
+ | Real: | ||
energy_adcgam(1,i) - energy in GeV | energy_adcgam(1,i) - energy in GeV | ||
x_adcgam(1,i) - log. method x-coord (angular correction applied) | x_adcgam(1,i) - log. method x-coord (angular correction applied) | ||
Line 49: | Line 50: | ||
chi2_adcgam(1,i) - chi2 | chi2_adcgam(1,i) - chi2 | ||
+ | Integer: | ||
type_adcgam(1,i) - type of the gamma: | type_adcgam(1,i) - type of the gamma: | ||
= 0 default value | = 0 default value | ||
Line 54: | Line 56: | ||
+= 10 for clusters selected from mutlimaxima island | += 10 for clusters selected from mutlimaxima island | ||
+= 100 for LG clusters in hybrid detector | += 100 for LG clusters in hybrid detector | ||
− | = 50 + type1 + type2 for merged together subclusters 1 and 2 (transition case) | + | = 50 + type1 + type2 for merged together subclusters 1 and 2 (transition case)<br> |
− | + | (hole size is hardcoded to 5x5 cells for PWO and 3x3 cells for LG) | |
− | hole size is hardcoded to 5x5 cells for PWO and 3x3 cells for LG | + | |
Revision as of 12:33, 5 September 2013
Short description of the reconstruction package developed for GAMs
(with changes and updates)
Package includes:
adcgam_bk.inc common with output reconstruction
phot_cell_cbk_mc.inc common with shower profile parameters
cphoto.inc common with control keys, flags and other aux. variables
addhits_lib.F generates hits based on shower library data
addhits_sim.F performs hits generation using fast simulation based on profile functions
isim.F main routine calling others to generate hits
recon_sim.F contains reconstruction package itself
guout.F call recon_sim() from here to get event reconstruction
gustep.F call isim(chdet,chsimopt) to simulate and store hits. Normally, after that particle stops in GUSTEP (ISTOP = 1) to avoid multiple counting
Files description:
adcgam_bk.inc Reconstruction common looks like:
Real:
energy_adcgam(1,i) - energy in GeV x_adcgam(1,i) - log. method x-coord (angular correction applied) y_adcgam(1,i) - log. method y-coord (angular correction applied) z_adcgam(1,i) - z-coord xc_adcgam(1,i) - center of gravity x-coord (no angular correction) yc_adcgam(1,i) - center of gravity y-coord (no angular correction) chi2_adcgam(1,i) - chi2
Integer:
type_adcgam(1,i) - type of the gamma: = 0 default value += 1 or 2 for clusters selected by splitting single maxima shape += 10 for clusters selected from mutlimaxima island += 100 for LG clusters in hybrid detector = 50 + type1 + type2 for merged together subclusters 1 and 2 (transition case)
(hole size is hardcoded to 5x5 cells for PWO and 3x3 cells for LG)
dime_adcgam(1,i) - number of counters in the gamma id_adcgam(1,i) - id of the cell with the max energy deposition: = 10*col+row
nadcgam - number of reconstructions in event (i = 1 ... nadcgam)