Higgsino Production

The example SUSY model

• We used softsusy-3.3.8 to generate a spectrum.
• We then used susyhit-1.3 to generate the decay table from the spectrum
• We then struggled a lot, until Alexis told us to use the madgraph web interface to translate the output from susyhit-1.3 into something that can be used as param_card.dat file in madgraph.

Here are the various files necessary to reproduce this chain:

softsusy input file

Block MODSEL                 # Select model
 1    0                        # Low energy MSSM
Block SMINPUTS               # Standard Model inputs
 1   1.27908953E+02       # alpha_em^(-1)(MZ) SM MSbar
 2   1.16637000E-05       # G_Fermi 
 3   1.18400000E-01       # alpha_s(MZ) SM MSbar
 5   4.19000000E+00       # mb(mb) SM MSbar
 6   1.73300000E+02       # mtop(pole)
 7   1.77700000E+00       #  Mtau     
Block MINPAR                 # Input parameters
 3   1.00000000E+01       # tanb
Block EXTPAR
 0   -1        # EWSB
 1   5.50000000E+02      # MG1
 2   5.13350000E+02      # MG2
 3   2.00000000E+03      # MG3
 11  2.70000000E+03      # At 
 12  0.00000000E+00      # Ab 
 13  0.00000000E+00      # Atau
 23  2.00000000E+02      # mu
 26  1.00000000E+03      # MA
 31  3.00000000E+03      # Ml1
 32  3.00000000E+03      # Ml2
 33  3.00000000E+03      # Ml3
 34  3.00000000E+03      # MR2
 35  3.00000000E+03      # MR2
 36  3.00000000E+03      # MR3
 41  3.00000000E+03      # Mq1
 42  3.00000000E+03      # Mq2
 43  9.00000000E+02      # Mq3
 44  3.00000000E+03      # Mu1
 45  3.00000000E+03      # Mu2
 46  1.35000000E+03      # Mu3
 47  3.00000000E+03      # Md1
 48  3.00000000E+03      # Md2
 49  3.00000000E+03      # Md3
Block SOFTSUSY               # SOFTSUSY specific inputs
         2   0                    # quark mixing option

softsusy output file which is also the susyhit input file

• naturalSusyHiggsinoLSP.slha: softsusy output file

susyhit output file

• naturalSusyHiggsinoLSP-out.slha: susyhit output file

param_card.dat file after madgraph conversion

• param_card_1.dat: madgraph param_card after conversion from susyhit output

Cross sections from prospino for this model

nn  1  1     0.0    0.0    1.0  190.5  190.5  0.000 0.118E-03 0.821E-03 0.149E-03 0.693E-03 1.2582 0.118E-03 0.149E-03
nn  1  2     0.0    0.0    1.0  190.5 -207.5  0.000 0.205     0.758E-03 0.263     0.605E-03 1.2872 0.205     0.264    
nn  2  2     0.0    0.0    1.0 -207.5 -207.5  0.000 0.169E-04 0.817E-03 0.213E-04 0.666E-03 1.2556 0.169E-04 0.213E-04
nn  1  5     0.0    0.0    1.0  190.5  196.5  0.000 0.297     0.754E-03 0.379     0.604E-03 1.2759 0.298     0.380    
nn  1  7     0.0    0.0    1.0  190.5  196.5  0.000 0.161     0.788E-03 0.208     0.628E-03 1.2911 0.162     0.209    
nn  2  5     0.0    0.0    1.0 -207.5  196.5  0.000 0.239     0.754E-03 0.304     0.604E-03 1.2737 0.239     0.305    
nn  2  7     0.0    0.0    1.0 -207.5  196.5  0.000 0.128     0.789E-03 0.165     0.629E-03 1.2891 0.128     0.166    
nn  5  7     0.0    0.0    1.0  196.5  196.5  0.000 0.244     0.758E-03 0.311     0.606E-03 1.2758 0.245     0.312    

    i1 i2  dummy0 dummy1 scafac  m1    m2      angle LO[pb]   rel_error NLO[pb]   rel_error   K    LO_ms[pb] NLO_ms[pb]  

Here’s how to read this:
1 = X10
2 = X20
5 = X1+
7 = X1-

then you go to the last column and read off the cross section in pb.
E.g. X20X10 = 0.264 pb, X10 X1+ = 0.380 pb etc.
As you can see, two of the cross sections (X10X10 and X20X20) are negligible, and you can ignore these when you generate
events. The rest are all relevant.

The total Xsection summed over all production channels is ~1.6 pb.

Also, the second to last column is a less accurate calculation. The cross section mad graph calculates should be somewhere between
the last and second to last column. If it isn’t then something is wrong.

Note that I did not bother to calculate in the original calculation the Xsections for production involving n3,n4, or x2.

-- FkW - 2015/09/09