# MSSM: MSSM.nb

File MSSM.nb, 13.0 KB (added by BenjF, 9 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[ 13172, 405] |

14 | NotebookOptionsPosition[ 11175, 330] |

15 | NotebookOutlinePosition[ 11533, 346] |

16 | CellTagsIndexPosition[ 11490, 343] |

17 | WindowFrame->Normal*) |

18 | |

19 | (* Beginning of Notebook Content *) |

20 | Notebook[{ |

21 | |

22 | Cell[CellGroupData[{ |

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

24 | |

25 | Cell[BoxData[{ |

26 | RowBox[{ |

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

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

29 | RowBox[{ |

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

31 | RowBox[{ |

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

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

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

35 | RowBox[{ |

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

37 | RowBox[{ |

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

39 | RowBox[{ |

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

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

42 | }, Open ]], |

43 | |

44 | Cell[CellGroupData[{ |

45 | |

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

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

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

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

50 | |

51 | Cell[CellGroupData[{ |

52 | |

53 | Cell["\<\ |

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

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

56 | \>", "Subsubsection", |

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

58 | |

59 | Cell[BoxData[ |

60 | RowBox[{ |

61 | RowBox[{"If", "[", |

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

63 | RowBox[{ |

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

65 | RowBox[{ |

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

67 | ",", |

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

69 | }, Open ]], |

70 | |

71 | Cell[CellGroupData[{ |

72 | |

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

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

75 | |

76 | Cell[BoxData[{ |

77 | RowBox[{ |

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

79 | RowBox[{ |

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

81 | RowBox[{ |

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

83 | RowBox[{ |

84 | RowBox[{"LSUSYBr", "=", "LSoft"}], ";"}], "\[IndentingNewLine]", |

85 | RowBox[{ |

86 | RowBox[{"LGaugeFixing", "=", "LFeynmanGFix"}], ";"}]}], "Input", |

87 | CellChangeTimes->{ |

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

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

90 | 3.506864211817175*^9}, {3.510511899307577*^9, 3.510511913390203*^9}}] |

91 | }, Open ]], |

92 | |

93 | Cell[CellGroupData[{ |

94 | |

95 | Cell["\<\ |

96 | Unitarity gauge (no ghosts and Goldstones, if the flag $FeynmanGauge was not \ |

97 | set correctly)\ |

98 | \>", "Subsubsection", |

99 | CellChangeTimes->{{3.508650830315112*^9, 3.508650832348791*^9}, { |

100 | 3.510511936427951*^9, 3.510511947862441*^9}}], |

101 | |

102 | Cell[BoxData[ |

103 | RowBox[{ |

104 | RowBox[{"LagFeynman", "=", |

105 | RowBox[{"lagr", "/.", |

106 | RowBox[{"{", |

107 | RowBox[{ |

108 | RowBox[{ |

109 | RowBox[{"ghG", "[", "__", "]"}], "\[Rule]", "0"}], ",", " ", |

110 | RowBox[{ |

111 | RowBox[{"ghGbar", "[", "__", "]"}], "\[Rule]", "0"}], ",", |

112 | RowBox[{"ghWp", "\[Rule]", "0"}], ",", |

113 | RowBox[{"ghWpbar", "\[Rule]", "0"}], ",", |

114 | RowBox[{"ghWmbar", "\[Rule]", "0"}], ",", |

115 | RowBox[{"ghWm", "\[Rule]", "0"}], ",", |

116 | RowBox[{"ghZ", "\[Rule]", "0"}], ",", |

117 | RowBox[{"ghZbar", "\[Rule]", "0"}], ",", |

118 | RowBox[{"ghA", "\[Rule]", "0"}], ",", |

119 | RowBox[{"ghAbar", "\[Rule]", "0"}], ",", " ", |

120 | RowBox[{"G0", "\[Rule]", "0"}], ",", |

121 | RowBox[{"GP", "\[Rule]", "0"}], ",", |

122 | RowBox[{"GPbar", "\[Rule]", "0"}]}], "}"}]}]}], ";"}]], "Input", |

123 | CellChangeTimes->{{3.507308829401106*^9, 3.5073089177178507`*^9}, |

124 | 3.507310603266782*^9}] |

125 | }, Open ]], |

126 | |

127 | Cell[CellGroupData[{ |

128 | |

129 | Cell["\<\ |

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

131 | compatible with the CalcHep interface)\ |

132 | \>", "Subsubsection", |

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

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

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

136 | |

137 | Cell[BoxData[{ |

138 | RowBox[{ |

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

140 | RowBox[{ |

141 | RowBox[{"ReadLHAFile", "[", |

142 | RowBox[{"Input", "\[Rule]", "\"\<sps1a.dat\>\""}], "]"}], ";"}]}], "Input",\ |

143 | |

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

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

146 | 3.510511927139632*^9}, 3.510511971795094*^9}] |

147 | }, Open ]], |

148 | |

149 | Cell[CellGroupData[{ |

150 | |

151 | Cell["\<\ |

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

153 | compatible with the Calchep interface)\ |

154 | \>", "Subsubsection", |

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

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

157 | 3.510514606763693*^9}, {3.5105147982031612`*^9, 3.510514828422709*^9}}], |

158 | |

159 | Cell[BoxData[{ |

160 | RowBox[{ |

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

162 | RowBox[{ |

163 | RowBox[{"ReadLHAFile", "[", |

164 | RowBox[{"Input", "\[Rule]", "\"\<sps1a_wo.dat\>\""}], "]"}], |

165 | ";"}]}], "Input", |

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

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

168 | 3.510511927139632*^9}, 3.510511971795094*^9, {3.510514852406598*^9, |

169 | 3.510514852821206*^9}}] |

170 | }, Open ]] |

171 | }, Open ]], |

172 | |

173 | Cell[CellGroupData[{ |

174 | |

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

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

177 | |

178 | Cell[BoxData[ |

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

180 | RowBox[{"FeynmanRules", "[", |

181 | RowBox[{"lagr", ",", |

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

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

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

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

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

187 | 3.510572566899229*^9}}] |

188 | }, Open ]], |

189 | |

190 | Cell[CellGroupData[{ |

191 | |

192 | Cell[TextData[{ |

193 | Cell[BoxData[ |

194 | FormBox[ |

195 | StyleBox[ |

196 | RowBox[{"L", |

197 | StyleBox[ |

198 | AdjustmentBox["A", |

199 | BoxBaselineShift->-0.4, |

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

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

202 | AdjustmentBox["E", |

203 | BoxBaselineShift->0.5, |

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

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

206 | " Output" |

207 | }], "Section", |

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

209 | |

210 | Cell[BoxData[ |

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

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

213 | }, Open ]], |

214 | |

215 | Cell[CellGroupData[{ |

216 | |

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

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

219 | |

220 | Cell[CellGroupData[{ |

221 | |

222 | Cell["Optimization", "Subsubsection", |

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

224 | |

225 | Cell[BoxData[ |

226 | RowBox[{ |

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

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

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

230 | |

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

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

233 | }, Open ]], |

234 | |

235 | Cell[CellGroupData[{ |

236 | |

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

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

239 | |

240 | Cell[BoxData[{ |

241 | RowBox[{ |

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

243 | RowBox[{ |

244 | RowBox[{"WriteUFO", "[", |

245 | RowBox[{"LagFeynman", ",", " ", |

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

247 | |

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

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

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

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

252 | }, Open ]], |

253 | |

254 | Cell[CellGroupData[{ |

255 | |

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

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

258 | |

259 | Cell[BoxData[{ |

260 | RowBox[{ |

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

262 | RowBox[{ |

263 | RowBox[{"WriteCHOutput", "[", |

264 | RowBox[{"lagr", ",", " ", |

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

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

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

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

269 | }, Open ]], |

270 | |

271 | Cell[CellGroupData[{ |

272 | |

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

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

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

276 | |

277 | Cell[BoxData[{ |

278 | RowBox[{ |

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

280 | RowBox[{ |

281 | RowBox[{"WriteMGOutput", "[", |

282 | RowBox[{"lagr", ",", " ", |

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

284 | |

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

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

287 | 3.510513955605171*^9}}] |

288 | }, Open ]], |

289 | |

290 | Cell[CellGroupData[{ |

291 | |

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

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

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

295 | 3.510514891766358*^9}}], |

296 | |

297 | Cell[BoxData[{ |

298 | RowBox[{ |

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

300 | RowBox[{ |

301 | RowBox[{"WriteWOOutput", "[", |

302 | RowBox[{"lagr", ",", " ", |

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

304 | |

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

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

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

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

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

310 | }, Open ]], |

311 | |

312 | Cell[CellGroupData[{ |

313 | |

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

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

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

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

318 | |

319 | Cell[BoxData[{ |

320 | RowBox[{ |

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

322 | RowBox[{ |

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

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

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

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

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

328 | }, Open ]] |

329 | }, Open ]] |

330 | }, |

331 | WindowSize->{1425, 913}, |

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

333 | ShowSelection->True, |

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

335 | StyleDefinitions->"Default.nb" |

336 | ] |

337 | (* End of Notebook Content *) |

338 | |

339 | (* Internal cache information *) |

340 | (*CellTagsOutline |

341 | CellTagsIndex->{} |

342 | *) |

343 | (*CellTagsIndex |

344 | CellTagsIndex->{} |

345 | *) |

346 | (*NotebookFileOutline |

347 | Notebook[{ |

348 | Cell[CellGroupData[{ |

349 | Cell[567, 22, 24, 0, 73, "Section"], |

350 | Cell[594, 24, 585, 16, 143, "Input"] |

351 | }, Open ]], |

352 | Cell[CellGroupData[{ |

353 | Cell[1216, 45, 260, 3, 73, "Section"], |

354 | Cell[CellGroupData[{ |

355 | Cell[1501, 52, 243, 4, 29, "Subsubsection"], |

356 | Cell[1747, 58, 316, 9, 55, "Input"] |

357 | }, Open ]], |

358 | Cell[CellGroupData[{ |

359 | Cell[2100, 72, 140, 1, 29, "Subsubsection"], |

360 | Cell[2243, 75, 631, 14, 121, "Input"] |

361 | }, Open ]], |

362 | Cell[CellGroupData[{ |

363 | Cell[2911, 94, 241, 5, 29, "Subsubsection"], |

364 | Cell[3155, 101, 911, 22, 32, "Input"] |

365 | }, Open ]], |

366 | Cell[CellGroupData[{ |

367 | Cell[4103, 128, 360, 6, 29, "Subsubsection"], |

368 | Cell[4466, 136, 414, 9, 55, "Input"] |

369 | }, Open ]], |

370 | Cell[CellGroupData[{ |

371 | Cell[4917, 150, 357, 6, 29, "Subsubsection"], |

372 | Cell[5277, 158, 468, 10, 55, "Input"] |

373 | }, Open ]] |

374 | }, Open ]], |

375 | Cell[CellGroupData[{ |

376 | Cell[5794, 174, 112, 1, 73, "Section"], |

377 | Cell[5909, 177, 463, 9, 32, "Input"] |

378 | }, Open ]], |

379 | Cell[CellGroupData[{ |

380 | Cell[6409, 191, 461, 16, 75, "Section"], |

381 | Cell[6873, 209, 130, 2, 32, "Input"] |

382 | }, Open ]], |

383 | Cell[CellGroupData[{ |

384 | Cell[7040, 216, 139, 1, 73, "Section"], |

385 | Cell[CellGroupData[{ |

386 | Cell[7204, 221, 103, 1, 29, "Subsubsection"], |

387 | Cell[7310, 224, 356, 7, 32, "Input"] |

388 | }, Open ]], |

389 | Cell[CellGroupData[{ |

390 | Cell[7703, 236, 104, 1, 29, "Subsubsection"], |

391 | Cell[7810, 239, 520, 11, 55, "Input"] |

392 | }, Open ]], |

393 | Cell[CellGroupData[{ |

394 | Cell[8367, 255, 110, 1, 29, "Subsubsection"], |

395 | Cell[8480, 258, 429, 9, 55, "Input"] |

396 | }, Open ]], |

397 | Cell[CellGroupData[{ |

398 | Cell[8946, 272, 160, 2, 29, "Subsubsection"], |

399 | Cell[9109, 276, 426, 10, 55, "Input"] |

400 | }, Open ]], |

401 | Cell[CellGroupData[{ |

402 | Cell[9572, 291, 208, 3, 29, "Subsubsection"], |

403 | Cell[9783, 296, 615, 12, 55, "Input"] |

404 | }, Open ]], |

405 | Cell[CellGroupData[{ |

406 | Cell[10435, 313, 259, 3, 29, "Subsubsection"], |

407 | Cell[10697, 318, 450, 8, 55, "Input"] |

408 | }, Open ]] |

409 | }, Open ]] |

410 | } |

411 | ] |

412 | *) |

413 | |

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