Make Digits() and SDigits() methods consts.
[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"
f29ba3e1 23#include "AliLog.h"
8c343c7c 24
d5bfadcc 25ClassImp(AliMUONResponseV0)
26
30178c30 27//__________________________________________________________________________
28AliMUONResponseV0::AliMUONResponseV0()
a713db22 29 : AliMUONResponse()
30178c30 30{
31// Default constructor
32
a713db22 33 fMathieson = new AliMUONMathieson();
30178c30 34 fChargeCorrel = 0;
35}
f29ba3e1 36
37 //_________________________________________________________________________
38AliMUONResponseV0::AliMUONResponseV0(const AliMUONResponseV0& rhs)
39 : AliMUONResponse(rhs)
40{
41// Protected copy constructor
42
43 AliFatal("Not implemented.");
44}
45
a713db22 46 //__________________________________________________________________________
47AliMUONResponseV0::~AliMUONResponseV0()
48{
49 delete fMathieson;
50}
f29ba3e1 51
52 //________________________________________________________________________
53AliMUONResponseV0& AliMUONResponseV0::operator = (const AliMUONResponseV0& rhs)
54{
55// Protected assignement operator
56
57 if (this == &rhs) return *this;
58
59 AliFatal("Not implemented.");
60
61 return *this;
62}
63
d5bfadcc 64 //__________________________________________________________________________
65void AliMUONResponseV0::SetSqrtKx3AndDeriveKx2Kx4(Float_t SqrtKx3)
66{
67 // Set to "SqrtKx3" the Mathieson parameter K3 ("fSqrtKx3")
68 // in the X direction, perpendicular to the wires,
69 // and derive the Mathieson parameters K2 ("fKx2") and K4 ("fKx4")
70 // in the same direction
a713db22 71 fMathieson->SetSqrtKx3AndDeriveKx2Kx4(SqrtKx3);
d5bfadcc 72}
73
74 //__________________________________________________________________________
75void AliMUONResponseV0::SetSqrtKy3AndDeriveKy2Ky4(Float_t SqrtKy3)
76{
77 // Set to "SqrtKy3" the Mathieson parameter K3 ("fSqrtKy3")
78 // in the Y direction, along the wires,
79 // and derive the Mathieson parameters K2 ("fKy2") and K4 ("fKy4")
80 // in the same direction
a713db22 81 fMathieson->SetSqrtKy3AndDeriveKy2Ky4(SqrtKy3);
d5bfadcc 82}
a713db22 83 //__________________________________________________________________________
a9e2aefa 84Float_t AliMUONResponseV0::IntPH(Float_t eloss)
85{
86 // Calculate charge from given ionization energy loss
87 Int_t nel;
4ac9d21e 88 nel= Int_t(eloss*1.e9/27.4);
a9e2aefa 89 Float_t charge=0;
90 if (nel == 0) nel=1;
91 for (Int_t i=1;i<=nel;i++) {
01997fa2 92 Float_t arg=0.;
93 while(!arg) arg = gRandom->Rndm();
94 charge -= fChargeSlope*TMath::Log(arg);
a9e2aefa 95 }
96 return charge;
97}
a713db22 98
a713db22 99 //-------------------------------------------
100Float_t AliMUONResponseV0::IntXY(Int_t idDE, AliMUONGeometrySegmentation* segmentation)
101{
102 // Calculate charge on current pad according to Mathieson distribution
a9e2aefa 103
a713db22 104 return fMathieson->IntXY(idDE, segmentation);
105}
106 //-------------------------------------------
f7b62f08 107Int_t AliMUONResponseV0::DigitResponse(Int_t digit, AliMUONTransientDigit* /*where*/)
a9e2aefa 108{
109 // add white noise and do zero-suppression and signal truncation
b64652f5 110// Float_t meanNoise = gRandom->Gaus(1, 0.2);
111 // correct noise for slat chambers;
112 // one more field to add to AliMUONResponseV0 to allow different noises ????
4ac9d21e 113 Float_t meanNoise = gRandom->Gaus(1., 0.2);
114 Float_t noise = gRandom->Gaus(0., meanNoise);
104b5ac2 115 digit += TMath::Nint(noise);
ddc10e24 116 if ( digit <= ZeroSuppression()) digit = 0;
4ac9d21e 117 // if ( digit > MaxAdc()) digit=MaxAdc();
118 if ( digit > Saturation()) digit=Saturation();
119
a9e2aefa 120 return digit;
121}
122
123
124
125
126
127
128
129
130