]>
Commit | Line | Data |
---|---|---|
61b8353b | 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 | ||
16 | //////////////////////////////////////////////////////////// | |
17 | // Factory for muon response // | |
18 | //////////////////////////////////////////////////////////// | |
19 | ||
20 | /* $Id$ */ | |
21 | ||
22 | #include "AliMUONResponseFactory.h" | |
23 | #include "AliRun.h" | |
24 | #include "AliLog.h" | |
25 | ||
26 | #include "AliMpPlaneType.h" | |
27 | ||
28 | #include "AliMUON.h" | |
29 | #include "AliMUONConstants.h" | |
30 | #include "AliMUONTriggerConstants.h" | |
31 | #include "AliMUONChamber.h" | |
32 | #include "AliMUONResponseV0.h" | |
33 | #include "AliMUONResponseTrigger.h" | |
34 | ||
35 | ClassImp(AliMUONResponseFactory) | |
36 | ||
37 | //__________________________________________________________________________ | |
38 | AliMUONResponseFactory::AliMUONResponseFactory(const char* name) | |
39 | : TNamed(name, ""), | |
40 | fMUON(0), | |
41 | fResponse0(0) | |
42 | { | |
43 | AliDebug(1,Form("ctor this=%p",this)); | |
44 | } | |
45 | ||
46 | //__________________________________________________________________________ | |
47 | AliMUONResponseFactory::AliMUONResponseFactory() | |
48 | : TNamed(), | |
49 | fMUON(0), | |
50 | fResponse0(0) | |
51 | { | |
52 | AliDebug(1,Form("default (empty) ctor this=%p",this)); | |
53 | // Default constructor | |
54 | } | |
55 | ||
56 | //__________________________________________________________________________ | |
57 | AliMUONResponseFactory::AliMUONResponseFactory(const AliMUONResponseFactory& rhs) | |
58 | : TNamed(rhs) | |
59 | { | |
60 | // Protected copy constructor | |
61 | ||
62 | AliFatal("Not implemented."); | |
63 | } | |
64 | ||
65 | //__________________________________________________________________________ | |
66 | ||
67 | AliMUONResponseFactory::~AliMUONResponseFactory() | |
68 | { | |
69 | // Destructor | |
70 | AliDebug(1,Form("dtor this=%p",this)); | |
71 | } | |
72 | ||
73 | //__________________________________________________________________________ | |
74 | AliMUONResponseFactory& AliMUONResponseFactory::operator=(const AliMUONResponseFactory& rhs) | |
75 | { | |
76 | // Protected assignement operator | |
77 | ||
78 | if (this == &rhs) return *this; | |
79 | ||
80 | AliFatal("Not implemented."); | |
81 | ||
82 | return *this; | |
83 | } | |
84 | ||
85 | //__________________________________________________________________________ | |
86 | void AliMUONResponseFactory::BuildCommon() | |
87 | { | |
88 | // | |
89 | // Construct the default response. | |
90 | // | |
91 | ||
92 | // Default response: 5 mm of gas | |
93 | fResponse0 = new AliMUONResponseV0; | |
94 | fResponse0->SetSqrtKx3AndDeriveKx2Kx4(0.7131); // sqrt(0.5085) | |
95 | fResponse0->SetSqrtKy3AndDeriveKy2Ky4(0.7642); // sqrt(0.5840) | |
96 | fResponse0->SetPitch(AliMUONConstants::Pitch()); // anode-cathode distance | |
97 | fResponse0->SetSigmaIntegration(10.); | |
98 | fResponse0->SetChargeSlope(10); | |
99 | fResponse0->SetChargeSpread(0.18, 0.18); | |
100 | fResponse0->SetMaxAdc(4096); | |
101 | fResponse0->SetSaturation(3000); | |
102 | fResponse0->SetZeroSuppression(6); | |
103 | } | |
104 | ||
105 | //__________________________________________________________________________ | |
106 | void AliMUONResponseFactory::BuildStation1() | |
107 | { | |
108 | /// Configuration for Chamber TC1/2 (Station 1) ---------- | |
109 | ||
110 | // Response for 4 mm of gas (station 1) | |
111 | // automatic consistency with width of sensitive medium in CreateGeometry ???? | |
112 | AliMUONResponseV0* responseSt1 = new AliMUONResponseV0; | |
113 | // Mathieson parameters from L.Kharmandarian's thesis, page 190 | |
114 | responseSt1->SetSqrtKx3AndDeriveKx2Kx4(0.7000); // sqrt(0.4900) | |
115 | responseSt1->SetSqrtKy3AndDeriveKy2Ky4(0.7550); // sqrt(0.5700) | |
116 | responseSt1->SetPitch(AliMUONConstants::PitchSt1()); // anode-cathode distance | |
117 | responseSt1->SetSigmaIntegration(10.); | |
118 | // ChargeSlope larger to compensate for the smaller anode-cathode distance | |
119 | // and keep the same most probable ADC channel for mip's | |
120 | responseSt1->SetChargeSlope(62.5); | |
121 | // assumed proportionality to anode-cathode distance for ChargeSpread | |
122 | responseSt1->SetChargeSpread(0.144, 0.144); | |
123 | responseSt1->SetMaxAdc(4096); | |
124 | responseSt1->SetSaturation(3000); | |
125 | responseSt1->SetZeroSuppression(6); | |
126 | ||
127 | for (Int_t chamber = 0; chamber < 2; chamber++) { | |
128 | fMUON->SetResponseModel(chamber, responseSt1); // special response | |
129 | fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread | |
130 | } | |
131 | } | |
132 | ||
133 | //__________________________________________________________________________ | |
134 | void AliMUONResponseFactory::BuildStation2() | |
135 | { | |
136 | /// Configuration for Chamber TC3/4 (Station 2) ----------- | |
137 | ||
138 | for (Int_t chamber = 2; chamber < 4; chamber++) { | |
139 | fMUON->SetResponseModel(chamber, fResponse0); // normal response | |
140 | fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread | |
141 | } | |
142 | } | |
143 | ||
144 | //__________________________________________________________________________ | |
145 | void AliMUONResponseFactory::BuildStation3() | |
146 | { | |
147 | /// Configuration for Chamber TC5/6 (Station 3) ---------- | |
148 | ||
149 | for (Int_t chamber = 4; chamber < 6; chamber++) { | |
150 | fMUON->SetResponseModel(chamber, fResponse0); // normal response | |
151 | fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread | |
152 | } | |
153 | } | |
154 | ||
155 | //__________________________________________________________________________ | |
156 | void AliMUONResponseFactory::BuildStation4() | |
157 | { | |
158 | /// Configuration for Chamber TC7/8 (Station 4) ---------- | |
159 | ||
160 | for (Int_t chamber = 6; chamber < 8; chamber++) { | |
161 | fMUON->SetResponseModel(chamber, fResponse0); // normal response | |
162 | fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread | |
163 | } | |
164 | } | |
165 | ||
166 | //__________________________________________________________________________ | |
167 | void AliMUONResponseFactory::BuildStation5() | |
168 | { | |
169 | /// Configuration for Chamber TC9/10 (Station 5) --------- | |
170 | ||
171 | for (Int_t chamber = 8; chamber < 10; chamber++) { | |
172 | fMUON->SetResponseModel(chamber, fResponse0); // normal response | |
173 | fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread | |
174 | } | |
175 | } | |
176 | ||
177 | //__________________________________________________________________________ | |
178 | void AliMUONResponseFactory::BuildStation6() | |
179 | { | |
180 | /// Configuration for Trigger Chambers (Station 6,7) --------- | |
181 | ||
182 | for (Int_t chamber = 10; chamber < 14; chamber++) { | |
183 | ||
184 | fMUON->SetResponseModel(chamber, new AliMUONResponseTrigger); | |
185 | fMUON->Chamber(chamber).SetChargeCorrel(0); // same charge on both cathodes | |
186 | } | |
187 | } | |
188 | ||
189 | //__________________________________________________________________________ | |
190 | void AliMUONResponseFactory::Build(AliMUON* where) | |
191 | { | |
192 | // Construct MUON responses | |
193 | ||
194 | fMUON = where; | |
195 | ||
196 | // Set default parameters | |
197 | fMUON->SetIshunt(0); | |
198 | fMUON->SetMaxStepGas(0.1); | |
199 | fMUON->SetMaxStepAlu(0.1); | |
200 | ||
201 | // Build stations | |
202 | BuildCommon(); | |
203 | BuildStation1(); | |
204 | BuildStation2(); | |
205 | BuildStation3(); | |
206 | BuildStation4(); | |
207 | BuildStation5(); | |
208 | BuildStation6(); | |
209 | } | |
210 | ||
211 | //__________________________________________________________________________ | |
212 | void AliMUONResponseFactory::BuildStation(AliMUON* where, Int_t stationNumber) | |
213 | { | |
214 | // Construct MUON responses | |
215 | ||
216 | fMUON = where; | |
217 | if (!fResponse0) BuildCommon(); | |
218 | ||
219 | switch (stationNumber) { | |
220 | case 1: BuildStation1(); break; | |
221 | case 2: BuildStation2(); break; | |
222 | case 3: BuildStation3(); break; | |
223 | case 4: BuildStation4(); break; | |
224 | case 5: BuildStation5(); break; | |
225 | case 6: BuildStation6(); break; | |
226 | ||
227 | default: AliFatal("Wrong station number"); | |
228 | } | |
229 | } |