// obj declaration:
//  gen_obj the_name(arg1, "arg2", "arg3");
//  with: 
//  arg 1 values, true->parameter, false->observable of interest, not a parameter
//  arg 2 is the working label
//  arg 3 is the latex label
//  the_name.declare(themodel.logfile, arg1, arg2, arg3, arg4);
//  themodel.logfile is the log file to which declaration is reported
//  arg1 is the number of bins to be considered
//  arg2 is the min of the range
//  arg3 is the max of the range
//  arg4 is the step in case of a parameter  

int NBinsL = 200, NBinsM = 80, NBinsS = 25;
obs CHI2_min = 0, CHI2_max = 15.0;  
obs DCHI2_min = 0, DCHI2_max = 5.0;  

obs log_BRmin = -4.0;

// Scalar sector
gen_obj log10_tanbeta_obj(true, "log10_tanbeta", "$\\\\log_{10}\\\\tan\\\\beta$");
log10_tanbeta_obj.declare(themodel.logfile, NBinsL, -0.50, 2.5, 5.0e-4);

gen_obj rot_ba_obj(true, "rot_ba", "$\\\\beta-\\\\alpha$");
rot_ba_obj.declare(themodel.logfile, NBinsM, 1.55, 1.59, 1.0e-5);

gen_obj R11_obj(false, "R11", "$R_{11}$");
R11_obj.declare(themodel.logfile, NBinsM, 0.75, 1.0);

gen_obj R21_obj(false, "R21", "$R_{21}(\\\\text{log}_{10},[-1;\\\\pm 10^{-4};+1])$");
R21_obj.declare(themodel.logfile, NBinsL, -4.0, 4.0);

gen_obj m_ChH_obj(true, "m_ChH", "$m_{H^\\\\pm}$ GeV");
m_ChH_obj.declare(themodel.logfile, NBinsL, 0.0, 2500.0, 2.0e-01);

gen_obj m_H_obj(true, "m_H", "$m_{H}$ GeV");
m_H_obj.declare(themodel.logfile, NBinsL, 0.0, 2500.0, 2.0e-01);

gen_obj m_A_obj(true, "m_A", "$m_{A}$ GeV");
m_A_obj.declare(themodel.logfile, NBinsL, 0.0, 2500.0, 2.0e-01);

gen_obj re_mu12_obj(true, "re_mu12", "Re$(\\\\mu_{12}^2)$ GeV$^2$");
re_mu12_obj.declare(themodel.logfile, NBinsL, -2.0e6, 2.0e6, 0.20e2);

gen_obj lambda_1_obj(false, "lambda_1", "$\\\\lambda_{1}$");
lambda_1_obj.declare(themodel.logfile, NBinsM, -5.0, 15.0);

gen_obj lambda_2_obj(false, "lambda_2", "$\\\\lambda_{2}$");
lambda_2_obj.declare(themodel.logfile, NBinsM, -5.0, 15.0);

gen_obj lambda_3_obj(false, "lambda_3", "$\\\\lambda_{3}$");
lambda_3_obj.declare(themodel.logfile, NBinsM, -15.0, 15.0);

gen_obj lambda_4_obj(false, "lambda_4", "$\\\\lambda_{4}$");
lambda_4_obj.declare(themodel.logfile, NBinsM, -15.0, 15.0);

gen_obj lambda_5_obj(false, "lambda_5", "$|\\\\lambda_{5}|$");
lambda_5_obj.declare(themodel.logfile, NBinsM, -5.0, 15.0);

gen_obj a_nle_obj(true, "a_nle", "$Re(n_{e})$");
//a_nle_obj.declare(themodel.logfile, NBinsL, 0.0, 50.0, 5.0e-2);
a_nle_obj.declare(themodel.logfile, NBinsL, 0.0, 60.0, 5.0e-2);

gen_obj b_nle_obj(true, "b_nle", "$Im(n_{e})$");
b_nle_obj.declare(themodel.logfile, NBinsL, -1.0e-1, 1.0e-1, 1.0e-5);

gen_obj a_nlm_obj(true, "a_nlm", "$Re(n_{\\\\mu})$");
//a_nlm_obj.declare(themodel.logfile, NBinsL, -100.0, 0.0, 5.0e-2);
a_nlm_obj.declare(themodel.logfile, NBinsL, -100.0, 100.0, 5.0e-2);

