anchor.diff:  #endif
     end if
   
  -!++  Initialize efficiency array
  +  ! Initialize efficiency array
     do i=1,iu
       sum_ax(i)   = zero
       sqsum_ax(i) = zero
  
positive.code:   write(lout,"(a,l1)")    'COLL> Info: DO_COLL             = ', do_coll
    write(lout,"(a,l1)")    'COLL> Info: DO_NSIG             = ', cdb_doNSig
    do i=1,cdb_nFam
      write(lout,"(a,a19,a3,f13.6)") "COLL> Info: ",chr_rPad("NSIG_"//trim(cdb_famName(i)),19)," = ",cdb_famNSig(i)
    end do
  
    write(lout,"(a)")
    write(lout,"(a)")       'COLL> INPUT PARAMETERS FOR THE SLICING:'
    write(lout,"(a)")
    write(lout,"(a,i0)")    'COLL> Info: N_SLICES            = ',n_slices
    write(lout,"(a,e15.8)") 'COLL> Info: SMIN_SLICES         = ',smin_slices
    write(lout,"(a,e15.8)") 'COLL> Info: SMAX_SLICES         = ',smax_slices
    write(lout,"(a,e15.8)") 'COLL> Info: RECENTER1           = ',recenter1
    write(lout,"(a,e15.8)") 'COLL> Info: RECENTER2           = ',recenter2
    write(lout,"(a)")
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(1,1)        = ',jaw_fit(1,1)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(1,2)        = ',jaw_fit(1,2)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(1,3)        = ',jaw_fit(1,3)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(1,4)        = ',jaw_fit(1,4)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(1,5)        = ',jaw_fit(1,5)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(1,6)        = ',jaw_fit(1,6)
    write(lout,"(a,e15.8)") 'COLL> Info: SCALING1            = ',jaw_ssf(1)
    write(lout,"(a)")
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(2,1)        = ',jaw_fit(2,1)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(2,2)        = ',jaw_fit(2,2)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(2,3)        = ',jaw_fit(2,3)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(2,4)        = ',jaw_fit(2,4)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(2,5)        = ',jaw_fit(2,5)
    write(lout,"(a,e15.8)") 'COLL> Info: jaw_fit(2,6)        = ',jaw_fit(2,6)
    write(lout,"(a,e15.8)") 'COLL> Info: SCALING2            = ',jaw_ssf(2)
    write(lout,"(a)")
  
  !SEPT2005
  !
  ! HERE WE CHECK IF THE NEW INPUT IS READ CORRECTLY
  !
    write(lout,"(a,e15.8)") 'COLL> Info: EMITXN0_DIST        = ', emitnx0_dist
    write(lout,"(a,e15.8)") 'COLL> Info: EMITYN0_DIST        = ', emitny0_dist
    write(lout,"(a,e15.8)") 'COLL> Info: EMITXN0_COLLGAP     = ', emitnx0_collgap
  
negative.code:   write(lout,"(a,e15.8)") 'COLL> Info: SIGSECUT2           = ', sigsecut2
    write(lout,"(a,e15.8)") 'COLL> Info: SIGSECUT3           = ', sigsecut3
    write(lout,"(a)")
    write(lout,"(a,i0)")    'COLL> Info: NAPX                = ', napx
    write(lout,"(a,e15.8)") 'COLL> Info: Sigma_x0            = ', sqrt(mybetax*myemitx0_dist)
    write(lout,"(a,e15.8)") 'COLL> Info: Sigma_y0            = ', sqrt(mybetay*myemity0_dist)
    write(lout,"(a)")
  
    ! Initialize random number generator
    if(rnd_seed == 0) rnd_seed = mclock_liar()
    if(rnd_seed <  0) rnd_seed = abs(rnd_seed)
    rnd_lux = 3
    rnd_k1  = 0
    rnd_k2  = 0
    call rluxgo(rnd_lux, rnd_seed, rnd_k1, rnd_k2)
    write(outlun,*) 'INFO>  rnd_seed: ', rnd_seed
  
    ! Call distribution routines only if collimation block is in fort.3
    cdist_logUnit = outlun
    cdist_energy  = myenom
    cdist_alphaX  = myalphax
  