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

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