gen_obj b_nlm_obj(true, "b_nlm", "$Im(n_{\\\\mu})$"); 
b_nlm_obj.declare(themodel.logfile, NBinsL, -100.0, 25.0, 1.0e-2);

gen_obj a_nlt_obj(true, "a_nlt", "$Re(n_{\\\\tau})$");
a_nlt_obj.declare(themodel.logfile, NBinsL, -100.0, 100.0, 5.0e-2);

gen_obj b_nlt_obj(true, "b_nlt", "$Im(n_{\\\\tau})$"); 
b_nlt_obj.declare(themodel.logfile, NBinsL, -100.0, 100.0, 1.0e-2);

gen_obj anom_le_obj(false, "anom_le", "$\\\\delta a_{e}\\\\times 10^{13}$");
anom_le_obj.declare(themodel.logfile, NBinsL, -20.0, 0.0);

gen_obj anom_lm_obj(false, "anom_lm", "$\\\\delta a_{\\\\mu}\\\\times 10^{9}$");
anom_lm_obj.declare(themodel.logfile, NBinsL, 0.0, 6.0);

gen_obj anom_lt_obj(false, "anom_lt", "log$_{10}|\\\\delta a_{\\\\tau}|$");
anom_lt_obj.declare(themodel.logfile, NBinsL, -10.0, -3.0);

gen_obj edm_le_obj(false, "edm_le", "$|d_{e}|\\\\times 10^{29}$ e.cm");
edm_le_obj.declare(themodel.logfile, NBinsL, 0.0, 2.0);

gen_obj edm_lm_obj(false, "edm_lm", "log$_{10}|d_{\\\\mu}|$ (e.cm)");
edm_lm_obj.declare(themodel.logfile, NBinsL, -26.0, -18.0);

gen_obj edm_lt_obj(false, "edm_lt", "$log_{10}|d_{\\\\tau}|$ e.cm");
edm_lt_obj.declare(themodel.logfile, NBinsL, -26.0, -18.0);

gen_obj G_h_obj(false, "G_h", "$\\\\Gamma_h$");
G_h_obj.declare(themodel.logfile, NBinsM, 0.0, 10.0);

gen_obj G_H_obj(false, "G_H", "$\\\\log_{10}\\\\Gamma_H/m_H$");
G_H_obj.declare(themodel.logfile, NBinsL, -2.0, 1.0);

gen_obj G_A_obj(false, "G_A", "$\\\\log_{10}\\\\Gamma_A/m_A$");
G_A_obj.declare(themodel.logfile, NBinsL, -2.0, 1.0);

gen_obj G_ChH_obj(false, "G_ChH", "$\\\\log_{10}\\\\Gamma_{H^\\\\pm}/m_{H^\\\\pm}$");
G_ChH_obj.declare(themodel.logfile, NBinsL, -2.0, 1.0);

