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 **************************************************************************/
16 ////////////////////////////////////////////////////////////
17 // Factory for muon chambers, segmentations and response //
18 ////////////////////////////////////////////////////////////
22 #include "AliMUONFactory.h"
24 #include "AliMUONChamber.h"
25 #include "AliMUONResponseV0.h"
26 #include "AliMUONResponseTrigger.h"
27 #include "AliMUONSegmentationV01.h"
28 #include "AliMUONSegmentationV02.h"
29 #include "AliMUONSegmentationSlat.h"
30 #include "AliMUONSegmentationSlatN.h"
31 #include "AliMUONSegmentationTriggerX.h"
32 #include "AliMUONSegmentationTriggerY.h"
34 ClassImp(AliMUONFactory)
36 //__________________________________________________________________________
37 AliMUONFactory::AliMUONFactory()
45 //__________________________________________________________________________
46 AliMUONFactory::AliMUONFactory(const AliMUONFactory& rhs)
49 // Protected copy constructor
51 Fatal("AliMUONFactory", "Not implemented.");
54 //__________________________________________________________________________
56 AliMUONFactory::~AliMUONFactory()
61 //__________________________________________________________________________
62 AliMUONFactory& AliMUONFactory::operator=(const AliMUONFactory& rhs)
64 // Protected assignement operator
66 if (this == &rhs) return *this;
68 Fatal("operator=", "Not implemented.");
73 //__________________________________________________________________________
74 void AliMUONFactory::BuildCommon()
77 // Construct the default response.
80 // Default response: 5 mm of gas
81 fResponse0 = new AliMUONResponseV0;
82 fResponse0->SetSqrtKx3AndDeriveKx2Kx4(0.7131); // sqrt(0.5085)
83 fResponse0->SetSqrtKy3AndDeriveKy2Ky4(0.7642); // sqrt(0.5840)
84 fResponse0->SetPitch(0.25); // anode-cathode distance
85 fResponse0->SetSigmaIntegration(10.);
86 fResponse0->SetChargeSlope(10);
87 fResponse0->SetChargeSpread(0.18, 0.18);
88 fResponse0->SetMaxAdc(4096);
89 fResponse0->SetSaturation(3000);
90 fResponse0->SetZeroSuppression(6);
93 //__________________________________________________________________________
94 void AliMUONFactory::BuildStation1()
96 //--------------------------------------------------------
97 // Configuration for Chamber TC1/2 (Station 1) ----------
98 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
100 // Response for 4 mm of gas (station 1)
101 // automatic consistency with width of sensitive medium in CreateGeometry ????
102 AliMUONResponseV0* responseSt1 = new AliMUONResponseV0;
103 // Mathieson parameters from L.Kharmandarian's thesis, page 190
104 responseSt1->SetSqrtKx3AndDeriveKx2Kx4(0.7000); // sqrt(0.4900)
105 responseSt1->SetSqrtKy3AndDeriveKy2Ky4(0.7550); // sqrt(0.5700)
106 responseSt1->SetPitch(0.20); // anode-cathode distance
107 responseSt1->SetSigmaIntegration(10.);
108 // ChargeSlope larger to compensate for the smaller anode-cathode distance
109 // and keep the same most probable ADC channel for mip's
110 responseSt1->SetChargeSlope(62.5);
111 // assumed proportionality to anode-cathode distance for ChargeSpread
112 responseSt1->SetChargeSpread(0.144, 0.144);
113 responseSt1->SetMaxAdc(4096);
114 responseSt1->SetSaturation(3000);
115 responseSt1->SetZeroSuppression(6);
117 //--------------------------------------------------------
118 // Configuration for Chamber TC1/2 (Station 1) ----------
120 Float_t rseg1[4]={17.5, 55.2, 71.3, 95.5};
121 Int_t nseg1[4]={4, 4, 2, 1};
125 fMUON->SetNsec(chamber-1,2);
127 AliMUONSegmentationV01 *seg11=new AliMUONSegmentationV01(4);
129 seg11->SetSegRadii(rseg1);
130 seg11->SetPadSize(2.4, 0.4); // smaller pad size
131 seg11->SetDAnod(0.20); // smaller distance between anode wires
132 seg11->SetPadDivision(nseg1);
134 fMUON->SetSegmentationModel(chamber-1, 1, seg11);
136 AliMUONSegmentationV02 *seg12=new AliMUONSegmentationV02(4);
137 seg12->SetSegRadii(rseg1);
138 seg12->SetPadSize(0.6, 1.6); // smaller pad size
139 seg12->SetDAnod(0.20); // smaller distance between anode wires
140 seg12->SetPadDivision(nseg1);
142 fMUON->SetSegmentationModel(chamber-1, 2, seg12);
144 fMUON->SetResponseModel(chamber-1, responseSt1); // special response
145 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
150 fMUON->SetNsec(chamber-1,2);
152 AliMUONSegmentationV01 *seg21=new AliMUONSegmentationV01(4);
153 seg21->SetSegRadii(rseg1);
154 seg21->SetPadSize(2.4, 0.4); // smaller pad size
155 seg21->SetDAnod(0.20); // smaller distance between anode wires
156 seg21->SetPadDivision(nseg1);
157 fMUON->SetSegmentationModel(chamber-1, 1, seg21);
159 AliMUONSegmentationV02 *seg22=new AliMUONSegmentationV02(4);
160 seg22->SetSegRadii(rseg1);
161 seg22->SetPadSize(0.6, 1.6); // smaller pad size
162 seg22->SetDAnod(0.20); // smaller distance between anode wires
163 seg22->SetPadDivision(nseg1);
164 fMUON->SetSegmentationModel(chamber-1, 2, seg22);
166 fMUON->SetResponseModel(chamber-1, responseSt1); // special response
167 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
171 //__________________________________________________________________________
172 void AliMUONFactory::BuildStation2()
175 //--------------------------------------------------------
176 // Configuration for Chamber TC3/4 (Station 2) -----------
177 ///^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
178 // Float_t rseg2[4]={23.5, 87.7, 122.4, 122.5};
180 Float_t rseg2[4]={23.5, 53.5, 90.5, 122.5};
181 Int_t nseg2[4]={4, 4, 2, 1};
185 fMUON->SetNsec(chamber-1,2);
187 AliMUONSegmentationV01 *seg31=new AliMUONSegmentationV01(4);
188 seg31->SetSegRadii(rseg2);
189 seg31->SetPadSize(3.0, 0.5);
190 seg31->SetDAnod(3.0/3./4);
191 seg31->SetPadDivision(nseg2);
192 fMUON->SetSegmentationModel(chamber-1, 1, seg31);
194 AliMUONSegmentationV02 *seg32=new AliMUONSegmentationV02(4);
195 seg32->SetSegRadii(rseg2);
196 seg32->SetPadSize(0.75, 2.0);
197 seg32->SetPadDivision(nseg2);
198 seg32->SetDAnod(3.0/3./4);
200 fMUON->SetSegmentationModel(chamber-1, 2, seg32);
202 fMUON->SetResponseModel(chamber-1, fResponse0);
203 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
208 fMUON->SetNsec(chamber-1,2);
210 AliMUONSegmentationV01 *seg41=new AliMUONSegmentationV01(4);
211 seg41->SetSegRadii(rseg2);
212 seg41->SetPadSize(3.0, 0.5);
213 seg41->SetDAnod(3.0/3./4);
214 seg41->SetPadDivision(nseg2);
215 fMUON->SetSegmentationModel(chamber-1, 1, seg41);
217 AliMUONSegmentationV02 *seg42=new AliMUONSegmentationV02(4);
218 seg42->SetSegRadii(rseg2);
219 seg42->SetPadSize(0.75, 2.0);
220 seg42->SetPadDivision(nseg2);
221 seg42->SetDAnod(3.0/3./4);
223 fMUON->SetSegmentationModel(chamber-1, 2, seg42);
225 fMUON->SetResponseModel(chamber-1, fResponse0);
226 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
230 //__________________________________________________________________________
231 void AliMUONFactory::BuildStation3()
233 //--------------------------------------------------------
234 // Configuration for Chamber TC5/6 (Station 3) ----------
235 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
237 Int_t nseg3[4]={4, 4, 2, 1};
238 Int_t npcb5[36] = {0,0,2,0,
248 Float_t shift = 1.5/2.;
249 Float_t xpos5[9] = {2., 2., 2., 2.,33., 2., 2., 2., 2.};
250 Float_t ypos5 = -(20.+4.*(40.-2.*shift));
253 fMUON->SetNsec(chamber-1,2);
254 AliMUONSegmentationSlat *seg51=new AliMUONSegmentationSlat(4);
256 seg51->SetShift(shift);
257 seg51->SetNPCBperSector(npcb5);
258 seg51->SetSlatXPositions(xpos5);
259 seg51->SetSlatYPosition(ypos5);
260 seg51->SetPadSize(10.,0.5);
261 seg51->SetDAnod(0.25);
262 seg51->SetPadDivision(nseg3);
263 fMUON->SetSegmentationModel(chamber-1, 1, seg51);
265 AliMUONSegmentationSlatN *seg52=new AliMUONSegmentationSlatN(4);
267 seg52->SetShift(shift);
268 seg52->SetNPCBperSector(npcb5);
269 seg52->SetSlatXPositions(xpos5);
270 seg52->SetSlatYPosition(ypos5);
271 seg52->SetPadSize(1., 10.); // DeltaX(non bending) = 2 * DeltaY(bending)
272 seg52->SetDAnod(0.25);
273 seg52->SetPadDivision(nseg3);
274 fMUON->SetSegmentationModel(chamber-1, 2, seg52);
275 fMUON->SetResponseModel(chamber-1, fResponse0);
276 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
279 fMUON->SetNsec(chamber-1,2);
280 AliMUONSegmentationSlat *seg61=new AliMUONSegmentationSlat(4);
282 seg61->SetShift(shift);
283 seg61->SetNPCBperSector(npcb5);
284 seg61->SetSlatXPositions(xpos5);
285 seg61->SetSlatYPosition(ypos5);
286 seg61->SetPadSize(10.,0.5);
287 seg61->SetDAnod(0.25);
288 seg61->SetPadDivision(nseg3);
289 fMUON->SetSegmentationModel(chamber-1, 1, seg61);
291 AliMUONSegmentationSlatN *seg62=new AliMUONSegmentationSlatN(4);
293 seg62->SetShift(shift);
294 seg62->SetNPCBperSector(npcb5);
295 seg62->SetSlatXPositions(xpos5);
296 seg62->SetSlatYPosition(ypos5);
297 seg62->SetPadSize(1., 10.); // DeltaX(non bending) = 2 * DeltaY(bending)
298 seg62->SetDAnod(0.25);
299 seg62->SetPadDivision(nseg3);
300 fMUON->SetSegmentationModel(chamber-1, 2, seg62);
301 fMUON->SetResponseModel(chamber-1, fResponse0);
302 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
306 //__________________________________________________________________________
307 void AliMUONFactory::BuildStation4()
309 //--------------------------------------------------------
310 // Configuration for Chamber TC7/8 (Station 4) ----------
311 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
313 Int_t nseg4[4]={4, 4, 2, 1};
318 fMUON->SetNsec(chamber-1,2);
320 AliMUONSegmentationSlat *seg71=new AliMUONSegmentationSlat(4);
321 Float_t shift = 1.5/2.;
322 Int_t npcb7[44] = {0,0,0,3,
333 Float_t xpos7[11] = {2., 2., 2., 2., 2., 40.5, 2., 2., 2., 2., 2.};
334 Float_t ypos7 = -(20.+5.*(40.-2.*shift));
336 seg71->SetNSlats(11);
337 seg71->SetShift(shift);
338 seg71->SetNPCBperSector(npcb7);
339 seg71->SetSlatXPositions(xpos7);
340 seg71->SetSlatYPosition(ypos7);
342 seg71->SetPadSize(10.,0.5);
343 seg71->SetDAnod(0.25);
344 seg71->SetPadDivision(nseg4);
345 fMUON->SetSegmentationModel(chamber-1, 1, seg71);
347 AliMUONSegmentationSlatN *seg72=new AliMUONSegmentationSlatN(4);
349 fMUON->SetSegmentationModel(chamber-1, 2, seg72);
350 seg72->SetNSlats(11);
351 seg72->SetShift(shift);
352 seg72->SetNPCBperSector(npcb7);
353 seg72->SetSlatXPositions(xpos7);
354 seg72->SetSlatYPosition(ypos7);
355 seg72->SetPadSize(1., 10.); // DeltaX(non bending) = 2 * DeltaY(bending)
356 seg72->SetDAnod(0.25);
357 seg72->SetPadDivision(nseg4);
359 fMUON->SetResponseModel(chamber-1, fResponse0);
360 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
364 fMUON->SetNsec(chamber-1,2);
366 AliMUONSegmentationSlat *seg81=new AliMUONSegmentationSlat(4);
368 seg81->SetNSlats(11);
369 seg81->SetShift(shift);
370 seg81->SetNPCBperSector(npcb7);
371 seg81->SetSlatXPositions(xpos7);
372 seg81->SetSlatYPosition(ypos7);
373 seg81->SetPadSize(10.,0.5);
374 seg81->SetDAnod(0.25);
375 seg81->SetPadDivision(nseg4);
376 fMUON->SetSegmentationModel(chamber-1, 1, seg81);
378 AliMUONSegmentationSlat *seg82=new AliMUONSegmentationSlatN(4);
380 fMUON->SetSegmentationModel(chamber-1, 2, seg82);
381 seg82->SetNSlats(11);
382 seg82->SetShift(shift);
383 seg82->SetNPCBperSector(npcb7);
384 seg82->SetSlatXPositions(xpos7);
385 seg82->SetSlatYPosition(ypos7);
386 seg82->SetPadSize(1., 10.); // DeltaX(non bending) = 2 * DeltaY(bending)
387 seg82->SetDAnod(0.25);
388 seg82->SetPadDivision(nseg4);
390 fMUON->SetResponseModel(chamber-1, fResponse0);
391 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
394 //__________________________________________________________________________
395 void AliMUONFactory::BuildStation5()
397 //--------------------------------------------------------
398 // Configuration for Chamber TC9/10 (Station 5) ---------
399 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
404 fMUON->SetNsec(chamber-1,2);
406 AliMUONSegmentationSlat *seg91=new AliMUONSegmentationSlat(4);
408 Int_t nseg4[4]={4, 4, 2, 1};
409 Float_t shift = 1.5/2.;
410 Int_t npcb9[52] = {0,0,0,3,
424 Float_t xpos9[13] = {2., 2., 2., 2., 2., 2., 40.5, 2., 2., 2., 2., 2., 2.};
425 Float_t ypos9 = -(20.+6.*(40.-2.*shift));
427 seg91->SetNSlats(13);
428 seg91->SetShift(shift);
429 seg91->SetNPCBperSector(npcb9);
430 seg91->SetSlatXPositions(xpos9);
431 seg91->SetSlatYPosition(ypos9);
432 seg91->SetPadSize(10.,0.5);
433 seg91->SetDAnod(0.25);
434 seg91->SetPadDivision(nseg4);
435 fMUON->SetSegmentationModel(chamber-1, 1, seg91);
437 AliMUONSegmentationSlatN *seg92=new AliMUONSegmentationSlatN(4);
439 fMUON->SetSegmentationModel(chamber-1, 2, seg92);
440 seg92->SetNSlats(13);
441 seg92->SetShift(shift);
442 seg92->SetNPCBperSector(npcb9);
443 seg92->SetSlatXPositions(xpos9);
444 seg92->SetSlatYPosition(ypos9);
445 seg92->SetPadSize(1., 10.); // DeltaX(non bending) = 2 * DeltaY(bending)
446 seg92->SetDAnod(0.25);
447 seg92->SetPadDivision(nseg4);
449 fMUON->SetResponseModel(chamber-1, fResponse0);
450 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
454 fMUON->SetNsec(chamber-1,2);
456 AliMUONSegmentationSlat *seg101=new AliMUONSegmentationSlat(4);
458 seg101->SetNSlats(13);
459 seg101->SetShift(shift);
460 seg101->SetNPCBperSector(npcb9);
461 seg101->SetSlatXPositions(xpos9);
462 seg101->SetSlatYPosition(ypos9);
463 seg101->SetPadSize(10.,0.5);
464 seg101->SetDAnod(0.25);
465 seg101->SetPadDivision(nseg4);
466 fMUON->SetSegmentationModel(chamber-1, 1, seg101);
468 AliMUONSegmentationSlatN *seg102=new AliMUONSegmentationSlatN(4);
470 fMUON->SetSegmentationModel(chamber-1, 2, seg102);
471 seg102->SetNSlats(13);
472 seg102->SetShift(shift);
473 seg102->SetNPCBperSector(npcb9);
474 seg102->SetSlatXPositions(xpos9);
475 seg102->SetSlatYPosition(ypos9);
476 seg102->SetPadSize(1., 10.); // DeltaX(non bending) = 2 * DeltaY(bending)
477 seg102->SetDAnod(0.25);
478 seg102->SetPadDivision(nseg4);
480 fMUON->SetResponseModel(chamber-1, fResponse0);
481 fMUON->Chamber(chamber-1).SetChargeCorrel(0.11); // 11% charge spread
484 //__________________________________________________________________________
485 void AliMUONFactory::BuildStation6()
487 //--------------------------------------------------------
488 // Configuration for Trigger Stations --------------------
489 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
492 AliMUONResponseTrigger* responseTrigger0 = new AliMUONResponseTrigger;
494 // AliMUONResponseTriggerV1* responseTrigger0 = new AliMUONResponseTriggerV1;
497 fMUON->SetNsec(chamber-1,2);
498 AliMUONSegmentationTriggerX *seg111=new AliMUONSegmentationTriggerX;
499 fMUON->SetSegmentationModel(chamber-1, 1, seg111);
500 AliMUONSegmentationTriggerY *seg112=new AliMUONSegmentationTriggerY;
501 fMUON->SetSegmentationModel(chamber-1, 2, seg112);
503 fMUON->SetResponseModel(chamber-1, responseTrigger0);
504 fMUON->Chamber(chamber-1).SetChargeCorrel(0); // same charge on cathodes
508 fMUON->SetNsec(chamber-1,2);
509 AliMUONSegmentationTriggerX *seg121=new AliMUONSegmentationTriggerX;
510 fMUON->SetSegmentationModel(chamber-1, 1, seg121);
511 AliMUONSegmentationTriggerY *seg122=new AliMUONSegmentationTriggerY;
512 fMUON->SetSegmentationModel(chamber-1, 2, seg122);
514 fMUON->SetResponseModel(chamber-1, responseTrigger0);
515 fMUON->Chamber(chamber-1).SetChargeCorrel(0); // same charge on cathodes
518 fMUON->SetNsec(chamber-1,2);
519 AliMUONSegmentationTriggerX *seg131=new AliMUONSegmentationTriggerX;
520 fMUON->SetSegmentationModel(chamber-1, 1, seg131);
521 AliMUONSegmentationTriggerY *seg132=new AliMUONSegmentationTriggerY;
522 fMUON->SetSegmentationModel(chamber-1, 2, seg132);
523 fMUON->SetResponseModel(chamber-1, responseTrigger0);
524 fMUON->Chamber(chamber-1).SetChargeCorrel(0); // same charge on cathodes
527 fMUON->SetNsec(chamber-1,2);
528 AliMUONSegmentationTriggerX *seg141=new AliMUONSegmentationTriggerX;
529 fMUON->SetSegmentationModel(chamber-1, 1, seg141);
530 AliMUONSegmentationTriggerY *seg142=new AliMUONSegmentationTriggerY;
531 fMUON->SetSegmentationModel(chamber-1, 2, seg142);
533 fMUON->SetResponseModel(chamber-1, responseTrigger0);
534 fMUON->Chamber(chamber-1).SetChargeCorrel(0); // same charge on cathodes
537 //__________________________________________________________________________
538 void AliMUONFactory::Build(AliMUON* where, const char* what)
541 // Construct MUON from chambers, segmentation and responses
548 if (strcmp(tmp, "default")==0) {
549 // Set default parameters
551 fMUON->SetMaxStepGas(0.1);
552 fMUON->SetMaxStepAlu(0.1);
554 // Build all stations
564 if(fMUON->GetDebug()) {
565 printf("\nAliMUONFactory: --------AliMUONFactory------------------------------");
566 printf("\nAliMUONFactory: Non default version of MUON selected ");
567 printf("\nAliMUONFactory: You have to construct yourself the MUON elements !!");
568 printf("\nAliMUONFactory: ----------------------------------------------------");
573 //__________________________________________________________________________
574 void AliMUONFactory::BuildStation(AliMUON* where, Int_t stationNumber)
577 // Construct MUON from chambers, segmentation and responses
581 // First define the number of planes that are segmented (1 or 2) by a call
583 // Then chose for each chamber (chamber plane) the segmentation
584 // and response model.
585 // They should be equal for the two chambers of each station. In a future
586 // version this will be enforced.
590 if (!fResponse0) BuildCommon();
592 switch (stationNumber) {
593 case 1: BuildStation1(); break;
594 case 2: BuildStation2(); break;
595 case 3: BuildStation3(); break;
596 case 4: BuildStation4(); break;
597 case 5: BuildStation5(); break;
598 case 6: BuildStation6(); break;
600 default: Fatal("Build", "Wrong station number");