# RMSSM: MRSSM.nb

File MRSSM.nb, 11.1 KB (added by BenjF, 5 years ago) |
---|

Line | |
---|---|

1 | (* Content-type: application/mathematica *) |

2 | |

3 | (*** Wolfram Notebook File ***) |

4 | (* http://www.wolfram.com/nb *) |

5 | |

6 | (* CreatedBy='Mathematica 6.0' *) |

7 | |

8 | (*CacheID: 234*) |

9 | (* Internal cache information: |

10 | NotebookFileLineBreakTest |

11 | NotebookFileLineBreakTest |

12 | NotebookDataPosition[ 145, 7] |

13 | NotebookDataLength[ 11173, 347] |

14 | NotebookOptionsPosition[ 9378, 279] |

15 | NotebookOutlinePosition[ 9736, 295] |

16 | CellTagsIndexPosition[ 9693, 292] |

17 | WindowFrame->Normal*) |

18 | |

19 | (* Beginning of Notebook Content *) |

20 | Notebook[{ |

21 | Cell[BoxData[ |

22 | RowBox[{"Quit", "[", "]"}]], "Input", |

23 | CellChangeTimes->{{3.511874428548209*^9, 3.511874430667358*^9}}], |

24 | |

25 | Cell[CellGroupData[{ |

26 | |

27 | Cell["Setup", "Section"], |

28 | |

29 | Cell[BoxData[{ |

30 | RowBox[{ |

31 | RowBox[{"$OldDir", "=", |

32 | RowBox[{"Directory", "[", "]"}]}], ";"}], "\n", |

33 | RowBox[{ |

34 | RowBox[{"$FeynRulesPath", "=", |

35 | RowBox[{ |

36 | "SetDirectory", "[", "\"\<~/FeynRules/trunk/feynrules-development\>\"", |

37 | "]"}]}], ";"}], "\n", |

38 | RowBox[{"<<", "FeynRules`"}], "\n", |

39 | RowBox[{ |

40 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\n", |

41 | RowBox[{ |

42 | RowBox[{"LoadModel", "[", "\"\<mrssm.fr\>\"", "]"}], ";"}], "\n", |

43 | RowBox[{ |

44 | RowBox[{"$CalcLag", "=", "False"}], ";"}]}], "Input", |

45 | CellChangeTimes->{{3.510511714359152*^9, 3.5105117232569723`*^9}, { |

46 | 3.511873979058074*^9, 3.511873982075786*^9}, {3.537793110077643*^9, |

47 | 3.537793110571989*^9}}] |

48 | }, Open ]], |

49 | |

50 | Cell[CellGroupData[{ |

51 | |

52 | Cell["Lagrangian computation and update of the parameters", "Section", |

53 | CellChangeTimes->{{3.4912194701600723`*^9, 3.4912194715325117`*^9}, { |

54 | 3.50686100488665*^9, 3.506861008141419*^9}, {3.506866546045188*^9, |

55 | 3.50686654675209*^9}, 3.510511732200982*^9}], |

56 | |

57 | Cell[CellGroupData[{ |

58 | |

59 | Cell["\<\ |

60 | Calculation of the Lagrangian itself. This can take ~10 minutes. Therefore, \ |

61 | it is saved in a file which can be read each time it is necessary.\ |

62 | \>", "Subsubsection", |

63 | CellChangeTimes->{{3.510511740551598*^9, 3.510511755990991*^9}}], |

64 | |

65 | Cell[BoxData[ |

66 | RowBox[{ |

67 | RowBox[{"If", "[", |

68 | RowBox[{"$CalcLag", ",", |

69 | RowBox[{ |

70 | RowBox[{"lagr", "=", "Lag"}], ";", "\[IndentingNewLine]", |

71 | RowBox[{ |

72 | RowBox[{"Definition", "[", "lagr", "]"}], ">>", "lagrangian.dat"}]}], |

73 | ",", |

74 | RowBox[{"<<", "lagrangian.dat"}]}], "]"}], ";"}]], "Input"] |

75 | }, Open ]], |

76 | |

77 | Cell[CellGroupData[{ |

78 | |

79 | Cell["Access to the different pieces of the Lagrangian", "Subsubsection", |

80 | CellChangeTimes->{{3.510511759479219*^9, 3.510511767623085*^9}}], |

81 | |

82 | Cell[BoxData[{ |

83 | RowBox[{ |

84 | RowBox[{"LVec", "=", "LVector"}], ";"}], "\[IndentingNewLine]", |

85 | RowBox[{ |

86 | RowBox[{"LC", "=", "LChiral"}], ";"}], "\[IndentingNewLine]", |

87 | RowBox[{ |

88 | RowBox[{"LW", "=", "LSuperW"}], ";"}], "\[IndentingNewLine]", |

89 | RowBox[{ |

90 | RowBox[{"LSUSYBr", "=", "LSoft"}], ";"}]}], "Input", |

91 | CellChangeTimes->{ |

92 | 3.506861831832621*^9, {3.50686396987472*^9, 3.50686397956723*^9}, { |

93 | 3.506864125221086*^9, 3.506864126560401*^9}, {3.506864198923688*^9, |

94 | 3.506864211817175*^9}, {3.510511899307577*^9, 3.510511913390203*^9}, |

95 | 3.5121224264764137`*^9, {3.512122465053269*^9, 3.512122467303591*^9}}] |

96 | }, Open ]], |

97 | |

98 | Cell[CellGroupData[{ |

99 | |

100 | Cell["\<\ |

101 | Update of the parameters (not compatible with the Whizard interface, but \ |

102 | compatible with the CalcHep interface)\ |

103 | \>", "Subsubsection", |

104 | CellChangeTimes->{{3.508650840502183*^9, 3.508650841544434*^9}, { |

105 | 3.510511953402149*^9, 3.5105119585185213`*^9}, {3.5105145827435713`*^9, |

106 | 3.510514606763693*^9}, {3.510514811212164*^9, 3.510514819639832*^9}}], |

107 | |

108 | Cell[BoxData[{ |

109 | RowBox[{ |

110 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]", |

111 | RowBox[{ |

112 | RowBox[{"ReadLHAFile", "[", |

113 | RowBox[{"Input", "\[Rule]", "\"\<spectrum.dat\>\""}], "]"}], |

114 | ";"}]}], "Input", |

115 | CellChangeTimes->{{3.507310458304277*^9, 3.507310473976326*^9}, { |

116 | 3.507310607848939*^9, 3.507310608185873*^9}, {3.510511925987886*^9, |

117 | 3.510511927139632*^9}, 3.510511971795094*^9, {3.5439086684386797`*^9, |

118 | 3.543908669966134*^9}}] |

119 | }, Open ]] |

120 | }, Open ]], |

121 | |

122 | Cell[CellGroupData[{ |

123 | |

124 | Cell["Feynman rules calculation", "Section", |

125 | CellChangeTimes->{{3.508650859336681*^9, 3.5086508623247147`*^9}}], |

126 | |

127 | Cell[BoxData[ |

128 | RowBox[{"rules", "=", |

129 | RowBox[{"FeynmanRules", "[", |

130 | RowBox[{"lagr", ",", |

131 | RowBox[{"ScreenOutput", "\[Rule]", "False"}], ",", |

132 | RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}]}]], "Input", |

133 | CellChangeTimes->{{3.508061383518343*^9, 3.508061413431903*^9}, { |

134 | 3.510571497855986*^9, 3.510571568813633*^9}, {3.510572173986182*^9, |

135 | 3.5105721743854856`*^9}, 3.510572210161304*^9, {3.510572541184834*^9, |

136 | 3.510572566899229*^9}}] |

137 | }, Open ]], |

138 | |

139 | Cell[CellGroupData[{ |

140 | |

141 | Cell[TextData[{ |

142 | Cell[BoxData[ |

143 | FormBox[ |

144 | StyleBox[ |

145 | RowBox[{"L", |

146 | StyleBox[ |

147 | AdjustmentBox["A", |

148 | BoxBaselineShift->-0.4, |

149 | BoxMargins->{{-0.5, -0.3}, {0, 0}}], |

150 | FontSize->Smaller], "T", |

151 | AdjustmentBox["E", |

152 | BoxBaselineShift->0.5, |

153 | BoxMargins->{{-0.3, 0}, {0, 0}}], "X"}], |

154 | SingleLetterItalics->False], TraditionalForm]]], |

155 | " Output" |

156 | }], "Section", |

157 | CellChangeTimes->{{3.510516585335*^9, 3.510516590423654*^9}}], |

158 | |

159 | Cell[BoxData[ |

160 | RowBox[{"WriteLaTeXOutput", "[", "]"}]], "Input", |

161 | CellChangeTimes->{{3.510516592624229*^9, 3.510516605419121*^9}}] |

162 | }, Open ]], |

163 | |

164 | Cell[CellGroupData[{ |

165 | |

166 | Cell["Generating MC model files with optimized restrictions", "Section", |

167 | CellChangeTimes->{{3.51051202409053*^9, 3.5105120509193287`*^9}}], |

168 | |

169 | Cell[CellGroupData[{ |

170 | |

171 | Cell["Optimization", "Subsubsection", |

172 | CellChangeTimes->{{3.508650844785081*^9, 3.508650846244998*^9}}], |

173 | |

174 | Cell[BoxData[ |

175 | RowBox[{ |

176 | RowBox[{"WriteRestrictionFile", "[", "]"}], ";", " ", |

177 | RowBox[{"LoadRestriction", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";", " ", |

178 | RowBox[{"DeleteFile", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";"}]], "Input",\ |

179 | |

180 | CellChangeTimes->{{3.5080579535629387`*^9, 3.508057961830206*^9}, { |

181 | 3.508058109308621*^9, 3.508058122229617*^9}}] |

182 | }, Open ]], |

183 | |

184 | Cell[CellGroupData[{ |

185 | |

186 | Cell["UFO interface", "Subsubsection", |

187 | CellChangeTimes->{{3.508650849632307*^9, 3.508650853575732*^9}}], |

188 | |

189 | Cell[BoxData[{ |

190 | RowBox[{ |

191 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]", |

192 | RowBox[{ |

193 | RowBox[{"WriteUFO", "[", |

194 | RowBox[{"LagFeynman", ",", " ", |

195 | RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\ |

196 | |

197 | CellChangeTimes->{{3.5073089286898727`*^9, 3.507308937458375*^9}, { |

198 | 3.507311910076086*^9, 3.507311939653351*^9}, {3.507398897853717*^9, |

199 | 3.507398901626676*^9}, {3.507399273736373*^9, 3.507399274557005*^9}, { |

200 | 3.510512090662549*^9, 3.510512091990087*^9}}] |

201 | }, Open ]], |

202 | |

203 | Cell[CellGroupData[{ |

204 | |

205 | Cell["CalcHEP model files", "Subsubsection", |

206 | CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}}], |

207 | |

208 | Cell[BoxData[{ |

209 | RowBox[{ |

210 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]", |

211 | RowBox[{ |

212 | RowBox[{"WriteCHOutput", "[", |

213 | RowBox[{"lagr", ",", " ", |

214 | RowBox[{"Exclude4Scalars", "\[Rule]", "True"}], ",", |

215 | RowBox[{"CHSimplify", "\[Rule]", "False"}]}], "]"}], ";"}]}], "Input", |

216 | CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, { |

217 | 3.5105127146047497`*^9, 3.510512717826407*^9}}] |

218 | }, Open ]], |

219 | |

220 | Cell[CellGroupData[{ |

221 | |

222 | Cell["MadGraph model files", "Subsubsection", |

223 | CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, { |

224 | 3.510513948444434*^9, 3.510513949301979*^9}}], |

225 | |

226 | Cell[BoxData[{ |

227 | RowBox[{ |

228 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]", |

229 | RowBox[{ |

230 | RowBox[{"WriteMGOutput", "[", |

231 | RowBox[{"lagr", ",", " ", |

232 | RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\ |

233 | |

234 | CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, { |

235 | 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9, |

236 | 3.510513955605171*^9}}] |

237 | }, Open ]], |

238 | |

239 | Cell[CellGroupData[{ |

240 | |

241 | Cell["Whizard model files", "Subsubsection", |

242 | CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, { |

243 | 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9, |

244 | 3.510514891766358*^9}}], |

245 | |

246 | Cell[BoxData[{ |

247 | RowBox[{ |

248 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]", |

249 | RowBox[{ |

250 | RowBox[{"WriteWOOutput", "[", |

251 | RowBox[{"lagr", ",", " ", |

252 | RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\ |

253 | |

254 | CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, { |

255 | 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9, |

256 | 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, { |

257 | 3.510572235119359*^9, 3.51057223551847*^9}, {3.510572377588068*^9, |

258 | 3.510572378543983*^9}, {3.510572531362316*^9, 3.510572531820734*^9}}] |

259 | }, Open ]], |

260 | |

261 | Cell[CellGroupData[{ |

262 | |

263 | Cell["FeynArts model files", "Subsubsection", |

264 | CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, { |

265 | 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9, |

266 | 3.510514891766358*^9}, {3.5105162477375813`*^9, 3.5105162494145308`*^9}}], |

267 | |

268 | Cell[BoxData[{ |

269 | RowBox[{ |

270 | RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]", |

271 | RowBox[{ |

272 | RowBox[{"WriteFeynArtsOutput", "[", "lagr", "]"}], ";"}]}], "Input", |

273 | CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, { |

274 | 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9, |

275 | 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, { |

276 | 3.5105162369063063`*^9, 3.510516243526369*^9}}] |

277 | }, Open ]] |

278 | }, Open ]] |

279 | }, |

280 | WindowSize->{1440, 825}, |

281 | WindowMargins->{{3, Automatic}, {Automatic, 23}}, |

282 | ShowSelection->True, |

283 | FrontEndVersion->"7.0 for Linux x86 (64-bit) (February 25, 2009)", |

284 | StyleDefinitions->"Default.nb" |

285 | ] |

286 | (* End of Notebook Content *) |

287 | |

288 | (* Internal cache information *) |

289 | (*CellTagsOutline |

290 | CellTagsIndex->{} |

291 | *) |

292 | (*CellTagsIndex |

293 | CellTagsIndex->{} |

294 | *) |

295 | (*NotebookFileOutline |

296 | Notebook[{ |

297 | Cell[545, 20, 118, 2, 32, "Input"], |

298 | Cell[CellGroupData[{ |

299 | Cell[688, 26, 24, 0, 73, "Section"], |

300 | Cell[715, 28, 684, 18, 143, "Input"] |

301 | }, Open ]], |

302 | Cell[CellGroupData[{ |

303 | Cell[1436, 51, 260, 3, 73, "Section"], |

304 | Cell[CellGroupData[{ |

305 | Cell[1721, 58, 243, 4, 29, "Subsubsection"], |

306 | Cell[1967, 64, 316, 9, 55, "Input"] |

307 | }, Open ]], |

308 | Cell[CellGroupData[{ |

309 | Cell[2320, 78, 140, 1, 29, "Subsubsection"], |

310 | Cell[2463, 81, 615, 13, 99, "Input"] |

311 | }, Open ]], |

312 | Cell[CellGroupData[{ |

313 | Cell[3115, 99, 360, 6, 29, "Subsubsection"], |

314 | Cell[3478, 107, 470, 10, 55, "Input"] |

315 | }, Open ]] |

316 | }, Open ]], |

317 | Cell[CellGroupData[{ |

318 | Cell[3997, 123, 112, 1, 73, "Section"], |

319 | Cell[4112, 126, 463, 9, 32, "Input"] |

320 | }, Open ]], |

321 | Cell[CellGroupData[{ |

322 | Cell[4612, 140, 461, 16, 75, "Section"], |

323 | Cell[5076, 158, 130, 2, 32, "Input"] |

324 | }, Open ]], |

325 | Cell[CellGroupData[{ |

326 | Cell[5243, 165, 139, 1, 73, "Section"], |

327 | Cell[CellGroupData[{ |

328 | Cell[5407, 170, 103, 1, 29, "Subsubsection"], |

329 | Cell[5513, 173, 356, 7, 32, "Input"] |

330 | }, Open ]], |

331 | Cell[CellGroupData[{ |

332 | Cell[5906, 185, 104, 1, 29, "Subsubsection"], |

333 | Cell[6013, 188, 520, 11, 55, "Input"] |

334 | }, Open ]], |

335 | Cell[CellGroupData[{ |

336 | Cell[6570, 204, 110, 1, 29, "Subsubsection"], |

337 | Cell[6683, 207, 429, 9, 55, "Input"] |

338 | }, Open ]], |

339 | Cell[CellGroupData[{ |

340 | Cell[7149, 221, 160, 2, 29, "Subsubsection"], |

341 | Cell[7312, 225, 426, 10, 55, "Input"] |

342 | }, Open ]], |

343 | Cell[CellGroupData[{ |

344 | Cell[7775, 240, 208, 3, 29, "Subsubsection"], |

345 | Cell[7986, 245, 615, 12, 55, "Input"] |

346 | }, Open ]], |

347 | Cell[CellGroupData[{ |

348 | Cell[8638, 262, 259, 3, 29, "Subsubsection"], |

349 | Cell[8900, 267, 450, 8, 55, "Input"] |

350 | }, Open ]] |

351 | }, Open ]] |

352 | } |

353 | ] |

354 | *) |

355 | |

356 | (* End of internal cache information *) |