gen_obj BR_H_lele_obj(false, "BR_H_lele", "$\\\\log_{10}$BR$(H\\\\to e^+e^-)$");
BR_H_lele_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_H_lmlm_obj(false, "BR_H_lmlm", "$\\\\log_{10}$BR$(H\\\\to\\\\mu^+\\\\mu^-)$");
BR_H_lmlm_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_H_ltlt_obj(false, "BR_H_ltlt", "$\\\\log_{10}$BR$(H\\\\to\\\\tau^+\\\\tau^-)$");
BR_H_ltlt_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_H_qtqt_obj(false, "BR_H_qtqt", "$\\\\log_{10}$BR$(H\\\\to\\\\bar tt)$");
BR_H_qtqt_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_H_qbqb_obj(false, "BR_H_qbqb", "$\\\\log_{10}$BR$(H\\\\to\\\\bar bb)$");
BR_H_qbqb_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_A_lele_obj(false, "BR_A_lele", "$\\\\log_{10}$BR$(A\\\\to e^+e^-)$");
BR_A_lele_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_A_lmlm_obj(false, "BR_A_lmlm", "$\\\\log_{10}$BR$(A\\\\to\\\\mu^+\\\\mu^-)$");
BR_A_lmlm_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_A_ltlt_obj(false, "BR_A_ltlt", "$\\\\log_{10}$BR$(A\\\\to\\\\tau^+\\\\tau^-)$");
BR_A_ltlt_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_A_qtqt_obj(false, "BR_A_qtqt", "$\\\\log_{10}$BR$(A\\\\to\\\\bar tt)$");
BR_A_qtqt_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_A_qbqb_obj(false, "BR_A_qbqb", "$\\\\log_{10}$BR$(A\\\\to\\\\bar bb)$");
BR_A_qbqb_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_A_HZ_obj(false, "BR_A_HZ", "$\\\\log_{10}$BR$(A\\\\to HZ)$");
BR_A_HZ_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj pp_H_lele_xSggF_BR_obj(false, "pp_H_lele_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H\\\\to e^+e^-)$ (fb)");
pp_H_lele_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj pp_H_lmlm_xSggF_BR_obj(false, "pp_H_lmlm_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H\\\\to\\\\mu^+\\\\mu^-)$ (fb)");
pp_H_lmlm_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj pp_H_ltlt_xSggF_BR_obj(false, "pp_H_ltlt_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H\\\\to\\\\tau^+\\\\tau^-)$ (fb)");
pp_H_ltlt_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj pp_A_lele_xSggF_BR_obj(false, "pp_A_lele_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to A\\\\to e^+e^-)$ (fb)");
pp_A_lele_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj pp_A_lmlm_xSggF_BR_obj(false, "pp_A_lmlm_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to A\\\\to\\\\mu^+\\\\mu^-)$ (fb)");
pp_A_lmlm_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj pp_A_ltlt_xSggF_BR_obj(false, "pp_A_ltlt_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to A\\\\to\\\\tau^+\\\\tau^-)$ (fb)");
pp_A_ltlt_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj pp_A_HZ_xSggF_BR_obj(false, "pp_A_HZ_xSggF_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to A\\\\to HZ)$ (fb)");
pp_A_HZ_xSggF_BR_obj.declare(themodel.logfile, NBinsL, -2.0, 4.0);

gen_obj BR_ChH_lenu_obj(false, "BR_ChH_lenu", "$\\\\log_{10}$BR$(H^\\\\pm\\\\to e^\\\\pm\\\\nu)$");
BR_ChH_lenu_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_ChH_lmnu_obj(false, "BR_ChH_lmnu", "$\\\\log_{10}$BR$(H^\\\\pm\\\\to\\\\mu^\\\\pm\\\\nu)$");
BR_ChH_lmnu_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_ChH_ltnu_obj(false, "BR_ChH_ltnu", "$\\\\log_{10}$BR$(H^\\\\pm\\\\to\\\\tau^\\\\pm\\\\nu)$");
BR_ChH_ltnu_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_ChH_qtqb_obj(false, "BR_ChH_qtqb", "$\\\\log_{10}$BR$(H^\\\\pm\\\\to tb)$");
BR_ChH_qtqb_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj BR_ChH_HW_obj(false, "BR_ChH_HW", "$\\\\log_{10}$BR$(H^\\\\pm\\\\to HW^\\\\pm)$");
BR_ChH_HW_obj.declare(themodel.logfile, NBinsL, log_BRmin, 0.0);

gen_obj pp_ChHtb_ltln_xS_BR_obj(false, "pp_ChHtb_ltln_xS_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H^\\\\pm tb)(H^\\\\pm\\\\to\\\\tau^\\\\pm\\\\nu)$ (fb)");
pp_ChHtb_ltln_xS_BR_obj.declare(themodel.logfile, NBinsL, -3.0, 3.0);

gen_obj pp_ChHtb_lmln_xS_BR_obj(false, "pp_ChHtb_lmln_xS_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H^\\\\pm tb)(H^\\\\pm\\\\to\\\\mu^\\\\pm\\\\nu)$ (fb)");
pp_ChHtb_lmln_xS_BR_obj.declare(themodel.logfile, NBinsL, -3.0, 3.0);

gen_obj pp_ChHtb_leln_xS_BR_obj(false, "pp_ChHtb_leln_xS_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H^\\\\pm tb)(H^\\\\pm\\\\to e^\\\\pm\\\\nu)$ (fb)");
pp_ChHtb_leln_xS_BR_obj.declare(themodel.logfile, NBinsL, -3.0, 3.0);

