+ static const Float_t lSy[6][24] = {
+ {75.7561, 0.0325, 0.0175, 0.0174, 0.0206, 0.0232,
+ 0.0253, 0.0262, 0.0265, 0.0264, 0.0266, 0.0257,
+ 0.0258, 0.0261, 0.0259, 0.0253, 0.0257, 0.0261,
+ 0.0255, 0.0250, 0.0259, 0.0266, 0.0278, 0.0319
+ },
+ {49.2252, 0.0371, 0.0204, 0.0189, 0.0230, 0.0261,
+ 0.0281, 0.0290, 0.0292, 0.0286, 0.0277, 0.0279,
+ 0.0285, 0.0281, 0.0291, 0.0281, 0.0281, 0.0282,
+ 0.0272, 0.0282, 0.0282, 0.0284, 0.0310, 0.0334
+ },
+ {55.1674, 0.0388, 0.0212, 0.0200, 0.0239, 0.0271,
+ 0.0288, 0.0299, 0.0306, 0.0300, 0.0296, 0.0303,
+ 0.0293, 0.0290, 0.0291, 0.0294, 0.0295, 0.0290,
+ 0.0293, 0.0292, 0.0292, 0.0293, 0.0316, 0.0358
+ },
+ {45.1004, 0.0411, 0.0225, 0.0215, 0.0249, 0.0281,
+ 0.0301, 0.0315, 0.0320, 0.0308, 0.0318, 0.0321,
+ 0.0312, 0.0311, 0.0316, 0.0315, 0.0310, 0.0308,
+ 0.0313, 0.0303, 0.0314, 0.0314, 0.0324, 0.0369
+ },
+ {43.8614, 0.0420, 0.0239, 0.0224, 0.0268, 0.0296,
+ 0.0322, 0.0336, 0.0333, 0.0326, 0.0321, 0.0325,
+ 0.0329, 0.0326, 0.0323, 0.0322, 0.0326, 0.0320,
+ 0.0329, 0.0319, 0.0314, 0.0329, 0.0341, 0.0373
+ },
+ {40.5440, 0.0434, 0.0246, 0.0236, 0.0275, 0.0311,
+ 0.0332, 0.0345, 0.0347, 0.0347, 0.0340, 0.0336,
+ 0.0339, 0.0344, 0.0339, 0.0341, 0.0341, 0.0342,
+ 0.0345, 0.0328, 0.0341, 0.0332, 0.0356, 0.0398
+ },
+ };
+ // adjusted ...
+ return TMath::Max(lSy[ly][tb]-0.0150, 0.0010);
+
+/* const Double_t sy[24][10]={