Ticket #61: Run.dat

(run){
  !-------------------------------------------------------
  !--- Run parameters ------------------------------------
  !-------------------------------------------------------
  EVENTS          = 5000      ! Number of Events
  NUM_ACCURACY   = 1.e-10       !
  OUTPUT          = 18       ! Output level:
  !                             ! 0=errors,1=events,2=run info,4=Tracking message
  !                             ! 3=events+information, 8=Debugging message
  !
  ANALYSIS        = 0                   ! Analysis: 0=off, 1=on
  !  ANALYSIS_OUTPUT = Analysis.MI
  !
  LOG_FILE      = sherpa_Ttop800_1.log        ! log file
  RANDOM_SEED = 77345 76772
  FILE_SIZE = 5000
  MASSIVE[5] 1
  !EVENT_MODE = Sherpa
  !EVENT_MODE = Hepevt
  EVENT_MODE = HepMC
  !HEPEVT_OUTPUT = Results1
  HEPMC2_GENEVENT_OUTPUT = Results1_HepMC_GENEVENT
  !HEPMC2_ASCII_OUTPUT = Results1_HepMC_ASCII
  !HEPMC2_OUTPUT = Results_HepMC
  !SHERPA_OUTPUT = Results1
  !ANALYSIS_OUTPUT = Results1
}(run)

(beam){
  !-------------------------------------------------------
  !-- Beam parameters ------------------------------------
  !-------------------------------------------------------
  BEAM_1          = 2212                 ! possible beam particles:  P+, P-, e+, e-
  BEAM_ENERGY_1   = 7000.                ! in GeV
  BEAM_POL_1      = 0.                   ! Polarization degree -1 ... 1
  BEAM_SPECTRUM_1 = Monochromatic        ! Monochromatic,Laser_Backscattering,Simple_Compton
  K_PERP_MEAN_1 = 0.2
  K_PERP_SIGMA_1        = 0.8
  !
  BEAM_2          = 2212                 ! possible beam particles:  P+, P-, e+, e- 
  BEAM_ENERGY_2   = 7000.                ! in GeV
  BEAM_POL_2      = 0.                   ! Polarization degree -1 ... 1
  BEAM_SPECTRUM_2 = Monochromatic        ! Monochromatic,Laser_Backscattering,Simple_Compton
  K_PERP_MEAN_2 = 0.2
  K_PERP_SIGMA_2        = 0.8
  !
  !BEAM_SMIN       = 1.e-10               ! Minimal fraction of nominal s after beam spectra
  !BEAM_SMAX       = 1.0                  ! Maximal fraction of nominal s after beam spectra
  !
  !E_LASER_1          = 1.17e-9           ! Laser energy in GeV
  !P_LASER_1          = 0.                ! Laser polarization +-1
  !E_LASER_2          = 1.17e-9           ! Laser energy in GeV
  !P_LASER_2          = 0.                ! Laser polarization +-1
  !LASER_MODE         = 0                 ! 0 = all, 1,2,3 = individual components
  !LASER_ANGLES       = Off               ! On/Off
  !LASER_NONLINEARITY = On                ! On/Off
}(beam)

(isr){
  !-------------------------------------------------------
  !-- ISR parameters ------------------------------------
  !-------------------------------------------------------
  BUNCH_1         = 2212      ! possible beam particles:  P+, P-, e+, e-
  ISR_1           = On        ! On/Off
  ! 
  BUNCH_2         = 2212      ! possible beam particles:  P+, P-, e+, e- 
  ISR_2           = On        ! On/Off
  !
  ISR_SMIN        = 1.e-10    ! Minimal fraction of nominal s for parton after ISR
  ISR_SMAX        = 1.0       ! Maximal fraction of nominal s for parton after ISR
  !
  !ISR_E_ORDER     = 1         ! Perturbative order of electron structure function
  !ISR_E_SCHEME    = 2         ! Beta-scheme : 0,1,2 , default = 2 
  !
  PDF_SET         = cteq6m.LHpdf    ! the CTEQ6 PDFs cteq6m,cteq6d,cteq6l,cteq6l1
                            ! can be accessed when PDF_GRID_PATH = CTEQ6Grid
                            !
                            ! to use LHAPDF list one of the PDF sets,
                            ! , e.g. Alekhin_100.LHpdf, that comes with your 
                            ! LHAPDF installation, the full list can be found in
                              ! 'lhapdf-config --pdfsets-path', in addition set
                            ! PDF_GRID_PATH = PDFsets
                            !
                              ! or use the interface to MRST99 through
                            ! PDF_GRID_PATH = MRST99Grid
                              ! 
  PDF_SET_VERSION = 1
  !
  PDF_GRID_PATH   = PDFsets ! CTEQ6Grid,PDFsets,MRST99Grid
}(isr)

(processes){
  Process 93 93 -> 6[a] -6[b] 93{1};
  DecayOS 6[a] -> 5 24[c];
  DecayOS 24[c] -> 94 94;
  DecayOS -6[b] -> -5 -24[d];
  DecayOS -24[d] -> 94 94;
  CKKW sqr(30/E_CMS);
  Order_EW 4;   
  Selector_File Run.dat|(sel0j){|}(sel0j) {6}
  Selector_File Run.dat|(sel1j){|}(sel1j) {7}
  Integration_Error 0.02 {3};
  End process;

}(processes)


