model = SM_CKM

include("sindarins/processes.sin") 
include("sindarins/cuts.sin") 
include("sindarins/process_input.sin") 
!  The following works: if preint is declared and set to true on the command-line
!  all is fine, because the second include over-rules the firs. But also if it is not 
!  declared at all: the if statement undefined, which means *no* branch will be
!  executed: all is still fine!
include("sindarins/integration_input.sin") 
if ( ?preint ) then
  include("sindarins/pre_integration_input.sin") 
endif  
include("sindarins/simulation_input.sin") 
include("sindarins/parameter_input.sin") 
include("sindarins/beams.sin")






# ====================== In 1.95: actually doing it =============== 
! ..... In 1.95, this could be controled on the
! command line with make prg; whizard --(whizard.in section)=...

!>> 1.95: make prg:
compile

!! example: set lumi, Ecms and polarisation (can be done before reading the file, but
!! note that they must NOT be set in the file as well, because the file (oddly) takes
!! precedence over the command line.
!! Addend (on the fly) what you want to do to the end of the file, for ~ the same
!! reason (it makes no sense before the file has been read, because the processes
!! are not yet defined 

!! whizard -e  'luminosity = 500  sqrts = 500 GeV beams_pol_density = @(1),@(-1) $sample = "ww_sl0muq"'\
!!   <( cat dbd-2f-transl-nophincut.sin ; 
!!      echo 'simulate ( ww_sl0muq ){ ?polarized_events = true }' )


! if not on the commad line:

!!>> !>> 1.95: whizard --process_input='luminosity=0 "process_id="ww_sl0muq"' :
!!>> 
!!>> luminosity = 500  
!!>> sqrts = 500 GeV 
!!>> beams_pol_density = @(-1),@(1)
!!>> 
!!>> integrate (ww_sl0muq){ ?polarized_events = true ?vis_channels = true}
!!>> 
!!>> 
!!>> 
!!>> !>> 1.95: whizard --process_input='process_id="ww_sl0muq"' --simulation_input='write_events_file="ww_sl0muq"':
!!>> 
!!>> ?rebuild_events = true
!!>> $sample = "ww_sl0muq"
!!>> simulate ( ww_sl0muq ){ ?polarized_events = true }
!!>>