NewIm[expression_] := (num = Numerator[expression];
   den = Denominator[expression];
   {reNum, imNum} = ComplexExpand[ReIm[num]]; {reDen, imDen} =
    ComplexExpand[ReIm[den]];
   imout = (imNum reDen - reNum imDen)/(reDen^2 + imDen^2));


NewRe[expression_] := (num = Numerator[expression];
   den = Denominator[expression]; {reNum, imNum} =
    ComplexExpand[ReIm[num]]; {reDen, imDen} =
    ComplexExpand[ReIm[den]];
   reout = (reNum reDen + imNum imDen)/(reDen^2 + imDen^2));



MaxFINDER[[Chi]_] :=
  NMaximize[
   Simplify[
    NewIm[[Chi]]], [CapitalOmega]c > .001 && [CapitalOmega]d >
     0, {[CapitalOmega]c, [CapitalOmega]d, [CapitalDelta]p,
[CapitalDelta]c, [CapitalDelta]s, [CapitalDelta]d, [Phi]},
   MaxIterations -> 10000];



pCs[Chi]2fun[[Beta]_]  := (I E^(I [Phi]) [CapitalOmega]c)/(
   8 (-([Gamma]a/2) + I [CapitalDelta]a) (([Beta] [Gamma]a)/2 + (
      I [Gamma]a [CapitalDelta]c)/2 - I [Beta] [CapitalDelta]p - (
      I [Gamma]a [CapitalDelta]p)/
      2 + [CapitalDelta]c [CapitalDelta]p - [CapitalDelta]p^2 +
[CapitalOmega]c^2/4)) /. {[Gamma]a ->
     1, [CapitalDelta]a -> [CapitalDelta]c + [CapitalDelta]p -
[CapitalDelta]s};

(*level5[Chi]2 = -([ExponentialE]^([ImaginaryI] [Phi])(50
[CapitalOmega]c [CapitalOmega]d)/(([Gamma]a^2+2 [ImaginaryI]
[Gamma]a (50 [CapitalDelta]c-51 [CapitalDelta]p)+50 (4
[CapitalDelta]c [CapitalDelta]p-4
[CapitalDelta]p^2+[CapitalOmega]c^2)) (([Gamma]a+2 [ImaginaryI] (
[CapitalDelta]c-[CapitalDelta]p-[CapitalDelta]s)) ([Gamma]a+2
[ImaginaryI] ([CapitalDelta]c+[CapitalDelta]d-[CapitalDelta]p-
[CapitalDelta]s))+[CapitalOmega]d^2)))/.{[Gamma]a[Rule]1,
dephasing[Rule]10^-4,[CapitalDelta]a->-1(-[CapitalDelta]p +
[CapitalDelta]c-[CapitalDelta]s+[CapitalDelta]d)};*)
(*
MaxFINDER[psC[Chi]2]
MaxFINDER[level5[Chi]2]*)
MaxFINDER[pCs[Chi]2fun[.1]]
Plot[Evaluate[MaxFINDER[pCs[Chi]2fun[bacon]]], {bacon, .1, 2},
 PlotPoints -> 15]