gen_obj pp_ChHtb_qtqb_xS_BR_obj(false, "pp_ChHtb_qtqb_xS_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H^\\\\pm tb)(H^\\\\pm\\\\to tb)$ (fb)");
pp_ChHtb_qtqb_xS_BR_obj.declare(themodel.logfile, NBinsL, -3.0, 3.0);

gen_obj pp_ChHtb_HW_xS_BR_obj(false, "pp_ChHtb_HW_xS_BR", "$\\\\log_{10}\\\\sigma(pp\\\\to H^\\\\pm tb)(H^\\\\pm\\\\to HW^\\\\pm)$ (fb)");
pp_ChHtb_HW_xS_BR_obj.declare(themodel.logfile, NBinsL, -3.0, 3.0);

gen_obj Oblique_T_obj(false, "Oblique_T", "$\\\\Delta T$");
Oblique_T_obj.declare(themodel.logfile, NBinsM, -0.50, 0.50);

gen_obj Oblique_S_obj(false, "Oblique_S", "$\\\\Delta S$");
Oblique_S_obj.declare(themodel.logfile, NBinsM, -0.50, 0.50);

gen_obj BR_bsg_obj(false, "BR_bsg", "BR$(bs\\\\gamma)$");
BR_bsg_obj.declare(themodel.logfile, NBinsM, 2.50, 4.50);

gen_obj chi2_2HDMScalar_pert_obj(false, "chi2_2HDMScalar_pert", "$\\\\chi^2$(2HDM-pert)");
chi2_2HDMScalar_pert_obj.declare(themodel.logfile, NBinsL, 0.0, 20.0);
gen_obj chi2_yukawa_pert_obj(false, "chi2_yukawa_pert", "$\\\\chi^2$(Yukawa-pert)");
chi2_yukawa_pert_obj.declare(themodel.logfile, NBinsL, 0.0, 20.0);
gen_obj chi2_2HDMScalar_pert_unit_obj(false, "chi2_2HDMScalar_pert_unit", "$\\\\chi^2$(2HDM-pert-unit)");
chi2_2HDMScalar_pert_unit_obj.declare(themodel.logfile, NBinsL, 0.0, 20.0);
gen_obj chi2_ScalarMasses_obj(false, "chi2_ScalarMasses", "$\\\\chi^2$(ScalarMasses)");
chi2_ScalarMasses_obj.declare(themodel.logfile, NBinsS, 0.0, 20.0);
gen_obj chi2_Oblique_obj(false, "chi2_Oblique", "$\\\\chi^2$(Oblique)");
chi2_Oblique_obj.declare(themodel.logfile, NBinsS, 0.0, 20.0);
gen_obj chi2_Higgs_obj(false, "chi2_Higgs", "$\\\\chi^2$(Higgs)");
chi2_Higgs_obj.declare(themodel.logfile, NBinsS, 0.0, 20.0);
gen_obj chi2_dipoles_obj(false, "chi2_dipoles", "$\\\\chi^2$(dipoles)");
chi2_dipoles_obj.declare(themodel.logfile, NBinsS, 0.0, 20.0);
gen_obj chi2_bsgamma_obj(false, "chi2_bsgamma", "$\\\\chi^2$(bsg)");
chi2_bsgamma_obj.declare(themodel.logfile, NBinsS, 0.0, 20.0);
gen_obj chi2_tree_Charged_obj(false, "chi2_tree_Charged", "$\\\\chi^2$(tree Charged)");
chi2_tree_Charged_obj.declare(themodel.logfile, NBinsM, 0.0, 20.0);
gen_obj chi2_LEP_eemm_obj(false, "chi2_LEP_eemm", "$\\\\chi^2$(LEP$(e^+e^-\\\\to\\\\mu^+\\\\mu^-)$)");
chi2_LEP_eemm_obj.declare(themodel.logfile, NBinsM, 0.0, 20.0);
gen_obj chi2_LEP_eett_obj(false, "chi2_LEP_eett", "$\\\\chi^2$(LEP$(e^+e^-\\\\to\\\\tau^+\\\\tau^-)$)");
chi2_LEP_eett_obj.declare(themodel.logfile, NBinsM, 0.0, 20.0);
gen_obj chi2_LHC_ll_obj(false, "chi2_LHC_ll", "$\\\\chi^2$(LHC$(pp\\\\to S\\\\to\\\\ell^+\\\\ell^-)$)");
chi2_LHC_ll_obj.declare(themodel.logfile, NBinsM, 0.0, 20.0);
gen_obj chi2_LHC_ChH_obj(false, "chi2_LHC_ChH", "$\\\\chi^2$(LHC$(pp\\\\to H^\\\\pm tb)(H^\\\\to\\\\tau^\\\\pm\\\\nu,tb)$)");
chi2_LHC_ChH_obj.declare(themodel.logfile, NBinsM, 0.0, 20.0);
gen_obj chi2_mesonmixing_obj(false, "chi2_mesonmixing", "$\\\\chi^2$(Meson Mixing)");
chi2_mesonmixing_obj.declare(themodel.logfile, NBinsL, 0.0, 20.0);

