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 #47: Run.dat

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

Run.dat

Line 
1(run){
2  !-------------------------------------------------------
3  !--- Run parameters ------------------------------------
4  !-------------------------------------------------------
5  EVENTS          = 0             ! Number of Events
6  NUM_ACCURACY   = 1.e-10        !
7  OUTPUT          = 2             ! Output level:
8  !                               ! 0=errors,1=events,2=run info,
9  !                               ! 3=events+information
10  !
11  ANALYSIS        = 1                    ! Analysis: 0=off, 1=on
12  !
13  !
14  !
15  !
16  LOG_FILE        = sherpa_log.log        ! log file
17 
18}(run)
19(me){
20  !-------------------------------------------------------
21  !-- ME generators --------------------------------------
22  !-------------------------------------------------------
23  ME_SIGNAL_GENERATOR   = Amegic      ! Internal or Amegic
24  EVENT_GENERATION_MODE = Weighted
25  !
26  SIGNAL_MODEL          = MSSM        ! subset of a physical model
27                                      ! or the whole SM, MSSM, ADD
28  !
29  COUPLING_SCHEME       = Running_alpha_S ! Fixed   (default value) or
30                                      ! Running (s for ME generators, pt for 2->2)
31  YUKAWA_MASSES         = Fixed       ! Fixed   (polemass) or
32                                      ! Running (higgs mass)
33  YUKAWA_MASSES_FACTOR  = 1.          ! Additional prefactor for yukawas
34  !
35  SCALE_SCHEME          = S_HAT       ! see 'http://www.sherpa-mc.de/scales.html'
36  KFACTOR_SCHEME        = 1           ! 1=as(pt2)/as(ecms)^nstrong,default=1
37  SUDAKOV_WEIGHT        = 0           ! apply sudakov weight on single events
38  SCALE_FACTOR          = 1.          ! factor the scale is multiplied with
39}(me)
40(shower){
41  !-------------------------------------------------------
42  !-- Parton showers -------------------------------------
43  !-------------------------------------------------------
44  SHOWER_GENERATOR  = Apacic   ! 
45  !
46  FSR_SHOWER        = 1        !  1=On,0=Off
47  ISR_SHOWER        = 1        !  1=On,0=Off
48  !
49  IS_PT2MIN         = 4.       !  IS Shower cutoff
50  FS_PT2MIN         = 1.       !  FS Shower cutoff
51}(shower)
52(fragmentation){
53  !-------------------------------------------------------
54  !-- Fragmentation parameters----------------------------
55  !-------------------------------------------------------
56 FRAGMENTATION = Off       ! Off, Lund (Pythia string fragmentation) or Ahadic
57 DECAYMODEL    = Hadrons      ! Lund or Hadrons
58
59  !-------------------------------------------------------
60  !-- Photons parameters ---------------------------------
61  !-------------------------------------------------------
62YFS_MODE       = 2          ! settings for Photons-module:
63                            ! 0        - off/no corrections
64                            ! 1        - soft only
65                            ! 2        - soft + order(alpha) hard correction (default)
66YFS_USE_ME     = 1          ! use exact MEs for hard correction if possible
67YFS_IR_CUTOFF  = 1E-3       ! IR cut-off
68
69}(fragmentation)
70(lund){
71!-------------------------------------------------------
72!--- Lund Run Parameters -------------------------------
73!-------------------------------------------------------
74PARJ(41) = 0.60 ! a
75PARJ(42) = 0.90 ! b
76PARJ(21) = 0.39 ! sigma
77}(lund)
78(model){
79  !-------------------------------------------------------
80  !-- Model parameters -----------------------------------
81  !-------------------------------------------------------
82  MODEL                 = MSSM       ! Model
83  !
84  GENERATOR_ON          = 0        ! use external spectrum generator
85  !
86  ! SM parameters
87  !
88  EW_SCHEME             = 0        ! which parameters define the ew sector.
89  ALPHAS(MZ)            = 0.118    ! strong coupling at scale M_Z
90  ALPHAS(default)       = 0.0800   ! strong coupling
91  ORDER_ALPHAS          = 1        ! NLO
92  1/ALPHAQED(0)         = 137.036  ! inverse of alpha QED in the Thomson limit
93  1/ALPHAQED(default)   = 132.51   ! inverse of alpha QED
94  SIN2THETAW            = 0.2222   ! Weinberg angle at scale M_Z
95  VEV                   = 246.     ! Higgs vev       
96  LAMBDA                = 0.47591  ! SM Higgs self coupling
97  YUKAWA_E              = 0.       ! Yukawa coupling of electron
98  YUKAWA_MU             = 0.105    ! Yukawa coupling of mu
99  YUKAWA_TAU            = 1.77     ! Yukawa coupling of tau
100  YUKAWA_U              = 0.       ! Yukawa coupling of u quark
101  YUKAWA_D              = 0.       ! Yukawa coupling of d quark
102  YUKAWA_S              = 0.       ! Yukawa coupling of s quark
103  YUKAWA_C              = 1.3      ! Yukawa coupling of c quark
104  YUKAWA_B              = 4.5      ! Yukawa coupling of b quark
105  YUKAWA_T              = 175.     ! Yukawa coupling of t quark
106  MSTRANGE_EFF          = 0.48     ! effective strange mass for higgs-loops.
107  CKMORDER              = 0        ! order of expansion of CKM matrix in Cabibbo angle
108  CABIBBO               = 0.22     ! Cabibbo angle (Wolfenstein parametrization)
109  A                     = 0.85     ! A (Wolfenstein parametrization)
110  RHO                   = 0.50     ! rho (Wolfenstein parametrization)
111  ETA                   = 0.50     ! eta (Wolfenstein parametrization)
112  !
113  ! SUSY parameters
114  !
115  SUSY_GENERATOR   = LesHouches
116  SLHA_INPUT       = modelA.out
117  !
118  ! ADD parameters
119  !
120  N_ED                  = 2         ! Number of extra dimensions in ADD model
121  G_NEWTON                    = 6.707e-39 ! Newton gravity constant
122  M_S                   = 2.5e3     ! string scale for ADD model, meaning depends on KK_CONVENTION
123  M_CUT                 = 2.5e3     ! cut-off scale for the c.m. energy
124  KK_CONVENTION         = 1         ! 0=const mass   
125                                    ! 1=simplified sum(HLZ)   2=exact sum(HLZ)
126                                    ! 3=Hewett +1    4=Hewett -1
127                                    ! 5=GRW (Lambda_T(GRW)=M_S for virtual Graviton exchange
128                                  !        or M_D=M_S for real Gravtion production)
129                                  ! only KK_CONVENTION's 1,2 or 5 are applicable for real Graviton production!
130 MASSIVE[5] = 1
131 MASSIVE[15]=1
132 STABLE[15]=1
133 
134 
135}(model)
136(isr){
137  !-------------------------------------------------------
138  !-- ISR parameters ------------------------------------
139  !-------------------------------------------------------
140  BUNCH_1         = 2212      ! possible beam particles:  P+, P-, e+, e-
141  ISR_1           = On        ! On/Off
142  !
143  BUNCH_2         = 2212      ! possible beam particles:  P+, P-, e+, e-
144  ISR_2           = On        ! On/Off
145  !
146  ISR_SMIN        = 1.e-10    ! Minimal fraction of nominal s for parton after ISR
147  ISR_SMAX        = 1.0       ! Maximal fraction of nominal s for parton after ISR
148  !
149  ISR_E_ORDER     = 1         ! Perturbative order of electron structure function
150  ISR_E_SCHEME    = 2         ! Beta-scheme : 0,1,2 , default = 2
151  !
152  PDF_SET         = cteq6l    ! For LHAPDF use one of the following PDFs in
153                              ! in combination with PDFsets
154                              ! a02_lo_v.LHgrid, a02_nlo_v.LHgrid,
155                              ! a02_nnlo_v.LHgrid
156                              ! Alekhin_1000.LHpdf, Alekhin_100.LHpdf
157                              ! Botje_1000.LHpdf, Botje_100.LHpdf
158                              ! cteq4d.LHgrid, cteq4l.LHgrid, cteq4m.LHgrid
159                              ! cteq5d.LHgrid, cteq5l.LHgrid, cteq5m1.LHgrid
160                              ! cteq5m.LHgrid, cteq61.LHgrid, cteq61.LHpdf
161                              ! cteq6.LHpdf, cteq6l.LHpdf, cteq6ll.LHpdf
162                              ! cteq6mE.LHgrid, cteq6m.LHpdf
163                              ! Fermi2002_1000.LHpdf, Fermi2002_100.LHpdf
164                              ! GRV98lo.LHgrid, GRV98nlo.LHgrid
165                              ! H12000disE.LHgrid, H12000dis.LHgrid
166                              ! H12000lo2E.LHgrid, H12000lo2.LHgrid
167                              ! H12000loE.LHgrid, H12000lo.LHgrid
168                              ! H12000msE.LHgrid, H12000ms.LHgrid
169                              ! MRST2001E.LHgrid, MRST2001E.LHpdf
170                              ! MRST2001.LHpdf, MRST2001lo.LHgrid
171                              ! MRST2001nlo.LHgrid. MRST2001nlo.LHpdf
172                              ! MRST2001nnlo.LHgrid, MRST2002nlo.LHgrid
173                              ! MRST2002nlo.LHpdf, MRST2002nnlo.LHgrid
174                              ! MRST2003cnlo.LHgrid, MRST2003cnlo.LHpdf
175                              ! MRST2003cnnlo.LHgrid, MRST98.LHpdf
176                              ! ZEUS2002_FF.LHpdf, ZEUS2002_TR.LHpdf
177                              ! ZEUS2002_ZM.LHpdf
178                              !
179                              ! alternatively use cteq6 via CTEQ6Grid
180                              ! cteq6m,cteq6d,cteq6l,cteq6l1
181                              !
182                              ! Or use interface to MRST99 via MRST99Grid
183                              ! MRST99
184  PDF_SET_VERSION = 1
185  PDF_GRID_PATH   = CTEQ6Grid ! PDFsets,MRST99Grid,CTEQ6Grid
186}(isr)
187(beam){
188  !-------------------------------------------------------
189  !-- Beam parameters ------------------------------------
190  !-------------------------------------------------------
191  BEAM_1          = 2212                 ! possible beam particles:  P+, P-, e+, e-
192  BEAM_ENERGY_1   = 7000.                ! in GeV
193  BEAM_POL_1      = 0.                   ! Polarization degree -1 ... 1
194  BEAM_SPECTRUM_1 = Monochromatic        ! Monochromatic,Laser_Backscattering,Simple_Compton
195  K_PERP_MEAN_1 = 0.2
196  K_PERP_SIGMA_1        = 0.8
197  !
198  BEAM_2          = 2212                 ! possible beam particles:  P+, P-, e+, e-
199  BEAM_ENERGY_2   = 7000.                ! in GeV
200  BEAM_POL_2      = 0.                   ! Polarization degree -1 ... 1
201  BEAM_SPECTRUM_2 = Monochromatic        ! Monochromatic,Laser_Backscattering,Simple_Compton
202  K_PERP_MEAN_2 = 0.2
203  K_PERP_SIGMA_2        = 0.8
204  !
205  BEAM_SMIN       = 1.e-10               ! Minimal fraction of nominal s after beam spectra
206  BEAM_SMAX       = 1.0                  ! Maximal fraction of nominal s after beam spectra
207  !
208  E_LASER_1          = 1.17e-9           ! Laser energy in GeV
209  P_LASER_1          = 0.                ! Laser polarization +-1
210  E_LASER_2          = 1.17e-9           ! Laser energy in GeV
211  P_LASER_2          = 0.                ! Laser polarization +-1
212  LASER_MODE         = 0                 ! 0 = all, 1,2,3 = individual components
213  LASER_ANGLES       = Off               ! On/Off
214  LASER_NONLINEARITY = On                ! On/Off
215}(beam)
216(mi){
217  !=================================================!
218  !          Underlying Event Setup file            !
219  !=================================================!
220  !
221  ! general parameters
222  !
223  MI_HANDLER    = None                  ! Amisic / None
224  !
225  ! hard underlying event parameters
226  !
227  CREATE_GRID 93 93 -> 93 93            ! processes to generate
228  PS_ERROR           = 1.0e-2           ! error for integration
229  REGULATE_XS        = 0                        ! regulate cross section
230  XS_REGULATION      = 2.45             ! regulation parameter
231  SCALE_MIN          = 2.45             ! minimum scale
232  RESCALE_EXPONENT   = 0.16             ! rescaling exponent
233  REFERENCE_SCALE    = 1800.0           ! reference energy scale
234  PROFILE_FUNCTION   = Gaussian         ! Gaussian / Double_Gaussian
235  PROFILE_PARAMETERS = 1.0 0.5 0.5      ! size (must be 1), coresize,
236                                        ! matter fraction
237}(mi)
238(integration){
239  !-------------------------------------------------------
240  !-- Phase space setup ----------------------------------
241  !-------------------------------------------------------
242  ERROR =  1.e-2   ! Error by calculating matrix-elements
243  !     
244  FINISH_OPTIMIZATION = On ! Integrate until optimization is completed     
245  !
246  INTEGRATOR = 6   ! Phasespace : Rambo=0, Sarge=1,
247                   !              Rambo + Sarge=2,
248                 !              Rambo + Multichannel=3,
249                 !              Multichannel:
250                 !              =4 (Original channel generator)
251                 !              =5 (Original+extra channels for competing peaks
252                 !                  like in ->Higgs->ZZ->q qb q qb)
253                 !              =6 (New channel generator)
254                 ! Choice 4-6 only affects channel generation!
255  VEGAS = On
256}(integration)
257(processes){
258%-------------------------------------------------------
259%-- Processes to calculate -----------------------------
260%-------------------------------------------------------
261%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
262%
263!Process : 1000011 -> 90 91 -1000016
264!Print_Graphs
265!End process
266
267
268Process : -1000011 -> 1000016 -15 -12
269Print_Graphs
270End process
271
272Process : 1000011 -> -1000016 15 12
273Print_Graphs
274End process
275
276!Process : 93 93 -> -1000012[a] 1000011[b]
277!Decay : -1000012[a] -> -1000016 90 90
278!Decay : 1000011[b] -> -1000016 90 91
279!Print_Graphs
280!End process
281
282!Process : 93 93 -> -1000014[a] 1000013[b]
283!Decay : -1000014[a] -> -1000016 15 -13   
284!Decay : 1000013[b] -> -1000016 15 14 
285!Print_Graphs
286!End process
287}(processes)
288(selector){
289!-------------------------------------------------------
290!-- Parton level selectors -----------------------------
291!-------------------------------------------------------
292!JetFinder sqr(20/E_CMS) 1.
293
294!PseudoRapidity 93 -6.0 6.0
295!PT 93 10.0 14000.
296
297!PseudoRapidity 15 -5.0 5.0
298!PT 15 0.1 14000.
299!PseudoRapidity -15 -5.0 5.0
300!PT -15 0.1 14000.
301
302!PseudoRapidity 90 -5.0 5.0
303!PT 90 0.1 14000.
304
305
306}(selector)