constrained parameter rd

rd is constrained variable which contains result of loading of output of
spectrum calculator(SPHeno).
For instance, if SPheno can not calculate particle spectrum for given
parameters and output file is absent, then rd can not be evaluated.
I'll try to check what happens for your input.

Best
    Alexander Pukhov

On 07/29/2016 11:52 PM, ali wrote:
> New question #312370 on CalcHEP:
> https://answers.launchpad.net/calchep/+question/312370
>
> I have problem with CalcHEP. I am studying with MSSM model file. The problem is that, when i entered values which i got from Spheno output into the variables in CalcHEP, i get same error message every time and every i tried process. The error message is "can not evaluate constrained parameter rd". I didn't change anything in constraints. I entered the values into variables as follows
>
> ============================================================================================
> MSSM(EWSB)
> Variables
> Name | Value <|> Comment <|
> alfSMZ |0.1184 |alpha_QCD(MZ)
> MZ |91.3 |Z-boson mass. PDG value 91.1876 +/- 0.0021
> MW |80.2 |W-boson mass. PDG valeu 80.385 +/- 0.015
> Mtp |173.07 |top quark pole mass
> MbMb |4.23 |Mb(Mb) running mass
> McMc |1.27 |Mc(Mc) running mass
> GG |1.117 |Drived by CalcHEP
> Q |100 |QCD scale for running quark masses
> EE |0.3123 |PDG value 4*pi/EE^2=128 corresponds to EE=0.3134
> Ml |1.777 |mass of tau-lepton
> tb |24.35296361 |Tangent beta
> MH3 |143.849602 |Odd Higgs Mass
> mu |-1357.8100 |mu
> MG1 |-421.51075 |
> MG2 |-615.74671 |
> MG3 |123.441571 |
> Ml1 |1352.42389 |ML1
> Ml2 |1352.07558 |ML2
> Ml3 |1254.33106 |ML3
> Mr1 |1302.66244 |MR1
> Mr2 |1301.93640 |MR2
> Mr3 |1087.68383 |MR3
> Mq1 |1694.45904 |
> Mq2 |1693.81237 |
> Mq3 |1134.94101 |
> Mu1 |1644.27234 |
> Mu2 |1644.25408 |
> Mu3 |612.868032 |
> Md1 |1632.41052 |
> Md2 |1632.28849 |
> Md3 |1339.13473 |
> At |-1707.575 |
> Ab |-3303.151 |
> Al |-2461.803 |
> Maux |1 |mass of aux particles
> ============================================================================================
>
>
> Thanks
>

Thank you for your answer. I didnt use spheno with calchep. I used spheno as external for low scale output. Then i gave parameters to calchep. I just changed the Variables in calchep, remaining all the same.

Here SuSpect would not like to work with given set of parameters.
File results/suspect2_lha.out contains

------------------------------------------
BLOCK SPINFO # Spectrum Program information
      1 SuSpect # RGE +Spectrum calculator
      2 2.43 # version number
# Caution: warning or error message follows
      3 Warning: MA^2(Q) <0 at a scale MZ<Q<EWSB !
      4 EWSB maybe unconsistent/not realized
      3 RG-improved V_eff has CCB or UFB problems
---------------------------------------------

Here record "4 EWSB ..." means error.

It is not a good idea to use parameters defined via SPheno
for SuSpect. Although for "standard" scenario a reasonable agreement
is expected, any exotic one can lead to a noticeable difference.
So, I propose to use CalcHEP MSSM model with SPheno. To realize it
you have to change Constraints.

  a) open |openAppend("suspect2_lha.in") ===> open
|openAppend("LesHouches.in")
        name if input file for SPheno

  b) delOld |System("rm -f suspect2_lha.out") ===> delOld |System("rm
-f SPheno.spc")
        instruction to delete old SPheno output.

  c) sys
|System("/home/pukhov/CALCHEP/calchep_3.6.28/MSSM_10.14/SuSpect/suspect.exe")

                                               ===>sys |System("
path to SPheno/SPheno")
                     call SPheno instead on SuSpect

   d) rd |slhaRead("suspect2_lha.out", 0) ===> rd |
slhaRead("SPheno.spc", 0)
             read SPheno output file.

