Initialization of Reconstruction Parameters

Subroutine CALOR_RECONST_INI
called once a run; makes an initialization of reconstruction parameters and counters. All reconstruction parameters are in COMMON / RECONPARAM /.

In the other environment as BRAHMS this subroutine can be called at any reasonable place.

Brief parameters description:

$\bullet$ energy_coeff1, energy_coeff2, energy_coeff - calibration coefficients are taken from muon run (10 GeV) to convert the hit amplitudes to physical scale [GEV] for two ECAL and one HCAL sampling structures;

$\bullet$ sph_cutE_sup, sph_cutE_clu, sph_cut_sup, sph_cut_clu - cuts for ``angular'' distances, used during super-cluster and cluster searching in ECAL and HCAL (see utility function angle_dist($\Theta$,$\Phi$) that calculates an angular distance);

$\bullet$ coeff_radE_sup, coeff_radE_clu, coeff_rad_sup, coeff_rad_clu - additional parameters for cluster ``angular'' radius; it makes sph_cuts to be proportional to Log(Energy);

$\bullet$ Cl_Dist_Merge - if the distance between centers of clusters less then this parameter - try to MERGE them;

$\bullet$ Cl_Dist_Merge_EM - allow to merge if distance between ECAL cluster and HCAL cluster less then this parameter;

$\bullet$ Dist_Merge - minimal distance from line or helix to cluster to MERGE ECAL and HCAL;

$\bullet$ EC_min_mu - minimal energy for ECAL clusters to be merged

$\bullet$ EC_max_mu - maximal energy for muon can deposit to be recognized as ECAL cluster;

$\bullet$ EM_Energy - there are special condition for EM clusters - if ECAL Gamma quanta or electron has an energy greater then this it can produce signal in HACL;

$\bullet$ x0_1ECAL, x0_2ECAL - number of X0 in one sampling in first and second ECAL structures, used in subroutine EClust_Shower_Hypo;

$\bullet$ R_moliere - Moliere radius of ECAL first structure

$\bullet$ theor_mu_densE, theor_mu_dens - the average theoretical density in [MeV/hit] for the muon tracks; used in subroutines Eclust_axsis, Eclust_property and Hclust_axsis;

$\bullet$ ec1_theor_dens, ec2_theor_dens, ec1_theor_RMS, ec2_theor_RMS - coefficients to calculate predicted, energy dependent theoretical density in [MeV/hit] and its error for the electromagnetic and hadronic showers. Special electron and gamma runs were made. Afterwards the energy density and its RMS energy dependence were fitted with the simplest $theor\_dens = C1 + C2*\sqrt{E_{clust}}$, $theor\_dens\_RMS = C3 + C4*\sqrt{E_{clust}}$ formulae. Coefficients are used in subroutines Eclust_axsis, Eclust_property to appropriate cluster particle identification hypothesis in ECAL;

$\bullet$ nindbeg_thr, nindend_thr - thresholds used for particle identification in ECAL;

$\bullet$ thresh_clust - threshold distance for ECAL cluster joining; used in subroutine EC_CLUST_JOIN

$\bullet$ thresh_hits - threshold distance for ECAL hit joining; used in subroutine EC_HIT_JOIN

$\bullet$ nhit_thrE, Esup_clustE_thr - thresholds to create cluster from super-cluster ECAL

$\bullet$ nhit_thr, Esup_clust_thr - thresholds to create cluster from super-cluster HCAL

$\bullet$ TPC_trk_Dist_Merge - the largest allowed distance between tracks at the calorimeter face to be a group

$\bullet$ TPC_Cl_Dist_Merge, TPC_Cl_Dist_MergeH - if the distance between center of cluster and track is larger then this parameter it will not be calculated at all (save CPU time)

$\bullet$ TPC_Dist_Merge, TPC_Dist_MergeH - maximal distance from TPC particle to cluster to be MERGED (from strait line along TPC track direction to the center of ECAL/HCAL Cluster)

$\bullet$ Energy_sigmaE - actually this value is near by sigma for ECAL 12%/sqrt(E).

$\bullet$ Energy_sigmaH - actually this value is near by sigma for HCAL 50%/sqrt(E).
Parameters to choose the subtraction or substitution for the charged tracks with TPC momentum.

$\bullet$ -