Guest User

Untitled

a guest
Jan 22nd, 2018
70
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.44 KB | None | 0 0
  1. Where Im@qpC32 != 0 do qpC32=0 && Plot[Re@qpC32, {q, 0, 2000}, {m5, 0, 2000}]
  2. Where Im@qpC34 != 0 do qpC34=0 && Plot[Re@qpC34, {q, 0, 2000}, {m5, 0, 2000}]
  3.  
  4. Plot[Re@qpC32+Re@qpC34, {q, 0, 2000}, {m5, 0, 2000}]
  5.  
  6. eqn = j -
  7. Sqrt[q^2 + qp^2 -
  8. 2 q qp Cos[[Theta]]] - [Sqrt](qp^2 +
  9. 1/2 (16 m5^2 + ma^2 + mp^2 -
  10. Sqrt[(-(16 m5^2) - ma^2 - mp^2)^2 -
  11. 4 (ma^2 mp^2 - 16 m5^2 qp^2)])) == 0;
  12. With[{gensol = Solve[eqn , qp]},
  13. Block[{[Theta] = Pi/12, m = 5.5, M = 300, Nc = 3, c = !(*
  14. TagBox[
  15. InterpretationBox[
  16. RowBox[{""<-4.46874>"", "*",
  17. SuperscriptBox["10", ""<4>""]}],
  18. -44687.3983417778,
  19. AutoDelete->True],
  20. ScientificForm]), b = !(*
  21. TagBox[
  22. InterpretationBox[
  23. RowBox[{""<1.61594>"", "*",
  24. SuperscriptBox["10", ""<5>""]}],
  25. 161593.81818181818`,
  26. AutoDelete->True],
  27. ScientificForm]), k1 = !(*
  28. TagBox[
  29. InterpretationBox[
  30. RowBox[{""<1.6485>"", "*",
  31. SuperscriptBox["10", ""<1>""]}],
  32. 16.485010961790245`,
  33. AutoDelete->True],
  34. ScientificForm]), k2 = !(*
  35. TagBox[
  36. InterpretationBox[
  37. RowBox[{""<-1.31313>"", "*",
  38. SuperscriptBox["10", ""<1>""]}],
  39. -13.131344420001051`,
  40. AutoDelete->True],
  41. ScientificForm]), ma, mp,
  42. j},(*subs vals when gensol is evaluated*){j = [Sqrt](q^2 +
  43. 1/2 (16 m5^2 + ma^2 + mp^2 +
  44. Sqrt[(-(16 m5^2) - ma^2 - mp^2)^2 -
  45. 4 (ma^2 mp^2 - 16 m5^2 q^2)])),
  46. ma = [Sqrt](-2 (M^2 -
  47. 2 (3 k1 + k2) (Sqrt[(c + M^2 + 2 m5^2)/(2 (k1 + k2))] +
  48. m b/(2 (c + M^2 + 2 m5^2)))^2 - c + 2 m5^2)),
  49. mp = Sqrt[
  50. 2 b m (Sqrt[(c + M^2 + 2 m5^2)/(2 (k1 + k2))] +
  51. m b/(2 (c + M^2 + 2 m5^2)))^-1]};
  52. sols = gensol]];
  53. qpC32 = Compile[{{q, _Complex}, {m5, _Complex}},
  54. Evaluate[qp /. sols[[2]]],
  55. RuntimeOptions -> "EvaluateSymbolically" -> False] ;
  56. qpC34 = Compile[{{q, _Complex}, {m5, _Complex}},
  57. Evaluate[qp /. sols[[4]]],
  58. RuntimeOptions -> "EvaluateSymbolically" -> False] ;
  59. Plot3D[ Re@qpC34[q, m5], {q, 0, 2000}, {m5, 0, 2000},
  60. PlotRange -> Full]
  61. Plot3D[ Im@qpC34[q, m5], {q, 0, 2000}, {m5, 0, 2000},
  62. PlotRange -> Full]
  63. Plot3D[ Re@qpC32[q, m5], {q, 0, 2000}, {m5, 0, 2000},
  64. PlotRange -> Full]
  65. Plot3D[ Im@qpC32[q, m5], {q, 0, 2000}, {m5, 0, 2000},
  66. PlotRange -> Full]
  67.  
  68. rg32 = Function[{q, m5}, Chop[Im[qpC32[q, m5]]] == 0];
  69. rg34 = Function[{q, m5}, Chop[Im[qpC34[q, m5]]] == 0];
  70. plt32 = Plot3D[Re@qpC32[q, m5], {q, 0, 2000}, {m5, 0, 2000}, RegionFunction -> rg32];
  71. plt34 = Plot3D[Re@qpC34[q, m5], {q, 0, 2000}, {m5, 0, 2000}, RegionFunction -> rg34];
Add Comment
Please, Sign In to add comment