1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 Revision 1.5 2000/11/21 13:47:55 gosset
19 All Mathieson parameters (Sqrt(K3), K2 and K4) set in one function,
20 SetSqrtKx3AndDeriveKx2Kx4 or SetSqrtKx3AndDeriveKx2Kx4,
21 for each cathode plane
23 Revision 1.4 2000/10/25 10:41:52 morsch
24 IntPH(..): Protec Log against random numbers equal to 0.
26 Revision 1.3 2000/07/03 11:54:57 morsch
27 AliMUONSegmentation and AliMUONHitMap have been replaced by AliSegmentation and AliHitMap in STEER
28 The methods GetPadIxy and GetPadXxy of AliMUONSegmentation have changed name to GetPadI and GetPadC.
30 Revision 1.2 2000/06/15 07:58:48 morsch
31 Code from MUON-dev joined
33 Revision 1.1.2.1 2000/06/09 21:33:35 morsch
34 AliMUONResponse code from AliMUONSegResV0.cxx
38 #include "AliMUONResponseV0.h"
39 #include "AliSegmentation.h"
44 ClassImp(AliMUONResponseV0)
46 //__________________________________________________________________________
47 void AliMUONResponseV0::SetSqrtKx3AndDeriveKx2Kx4(Float_t SqrtKx3)
49 // Set to "SqrtKx3" the Mathieson parameter K3 ("fSqrtKx3")
50 // in the X direction, perpendicular to the wires,
51 // and derive the Mathieson parameters K2 ("fKx2") and K4 ("fKx4")
52 // in the same direction
54 fKx2 = TMath::Pi() / 2. * (1. - 0.5 * fSqrtKx3);
55 Float_t cx1 = fKx2 * fSqrtKx3 / 4. / TMath::ATan(Double_t(fSqrtKx3));
56 fKx4 = cx1 / fKx2 / fSqrtKx3;
59 //__________________________________________________________________________
60 void AliMUONResponseV0::SetSqrtKy3AndDeriveKy2Ky4(Float_t SqrtKy3)
62 // Set to "SqrtKy3" the Mathieson parameter K3 ("fSqrtKy3")
63 // in the Y direction, along the wires,
64 // and derive the Mathieson parameters K2 ("fKy2") and K4 ("fKy4")
65 // in the same direction
67 fKy2 = TMath::Pi() / 2. * (1. - 0.5 * fSqrtKy3);
68 Float_t cy1 = fKy2 * fSqrtKy3 / 4. / TMath::ATan(Double_t(fSqrtKy3));
69 fKy4 = cy1 / fKy2 / fSqrtKy3;
72 Float_t AliMUONResponseV0::IntPH(Float_t eloss)
74 // Calculate charge from given ionization energy loss
76 nel= Int_t(eloss*1.e9/32.);
79 for (Int_t i=1;i<=nel;i++) {
81 while(!arg) arg = gRandom->Rndm();
82 charge -= fChargeSlope*TMath::Log(arg);
86 // -------------------------------------------
88 Float_t AliMUONResponseV0::IntXY(AliSegmentation * segmentation)
90 // Calculate charge on current pad according to Mathieson distribution
92 const Float_t kInversePitch = 1/fPitch;
94 // Integration limits defined by segmentation model
96 Float_t xi1, xi2, yi1, yi2;
97 segmentation->IntegrationLimits(xi1,xi2,yi1,yi2);
98 xi1=xi1*kInversePitch;
99 xi2=xi2*kInversePitch;
100 yi1=yi1*kInversePitch;
101 yi2=yi2*kInversePitch;
103 // The Mathieson function
104 Double_t ux1=fSqrtKx3*TMath::TanH(fKx2*xi1);
105 Double_t ux2=fSqrtKx3*TMath::TanH(fKx2*xi2);
107 Double_t uy1=fSqrtKy3*TMath::TanH(fKy2*yi1);
108 Double_t uy2=fSqrtKy3*TMath::TanH(fKy2*yi2);
111 return Float_t(4.*fKx4*(TMath::ATan(ux2)-TMath::ATan(ux1))*
112 fKy4*(TMath::ATan(uy2)-TMath::ATan(uy1)));
115 Int_t AliMUONResponseV0::DigitResponse(Int_t digit)
117 // add white noise and do zero-suppression and signal truncation
118 // Float_t meanNoise = gRandom->Gaus(1, 0.2);
119 // correct noise for slat chambers;
120 // one more field to add to AliMUONResponseV0 to allow different noises ????
121 Float_t meanNoise = gRandom->Gaus(1.5, 0.2);
122 Float_t noise = gRandom->Gaus(0, meanNoise);
124 if ( digit <= ZeroSuppression()) digit = 0;
125 if ( digit > MaxAdc()) digit=MaxAdc();