Skip to content
Snippets Groups Projects
Commit 9aa9fc88 authored by Sylvester Joosten's avatar Sylvester Joosten
Browse files

streamlined generator settings

parent e0138cde
No related branches found
No related tags found
No related merge requests found
...@@ -4,31 +4,29 @@ ...@@ -4,31 +4,29 @@
#ifdef __cplusplus #ifdef __cplusplus
extern "C" { extern "C" {
#endif #endif
struct GeneratorSettings { struct GeneratorSettings {
int photoStruc; // photon structure parameter (MSTP(14)) int PhotoStruc; // photon structure parameter (MSTP(14))
int photoType; // user control switch (MSEL, 9.2) int PhotoType; // user control switch (MSEL, 9.2)
int cosNPhiModel; // disable with 0, 1 for HERMES hadron parameterization int PyMaxFl; // PYTHIA max flavor
int pyMaxFlavor; // PYTHIA max flavor int PartonSet; // PDF set for PYTHIA/JETSET
int PDFSet; // PDF set for PYTHIA/JETSET int DecayPiK; // Have LUND decay pi+-, K+- and K0_L
int decayPiK; // Have LUND decay pi+-, K+- and K0_L int DecayLamKs; // Have LUND decay Lambda, Sigma, Xi, Omega, K0_S
int decayLamKs; // Have LUND decay Lambda, Sigma, Xi, Omega, K0_S int radgenGenLUT; // generate LUT for pol. rad. corr.
int radgenGenLUT; // generate LUT for pol. rad. corr. int radgenUseLUT; // use LUT for pol. rad. corr.
int radgenUseLUT; // use LUT for pol. rad. corr. int enableFM; // activate or not the FM
int enableFM; // activate or not the FM int FMNBins; // Number of bin for FM discretization
int FMNBins; // Number of bin for FM discretization float FMCutOff; // Max FM (in GeV)
float FMCutOff; // Max FM (in GeV) char PyModel[4]; // l/N (DIS) or gamma/l pN (GVMD,REAL,DIS,RAD)
char pyModel[4]; // l/N (DIS) or gamma/l pN (GVMD,REAL,DIS,RAD) char GenRad[4]; // radiative effects (NO,POL)
char genRad[4]; // radiative effects (NO,POL) char FStruct[4]; // structure function parameterization
char paramF[4]; // structure function parameterization char R[4]; // param. of R=sL/sT (1990, 1998; or 0 for R=0)
char paramR[4]; // param. of R=sL/sT (1990, 1998; or 0 for R=0) };
}; extern struct GeneratorSettings common_gen_set_;
extern struct GeneratorSettings common_gen_set_;
#define glb_gen common_gen_set_ #define glb_gen common_gen_set_
#ifdef __cplusplus #ifdef __cplusplus
} }
#endif #endif
#endif #endif
common /common_gen_set/ common /common_gen_set/
+ genSet_PhotoStruc,
+ genSet_PhotoType, + genSet_PhotoType,
+ genSet_CosNPhiParam,
+ genSet_PyMaxFl,
+ genSet_PartonSet,
+ genSet_DecayPiK, + genSet_DecayPiK,
+ genSet_DecayLamKs, + genSet_DecayLamKs,
+ genSet_GenLUT, + genSet_GenLUT,
...@@ -18,11 +14,7 @@ ...@@ -18,11 +14,7 @@
save /common_gen_set/ save /common_gen_set/
integer integer
+ genSet_PhotoStruc,
+ genSet_PhotoType, + genSet_PhotoType,
+ genSet_CosNPhiParam,
+ genSet_PyMaxFl,
+ genSet_PartonSet,
+ genSet_fastMC, + genSet_fastMC,
+ genSet_enableFM, + genSet_enableFM,
+ genSet_FMNBins + genSet_FMNBins
......
...@@ -43,31 +43,12 @@ ...@@ -43,31 +43,12 @@
! Other parameters ! Other parameters
! ... structue of incoming photon (resolved, VMD) (sec 9.3)
mstp(14) = genSet_PhotoStruc
! ... select nucleon parton distribution set
! ... GNOME? mstp(56) is photo PDF
! TODO move to after gmc_loadtune
mstp(51) = genSet_PartonSet
if (genSet_PartonSet.lt.1000) then
mstp(52) = 1
else
mstp(52) = 2
endif
! ... number of allowed flavours (sec 6.3.3 ... e+e- only?)
mstp(58) = genSet_PyMaxFl
! ......... (e+e- only?) ! ......... (e+e- only?)
mstj(104) = genSet_PyMaxFl C mstj(104) = genSet_PyMaxFl
mstj(104) = 4
! ... number of quark generations ! ... number of quark generations
mstp(1) = genSet_PyMaxFl C mstp(1) = genSet_PyMaxFl
mstp(1) = 4
! ... QCD and QED radiative corrections (sec 10.4; 1=on=default, 0=off)
! .... according to me (E.C.A) we should run with mstp(61) = 0 as that turns
! .... of parton showering
mstp(61) = 0
mstp(71) = 0
! ... parameters describing fragmentation functions ! ... parameters describing fragmentation functions
call gmc_setfrag (0) call gmc_setfrag (0)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment