CalcHEP

The batch-file edited by myself don't work in the model file of NMSSMTools

Asked by tang yilei

I just downloaded a NMSSMTools and manipulate with the model files listed in microRun/work/*.mdi. The batch file looks like this,
#Begin batch file
Model: NMSSM(micro)
Model changed: True
Gauge: Feynman

Process: ~o1,~o1->AllEven,AllEven
Composite: AllEven=A,Z,e,E,m,M

p1: 0.91965
p2: 0.91965

Parameter: EE=0.312853
Parameter: SW=0.482000
Parameter: MZ=91.200000
Parameter: GG=1.117
Parameter: alfSMZ=0.117200
Parameter: MbMb=4.214000
Parameter: McMc=1.300000
Parameter: Mtp=173.100000
Parameter: Ml=1.777000
Parameter: wt=1.442000
Parameter: Mq=0.050000
Parameter: MG1=500.000000
Parameter: MG2=1000.000000
Parameter: MG3=3000.000000
Parameter: tb=0.000000
Parameter: mu=345.113567
Parameter: Lambda=0.271694
Parameter: Kappa=0.035374
Parameter: aLambda=2000.000000
Parameter: aKappa=0.000000
Parameter: At=1500.000000
Parameter: Ab=1500.000000
Parameter: Al=1500.000000
Parameter: Ml2=200.000000
Parameter: Ml3=200.000000
Parameter: Mr2=200.000000
Parameter: Mr3=200.000000
Parameter: Mq2=1000.000000
Parameter: Mq3=1000.000000
Parameter: Mu2=1000.000000
Parameter: Mu3=1000.000000
Parameter: Md2=1000.000000
Parameter: Md3=1000.000000
Parameter: Mh1=83.501298
Parameter: Mh2=127.639686
Parameter: Mh3=1929.035320
Parameter: Mha=21.552419
Parameter: Mhb=1928.774487
Parameter: MHc=1929.187951
Parameter: MNE1=91.965364
Parameter: MNE2=339.275457
Parameter: MNE3=-356.744288
Parameter: MNE4=500.358064
Parameter: MNE5=1010.191822
Parameter: MC1=345.898075
Parameter: MC2=1010.175296
Parameter: MSuL=824.129765
Parameter: MSuR=825.036208
Parameter: MSdL=827.797344
Parameter: MSdR=826.151397
Parameter: MScL=824.129765
Parameter: MScR=825.036208
Parameter: MSsL=827.797344
Parameter: MSsR=826.151397
Parameter: MSt1=645.382799
Parameter: MSt2=993.263600
Parameter: MSb1=826.168454
Parameter: MSb2=827.787660
Parameter: MSeL=204.875686
Parameter: MSeR=204.264199
Parameter: MSmL=204.875686
Parameter: MSmR=204.264199
Parameter: MSl1=203.259866
Parameter: MSl2=205.887474
Parameter: MSne=190.531074
Parameter: MSnm=190.531074
Parameter: MSnl=190.531074
Parameter: Zn11=0.014638
Parameter: Zn12=-0.012030
Parameter: Zn13=0.014820
Parameter: Zn14=-0.135764
Parameter: Zn15=0.990449
Parameter: Zn21=0.222012
Parameter: Zn22=-0.095529
Parameter: Zn23=0.696021
Parameter: Zn24=-0.667700
Parameter: Zn25=-0.106380
Parameter: Zn31=-0.028332
Parameter: Zn32=0.033100
Parameter: Zn33=0.702851
Parameter: Zn34=0.704664
Parameter: Zn35=0.086895
Parameter: Zn41=0.974485
Parameter: Zn42=0.031546
Parameter: Zn43=-0.137942
Parameter: Zn44=0.173822
Parameter: Zn45=0.011871
Parameter: Zn51=-0.008467
Parameter: Zn52=0.994303
Parameter: Zn53=0.048029
Parameter: Zn54=-0.094765
Parameter: Zn55=-0.001507
Parameter: Zu11=0.067584
Parameter: Zu12=-0.997714
Parameter: Zu21=0.997714
Parameter: Zu22=0.067584
Parameter: Zv11=0.133798
Parameter: Zv12=-0.991009
Parameter: Zv21=0.991009
Parameter: Zv22=0.133798
Parameter: Zt11=-0.708100
Parameter: Zt12=0.706112
Parameter: Zt21=-0.706112
Parameter: Zt22=-0.708100
Parameter: Zb11=0.043687
Parameter: Zb12=0.999045
Parameter: Zb21=-0.999045
Parameter: Zb22=0.043687
Parameter: Zl11=0.619390
Parameter: Zl12=0.785083
Parameter: Zl21=-0.785083
Parameter: Zl22=0.619390
Parameter: Zh11=0.113950
Parameter: Zh12=0.478659
Parameter: Zh13=0.870575
Parameter: Zh21=0.153624
Parameter: Zh22=0.857254
Parameter: Zh23=-0.491443
Parameter: Zh31=0.981537
Parameter: Zh32=-0.189741
Parameter: Zh33=-0.024150
Parameter: Za11=-0.023287
Parameter: Za12=-0.004479
Parameter: Za13=0.999719
Parameter: Za21=0.981723
Parameter: Za22=0.188835
Parameter: Za23=0.023714
Parameter: la1=0.117365
Parameter: la2=0.247896
Parameter: la3=0.057439
Parameter: la4=-0.121916
Parameter: la5=-0.000623
Parameter: la6=0.000435
Parameter: la7=0.002597
Parameter: la1s=0.072002
Parameter: la2s=0.072004
Parameter: la3s=0.001393
Parameter: la4s=-0.000003
Parameter: la5s=0.000000
Parameter: la6s=0.009392
Parameter: la7s=0.000000
Parameter: la8s=0.000000
Parameter: aa1=-0.000006
Parameter: aa2=-0.000006
Parameter: aa3=0.000000
Parameter: aa4=0.000000
Parameter: aa5=529.523208
Parameter: aa6=0.000000
Parameter: B1=0.000000
Parameter: B2=0.000000
Parameter: X=0.000000
Parameter: dMb=0.036955
Parameter: Q=183.930728
Parameter: Au=0.000000
Parameter: Ad=0.000000
Parameter: Maux=1.000000
Parameter: tB=5.198849
Parameter: MSG=2882.240418
Parameter: wNE1=0.919654

Number of events (per run step): 10000
Filename: test
NTuple: False
Cleanup: False

Parallelization method: local
Max number of cpus: 2
sleep time: 3
nice level : 19

Dist parameter: cos(p1,p3)
Dist min: 0
Dist max: 1
Dist n bins: 10
#Dist title:
#Dist x-title:

nSess_1: 5
nCalls_1: 100000
nSess_2: 5
nCalls_2: 100000
#end batch file

Then it always report error that xvev cannot be calculated. I think that's because the default value in vars1.mdl causes 0/0 when n_calchep is launched, so I modified the vars1.mdl into

#begin vars1.mdl
NMSSM(micro)
 Variables
 Name | Value |> Comment <|
EE |0.312853 |Electromagnetic coupling constant
SW |0.482 |sin of the Weinberg angle
MZ |91.2 |Z mass
GG |1.117 |Strong coupling constant
alfSMZ |0.1184 |
MbMb |4.5 |Mb(Mb)
McMc |1.3 |Mc(Mc)
Mtp |174.5 |pole mass of t-quark
Ml |1.777 |mass of tau-lepton
wt |1.442 |width of top quark
Mq |1.15 |mass for light quarks
MG1 |1 |
MG2 |1 |
MG3 |1 |
tb |0 |
mu |1 |
Lambda |1 |
Kappa |1 |
aLambda|1 |
aKappa |1 |
At |1 |
Ab |1 |
Al |1 |
Ml2 |211 |ML2
Ml3 |211 |ML3
Mr2 |212 |MR2
Mr3 |211 |MR3
Mq2 |1111 |
Mq3 |1111 |
Mu2 |1111 |
Mu3 |1111 |
Md2 |1111 |
Md3 |1111 |
Mh1 |11 |
Mh2 |12 |
Mh3 |13 |
Mha |4 |
Mhb |5 |
MHc |6 |
MNE1 |7 |
MNE2 |8 |
MNE3 |9 |
MNE4 |0.1 |
MNE5 |0.2 |
MC1 |1 |
MC2 |1 |
MSuL |1 |
MSuR |2 |
MSdL |3 |
MSdR |4 |
MScL |5 |
MScR |6 |
MSsL |7 |
MSsR |8 |
MSt1 |2 |
MSt2 |3 |
MSb1 |4 |
MSb2 |5 |
MSeL |6 |
MSeR |7 |
MSmL |8 |
MSmR |9 |
MSl1 |1.1 |
MSl2 |1.2 |
MSne |5 |
MSnm |7 |
MSnl |1 |
Zn11 |1 |
Zn12 |0 |
Zn13 |0 |
Zn14 |0 |
Zn15 |0 |
Zn21 |0 |
Zn22 |1 |
Zn23 |0 |
Zn24 |0 |
Zn25 |0 |
Zn31 |0 |
Zn32 |0 |
Zn33 |1 |
Zn34 |0 |
Zn35 |0 |
Zn41 |0 |
Zn42 |0 |
Zn43 |0 |
Zn44 |1 |
Zn45 |0 |
Zn51 |0 |
Zn52 |0 |
Zn53 |0 |
Zn54 |0 |
Zn55 |1 |
Zu11 |1 |
Zu12 |0 |
Zu21 |0 |
Zu22 |1 |
Zv11 |1 |
Zv12 |0 |
Zv21 |0 |
Zv22 |1 |
Zt11 |1 |
Zt12 |0 |
Zt21 |0 |
Zt22 |1 |
Zb11 |1 |
Zb12 |0 |
Zb21 |0 |
Zb22 |1 |
Zl11 |1 |
Zl12 |0 |
Zl21 |0 |
Zl22 |1 |
Zh11 |1 |
Zh12 |0 |
Zh13 |0 |
Zh21 |0 |
Zh22 |1 |
Zh23 |0 |
Zh31 |0 |
Zh32 |0 |
Zh33 |1 |
Za11 |1 |
Za12 |0 |
Za13 |0 |
Za21 |0 |
Za22 |1 |
Za23 |0 |
la1 |1 |
la2 |1 |
la3 |1 |
la4 |1 |
la5 |1 |
la6 |1 |
la7 |1 |
la1s |1 |
la2s |1 |
la3s |1 |
la4s |1 |
la5s |1 |
la6s |1 |
la7s |1 |
la8s |1 |
aa1 |1 |
aa2 |1 |
aa3 |1 |
aa4 |1 |
aa5 |1 |
aa6 |1 |
B1 |1 |
B2 |1 |
X |1 |
dMb |1 |
Q |111 |QCD scale for quark masses
Au |1 |
Ad |1 |
Maux |1 |mass of aux particles
tB |2 |
MSG |1 |
wNE1 |1 | neutralino width
#end vars1.mdl

  I guess that the calchep_batch will calculate the functions in fun1 before it loads the parameter values in the batch file, and then calculate the functions again(am I right?). However, new errors arised. Here is the error message:
calchep_batch version 1.6

Processing batch:
Progress information can be found in the html directory.
Simply open the following link in your browser:
file:///home/tangyilei/OrigNMSSMTools/NMSSMTools_4.3.0/microRun/work/html/index.html
You can also view textual progress reports in /home/tangyilei/OrigNMSSMTools/NMSSMTools_4.3.0/microRun/work/html/index.txt
 and the other .txt files in the html directory.
Events will be stored in the Events directory.

Segmentation fault (core dumped)
mv: 无法获取"/home/tangyilei/OrigNMSSMTools/NMSSMTools_4.3.0/microRun/work/Processes/m1/Feynman/p3/single/session.dat" 的文件状态(stat): 没有那个文件或目录
Error in calculation of cross section for ~o1,~o1->Z,Z
Please have a look at cs.o and cs.out in Processes/m1/Feynman/p4/single.
Killing all jobs and quitting.

The Chinese sentences probably mean, mv: cannot acquire the status(stat) of blablabla, no such file or directory.
What's the matter here?

Question information

Language:
English Edit question
Status:
Solved
For:
CalcHEP Edit question
Assignee:
No assignee Edit question
Solved by:
Alexander Belyaev
Solved:
Last query:
Last reply:
Revision history for this message
Alexander Pukhov (pukhov) said : 		#1

Sorry for delay.
Models in microRun/work/models are not intended for CalcHEP batch and interactive sessions. In version NMSSM(micro) all parameters are assigned by other program. It allows to microRun/main to work with different scenarios. But it is difficult to assign them 'by hand'.

The models which one should use for CalcHEP sessions ( both batch and interactive) are presented in directory microRun/calchep. Just go to this directory and launch ./calchep or ./calchep_batch.

If you would like to organize some cycles for calculation of relic and other Dark Matter observables, you don't need calchep_batch. microRun/main.c allows you to do it.

I just put on CalcHEP web site new version of NMSSM working with NMSSMTools_4.4.0. You can use it for calchep_batch.

The NMSSM model included in NMSSMTools_4.4.0 is not absolutely correct. It does not take into account Z^3-non-invariant terms which was included in spectrum calculator of NMSSMTools_4.4.0. In NMSSM_10.14 disposed on CalcHEP web cite this bug is fixed.

      New version of NMSSMTools should appear on WEB this week.

   Best
      Alexander Pukhov

Revision history for this message
Best Alexander Belyaev (alexander.belyaev) said : 		#2

Do you still have problem?

Revision history for this message
tang yilei (tangyilei10) said : 		#3

Thanks Alexander Belyaev, that solved my question.

Revision history for this message
tang yilei (tangyilei10) said : 		#4

Thanks for your answer. I was trying to use calchep to run a process in NMSSM and later link it into pythia, so what I need is the event file. I just modified the model file and assign all parameters with a nonzero value to solve this problem.