Gams reconstruction
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)