Now lets get HDECAY running to generate a MSSM Higgs spectrum to test with CheckMATE.
> mkdir HDecay
> cd HDecay/
> wget http://tiger.web.psi.ch/hdecay/hdecay.tar.gz
> tar -xvzf hdecay.tar.gz
> make
The executable is calles "run", and it reads the input file "hdecay.in". Take some time to familiarize yourself with this file. Now modify hdecay.in to produce an SLHA (SUSY Les Houches Accord) file to feed to Pythia8. Change the line "SLHAOUT = 0" to "SLHAOUT = 1", and run:
> ./run
The SLHA file is called "slha.out".
> cp slha.out /home/potter/pythia8235/lhc/
> cd /home/potter/pythia8235/lhc/
Modify mainLHC.cc to produce all BSM Higgs processes at 13 TeV, using your slha.out spectrum:
...
HepMC::IO_GenEvent ascii_io("bsm_higgs_lhc.hepmc", std::ios::out);
...
pythia.readString("Beams:eCM = 13000.");
pythia.readString("Higgs:useBSM=on");
pythia.readString("HiggsBSM:all=on");
pythia.readString("SLHA:readFrom=2");
pythia.readString("SLHA:file=slha.out");
...
for (int iEvent = 0; iEvent < 10000; ++iEvent) {s
Excellent! Now compile and run, then copy the HepMC file to CheckMATE so we can check its viability. Note the total cross section and error so we can tell this to CheckMATE. Pythia8 gives it in millibarns, CheckMATE want femtobarns.
> make mainLHC
> ./mainLHC | tee bsm.log
> cp bsm_higgs_lhc.hepmc /home/potter/CheckMATE-2.0.26
Now let's try it on CheckMATE. Run "./bin/CheckMATE" without arguments to get help, or "./bin/CheckMATE -h" for more help. Here we give this run a name (bsm_higgs), tell CheckMATE which analyses to use (all implemented ATLAS searches for 13 TeV), give the process a name (bsm_hall), give it the cross section and error reported by Pythia8, and tell it to use our Pythia8 output file:
> cd /home/potter/CheckMATE-2.0.26
> ./bin/CheckMATE -n bsm_higgs -a "atlas_1604_01306,atlas_1605_09318,atlas_1609_01599,atlas_1704_03848,atlas_conf_2015_082,atlas_conf_2016_013,atlas_conf_2016_050,atlas_conf_2016_054,atlas_conf_2016_076,atlas_conf_2016_096,atlas_conf_2017_060,atlas_conf_2016_066,atlas_1712_08119,atlas_1712_02332,atlas_1709_04183,atlas_1802_03158,atlas_1708_07875,atlas_1706_03731" -p "bsm_hall" -xs "231400*FB" -xse "48170*FB" -ev "bsm_higgs_lhc.hepmc"