Pythia8_Base_Fragment.py

Go to the documentation of this file.

from Pythia8_i.Pythia8_iConf import Pythia8_i
genSeq += Pythia8_i("Pythia8")
evgenConfig.generators += ["Pythia8"]
 
if evgenConfig.PDGparams:

   genSeq.Pythia8.Commands += [
       "Main:timesAllowErrors = 500",
       "ParticleDecays:limitTau0 = on",
       "ParticleDecays:tau0Max = 10.0"]
 

   from EvgenProdTools.offline_dict import parameters
 

   particle_params = parameters.get("particles")
   if particle_params:
       for key, value in particle_params.items():
           
           if int(key) in range(6,26):
               genSeq.Pythia8.Commands += [
                   "{:d}:m0 = {}".format(int(key), value['mass']),
                   "{:d}:mWidth = {}".format(int(key), value['width'])
               ]
   else:
       print ("Could not retrieve standard ATLAS particle parameters")
 

   ew_params = parameters.get("EW_parameters")
   if ew_params:
       
       for key,value in ew_params.items():
           if key[1] in ('sin2thetaW', 'sin2thetaWbar'):
               genSeq.Pythia8.Commands += [
                   'StandardModel:{} = {}'.format(key[1], value)
               ]
   else:
       print ("Could not retrieve standard ATLAS EW parameters")
else:

   genSeq.Pythia8.Commands += [
       "Main:timesAllowErrors = 500",
        "6:m0 = 172.5",
        "23:m0 = 91.1876",
        "23:mWidth = 2.4952",
        "24:m0 = 80.399",
        "24:mWidth = 2.085",
        "StandardModel:sin2thetaW = 0.23113",
        "StandardModel:sin2thetaWbar = 0.23146",
        "ParticleDecays:limitTau0 = on",
        "ParticleDecays:tau0Max = 10.0"]
 

from MCJobOptionUtils.JOsupport import check_reset_proc_number
check_reset_proc_number(opts)
 

 
if "PYTHIA8VER" in os.environ:
    verph8str = str(os.environ['PYTHIA8VER'])[:3]
    verph8 = int(verph8str)
    if verph8 >= 314:
        genSeq.Pythia8.Commands += ["TimeShower:recoilStrategyRF = 1"]
 

 
    if verph8 >= 316:    
        genSeq.Pythia8.Commands += ["StringFragmentation:stopMass = 1.0"]