Page MenuHomeHEPForge
Files F8308543

HdecayWrapCppFuncs.cpp
No OneTemporary
Actions

HdecayWrapCppFuncs.cpp
View Options

#include "HdecayWrap.h"
using namespace std;
void HdecayCalc(void){
double dummy;
hdecaycallplain_(dummy);
}
void HdecayCalc2HDM(std::map<std::string,double> & HdecayA,std::map<std::string,double> & HdecayHlight,std::map<std::string,double> & HdecayHheavy,std::map<std::string,double> & HdecayHcharged,std::map<std::string,double> & HdecayTop){
double ABRB_out,ABRL_out,ABRM_out,ABRS_out,ABRC_out,ABRT_out,ABRG_out,ABRGA_out, ABRZGA_out,ABRZ_out,AWDTH_out,HLBRB_out,HLBRL_out,HLBRM_out,HLBRS_out,HLBRC_out,HLBRT_out,HLBRG_out,HLBRGA_out,HLBRZGA_out,HLBRW_out,HLBRZ_out,HLBRA_out,HLBRAZ_out,HLBRHW_out,HLWDTH_out,HHBRB_out,HHBRL_out,HHBRM_out,HHBRS_out,HHBRC_out,HHBRT_out,HHBRG_out,HHBRGA_out,HHBRZGA_out,HHBRW_out,HHBRZ_out,HHBRH_out,HHBRA_out,HHBRAZ_out,HHBRHW_out,HHWDTH_out,HCBRB_out,HCBRL_out,HCBRM_out,HCBRBU_out,HCBRS_out,HCBRC_out,HCBRT_out,HCBRW_out,HCBRA_out,HCWDTH_out,hcbrcd_out,hcbrts_out,hcbrtd_out,hcbrwhh_out,hhbrchch_out,hlbrchch_out,abrhhaz_out,abrhawphm_out,gamt0_out,gamt1_out;
hdecaycall2hdm_(ABRB_out,ABRL_out,ABRM_out,ABRS_out,ABRC_out,ABRT_out,ABRG_out,ABRGA_out, ABRZGA_out,ABRZ_out,AWDTH_out,HLBRB_out,HLBRL_out,HLBRM_out,HLBRS_out,HLBRC_out,HLBRT_out,HLBRG_out,HLBRGA_out,HLBRZGA_out,HLBRW_out,HLBRZ_out,HLBRA_out,HLBRAZ_out,HLBRHW_out,HLWDTH_out,HHBRB_out,HHBRL_out,HHBRM_out,HHBRS_out,HHBRC_out,HHBRT_out,HHBRG_out,HHBRGA_out,HHBRZGA_out,HHBRW_out,HHBRZ_out,HHBRH_out,HHBRA_out,HHBRAZ_out,HHBRHW_out,HHWDTH_out,HCBRB_out,HCBRL_out,HCBRM_out,HCBRBU_out,HCBRS_out,HCBRC_out,HCBRT_out,HCBRW_out,HCBRA_out,HCWDTH_out,hcbrcd_out,hcbrts_out,hcbrtd_out,hcbrwhh_out,hhbrchch_out,hlbrchch_out,abrhhaz_out,abrhawphm_out,gamt0_out,gamt1_out);
HdecayA["BR(A -> b bbar)"]=ABRB_out;
HdecayA["BR(A -> tau+ tau-)"]=ABRL_out;
HdecayA["BR(A -> mu+ mu-)"]=ABRM_out;
HdecayA["BR(A -> s sbar)"]=ABRS_out;
HdecayA["BR(A -> c cbar)"]=ABRC_out;
HdecayA["BR(A -> t tbar)"]=ABRT_out;
HdecayA["BR(A -> g g)"]=ABRG_out;
HdecayA["BR(A -> gamma gamma)"]=ABRGA_out;
HdecayA["BR(A -> Z gamma)"]=ABRZGA_out;
HdecayA["BR(A -> Z h)"]=ABRZ_out;
HdecayA["BR(A -> A Z)"]=abrhhaz_out;
HdecayA["BR(A -> A W+ H-)"]=abrhawphm_out;
HdecayA["Width"]=AWDTH_out;
HdecayHlight["BR(h -> b bbar)"]=HLBRB_out;
HdecayHlight["BR(h -> tau+ tau-)"]=HLBRL_out;
HdecayHlight["BR(h -> mu+ mu-)"]=HLBRM_out;
HdecayHlight["BR(h -> s sbar)"]=HLBRS_out;
HdecayHlight["BR(h -> c cbar)"]=HLBRC_out;
HdecayHlight["BR(h -> t tbar)"]=HLBRT_out;
HdecayHlight["BR(h -> g g)"]=HLBRG_out;
HdecayHlight["BR(h -> gamma gamma)"]=HLBRGA_out;
HdecayHlight["BR(h -> Z gamma)"]=HLBRZGA_out;
HdecayHlight["BR(h -> W+ W-)"]=HLBRW_out;
HdecayHlight["BR(h -> Z Z)"]=HLBRZ_out;
HdecayHlight["BR(h -> A A)"]=HLBRA_out;
HdecayHlight["BR(h -> Z A)"]=HLBRAZ_out;
HdecayHlight["BR(h -> H+ H-)"]=hlbrchch_out;
HdecayHlight["BR(h -> W+ H-)+BR(h -> W- H+)"]=HLBRHW_out;
HdecayHlight["Width"]=HLWDTH_out;
HdecayHheavy["BR(H -> b bbar)"]=HHBRB_out;
HdecayHheavy["BR(H -> tau+ tau-)"]=HHBRL_out;
HdecayHheavy["BR(H -> mu+ mu-)"]=HHBRM_out;
HdecayHheavy["BR(H -> s sbar)"]=HHBRS_out;
HdecayHheavy["BR(H -> c cbar)"]=HHBRC_out;
HdecayHheavy["BR(H -> t tbar)"]=HHBRT_out;
HdecayHheavy["BR(H -> g g)"]=HHBRG_out;
HdecayHheavy["BR(H -> gamma gamma)"]=HHBRGA_out;
HdecayHheavy["BR(H -> Z gamma)"]=HHBRZGA_out;
HdecayHheavy["BR(H -> W+ W-)"]=HHBRW_out;
HdecayHheavy["BR(H -> Z Z)"]=HHBRZ_out;
HdecayHheavy["BR(H -> h h)"]=HHBRH_out;
HdecayHheavy["BR(H -> A A)"]=HHBRA_out;
HdecayHheavy["BR(H -> Z A)"]=HHBRAZ_out;
HdecayHheavy["BR(H -> W+ H-)+BR(H -> W- H+)"]=HHBRHW_out;
HdecayHheavy["BR(H -> H+ H-)"]=hhbrchch_out;
HdecayHheavy["Width"]=HHWDTH_out;
HdecayHcharged["BR(H+ -> c bbar)"]=HCBRB_out;
HdecayHcharged["BR(H+ -> tau+ nu_tau)"]=HCBRL_out;
HdecayHcharged["BR(H+ -> mu+ nu_mu)"]=HCBRM_out;
HdecayHcharged["BR(H+ -> u bbar)"]=HCBRBU_out;
HdecayHcharged["BR(H+ -> u sbar)"]=HCBRS_out;
HdecayHcharged["BR(H+ -> c dbar)"]=hcbrcd_out;
HdecayHcharged["BR(H+ -> c sbar)"]=HCBRC_out;
HdecayHcharged["BR(H+ -> t bbar)"]=HCBRT_out;
HdecayHcharged["BR(H+ -> t sbar)"]=hcbrts_out;
HdecayHcharged["BR(H+ -> t dbar)"]=hcbrtd_out;
HdecayHcharged["BR(H+ -> W+ h)"]=HCBRW_out;
HdecayHcharged["BR(H+ -> W+ H)"]=hcbrwhh_out;
HdecayHcharged["BR(H+ -> W+ A)"]=HCBRA_out;
HdecayHcharged["Width"]=HCWDTH_out;
HdecayTop["BR(t -> b W+)"]=gamt0_out/gamt1_out;
HdecayTop["BR(t -> b H+)"]=(gamt1_out-gamt0_out)/gamt1_out;
HdecayTop["Width"]=gamt1_out;
}
void HdecayCalcSM(std::map<std::string,double> & HdecayHsm){
double SMBRB_out,SMBRL_out,SMBRM_out,SMBRS_out,SMBRC_out,SMBRT_out,SMBRG_out,SMBRGA_out,SMBRZGA_out,SMBRW_out,SMBRZ_out,SMWDTH_out;
hdecaycallsm_(SMBRB_out,SMBRL_out,SMBRM_out,SMBRS_out,SMBRC_out,SMBRT_out,SMBRG_out,SMBRGA_out,SMBRZGA_out,SMBRW_out,SMBRZ_out,SMWDTH_out);
HdecayHsm["BR(Hsm -> b bbar)"]=SMBRB_out;
HdecayHsm["BR(Hsm -> tau+ tau-)"]=SMBRL_out;
HdecayHsm["BR(Hsm -> mu+ mu-)"]=SMBRM_out;
HdecayHsm["BR(Hsm -> s sbar)"]=SMBRS_out;
HdecayHsm["BR(Hsm -> c cbar)"]=SMBRC_out;
HdecayHsm["BR(Hsm -> t tbar)"]=SMBRT_out;
HdecayHsm["BR(Hsm -> g g)"]=SMBRG_out;
HdecayHsm["BR(Hsm -> gamma gamma)"]=SMBRGA_out;
HdecayHsm["BR(Hsm -> Z gamma)"]=SMBRZGA_out;
HdecayHsm["BR(Hsm -> W+ W-)"]=SMBRW_out;
HdecayHsm["BR(Hsm -> Z Z)"]=SMBRZ_out;
HdecayHsm["Width"]=SMWDTH_out;
}
void CreateInputFileHdecay2HDM(int Type, double tanbeta,double alpha, double mHlight,double mHheavy, double mA, double mCharged,double m12sq){
ofstream outfile ; // ofstream output file object to write file;
outfile.open("hdecay.in"); //open file named "tempsuperiso.lha" in the current directory
outfile<<"SLHAIN = 0"<<endl;
outfile<<"SLHAOUT = 0"<<endl;
outfile<<"COUPVAR = 0"<<endl;
outfile<<"HIGGS = 5"<<endl;
outfile<<"SM4 = 0"<<endl;
outfile<<"FERMPHOB = 0"<<endl;
outfile<<"2HDM = 1"<<endl;
outfile<<"MODEL = 1"<<endl;
outfile<<fixed<<setprecision(3);
int exponent=log10(abs(tanbeta));
outfile<<"TGBET = "<<tanbeta/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"MABEG = 125.D0"<<endl;
outfile<<"MAEND = 1000.D0"<<endl;
outfile<<"NMA = 1"<<endl;
outfile<<"********************* hMSSM (MODEL = 10) *********************************"<<endl;
outfile<<"MHL = 125.D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"ALS(MZ) = 0.119D0"<<endl;
outfile<<"MSBAR(2) = 0.100D0"<<endl;
outfile<<"MC = 1.42D0"<<endl;
outfile<<"MB = 4.75D0"<<endl;
outfile<<"MT = 172.5D0"<<endl;
outfile<<"MTAU = 1.77684D0"<<endl;
outfile<<"MMUON = 0.105658367D0"<<endl;
outfile<<"1/ALPHA = 137.0359997D0"<<endl;
outfile<<"GF = 1.16637D-5"<<endl;
outfile<<"GAMW = 2.08856D0"<<endl;
outfile<<"GAMZ = 2.49581D0"<<endl;
outfile<<"MZ = 91.15349D0"<<endl;
outfile<<"MW = 80.36951D0"<<endl;
outfile<<"VTB = 0.9991D0"<<endl;
outfile<<"VTS = 0.0404D0"<<endl;
outfile<<"VTD = 0.00867D0"<<endl;
outfile<<"VCB = 0.0412D0"<<endl;
outfile<<"VCS = 0.97344D0"<<endl;
outfile<<"VCD = 0.22520D0"<<endl;
outfile<<"VUB = 0.00351D0"<<endl;
outfile<<"VUS = 0.22534D0"<<endl;
outfile<<"VUD = 0.97427D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<" SCENARIO FOR ELW. CORRECTIONS TO H -> GG (EVERYTHING IN GEV):"<<endl;
outfile<<" GG_ELW = 1: MTP = 500 MBP = 450 MNUP = 375 MEP = 450"<<endl;
outfile<<" GG_ELW = 2: MBP = MNUP = MEP = 600 MTP = MBP+50*(1+LOG(M_H/115)/5)"<<endl<<endl;
outfile<<"GG_ELW = 1"<<endl;
outfile<<"MTP = 500.D0"<<endl;
outfile<<"MBP = 450.D0"<<endl;
outfile<<"MNUP = 375.D0"<<endl;
outfile<<"MEP = 450.D0"<<endl;
outfile<<"************************** 2 Higgs Doublet Model *************************"<<endl;
outfile<<" TYPE: 1 (I), 2 (II), 3 (lepton-specific), 4 (flipped)"<<endl;
outfile<<" PARAM: 1 (masses), 2 (lambda_i)"<<endl<<endl;
outfile<<"PARAM = 1"<<endl;
outfile<<"TYPE = "<<fixed<<Type<<setprecision(3)<<endl;
outfile<<"********************"<<endl;
exponent=log10(abs(tanbeta));
outfile<<"TGBET2HDM= "<<tanbeta/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(abs(m12sq));
outfile<<"M_12^2 = "<<m12sq/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"******************** PARAM=1:"<<endl;
exponent=log10(abs(alpha));
outfile<<"ALPHA_H = "<<alpha/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mHlight);
outfile<<"MHL = "<<mHlight/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mHheavy);
outfile<<"MHH = "<<mHheavy/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mA);
outfile<<"MHA = "<<mA/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mCharged);
outfile<<"MH+- = "<<mCharged/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"******************** PARAM=2:"<<endl;
outfile<<"LAMBDA1 = 2.6885665050462264D0"<<endl;
outfile<<"LAMBDA2 = 0.000156876030254505681D0"<<endl;
outfile<<"LAMBDA3 = 0.46295674052962260D0"<<endl;
outfile<<"LAMBDA4 = 0.96605498373771792D0"<<endl;
outfile<<"LAMBDA5 = -0.88138084173680198D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"SUSYSCALE= 1000.D0"<<endl;
outfile<<"MU = 1000.D0"<<endl;
outfile<<"M2 = 1000.D0"<<endl;
outfile<<"MGLUINO = 1000.D0"<<endl;
outfile<<"MSL1 = 1000.D0"<<endl;
outfile<<"MER1 = 1000.D0"<<endl;
outfile<<"MQL1 = 1000.D0"<<endl;
outfile<<"MUR1 = 1000.D0"<<endl;
outfile<<"MDR1 = 1000.D0"<<endl;
outfile<<"MSL = 1000.D0"<<endl;
outfile<<"MER = 1000.D0"<<endl;
outfile<<"MSQ = 1000.D0"<<endl;
outfile<<"MUR = 1000.D0"<<endl;
outfile<<"MDR = 1000.D0"<<endl;
outfile<<"AL = 1000.D0"<<endl;
outfile<<"AU = 1000.D0"<<endl;
outfile<<"AD = 1000.D0"<<endl;
outfile<<"NNLO (M) = 0"<<endl;
outfile<<"ON-SHELL = 0"<<endl;
outfile<<"ON-SH-WZ = 0"<<endl;
outfile<<"IPOLE = 0"<<endl;
outfile<<"OFF-SUSY = 0"<<endl;
outfile<<"INDIDEC = 0"<<endl;
outfile<<"NF-GG = 5"<<endl;
outfile<<"IGOLD = 0"<<endl;
outfile<<"MPLANCK = 2.4D18"<<endl;
outfile<<"MGOLD = 1.D-13"<<endl;
outfile<<"******************* VARIATION OF HIGGS COUPLINGS *************************"<<endl;
outfile<<"ELWK = 0"<<endl;
outfile<<"CW = 1.D0"<<endl;
outfile<<"CZ = 1.D0"<<endl;
outfile<<"Ctau = 1.D0"<<endl;
outfile<<"Cmu = 1.D0"<<endl;
outfile<<"Ct = 1.D0"<<endl;
outfile<<"Cb = 1.D0"<<endl;
outfile<<"Cc = 1.D0"<<endl;
outfile<<"Cs = 1.D0"<<endl;
outfile<<"Cgaga = 0.D0"<<endl;
outfile<<"Cgg = 0.D0"<<endl;
outfile<<"CZga = 0.D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<"Ctp = 1.D0"<<endl;
outfile<<"Cbp = 1.D0"<<endl;
outfile<<"Cnup = 1.D0"<<endl;
outfile<<"Cep = 1.D0"<<endl;
#ifdef SHDECAYMODE
//singlet stuff
outfile<<"********************** real or complex singlet Model *********************"<<endl;
outfile<<"Model: 1 - real broken phase, 2 - real dark matter phase"<<endl;
outfile<<" 3 - complex broken phase, 4 - complex dark matter phase"<<endl;
outfile<<"icxSM = 4"<<endl;
outfile<<"*** real singlet broken phase ***"<<endl;
outfile<<"alph1 = -0.118574"<<endl;
outfile<<"mH1 = 125.1D0"<<endl;
outfile<<"mH2 = 306.361D0"<<endl;
outfile<<"vs = 293.222D0"<<endl;
outfile<<"*** real singlet dark matter phase ***"<<endl;
outfile<<"mH1 = 125.1D0"<<endl;
outfile<<"mD = 48.0215D0"<<endl;
outfile<<"m2s = -463128.D0"<<endl;
outfile<<"lambdas = 3.56328D0"<<endl;
outfile<<"*** complex singlet broken phase ***"<<endl;
outfile<<"alph1 = 0.160424D0"<<endl;
outfile<<"alph2 = -0.362128D0"<<endl;
outfile<<"alph3 = -0.552533D0"<<endl;
outfile<<"m1 = 125.518D0"<<endl;
outfile<<"m3 = 500.705D0"<<endl;
outfile<<"vs = 510.922D0"<<endl;
outfile<<"*** complex singlet dark matter phase ***"<<endl;
outfile<<"alph1 = -0.317120D0"<<endl;
outfile<<"m1 = 125.3D0"<<endl;
outfile<<"m2 = 400.D0"<<endl;
outfile<<"m3 = 731.205D0"<<endl;
outfile<<"vs = 522.181D0"<<endl;
outfile<<"a1 = -3.12115D07"<<endl;
#endif
outfile.close();
}
void CreateInputFileHdecaySM(double mH){
ofstream outfile ; // ofstream output file object to write file;
outfile.open("hdecay.in"); //open file named "tempsuperiso.lha" in the current directory
outfile<<"SLHAIN = 0"<<endl;
outfile<<"SLHAOUT = 0"<<endl;
outfile<<"COUPVAR = 0"<<endl;
outfile<<"HIGGS = 0"<<endl;
outfile<<"SM4 = 0"<<endl;
outfile<<"FERMPHOB = 0"<<endl;
outfile<<"2HDM = 0"<<endl;
outfile<<"MODEL = 1"<<endl;
outfile<<fixed<<setprecision(3);
outfile<<"TGBET = 30.D0"<<endl;
int exponent=log10(mH);
outfile<<"MABEG = "<<mH/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"MAEND = "<<mH/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"NMA = 1"<<endl;
outfile<<"********************* hMSSM (MODEL = 10) *********************************"<<endl;
outfile<<"MHL = 125.D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"ALS(MZ) = 0.119D0"<<endl;
outfile<<"MSBAR(2) = 0.100D0"<<endl;
outfile<<"MC = 1.42D0"<<endl;
outfile<<"MB = 4.75D0"<<endl;
outfile<<"MT = 172.5D0"<<endl;
outfile<<"MTAU = 1.77684D0"<<endl;
outfile<<"MMUON = 0.105658367D0"<<endl;
outfile<<"1/ALPHA = 137.0359997D0"<<endl;
outfile<<"GF = 1.16637D-5"<<endl;
outfile<<"GAMW = 2.08856D0"<<endl;
outfile<<"GAMZ = 2.49581D0"<<endl;
outfile<<"MZ = 91.15349D0"<<endl;
outfile<<"MW = 80.36951D0"<<endl;
outfile<<"VTB = 0.9991D0"<<endl;
outfile<<"VTS = 0.0404D0"<<endl;
outfile<<"VTD = 0.00867D0"<<endl;
outfile<<"VCB = 0.0412D0"<<endl;
outfile<<"VCS = 0.97344D0"<<endl;
outfile<<"VCD = 0.22520D0"<<endl;
outfile<<"VUB = 0.00351D0"<<endl;
outfile<<"VUS = 0.22534D0"<<endl;
outfile<<"VUD = 0.97427D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<" SCENARIO FOR ELW. CORRECTIONS TO H -> GG (EVERYTHING IN GEV):"<<endl;
outfile<<" GG_ELW = 1: MTP = 500 MBP = 450 MNUP = 375 MEP = 450"<<endl;
outfile<<" GG_ELW = 2: MBP = MNUP = MEP = 600 MTP = MBP+50*(1+LOG(M_H/115)/5)"<<endl<<endl;
outfile<<"GG_ELW = 1"<<endl;
outfile<<"MTP = 500.D0"<<endl;
outfile<<"MBP = 450.D0"<<endl;
outfile<<"MNUP = 375.D0"<<endl;
outfile<<"MEP = 450.D0"<<endl;
outfile<<"************************** 2 Higgs Doublet Model *************************"<<endl;
outfile<<" TYPE: 1 (I), 2 (II), 3 (lepton-specific), 4 (flipped)"<<endl;
outfile<<" PARAM: 1 (masses), 2 (lambda_i)"<<endl<<endl;
outfile<<"PARAM = 1"<<endl;
outfile<<"TYPE = 1"<<endl;
outfile<<"********************"<<endl;
outfile<<"TGBET2HDM= 2.275D0"<<endl;
outfile<<"M_12^2 = 3.045D3"<<endl;
outfile<<"******************** PARAM=1:"<<endl;
outfile<<"ALPHA_H = -1.089D0"<<endl;
exponent=log10(mH);
outfile<<"MHL = "<<mH/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"MHH = 5.517D2"<<endl;
outfile<<"MHA = 2.915D2"<<endl;
outfile<<"MH+- = 4.682D2"<<endl;
outfile<<"******************** PARAM=2:"<<endl;
outfile<<"LAMBDA1 = 2.6885665050462264D0"<<endl;
outfile<<"LAMBDA2 = 0.000156876030254505681D0"<<endl;
outfile<<"LAMBDA3 = 0.46295674052962260D0"<<endl;
outfile<<"LAMBDA4 = 0.96605498373771792D0"<<endl;
outfile<<"LAMBDA5 = -0.88138084173680198D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"SUSYSCALE= 1000.D0"<<endl;
outfile<<"MU = 1000.D0"<<endl;
outfile<<"M2 = 1000.D0"<<endl;
outfile<<"MGLUINO = 1000.D0"<<endl;
outfile<<"MSL1 = 1000.D0"<<endl;
outfile<<"MER1 = 1000.D0"<<endl;
outfile<<"MQL1 = 1000.D0"<<endl;
outfile<<"MUR1 = 1000.D0"<<endl;
outfile<<"MDR1 = 1000.D0"<<endl;
outfile<<"MSL = 1000.D0"<<endl;
outfile<<"MER = 1000.D0"<<endl;
outfile<<"MSQ = 1000.D0"<<endl;
outfile<<"MUR = 1000.D0"<<endl;
outfile<<"MDR = 1000.D0"<<endl;
outfile<<"AL = 1000.D0"<<endl;
outfile<<"AU = 1000.D0"<<endl;
outfile<<"AD = 1000.D0"<<endl;
outfile<<"NNLO (M) = 0"<<endl;
outfile<<"ON-SHELL = 0"<<endl;
outfile<<"ON-SH-WZ = 0"<<endl;
outfile<<"IPOLE = 0"<<endl;
outfile<<"OFF-SUSY = 0"<<endl;
outfile<<"INDIDEC = 0"<<endl;
outfile<<"NF-GG = 5"<<endl;
outfile<<"IGOLD = 0"<<endl;
outfile<<"MPLANCK = 2.4D18"<<endl;
outfile<<"MGOLD = 1.D-13"<<endl;
outfile<<"******************* VARIATION OF HIGGS COUPLINGS *************************"<<endl;
outfile<<"ELWK = 0"<<endl;
outfile<<"CW = 1.D0"<<endl;
outfile<<"CZ = 1.D0"<<endl;
outfile<<"Ctau = 1.D0"<<endl;
outfile<<"Cmu = 1.D0"<<endl;
outfile<<"Ct = 1.D0"<<endl;
outfile<<"Cb = 1.D0"<<endl;
outfile<<"Cc = 1.D0"<<endl;
outfile<<"Cs = 1.D0"<<endl;
outfile<<"Cgaga = 0.D0"<<endl;
outfile<<"Cgg = 0.D0"<<endl;
outfile<<"CZga = 0.D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<"Ctp = 1.D0"<<endl;
outfile<<"Cbp = 1.D0"<<endl;
outfile<<"Cnup = 1.D0"<<endl;
outfile<<"Cep = 1.D0"<<endl;
#ifdef SHDECAYMODE
//singlet stuff
outfile<<"********************** real or complex singlet Model *********************"<<endl;
outfile<<"Model: 1 - real broken phase, 2 - real dark matter phase"<<endl;
outfile<<" 3 - complex broken phase, 4 - complex dark matter phase"<<endl;
outfile<<"icxSM = 4"<<endl;
outfile<<"*** real singlet broken phase ***"<<endl;
outfile<<"alph1 = -0.118574"<<endl;
outfile<<"mH1 = 125.1D0"<<endl;
outfile<<"mH2 = 306.361D0"<<endl;
outfile<<"vs = 293.222D0"<<endl;
outfile<<"*** real singlet dark matter phase ***"<<endl;
outfile<<"mH1 = 125.1D0"<<endl;
outfile<<"mD = 48.0215D0"<<endl;
outfile<<"m2s = -463128.D0"<<endl;
outfile<<"lambdas = 3.56328D0"<<endl;
outfile<<"*** complex singlet broken phase ***"<<endl;
outfile<<"alph1 = 0.160424D0"<<endl;
outfile<<"alph2 = -0.362128D0"<<endl;
outfile<<"alph3 = -0.552533D0"<<endl;
outfile<<"m1 = 125.518D0"<<endl;
outfile<<"m3 = 500.705D0"<<endl;
outfile<<"vs = 510.922D0"<<endl;
outfile<<"*** complex singlet dark matter phase ***"<<endl;
outfile<<"alph1 = -0.317120D0"<<endl;
outfile<<"m1 = 125.3D0"<<endl;
outfile<<"m2 = 400.D0"<<endl;
outfile<<"m3 = 731.205D0"<<endl;
outfile<<"vs = 522.181D0"<<endl;
outfile<<"a1 = -3.12115D07"<<endl;
#endif
outfile.close();
}
#ifdef SHDECAYMODE
void CreateInputFileHdecayRsinglet(double mH1,double mH2,double alph1,double vs,double m2s,double lambdas,int type){
ofstream outfile ; // ofstream output file object to write file;
outfile.open("shdecay.in"); //open file named "tempsuperiso.lha" in the current directory
clog<<"CALLED"<<endl;
outfile<<"SLHAIN = 0"<<endl;
outfile<<"SLHAOUT = 0"<<endl;
outfile<<"COUPVAR = 0"<<endl;
outfile<<"HIGGS = 0"<<endl;
outfile<<"SM4 = 0"<<endl;
outfile<<"FERMPHOB = 0"<<endl;
outfile<<"2HDM = 0"<<endl;
outfile<<"MODEL = 1"<<endl;
outfile<<fixed<<setprecision(3);
outfile<<"TGBET = 30.D0"<<endl;
outfile<<"MABEG = 125.D0"<<endl;
outfile<<"MAEND = 1000.D0"<<endl;
outfile<<"NMA = 1"<<endl;
outfile<<"********************* hMSSM (MODEL = 10) *********************************"<<endl;
outfile<<"MHL = 125.D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"ALS(MZ) = 0.119D0"<<endl;
outfile<<"MSBAR(2) = 0.100D0"<<endl;
outfile<<"MC = 1.42D0"<<endl;
outfile<<"MB = 4.75D0"<<endl;
outfile<<"MT = 172.5D0"<<endl;
outfile<<"MTAU = 1.77684D0"<<endl;
outfile<<"MMUON = 0.105658367D0"<<endl;
outfile<<"1/ALPHA = 137.0359997D0"<<endl;
outfile<<"GF = 1.16637D-5"<<endl;
outfile<<"GAMW = 2.08856D0"<<endl;
outfile<<"GAMZ = 2.49581D0"<<endl;
outfile<<"MZ = 91.15349D0"<<endl;
outfile<<"MW = 80.36951D0"<<endl;
outfile<<"VTB = 0.9991D0"<<endl;
outfile<<"VTS = 0.0404D0"<<endl;
outfile<<"VTD = 0.00867D0"<<endl;
outfile<<"VCB = 0.0412D0"<<endl;
outfile<<"VCS = 0.97344D0"<<endl;
outfile<<"VCD = 0.22520D0"<<endl;
outfile<<"VUB = 0.00351D0"<<endl;
outfile<<"VUS = 0.22534D0"<<endl;
outfile<<"VUD = 0.97427D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<" SCENARIO FOR ELW. CORRECTIONS TO H -> GG (EVERYTHING IN GEV):"<<endl;
outfile<<" GG_ELW = 1: MTP = 500 MBP = 450 MNUP = 375 MEP = 450"<<endl;
outfile<<" GG_ELW = 2: MBP = MNUP = MEP = 600 MTP = MBP+50*(1+LOG(M_H/115)/5)"<<endl<<endl;
outfile<<"GG_ELW = 1"<<endl;
outfile<<"MTP = 500.D0"<<endl;
outfile<<"MBP = 450.D0"<<endl;
outfile<<"MNUP = 375.D0"<<endl;
outfile<<"MEP = 450.D0"<<endl;
outfile<<"************************** 2 Higgs Doublet Model *************************"<<endl;
outfile<<" TYPE: 1 (I), 2 (II), 3 (lepton-specific), 4 (flipped)"<<endl;
outfile<<" PARAM: 1 (masses), 2 (lambda_i)"<<endl<<endl;
outfile<<"PARAM = 1"<<endl;
outfile<<"TYPE = 1"<<endl;
outfile<<"********************"<<endl;
outfile<<"TGBET2HDM= 2.275D0"<<endl;
outfile<<"M_12^2 = 3.045D3"<<endl;
outfile<<"******************** PARAM=1:"<<endl;
outfile<<"ALPHA_H = -1.089D0"<<endl;
outfile<<"MHL = 5.517D2"<<endl;
outfile<<"MHH = 5.517D2"<<endl;
outfile<<"MHA = 2.915D2"<<endl;
outfile<<"MH+- = 4.682D2"<<endl;
outfile<<"******************** PARAM=2:"<<endl;
outfile<<"LAMBDA1 = 2.6885665050462264D0"<<endl;
outfile<<"LAMBDA2 = 0.000156876030254505681D0"<<endl;
outfile<<"LAMBDA3 = 0.46295674052962260D0"<<endl;
outfile<<"LAMBDA4 = 0.96605498373771792D0"<<endl;
outfile<<"LAMBDA5 = -0.88138084173680198D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"SUSYSCALE= 1000.D0"<<endl;
outfile<<"MU = 1000.D0"<<endl;
outfile<<"M2 = 1000.D0"<<endl;
outfile<<"MGLUINO = 1000.D0"<<endl;
outfile<<"MSL1 = 1000.D0"<<endl;
outfile<<"MER1 = 1000.D0"<<endl;
outfile<<"MQL1 = 1000.D0"<<endl;
outfile<<"MUR1 = 1000.D0"<<endl;
outfile<<"MDR1 = 1000.D0"<<endl;
outfile<<"MSL = 1000.D0"<<endl;
outfile<<"MER = 1000.D0"<<endl;
outfile<<"MSQ = 1000.D0"<<endl;
outfile<<"MUR = 1000.D0"<<endl;
outfile<<"MDR = 1000.D0"<<endl;
outfile<<"AL = 1000.D0"<<endl;
outfile<<"AU = 1000.D0"<<endl;
outfile<<"AD = 1000.D0"<<endl;
outfile<<"NNLO (M) = 0"<<endl;
outfile<<"ON-SHELL = 0"<<endl;
outfile<<"ON-SH-WZ = 0"<<endl;
outfile<<"IPOLE = 0"<<endl;
outfile<<"OFF-SUSY = 0"<<endl;
outfile<<"INDIDEC = 0"<<endl;
outfile<<"NF-GG = 5"<<endl;
outfile<<"IGOLD = 0"<<endl;
outfile<<"MPLANCK = 2.4D18"<<endl;
outfile<<"MGOLD = 1.D-13"<<endl;
outfile<<"******************* VARIATION OF HIGGS COUPLINGS *************************"<<endl;
outfile<<"ELWK = 0"<<endl;
outfile<<"CW = 1.D0"<<endl;
outfile<<"CZ = 1.D0"<<endl;
outfile<<"Ctau = 1.D0"<<endl;
outfile<<"Cmu = 1.D0"<<endl;
outfile<<"Ct = 1.D0"<<endl;
outfile<<"Cb = 1.D0"<<endl;
outfile<<"Cc = 1.D0"<<endl;
outfile<<"Cs = 1.D0"<<endl;
outfile<<"Cgaga = 0.D0"<<endl;
outfile<<"Cgg = 0.D0"<<endl;
outfile<<"CZga = 0.D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<"Ctp = 1.D0"<<endl;
outfile<<"Cbp = 1.D0"<<endl;
outfile<<"Cnup = 1.D0"<<endl;
outfile<<"Cep = 1.D0"<<endl;
//singlet stuff
outfile<<"********************** real or complex singlet Model *********************"<<endl;
outfile<<"Model: 1 - real broken phase, 2 - real dark matter phase"<<endl;
outfile<<" 3 - complex broken phase, 4 - complex dark matter phase"<<endl;
outfile<<"icxSM = "<<type<<endl;
outfile<<"*** real singlet broken phase ***"<<endl;
int exponent=log10(alph1);
outfile<<"alph1 = "<<alph1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH1);
outfile<<"mH1 = "<<mH1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH2);
outfile<<"mH2 = "<<mH2/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(vs);
outfile<<"vs = "<<vs/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"*** real singlet dark matter phase ***"<<endl;
exponent=log10(mH1);
outfile<<"mH1 = "<<mH1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH2);
outfile<<"mD = "<<mH2/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(m2s);
outfile<<"m2s = "<<m2s/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(lambdas);
outfile<<"lambdas = "<<lambdas/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"*** complex singlet broken phase ***"<<endl;
outfile<<"alph1 = 0.160424D0"<<endl;
outfile<<"alph2 = -0.362128D0"<<endl;
outfile<<"alph3 = -0.552533D0"<<endl;
outfile<<"m1 = 125.518D0"<<endl;
outfile<<"m3 = 500.705D0"<<endl;
outfile<<"vs = 510.922D0"<<endl;
outfile<<"*** complex singlet dark matter phase ***"<<endl;
outfile<<"alph1 = -0.317120D0"<<endl;
outfile<<"m1 = 125.3D0"<<endl;
outfile<<"m2 = 400.D0"<<endl;
outfile<<"m3 = 731.205D0"<<endl;
outfile<<"vs = 522.181D0"<<endl;
outfile<<"a1 = -3.12115D07"<<endl;
outfile.close();
}
void CreateInputFileHdecayCsinglet(double mH1,double mH2,double mH3,double alph1,double alph2,double alph3,double vs,double a1,int type){
ofstream outfile ; // ofstream output file object to write file;
outfile.open("shdecay.in"); //open file named "tempsuperiso.lha" in the current directory
outfile<<"SLHAIN = 0"<<endl;
outfile<<"SLHAOUT = 0"<<endl;
outfile<<"COUPVAR = 0"<<endl;
outfile<<"HIGGS = 0"<<endl;
outfile<<"SM4 = 0"<<endl;
outfile<<"FERMPHOB = 0"<<endl;
outfile<<"2HDM = 0"<<endl;
outfile<<"MODEL = 1"<<endl;
outfile<<fixed<<setprecision(3);
outfile<<"TGBET = 30.D0"<<endl;
outfile<<"MABEG = 125.D0"<<endl;
outfile<<"MAEND = 1000.D0"<<endl;
outfile<<"NMA = 1"<<endl;
outfile<<"********************* hMSSM (MODEL = 10) *********************************"<<endl;
outfile<<"MHL = 125.D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"ALS(MZ) = 0.119D0"<<endl;
outfile<<"MSBAR(2) = 0.100D0"<<endl;
outfile<<"MC = 1.42D0"<<endl;
outfile<<"MB = 4.75D0"<<endl;
outfile<<"MT = 172.5D0"<<endl;
outfile<<"MTAU = 1.77684D0"<<endl;
outfile<<"MMUON = 0.105658367D0"<<endl;
outfile<<"1/ALPHA = 137.0359997D0"<<endl;
outfile<<"GF = 1.16637D-5"<<endl;
outfile<<"GAMW = 2.08856D0"<<endl;
outfile<<"GAMZ = 2.49581D0"<<endl;
outfile<<"MZ = 91.15349D0"<<endl;
outfile<<"MW = 80.36951D0"<<endl;
outfile<<"VTB = 0.9991D0"<<endl;
outfile<<"VTS = 0.0404D0"<<endl;
outfile<<"VTD = 0.00867D0"<<endl;
outfile<<"VCB = 0.0412D0"<<endl;
outfile<<"VCS = 0.97344D0"<<endl;
outfile<<"VCD = 0.22520D0"<<endl;
outfile<<"VUB = 0.00351D0"<<endl;
outfile<<"VUS = 0.22534D0"<<endl;
outfile<<"VUD = 0.97427D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<" SCENARIO FOR ELW. CORRECTIONS TO H -> GG (EVERYTHING IN GEV):"<<endl;
outfile<<" GG_ELW = 1: MTP = 500 MBP = 450 MNUP = 375 MEP = 450"<<endl;
outfile<<" GG_ELW = 2: MBP = MNUP = MEP = 600 MTP = MBP+50*(1+LOG(M_H/115)/5)"<<endl<<endl;
outfile<<"GG_ELW = 1"<<endl;
outfile<<"MTP = 500.D0"<<endl;
outfile<<"MBP = 450.D0"<<endl;
outfile<<"MNUP = 375.D0"<<endl;
outfile<<"MEP = 450.D0"<<endl;
outfile<<"************************** 2 Higgs Doublet Model *************************"<<endl;
outfile<<" TYPE: 1 (I), 2 (II), 3 (lepton-specific), 4 (flipped)"<<endl;
outfile<<" PARAM: 1 (masses), 2 (lambda_i)"<<endl<<endl;
outfile<<"PARAM = 1"<<endl;
outfile<<"TYPE = 1"<<endl;
outfile<<"********************"<<endl;
outfile<<"TGBET2HDM= 2.275D0"<<endl;
outfile<<"M_12^2 = 3.045D3"<<endl;
outfile<<"******************** PARAM=1:"<<endl;
outfile<<"ALPHA_H = -1.089D0"<<endl;
outfile<<"MHL = 5.517D2"<<endl;
outfile<<"MHH = 5.517D2"<<endl;
outfile<<"MHA = 2.915D2"<<endl;
outfile<<"MH+- = 4.682D2"<<endl;
outfile<<"******************** PARAM=2:"<<endl;
outfile<<"LAMBDA1 = 2.6885665050462264D0"<<endl;
outfile<<"LAMBDA2 = 0.000156876030254505681D0"<<endl;
outfile<<"LAMBDA3 = 0.46295674052962260D0"<<endl;
outfile<<"LAMBDA4 = 0.96605498373771792D0"<<endl;
outfile<<"LAMBDA5 = -0.88138084173680198D0"<<endl;
outfile<<"**************************************************************************"<<endl;
outfile<<"SUSYSCALE= 1000.D0"<<endl;
outfile<<"MU = 1000.D0"<<endl;
outfile<<"M2 = 1000.D0"<<endl;
outfile<<"MGLUINO = 1000.D0"<<endl;
outfile<<"MSL1 = 1000.D0"<<endl;
outfile<<"MER1 = 1000.D0"<<endl;
outfile<<"MQL1 = 1000.D0"<<endl;
outfile<<"MUR1 = 1000.D0"<<endl;
outfile<<"MDR1 = 1000.D0"<<endl;
outfile<<"MSL = 1000.D0"<<endl;
outfile<<"MER = 1000.D0"<<endl;
outfile<<"MSQ = 1000.D0"<<endl;
outfile<<"MUR = 1000.D0"<<endl;
outfile<<"MDR = 1000.D0"<<endl;
outfile<<"AL = 1000.D0"<<endl;
outfile<<"AU = 1000.D0"<<endl;
outfile<<"AD = 1000.D0"<<endl;
outfile<<"NNLO (M) = 0"<<endl;
outfile<<"ON-SHELL = 0"<<endl;
outfile<<"ON-SH-WZ = 0"<<endl;
outfile<<"IPOLE = 0"<<endl;
outfile<<"OFF-SUSY = 0"<<endl;
outfile<<"INDIDEC = 0"<<endl;
outfile<<"NF-GG = 5"<<endl;
outfile<<"IGOLD = 0"<<endl;
outfile<<"MPLANCK = 2.4D18"<<endl;
outfile<<"MGOLD = 1.D-13"<<endl;
outfile<<"******************* VARIATION OF HIGGS COUPLINGS *************************"<<endl;
outfile<<"ELWK = 0"<<endl;
outfile<<"CW = 1.D0"<<endl;
outfile<<"CZ = 1.D0"<<endl;
outfile<<"Ctau = 1.D0"<<endl;
outfile<<"Cmu = 1.D0"<<endl;
outfile<<"Ct = 1.D0"<<endl;
outfile<<"Cb = 1.D0"<<endl;
outfile<<"Cc = 1.D0"<<endl;
outfile<<"Cs = 1.D0"<<endl;
outfile<<"Cgaga = 0.D0"<<endl;
outfile<<"Cgg = 0.D0"<<endl;
outfile<<"CZga = 0.D0"<<endl;
outfile<<"********************* 4TH GENERATION *************************************"<<endl;
outfile<<"Ctp = 1.D0"<<endl;
outfile<<"Cbp = 1.D0"<<endl;
outfile<<"Cnup = 1.D0"<<endl;
outfile<<"Cep = 1.D0"<<endl;
//singlet stuff
outfile<<"********************** real or complex singlet Model *********************"<<endl;
outfile<<"Model: 1 - real broken phase, 2 - real dark matter phase"<<endl;
outfile<<" 3 - complex broken phase, 4 - complex dark matter phase"<<endl;
outfile<<"icxSM = "<<type<<endl;
outfile<<"*** real singlet broken phase ***"<<endl;
outfile<<"alph1 = -0.118574"<<endl;
outfile<<"mH1 = 125.1D0"<<endl;
outfile<<"mH2 = 306.361D0"<<endl;
outfile<<"vs = 293.222D0"<<endl;
outfile<<"*** real singlet dark matter phase ***"<<endl;
outfile<<"mH1 = 125.1D0"<<endl;
outfile<<"mD = 48.0215D0"<<endl;
outfile<<"m2s = -463128.D0"<<endl;
outfile<<"lambdas = 3.56328D0"<<endl;
outfile<<"*** complex singlet broken phase ***"<<endl;
int exponent=log10(alph1);
outfile<<"alph1 = "<<alph1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(alph2);
outfile<<"alph2 = "<<alph2/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(alph3);
outfile<<"alph3 = "<<alph3/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH1);
outfile<<"m1 = "<<mH1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH2);
outfile<<"m3 = "<<mH3/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(vs);
outfile<<"vs = "<<vs/pow(10,exponent)<<"D"<<exponent<<endl;
outfile<<"*** complex singlet dark matter phase ***"<<endl;
exponent=log10(alph1);
outfile<<"alph1 = "<<alph1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH1);
outfile<<"m1 = "<<mH1/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH2);
outfile<<"m2 = "<<mH2/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(mH3);
outfile<<"m3 = "<<mH3/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(vs);
outfile<<"vs = "<<vs/pow(10,exponent)<<"D"<<exponent<<endl;
exponent=log10(a1);
outfile<<"a1 = "<<a1/pow(10,exponent)<<"D"<<exponent<<endl;
outfile.close();
}
void HdecayCalcRxSM(std::map<std::string,double> & HdecayHsing){
//NOT DONE
double SING1_BRB_out,SING1_BRL_out,SING1_BRM_out,SING1_BRS_out,SING1_BRC_out,SING1_BRT_out,SING1_BRG_out,SING1_BRGA_out,SING1_BRZGA_out,SING1_BRW_out,SING1_BRZ_out,SING1_WDTH_out,SING2_BRB_out,SING2_BRL_out,SING2_BRM_out,SING2_BRS_out,SING2_BRC_out,SING2_BRT_out,SING2_BRG_out,SING2_BRGA_out,SING2_BRZGA_out,SING2_BRW_out,SING2_BRZ_out,SING2_WDTH_out,SING3_BRB_out,SING3_BRL_out,SING3_BRM_out,SING3_BRS_out,SING3_BRC_out,SING3_BRT_out,SING3_BRG_out,SING3_BRGA_out,SING3_BRZGA_out,SING3_BRW_out,SING3_BRZ_out,SING3_WDTH_out,SING_BRH1H2H2_out,SING_BRH1H3H3_out,SING_BRH1H2H3_out,SING_BRH2H1H1_out,SING_BRH2H3H3_out,SING_BRH2H1H3_out,SING_BRH3H1H1_out,SING_BRH3H2H2_out,SING_BRH3H1H2_out;
hdecaycallsinglets_(SING1_BRB_out,SING1_BRL_out,SING1_BRM_out,SING1_BRS_out,SING1_BRC_out,SING1_BRT_out,SING1_BRG_out,SING1_BRGA_out,SING1_BRZGA_out,SING1_BRW_out,SING1_BRZ_out,SING1_WDTH_out,SING2_BRB_out,SING2_BRL_out,SING2_BRM_out,SING2_BRS_out,SING2_BRC_out,SING2_BRT_out,SING2_BRG_out,SING2_BRGA_out,SING2_BRZGA_out,SING2_BRW_out,SING2_BRZ_out,SING2_WDTH_out,SING3_BRB_out,SING3_BRL_out,SING3_BRM_out,SING3_BRS_out,SING3_BRC_out,SING3_BRT_out,SING3_BRG_out,SING3_BRGA_out,SING3_BRZGA_out,SING3_BRW_out,SING3_BRZ_out,SING3_WDTH_out,SING_BRH1H2H2_out,SING_BRH1H3H3_out,SING_BRH1H2H3_out,SING_BRH2H1H1_out,SING_BRH2H3H3_out,SING_BRH2H1H3_out,SING_BRH3H1H1_out,SING_BRH3H2H2_out,SING_BRH3H1H2_out);
HdecayHsing["BR(H1 -> b bbar)"]=SING1_BRB_out;
HdecayHsing["BR(H1 -> tau+ tau-)"]=SING1_BRL_out;
HdecayHsing["BR(H1 -> mu+ mu-)"]=SING1_BRM_out;
HdecayHsing["BR(H1 -> s sbar)"]=SING1_BRS_out;
HdecayHsing["BR(H1 -> c cbar)"]=SING1_BRC_out;
HdecayHsing["BR(H1 -> t tbar)"]=SING1_BRT_out;
HdecayHsing["BR(H1 -> g g)"]=SING1_BRG_out;
HdecayHsing["BR(H1 -> gamma gamma)"]=SING1_BRGA_out;
HdecayHsing["BR(H1 -> Z gamma)"]=SING1_BRZGA_out;
HdecayHsing["BR(H1 -> W+ W-)"]=SING1_BRW_out;
HdecayHsing["BR(H1 -> Z Z)"]=SING1_BRZ_out;
HdecayHsing["H1Width"]=SING1_WDTH_out;
HdecayHsing["BR(H2 -> b bbar)"]=SING2_BRB_out;
HdecayHsing["BR(H2 -> tau+ tau-)"]=SING2_BRL_out;
HdecayHsing["BR(H2 -> mu+ mu-)"]=SING2_BRM_out;
HdecayHsing["BR(H2 -> s sbar)"]=SING2_BRS_out;
HdecayHsing["BR(H2 -> c cbar)"]=SING2_BRC_out;
HdecayHsing["BR(H2 -> t tbar)"]=SING2_BRT_out;
HdecayHsing["BR(H2 -> g g)"]=SING2_BRG_out;
HdecayHsing["BR(H2 -> gamma gamma)"]=SING2_BRGA_out;
HdecayHsing["BR(H2 -> Z gamma)"]=SING2_BRZGA_out;
HdecayHsing["BR(H2 -> W+ W-)"]=SING2_BRW_out;
HdecayHsing["BR(H2 -> Z Z)"]=SING2_BRZ_out;
HdecayHsing["H2Width"]=SING2_WDTH_out;
HdecayHsing["BR(H3 -> b bbar)"]=SING3_BRB_out;
HdecayHsing["BR(H3 -> tau+ tau-)"]=SING3_BRL_out;
HdecayHsing["BR(H3 -> mu+ mu-)"]=SING3_BRM_out;
HdecayHsing["BR(H3 -> s sbar)"]=SING3_BRS_out;
HdecayHsing["BR(H3 -> c cbar)"]=SING3_BRC_out;
HdecayHsing["BR(H3 -> t tbar)"]=SING3_BRT_out;
HdecayHsing["BR(H3 -> g g)"]=SING3_BRG_out;
HdecayHsing["BR(H3 -> gamma gamma)"]=SING3_BRGA_out;
HdecayHsing["BR(H3 -> Z gamma)"]=SING3_BRZGA_out;
HdecayHsing["BR(H3 -> W+ W-)"]=SING3_BRW_out;
HdecayHsing["BR(H3 -> Z Z)"]=SING3_BRZ_out;
HdecayHsing["H3Width"]=SING3_WDTH_out;
HdecayHsing["BR(H1 -> H2 H2)"]=SING_BRH1H2H2_out;
HdecayHsing["BR(H1 -> H3 H3)"]=SING_BRH1H3H3_out;
HdecayHsing["BR(H1 -> H2 H3)"]=SING_BRH1H2H3_out;
HdecayHsing["BR(H2 -> H1 H1)"]=SING_BRH2H1H1_out;
HdecayHsing["BR(H2 -> H3 H3)"]=SING_BRH2H3H3_out;
HdecayHsing["BR(H2 -> H1 H3)"]=SING_BRH2H1H3_out;
HdecayHsing["BR(H3 -> H1 H1)"]=SING_BRH3H1H1_out;
HdecayHsing["BR(H3 -> H2 H2)"]=SING_BRH3H2H2_out;
HdecayHsing["BR(H3 -> H1 H2)"]=SING_BRH3H1H2_out;
}
#endif

File Metadata

Mime Type
text/x-c
Expires
Sat, Dec 21, 12:22 PM (1 d, 17 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4022752
Default Alt Text
HdecayWrapCppFuncs.cpp (36 KB)

Event Timeline