In the current moment I have not SPheno on my computer and can not check
how it works. Indeed with way was tested in micrOMEGAs package and has
to work. The only expected problem is my mistyping in instructions
presented above. Let me know if you still need my help.

Best
     Alexander Pukhov

On 07/30/2016 05:02 PM, ali wrote:
> Question #312370 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/312370
>
> ali posted a new comment:
> Thank you for your answer. I didnt use spheno with calchep. I used
> spheno as external for low scale output. Then i gave parameters to
> calchep. I just changed the Variables in calchep, remaining all the
> same.
>

Thanks you again. I did what you said but i got same error and with some terminal messages, these are;

==================================================

 Warning from routine CalculateMSd
 in the calculation of the masses
 occurred a negative mass squared!
           1 -333581.51773008227
  Loop corrections to Yukawa couplings at EW scale too large!
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP Subroutine TerminateProgram

==================================================

and my constraints file or func4.mdl looks like;

=================================================

..
....
......
open |openAppend("LesHouches.in.MSSM")
input0 |aPrintF("Block MODSEL # Select model\n 1 0 # SUSY scale input\n")
input1 |aPrintF("Block SMINPUTS \n 3 % E # alpha_s(MZ)\n 5 %E # mb(mb)SM MSbar\n 6 %E # Mtop(Pole)\n", alfSMZ, MbMb, Mtp)
input2 |aPrintF("BLOCK MINPAR\n 3 %E # tb\n", tb)
input3 |aPrintF("Block EXTPAR\n 0 -1 # EWSB scale\n 1 %E # MG1\n 2 %E # MG2\n 3 %E # MG3\n 11 %E # At\n 12 %E # Ab\n 13 %E # AL\n", MG1, MG2, MG3, At, Ab, Al)
input4 |aPrintF(" 23 %E # mu \n 26 %E # MH3\n 31 %E # Ml1\n 32 %E #Ml2 \n 33 %E # Ml3\n 34 %E # Mr1\n 35 %E # Mr2\n 36 %E # Mr3\n", mu, MH3, Ml1, Ml2, Ml3, Mr1, Mr2, Mr3)
input5 |aPrintF(" 41 %E # Mq1\n 42 %E # Mq2\n 43 %E # Mq3\n 44 %E # Mu1\n 45 %E #Mu2 \n 46 %E # Mu3\n 47 %E # Md1\n 48 %E # Md2\n 49 %E # Md3\n", Mq1, Mq2, Mq3, Mu1, Mu2, Mu3, Md1, Md2, Md3)
delOld |System("rm -f SPheno.spc.MSSM")
sys |System("/home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM")
rd |slhaRead("SPheno.spc.MSSM", 0)
tB |slhaVal("HMIX", 0, 1, 2)
sb |tB/sqrt(1+tB^2) % Sine beta
cb |sqrt(1-sb^2) % Cosine beta
c2b |cb^2-sb^2 % cos(2b)
MSne |slhaVal("MASS", MZ, 1, 1000012)
MSnm |slhaVal("MASS", MZ, 1, 1000014)
MSeL |slhaVal("MASS", MZ, 1, 1000011)
MSeR |slhaVal("MASS", MZ, 1, 2000011)
MSmL |slhaVal("MASS", MZ, 1, 1000013)
MSmR |slhaVal("MASS", MZ, 1, 2000013)
MSdL |slhaVal("MASS", MZ, 1, 1000001)
...
.....
.......

==========================================

should i change another thing in this file ?

We have input file

        LesHouches.in.MSSM

in directory results. Then we call SPhenoMSSM
from the same directory. I hope SPhenoMSSM reads namely
LesHouches.in.MSSM. In the old version that I tested about two years ago
  SPheno reads LesHouches.in. It is quite reasonable that modern
SPhenoMSSM reads LesHouches.in.MSSM.
So, it should work.

Please, check 2 points. a) Your version CalcHEP+SPheno works fine for
default MSSM point.
                                       b) LesHouches.in.MSSM agrees
with SPheno call you used to find your parameter set.

If you will not find an answer, let me know and I will download last
version of SPheno and check what is a problem.

Best
    Alexander Pukhov

