Quantum efficiency of 10-inch Hamamatsu PMT.
819 tuple(
const double __QE,
820 const double __l_gel,
821 const double __l_glass) :
832 static const tuple ntuple[] = {
833 tuple(0.000e-2, 0.00, 0.00),
834 tuple(1.988e-2, 100.81, 148.37),
835 tuple(2.714e-2, 99.94, 142.87),
836 tuple(3.496e-2, 99.89, 135.64),
837 tuple(4.347e-2, 96.90, 134.58),
838 tuple(5.166e-2, 96.42, 138.27),
839 tuple(6.004e-2, 94.36, 142.40),
840 tuple(6.885e-2, 89.09, 147.16),
841 tuple(8.105e-2, 90.10, 151.80),
842 tuple(10.13e-2, 86.95, 150.88),
843 tuple(13.03e-2, 85.88, 145.68),
844 tuple(15.29e-2, 84.49, 139.70),
845 tuple(16.37e-2, 81.08, 126.55),
846 tuple(17.11e-2, 78.18, 118.86),
847 tuple(17.86e-2, 76.48, 113.90),
848 tuple(18.95e-2, 74.55, 116.08),
849 tuple(20.22e-2, 72.31, 109.23),
850 tuple(21.26e-2, 68.05, 81.63),
851 tuple(22.10e-2, 66.91, 65.66),
852 tuple(22.65e-2, 64.48, 77.30),
853 tuple(23.07e-2, 62.53, 73.02),
854 tuple(23.14e-2, 59.38, 81.25),
855 tuple(23.34e-2, 56.64, 128.04),
856 tuple(22.95e-2, 53.29, 61.84),
857 tuple(22.95e-2, 48.96, 19.23),
858 tuple(22.74e-2, 45.71, 27.21),
859 tuple(23.48e-2, 41.88, 18.09),
860 tuple(22.59e-2, 37.14, 8.41),
861 tuple(20.61e-2, 30.49, 3.92),
862 tuple(17.68e-2, 23.08, 1.82),
863 tuple(13.18e-2, 15.60, 0.84),
864 tuple(7.443e-2, 8.00, 0.39),
865 tuple(2.526e-2, 0.00, 0.17),
866 tuple(0.000e-2, 0.00, 0.00)
869 static const double cola = 0.9;
870 static const double x_glass = 1.5;
871 static const double x_gel = 1.0;
874 static const int N =
sizeof(ntuple) /
sizeof(ntuple[0]) - 1;
876 static const double xmax = 620.0;
877 static const double xmin = 290.0;
879 const double x = lambda;
883 if (x > xmin && x < xmax) {
885 const int i = (int) (N * (x - xmax) / (xmin - xmax));
886 const int j = (i == N ? i - 1 : i + 1);
888 const double x1 = xmax + i * (xmin - xmax) / N;
889 const double x2 = xmax + j * (xmin - xmax) / N;
891 const double dx = (x - x1) / (x1 - x2);
893 const double QE = ntuple[i].QE + (ntuple[i].QE - ntuple[
j].QE ) * dx;
894 const double l_gel = ntuple[i].l_gel + (ntuple[i].l_gel - ntuple[
j].l_gel ) * dx;
895 const double l_glass = ntuple[i].l_glass + (ntuple[i].l_glass - ntuple[
j].l_glass) * dx;
901 if (l_glass > 0.0 && l_gel > 0.0)
902 y *=
exp(-x_glass/l_glass) *
exp(-x_gel/l_gel);
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` typeset -Z 4 STRING JOpera1D -f hydrophone.root
then usage $script[input file[working directory[option]]] nWhere option can be N