main_thdm.cpp

Go to the documentation of this file.

#include <iostream>
 
#include "MemoryManager.h"
#include "Parameters.h"
#include "WilsonInterface.h"
#include "BlockProxy.h"
 
int main() {
    auto hyp = HyperisoMaster();
    HyperisoConfig config_hyp;
    // config_hyp.flags[ExternalFlag::USE_MARTY] = false;
    config_hyp.model = Model::THDM;
    config_hyp.mty_model_name = "THDM";
    config_hyp.mty_model_path = std::string(project_tp_root.data()) + "MARTY/src/MARTY/src/marty/models/thdm.h";
 
    hyp.init("lha/testinput_thdm.lha", config_hyp); // Initialize program manager with LHA file
 
    auto wi = WilsonInterface(); // Initialize interface and build the required groups
    std::cout << "here " << std::endl;
    WilsonBuildConfig config({GroupMapper::to_id(WGroup::BScalar), GroupMapper::to_id(WGroup::B), GroupMapper::to_id(WGroup::BPrime)}, 81, 42, QCDOrder::LO);
    std::cout << "here " << std::endl;
    wi.build(config);
    std::cout << "here " << std::endl;
    // BlockProxy().log_block(ParameterType::WILSON, "BCoefficients_B_SCALE_STANDARD");
    BlockProxy().log_all_blocks(ParameterType::WILSON);
    BlockProxy().log_all_blocks(ParameterType::SM);
    BlockProxy().log_all_blocks(ParameterType::BSM);
    for (auto elem : {WCoef::CQ1_MU, WCoef::CQ2_MU}) {
        LOG_INFO(WCoefMapper::str(elem), "(MW) at LO =", wi.getM(WGroup::BScalar, elem, QCDOrder::LO, ContributionType::BSM));
        LOG_INFO(WCoefMapper::str(elem), "(MW) at NLO =", wi.getM(WGroup::BScalar, elem, QCDOrder::NLO, ContributionType::BSM));
        LOG_INFO("\n");
    }
 
    for (auto elem : {WCoef::C1, WCoef::C2, WCoef::C3, WCoef::C4, WCoef::C5, WCoef::C6, WCoef::C7, WCoef::C8, WCoef::C9, WCoef::C10}) {
        LOG_INFO(WCoefMapper::str(elem), "(MW) at LO =", wi.getM(WGroup::B, elem, QCDOrder::LO, ContributionType::BSM));
        LOG_INFO(WCoefMapper::str(elem), "(MW) at NLO =", wi.getM(WGroup::B, elem, QCDOrder::NLO, ContributionType::BSM));
        LOG_INFO(WCoefMapper::str(elem), "(MW) at NNLO =", wi.getM(WGroup::B, elem, QCDOrder::NNLO, ContributionType::BSM));
        LOG_INFO("\n");
    }
 
    for (auto elem : {WCoef::CP1, WCoef::CP2, WCoef::CP3, WCoef::CP4, WCoef::CP5, WCoef::CP6, WCoef::CP7, WCoef::CP8, WCoef::CP9, WCoef::CP10, WCoef::CPQ1_MU, WCoef::CPQ2_MU}) {
        LOG_INFO(WCoefMapper::str(elem), "(MW) at LO =", wi.getM(WGroup::BPrime, elem, QCDOrder::LO, ContributionType::BSM));
        LOG_INFO(WCoefMapper::str(elem), "(MW) at NLO =", wi.getM(WGroup::BPrime, elem, QCDOrder::NLO, ContributionType::BSM));
        LOG_INFO("\n");
    }
 
    LOG_INFO("Parameters created");
    // LOG_INFO("C9(mu_h) at LO =", wi.getR(WGroup::B, WCoef::C9, QCDOrder::LO, ContributionType::TOTAL));
    // LOG_INFO("AGAIN");
    // LOG_INFO("C9(mu_h) at NLO =", wi.getR(WGroup::B, WCoef::C9, QCDOrder::NLO, ContributionType::TOTAL));
    // LOG_INFO("C9(mu_h) at NNLO =", wi.getR(WGroup::B, WCoef::C9, QCDOrder::NNLO, ContributionType::TOTAL));
    
    // LOG_INFO("C9(mu_h) full =", wi.getFR(WGroup::B, WCoef::C9, QCDOrder::NNLO, ContributionType::TOTAL));
    return 0;
}