On 07/31/2016 05:02 PM, ali wrote:
> Question #312370 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/312370
>
> ali posted a new comment:
> Thanks you again. I did what you said but i got same error and with some
> terminal messages, these are;
>
> ==================================================
>
> Warning from routine CalculateMSd
> in the calculation of the masses
> occurred a negative mass squared!
> 1 -333581.51773008227
> Loop corrections to Yukawa couplings at EW scale too large!
> Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
> STOP Subroutine TerminateProgram
>
> ==================================================
>
> and my constraints file or func4.mdl looks like;
>
>
> =================================================
>
> ..
> ....
> ......
> open |openAppend("LesHouches.in.MSSM")
> input0 |aPrintF("Block MODSEL # Select model\n 1 0 # SUSY scale input\n")
> input1 |aPrintF("Block SMINPUTS \n 3 % E # alpha_s(MZ)\n 5 %E # mb(mb)SM MSbar\n 6 %E # Mtop(Pole)\n", alfSMZ, MbMb, Mtp)
> input2 |aPrintF("BLOCK MINPAR\n 3 %E # tb\n", tb)
> input3 |aPrintF("Block EXTPAR\n 0 -1 # EWSB scale\n 1 %E # MG1\n 2 %E # MG2\n 3 %E # MG3\n 11 %E # At\n 12 %E # Ab\n 13 %E # AL\n", MG1, MG2, MG3, At, Ab, Al)
> input4 |aPrintF(" 23 %E # mu \n 26 %E # MH3\n 31 %E # Ml1\n 32 %E #Ml2 \n 33 %E # Ml3\n 34 %E # Mr1\n 35 %E # Mr2\n 36 %E # Mr3\n", mu, MH3, Ml1, Ml2, Ml3, Mr1, Mr2, Mr3)
> input5 |aPrintF(" 41 %E # Mq1\n 42 %E # Mq2\n 43 %E # Mq3\n 44 %E # Mu1\n 45 %E #Mu2 \n 46 %E # Mu3\n 47 %E # Md1\n 48 %E # Md2\n 49 %E # Md3\n", Mq1, Mq2, Mq3, Mu1, Mu2, Mu3, Md1, Md2, Md3)
> delOld |System("rm -f SPheno.spc.MSSM")
> sys |System("/home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM")
> rd |slhaRead("SPheno.spc.MSSM", 0)
> tB |slhaVal("HMIX", 0, 1, 2)
> sb |tB/sqrt(1+tB^2) % Sine beta
> cb |sqrt(1-sb^2) % Cosine beta
> c2b |cb^2-sb^2 % cos(2b)
> MSne |slhaVal("MASS", MZ, 1, 1000012)
> MSnm |slhaVal("MASS", MZ, 1, 1000014)
> MSeL |slhaVal("MASS", MZ, 1, 1000011)
> MSeR |slhaVal("MASS", MZ, 1, 2000011)
> MSmL |slhaVal("MASS", MZ, 1, 1000013)
> MSmR |slhaVal("MASS", MZ, 1, 2000013)
> MSdL |slhaVal("MASS", MZ, 1, 1000001)
> ...
> .....
> .......
>
> ==========================================
>
> should i change another thing in this file ?
>

My spheno directory is in : /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8

spheno executable file is : /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM

input file : /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM/LesHouches.in.MSSM

so far the spheno works very well. I can get an output.

my calchep work directory : /home/ali/Desktop/toolbox-2.0.2/calchep_3.6.25/mssm

But then in calchep i changed func4.mdl like this
===========================================
open |openAppend("LesHouches.in.MSSM")
..
..
delOld |System("rm -f SPheno.spc.MSSM")
sys |System("/home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM")
rd |slhaRead("SPheno.spc.MSSM", 0)
...

Zt11 |slhaVal("USQMIX", QSUSY, 2, 6, 3)
...
...
Zb12 |slhaVal("DSQMIX", QSUSY, 2, 6, 4)
...
Zl12 |slhaVal("SELMIX", QSUSY, 2, 6, 4)

here USQMIX, DSQMIX and SELMIX are same in spheno output.

================================================

vars4.mdl file is same as in the first comment.

When i try to use default MSSM values which is in vars4.mdl file with spheno i get same error messages.

