!-----------------------------------------------
!        Electron-positron annihilation
!-----------------------------------------------

!-----------------------------
!        Define run
!-----------------------------

application ee
set engy 91.2
set nevent 10000
set modsho 100                                                
set ndecay 1000000
set irescl 0 !otherwise crash in utresc

!------------------------------------------------
!            Printout particle list 
!------------------------------------------------

!print * 2   !uncomment to get list in z-expl2.check 

!----------------------------------------------------
!     Define analysis (multiplicity distribution)
!----------------------------------------------------

beginanalysis 
  histogram 
    mulevt !x = multiplicity
    numevt !y = number of events 
    1      !normalisation (per event)
    -0.5   !x-min
    61.5   !x-max
    31     !number of bins
  idcode  9970 !charged particles
endanalysis

!---------------------------------------------------
!      Write out final results to output file
!---------------------------------------------------

write "multiplicity distribution"
histoweight 
writearray 3 

!----------------------------------------------------
!     Define analysis (rapidity distribution)
!----------------------------------------------------

beginanalysis 
  histogram absrap numptl 11 0 6 30 !absrap = absolute value of rapidity
  frame thrust !particular frame used in e+e-
  idcode 9970  
endanalysis

!---------------------------------------------------
!      Write out final results to output file
!---------------------------------------------------

write "rapidity distribution"
histoweight 
writearray 3

!----------------------------------------------------
!     Define analysis (xp distribution)
!----------------------------------------------------

beginanalysis 
  frame total
  binning log
  histogram xp numptl 11 0.001 1 30 
  idcode 9970  
endanalysis

!---------------------------------------------------
!      Write out final results to output file
!---------------------------------------------------

write "xp distribution"
histoweight 
writearray 3

!----------------------------------------------------
!     Define analysis (xi distribution)
!----------------------------------------------------

beginanalysis 
  frame total
  binning lin
  histogram xi numptl 11 0.1 6 30 
  idcode 9970  
endanalysis

!---------------------------------------------------
!      Write out final results to output file
!---------------------------------------------------

write "xi distribution"
histoweight 
writearray 3  

!----------------------------------------------------
!     Define analysis (ptr-distr)
!----------------------------------------------------

beginanalysis 
  frame thrust
  histogram p1a numptl 11   0.0001  10.0   20
  trigger rap   -10.0    10.0
  idcode       9970
endanalysis

!---------------------------------------------------
!      Write out final results to output file
!---------------------------------------------------

write "ptr distribution"
histoweight  
writearray 3 

echo off

!##########################################
!  the normalization is a 
!    number hijk with the following options
!##########################################
!
!  k  0:  * 1
!     1:  / number of events
!     2:  / number of triggered events
!     4:  / bin-counts
!     5:  / bin sum
!     6:  / number of summed bin-counts (yield=1.)
!     7:  uses same normalization as one histo before
!
!  j  0:  * 1
!     1:  / bin-width
!     2:  * sigma_total / bin-width
!     3:  * sigma_diff / bin-width
!
!  i  0:  * 1
!     1:  y => y*x
!     2:  y => y/x/2/pi (modified for mt0)
!     3:  kno-scaling
!     4:  y => y/x**1.5
!     5:  y => y/x
!     6:  y => y*xi (for conex, xi=x of the bin)
!     7:  y => y/x/(x-m)
!
!  h  0: normal
!     1: accumulated