namespace Gambit::ColliderBit

[No description available]

Classes

Name
structGambit::ColliderBit::_gsl_target_func_params
A struct to contain parameters for the GSL optimiser target function.
classGambit::ColliderBit::Analysis_ATLAS_13TeV_0LEP_139invfb
ATLAS Run 2 0-lepton jet+MET SUSY analysis, with 139/fb of data.
classGambit::ColliderBit::Analysis_ATLAS_13TeV_0LEP_13invfb
ATLAS Run 2 0-lepton jet+MET SUSY analysis, with 13/fb of data.
classGambit::ColliderBit::Analysis_ATLAS_13TeV_0LEP_36invfb
ATLAS Run 2 0-lepton jet+MET SUSY analysis, with 36/fb of data.
classGambit::ColliderBit::Analysis_ATLAS_13TeV_0LEPStop_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_1LEPStop_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2bMET_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2LEPStop_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2LEPStop_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2LEPStop_exclusive_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2LEPStop_inclusive_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_chargino_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_chargino_80invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_chargino_binned_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_chargino_binned_80invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_chargino_inclusive_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_chargino_inclusive_80invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_2OSLEP_Z_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_3b_24invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_3b_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_3b_discoverySR_24invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_3b_discoverySR_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_4LEP_139invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_4LEP_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_2Lep0Jets_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_2LepPlusJets_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_3Lep_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_strong_139invfb
ATLAS Run 2 search for same-sign leptons and jets, with 139/fb of data.
classGambit::ColliderBit::Analysis_ATLAS_13TeV_PhotonGGM_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_RJ3L_2Lep2Jets_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_RJ3L_3Lep_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_RJ3L_lowmass_36invfb
classGambit::ColliderBit::Analysis_ATLAS_13TeV_ZGammaGrav_CONFNOTE_80invfb
ATLAS ZH(->photon+gravitino) (79.8 fb^-1)
classGambit::ColliderBit::Analysis_ATLAS_7TeV_1OR2LEPStop_4_7invfb
classGambit::ColliderBit::Analysis_ATLAS_7TeV_2LEPStop_4_7invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_0LEP_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_0LEPStop_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_1LEPbb_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_1LEPStop_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_2bStop_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_2LEPEW_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_2LEPStop_20invfb
classGambit::ColliderBit::Analysis_ATLAS_8TeV_3LEPEW_20invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_0LEP_137invfb
CMS Run 2 0-lepton jet+MET SUSY analysis, with 137/fb of data.
classGambit::ColliderBit::Analysis_CMS_13TeV_0LEP_13invfb
CMS Run 2 0-lepton jet+MET SUSY analysis, with 13/fb of data.
classGambit::ColliderBit::Analysis_CMS_13TeV_0LEP_36invfb
CMS Run 2 0-lepton jet+MET SUSY analysis, with 36/fb of data.
classGambit::ColliderBit::Analysis_CMS_13TeV_1LEPbb_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_1LEPStop_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_1Photon1Lepton_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_1Photon1Lepton_emu_combined_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2LEPsoft_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2LEPsoft_36invfb_nocovar
classGambit::ColliderBit::Analysis_CMS_13TeV_2LEPsoft_stop_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2LEPsoft_stop_36invfb_nocovar
classGambit::ColliderBit::Analysis_CMS_13TeV_2LEPStop_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2OSLEP_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2OSLEP_36invfb_nocovar
classGambit::ColliderBit::Analysis_CMS_13TeV_2OSLEP_chargino_stop_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2OSLEP_confnote_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2OSLEP_for_chargino_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2OSLEP_for_stop_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2Photon_GMSB_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2SSLEP_Stop_137invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2SSLEP_Stop_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2SSLEP_Stop_exclusive_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_2SSLEP_Stop_inclusive_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MONOJET_36invfb
CMS Run 2 monojet analysis (no W/Z region) with 36/fb of data.
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_2SSLep_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_3Lep_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_Full_2SSLep_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_Full_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_Full_3Lep_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_MultiLEP_Full_3Lep_rebinned_36invfb
classGambit::ColliderBit::Analysis_CMS_13TeV_Photon_GMSB_36invfb
classGambit::ColliderBit::Analysis_CMS_8TeV_1LEPDMTOP_20invfb
classGambit::ColliderBit::Analysis_CMS_8TeV_2LEPDMTOP_20invfb
classGambit::ColliderBit::Analysis_CMS_8TeV_MONOJET_20invfb
classGambit::ColliderBit::Analysis_CMS_8TeV_MultiLEP_20invfb
classGambit::ColliderBit::Analysis_CMS_8TeV_MultiLEP_3Lep_20invfb
classGambit::ColliderBit::Analysis_CMS_8TeV_MultiLEP_4Lep_20invfb
classGambit::ColliderBit::Analysis_Covariance
Dummy analysis code with a hard-coded return including a SR covariance matrix.
classGambit::ColliderBit::Analysis_Minimum
Basic analysis code for copying.
classGambit::ColliderBit::DMEFT_analysis_info
classGambit::ColliderBit::MT2
A useful MT2 class for this module.

Functions

