hec_his_ana1.f        [#]

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

    ( is called by:   hec_ped_pass1.f  , hec_online_ana1.f  )
 

    includes the following commons:

    hec_par.inc
    hec_datacard.inc
    hec_evth.inc
    hec_geo.inc
   hec_his_sys.inc
    hec_ped.inc
    hec_ped_sys.inc
    hec_runh.inc
    hec_tdc.inc
    hec_hvpar.inc
    hec_slow_pc.inc
    hec_imp.inc
 
 

    This is quite a lengthy routine.
    First it vzero-s some arrays.

    --- Then start looping over all adc channels here

    If histogram switch:    histo_switch(1)=1     then fill pedestal histograms and
                                                                         histograms of the adc-rms of the run.

    If histogram switch:    histo_switch(1)=1     then use the moments (see: hec_his_evt_hfill_ch.f)
                                                                        and fill the arrays for the average maximum
                                                                        and the error on it.
                                                                        Also fill a 2-dim histogram (phi,eta) with
                                                                        weighted hits in each layer

    If histogram switch:    histo_switch(3)=1     then loop over all used time slices and calculate:
                                                                        average = adcmom(ic,i_t,2) / adcmom(ic,i_t,1)
                                                                        use this average as weight in filling time-slice
                                                                        histogram.
 

    If histogram switch:    histo_switch(4)=1     then pack the average tdc and its error into
                                                                         a histogram. (use moments again)

    --  Here ends the loop over all channels
 
 

    If histogram switch:    histo_switch(1)=1     then pack adc-rms for run versus channel number
                                                                        into histograms. Also pedestals vs. channel.
                                                                        This is for the 4 channels chosen in:
                                                                       hec_his_evt_fill.f

    If histogram switch:    histo_switch(2)=1     Pack the average adc maximum vs. channel number.
                                                                        Fill signal sums as a function of z coordinate (depth)
                                                                        into histogram.
 

    If histogram switch:    histo_switch(7)=1     for run period 7, fill slow control histograms.
 
 

    Loop over   clusters. If particle id matches chosen cluster title: copy histograms around.