sherpa is hosted by Hepforge, IPPP Durham
close Warning: Can't synchronize with repository "(default)" (/hepforge/svn/sherpa does not appear to be a Subversion repository.). Look in the Trac log for more information.

Ticket #48: Run-LHC.dat

File Run-LHC.dat, 9.4 KB (added by Generic User (don't modify these fields), 16 years ago)

This RUNDATA crashes as is, but works when changing the process back to 11 -11 93(1)

Line 
1(beam){
2  !-------------------------------------------------------
3  !-- Beam parameters ------------------------------------
4  !-------------------------------------------------------
5  BEAM_1          = 2212                 ! possible beam particles:  P+, P-, e+, e-
6  BEAM_ENERGY_1   = 7000.                ! in GeV
7  BEAM_POL_1      = 0.                   ! Polarization degree -1 ... 1
8  BEAM_SPECTRUM_1 = Monochromatic        ! Monochromatic,Laser_Backscattering,Simple_Compton
9  K_PERP_MEAN_1 = 0.2
10  K_PERP_SIGMA_1        = 0.8
11  !
12  BEAM_2          = 2212                 ! possible beam particles:  P+, P-, e+, e-
13  BEAM_ENERGY_2   = 7000.                ! in GeV
14  BEAM_POL_2      = 0.                   ! Polarization degree -1 ... 1
15  BEAM_SPECTRUM_2 = Monochromatic        ! Monochromatic,Laser_Backscattering,Simple_Compton
16  K_PERP_MEAN_2 = 0.2
17  K_PERP_SIGMA_2        = 0.8
18  !
19  BEAM_SMIN       = 1.e-10               ! Minimal fraction of nominal s after beam spectra
20  BEAM_SMAX       = 1.0                  ! Maximal fraction of nominal s after beam spectra
21  !
22  E_LASER_1          = 1.17e-9           ! Laser energy in GeV
23  P_LASER_1          = 0.                ! Laser polarization +-1
24  E_LASER_2          = 1.17e-9           ! Laser energy in GeV
25  P_LASER_2          = 0.                ! Laser polarization +-1
26  LASER_MODE         = 0                 ! 0 = all, 1,2,3 = individual components
27  LASER_ANGLES       = Off               ! On/Off
28  LASER_NONLINEARITY = On                ! On/Off
29}(beam)
30(isr){
31  !-------------------------------------------------------
32  !-- ISR parameters ------------------------------------
33  !-------------------------------------------------------
34  BUNCH_1         = 2212      ! possible beam particles:  P+, P-, e+, e-
35  ISR_1           = On        ! On/Off
36  !
37  BUNCH_2         = 2212      ! possible beam particles:  P+, P-, e+, e-
38  ISR_2           = On        ! On/Off
39  !
40  ISR_SMIN        = 1.e-10    ! Minimal fraction of nominal s for parton after ISR
41  ISR_SMAX        = 1.0       ! Maximal fraction of nominal s for parton after ISR
42  !
43  ISR_E_ORDER     = 1         ! Perturbative order of electron structure function
44  ISR_E_SCHEME    = 2         ! Beta-scheme : 0,1,2 , default = 2
45  !
46  PDF_SET         = cteq6l    ! the CTEQ6 PDFs cteq6m,cteq6d,cteq6l,cteq6l1
47                            ! can be accessed when PDF_GRID_PATH = CTEQ6Grid
48                            !
49                            ! to use LHAPDF list one of the PDF sets,
50                            ! , e.g. Alekhin_100.LHpdf, that comes with your
51                            ! LHAPDF installation, the full list can be found in
52                              ! 'lhapdf-config --pdfsets-path', in addition set
53                            ! PDF_GRID_PATH = PDFsets
54                            !
55                              ! or use the interface to MRST99 through
56                            ! PDF_GRID_PATH = MRST99Grid
57                              !
58  PDF_SET_VERSION = 1
59  !
60  PDF_GRID_PATH   = CTEQ6Grid ! CTEQ6Grid,PDFsets,MRST99Grid
61}(isr)
62(model){
63  !-------------------------------------------------------
64  !-- Model parameters -----------------------------------
65  !-------------------------------------------------------
66  MODEL                 = SM       ! Model
67  !
68  ! SM parameters
69  !
70  EW_SCHEME             = 0        ! which parameters define the ew sector.
71  ALPHAS(MZ)            = 0.118    ! strong coupling at scale M_Z
72  ALPHAS(default)       = 0.0800   ! strong coupling
73  ORDER_ALPHAS          = 1        ! NLO
74  1/ALPHAQED(0)         = 137.036  ! inverse of alpha QED in the Thomson limit
75  1/ALPHAQED(default)   = 132.51   ! inverse of alpha QED
76  SIN2THETAW            = 0.2222   ! Weinberg angle at scale M_Z
77  VEV                   = 246.     ! Higgs vev       
78  LAMBDA                = 0.47591  ! SM Higgs self coupling
79  CKMORDER              = 0        ! order of expansion of CKM matrix in Cabibbo angle
80  CABIBBO               = 0.22     ! Cabibbo angle (Wolfenstein parametrization)
81  A                     = 0.85     ! A (Wolfenstein parametrization)
82  RHO                   = 0.50     ! rho (Wolfenstein parametrization)
83  ETA                   = 0.50     ! eta (Wolfenstein parametrization)
84  !
85  ! SUSY parameters
86  !
87  SLHA_INPUT            = LesHouches_SPS1A.dat
88  !
89  ! ADD parameters
90  !
91  N_ED                  = 2         ! Number of extra dimensions in ADD model
92  G_NEWTON                    = 6.707e-39 ! Newton gravity constant
93  M_S                   = 2.5e3     ! string scale for ADD model, meaning depends on KK_CONVENTION
94  M_CUT                 = 2.5e3     ! cut-off scale for the c.m. energy
95  KK_CONVENTION         = 1         ! 0=const mass   
96                                    ! 1=simplified sum(HLZ)   2=exact sum(HLZ)
97                                    ! 3=Hewett +1    4=Hewett -1
98                                    ! 5=GRW (Lambda_T(GRW)=M_S for virtual Graviton exchange
99                                  !        or M_D=M_S for real Gravtion production)
100                                  ! only KK_CONVENTION's 1,2 or 5 are applicable for real Graviton production!
101  ! set tau unstable
102  STABLE[15]=0
103}(model)
104(me){
105  !-------------------------------------------------------
106  !-- ME generators --------------------------------------
107  !-------------------------------------------------------
108  ME_SIGNAL_GENERATOR   = Amegic      ! Internal or Amegic
109  EVENT_GENERATION_MODE = Unweighted
110  !
111  COUPLING_SCHEME       = Running_alpha_S ! Fixed   (default value) or
112                                      ! Running (s for ME generators, pt for 2->2)
113  YUKAWA_MASSES         = Fixed       ! Fixed   (polemass) or
114                                      ! Running (higgs mass)
115  YUKAWA_MASSES_FACTOR  = 1.          ! Additional prefactor for yukawas
116  !
117  SCALE_SCHEME          = CKKW        ! see 'http://www.sherpa-mc.de/scales.html'
118  KFACTOR_SCHEME        = 1           ! 1=as(pt2)/as(ecms)^nstrong,default=1
119  SUDAKOV_WEIGHT        = 1           ! apply sudakov weight on single events
120  SCALE_FACTOR          = 1.          ! factor the scale is multiplied with
121}(me)
122(integration){
123  !-------------------------------------------------------
124  !-- Phase space setup ----------------------------------
125  !-------------------------------------------------------
126  ERROR =  1.e-2   ! Error by calculating matrix-elements
127  !     
128  FINISH_OPTIMIZATION = On ! Integrate until optimization is completed     
129  !
130  INTEGRATOR = 6   ! Phasespace : Rambo=0, Sarge=1,
131                   !              Rambo + Sarge=2,
132                 !              Rambo + Multichannel=3,
133                 !              Multichannel:
134                 !              =4 (Original channel generator)
135                 !              =5 (Original+extra channels for competing peaks
136                 !                  like in ->Higgs->ZZ->q qb q qb)
137                 !              =6 (New channel generator)
138                 ! Choice 4-6 only affects channel generation!
139  VEGAS = On
140}(integration)
141(mi){
142  !=================================================!
143  !          Underlying Event Setup file            !
144  !=================================================!
145  !
146  ! general parameters
147  !
148  MI_HANDLER    = None                  ! Amisic / None
149  !
150  ! hard underlying event parameters
151  !
152  CREATE_GRID 93 93 -> 93 93            ! processes to generate
153  PS_ERROR          = 1.0e-2            ! error for integration
154  REGULATE_XS       = 0                 ! regulate cross section
155  XS_REGULATION     = 2.225             ! regulation parameter
156  SCALE_MIN         = 2.225             ! minimum scale
157  RESCALE_EXPONENT    = 0.16            ! rescaling exponent
158  REFERENCE_SCALE           = 1800.0            ! reference energy scale
159  PROFILE_FUNCTION    = Gaussian                ! Gaussian / Double_Gaussian
160  PROFILE_PARAMETERS  = 1.0 0.5 0.5     ! size (must be 1), coresize,
161                                        ! matter fraction
162}(mi)
163(shower){
164  !-------------------------------------------------------
165  !-- Parton showers -------------------------------------
166  !-------------------------------------------------------
167  SHOWER_GENERATOR  = Apacic   ! 
168  !
169  FSR_SHOWER        = 1        !  1=On,0=Off
170  ISR_SHOWER        = 1        !  1=On,0=Off
171  !
172  IS_PT2MIN         = 4.       !  IS Shower cutoff
173  FS_PT2MIN         = 1.       !  FS Shower cutoff
174}(shower)
175(fragmentation){
176  !-------------------------------------------------------
177  !-- Fragmentation parameters----------------------------
178  !-------------------------------------------------------
179  FRAGMENTATION = Ahadic       ! Off, Lund (Pythia string fragmentation) or Ahadic
180  DECAYMODEL    = Hadrons      ! Lund or Hadrons
181 
182  !-------------------------------------------------------
183  !-- Photons parameters ---------------------------------
184  !-------------------------------------------------------
185  YFS_MODE       = 2          ! settings for Photons-module:
186                              ! 0        - off/no corrections
187                              ! 1        - soft only
188                              ! 2        - soft + order(alpha) hard correction (default)
189  YFS_USE_ME     = 1          ! use exact MEs for hard correction if possible
190  YFS_IR_CUTOFF  = 1E-3       ! IR cut-off
191}(fragmentation)
192(processes){
193  %-------------------------------------------------------
194  %-- Processes to calculate -----------------------------
195  %-------------------------------------------------------
196  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
197  %
198  % jet jet -> electron positron + n-jets
199  %
200  % Process : 93 93 -> 11 -11 93{1}
201  Process : 93 93 -> 22 23 93{1}
202  Order electroweak : 2
203  End process
204  %
205}(processes)
206(selector){
207  !-------------------------------------------------------
208  !-- Parton level selectors -----------------------------
209  !-------------------------------------------------------
210  JetFinder   sqr(20/E_CMS) 1.
211  Mass 11 -11 66 116
212}(selector)