Name
voidLEP208_SLHA1_convention_xsec_selselbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_selserbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_serserbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_serselbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_se1se1bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_se1se2bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_se2se2bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_se2se1bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smulsmulbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smulsmurbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smursmurbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smursmulbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smu1smu1bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smu1smu2bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smu2smu2bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_smu2smu1bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_staulstaulbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_staulstaurbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_staurstaurbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_staurstaulbar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_stau1stau1bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_stau1stau2bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_stau2stau2bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_stau2stau1bar(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_11(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_12(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_13(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_14(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_22(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_23(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_24(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_33(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_34(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chi00_44(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chipm_11(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chipm_12(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chipm_22(triplet< double > & result)
voidLEP208_SLHA1_convention_xsec_chipm_21(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_selselbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_selserbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_serserbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_serselbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_se1se1bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_se1se2bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_se2se2bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_se2se1bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smulsmulbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smulsmurbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smursmurbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smursmulbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smu1smu1bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smu1smu2bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smu2smu2bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_smu2smu1bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_staulstaulbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_staulstaurbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_staurstaurbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_staurstaulbar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_stau1stau1bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_stau1stau2bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_stau2stau2bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_stau2stau1bar(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_11(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_12(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_13(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_14(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_22(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_23(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_24(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_33(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_34(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chi00_44(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chipm_11(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chipm_12(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chipm_22(triplet< double > & result)
voidLEP205_SLHA1_convention_xsec_chipm_21(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_selselbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_selserbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_serserbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_serselbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_se1se1bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_se1se2bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_se2se2bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_se2se1bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smulsmulbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smulsmurbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smursmurbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smursmulbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smu1smu1bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smu1smu2bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smu2smu2bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_smu2smu1bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_staulstaulbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_staulstaurbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_staurstaurbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_staurstaulbar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_stau1stau1bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_stau1stau2bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_stau2stau2bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_stau2stau1bar(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_11(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_12(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_13(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_14(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_22(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_23(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_24(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_33(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_34(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chi00_44(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chipm_11(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chipm_12(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chipm_22(triplet< double > & result)
voidLEP188_SLHA1_convention_xsec_chipm_21(triplet< double > & result)
voidALEPH_Selectron_Conservative_LLike(double & result)
voidALEPH_Smuon_Conservative_LLike(double & result)
voidALEPH_Stau_Conservative_LLike(double & result)
voidL3_Selectron_Conservative_LLike(double & result)
voidL3_Smuon_Conservative_LLike(double & result)
voidL3_Stau_Conservative_LLike(double & result)
voidL3_Neutralino_All_Channels_Conservative_LLike(double & result)
voidL3_Neutralino_Leptonic_Conservative_LLike(double & result)
voidL3_Chargino_All_Channels_Conservative_LLike(double & result)
voidL3_Chargino_Leptonic_Conservative_LLike(double & result)
voidOPAL_Chargino_Hadronic_Conservative_LLike(double & result)
voidOPAL_Chargino_SemiLeptonic_Conservative_LLike(double & result)
voidOPAL_Chargino_Leptonic_Conservative_LLike(double & result)
voidOPAL_Degenerate_Chargino_Conservative_LLike(double & result)
voidOPAL_Chargino_All_Channels_Conservative_LLike(double & result)
voidOPAL_Neutralino_Hadronic_Conservative_LLike(double & result)
doubleI1(double s, double m1, double m2, double mk, double ml)
doubleI2(double s, double m1, double m2, double mk, double ml)
doubleI3(double s, double m1, double m2, double mk)
voidSLHA2BFM_NN(MixMatrix & NN, double tanb, double sin2thetaW)
Functions to convert mass matrices between SLHA and BFM conventions ////////////////////////////////////////////////////////////////////.
voidSLHA2BFM_VV(MixMatrix & VV)
Converts the chargino mixing matrix V in SLHA conventions to BFM conventions.
voidBFM2SLHA_NN(MixMatrix & NN, double tanb, double sin2thetaW)
Converts a neutralino mixing matrix in BFM conventions to SLHA conventions, (\tan\beta) is as defined in SLHA.
voidBFM2SLHA_VV(MixMatrix & VV)
Converts the chargino mixing matrix V in BFM conventions to SLHA conventions.
MixMatrixmultiply(MixMatrix A, MixMatrix B)
Helper function to multiply matrices.
MixMatrixtranspose(MixMatrix A)
Helper function to find matrix transpose.
voidprint(MixMatrix A)
Helper function to print a matrix.
boolsortByPT13(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByPT13_sharedptr(std::shared_ptr< HEPUtils::Jet > jet1, std::shared_ptr< HEPUtils::Jet > jet2)
boolsortByMass(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByMass_sharedptr(std::shared_ptr< HEPUtils::Jet > jet1, std::shared_ptr< HEPUtils::Jet > jet2)
doublecalcMT(HEPUtils::P4 jetMom, HEPUtils::P4 metMom)
boolsortByPT_1l(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByPT_1l_sharedptr(std::shared_ptr< HEPUtils::Jet > jet1, std::shared_ptr< HEPUtils::Jet > jet2)
boolsortByMass_1l(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByMass_1l_sharedptr(std::shared_ptr< HEPUtils::Jet > jet1, std::shared_ptr< HEPUtils::Jet > jet2)
doublecalcMT_1l(HEPUtils::P4 jetMom, HEPUtils::P4 metMom)
boolsortByPT_j(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByPT_l(const HEPUtils::Particle * lep1, const HEPUtils::Particle * lep2)
boolsortByPT_jet(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByPT_lep(const HEPUtils::Particle * lep1, const HEPUtils::Particle * lep2)
boolsortByPT_RJ3L(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortLepByPT_RJ3L(const HEPUtils::Particle * lep1, const HEPUtils::Particle * lep2)
boolSortLeptons(const pair< TLorentzVector, int > lv1, const pair< TLorentzVector, int > lv2)
boolSortJets(const TLorentzVector jv1, const TLorentzVector jv2)
DEFINE_ANALYSIS_FACTORY(ATLAS_13TeV_ZGammaGrav_CONFNOTE_80invfb )
boolsortByPT(const HEPUtils::Jet * jet1, const HEPUtils::Jet * jet2)
boolsortByPT_2lep(const HEPUtils::Particle * lep1, const HEPUtils::Particle * lep2)
doublePhi_mpi_pi(double x)
double_Phi_mpi_pi(double x)
Analysis *mkAnalysis(const str & name)
Factory definition.
strgetDetector(const str & name)
Return the detector to be used for a given analysis name (and check that the analysis exists).
void_setEventWeight_unity(HEPUtils::Event & event, const BaseCollider * )
A function that sets the event weight to unity, with zero uncertainty.
voidsetEventWeight_unity(EventWeighterFunctionType & result)
void_setEventWeight_fromCrossSection(HEPUtils::Event & event, const BaseCollider * HardScatteringSim_ptr, const map_int_process_xsec & ProcessCrossSectionsMap, const int use_trust_level)
A function that sets the event weight based on the process cross-sections.
voidsetEventWeight_fromCrossSection(EventWeighterFunctionType & result)
voidgetDummyColliderObservable(double & result)
Dummy observable that creates a dependency on TestModel1D.
voidoperateLHCLoop(MCLoopInfo & result)
LHC Loop Manager.
voidgetLHCEventLoopInfo(map_str_dbl & result)
Store some information about the event generation.
voidCollectAnalyses(AnalysisDataPointers & result)
Loop over all analyses and collect them in one place.
voidset_CS(hb_ModelParameters & result, const HiggsCouplingsTable & couplings, int n_neutral_higgses)
Helper function to set HiggsBounds/Signals parameters cross-section ratios from a GAMBIT HiggsCouplingsTable.
voidset_SMLikeHiggs_ModelParameters(const SubSpectrum & spec, const HiggsCouplingsTable & couplings, hb_ModelParameters & result)
Helper function for populating a HiggsBounds/Signals ModelParameters object for SM-like Higgs.
voidSMLikeHiggs_ModelParameters(hb_ModelParameters & result)
SM-like (SM + possible invisibles) Higgs model parameters for HiggsBounds/Signals.
voidMSSMLikeHiggs_ModelParameters(hb_ModelParameters & result)
MSSM-like (MSSM + NMSSM + …) Higgs model parameters for HiggsBounds/Signals.
voidcalc_HB_LEP_LogLike(double & result)
Get a LEP chisq from HiggsBounds.
voidcalc_HS_LHC_LogLike(double & result)
Get an LHC chisq from HiggsSignals.
voidFeynHiggs_HiggsProd(fh_HiggsProd_container & result)
Higgs production cross-sections from FeynHiggs.
AnalysisLogLikescalc_loglikes_for_analysis(const AnalysisData & adata, bool USE_COVAR, bool USE_MARG, bool combine_nocovar_SRs, bool set_zero_loglike)
Forward declaration of funtion in LHC_likelihoods.
voidfill_analysis_info_map()
Forward declarations of functions in this file.
voidDMEFT_results(AnalysisDataPointers & result)
Results from DMEFT analyses before any modification of the MET spectrum.
voidget_all_DMEFT_signal_yields(std::vector< double > & sr_nums, const DMEFT_analysis_info & analysis_info, const Spectrum & spec)
Fill the input vector with the total DMEFT signal prediction for each SR in the given LHC analysis.
voidget_DMEFT_signal_yields_dim6_operator(std::vector< double > & signal_yields, const str, const DMEFT_analysis_info & analysis_info, double m, double O1, double O2, double lambda)
Fill the input vector with the DMEFT signal prediction for a given set of dim-6 operators.
voidget_DMEFT_signal_yields_dim7_operator(std::vector< double > & signal_yields, const str, const DMEFT_analysis_info & analysis_info, double m, double O, double lambda)
Fill the input vector with the DMEFT signal prediction for a given dim-7 operator.
voidDMEFT_results_profiled(AnalysisDataPointers & result)
Results from DMEFT analyses after profiling over the ‘a’ parameter in the smooth cut-off of the MET spectrum.
voidDMEFT_results_cutoff(AnalysisDataPointers & result)
Results from DMEFT analyses after imposing a hard cut-off of the MET spectrum.
voidsignal_modifier_function(AnalysisData & adata, double lambda, double a)
voidsignal_cutoff_function(AnalysisData & adata, double lambda)
void_gsl_target_func(const size_t, const double * a, void * fparams, double * fval)
A target function for the GSL optimiser.
voidcalc_DMEFT_profiled_LHC_nuisance_params(map_str_dbl & result)
voidInterpolatedMCInfo(MCLoopInfo & result)
doublelimit_LLike(double x, double x95, double sigma)
LEP limit likelihood function.
boolis_xsec_sane(const triplet< double > & xsecWithError)
LEP limit debugging function.
voidLEP207_SLHA1_convention_xsec_chi00_11(triplet< double > & result)
voidL3_Gravitino_LLike(double & result)
const multimap_PID_pair_int &all_PID_pairs_to_process_codes()
voidgetActiveProcessCodes(std::vector< int > & result)
Get the list of active collider process codes.
voidgetAnalysisContainer(AnalysisContainer & result, const str & detname, const MCLoopInfo & RunMC, const xsec_container & TotalCrossSection, int iteration)
Retrieve an analysis container for a specific detector.
voidgetBuckFastATLAS(BaseDetector *& result)
Retrieve a BuckFast sim of ATLAS.
voidgetBuckFastCMS(BaseDetector *& result)
Retrieve a BuckFast sim of CMS.
voidgetBuckFastIdentity(BaseDetector *& result)
Retrieve an Identity BuckFast sim (no sim)
voidgetHepMCEvent(HEPUtils::Event & result)
A nested function that reads in HepMC event files and converts them to HEPUtils::Event format.
voidgetLHEvent(HEPUtils::Event & result)
A nested function that reads in Les Houches Event files and converts them to HEPUtils::Event format.
std::pair< double, double >convert_xsecs_to_fb(double input_xsec, double input_xsec_uncert, str input_unit, bool input_fractional_uncert)
voidgetPIDPairCrossSectionsMap_prospino(map_PID_pair_PID_pair_xsec & result)
PID_pair_xsec_containersilly_pid_xsec_constructor(PID_pair pid_pair, double xsec_val)
Test functions for provding PIDPairCrossSectionsMap (cross-sections in fb)
voidgetPIDPairCrossSectionsMap_testing(map_PID_pair_PID_pair_xsec & result)
voidgetProcessCrossSectionsMap(map_int_process_xsec & result)
Get a map between Pythia process codes and cross-sections.
voidgetEvGenCrossSection(MC_xsec_container & result)
Compute a cross-section from Monte Carlo.
voidgetEvGenCrossSection_as_base(xsec_container & result)
Return MC_xsec_container as the base xsec_container.
voidgetNLLFastCrossSection(xsec_container & result)
Get a cross-section from NLL-FAST.
boolcheckOptions_getYAMLCrossSection(const Options & runOptions, const str calling_function, std::pair< str, str > & xsec_pnames, str & input_unit, bool & input_fractional_uncert, str & errmsg)
A helper function to check the YAML options for getYAMLCrossSection and getYAMLCrossSection_SLHA.
voidgetYAMLCrossSection(xsec_container & result)
A function that reads the total cross-section from the input file, but builds up the number of events from the event loop.
voidgetYAMLCrossSection_SLHA(xsec_container & result)
A function that reads a list of (SLHA file, total cross-section) pairs from the input YAML file.
voidgetYAMLCrossSection_param(xsec_container & result)
voidgetTotalCrossSectionAsMap(map_str_dbl & result)
Get cross-section info as map_str_dbl (for simple printing)
voidgetPIDPairCrossSectionsInfo(map_str_dbl & result)
Output PID pair cross-sections as a str-dbl map, for easy printing.
voiddoCrossSectionConsistencyCheck(bool & result)
voidget_sigma_ee_ll(triplet< double > & result, const double sqrts, const int generation, const int l_chirality, const int lbar_chirality, const double gtol, const double ftol, const bool gpt_error, const bool fpt_error, const Spectrum & spec, const double gammaZ, const bool l_are_gauge_es)
voidget_sigma_ee_chi00(triplet< double > & result, const double sqrts, const int chi_first, const int chi_second, const double tol, const bool pt_error, const Spectrum & spec, const double gammaZ)
Retrieve the production cross-section at an e+e- collider for neutralino pairs.
voidget_sigma_ee_chipm(triplet< double > & result, const double sqrts, const int chi_plus, const int chi_minus, const double tol, const bool pt_error, const Spectrum & spec, const double gammaZ)
Retrieve the production cross-section at an e+e- collider for chargino pairs.
doublexsec_sleislej(int pid1, int pid2, double sqrts, double m1, double m2, MixMatrix F, MixMatrix N, const double mN[4], double alpha, double mZ, double gZ, double sin2thetaW, bool CP_lock)
doublexsec_neuineuj(int pid1, int pid2, double sqrts, double mi, double mj, MixMatrix N, const double mS[2], double tanb, double alpha, double mZ, double gZ, double sin2thetaW)
doublexsec_chaichaj(int pid1, int pid2, double sqrts, double mi, double mj, MixMatrix V, MixMatrix U, double ms, double alpha, double mZ, double gZ, double sin2thetaW)
voidcalc_LHC_signals(map_str_dbl & result)
Loop over all analyses and fill a map of predicted counts.
void_gsl_calc_Analysis_MinusLogLike(const size_t n, const double * unit_nuisances_dbl, void * fixedparamspack, double * fval)
void_gsl_calc_Analysis_MinusLogLikeGrad(const size_t n, const double * unit_nuisances_dbl, void * fixedparamspack, double * fgrad)
Loglike gradient-function wrapper to provide the signature for GSL multimin.
void_gsl_calc_Analysis_MinusLogLikeAndGrad(const size_t n, const double * unit_nuisances_dbl, void * fixedparamspack, double * fval, double * fgrad)
std::vector< double >_gsl_mkpackedarray(const Eigen::ArrayXd & n_preds, const Eigen::ArrayXd & n_obss, const Eigen::ArrayXd & sqrtevals, const Eigen::MatrixXd & evecs)
doubleprofile_loglike_cov(const Eigen::ArrayXd & n_preds, const Eigen::ArrayXd & n_obss, const Eigen::ArrayXd & sqrtevals, const Eigen::MatrixXd & evecs)
doublemarg_loglike_nulike1sr(const Eigen::ArrayXd & n_preds, const Eigen::ArrayXd & n_obss, const Eigen::ArrayXd & sqrtevals)
doublemarg_loglike_cov(const Eigen::ArrayXd & n_preds, const Eigen::ArrayXd & n_obss, const Eigen::ArrayXd & sqrtevals, const Eigen::MatrixXd & evecs)
voidcalc_LHC_LogLikes(map_str_AnalysisLogLikes & result)
Loop over all analyses and fill a map of AnalysisLogLikes objects.
voidget_LHC_LogLike_per_analysis(map_str_dbl & result)
Extract the combined log likelihood for each analysis.
voidget_LHC_LogLike_per_SR(map_str_dbl & result)
Extract the log likelihood for each SR.
voidget_LHC_LogLike_SR_labels(map_str_str & result)
Extract the labels for the SRs used in the analysis loglikes.
voidget_LHC_LogLike_SR_indices(map_str_dbl & result)
voidcalc_combined_LHC_LogLike(double & result)
Compute the total likelihood combining all analyses.
voidcalc_LHC_LogLike_scan_guide(double & result)
voidgetNextSLHAFileNameAndContent(pair_str_SLHAstruct & result)
voidgetAndReplaceSLHAContent(pair_str_SLHAstruct & result)
voidgetSLHAFileElements(map_str_dbl & result)
voidgetSLHA1Spectrum(SLHAstruct & result)
voidgetSLHA2Spectrum(SLHAstruct & result)
voidcalc_susy_spectrum_scan_guide(double & result)
voidgetActiveProcessCodeToPIDPairsMap(multimap_int_PID_pair & result)
voidgetActivePIDPairs(vec_PID_pair & result)
Get a list of all the PID pairs related to active process codes.
voidrunAnalyses(AnalysisDataPointers & result, const str & , const MCLoopInfo & RunMC, const AnalysisContainer & Container, const HEPUtils::Event & SmearedEvent, int iteration, void(*)() wrapup)
Run all the analyses in a given container.
voidsmearEvent(HEPUtils::Event & result, const HEPUtils::Event & HardScatteringEvent, const BaseDetector & detector, const MCLoopInfo & RunMC, const int iteration, const str & detname)
Smear an event.
boolrandom_bool(double eff)
voidfiltereff(std::vector< const HEPUtils::Particle * > & particles, double eff, bool do_delete)
voidfiltereff(std::vector< const HEPUtils::Particle * > & particles, std::function< double(const HEPUtils::Particle *)> eff_fn, bool do_delete)
Utility function for filtering a supplied particle vector by sampling wrt a binned 1D efficiency map in pT.
voidfiltereff_pt(std::vector< const HEPUtils::Particle * > & particles, const HEPUtils::BinnedFn1D< double > & eff_pt, bool do_delete)
voidfiltereff_etapt(std::vector< const HEPUtils::Particle * > & particles, const HEPUtils::BinnedFn2D< double > & eff_etapt, bool do_delete)
std::vector< std::vector< const HEPUtils::Particle * > >getSFOSpairs(std::vector< const HEPUtils::Particle * > particles)
std::vector< std::vector< const HEPUtils::Particle * > >getOSpairs(std::vector< const HEPUtils::Particle * > particles)
std::vector< std::vector< const HEPUtils::Particle * > >getSSpairs(std::vector< const HEPUtils::Particle * > particles)

Attributes

Name
std::map< str, DMEFT_analysis_info >analysis_info_map

Functions Documentation

function LEP208_SLHA1_convention_xsec_selselbar

void LEP208_SLHA1_convention_xsec_selselbar(
    triplet< double > & result
)

ee –> selectron pair production cross-sections at 208 GeV

function LEP208_SLHA1_convention_xsec_selserbar

void LEP208_SLHA1_convention_xsec_selserbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_serserbar

void LEP208_SLHA1_convention_xsec_serserbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_serselbar

void LEP208_SLHA1_convention_xsec_serselbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_se1se1bar

void LEP208_SLHA1_convention_xsec_se1se1bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_se1se2bar

void LEP208_SLHA1_convention_xsec_se1se2bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_se2se2bar

void LEP208_SLHA1_convention_xsec_se2se2bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_se2se1bar

void LEP208_SLHA1_convention_xsec_se2se1bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smulsmulbar

void LEP208_SLHA1_convention_xsec_smulsmulbar(
    triplet< double > & result
)

ee –> smuon pair production cross-sections at 208 GeV

function LEP208_SLHA1_convention_xsec_smulsmurbar

void LEP208_SLHA1_convention_xsec_smulsmurbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smursmurbar

void LEP208_SLHA1_convention_xsec_smursmurbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smursmulbar

void LEP208_SLHA1_convention_xsec_smursmulbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smu1smu1bar

void LEP208_SLHA1_convention_xsec_smu1smu1bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smu1smu2bar

void LEP208_SLHA1_convention_xsec_smu1smu2bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smu2smu2bar

void LEP208_SLHA1_convention_xsec_smu2smu2bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_smu2smu1bar

void LEP208_SLHA1_convention_xsec_smu2smu1bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_staulstaulbar

void LEP208_SLHA1_convention_xsec_staulstaulbar(
    triplet< double > & result
)

ee –> stau pair production cross-sections at 208 GeV

function LEP208_SLHA1_convention_xsec_staulstaurbar

void LEP208_SLHA1_convention_xsec_staulstaurbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_staurstaurbar

void LEP208_SLHA1_convention_xsec_staurstaurbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_staurstaulbar

void LEP208_SLHA1_convention_xsec_staurstaulbar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_stau1stau1bar

void LEP208_SLHA1_convention_xsec_stau1stau1bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_stau1stau2bar

void LEP208_SLHA1_convention_xsec_stau1stau2bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_stau2stau2bar

void LEP208_SLHA1_convention_xsec_stau2stau2bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_stau2stau1bar

void LEP208_SLHA1_convention_xsec_stau2stau1bar(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_11

void LEP208_SLHA1_convention_xsec_chi00_11(
    triplet< double > & result
)

ee –> neutralino pair production cross-sections at 208 GeV

function LEP208_SLHA1_convention_xsec_chi00_12

void LEP208_SLHA1_convention_xsec_chi00_12(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_13

void LEP208_SLHA1_convention_xsec_chi00_13(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_14

void LEP208_SLHA1_convention_xsec_chi00_14(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_22

void LEP208_SLHA1_convention_xsec_chi00_22(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_23

void LEP208_SLHA1_convention_xsec_chi00_23(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_24

void LEP208_SLHA1_convention_xsec_chi00_24(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_33

void LEP208_SLHA1_convention_xsec_chi00_33(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_34

void LEP208_SLHA1_convention_xsec_chi00_34(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chi00_44

void LEP208_SLHA1_convention_xsec_chi00_44(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chipm_11

void LEP208_SLHA1_convention_xsec_chipm_11(
    triplet< double > & result
)

ee –> chargino pair production cross-sections at 208 GeV

function LEP208_SLHA1_convention_xsec_chipm_12

void LEP208_SLHA1_convention_xsec_chipm_12(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chipm_22

void LEP208_SLHA1_convention_xsec_chipm_22(
    triplet< double > & result
)

function LEP208_SLHA1_convention_xsec_chipm_21

void LEP208_SLHA1_convention_xsec_chipm_21(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_selselbar

void LEP205_SLHA1_convention_xsec_selselbar(
    triplet< double > & result
)

ee –> selectron pair production cross-sections at 205 GeV

function LEP205_SLHA1_convention_xsec_selserbar

void LEP205_SLHA1_convention_xsec_selserbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_serserbar

void LEP205_SLHA1_convention_xsec_serserbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_serselbar

void LEP205_SLHA1_convention_xsec_serselbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_se1se1bar

void LEP205_SLHA1_convention_xsec_se1se1bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_se1se2bar

void LEP205_SLHA1_convention_xsec_se1se2bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_se2se2bar

void LEP205_SLHA1_convention_xsec_se2se2bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_se2se1bar

void LEP205_SLHA1_convention_xsec_se2se1bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smulsmulbar

void LEP205_SLHA1_convention_xsec_smulsmulbar(
    triplet< double > & result
)

ee –> smuon pair production cross-sections at 205 GeV

function LEP205_SLHA1_convention_xsec_smulsmurbar

void LEP205_SLHA1_convention_xsec_smulsmurbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smursmurbar

void LEP205_SLHA1_convention_xsec_smursmurbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smursmulbar

void LEP205_SLHA1_convention_xsec_smursmulbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smu1smu1bar

void LEP205_SLHA1_convention_xsec_smu1smu1bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smu1smu2bar

void LEP205_SLHA1_convention_xsec_smu1smu2bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smu2smu2bar

void LEP205_SLHA1_convention_xsec_smu2smu2bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_smu2smu1bar

void LEP205_SLHA1_convention_xsec_smu2smu1bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_staulstaulbar

void LEP205_SLHA1_convention_xsec_staulstaulbar(
    triplet< double > & result
)

ee –> stau pair production cross-sections at 205 GeV

function LEP205_SLHA1_convention_xsec_staulstaurbar

void LEP205_SLHA1_convention_xsec_staulstaurbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_staurstaurbar

void LEP205_SLHA1_convention_xsec_staurstaurbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_staurstaulbar

void LEP205_SLHA1_convention_xsec_staurstaulbar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_stau1stau1bar

void LEP205_SLHA1_convention_xsec_stau1stau1bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_stau1stau2bar

void LEP205_SLHA1_convention_xsec_stau1stau2bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_stau2stau2bar

void LEP205_SLHA1_convention_xsec_stau2stau2bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_stau2stau1bar

void LEP205_SLHA1_convention_xsec_stau2stau1bar(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_11

void LEP205_SLHA1_convention_xsec_chi00_11(
    triplet< double > & result
)

ee –> neutralino pair production cross-sections at 205 GeV

function LEP205_SLHA1_convention_xsec_chi00_12

void LEP205_SLHA1_convention_xsec_chi00_12(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_13

void LEP205_SLHA1_convention_xsec_chi00_13(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_14

void LEP205_SLHA1_convention_xsec_chi00_14(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_22

void LEP205_SLHA1_convention_xsec_chi00_22(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_23

void LEP205_SLHA1_convention_xsec_chi00_23(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_24

void LEP205_SLHA1_convention_xsec_chi00_24(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_33

void LEP205_SLHA1_convention_xsec_chi00_33(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_34

void LEP205_SLHA1_convention_xsec_chi00_34(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chi00_44

void LEP205_SLHA1_convention_xsec_chi00_44(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chipm_11

void LEP205_SLHA1_convention_xsec_chipm_11(
    triplet< double > & result
)

ee –> chargino pair production cross-sections at 205 GeV

function LEP205_SLHA1_convention_xsec_chipm_12

void LEP205_SLHA1_convention_xsec_chipm_12(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chipm_22

void LEP205_SLHA1_convention_xsec_chipm_22(
    triplet< double > & result
)

function LEP205_SLHA1_convention_xsec_chipm_21

void LEP205_SLHA1_convention_xsec_chipm_21(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_selselbar

void LEP188_SLHA1_convention_xsec_selselbar(
    triplet< double > & result
)

ee –> selectron pair production cross-sections at 188.6 GeV

function LEP188_SLHA1_convention_xsec_selserbar

void LEP188_SLHA1_convention_xsec_selserbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_serserbar

void LEP188_SLHA1_convention_xsec_serserbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_serselbar

void LEP188_SLHA1_convention_xsec_serselbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_se1se1bar

void LEP188_SLHA1_convention_xsec_se1se1bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_se1se2bar

void LEP188_SLHA1_convention_xsec_se1se2bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_se2se2bar

void LEP188_SLHA1_convention_xsec_se2se2bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_se2se1bar

void LEP188_SLHA1_convention_xsec_se2se1bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smulsmulbar

void LEP188_SLHA1_convention_xsec_smulsmulbar(
    triplet< double > & result
)

ee –> smuon pair production cross-sections at 188.6 GeV

function LEP188_SLHA1_convention_xsec_smulsmurbar

void LEP188_SLHA1_convention_xsec_smulsmurbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smursmurbar

void LEP188_SLHA1_convention_xsec_smursmurbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smursmulbar

void LEP188_SLHA1_convention_xsec_smursmulbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smu1smu1bar

void LEP188_SLHA1_convention_xsec_smu1smu1bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smu1smu2bar

void LEP188_SLHA1_convention_xsec_smu1smu2bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smu2smu2bar

void LEP188_SLHA1_convention_xsec_smu2smu2bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_smu2smu1bar

void LEP188_SLHA1_convention_xsec_smu2smu1bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_staulstaulbar

void LEP188_SLHA1_convention_xsec_staulstaulbar(
    triplet< double > & result
)

ee –> stau pair production cross-sections at 188.6 GeV

function LEP188_SLHA1_convention_xsec_staulstaurbar

void LEP188_SLHA1_convention_xsec_staulstaurbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_staurstaurbar

void LEP188_SLHA1_convention_xsec_staurstaurbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_staurstaulbar

void LEP188_SLHA1_convention_xsec_staurstaulbar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_stau1stau1bar

void LEP188_SLHA1_convention_xsec_stau1stau1bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_stau1stau2bar

void LEP188_SLHA1_convention_xsec_stau1stau2bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_stau2stau2bar

void LEP188_SLHA1_convention_xsec_stau2stau2bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_stau2stau1bar

void LEP188_SLHA1_convention_xsec_stau2stau1bar(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_11

void LEP188_SLHA1_convention_xsec_chi00_11(
    triplet< double > & result
)

ee –> neutralino pair production cross-sections at 188.6 GeV

function LEP188_SLHA1_convention_xsec_chi00_12

void LEP188_SLHA1_convention_xsec_chi00_12(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_13

void LEP188_SLHA1_convention_xsec_chi00_13(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_14

void LEP188_SLHA1_convention_xsec_chi00_14(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_22

void LEP188_SLHA1_convention_xsec_chi00_22(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_23

void LEP188_SLHA1_convention_xsec_chi00_23(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_24

void LEP188_SLHA1_convention_xsec_chi00_24(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_33

void LEP188_SLHA1_convention_xsec_chi00_33(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_34

void LEP188_SLHA1_convention_xsec_chi00_34(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chi00_44

void LEP188_SLHA1_convention_xsec_chi00_44(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chipm_11

void LEP188_SLHA1_convention_xsec_chipm_11(
    triplet< double > & result
)

ee –> chargino pair production cross-sections at 188.6 GeV

function LEP188_SLHA1_convention_xsec_chipm_12

void LEP188_SLHA1_convention_xsec_chipm_12(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chipm_22

void LEP188_SLHA1_convention_xsec_chipm_22(
    triplet< double > & result
)

function LEP188_SLHA1_convention_xsec_chipm_21

void LEP188_SLHA1_convention_xsec_chipm_21(
    triplet< double > & result
)

function ALEPH_Selectron_Conservative_LLike

void ALEPH_Selectron_Conservative_LLike(
    double & result
)

LEP Slepton Log-Likelihoods

function ALEPH_Smuon_Conservative_LLike

void ALEPH_Smuon_Conservative_LLike(
    double & result
)

function ALEPH_Stau_Conservative_LLike

void ALEPH_Stau_Conservative_LLike(
    double & result
)

function L3_Selectron_Conservative_LLike

void L3_Selectron_Conservative_LLike(
    double & result
)

function L3_Smuon_Conservative_LLike

void L3_Smuon_Conservative_LLike(
    double & result
)

function L3_Stau_Conservative_LLike

void L3_Stau_Conservative_LLike(
    double & result
)

function L3_Neutralino_All_Channels_Conservative_LLike

void L3_Neutralino_All_Channels_Conservative_LLike(
    double & result
)

LEP Gaugino Log-Likelihoods

function L3_Neutralino_Leptonic_Conservative_LLike

void L3_Neutralino_Leptonic_Conservative_LLike(
    double & result
)

function L3_Chargino_All_Channels_Conservative_LLike

void L3_Chargino_All_Channels_Conservative_LLike(
    double & result
)

function L3_Chargino_Leptonic_Conservative_LLike

void L3_Chargino_Leptonic_Conservative_LLike(
    double & result
)

function OPAL_Chargino_Hadronic_Conservative_LLike

void OPAL_Chargino_Hadronic_Conservative_LLike(
    double & result
)

function OPAL_Chargino_SemiLeptonic_Conservative_LLike

void OPAL_Chargino_SemiLeptonic_Conservative_LLike(
    double & result
)

function OPAL_Chargino_Leptonic_Conservative_LLike

void OPAL_Chargino_Leptonic_Conservative_LLike(
    double & result
)

function OPAL_Degenerate_Chargino_Conservative_LLike

void OPAL_Degenerate_Chargino_Conservative_LLike(
    double & result
)

function OPAL_Chargino_All_Channels_Conservative_LLike

void OPAL_Chargino_All_Channels_Conservative_LLike(
    double & result
)

function OPAL_Neutralino_Hadronic_Conservative_LLike

void OPAL_Neutralino_Hadronic_Conservative_LLike(
    double & result
)

function I1

double I1(
    double s,
    double m1,
    double m2,
    double mk,
    double ml
)

Integrals for t-channel neutralino diagrams m1 and m2 are masses of final state sleptons mk and ml are neutralino masses

function I2

double I2(
    double s,
    double m1,
    double m2,
    double mk,
    double ml
)

function I3

double I3(
    double s,
    double m1,
    double m2,
    double mk
)

function SLHA2BFM_NN

void SLHA2BFM_NN(
    MixMatrix & NN,
    double tanb,
    double sin2thetaW
)

Functions to convert mass matrices between SLHA and BFM conventions ////////////////////////////////////////////////////////////////////.

Converts a neutralino mixing matrix in SLHA conventions to BFM conventions, (\tan\beta) is as defined in SLHA

function SLHA2BFM_VV

void SLHA2BFM_VV(
    MixMatrix & VV
)

Converts the chargino mixing matrix V in SLHA conventions to BFM conventions.

function BFM2SLHA_NN

void BFM2SLHA_NN(
    MixMatrix & NN,
    double tanb,
    double sin2thetaW
)

Converts a neutralino mixing matrix in BFM conventions to SLHA conventions, (\tan\beta) is as defined in SLHA.

function BFM2SLHA_VV

void BFM2SLHA_VV(
    MixMatrix & VV
)

Converts the chargino mixing matrix V in BFM conventions to SLHA conventions.

function multiply

MixMatrix multiply(
    MixMatrix A,
    MixMatrix B
)

Helper function to multiply matrices.

function transpose

MixMatrix transpose(
    MixMatrix A
)

Helper function to find matrix transpose.

function print

void print(
    MixMatrix A
)

Helper function to print a matrix.

function sortByPT13

bool sortByPT13(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByPT13_sharedptr

bool sortByPT13_sharedptr(
    std::shared_ptr< HEPUtils::Jet > jet1,
    std::shared_ptr< HEPUtils::Jet > jet2
)

function sortByMass

bool sortByMass(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByMass_sharedptr

bool sortByMass_sharedptr(
    std::shared_ptr< HEPUtils::Jet > jet1,
    std::shared_ptr< HEPUtils::Jet > jet2
)

function calcMT

double calcMT(
    HEPUtils::P4 jetMom,
    HEPUtils::P4 metMom
)

function sortByPT_1l

bool sortByPT_1l(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByPT_1l_sharedptr

bool sortByPT_1l_sharedptr(
    std::shared_ptr< HEPUtils::Jet > jet1,
    std::shared_ptr< HEPUtils::Jet > jet2
)

function sortByMass_1l

bool sortByMass_1l(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByMass_1l_sharedptr

bool sortByMass_1l_sharedptr(
    std::shared_ptr< HEPUtils::Jet > jet1,
    std::shared_ptr< HEPUtils::Jet > jet2
)

function calcMT_1l

double calcMT_1l(
    HEPUtils::P4 jetMom,
    HEPUtils::P4 metMom
)

function sortByPT_j

bool sortByPT_j(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByPT_l

bool sortByPT_l(
    const HEPUtils::Particle * lep1,
    const HEPUtils::Particle * lep2
)

function sortByPT_jet

bool sortByPT_jet(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByPT_lep

bool sortByPT_lep(
    const HEPUtils::Particle * lep1,
    const HEPUtils::Particle * lep2
)

function sortByPT_RJ3L

bool sortByPT_RJ3L(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortLepByPT_RJ3L

bool sortLepByPT_RJ3L(
    const HEPUtils::Particle * lep1,
    const HEPUtils::Particle * lep2
)

function SortLeptons

bool SortLeptons(
    const pair< TLorentzVector, int > lv1,
    const pair< TLorentzVector, int > lv2
)

function SortJets

bool SortJets(
    const TLorentzVector jv1,
    const TLorentzVector jv2
)

function DEFINE_ANALYSIS_FACTORY

DEFINE_ANALYSIS_FACTORY(
    ATLAS_13TeV_ZGammaGrav_CONFNOTE_80invfb 
)

function sortByPT

bool sortByPT(
    const HEPUtils::Jet * jet1,
    const HEPUtils::Jet * jet2
)

function sortByPT_2lep

bool sortByPT_2lep(
    const HEPUtils::Particle * lep1,
    const HEPUtils::Particle * lep2
)

function Phi_mpi_pi

double Phi_mpi_pi(
    double x
)

function _Phi_mpi_pi

double _Phi_mpi_pi(
    double x
)

function mkAnalysis

Analysis * mkAnalysis(
    const str & name
)

Factory definition.

function getDetector

str getDetector(
    const str & name
)

Return the detector to be used for a given analysis name (and check that the analysis exists).

function _setEventWeight_unity

void _setEventWeight_unity(
    HEPUtils::Event & event,
    const BaseCollider * 
)

A function that sets the event weight to unity, with zero uncertainty.

function setEventWeight_unity

void setEventWeight_unity(
    EventWeighterFunctionType & result
)

Module function providing an instance of EventWeighterFunctionType pointing to _setEventWeight_unity

function _setEventWeight_fromCrossSection

void _setEventWeight_fromCrossSection(
    HEPUtils::Event & event,
    const BaseCollider * HardScatteringSim_ptr,
    const map_int_process_xsec & ProcessCrossSectionsMap,
    const int use_trust_level
)

A function that sets the event weight based on the process cross-sections.

function setEventWeight_fromCrossSection

void setEventWeight_fromCrossSection(
    EventWeighterFunctionType & result
)

Module function providing an instance of EventWeighterFunctionType pointing to _setEventWeight_fromCrossSection

function getDummyColliderObservable

void getDummyColliderObservable(
    double & result
)

Dummy observable that creates a dependency on TestModel1D.

This is used to satisfy the normal GAMBIT model requrements in a minimal way. This is useful in the case where we just want to run ColliderBit on a single point with a custom Pythia version, using Pythia’s SLHA interface.

function operateLHCLoop

void operateLHCLoop(
    MCLoopInfo & result
)

LHC Loop Manager.

function getLHCEventLoopInfo

void getLHCEventLoopInfo(
    map_str_dbl & result
)

Store some information about the event generation.

function CollectAnalyses

void CollectAnalyses(
    AnalysisDataPointers & result
)

Loop over all analyses and collect them in one place.

function set_CS

void set_CS(
    hb_ModelParameters & result,
    const HiggsCouplingsTable & couplings,
    int n_neutral_higgses
)

Helper function to set HiggsBounds/Signals parameters cross-section ratios from a GAMBIT HiggsCouplingsTable.

function set_SMLikeHiggs_ModelParameters

void set_SMLikeHiggs_ModelParameters(
    const SubSpectrum & spec,
    const HiggsCouplingsTable & couplings,
    hb_ModelParameters & result
)

Helper function for populating a HiggsBounds/Signals ModelParameters object for SM-like Higgs.

function SMLikeHiggs_ModelParameters

void SMLikeHiggs_ModelParameters(
    hb_ModelParameters & result
)

SM-like (SM + possible invisibles) Higgs model parameters for HiggsBounds/Signals.

function MSSMLikeHiggs_ModelParameters

void MSSMLikeHiggs_ModelParameters(
    hb_ModelParameters & result
)

MSSM-like (MSSM + NMSSM + …) Higgs model parameters for HiggsBounds/Signals.

function calc_HB_LEP_LogLike

void calc_HB_LEP_LogLike(
    double & result
)

Get a LEP chisq from HiggsBounds.

function calc_HS_LHC_LogLike

void calc_HS_LHC_LogLike(
    double & result
)

Get an LHC chisq from HiggsSignals.

function FeynHiggs_HiggsProd

void FeynHiggs_HiggsProd(
    fh_HiggsProd_container & result
)

Higgs production cross-sections from FeynHiggs.

function calc_loglikes_for_analysis

AnalysisLogLikes calc_loglikes_for_analysis(
    const AnalysisData & adata,
    bool USE_COVAR,
    bool USE_MARG,
    bool combine_nocovar_SRs,
    bool set_zero_loglike
)

Forward declaration of funtion in LHC_likelihoods.

For a given analysis, calculate per-SR loglikes and the overall analysis loglike. Return the results as an AnalysLogLikes object.

If (simplified) SR-correlation info is available, so use the covariance matrix to construct composite marginalised likelihood Despite initial thoughts, we can’t just do independent LL calculations in a rotated basis, but have to sample from the covariance matrix.

This means we can’t use the nulike LL functions, which operate in 1D only. Also, log-normal sampling in the diagonal basis is not helpful, since the rotation will re-generate negative rates.

TodoSupport NSL, i.e. skewness correction

TodoUnify this for both cov and no-cov, feeding in one-element Eigen blocks as Ref<>s for the latter?

TodoCompute the background-only covariance decomposition and likelihood only once

TodoOnly compute this once per run

TodoUse newer (?) one-step Eigen constructors for (const) single-element arrays

TodoOr compute all the exp DLLs first, then only the best-expected SR’s obs DLL?

TodoTodoOnly compute this once per run

TodoOnly compute this once per run

If (simplified) SR-correlation info is available, so use the covariance matrix to construct composite marginalised likelihood Despite initial thoughts, we can’t just do independent LL calculations in a rotated basis, but have to sample from the covariance matrix.

This means we can’t use the nulike LL functions, which operate in 1D only. Also, log-normal sampling in the diagonal basis is not helpful, since the rotation will re-generate negative rates.

TodoSupport NSL, i.e. skewness correction

TodoUnify this for both cov and no-cov, feeding in one-element Eigen blocks as Ref<>s for the latter?

TodoCompute the background-only covariance decomposition and likelihood only once

TodoOnly compute this once per run

TodoUse newer (?) one-step Eigen constructors for (const) single-element arrays

TodoOr compute all the exp DLLs first, then only the best-expected SR’s obs DLL?

TodoTodoOnly compute this once per run

TodoOnly compute this once per run

function fill_analysis_info_map

void fill_analysis_info_map()

Forward declarations of functions in this file.

A function for filling the analysis_info_map. This is where all the analysis-specific numbers and file names go.

function DMEFT_results

void DMEFT_results(
    AnalysisDataPointers & result
)

Results from DMEFT analyses before any modification of the MET spectrum.

function get_all_DMEFT_signal_yields

void get_all_DMEFT_signal_yields(
    std::vector< double > & sr_nums,
    const DMEFT_analysis_info & analysis_info,
    const Spectrum & spec
)

Fill the input vector with the total DMEFT signal prediction for each SR in the given LHC analysis.

function get_DMEFT_signal_yields_dim6_operator

void get_DMEFT_signal_yields_dim6_operator(
    std::vector< double > & signal_yields,
    const str,
    const DMEFT_analysis_info & analysis_info,
    double m,
    double O1,
    double O2,
    double lambda
)

Fill the input vector with the DMEFT signal prediction for a given set of dim-6 operators.

function get_DMEFT_signal_yields_dim7_operator

void get_DMEFT_signal_yields_dim7_operator(
    std::vector< double > & signal_yields,
    const str,
    const DMEFT_analysis_info & analysis_info,
    double m,
    double O,
    double lambda
)

Fill the input vector with the DMEFT signal prediction for a given dim-7 operator.

function DMEFT_results_profiled

void DMEFT_results_profiled(
    AnalysisDataPointers & result
)

Results from DMEFT analyses after profiling over the ‘a’ parameter in the smooth cut-off of the MET spectrum.

function DMEFT_results_cutoff

void DMEFT_results_cutoff(
    AnalysisDataPointers & result
)

Results from DMEFT analyses after imposing a hard cut-off of the MET spectrum.

function signal_modifier_function

void signal_modifier_function(
    AnalysisData & adata,
    double lambda,
    double a
)

Function to modify the DMEFT LHC signal prediction for ETmiss bins where ETmiss > Lambda. Alt 1: Gradually turn off the ETmiss spectrum above Lambda by multiplying the spectrum with (ETmiss/Lambda)^-a

function signal_cutoff_function

void signal_cutoff_function(
    AnalysisData & adata,
    double lambda
)

Function to modify the DMEFT LHC signal prediction for ETmiss bins where ETmiss > Lambda. Alt 2: Simply put a hard cut-off in the ETmiss spectrum for ETmiss > Lambda

function _gsl_target_func

void _gsl_target_func(
    const size_t,
    const double * a,
    void * fparams,
    double * fval
)

A target function for the GSL optimiser.

function calc_DMEFT_profiled_LHC_nuisance_params

void calc_DMEFT_profiled_LHC_nuisance_params(
    map_str_dbl & result
)

function InterpolatedMCInfo

void InterpolatedMCInfo(
    MCLoopInfo & result
)

This makes an MCLoopInfo object for satisfying the ColliderBit dependency chain (This will not be needed once we have a general system for simulation-less analyses.)

function limit_LLike

double limit_LLike(
    double x,
    double x95,
    double sigma
)

LEP limit likelihood function.

Incorporate theoretical uncertainty in a 95% limit xPredicted cross section

x95Experimental 95% upper limit on cross section

sigmaTheoretical uncertainty on predicted cross section

Log-likelihood

function is_xsec_sane

bool is_xsec_sane(
    const triplet< double > & xsecWithError
)

LEP limit debugging function.

function LEP207_SLHA1_convention_xsec_chi00_11

void LEP207_SLHA1_convention_xsec_chi00_11(
    triplet< double > & result
)

function L3_Gravitino_LLike

void L3_Gravitino_LLike(
    double & result
)

L3 search for gravitinos at 207 GeV

We use a limit from Fig. 6c of https://doi.org/10.1016/j.physletb.2004.01.010.

We use the 95% upper limit on [ \sigma(ee \to \chi^0_1\chi^0_1) \textrm{BR}(\chi^0_1 \to \tilde{G}\gamma)^2 ]

function all_PID_pairs_to_process_codes

const multimap_PID_pair_int & all_PID_pairs_to_process_codes()

function getActiveProcessCodes

void getActiveProcessCodes(
    std::vector< int > & result
)

Get the list of active collider process codes.

function getAnalysisContainer

void getAnalysisContainer(
    AnalysisContainer & result,
    const str & detname,
    const MCLoopInfo & RunMC,
    const xsec_container & TotalCrossSection,
    int iteration
)

Retrieve an analysis container for a specific detector.

function getBuckFastATLAS

void getBuckFastATLAS(
    BaseDetector *& result
)

Retrieve a BuckFast sim of ATLAS.

function getBuckFastCMS

void getBuckFastCMS(
    BaseDetector *& result
)

Retrieve a BuckFast sim of CMS.

function getBuckFastIdentity

void getBuckFastIdentity(
    BaseDetector *& result
)

Retrieve an Identity BuckFast sim (no sim)

function getHepMCEvent

void getHepMCEvent(
    HEPUtils::Event & result
)

A nested function that reads in HepMC event files and converts them to HEPUtils::Event format.

function getLHEvent

void getLHEvent(
    HEPUtils::Event & result
)

A nested function that reads in Les Houches Event files and converts them to HEPUtils::Event format.

function convert_xsecs_to_fb

std::pair< double, double > convert_xsecs_to_fb(
    double input_xsec,
    double input_xsec_uncert,
    str input_unit,
    bool input_fractional_uncert
)

Helper function that takes a cross-section value in fb or pb, along with an absolute or relative uncertainty, and returns the xsec and absolute uncertainty in fb.

function getPIDPairCrossSectionsMap_prospino

void getPIDPairCrossSectionsMap_prospino(
    map_PID_pair_PID_pair_xsec & result
)

Get a cross-section from Prospino WORK IN PROGRESS

function silly_pid_xsec_constructor

PID_pair_xsec_container silly_pid_xsec_constructor(
    PID_pair pid_pair,
    double xsec_val
)

Test functions for provding PIDPairCrossSectionsMap (cross-sections in fb)

function getPIDPairCrossSectionsMap_testing

void getPIDPairCrossSectionsMap_testing(
    map_PID_pair_PID_pair_xsec & result
)

function getProcessCrossSectionsMap

void getProcessCrossSectionsMap(
    map_int_process_xsec & result
)

Get a map between Pythia process codes and cross-sections.

function getEvGenCrossSection

void getEvGenCrossSection(
    MC_xsec_container & result
)

Compute a cross-section from Monte Carlo.

function getEvGenCrossSection_as_base

void getEvGenCrossSection_as_base(
    xsec_container & result
)

Return MC_xsec_container as the base xsec_container.

function getNLLFastCrossSection

void getNLLFastCrossSection(
    xsec_container & result
)

Get a cross-section from NLL-FAST.

function checkOptions_getYAMLCrossSection

bool checkOptions_getYAMLCrossSection(
    const Options & runOptions,
    const str calling_function,
    std::pair< str, str > & xsec_pnames,
    str & input_unit,
    bool & input_fractional_uncert,
    str & errmsg
)

A helper function to check the YAML options for getYAMLCrossSection and getYAMLCrossSection_SLHA.

function getYAMLCrossSection

void getYAMLCrossSection(
    xsec_container & result
)

A function that reads the total cross-section from the input file, but builds up the number of events from the event loop.

function getYAMLCrossSection_SLHA

void getYAMLCrossSection_SLHA(
    xsec_container & result
)

A function that reads a list of (SLHA file, total cross-section) pairs from the input YAML file.

function getYAMLCrossSection_param

void getYAMLCrossSection_param(
    xsec_container & result
)

A function that assigns a total cross-sections directly from the scan parameters (for model ColliderBit_SLHA_scan_model)

function getTotalCrossSectionAsMap

void getTotalCrossSectionAsMap(
    map_str_dbl & result
)

Get cross-section info as map_str_dbl (for simple printing)

function getPIDPairCrossSectionsInfo

void getPIDPairCrossSectionsInfo(
    map_str_dbl & result
)

Output PID pair cross-sections as a str-dbl map, for easy printing.

function doCrossSectionConsistencyCheck

void doCrossSectionConsistencyCheck(
    bool & result
)

A consistency check that ensures that if each event is weighted by a process-level cross-section from an external calculator, then the total cross-section is taken from the event generator

function get_sigma_ee_ll

void get_sigma_ee_ll(
    triplet< double > & result,
    const double sqrts,
    const int generation,
    const int l_chirality,
    const int lbar_chirality,
    const double gtol,
    const double ftol,
    const bool gpt_error,
    const bool fpt_error,
    const Spectrum & spec,
    const double gammaZ,
    const bool l_are_gauge_es
)

Retrieve the production cross-section at an e+e- collider for slepton pairs. If l_are_gauge_es = T, then l(bar)_chirality = 1 => (anti-)left-type slepton = 2 => (anti-)right-type slepton If l_are_gauge_es = F, then l(bar)_chirality = 1 => (anti-)slepton is lightest family state = 2 => (anti-)slepton is heaviest family state

function get_sigma_ee_chi00

void get_sigma_ee_chi00(
    triplet< double > & result,
    const double sqrts,
    const int chi_first,
    const int chi_second,
    const double tol,
    const bool pt_error,
    const Spectrum & spec,
    const double gammaZ
)

Retrieve the production cross-section at an e+e- collider for neutralino pairs.

function get_sigma_ee_chipm

void get_sigma_ee_chipm(
    triplet< double > & result,
    const double sqrts,
    const int chi_plus,
    const int chi_minus,
    const double tol,
    const bool pt_error,
    const Spectrum & spec,
    const double gammaZ
)

Retrieve the production cross-section at an e+e- collider for chargino pairs.

function xsec_sleislej

double xsec_sleislej(
    int pid1,
    int pid2,
    double sqrts,
    double m1,
    double m2,
    MixMatrix F,
    MixMatrix N,
    const double mN[4],
    double alpha,
    double mZ,
    double gZ,
    double sin2thetaW,
    bool CP_lock
)

Cross section [pb] for (e^+e^- -> \tilde l_i \tilde l_j^*) To use, call SLHA2BFM first on SLHA mixing matrices constructed as a vector of vectors

function xsec_neuineuj

double xsec_neuineuj(
    int pid1,
    int pid2,
    double sqrts,
    double mi,
    double mj,
    MixMatrix N,
    const double mS[2],
    double tanb,
    double alpha,
    double mZ,
    double gZ,
    double sin2thetaW
)

Cross section [pb] for (e^+e^- -> \tilde\chi^0_i \tilde\chi^0_j) Masses mi and mj for the neutralinos are signed. mS are the selectron masses (left = 0, right = 1). Warning! BFM uses inverted (\tan\beta)! Use tanb = 1 / tanb in converting from SLHA.

function xsec_chaichaj

double xsec_chaichaj(
    int pid1,
    int pid2,
    double sqrts,
    double mi,
    double mj,
    MixMatrix V,
    MixMatrix U,
    double ms,
    double alpha,
    double mZ,
    double gZ,
    double sin2thetaW
)

Cross section [pb] for (e^+e^- -> \tilde\chi^+_i \tilde\chi^-_j) Masses mi and mj for the charginos are signed. msn is electron sneutrino mass.

function calc_LHC_signals

void calc_LHC_signals(
    map_str_dbl & result
)

Loop over all analyses and fill a map of predicted counts.

function _gsl_calc_Analysis_MinusLogLike

void _gsl_calc_Analysis_MinusLogLike(
    const size_t n,
    const double * unit_nuisances_dbl,
    void * fixedparamspack,
    double * fval
)

Note: Doesn’t return a full log-like: the factorial term is missing since it’s expensive, fixed and cancels in DLLs

Loglike objective-function wrapper to provide the signature for GSL multimin

We’ve dropped the log(n_obs!) terms, since they’re expensive and cancel in computing DLL

function _gsl_calc_Analysis_MinusLogLikeGrad

void _gsl_calc_Analysis_MinusLogLikeGrad(
    const size_t n,
    const double * unit_nuisances_dbl,
    void * fixedparamspack,
    double * fgrad
)

Loglike gradient-function wrapper to provide the signature for GSL multimin.

function _gsl_calc_Analysis_MinusLogLikeAndGrad

void _gsl_calc_Analysis_MinusLogLikeAndGrad(
    const size_t n,
    const double * unit_nuisances_dbl,
    void * fixedparamspack,
    double * fval,
    double * fgrad
)

function _gsl_mkpackedarray

std::vector< double > _gsl_mkpackedarray(
    const Eigen::ArrayXd & n_preds,
    const Eigen::ArrayXd & n_obss,
    const Eigen::ArrayXd & sqrtevals,
    const Eigen::MatrixXd & evecs
)

function profile_loglike_cov

double profile_loglike_cov(
    const Eigen::ArrayXd & n_preds,
    const Eigen::ArrayXd & n_obss,
    const Eigen::ArrayXd & sqrtevals,
    const Eigen::MatrixXd & evecs
)

Note: Return value is missing the log(n_obs!) terms (n_SR of them) which cancel in LLR calculation

Todo: Pass in the cov, and compute the fixed evals, evecs, and corr matrix as fixed params in here? Via a helper function to reduce duplication

Return the best log likelihood

function marg_loglike_nulike1sr

double marg_loglike_nulike1sr(
    const Eigen::ArrayXd & n_preds,
    const Eigen::ArrayXd & n_obss,
    const Eigen::ArrayXd & sqrtevals
)

function marg_loglike_cov

double marg_loglike_cov(
    const Eigen::ArrayXd & n_preds,
    const Eigen::ArrayXd & n_obss,
    const Eigen::ArrayXd & sqrtevals,
    const Eigen::MatrixXd & evecs
)

TodoShould also implement a check of relative difference

How to correct negative rates? Discard (scales badly), set to epsilon (= discontinuous & unphysical pdf), transform to log-space (distorts the pdf quite badly), or something else (skew term)? We’re using the “set to epsilon” version for now. Ben: I would vote for ‘discard’. It can’t be that inefficient, surely? Andy: For a lot of signal regions, the probability of none having a negative sample is Prod_SR p(non-negative)_SR… which can get bad.

function calc_LHC_LogLikes

void calc_LHC_LogLikes(
    map_str_AnalysisLogLikes & result
)

Loop over all analyses and fill a map of AnalysisLogLikes objects.

function get_LHC_LogLike_per_analysis

void get_LHC_LogLike_per_analysis(
    map_str_dbl & result
)

Extract the combined log likelihood for each analysis.

function get_LHC_LogLike_per_SR

void get_LHC_LogLike_per_SR(
    map_str_dbl & result
)

Extract the log likelihood for each SR.

function get_LHC_LogLike_SR_labels

void get_LHC_LogLike_SR_labels(
    map_str_str & result
)

Extract the labels for the SRs used in the analysis loglikes.

function get_LHC_LogLike_SR_indices

void get_LHC_LogLike_SR_indices(
    map_str_dbl & result
)

Todo: Switch result type to map_str_int once we have implemented a printer for this type

Extract the indices for the SRs used in the analysis loglikes

function calc_combined_LHC_LogLike

void calc_combined_LHC_LogLike(
    double & result
)

Compute the total likelihood combining all analyses.

function calc_LHC_LogLike_scan_guide

void calc_LHC_LogLike_scan_guide(
    double & result
)

A dummy log-likelihood that helps the scanner track a given range of collider log-likelihood values

function getNextSLHAFileNameAndContent

void getNextSLHAFileNameAndContent(
    pair_str_SLHAstruct & result
)

function getAndReplaceSLHAContent

void getAndReplaceSLHAContent(
    pair_str_SLHAstruct & result
)

TodoAdd option to save the new SLHA content to file

function getSLHAFileElements

void getSLHAFileElements(
    map_str_dbl & result
)

function getSLHA1Spectrum

void getSLHA1Spectrum(
    SLHAstruct & result
)

function getSLHA2Spectrum

void getSLHA2Spectrum(
    SLHAstruct & result
)

function calc_susy_spectrum_scan_guide

void calc_susy_spectrum_scan_guide(
    double & result
)

function getActiveProcessCodeToPIDPairsMap

void getActiveProcessCodeToPIDPairsMap(
    multimap_int_PID_pair & result
)

Get a multimap between the active Pythia process codes and the corresponding PID pair for the two final state particles

function getActivePIDPairs

void getActivePIDPairs(
    vec_PID_pair & result
)

Get a list of all the PID pairs related to active process codes.

function runAnalyses

void runAnalyses(
    AnalysisDataPointers & result,
    const str & ,
    const MCLoopInfo & RunMC,
    const AnalysisContainer & Container,
    const HEPUtils::Event & SmearedEvent,
    int iteration,
    void(*)() wrapup
)

Run all the analyses in a given container.

function smearEvent

void smearEvent(
    HEPUtils::Event & result,
    const HEPUtils::Event & HardScatteringEvent,
    const BaseDetector & detector,
    const MCLoopInfo & RunMC,
    const int iteration,
    const str & detname
)

Smear an event.

function random_bool

bool random_bool(
    double eff
)

TodoHandle out-of-range eff values

function filtereff

void filtereff(
    std::vector< const HEPUtils::Particle * > & particles,
    double eff,
    bool do_delete
)

function filtereff

void filtereff(
    std::vector< const HEPUtils::Particle * > & particles,
    std::function< double(const HEPUtils::Particle *)> eff_fn,
    bool do_delete
)

Utility function for filtering a supplied particle vector by sampling wrt a binned 1D efficiency map in pT.

function filtereff_pt

void filtereff_pt(
    std::vector< const HEPUtils::Particle * > & particles,
    const HEPUtils::BinnedFn1D< double > & eff_pt,
    bool do_delete
)

function filtereff_etapt

void filtereff_etapt(
    std::vector< const HEPUtils::Particle * > & particles,
    const HEPUtils::BinnedFn2D< double > & eff_etapt,
    bool do_delete
)

function getSFOSpairs

std::vector< std::vector< const HEPUtils::Particle * > > getSFOSpairs(
    std::vector< const HEPUtils::Particle * > particles
)

function getOSpairs

std::vector< std::vector< const HEPUtils::Particle * > > getOSpairs(
    std::vector< const HEPUtils::Particle * > particles
)

function getSSpairs

std::vector< std::vector< const HEPUtils::Particle * > > getSSpairs(
    std::vector< const HEPUtils::Particle * > particles
)

Attributes Documentation

variable analysis_info_map

std::map< str, DMEFT_analysis_info > analysis_info_map;

A global map from analysis name to DMEFT_analysis_info instance. This map is initialized by the function fill_analysis_info_map, which is called the first time DMEFT_results run.


Updated on 2022-08-03 at 12:58:17 +0000