gen_obj *obs_grp[] = {&log10_tanbeta_obj, &rot_ba_obj, &m_ChH_obj, &m_H_obj, &m_A_obj, &G_h_obj, &G_H_obj, &G_A_obj, &G_ChH_obj, &re_mu12_obj, &lambda_1_obj, &lambda_2_obj, &lambda_3_obj, &lambda_4_obj, &lambda_5_obj, &R11_obj, &R21_obj, &a_nle_obj, &b_nle_obj, &a_nlm_obj, &b_nlm_obj, &a_nlt_obj, &b_nlt_obj, &anom_le_obj, &anom_lm_obj, &anom_lt_obj, &edm_le_obj, &edm_lm_obj, &edm_lt_obj, &BR_H_lele_obj, &BR_H_lmlm_obj, &BR_H_ltlt_obj, &BR_H_qtqt_obj, &BR_H_qbqb_obj, &BR_A_lele_obj, &BR_A_lmlm_obj, &BR_A_ltlt_obj, &BR_A_qtqt_obj, &BR_A_qbqb_obj, &BR_A_HZ_obj, &pp_H_lele_xSggF_BR_obj, &pp_H_lmlm_xSggF_BR_obj, &pp_H_ltlt_xSggF_BR_obj, &pp_A_lele_xSggF_BR_obj, &pp_A_lmlm_xSggF_BR_obj, &pp_A_ltlt_xSggF_BR_obj, &pp_A_HZ_xSggF_BR_obj, &BR_ChH_lenu_obj, &BR_ChH_lmnu_obj, &BR_ChH_ltnu_obj, &BR_ChH_qtqb_obj, &BR_ChH_HW_obj, &pp_ChHtb_leln_xS_BR_obj, &pp_ChHtb_lmln_xS_BR_obj, &pp_ChHtb_ltln_xS_BR_obj, &pp_ChHtb_qtqb_xS_BR_obj, &pp_ChHtb_HW_xS_BR_obj, &Oblique_T_obj, &Oblique_S_obj, &BR_bsg_obj, &chi2_2HDMScalar_pert_obj, &chi2_2HDMScalar_pert_unit_obj, &chi2_ScalarMasses_obj, &chi2_Oblique_obj, &chi2_Higgs_obj, &chi2_dipoles_obj, &chi2_yukawa_pert_obj, &chi2_bsgamma_obj, &chi2_tree_Charged_obj, &chi2_LEP_eemm_obj, &chi2_LEP_eett_obj, &chi2_LHC_ll_obj, &chi2_LHC_ChH_obj, &chi2_mesonmixing_obj};

// Set the number of observables and of groups of observables
themodel.n_obs = arraysize(obs_grp);
themodel.n_grps = 6;
themodel.set_obs_pt(obs_grp, themodel.n_obs);

// For each observable individual declaration observable_object.add_to_obs_grp(n) if it belongs to the observable group n

log10_tanbeta_obj.add_to_obs_grp(1);
R21_obj.add_to_obs_grp(1);
m_ChH_obj.add_to_obs_grp(1);
m_H_obj.add_to_obs_grp(1);
m_A_obj.add_to_obs_grp(1);

a_nle_obj.add_to_obs_grp(1);
//  b_nle_obj.add_to_obs_grp(1);
a_nlm_obj.add_to_obs_grp(1);
//  b_nlm_obj.add_to_obs_grp(1);
a_nlt_obj.add_to_obs_grp(1);
//  b_nlt_obj.add_to_obs_grp(1);

