Added <assert.h> include
[u/mrichter/AliRoot.git] / MUON / AliMUONResponseV0.cxx
CommitLineData
a9e2aefa 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
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 **************************************************************************/
15
88cb7938 16/* $Id$ */
a9e2aefa 17
a9e2aefa 18#include <TMath.h>
19#include <TRandom.h>
20
30178c30 21#include "AliMUONResponseV0.h"
a713db22 22#include "AliMUONGeometrySegmentation.h"
8c343c7c 23
d5bfadcc 24ClassImp(AliMUONResponseV0)
25
30178c30 26//__________________________________________________________________________
27AliMUONResponseV0::AliMUONResponseV0()
a713db22 28 : AliMUONResponse()
30178c30 29{
30// Default constructor
31
a713db22 32 fMathieson = new AliMUONMathieson();
30178c30 33 fChargeCorrel = 0;
34}
a713db22 35 //__________________________________________________________________________
36AliMUONResponseV0::~AliMUONResponseV0()
37{
38 delete fMathieson;
39}
d5bfadcc 40 //__________________________________________________________________________
41void AliMUONResponseV0::SetSqrtKx3AndDeriveKx2Kx4(Float_t SqrtKx3)
42{
43 // Set to "SqrtKx3" the Mathieson parameter K3 ("fSqrtKx3")
44 // in the X direction, perpendicular to the wires,
45 // and derive the Mathieson parameters K2 ("fKx2") and K4 ("fKx4")
46 // in the same direction
a713db22 47 fMathieson->SetSqrtKx3AndDeriveKx2Kx4(SqrtKx3);
d5bfadcc 48}
49
50 //__________________________________________________________________________
51void AliMUONResponseV0::SetSqrtKy3AndDeriveKy2Ky4(Float_t SqrtKy3)
52{
53 // Set to "SqrtKy3" the Mathieson parameter K3 ("fSqrtKy3")
54 // in the Y direction, along the wires,
55 // and derive the Mathieson parameters K2 ("fKy2") and K4 ("fKy4")
56 // in the same direction
a713db22 57 fMathieson->SetSqrtKy3AndDeriveKy2Ky4(SqrtKy3);
d5bfadcc 58}
a713db22 59 //__________________________________________________________________________
a9e2aefa 60Float_t AliMUONResponseV0::IntPH(Float_t eloss)
61{
62 // Calculate charge from given ionization energy loss
63 Int_t nel;
4ac9d21e 64 nel= Int_t(eloss*1.e9/27.4);
a9e2aefa 65 Float_t charge=0;
66 if (nel == 0) nel=1;
67 for (Int_t i=1;i<=nel;i++) {
01997fa2 68 Float_t arg=0.;
69 while(!arg) arg = gRandom->Rndm();
70 charge -= fChargeSlope*TMath::Log(arg);
a9e2aefa 71 }
72 return charge;
73}
a713db22 74
a713db22 75 //-------------------------------------------
76Float_t AliMUONResponseV0::IntXY(Int_t idDE, AliMUONGeometrySegmentation* segmentation)
77{
78 // Calculate charge on current pad according to Mathieson distribution
a9e2aefa 79
a713db22 80 return fMathieson->IntXY(idDE, segmentation);
81}
82 //-------------------------------------------
f7b62f08 83Int_t AliMUONResponseV0::DigitResponse(Int_t digit, AliMUONTransientDigit* /*where*/)
a9e2aefa 84{
85 // add white noise and do zero-suppression and signal truncation
b64652f5 86// Float_t meanNoise = gRandom->Gaus(1, 0.2);
87 // correct noise for slat chambers;
88 // one more field to add to AliMUONResponseV0 to allow different noises ????
4ac9d21e 89 Float_t meanNoise = gRandom->Gaus(1., 0.2);
90 Float_t noise = gRandom->Gaus(0., meanNoise);
104b5ac2 91 digit += TMath::Nint(noise);
ddc10e24 92 if ( digit <= ZeroSuppression()) digit = 0;
4ac9d21e 93 // if ( digit > MaxAdc()) digit=MaxAdc();
94 if ( digit > Saturation()) digit=Saturation();
95
a9e2aefa 96 return digit;
97}
98
99
100
101
102
103
104
105
106