



double LZEff(double E)  
{
   double Egrid[200]=  {2.500e-01,7.500e-01,1.250e+00,1.750e+00,2.250e+00,2.750e+00,3.250e+00,3.750e+00,4.250e+00,4.750e+00,5.250e+00,5.750e+00,6.250e+00,6.750e+00,7.250e+00,7.750e+00,8.250e+00,8.750e+00,9.250e+00,9.750e+00,1.025e+01,1.075e+01,1.125e+01,1.175e+01,1.225e+01,1.275e+01,1.325e+01,1.375e+01,1.425e+01,1.475e+01,1.525e+01,1.575e+01,1.625e+01,1.675e+01,1.725e+01,1.775e+01,1.825e+01,1.875e+01,1.925e+01,1.975e+01,2.025e+01,2.075e+01,2.125e+01,2.175e+01,2.225e+01,2.275e+01,2.325e+01,2.375e+01,2.425e+01,2.475e+01,2.525e+01,2.575e+01,2.625e+01,2.675e+01,2.725e+01,2.775e+01,2.825e+01,2.875e+01,2.925e+01,2.975e+01,3.025e+01,3.075e+01,3.125e+01,3.175e+01,3.225e+01,3.275e+01,3.325e+01,3.375e+01,3.425e+01,3.475e+01,3.525e+01,3.575e+01,3.625e+01,3.675e+01,3.725e+01,3.775e+01,3.825e+01,3.875e+01,3.925e+01,3.975e+01,4.025e+01,4.075e+01,4.125e+01,4.175e+01,4.225e+01,4.275e+01,4.325e+01,4.375e+01,4.425e+01,4.475e+01,4.525e+01,4.575e+01,4.625e+01,4.675e+01,4.725e+01,4.775e+01,4.825e+01,4.875e+01,4.925e+01,4.975e+01,5.025e+01,5.075e+01,5.125e+01,5.175e+01,5.225e+01,5.275e+01,5.325e+01,5.375e+01,5.425e+01,5.475e+01,5.525e+01,5.575e+01,5.625e+01,5.675e+01,5.725e+01,5.775e+01,5.825e+01,5.875e+01,5.925e+01,5.975e+01,6.025e+01,6.075e+01,6.125e+01,6.175e+01,6.225e+01,6.275e+01,6.325e+01,6.375e+01,6.425e+01,6.475e+01,6.525e+01,6.575e+01,6.625e+01,6.675e+01,6.725e+01,6.775e+01,6.825e+01,6.875e+01,6.925e+01,6.975e+01,7.025e+01,7.075e+01,7.125e+01,7.175e+01,7.225e+01,7.275e+01,7.325e+01,7.375e+01,7.425e+01,7.475e+01,7.525e+01,7.575e+01,7.625e+01,7.675e+01,7.725e+01,7.775e+01,7.825e+01,7.875e+01,7.925e+01,7.975e+01,8.025e+01,8.075e+01,8.125e+01,8.175e+01,8.225e+01,8.275e+01,8.325e+01,8.375e+01,8.425e+01,8.475e+01,8.525e+01,8.575e+01,8.625e+01,8.675e+01,8.725e+01,8.775e+01,8.825e+01,8.875e+01,8.925e+01,8.975e+01,9.025e+01,9.075e+01,9.125e+01,9.175e+01,9.225e+01,9.275e+01,9.325e+01,9.375e+01,9.425e+01,9.475e+01,9.525e+01,9.575e+01,9.625e+01,9.675e+01,9.725e+01,9.775e+01,9.825e+01,9.875e+01,9.925e+01,9.975e+01};
   double   Eff[200]={0.000e+00,0.000e+00,1.647e-04,3.641e-03,2.408e-02,7.633e-02,1.509e-01,2.416e-01,3.377e-01,4.239e-01,5.062e-01,5.753e-01,6.401e-01,6.913e-01,7.328e-01,7.682e-01,7.929e-01,8.118e-01,8.269e-01,8.377e-01,8.455e-01,8.516e-01,8.566e-01,8.616e-01,8.652e-01,8.683e-01,8.706e-01,8.729e-01,8.751e-01,8.770e-01,8.790e-01,8.806e-01,8.822e-01,8.837e-01,8.852e-01,8.866e-01,8.878e-01,8.891e-01,8.902e-01,8.913e-01,8.922e-01,8.931e-01,8.940e-01,8.948e-01,8.955e-01,8.961e-01,8.967e-01,8.973e-01,8.978e-01,8.983e-01,8.987e-01,8.991e-01,8.995e-01,8.999e-01,9.002e-01,9.006e-01,9.009e-01,9.012e-01,9.015e-01,9.018e-01,9.020e-01,9.023e-01,9.026e-01,9.028e-01,9.031e-01,9.033e-01,9.035e-01,9.037e-01,9.040e-01,9.042e-01,9.044e-01,9.046e-01,9.047e-01,9.049e-01,9.049e-01,9.052e-01,9.049e-01,9.047e-01,9.049e-01,9.046e-01,9.044e-01,9.035e-01,9.021e-01,9.003e-01,8.993e-01,8.960e-01,8.923e-01,8.891e-01,8.826e-01,8.759e-01,8.701e-01,8.586e-01,8.466e-01,8.360e-01,8.196e-01,8.005e-01,7.842e-01,7.631e-01,7.353e-01,7.123e-01,6.866e-01,6.509e-01,6.188e-01,5.875e-01,5.576e-01,5.213e-01,4.813e-01,4.547e-01,4.123e-01,3.818e-01,3.464e-01,3.156e-01,2.807e-01,2.494e-01,2.228e-01,2.000e-01,1.765e-01,1.520e-01,1.368e-01,1.183e-01,9.838e-02,8.288e-02,7.037e-02,6.165e-02,5.069e-02,4.270e-02,3.620e-02,2.972e-02,2.453e-02,1.967e-02,1.581e-02,1.143e-02,9.893e-03,8.031e-03,6.912e-03,4.918e-03,3.794e-03,2.988e-03,2.572e-03,1.608e-03,1.268e-03,1.484e-03,8.772e-04,4.203e-04,5.077e-04,2.819e-04,2.222e-04,1.406e-04,2.513e-04,1.121e-04,5.593e-05,5.565e-05,0.000e+00,2.821e-05,2.810e-05,5.627e-05,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00,0.000e+00};
   if(E<Egrid[0] || E>Eff[199]) return 0;
   return polint3(E,200,Egrid,Eff);
}