# NMSSM: NMSSM.nb

File NMSSM.nb, 13.0 KB (added by BenjF, 8 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[ 13142, 404] |

14 | NotebookOptionsPosition[ 11144, 329] |

15 | NotebookOutlinePosition[ 11503, 345] |

16 | CellTagsIndexPosition[ 11460, 342] |

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[{"SetDirectory", "[", "\"\<<path to FeynRules>\>\"", "]"}]}], |

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

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

34 | RowBox[{ |

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

36 | RowBox[{ |

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

38 | RowBox[{ |

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

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

41 | 3.5105947682179747`*^9, 3.510594774165135*^9}, {3.511517628737059*^9, |

42 | 3.51151762923162*^9}, {3.5115238780058804`*^9, 3.5115238807474947`*^9}}] |

43 | }, Open ]], |

44 | |

45 | Cell[CellGroupData[{ |

46 | |

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

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

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

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

51 | |

52 | Cell[CellGroupData[{ |

53 | |

54 | Cell["\<\ |

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

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

57 | \>", "Subsubsection", |

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

59 | |

60 | Cell[BoxData[ |

61 | RowBox[{ |

62 | RowBox[{"If", "[", |

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

64 | RowBox[{ |

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

66 | RowBox[{ |

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

68 | ",", |

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

70 | }, Open ]], |

71 | |

72 | Cell[CellGroupData[{ |

73 | |

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

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

76 | |

77 | Cell[BoxData[{ |

78 | RowBox[{ |

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

80 | RowBox[{ |

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

82 | RowBox[{ |

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

84 | RowBox[{ |

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

86 | RowBox[{ |

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

88 | CellChangeTimes->{ |

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

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

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

92 | }, Open ]], |

93 | |

94 | Cell[CellGroupData[{ |

95 | |

96 | Cell["\<\ |

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

98 | set correctly)\ |

99 | \>", "Subsubsection", |

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

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

102 | |

103 | Cell[BoxData[ |

104 | RowBox[{ |

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

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

107 | RowBox[{"{", |

108 | RowBox[{ |

109 | RowBox[{ |

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

111 | RowBox[{ |

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

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

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

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

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

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

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

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

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

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

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

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

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

125 | 3.507310603266782*^9}] |

126 | }, Open ]], |

127 | |

128 | Cell[CellGroupData[{ |

129 | |

130 | Cell["\<\ |

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

132 | compatible with the CalcHep interface)\ |

133 | \>", "Subsubsection", |

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

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

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

137 | |

138 | Cell[BoxData[{ |

139 | RowBox[{ |

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

141 | RowBox[{ |

142 | RowBox[{"ReadLHAFile", "[", |

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

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, {3.511520172569066*^9, |

147 | 3.511520173343334*^9}}] |

148 | }, Open ]], |

149 | |

150 | Cell[CellGroupData[{ |

151 | |

152 | Cell["\<\ |

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

154 | compatible with the Calchep interface)\ |

155 | \>", "Subsubsection", |

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

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

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

159 | |

160 | Cell[BoxData[{ |

161 | RowBox[{ |

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

163 | RowBox[{ |

164 | RowBox[{"ReadLHAFile", "[", |

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

166 | |

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

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

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

170 | 3.510514852821206*^9}, {3.51152299692323*^9, 3.511522998284626*^9}}] |

171 | }, Open ]] |

172 | }, Open ]], |

173 | |

174 | Cell[CellGroupData[{ |

175 | |

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

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

178 | |

179 | Cell[BoxData[ |

180 | RowBox[{ |

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

182 | RowBox[{"FeynmanRules", "[", |

183 | RowBox[{"lagr", ",", |

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

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

186 | ";"}]], "Input", |

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

188 | 3.510594846536268*^9, 3.5105948467914553`*^9}}] |

189 | }, Open ]], |

190 | |

