Commit 7d71035f authored by Franziska Iltzsche's avatar Franziska Iltzsche
Browse files

309807 - 309809


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@813551 4525493e-7705-40b1-a816-d608a930855b
parent 9abde7fa
2018-03-28 franziska iltzsche <franziska.speiser@cern.ch>
* add 309807 - 309809
* add common/MadGraph/MadGraphControl_DarkEnergy_jetphiphi_c2_M2000.py
* add common/MadGraph/MadGraphControl_DarkEnergy_jetphiphi_c1_M200.py
* add common/MadGraph/MadGraphControl_DarkEnergy_jetphiphi_c1_M400.py
* add exocontrol/aMcAtNlo_param_card_c1_M200.dat
* add exocontrol/aMcAtNlo_param_card_c1_M400.dat
* add exocontrol/aMcAtNlo_param_card_c2_M2000.dat
* update 399000 - 399057
2018-03-27 Ewelina Maria Lobodzinska <ewelina.maria.lobodzinska@cern.ch>
......
from MadGraphControl.MadGraphUtils import *
#---------------------------------------------------------------------------------------------------
# Set parameters
#---------------------------------------------------------------------------------------------------
lhe_version = 3.0
extras = {'dynamical_scale_choice':3,
'auto_ptj_mjj':'False',
'event_norm':'sum',
'maxjetflavor':'5',
'pdlabel':'lhapdf',
'lhaid':263000,
'use_syst' :'T',
'sys_scalefact' :'0.5 1 2',
'sys_alpsfact' : 'None',
'sys_pdf' : 'NNPDF30_lo_as_0130',
}
parameters = {'YUKAWA':{'invMscale':invMscale,'c1':c1,'c2':c2}}
#---------------------------------------------------------------------------------------------------
# MadGraph process
#---------------------------------------------------------------------------------------------------
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model Standard_Model_cosmo_no_c10_UFO --modelname
generate p p > j phinew phinew NEW=1 QED=0
output -f""")
fcard.close()
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
# evt_multiplier ensure sufficient events generated to be passed to Pythia8 for showering
if evt_multiplier>0:
if runArgs.maxEvents>0:
nevents=runArgs.maxEvents*evt_multiplier
else:
nevents=5000*evt_multiplier
process_dir = new_process()
#---------------------------------------------------------------------------------------------------
# Build run card
#---------------------------------------------------------------------------------------------------
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
#---------------------------------------------------------------------------------------------------
# Build parameter card
#---------------------------------------------------------------------------------------------------
build_param_card(param_card_old='MadGraph_param_card_DarkEnergy.dat',param_card_new='aMcAtNlo_param_card_c1_M200.dat',params=parameters)
print_cards()
mode=0
runName='run_01'
generate(run_card_loc='run_card.dat',param_card_loc='aMcAtNlo_param_card_c1_M200.dat',mode=mode,proc_dir=process_dir,run_name=runName)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=lhe_version)
#---------------------------------------------------------------------------------------------------
# Shower with Pythia8
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8", "EvtGen"]
evgenConfig.description = "jet+MET"
evgenConfig.keywords = ["BSM", "jets"]
evgenConfig.contact = ['Michaela Queitsch-Maitland <michaela.queitsch-maitland@cern.ch>']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
# Phinew is non-interacting
bonus_file = open('pdg_extras.dat','w')
bonus_file.write('1000022 phinew 100.0 (MeV/c) boson phinew 0\n')
bonus_file.close()
testSeq.TestHepMC.G4ExtraWhiteFile='pdg_extras.dat'
# Pythia decays
pythiachans=[]
pythiachans.append("1000022:all phinew phinew 1 0 0 0.1")
pythiachans.append("1000022:mayDecay off")
genSeq.Pythia8.Commands += pythiachans
from MadGraphControl.MadGraphUtils import *
#---------------------------------------------------------------------------------------------------
# Set parameters
#---------------------------------------------------------------------------------------------------
lhe_version = 3.0
extras = {'dynamical_scale_choice':3,
'auto_ptj_mjj':'False',
'event_norm':'sum',
'maxjetflavor':'5',
'pdlabel':'lhapdf',
'lhaid':263000,
'use_syst' :'T',
'sys_scalefact' :'0.5 1 2',
'sys_alpsfact' : 'None',
'sys_pdf' : 'NNPDF30_lo_as_0130',
}
parameters = {'YUKAWA':{'invMscale':invMscale,'c1':c1,'c2':c2}}
#---------------------------------------------------------------------------------------------------
# MadGraph process
#---------------------------------------------------------------------------------------------------
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model Standard_Model_cosmo_no_c10_UFO --modelname
generate p p > j phinew phinew NEW=1 QED=0
output -f""")
fcard.close()
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
# evt_multiplier ensure sufficient events generated to be passed to Pythia8 for showering
if evt_multiplier>0:
if runArgs.maxEvents>0:
nevents=runArgs.maxEvents*evt_multiplier
else:
nevents=5000*evt_multiplier
process_dir = new_process()
#---------------------------------------------------------------------------------------------------
# Build run card
#---------------------------------------------------------------------------------------------------
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
#---------------------------------------------------------------------------------------------------
# Build parameter card
#---------------------------------------------------------------------------------------------------
build_param_card(param_card_old='MadGraph_param_card_DarkEnergy.dat',param_card_new='aMcAtNlo_param_card_c1_M400.dat',params=parameters)
print_cards()
mode=0
runName='run_01'
generate(run_card_loc='run_card.dat',param_card_loc='aMcAtNlo_param_card_c1_M400.dat',mode=mode,proc_dir=process_dir,run_name=runName)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=lhe_version)
#---------------------------------------------------------------------------------------------------
# Shower with Pythia8
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8", "EvtGen"]
evgenConfig.description = "jet+MET"
evgenConfig.keywords = ["BSM", "jets"]
evgenConfig.contact = ['Michaela Queitsch-Maitland <michaela.queitsch-maitland@cern.ch>']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
# Phinew is non-interacting
bonus_file = open('pdg_extras.dat','w')
bonus_file.write('1000022 phinew 100.0 (MeV/c) boson phinew 0\n')
bonus_file.close()
testSeq.TestHepMC.G4ExtraWhiteFile='pdg_extras.dat'
# Pythia decays
pythiachans=[]
pythiachans.append("1000022:all phinew phinew 1 0 0 0.1")
pythiachans.append("1000022:mayDecay off")
genSeq.Pythia8.Commands += pythiachans
from MadGraphControl.MadGraphUtils import *
#---------------------------------------------------------------------------------------------------
# Set parameters
#---------------------------------------------------------------------------------------------------
lhe_version = 3.0
extras = {'dynamical_scale_choice':3,
'auto_ptj_mjj':'False',
'event_norm':'sum',
'maxjetflavor':'5',
'pdlabel':'lhapdf',
'lhaid':263000,
'use_syst' :'T',
'sys_scalefact' :'0.5 1 2',
'sys_alpsfact' : 'None',
'sys_pdf' : 'NNPDF30_lo_as_0130',
}
parameters = {'YUKAWA':{'invMscale':invMscale,'c1':c1,'c2':c2}}
#---------------------------------------------------------------------------------------------------
# MadGraph process
#---------------------------------------------------------------------------------------------------
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model Standard_Model_cosmo_no_c10_UFO --modelname
generate p p > j phinew phinew NEW=1 QED=0
output -f""")
fcard.close()
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
# evt_multiplier ensure sufficient events generated to be passed to Pythia8 for showering
if evt_multiplier>0:
if runArgs.maxEvents>0:
nevents=runArgs.maxEvents*evt_multiplier
else:
nevents=5000*evt_multiplier
process_dir = new_process()
#---------------------------------------------------------------------------------------------------
# Build run card
#---------------------------------------------------------------------------------------------------
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
#---------------------------------------------------------------------------------------------------
# Build parameter card
#---------------------------------------------------------------------------------------------------
build_param_card(param_card_old='MadGraph_param_card_DarkEnergy.dat',param_card_new='aMcAtNlo_param_card_c2_M2000.dat',params=parameters)
print_cards()
mode=0
runName='run_01'
generate(run_card_loc='run_card.dat',param_card_loc='aMcAtNlo_param_card_c2_M2000.dat',mode=mode,proc_dir=process_dir,run_name=runName)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=lhe_version)
#---------------------------------------------------------------------------------------------------
# Shower with Pythia8
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8", "EvtGen"]
evgenConfig.description = "jet+MET"
evgenConfig.keywords = ["BSM", "jets"]
evgenConfig.contact = ['Michaela Queitsch-Maitland <michaela.queitsch-maitland@cern.ch>']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
# Phinew is non-interacting
bonus_file = open('pdg_extras.dat','w')
bonus_file.write('1000022 phinew 100.0 (MeV/c) boson phinew 0\n')
bonus_file.close()
testSeq.TestHepMC.G4ExtraWhiteFile='pdg_extras.dat'
# Pythia decays
pythiachans=[]
pythiachans.append("1000022:all phinew phinew 1 0 0 0.1")
pythiachans.append("1000022:mayDecay off")
genSeq.Pythia8.Commands += pythiachans
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO #####################
######################################################################
###################################
## INFORMATION FOR SMINPUTS
###################################
Block SMINPUTS
1 1.32348905E+02 # aEWM1
2 1.16637000E-05 # Gf
3 1.184000e-01 # aS
###################################
## INFORMATION FOR MASS
###################################
Block MASS
1 5.040000e-03 # MD
2 2.550000e-03 # MU
3 1.010000e-01 # MS
4 1.270000e+00 # MC
5 4.700000e+00 # MB
6 1.725000e+02 # MT
11 5.110000e-04 # Me
13 1.056600e-01 # MMU
15 1.777000e+00 # MTA
23 9.118760e+01 # MZ
25 1.250000e+02 # MH
1000022 1.000000e-01 # mphinew
## Not dependent paramater.
## Those values should be edited following analytical the
## analytical expression. Some generator could simply ignore
## those values and use the analytical expression
22 0.000000 # a : 0.0
24 80.399 # W+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
21 0.000000 # g : 0.0
9000001 0.000000 # ghA : 0.0
9000003 79.824360 # ghWp : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
9000004 79.824360 # ghWm : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
9000005 0.000000 # ghG : 0.0
12 0.000000 # ve : 0.0
14 0.000000 # vm : 0.0
16 0.000000 # vt : 0.0
251 79.824360 # G+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
9000002 91.187600 # ghZ : MZ
250 91.187600 # G0 : MZ
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00
DECAY 23 2.495200e+00
DECAY 24 2.085000e+00
DECAY 25 4.070000e-03
DECAY 1000022 4.070000e-03
## Not dependent paramater.
## Those values should be edited following analytical the
## analytical expression. Some generator could simply ignore
## those values and use the analytical expression
DECAY 22 0.000000 # a : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 9000001 0.000000 # ghA : 0.0
DECAY 9000005 0.000000 # ghG : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 1 0.000000 # d : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 9000002 2.495200 # ghZ : WZ
DECAY 9000003 2.085000 # ghWp : WW
DECAY 9000004 2.085000 # ghWm : WW
DECAY 250 2.495200 # G0 : WZ
DECAY 251 2.085000 # G+ : WW
###################################
## INFORMATION FOR CKMBLOCK
###################################
Block CKMBLOCK
1 2.277360e-01 # cabi
###################################
## INFORMATION FOR YUKAWA
###################################
Block YUKAWA
1 5.040000e-03 # ymdo
2 2.550000e-03 # ymup
3 1.010000e-01 # yms
4 1.270000e+00 # ymc
5 4.700000e+00 # ymb
6 1.725000e+02 # ymt
11 5.110000e-04 # yme
13 1.056600e-01 # ymm
15 1.777000e+00 # ymtau
16 0.005 # invMscale
17 1.000000e-03 # invNscale
18 1.0 # c1
19 0.0 # c2
20 0.000000e+00 # c3a
21 0.000000e+00 # c4a
22 0.000000e+00 # c5a
23 0.000000e+00 # c6a
24 0.000000e+00 # c6b
25 0.000000e+00 # c7
26 0.000000e+00 # c8
27 0.000000e+00 # c9
28 0.000000e+00 # c10a
29 0.000000e+00 # c10b
30 0.000000e+00 # c10c
31 0.000000e+00 # c10d
32 0.000000e+00 # c11a
33 0.000000e+00 # c11b
34 0.000000e+00 # c11c
35 0.000000e+00 # c11d
36 0.000000e+00 # c12a
37 0.000000e+00 # c12b
38 0.000000e+00 # c12c
39 0.000000e+00 # c12d
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================
Block QNUMBERS 9000001 # ghA
1 0 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000002 # ghZ
1 0 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000003 # ghWp
1 3 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000004 # ghWm
1 -3 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000005 # ghG
1 0 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 250 # G0
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 251 # G+
1 3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 1000022 # phinew
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/Antiparticle distinction (0=own anti)
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO #####################
######################################################################
###################################
## INFORMATION FOR SMINPUTS
###################################
Block SMINPUTS
1 1.32348905E+02 # aEWM1
2 1.16637000E-05 # Gf
3 1.184000e-01 # aS
###################################
## INFORMATION FOR MASS
###################################
Block MASS
1 5.040000e-03 # MD
2 2.550000e-03 # MU
3 1.010000e-01 # MS
4 1.270000e+00 # MC
5 4.700000e+00 # MB
6 1.725000e+02 # MT
11 5.110000e-04 # Me
13 1.056600e-01 # MMU
15 1.777000e+00 # MTA
23 9.118760e+01 # MZ
25 1.250000e+02 # MH
1000022 1.000000e-01 # mphinew
## Not dependent paramater.
## Those values should be edited following analytical the
## analytical expression. Some generator could simply ignore
## those values and use the analytical expression
22 0.000000 # a : 0.0
24 80.399 # W+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
21 0.000000 # g : 0.0
9000001 0.000000 # ghA : 0.0
9000003 79.824360 # ghWp : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
9000004 79.824360 # ghWm : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
9000005 0.000000 # ghG : 0.0
12 0.000000 # ve : 0.0
14 0.000000 # vm : 0.0
16 0.000000 # vt : 0.0
251 79.824360 # G+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))
9000002 91.187600 # ghZ : MZ
250 91.187600 # G0 : MZ
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00
DECAY 23 2.495200e+00
DECAY 24 2.085000e+00
DECAY 25 4.070000e-03
DECAY 1000022 4.070000e-03
## Not dependent paramater.
## Those values should be edited following analytical the
## analytical expression. Some generator could simply ignore
## those values and use the analytical expression
DECAY 22 0.000000 # a : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 9000001 0.000000 # ghA : 0.0
DECAY 9000005 0.000000 # ghG : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 1 0.000000 # d : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 9000002 2.495200 # ghZ : WZ
DECAY 9000003 2.085000 # ghWp : WW
DECAY 9000004 2.085000 # ghWm : WW
DECAY 250 2.495200 # G0 : WZ
DECAY 251 2.085000 # G+ : WW
###################################
## INFORMATION FOR CKMBLOCK
###################################
Block CKMBLOCK
1 2.277360e-01 # cabi
###################################
## INFORMATION FOR YUKAWA
###################################
Block YUKAWA
1 5.040000e-03 # ymdo
2 2.550000e-03 # ymup
3 1.010000e-01 # yms
4 1.270000e+00 # ymc
5 4.700000e+00 # ymb
6 1.725000e+02 # ymt
11 5.110000e-04 # yme
13 1.056600e-01 # ymm
15 1.777000e+00 # ymtau
16 0.0025 # invMscale
17 1.000000e-03 # invNscale
18 1.0 # c1
19 0.0 # c2
20 0.000000e+00 # c3a
21 0.000000e+00 # c4a
22 0.000000e+00 # c5a