(sel0j){
!selector for the 6-particle final state
 ! cut on pT of 1st combination: top
 "Calc(PPperp(p[0]+p[1]+p[2]))" 5,93,93 800,10000;
 ! cut on pT of 6th combination: anti-top
 "Calc(PPperp(p[0]+p[1]+p[2]))" -5,93,93 0,1e12:0,1e12:0,1e12:0,1e12:0,1e12:800,1e12;
}(sel0j)


(sel1j){
!selector for the 7-particle final state
 ! cut on pT of 1st combination: top
 "Calc(PPperp(p[0]+p[1]+p[2]))" 5,93,93 800,10000;
 ! cut on pT of 8th combination: anti-top
 "Calc(PPperp(p[0]+p[1]+p[2]))" -5,93,93 0,1e12:0,1e12:0,1e12:0,1e12:0,1e12:0,1e12:0,1e12:800,1e12;
}(sel1j)

(selector){

!  PT 6 850 14000
!  PT -6 850 14000
!  PseudoRapidity 6 -2.5 2.5
!  PseudoRapidity -6 -2.5 2.5

}(selector)

(model){
  !-------------------------------------------------------
  !-- Model parameters -----------------------------------
  !-------------------------------------------------------
  MODEL                 = SM       ! Model 
  !
  ! SM parameters
  !
  EW_SCHEME             = 0        ! which parameters define the ew sector.
  ALPHAS(MZ)            = 0.118    ! strong coupling at scale M_Z 
  ALPHAS(default)       = 0.0800   ! strong coupling
  ORDER_ALPHAS          = 1        ! NLO
  1/ALPHAQED(0)         = 137.036  ! inverse of alpha QED in the Thomson limit
  1/ALPHAQED(default)   = 132.51   ! inverse of alpha QED
  SIN2THETAW            = 0.2222   ! Weinberg angle at scale M_Z
  VEV                   = 246.     ! Higgs vev        
  LAMBDA                = 0.47591  ! SM Higgs self coupling
  CKMORDER              = 0        ! order of expansion of CKM matrix in Cabibbo angle
  !CABIBBO               = 0.2272   ! Cabibbo angle (Wolfenstein parametrization)
  !A                     = 0.85     ! A (Wolfenstein parametrization)
  !RHO                   = 0.50     ! rho (Wolfenstein parametrization)
  !ETA                   = 0.50     ! eta (Wolfenstein parametrization)
  !
  ! SUSY parameters
  !
  !SLHA_INPUT            = LesHouches_SPS1A.dat
  !
  ! ADD parameters
  !
  !N_ED                  = 2         ! Number of extra dimensions in ADD model
  !G_NEWTON                   = 6.707e-39 ! Newton gravity constant
  !M_S                   = 2.5e3     ! string scale for ADD model, meaning depends on KK_CONVENTION
  !M_CUT                 = 2.5e3     ! cut-off scale for the c.m. energy
  !KK_CONVENTION         = 1         ! 0=const mass   
                                    ! 1=simplified sum(HLZ)   2=exact sum(HLZ)
                                    ! 3=Hewett +1    4=Hewett -1 
                                    ! 5=GRW (Lambda_T(GRW)=M_S for virtual Graviton exchange
                                  !        or M_D=M_S for real Gravtion production)
                                  ! only KK_CONVENTION's 1,2 or 5 are applicable for real Graviton production!
  ! set tau unstable
  STABLE[15]=0
}(model)

(me){
  ME_SIGNAL_GENERATOR = Comix
  EVENT_GENERATION_MODE = Unweighted
}(me)

(mi){
  !=================================================!
  !          Underlying Event Setup file            !
  !=================================================!
  !
  ! general parameters
  !
  MI_HANDLER    = None                  ! Amisic / None
  !
  ! hard underlying event parameters
  !
  !CREATE_GRID 93 93 -> 93 93           ! processes to generate
  !PS_ERROR         = 1.0e-2            ! error for integration
  !REGULATE_XS      = 0                 ! regulate cross section
  !XS_REGULATION            = 2.225             ! regulation parameter
  !SCALE_MIN        = 2.225             ! minimum scale
  !RESCALE_EXPONENT    = 0.16           ! rescaling exponent
  !REFERENCE_SCALE          = 1800.0            ! reference energy scale
  !PROFILE_FUNCTION    = Gaussian               ! Gaussian / Double_Gaussian
  !PROFILE_PARAMETERS  = 1.0 0.5 0.5    ! size (must be 1), coresize, 
                                        ! matter fraction
}(mi)

(fragmentation){
  !-------------------------------------------------------
  !-- Fragmentation parameters----------------------------
  !-------------------------------------------------------
  FRAGMENTATION = Ahadic       ! Off, Lund (Pythia string fragmentation) or Ahadic
  DECAYMODEL    = Hadrons      ! Lund or Hadrons
  
  !-------------------------------------------------------
  !-- Photons parameters ---------------------------------
  !-------------------------------------------------------
  YFS_MODE       = 2          ! settings for Photons-module:
                              ! 0        - off/no corrections
                              ! 1        - soft only
                              ! 2        - soft + order(alpha) hard correction (default)
  YFS_USE_ME     = 1          ! use exact MEs for hard correction if possible
  YFS_IR_CUTOFF  = 1E-3       ! IR cut-off
}(fragmentation)