(* Here I look at the 1 parameter family and lambda = 1.  I make plots for the paper for T = 20, 40, 60 *)

(* I'll use the parameter lambda = (kp^2+m^2)/(2*|eE|)  *)

omplus = Sqrt[(kz+eE*T)^2 + 2*eE*lambda]

omminus = Sqrt[(kz-eE*T)^2 + 2*eE*lambda]

(* This is for the 1-parameter not the 2-parameter case.  See (4.6).   *)

B2 = (Cosh[(Pi/2)*Sqrt[(2*eE*T^2)^2-1]]^2+Sinh[(Pi/2)*(omplus-omminus)*T]^2)/(Sinh[Pi*omplus*T]*Sinh[Pi*omminus*T])

(* Let's fix eE = 1 without loss of generality  *)

eE = 1

(* For the plots we'll now fix lambda = 1.  We can change to other values later if necessary. *)

B21 = B2/.{lambda->1}

A = Pi*T*(omplus+omminus)

B = A - 2*Pi*T^2

T1 = Table[{N[i/50],N[B2/.{lambda->1,T->20,kz->i/50}]},{i,-1000,1000}];

T2 = Table[{N[i/50],N[B2/.{lambda->1,T->40,kz->i/50}]},{i,-1000,1000}];

T3 = Table[{N[i/50],N[B2/.{lambda->1,T->60,kz->i/50}]},{i,-2000,2000}];

ListPlot[{T1,T2,T3},PlotRange->{{-10,10},{0,0.0025}}]

(* I already have T = 20 from a previous version of the program *)

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_1_parameter_model\\B2_1D_T_40.dat",T2]

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_1_parameter_model\\B2_1D_T_60.dat",T3]