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-gammabb-z-analysis.dat

File Run-gammabb-z-analysis.dat, 9.9 KB (added by Generic User (don't modify these fields), 15 years ago)

This RUNDATA used to work for us in 1.1.2

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