def getWhizardModel(self, nbevts, energy, model):
""" Create a default whizard
"""
from ILCDIRAC.Interfaces.API.NewInterface.Applications import Whizard
proddict = "e2e2_o"
if model != "sm":
proddict = "se2se2_r"
whiz = Whizard(processlist=self.diracInstance.getProcessList())
whiz.setModel(model)
pdict = {}
pdict['process_input'] = {}
pdict['process_input']['process_id'] = proddict
pdict['process_input']['sqrts'] = energy
if model != 'sm':
pdict['process_input']['input_file'] = "LesHouches.msugra_1"
pdict['process_input']['input_slha_format'] = 'T'
pdict['process_input']['beam_recoil'] = 'T'
pdict['integration_input'] = {}
pdict['integration_input']['calls'] = '1 50000 10 5000 1 15000'
pdict['simulation_input'] = {}
pdict['simulation_input']['normalize_weight'] = 'F'
pdict['simulation_input']['n_events'] = nbevts
pdict['simulation_input']['keep_initials'] = 'T'
pdict['simulation_input']['events_per_file'] = 500000
if model != 'sm':
pdict['simulation_input'][
'pythia_parameters'] = "PMAS(25,1)=125; PMAS(25,2)=0.3605E-02; MSTU(22)=20 ;PARJ(21)=0.40000;PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000;PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;IMSS(1)=11; IMSS(21)=71; IMSS(22)=71"
else:
pdict['simulation_input'][
'pythia_parameters'] = "PMAS(25,1)=125; PMAS(25,2)=0.3605E-02; MSTU(22)=20 ; MSTJ(28)=2 ;PARJ(21)=0.40000;PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000;PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1"
pdict['parameter_input'] = {}
# pdict['parameter_input']['mmu']=mmu
# pdict['parameter_input']['mtau']=mtau
# pdict['parameter_input']['mb']=mb
# pdict['parameter_input']['mc']=mc
pdict['parameter_input']['mH'] = 125
pdict['beam_input_1'] = {}
pdict['beam_input_1']['particle_name'] = "e1"
pdict['beam_input_1']['polarization'] = "0.0 0.0"
pdict['beam_input_1']['USER_spectrum_on'] = 'T'
pdict['beam_input_1'][
'USER_spectrum_mode'] = 19 if energy == 1400 else 11
pdict['beam_input_1']['ISR_on'] = 'T'
pdict['beam_input_1']['EPA_on'] = "F"
pdict['beam_input_2'] = {}
pdict['beam_input_2']['particle_name'] = "E1"
pdict['beam_input_2']['polarization'] = "0.0 0.0"
pdict['beam_input_2']['USER_spectrum_on'] = 'T'
pdict['beam_input_2']['ISR_on'] = 'T'
pdict['beam_input_2'][
'USER_spectrum_mode'] = 19 if energy == 1400 else 11
pdict['beam_input_2']['EPA_on'] = 'F'
whiz.setFullParameterDict(pdict)
whiz.setOutputFile("testgen.stdhep")
return whiz
def getWhizardModel(self, nbevts, energy, model):
""" Create a default whizard
"""
from ILCDIRAC.Interfaces.API.NewInterface.Applications import Whizard
proddict = "e2e2_o"
if model != "sm":
proddict = "se2se2_r"
whiz = Whizard(processlist = self.diracInstance.getProcessList())
whiz.setModel(model)
pdict = {}
pdict['process_input'] = {}
pdict['process_input']['process_id'] = proddict
pdict['process_input']['sqrts'] = energy
if model != 'sm':
pdict['process_input']['input_file'] = "LesHouches.msugra_1"
pdict['process_input']['input_slha_format'] = 'T'
pdict['process_input']['beam_recoil'] = 'T'
pdict['integration_input'] = {}
pdict['integration_input']['calls'] = '1 50000 10 5000 1 15000'
pdict['simulation_input'] = {}
pdict['simulation_input']['normalize_weight'] = 'F'
pdict['simulation_input']['n_events'] = nbevts
pdict['simulation_input']['keep_initials'] = 'T'
pdict['simulation_input']['events_per_file'] = 500000
if model != 'sm':
pdict['simulation_input']['pythia_parameters'] = "PMAS(25,1)=125; PMAS(25,2)=0.3605E-02; MSTU(22)=20 ;PARJ(21)=0.40000;PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000;PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;IMSS(1)=11; IMSS(21)=71; IMSS(22)=71"
else:
pdict['simulation_input']['pythia_parameters'] = "PMAS(25,1)=125; PMAS(25,2)=0.3605E-02; MSTU(22)=20 ; MSTJ(28)=2 ;PARJ(21)=0.40000;PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000;PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1"
pdict['parameter_input'] = {}
# pdict['parameter_input']['mmu']=mmu
# pdict['parameter_input']['mtau']=mtau
# pdict['parameter_input']['mb']=mb
# pdict['parameter_input']['mc']=mc
pdict['parameter_input']['mH'] = 125
pdict['beam_input_1'] = {}
pdict['beam_input_1']['particle_name'] = "e1"
pdict['beam_input_1']['polarization'] = "0.0 0.0"
pdict['beam_input_1']['USER_spectrum_on'] = 'T'
pdict['beam_input_1']['USER_spectrum_mode'] = 19 if energy == 1400 else 11
pdict['beam_input_1']['ISR_on'] = 'T'
pdict['beam_input_1']['EPA_on'] = "F"
pdict['beam_input_2'] = {}
pdict['beam_input_2']['particle_name'] = "E1"
pdict['beam_input_2']['polarization'] = "0.0 0.0"
pdict['beam_input_2']['USER_spectrum_on'] = 'T'
pdict['beam_input_2']['ISR_on'] = 'T'
pdict['beam_input_2']['USER_spectrum_mode'] = 19 if energy == 1400 else 11
pdict['beam_input_2']['EPA_on'] = 'F'
whiz.setFullParameterDict(pdict)
whiz.setOutputFile("testgen.stdhep")
return whiz