//anom_le_obj.add_to_obs_grp(1);
//anom_lm_obj.add_to_obs_grp(1);
//anom_lt_obj.add_to_obs_grp(1);

//  edm_le_obj.add_to_obs_grp(1);
//  edm_lm_obj.add_to_obs_grp(1);
//  edm_lt_obj.add_to_obs_grp(1);

chi2_tree_Charged_obj.add_to_obs_grp(1);

m_H_obj.add_to_obs_grp(2);
m_A_obj.add_to_obs_grp(2);
log10_tanbeta_obj.add_to_obs_grp(2);
pp_H_lele_xSggF_BR_obj.add_to_obs_grp(2);
pp_H_lmlm_xSggF_BR_obj.add_to_obs_grp(2);
pp_H_ltlt_xSggF_BR_obj.add_to_obs_grp(2);
pp_A_lele_xSggF_BR_obj.add_to_obs_grp(2);
pp_A_lmlm_xSggF_BR_obj.add_to_obs_grp(2);
pp_A_ltlt_xSggF_BR_obj.add_to_obs_grp(2);
pp_A_HZ_xSggF_BR_obj.add_to_obs_grp(2);

m_ChH_obj.add_to_obs_grp(3);
log10_tanbeta_obj.add_to_obs_grp(3);
pp_ChHtb_leln_xS_BR_obj.add_to_obs_grp(3);
pp_ChHtb_lmln_xS_BR_obj.add_to_obs_grp(3);
pp_ChHtb_ltln_xS_BR_obj.add_to_obs_grp(3);
pp_ChHtb_qtqb_xS_BR_obj.add_to_obs_grp(3);
pp_ChHtb_HW_xS_BR_obj.add_to_obs_grp(3);

// chi2_tree_Charged_obj.add_to_obs_grp(2);
// chi2_LEP_eemm_obj.add_to_obs_grp(2);
// chi2_LEP_eett_obj.add_to_obs_grp(2);
// chi2_LHC_ll_obj.add_to_obs_grp(2);

// G_ChH_obj.add_to_obs_grp(3);
// G_H_obj.add_to_obs_grp(3);
// G_A_obj.add_to_obs_grp(3);

m_ChH_obj.add_to_obs_grp(4);
BR_ChH_lenu_obj.add_to_obs_grp(4);
BR_ChH_lmnu_obj.add_to_obs_grp(4);
BR_ChH_ltnu_obj.add_to_obs_grp(4);
BR_ChH_qtqb_obj.add_to_obs_grp(4);
BR_ChH_HW_obj.add_to_obs_grp(4);
G_ChH_obj.add_to_obs_grp(4);

m_H_obj.add_to_obs_grp(5);
BR_H_lele_obj.add_to_obs_grp(5);
BR_H_lmlm_obj.add_to_obs_grp(5);
BR_H_ltlt_obj.add_to_obs_grp(5);
BR_H_qtqt_obj.add_to_obs_grp(5);
BR_H_qbqb_obj.add_to_obs_grp(5);
G_H_obj.add_to_obs_grp(5);

m_A_obj.add_to_obs_grp(6);
BR_A_lele_obj.add_to_obs_grp(6);
BR_A_lmlm_obj.add_to_obs_grp(6);
BR_A_ltlt_obj.add_to_obs_grp(6);
BR_A_qtqt_obj.add_to_obs_grp(6);
BR_A_qbqb_obj.add_to_obs_grp(6);
BR_A_HZ_obj.add_to_obs_grp(6);
G_A_obj.add_to_obs_grp(6);

gen_obj *par_grp[] = {&log10_tanbeta_obj, &rot_ba_obj, &m_ChH_obj, &m_H_obj, &m_A_obj, &re_mu12_obj, &a_nle_obj, &b_nle_obj, &a_nlm_obj, &b_nlm_obj, &a_nlt_obj, &b_nlt_obj};
themodel.n_pars = arraysize(par_grp);
themodel.set_par_pt(par_grp, themodel.n_pars);

/*
  chi2_obj *chi2_grp[] = {&chi2_Smps};
  chi2_obj *dchi2_grp[] = {&dchi2_Smps};
  themodel.n_chi2 = arraysize(chi2_grp);
  themodel.set_chi2_pt(chi2_grp, themodel.n_chi2);
  themodel.set_dchi2_pt(dchi2_grp, themodel.n_chi2);
*/