an addition may be i can help. ../.../calchep/mssm/results/Messages.out looks like

....
......

 Error in routine LesHouches_Input
 Unknown entry for Block EXTPAR 0
 Error in routine LesHouches_Input
 Unknown entry for Block EXTPAR 11
 Error in routine LesHouches_Input
 Unknown entry for Block EXTPAR 12
 Error in routine LesHouches_Input
 Unknown entry for Block EXTPAR 13
....
......
 Warning from routine CalculateMSd
 in the calculation of the masses
 occurred a negative mass squared!
           1 -6006.3990449243065
 Warning from routine CalculateMSv
 in the calculation of the masses
 occurred a negative mass squared!
           1 -4156.6393339906181
 Warning from routine CalculateMSv
 in the calculation of the masses
 occurred a negative mass squared!
....
.......

Thank you again for your interest

I typed wrong above

input file : /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM/LesHouches.in.MSSM

will be

input file : /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/LesHouches.in.MSSM

they are in same directory executable and input file

I have downloaded SPheno-3.3.8 from https://spheno.hepforge.org/ but
don't see LesHouches.in.MSSM inside.

Best
    Alexander Pukhov

On 07/31/2016 09:52 PM, ali wrote:
> Question #312370 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/312370
>
> ali posted a new comment:
> I typed wrong above
>
> input file :
> /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/SPhenoMSSM/LesHouches.in.MSSM
>
> will be
>
> input file :
> /home/ali/Desktop/toolbox-2.0.2/SPheno-3.3.8/bin/LesHouches.in.MSSM
>
> they are in same directory executable and input file
>

Firstly, I downloaded sarah toolbox "https://www.hepforge.org/downloads/sarah" and i installed it.
After extracted in the sarah directory i typed ./configure then ./make. It will create some directories about some programs like calchep, madgraph, spheno. Then, download this file "https://sarah.hepforge.org/trac/attachment/wiki/MSSM/SPheno.tar.gz". Extract it in the "../toolbox-2.0.2/SPheno-3.3.8/", after this step, i changed the folder name which is SPheno to MSSM. While you are in the ../SPheno-3.3.8/ folder you should type on the terminal "make Model=MSSM" then spheno will create "SPhenoMSSM" executable file into the "/bin" folder and you can find the "LesHouches.in.MSSM" in MSSM folder(which was name SPheno at the beginning). After that, you can copy "LesHouches.in.MSSM" into the "/bin" folder. Finaly when you type "./SPhenoMSSM LesHouches.in.MSSM" on the terminal you can get an output like .spc.MSSM.

OK, I know about Sarah but never tested this way.
I'll look today.
Best
     Alexander Pukhov

On 08/01/2016 01:52 AM, ali wrote:
> Question #312370 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/312370
>
> ali posted a new comment:
> Firstly, I downloaded sarah toolbox "https://www.hepforge.org/downloads/sarah" and i installed it.
> After extracted in the sarah directory i typed ./configure then ./make. It will create some directories about some programs like calchep, madgraph, spheno. Then, download this file "https://sarah.hepforge.org/trac/attachment/wiki/MSSM/SPheno.tar.gz". Extract it in the "../toolbox-2.0.2/SPheno-3.3.8/", after this step, i changed the folder name which is SPheno to MSSM. While you are in the ../SPheno-3.3.8/ folder you should type on the terminal "make Model=MSSM" then spheno will create "SPhenoMSSM" executable file into the "/bin" folder and you can find the "LesHouches.in.MSSM" in MSSM folder(which was name SPheno at the beginning). After that, you can copy "LesHouches.in.MSSM" into the "/bin" folder. Finaly when you type "./SPhenoMSSM LesHouches.in.MSSM" on the terminal you can get an output like .spc.MSSM.
>

Hi again, did you have chance to be able to test spheno with calchep.

Sorry,
  I had not time to look at your problem. My colleagues asked me about
some urgent job.
  I remember your request and will try to find a time to answer you.
Best
    Alexander Pukhov

On 08/03/2016 11:08 AM, ali wrote:
> Question #312370 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/312370
>
> ali posted a new comment:
> Hi again, did you have chance to be able to test spheno with calchep.
>

ok no problem, thanks a lot.