- Double_t energy;
- Double_t dpos;
- Double_t dnde;
-
- Int_t pos1, pos2;
-
- const Int_t nV = 31;
-
- Float_t vxe[nV] = { 2.3026, 2.9957, 3.4012, 3.6889, 3.9120
- , 4.0943, 4.2485, 4.3820, 4.4998, 4.6052
- , 4.7005, 5.0752, 5.2983, 5.7038, 5.9915
- , 6.2146, 6.5221, 6.9078, 7.3132, 7.6009
- , 8.0064, 8.5172, 8.6995, 8.9872, 9.2103
- , 9.4727, 9.9035,10.3735,10.5966,10.8198
- ,11.5129 };
-
- Float_t vye[nV] = { 80.0 , 31.0 , 23.3 , 21.1 , 21.0
- , 20.9 , 20.8 , 20.0 , 16.0 , 11.0
- , 8.0 , 6.0 , 5.2 , 4.6 , 4.0
- , 3.5 , 3.0 , 1.4 , 0.67 , 0.44
- , 0.3 , 0.18 , 0.12 , 0.08 , 0.056
- , 0.04 , 0.023, 0.015, 0.011, 0.01
- , 0.004 };
-
- energy = x[0];
-
- // Find the position
- pos1 = pos2 = 0;
- dpos = 0;
- do {
- dpos = energy - vxe[pos2++];
- }
- while (dpos > 0);
- pos2--;
- if (pos2 > nV) pos2 = nV;
- pos1 = pos2 - 1;
-
- // Differentiate between the sampling points
- dnde = (vye[pos1] - vye[pos2]) / (vxe[pos2] - vxe[pos1]);
-
- return dnde;