MSSMSimpleSpec.cpp
Go to the documentation of this file.
36 bool orderer (std::pair<int, double> a, std::pair<int, double> b) { return a.second < b.second; } 41 MSSMea::MSSMea() 79 errmsg << "Could not find block \"GAUGE\" in SLHAea object (required to retrieve scale Q). Received out_of_range error with message: " << e.what(); 88 pdg[0] = initVector<int>(1000001, 1000003, 1000005, 2000001, 2000003, 2000005); // d-type squarks 89 pdg[1] = initVector<int>(1000002, 1000004, 1000006, 2000002, 2000004, 2000006); // u-type squarks 95 for (int i = 0; i < lengths[j]; i++) masses[j].push_back(std::pair<int, double>(pdg[j][i], getdata("MASS",pdg[j][i]))); 140 //data[blocks[j]][""] << i+1 << k+1 << datum << "# "+blocks[j]+boost::lexical_cast<str>(i*10+k+11); 141 SLHAea_add(data, blocks[j], i+1, k+1, datum, "# "+blocks[j]+boost::lexical_cast<str>(i*10+k+11), true); 147 sspair M[5] = {sspair("MSL2","ml2"), sspair("MSE2","me2"), sspair("MSQ2","mq2"), sspair("MSU2","mu2"), sspair("MSD2","md2")}; 151 if(not SLHAea_block_exists(data, block)) SLHAea_add_block(data, block, scale); //TODO: maybe just always delete and replace 159 entry = getdata("MSOFT",30+3*k+i+(k>1?4:0)); // black magic to get correct index in MSOFT matching diagonal elements 171 // Yukawa and trilinear blocks. YU, YD and YE, plus [YU, YD and YE; SLHA1 only], or [TU, TD and TE; SLHA2 only]. 227 else logger() << "Input SLHA data for setting up simple spectrum is SLHA2. *living in the future*" << EOM; 228 // In the end, we should always have converted to SLHA2 by now, so set this as the internally-tracked version. 237 double MSSMea::get_mA2() const{ return getdata("HMIX",4); } // m^2_A=[m3^2/cosBsinB](Q) DRbar, tree 288 double MSSMea::get_g1() const { return getdata("GAUGE",1)/sqrt(3./5.); } // Convert from gy (in SLHAea object) to g1 303 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_Mhh_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 307 double MSSMea::get_MW_pole() const { return getdata("MASS",24); } // REQUIRED output of MSSM-compatible subspectrum 313 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_MCha_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 325 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_MSd_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 335 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_MSd_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 345 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_MSd_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 352 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_MSd_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 360 else { utils_error().raise(LOCAL_INFO,"Invalid index input to get_MChi_pole_slha! Please check index range limits in wrapper SubSpectrum class!"); return -1; } // Should not return. 409 int MSSMSimpleSpec::get_index_offset() const {return 0.;} // we use indices starting from 1 in this file, matching user assumptions. (because Peter is god, he knows user assumptions before they do.) 412 // NOTE! No need to write this function, SubSpectrum base class can handle it if add_to_SLHAea exists. 428 const std::map<int, int>& MSSMSimpleSpec::PDG_translator() const { return slhawrap.PDG_translator(); } 519 tmp_map["tanbeta(mZ)"]= &Model::get_tanbeta_mZ; // Special entry for reproducing MINPAR entry in SLHA double get_UP_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:377 virtual const std::map< int, int > & PDG_translator() const Retrieve the PDG translation map. Definition: MSSMSimpleSpec.cpp:428 Definition: spectrum_helpers.hpp:55 Definition: spectrum_helpers.hpp:54 Definition: log_tags.hpp:54 void set_pole_mass_uncertainties(double) Set pole mass uncertainties. Definition: MSSMSimpleSpec.cpp:431 virtual int get_index_offset() const Ofset from user-input indices (user assumes 1,2,3 indexed, e.g. use offset=-1 for zero-indexing) ... Definition: MSSMSimpleSpec.cpp:409 double get_Mhh_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:299 double get_ZA_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:371 double get_MChi_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:354 void SLHAea_add_block(SLHAstruct &, const str &name, const double scale=-1) Add a new block to an SLHAea object, with or without a scale. Definition: slhaea_helpers.cpp:85 double get_MSu_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:327 double get_MSv_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:347 double get_ZV_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:367 Logging access header for GAMBIT. void SLHAea_add_GAMBIT_SPINFO(SLHAstruct &slha) Write the SPINFO block with GAMBIT name and version number. Definition: slhaea_helpers.cpp:174 void SLHAea_add(SLHAstruct &slha, const str &block, const int index, const double value, const str &comment="", const bool overwrite=false) Add an entry to an SLHAea object (if overwrite=false, only if it doesn't already exist) ... Definition: slhaea_helpers.cpp:188 Routines to help translate between SLHA2 sfermions and SLHA1 (or similar) sfermions. EXPORT_SYMBOLS double sqr(double a) returns square of double - saves tedious repetition Definition: util_functions.cpp:391 General small utility functions. double getdata(const std::string &block, int index) const Helper functions to do error checking for SLHAea object contents. Definition: SLHASimpleSpec.cpp:79 void set_override_vector(const Par::Tags, const double, const std::vector< str > &, const bool allow_new=false, const bool decouple=false) Vector override functions. Definition: subspectrum.cpp:516 Definition: spectrum_helpers.hpp:50 void add_MSSM_spectrum_to_SLHAea(const SubSpectrum &mssmspec, SLHAstruct &slha, int slha_version) Add an entire MSSM spectrum to an SLHAea object. Definition: mssm_slhahelp.cpp:772 bool SLHAea_block_exists(SLHAstruct &slha, const str &block) Check if a block exists in an SLHAea object. Definition: slhaea_helpers.cpp:97 double get_ZD_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:364 double get_ZP_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:373 const Logging::endofmessage EOM Explicit const instance of the end of message struct in Gambit namespace. Definition: logger.hpp:100 EXPORT_SYMBOLS Logging::LogMaster & logger() Function to retrieve a reference to the Gambit global log object. Definition: logger.cpp:95 Variadic utilty functions. double get_sinthW2_DRbar(const Model &model) Definition: MSSMSpec.hpp:234 std::vector< T > initVector(std::vector< T > vector) Definition: variadic_functions.hpp:52 virtual void add_to_SLHAea(int, SLHAea::Coll &) const Retrieve SLHAea object. Definition: MSSMSimpleSpec.cpp:416 double get_ZH_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:370 Definition: SLHASimpleSpec.hpp:64 double get_MSd_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:315 static GetterMaps fill_getter_maps() Map fillers Used to initialise maps in Spec class, accessed via SubSpectrum interface class (speciali... Definition: MSSMSimpleSpec.cpp:455 double get_ZE_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:368 double get_ZU_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:365 MSSM specialisation of SLHAea object wrapper version of SubSpectrum class. Definition: MSSMSimpleSpec.hpp:121 bool orderer(std::pair< int, double > a, std::pair< int, double > b) Helper function for sorting int, double pairs according to the double. Definition: MSSMSimpleSpec.cpp:36 double get_ZN_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:374 double get_UM_pole_slha(int i, int j) const Definition: MSSMSimpleSpec.cpp:376 std::map< int, int > PDG_translation_map PDG translation map (e.g. from SLHA1 to SLHA2 for MSSMskeleton) Definition: SLHASimpleSpec.hpp:40 double get_MCha_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:309 Skeleton "model" class which interacts with an SLHAea object. Definition: SLHASimpleSpec.hpp:32 bool SLHAea_check_block(SLHAstruct &slha, const str &block) Check if a block exists in an SLHAea object, add it if not. Definition: slhaea_helpers.cpp:105 std::map< Par::Tags, MapCollection< MTget > > GetterMaps Will need a map of map collections for both the getters and setters, containing the map collections f... Definition: spec_head.hpp:248 double get_MSe_pole_slha(int i) const Definition: MSSMSimpleSpec.cpp:337 |