191 | Cell[CellGroupData[{ |

192 | |

193 | Cell[TextData[{ |

194 | Cell[BoxData[ |

195 | FormBox[ |

196 | StyleBox[ |

197 | RowBox[{"L", |

198 | StyleBox[ |

199 | AdjustmentBox["A", |

200 | BoxBaselineShift->-0.4, |

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

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

203 | AdjustmentBox["E", |

204 | BoxBaselineShift->0.5, |

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

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

207 | " Output" |

208 | }], "Section", |

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

210 | |

211 | Cell[BoxData[ |

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

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

214 | }, Open ]], |

215 | |

216 | Cell[CellGroupData[{ |

217 | |

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

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

220 | |

221 | Cell[CellGroupData[{ |

222 | |

223 | Cell["Optimization", "Subsubsection", |

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

225 | |

226 | Cell[BoxData[ |

227 | RowBox[{ |

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

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

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

231 | |

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

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

234 | }, Open ]], |

235 | |

236 | Cell[CellGroupData[{ |

237 | |

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

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

240 | |

241 | Cell[BoxData[{ |

242 | RowBox[{ |

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

244 | RowBox[{ |

245 | RowBox[{"WriteUFO", "[", |

246 | RowBox[{"LagFeynman", ",", " ", |

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

248 | |

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

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

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

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

253 | }, Open ]], |

254 | |

255 | Cell[CellGroupData[{ |

256 | |

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

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

259 | |

260 | Cell[BoxData[{ |

261 | RowBox[{ |

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

263 | RowBox[{ |

264 | RowBox[{"WriteCHOutput", "[", |

265 | RowBox[{"lagr", ",", " ", |

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

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

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

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

270 | }, Open ]], |

271 | |

272 | Cell[CellGroupData[{ |

273 | |

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

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

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

277 | |

278 | Cell[BoxData[{ |

279 | RowBox[{ |

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

281 | RowBox[{ |

282 | RowBox[{"WriteMGOutput", "[", |

283 | RowBox[{"lagr", ",", " ", |

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

285 | |

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

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

288 | 3.510513955605171*^9}}] |

289 | }, Open ]], |

290 | |

291 | Cell[CellGroupData[{ |

292 | |

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

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

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

296 | 3.510514891766358*^9}}], |

297 | |

298 | Cell[BoxData[{ |

299 | RowBox[{ |

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

301 | RowBox[{ |

302 | RowBox[{"WriteWOOutput", "[", |

303 | RowBox[{"lagr", ",", " ", |

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

305 | |

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

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

308 | 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}}] |

309 | }, Open ]], |

310 | |

311 | Cell[CellGroupData[{ |

312 | |

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

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

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

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

317 | |

318 | Cell[BoxData[{ |

319 | RowBox[{ |

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

321 | RowBox[{ |

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

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

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

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

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

327 | }, Open ]] |

328 | }, Open ]] |

329 | }, |

330 | WindowSize->{1914, 1123}, |

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

332 | ShowSelection->True, |

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

334 | StyleDefinitions->"Default.nb" |

335 | ] |

336 | (* End of Notebook Content *) |

337 | |

338 | (* Internal cache information *) |

339 | (*CellTagsOutline |

340 | CellTagsIndex->{} |

341 | *) |

342 | (*CellTagsIndex |

343 | CellTagsIndex->{} |

344 | *) |

345 | (*NotebookFileOutline |

346 | Notebook[{ |

347 | Cell[CellGroupData[{ |

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

349 | Cell[594, 24, 709, 17, 143, "Input"] |

350 | }, Open ]], |

351 | Cell[CellGroupData[{ |

352 | Cell[1340, 46, 260, 3, 73, "Section"], |

353 | Cell[CellGroupData[{ |

354 | Cell[1625, 53, 243, 4, 29, "Subsubsection"], |

355 | Cell[1871, 59, 316, 9, 55, "Input"] |

356 | }, Open ]], |

357 | Cell[CellGroupData[{ |

358 | Cell[2224, 73, 140, 1, 29, "Subsubsection"], |

359 | Cell[2367, 76, 631, 14, 121, "Input"] |

360 | }, Open ]], |

361 | Cell[CellGroupData[{ |

362 | Cell[3035, 95, 241, 5, 29, "Subsubsection"], |

363 | Cell[3279, 102, 911, 22, 32, "Input"] |

364 | }, Open ]], |

365 | Cell[CellGroupData[{ |

366 | Cell[4227, 129, 360, 6, 29, "Subsubsection"], |

367 | Cell[4590, 137, 459, 9, 55, "Input"] |

368 | }, Open ]], |

369 | Cell[CellGroupData[{ |

370 | Cell[5086, 151, 357, 6, 29, "Subsubsection"], |

371 | Cell[5446, 159, 509, 10, 55, "Input"] |

372 | }, Open ]] |

373 | }, Open ]], |

374 | Cell[CellGroupData[{ |

375 | Cell[6004, 175, 112, 1, 73, "Section"], |

376 | Cell[6119, 178, 365, 9, 32, "Input"] |

377 | }, Open ]], |

378 | Cell[CellGroupData[{ |

379 | Cell[6521, 192, 461, 16, 75, "Section"], |

380 | Cell[6985, 210, 130, 2, 32, "Input"] |

381 | }, Open ]], |

382 | Cell[CellGroupData[{ |

383 | Cell[7152, 217, 139, 1, 73, "Section"], |

384 | Cell[CellGroupData[{ |

385 | Cell[7316, 222, 103, 1, 29, "Subsubsection"], |

386 | Cell[7422, 225, 356, 7, 32, "Input"] |

387 | }, Open ]], |

388 | Cell[CellGroupData[{ |

389 | Cell[7815, 237, 104, 1, 29, "Subsubsection"], |

390 | Cell[7922, 240, 520, 11, 55, "Input"] |

391 | }, Open ]], |

392 | Cell[CellGroupData[{ |

393 | Cell[8479, 256, 110, 1, 29, "Subsubsection"], |

394 | Cell[8592, 259, 429, 9, 55, "Input"] |

395 | }, Open ]], |

396 | Cell[CellGroupData[{ |

397 | Cell[9058, 273, 160, 2, 29, "Subsubsection"], |

398 | Cell[9221, 277, 426, 10, 55, "Input"] |

399 | }, Open ]], |

400 | Cell[CellGroupData[{ |

401 | Cell[9684, 292, 208, 3, 29, "Subsubsection"], |

402 | Cell[9895, 297, 472, 10, 55, "Input"] |

403 | }, Open ]], |

404 | Cell[CellGroupData[{ |

405 | Cell[10404, 312, 259, 3, 29, "Subsubsection"], |

406 | Cell[10666, 317, 450, 8, 55, "Input"] |

407 | }, Open ]] |

408 | }, Open ]] |

409 | } |

410 | ] |

411 | *) |

412 | |

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