- Float_t kK1=0.28278795,kK2=0.96242952, kSqrtK3 =0.77459667, kD=0.445;
- Float_t lambda = x/kD;
- Float_t a=1-TMath::TanH(kK2*lambda)*TMath::TanH(kK2*lambda);
- Float_t b=1+kSqrtK3*kSqrtK3*TMath::TanH(kK2*lambda)*TMath::TanH(kK2*lambda);
- Float_t mathi = kK1*a/b;
+
+ Double_t lambda = x/AliHMPIDParam::PitchAnodeCathode();
+ Double_t tanh = TMath::TanH(AliHMPIDParam::K2x()*lambda);
+ Double_t a=1-tanh*tanh;
+ Double_t b=1+AliHMPIDParam::SqrtK3x()*AliHMPIDParam::SqrtK3x()*tanh*tanh;
+ Double_t mathi = AliHMPIDParam::K1x()*a/b;
+ return mathi;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+
+Double_t AliHMPIDDigit::MathiesonY(Double_t y)const
+{
+// Mathieson function.
+// This is the answer to electrostatic problem of charge distrubution in MWPC described elsewhere. (NIM A370(1988)602-603)
+// Arguments: x- position of the center of Mathieson distribution
+// Returns: value of the Mathieson function
+
+ Double_t lambda = y/AliHMPIDParam::PitchAnodeCathode();
+ Double_t tanh = TMath::TanH(AliHMPIDParam::K2y()*lambda);
+ Double_t a=1-tanh*tanh;
+ Double_t b=1+AliHMPIDParam::SqrtK3y()*AliHMPIDParam::SqrtK3y()*tanh*tanh;
+ Double_t mathi = AliHMPIDParam::K1y()*a/b;