Difference between revisions of "Gams reconstruction"

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)