P3D Contour Plot by varying variables from different blocks

[akdatta]

Hi,

My goal is to generate a P3D contour plot of some observable by varying two variables that belong to two different BLOCKS.

For such a plot, when the two variables belong to the same BLOCK, say HMIXIN, the necessary variations are done
using a script like

-------------------------------------------------------------------------------------------------
DEFINITION[ScanMuMup][HMIXIN]={
{{101}, {Value->4*10^4}},
{{1}, {Min->-300,Max->300, Steps->7,Distribution->LINEAR}},
{{20}, {Min->-1000,Max->1000, Steps->9,Distribution->LINEAR}}
};
--------------------------------------------------------------------------------------------------

What would be the corresponding script if two variables from two different blocks are to be varied for the purpose?

Could you please comment.
AD

[FStaub]

Just write the ranges in the two blocks.

[akdatta]

Hi,
Thanks a lot for your reply.
I am afraid, I perhaps could not get your suggestion correctly.

A concrete example is as follows.

I am trying to scan over two input parameters:
M1 from the MSOFTIN block and
HMIXIN[20] from the block HMIXIN

in order to finally draw Higgs mass contours in the plane of
M1 and HMIX[20].

The following is setup is leading to the message

total number of parameter points: 0

finished!

Would you please help.
AD

--------------------------------------------------------------------------------------------
DEFINITION[ScanTest][MSOFTIN]={
{{1},{Min->150,Max->550,Steps->5,Distribution->Linear}} (*M1*) ,
{{2},{Value->400}} (*M2*) ,
{{3},{Value->2000}} (*M3*)
};

DEFINITION[ScanTest][HMIXIN]={
{{101}, {Value->4*10^4}},
{{1}, {Value->250}},
{{20}, {Min->-1000,Max->1000, Steps->9,Distribution->LINEAR}}
};

DEFINITION[ScanTest][Plots]={
{P3D, {MSOFT[1],HMIX[20],MASS[25]},Style5,"M1_mup_Mass25.eps"}
};
--------------------------------------------------------------------------------------------

[FStaub]

It should work like this. Did you run SPheno without SSP to see if the parameter combination produces a valid spectrum?

[akdatta]

Thanks a lot, Florian, for your comment.

I haven't run SPheno separately for the scan-ranges I mentioned in my last post.
However, using SSP, if I scan one of the parameters from one block (say, HMIXIN[20])
over the range I mentioned and keep M1, i.e. MSOFTIN[1] fixed at any particular
value in the mentioned range, I get an apparently valid output spectrum in the file
"SpectrumFiles.spc". I am attaching this spectrum file and my Input file I used in the
SSP scan for your perusal. The notebook message for such a one-variable scan is as
follows (I do not expect a valid plot from this 1D scan though, given the Input file I used).
-------------------------------------------------------------------------------------------------------------------------
Start["gridscan-nhssm.m"]

Starting parameter scan ScanTest
total number of parameter points: 9
Reading all spectrum files
Creating plots

finished!
-------------------------------------------------------------------------------------------------------------------------

It seems that in my reported attempt to scan over two parameters each from a separate
block, the run abruptly prompted

"total number of parameter points: 0"

To begin with, the run perhaps is not registering
5 MSOFTIN[1] points (150 GeV --- 550 GeV in 5 steps)
times 9 HMIXIN[20] points (-1000 GeV to 1000 GeV in 9 steps)

i.e., 45 parameter points.

Could you please comment.

Regards.
AD

[akdatta]

Sorry, Florian.
The output I posted is for M1=150 GeV.

However, the input file shows M1=550 GeV.
This is just because I changed the value of M1
from 150 to 550 GeV in this file for a later run.

AD

[FStaub]

Hi,

you write "Linear" but it must be "LINEAR"

Cheers,
Florian

[akdatta]

Extremely sorry, Florian for the sloppiness on my part.

It's working! Thank you so much for your time.

Cheers.
Asesh