MUON/AliMUONFactoryV2.cxx

   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 chambers, segmentations and response //
  18 ////////////////////////////////////////////////////////////
  19
  20 /* $Id$ */
  21
  22 #include "AliMUONFactoryV2.h"
  23 #include "AliRun.h"
  24 #include "AliLog.h"
  25
  26 #include "AliMUON.h"
  27 #include "AliMUONConstants.h"
  28 #include "AliMUONTriggerConstants.h"
  29 #include "AliMUONChamber.h"
  30 #include "AliMUONResponseV0.h"
  31 #include "AliMUONGeometryModule.h"
  32 #include "AliMUONGeometryStore.h"
  33 #include "AliMUONGeometrySegmentation.h"
  34 #include "AliMUONVGeometryDEIndexing.h"
  35 #include "AliMUONSt12QuadrantSegmentation.h"
  36 #include "AliMUONSt345SlatSegmentation.h"
  37 #include "AliMUONTriggerSegmentation.h"
  38 #include "AliMUONResponseTrigger.h"
  39
  40 ClassImp(AliMUONFactoryV2)
  41
  42 //__________________________________________________________________________
  43   AliMUONFactoryV2::AliMUONFactoryV2(const char* name)
  44     : TNamed(name, ""),
  45       fMUON(0),
  46       fResponse0(0),
  47       fDESegmentations(0)
  48 {
  49   // FactoryV2 inherite from Factory for switching in AliMUONv1::Init()
  50   // to be changed when old segmentation will be removed.
  51
  52   fDESegmentations = new TObjArray();
  53   fDESegmentations->SetOwner(kTRUE);
  54 }
  55
  56 //__________________________________________________________________________
  57   AliMUONFactoryV2::AliMUONFactoryV2()
  58     : TNamed(),
  59       fMUON(0),
  60       fResponse0(0),
  61       fDESegmentations(0)
  62 {
  63 // Default constructor
  64 }
  65
  66 //__________________________________________________________________________
  67 AliMUONFactoryV2::AliMUONFactoryV2(const AliMUONFactoryV2& rhs)
  68  : TNamed(rhs)
  69 {
  70   // Protected copy constructor
  71
  72   AliFatal("Not implemented.");
  73 }
  74
  75 //__________________________________________________________________________
  76
  77 AliMUONFactoryV2::~AliMUONFactoryV2()
  78 {
  79 // Destructor
  80
  81   delete fDESegmentations;
  82 }
  83
  84 //__________________________________________________________________________
  85 AliMUONFactoryV2&  AliMUONFactoryV2::operator=(const AliMUONFactoryV2& rhs)
  86 {
  87   // Protected assignement operator
  88
  89   if (this == &rhs) return *this;
  90
  91   AliFatal("Not implemented.");
  92
  93   return *this;
  94 }
  95
  96 //__________________________________________________________________________
  97 Bool_t AliMUONFactoryV2::IsGeometryDefined(Int_t ichamber)
  98 {
  99 // Return true, if det elements for the chamber with the given ichamber Id
 100 // are defined in geometry (the geometry builder for this chamber was activated)
 101
 102   if ( ! fMUON ||
 103        ! fMUON->Chamber(ichamber).GetGeometry() ||
 104        ! fMUON->Chamber(ichamber).GetGeometry()->GetDEIndexing() ||
 105        ! fMUON->Chamber(ichamber).GetGeometry()->GetDEIndexing()->GetNofDetElements() )
 106
 107     return kFALSE;
 108
 109   return kTRUE;
 110 }
 111
 112 //__________________________________________________________________________
 113 void AliMUONFactoryV2::BuildCommon()
 114 {
 115   //
 116   // Construct the default response.
 117   //
 118
 119   // Default response: 5 mm of gas
 120   fResponse0 = new AliMUONResponseV0;
 121   fResponse0->SetSqrtKx3AndDeriveKx2Kx4(0.7131); // sqrt(0.5085)
 122   fResponse0->SetSqrtKy3AndDeriveKy2Ky4(0.7642); // sqrt(0.5840)
 123   fResponse0->SetPitch(AliMUONConstants::Pitch()); // anode-cathode distance
 124   fResponse0->SetSigmaIntegration(10.);
 125   fResponse0->SetChargeSlope(10);
 126   fResponse0->SetChargeSpread(0.18, 0.18);
 127   fResponse0->SetMaxAdc(4096);
 128   fResponse0->SetSaturation(3000);
 129   fResponse0->SetZeroSuppression(6);
 130 }
 131
 132 //__________________________________________________________________________
 133 void AliMUONFactoryV2::BuildStation1()
 134 {
 135   //--------------------------------------------------------
 136   // Configuration for Chamber TC1/2  (Station 1) ----------
 137   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 138
 139
 140   // Response for 4 mm of gas (station 1)
 141   // automatic consistency with width of sensitive medium in CreateGeometry ????
 142   AliMUONResponseV0* responseSt1 = new AliMUONResponseV0;
 143   // Mathieson parameters from L.Kharmandarian's thesis, page 190
 144   responseSt1->SetSqrtKx3AndDeriveKx2Kx4(0.7000); // sqrt(0.4900)
 145   responseSt1->SetSqrtKy3AndDeriveKy2Ky4(0.7550); // sqrt(0.5700)
 146   responseSt1->SetPitch(AliMUONConstants::PitchSt1()); // anode-cathode distance
 147   responseSt1->SetSigmaIntegration(10.);
 148   // ChargeSlope larger to compensate for the smaller anode-cathode distance
 149   // and keep the same most probable ADC channel for mip's
 150   responseSt1->SetChargeSlope(62.5);
 151   // assumed proportionality to anode-cathode distance for ChargeSpread
 152   responseSt1->SetChargeSpread(0.144, 0.144);
 153   responseSt1->SetMaxAdc(4096);
 154   responseSt1->SetSaturation(3000);
 155   responseSt1->SetZeroSuppression(6);
 156
 157   // Quadrant segmentations:
 158   AliMUONSt12QuadrantSegmentation* bendSt1
 159     = new AliMUONSt12QuadrantSegmentation(kStation1, kBendingPlane);
 160   AliMUONSt12QuadrantSegmentation* nonbendSt1
 161     = new AliMUONSt12QuadrantSegmentation(kStation1, kNonBendingPlane);
 162
 163   // Add in the array (for safe deleting)
 164   fDESegmentations->Add(bendSt1);
 165   fDESegmentations->Add(nonbendSt1);
 166
 167   AliMUONGeometrySegmentation* segmentation[2];
 168
 169   for (Int_t chamber = 0; chamber < 2; chamber++) {
 170
 171     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 172     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 173
 174     // id detection elt for chamber 1
 175     Int_t id0 = (chamber+1)*100;
 176
 177     //--------------------------------------------------------
 178     // Configuration for Chamber TC1/2  (Station 1) ----------
 179
 180
 181     fMUON->SetNsec(chamber,2);
 182
 183     // cathode 0
 184     segmentation[0]->Add(id0,      bendSt1);
 185     segmentation[0]->Add(id0 +  1, nonbendSt1);
 186     segmentation[0]->Add(id0 + 50, bendSt1);
 187     segmentation[0]->Add(id0 + 51, nonbendSt1);
 188     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
 189
 190     // cathode 1
 191     segmentation[1]->Add(id0,      nonbendSt1);
 192     segmentation[1]->Add(id0 +  1, bendSt1);
 193     segmentation[1]->Add(id0 + 50, nonbendSt1);
 194     segmentation[1]->Add(id0 + 51, bendSt1);
 195     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
 196
 197     fMUON->SetResponseModel(chamber, responseSt1); // special response
 198     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
 199
 200   }
 201 }
 202
 203 //__________________________________________________________________________
 204 void AliMUONFactoryV2::BuildStation2()
 205 {
 206   //
 207   //--------------------------------------------------------
 208   // Configuration for Chamber TC3/4 (Station 2) -----------
 209   ///^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 210
 211
 212   // Quadrant segmentations:
 213   AliMUONSt12QuadrantSegmentation* bendSt2
 214     = new AliMUONSt12QuadrantSegmentation(kStation2, kBendingPlane);
 215   AliMUONSt12QuadrantSegmentation* nonbendSt2
 216     = new AliMUONSt12QuadrantSegmentation(kStation2, kNonBendingPlane);
 217
 218   // Add in the array (for safe deleting)
 219   fDESegmentations->Add(bendSt2);
 220   fDESegmentations->Add(nonbendSt2);
 221
 222   AliMUONGeometrySegmentation* segmentation[2];
 223
 224   for (Int_t chamber = 2; chamber < 4; chamber++) {
 225
 226     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 227     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 228
 229     // id detection elt for chamber 1
 230     Int_t id0 = (chamber+1)*100;
 231
 232     //--------------------------------------------------------
 233     // Configuration for Chamber TC3/4  (Station 2) ----------
 234
 235
 236     fMUON->SetNsec(chamber,2);
 237
 238     // cathode 0
 239     segmentation[0]->Add(id0,      bendSt2);
 240     segmentation[0]->Add(id0 +  1, nonbendSt2);
 241     segmentation[0]->Add(id0 + 50, bendSt2);
 242     segmentation[0]->Add(id0 + 51, nonbendSt2);
 243     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
 244
 245     // cathode 1
 246     segmentation[1]->Add(id0,      nonbendSt2);
 247     segmentation[1]->Add(id0 +  1, bendSt2);
 248     segmentation[1]->Add(id0 + 50, nonbendSt2);
 249     segmentation[1]->Add(id0 + 51, bendSt2);
 250     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
 251
 252     fMUON->SetResponseModel(chamber, fResponse0); // normal response
 253     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
 254
 255   }
 256 }
 257
 258
 259 //__________________________________________________________________________
 260 void AliMUONFactoryV2::BuildStation3()
 261 {
 262   //--------------------------------------------------------
 263   // Configuration for Chamber TC5/6  (Station 3) ----------
 264   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 265
 266   AliMUONGeometrySegmentation* segmentation[2];
 267
 268   //Slats Segmentations
 269   AliMUONSt345SlatSegmentation *slatsegB[4]; // Types of segmentation for St3
 270   AliMUONSt345SlatSegmentation *slatsegNB[4];
 271   // Bending
 272
 273   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones
 274   for(Int_t i=0; i<4; i++) {
 275     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
 276     fDESegmentations->Add(slatsegB[i]);
 277     slatsegB[i]->SetPadSize(10.,0.5);
 278     slatsegB[i]->SetPadDivision(ndiv);
 279     slatsegB[i]->SetId(1); // Id elt ????
 280     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
 281     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
 282     fDESegmentations->Add(slatsegNB[i]);
 283     slatsegNB[i]->SetPadSize(1./1.4,10.); // Nbending
 284     slatsegNB[i]->SetPadDivision(ndiv);
 285     slatsegNB[i]->SetId(1);
 286     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
 287   }
 288
 289   // Type 112200 for 500, 501, 508, 509, 510, 517
 290   // in Ch5 (similar for Ch6) for the futur official numbering
 291   // Type 112200 for 503, 504, 505, 555, 554, 553
 292   // in Ch5 (similar for Ch6) actual numbering in the code to be changed in jan05
 293   Int_t n0[4] = { 0, 2, 2, 0 };
 294   slatsegB[0]->SetPcbBoards(n0);
 295   slatsegB[0]->Init(0);
 296   slatsegNB[0]->SetPcbBoards(n0);
 297   slatsegNB[0]->Init(0);
 298
 299   // Type 122200 for 502, 507, 511, 516 (similar in Ch6)
 300   // for future official numbering of ALICE
 301   // Type 122200 for 502, 506, 556, 552 (similiarin Ch6)
 302   // for actual numbering in muon code to be changed in jan05
 303   Int_t n1[4] = { 0, 1, 3, 0 };
 304   slatsegB[1]->SetPcbBoards(n1);
 305   slatsegB[1]->Init(0);
 306   slatsegNB[1]->SetPcbBoards(n1);
 307   slatsegNB[1]->Init(0);
 308
 309   // Type 222000 for 503, 506, 512, 515 (similar in Ch6)
 310   // for future official numbering of ALICE
 311   // Type 222000 for 501, 507, 557, 551 (similiarin Ch6)
 312   // for actual numbering in muon code to be changed in jan05
 313   Int_t n2[4] = { 0, 0, 3, 0 };
 314   slatsegB[2]->SetPcbBoards(n2);
 315   slatsegB[2]->Init(0);
 316   slatsegNB[2]->SetPcbBoards(n2);
 317   slatsegNB[2]->Init(0);
 318
 319   // Type 220000 for 504, 505, 513, 514 (similar in Ch6)
 320   // for future official numbering of ALICE
 321   // Type 220000 for 500, 508, 558, 550 (similiarin Ch6)
 322   // for actual numbering in muon code to be changed in jan05
 323   Int_t n3[4] = { 0, 0, 2, 0 };
 324   slatsegB[3]->SetPcbBoards(n3);
 325   slatsegB[3]->Init(0);
 326   slatsegNB[3]->SetPcbBoards(n3);
 327   slatsegNB[3]->Init(0);
 328
 329   for (Int_t chamber = 4; chamber < 6; chamber++) {
 330
 331     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 332     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 333
 334     // id detection elt for chamber 1
 335     Int_t id0 = (chamber+1)*100;
 336
 337     //--------------------------------------------------------
 338     // Configuration for Chamber TC3/4  (Station 2) ----------
 339
 340
 341     fMUON->SetNsec(chamber,2);
 342
 343     // cathode 0
 344     // type 220000
 345     segmentation[0]->Add(id0, slatsegB[3]);
 346     segmentation[0]->Add(id0+ 8, slatsegB[3]);
 347     segmentation[0]->Add(id0+50, slatsegB[3]);
 348     segmentation[0]->Add(id0+58, slatsegB[3]);
 349     // type 222000
 350     segmentation[0]->Add(id0+ 1, slatsegB[2]);
 351     segmentation[0]->Add(id0+ 7, slatsegB[2]);
 352     segmentation[0]->Add(id0+51, slatsegB[2]);
 353     segmentation[0]->Add(id0+57, slatsegB[2]);
 354     // type 122200
 355     segmentation[0]->Add(id0+ 2, slatsegB[1]);
 356     segmentation[0]->Add(id0+ 6, slatsegB[1]);
 357     segmentation[0]->Add(id0+52, slatsegB[1]);
 358     segmentation[0]->Add(id0+56, slatsegB[1]);
 359     // type 112200
 360     segmentation[0]->Add(id0+ 3, slatsegB[0]);
 361     segmentation[0]->Add(id0+ 4, slatsegB[0]);
 362     segmentation[0]->Add(id0+ 5, slatsegB[0]);
 363     segmentation[0]->Add(id0+53, slatsegB[0]);
 364     segmentation[0]->Add(id0+54, slatsegB[0]);
 365     segmentation[0]->Add(id0+55, slatsegB[0]);
 366     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
 367
 368     // cathode 1
 369     // type 220000
 370     segmentation[1]->Add(id0, slatsegNB[3]);
 371     segmentation[1]->Add(id0+ 8, slatsegNB[3]);
 372     segmentation[1]->Add(id0+50, slatsegNB[3]);
 373     segmentation[1]->Add(id0+58, slatsegNB[3]);
 374     // type 222000
 375     segmentation[1]->Add(id0+ 1, slatsegNB[2]);
 376     segmentation[1]->Add(id0+ 7, slatsegNB[2]);
 377     segmentation[1]->Add(id0+51, slatsegNB[2]);
 378     segmentation[1]->Add(id0+57, slatsegNB[2]);
 379     // type 122200
 380     segmentation[1]->Add(id0+ 2, slatsegNB[1]);
 381     segmentation[1]->Add(id0+ 6, slatsegNB[1]);
 382     segmentation[1]->Add(id0+52, slatsegNB[1]);
 383     segmentation[1]->Add(id0+56, slatsegNB[1]);
 384     // type 112200
 385     segmentation[1]->Add(id0+ 3, slatsegNB[0]);
 386     segmentation[1]->Add(id0+ 4, slatsegNB[0]);
 387     segmentation[1]->Add(id0+ 5, slatsegNB[0]);
 388     segmentation[1]->Add(id0+53, slatsegNB[0]);
 389     segmentation[1]->Add(id0+54, slatsegNB[0]);
 390     segmentation[1]->Add(id0+55, slatsegNB[0]);
 391     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
 392
 393     fMUON->SetResponseModel(chamber, fResponse0);
 394     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
 395   }
 396 }
 397
 398
 399 //__________________________________________________________________________
 400 void AliMUONFactoryV2::BuildStation4()
 401 {
 402   //--------------------------------------------------------
 403   // Configuration for Chamber TC7/8  (Station 4) ----------
 404   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 405
 406
 407   AliMUONGeometrySegmentation* segmentation[2];
 408
 409   //Slats Segmentations
 410   AliMUONSt345SlatSegmentation *slatsegB[7]; // Types of segmentation for St4
 411   AliMUONSt345SlatSegmentation *slatsegNB[7];
 412   // Bending
 413
 414   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones
 415   for(Int_t i = 0; i < 7; i++) {
 416     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
 417     fDESegmentations->Add(slatsegB[i]);
 418     slatsegB[i]->SetPadSize(10.,0.5);
 419     slatsegB[i]->SetPadDivision(ndiv);
 420     slatsegB[i]->SetId(1);
 421     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
 422     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
 423     fDESegmentations->Add(slatsegNB[i]);
 424     slatsegNB[i]->SetPadSize(1./1.4,10.);
 425     slatsegNB[i]->SetPadDivision(ndiv);
 426     slatsegNB[i]->SetId(1);
 427     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
 428   }
 429
 430   Int_t n4[4] = { 0, 1, 2, 2 };
 431   slatsegB[0]->SetPcbBoards(n4);
 432   slatsegB[0]->Init(0); // 0 detection element id
 433   slatsegNB[0]->SetPcbBoards(n4);
 434   slatsegNB[0]->Init(0); // 0 detection element id
 435
 436   // Type 112233 for 701, 712, 714, 725 in Ch7 (similar for Ch8)
 437   // for the futur official numbering
 438   // Type 112233 for 705, 707, 755, 757 in Ch7 (similar for Ch8)
 439   // actual numbering in the code to be changed in jan05
 440   // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9
 441   // (similar for Ch10) for the futur official numbering
 442   // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9
 443   // (similar for Ch10) actual numbering in the code to be changed in jan05
 444   Int_t n5[4] = { 0, 2, 2, 2 };
 445   slatsegB[1]->SetPcbBoards(n5);
 446   slatsegB[1]->Init(0); // 0 detection element id
 447   slatsegNB[1]->SetPcbBoards(n5);
 448   slatsegNB[1]->Init(0); // 0 detection element id
 449
 450   // Type 112230 for 702, 711, 715, 724 in Ch7 (similar for Ch8)
 451   // for the futur official numbering
 452   // Type 112230 for 704, 708, 754, 758 in Ch7 (similar for Ch8)
 453   // actual numbering in the code to be changed in jan05
 454   Int_t n6[4] = { 0, 2, 2, 1 };
 455   slatsegB[2]->SetPcbBoards(n6);
 456   slatsegB[2]->Init(0); // 0 detection element id
 457   slatsegNB[2]->SetPcbBoards(n6);
 458   slatsegNB[2]->Init(0); // 0 detection element id
 459
 460   // Type 222330 for 703, 710, 716, 723 in Ch7 (similar for Ch8)
 461   // for the futur official numbering
 462   // Type 222330 for 703, 709, 753, 759 in Ch7 (similar for Ch8)
 463   // actual numbering in the code to be changed in jan05
 464   Int_t n7[4] = { 0, 0, 3, 2 };
 465   slatsegB[3]->SetPcbBoards(n7);
 466   slatsegB[3]->Init(0); // 0 detection element id
 467   slatsegNB[3]->SetPcbBoards(n7);
 468   slatsegNB[3]->Init(0); // 0 detection element id
 469
 470   // Type 223300 for 704, 709, 717, 722 in Ch7 (similar for Ch8)
 471   // for the futur official numbering
 472   // Type 223300 for 702, 710, 752, 760 in Ch7 (similar for Ch8)
 473   // actual numbering in the code to be changed in jan05
 474   Int_t n8[4] = { 0, 0, 2, 2 };
 475   slatsegB[4]->SetPcbBoards(n8);
 476   slatsegB[4]->Init(0); // 0 detection element id
 477   slatsegNB[4]->SetPcbBoards(n8);
 478   slatsegNB[4]->Init(0); // 0 detection element id
 479
 480   // Type 333000 for 705, 708, 718, 721 in Ch7 (similar for Ch8)
 481   // for the futur official numbering
 482   // Type 333000 for 701, 711, 751, 761 in Ch7 (similar for Ch8)
 483   // actual numbering in the code to be changed in jan05
 484   // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10)
 485   // for the futur official numbering
 486   // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10)
 487   // actual numbering in the code to be changed in jan05
 488   Int_t n9[4] = { 0, 0, 0, 3 };
 489   slatsegB[5]->SetPcbBoards(n9);
 490   slatsegB[5]->Init(0); // 0 detection element id
 491   slatsegNB[5]->SetPcbBoards(n9);
 492   slatsegNB[5]->Init(0); // 0 detection element id
 493
 494   // Type 330000 for 706, 707, 719, 720 in Ch7 (similar for Ch8)
 495   // for the futur official numbering
 496   // Type 330000 for 700, 712, 750, 762 in Ch7 (similar for Ch8)
 497   // actual numbering in the code to be changed in jan05
 498   Int_t n10[4] = { 0, 0, 0, 2 };
 499   slatsegB[6]->SetPcbBoards(n10);
 500   slatsegB[6]->Init(0); // 0 detection element id
 501   slatsegNB[6]->SetPcbBoards(n10);
 502   slatsegNB[6]->Init(0); // 0 detection element id
 503
 504
 505   for (Int_t chamber = 6; chamber < 8; chamber++) {
 506
 507     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 508     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 509
 510     // id detection elt for chamber 1
 511     Int_t id0 = (chamber+1)*100;
 512
 513     //--------------------------------------------------------
 514     // Configuration for Chamber TC6/7  (Station 4) ----------
 515
 516
 517     fMUON->SetNsec(chamber,2);
 518
 519     // cathode 0
 520     // type 122330
 521     segmentation[0]->Add(id0+ 6, slatsegB[0]);
 522     segmentation[0]->Add(id0+56, slatsegB[0]);
 523
 524     // type 112233
 525     segmentation[0]->Add(id0+ 5, slatsegB[1]);
 526     segmentation[0]->Add(id0+ 7, slatsegB[1]);
 527     segmentation[0]->Add(id0+55, slatsegB[1]);
 528     segmentation[0]->Add(id0+57, slatsegB[1]);
 529
 530     // type 112230
 531     segmentation[0]->Add(id0+ 4, slatsegB[2]);
 532     segmentation[0]->Add(id0+ 8, slatsegB[2]);
 533     segmentation[0]->Add(id0+54, slatsegB[2]);
 534     segmentation[0]->Add(id0+58, slatsegB[2]);
 535
 536     // type 222330
 537     segmentation[0]->Add(id0+ 3, slatsegB[3]);
 538     segmentation[0]->Add(id0+ 9, slatsegB[3]);
 539     segmentation[0]->Add(id0+53, slatsegB[3]);
 540     segmentation[0]->Add(id0+59, slatsegB[3]);
 541
 542     // type 223300
 543     segmentation[0]->Add(id0+ 2, slatsegB[4]);
 544     segmentation[0]->Add(id0+10, slatsegB[4]);
 545     segmentation[0]->Add(id0+52, slatsegB[4]);
 546     segmentation[0]->Add(id0+60, slatsegB[4]);
 547
 548     // type 333000
 549     segmentation[0]->Add(id0+ 1, slatsegB[5]);
 550     segmentation[0]->Add(id0+11, slatsegB[5]);
 551     segmentation[0]->Add(id0+51, slatsegB[5]);
 552     segmentation[0]->Add(id0+61, slatsegB[5]);
 553
 554     // type 330000
 555     segmentation[0]->Add(id0   , slatsegB[6]);
 556     segmentation[0]->Add(id0+12, slatsegB[6]);
 557     segmentation[0]->Add(id0+50, slatsegB[6]);
 558     segmentation[0]->Add(id0+62, slatsegB[6]);
 559     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
 560
 561     // cathode 1
 562     // type 122330
 563     segmentation[1]->Add(id0+ 6, slatsegNB[0]);
 564     segmentation[1]->Add(id0+56, slatsegNB[0]);
 565
 566     // type 112233
 567     segmentation[1]->Add(id0+ 5, slatsegNB[1]);
 568     segmentation[1]->Add(id0+ 7, slatsegNB[1]);
 569     segmentation[1]->Add(id0+55, slatsegNB[1]);
 570     segmentation[1]->Add(id0+57, slatsegNB[1]);
 571
 572     // type 112230
 573     segmentation[1]->Add(id0+ 4, slatsegNB[2]);
 574     segmentation[1]->Add(id0+ 8, slatsegNB[2]);
 575     segmentation[1]->Add(id0+54, slatsegNB[2]);
 576     segmentation[1]->Add(id0+58, slatsegNB[2]);
 577
 578     // type 222330
 579     segmentation[1]->Add(id0+ 3, slatsegNB[3]);
 580     segmentation[1]->Add(id0+ 9, slatsegNB[3]);
 581     segmentation[1]->Add(id0+53, slatsegNB[3]);
 582     segmentation[1]->Add(id0+59, slatsegNB[3]);
 583
 584     // type 223300
 585     segmentation[1]->Add(id0+ 2, slatsegNB[4]);
 586     segmentation[1]->Add(id0+10, slatsegNB[4]);
 587     segmentation[1]->Add(id0+52, slatsegNB[4]);
 588     segmentation[1]->Add(id0+60, slatsegNB[4]);
 589
 590     // type 333000
 591     segmentation[1]->Add(id0+ 1, slatsegNB[5]);
 592     segmentation[1]->Add(id0+11, slatsegNB[5]);
 593     segmentation[1]->Add(id0+51, slatsegNB[5]);
 594     segmentation[1]->Add(id0+61, slatsegNB[5]);
 595
 596     // type 330000
 597     segmentation[1]->Add(id0   , slatsegNB[6]);
 598     segmentation[1]->Add(id0+12, slatsegNB[6]);
 599     segmentation[1]->Add(id0+50, slatsegNB[6]);
 600     segmentation[1]->Add(id0+62, slatsegNB[6]);
 601     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
 602
 603     fMUON->SetResponseModel(chamber, fResponse0);
 604     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
 605   }
 606 }
 607
 608 //__________________________________________________________________________
 609 void AliMUONFactoryV2::BuildStation5()
 610 {
 611   //--------------------------------------------------------
 612   // Configuration for Chamber TC9/10  (Station 5) ---------
 613   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 614
 615   AliMUONGeometrySegmentation* segmentation[2];
 616
 617   //Slats Segmentations
 618   AliMUONSt345SlatSegmentation *slatsegB[6]; // Types of segmentation for St5
 619   AliMUONSt345SlatSegmentation *slatsegNB[6];
 620   // Bending
 621
 622   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones
 623   for(Int_t i = 0; i < 6; i++) {
 624     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
 625     fDESegmentations->Add(slatsegB[i]);
 626     slatsegB[i]->SetPadSize(10.,0.5);
 627     slatsegB[i]->SetPadDivision(ndiv);
 628     slatsegB[i]->SetId(1);
 629     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
 630     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
 631     fDESegmentations->Add(slatsegNB[i]);
 632     slatsegNB[i]->SetPadSize(1./1.4,10.);
 633     slatsegNB[i]->SetPadDivision(ndiv);
 634     slatsegNB[i]->SetId(1);
 635     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
 636   }
 637
 638   // Type 122330 for 900, 913 in Ch9 (similar for Ch10)
 639   // for the futur official numbering
 640   // Type 122330 for 906, 956 in Ch9 (similar for Ch10)
 641   // actual numbering in the code to be changed in jan05
 642   Int_t n4[4] = { 0, 1, 2, 2 };
 643   slatsegB[0]->SetPcbBoards(n4);
 644   slatsegB[0]->Init(0); // 0 detection element id
 645   slatsegNB[0]->SetPcbBoards(n4);
 646   slatsegNB[0]->Init(0); // 0 detection element id
 647
 648   // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9
 649   // (similar for Ch10) for the futur official numbering
 650   // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9
 651   // (similar for Ch10) actual numbering in the code to be changed in jan05
 652   Int_t n5[4] = { 0, 2, 2, 2 };
 653   slatsegB[1]->SetPcbBoards(n5);
 654   slatsegB[1]->Init(0); // 0 detection element id
 655   slatsegNB[1]->SetPcbBoards(n5);
 656   slatsegNB[1]->Init(0); // 0 detection element id
 657
 658   // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10)
 659   // for the futur official numbering
 660   // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10)
 661   // actual numbering in the code to be changed in jan05
 662   Int_t n9[4] = { 0, 0, 0, 3 };
 663   slatsegB[2]->SetPcbBoards(n9);
 664   slatsegB[2]->Init(0); // 0 detection element id
 665   slatsegNB[2]->SetPcbBoards(n9);
 666   slatsegNB[2]->Init(0); // 0 detection element id
 667
 668   // Type 222333 for 903, 910, 916, 923 in Ch9 (similar for Ch10)
 669   // for the futur official numbering
 670   // Type 222333 for 903, 909, 953, 959 in Ch9 (similar for Ch10)
 671   // actual numbering in the code to be changed in jan05
 672   Int_t n11[4] = { 0, 0, 3, 3 };
 673   slatsegB[3]->SetPcbBoards(n11);
 674   slatsegB[3]->Init(0); // 0 detection element id
 675   slatsegNB[3]->SetPcbBoards(n11);
 676   slatsegNB[3]->Init(0); // 0 detection element id
 677
 678   // Type 223330 for 904, 909, 917, 922 in Ch9 (similar for Ch10)
 679   // for the futur official numbering
 680   // Type 223330 for 902, 910, 952, 960 in Ch9 (similar for Ch10)
 681   // actual numbering in the code to be changed in jan05
 682   Int_t n12[4] = { 0, 0, 2, 3 };
 683   slatsegB[4]->SetPcbBoards(n12);
 684   slatsegB[4]->Init(0); // 0 detection element id
 685   slatsegNB[4]->SetPcbBoards(n12);
 686   slatsegNB[4]->Init(0); // 0 detection element id
 687
 688   // Type 333300 for 905, 908, 918, 921 in Ch9 (similar for Ch10)
 689   // for the futur official numbering
 690   // Type 333300 for 901, 911, 951, 961 in Ch9 (similar for Ch10)
 691   // actual numbering in the code to be changed in jan05
 692   Int_t n13[4] = { 0, 0, 0, 4 };
 693   slatsegB[5]->SetPcbBoards(n13);
 694   slatsegB[5]->Init(0); // 0 detection element id
 695   slatsegNB[5]->SetPcbBoards(n13);
 696   slatsegNB[5]->Init(0); // 0 detection element id
 697
 698   for (Int_t chamber = 8; chamber < 10; chamber++) {
 699
 700     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 701     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 702
 703     // id detection elt for chamber 1
 704     Int_t id0 = (chamber+1)*100;
 705
 706     //--------------------------------------------------------
 707     // Configuration for Chamber TC8/9  (Station 5) ----------
 708
 709
 710     fMUON->SetNsec(chamber,2);
 711
 712     // cathode 0
 713     // type 122330
 714     segmentation[0]->Add(id0+ 6, slatsegB[0]);
 715     segmentation[0]->Add(id0+56, slatsegB[0]);
 716
 717     // type 112233
 718     segmentation[0]->Add(id0+ 4, slatsegB[1]);
 719     segmentation[0]->Add(id0+ 5, slatsegB[1]);
 720     segmentation[0]->Add(id0+ 7, slatsegB[1]);
 721     segmentation[0]->Add(id0+ 8, slatsegB[1]);
 722     segmentation[0]->Add(id0+54, slatsegB[1]);
 723     segmentation[0]->Add(id0+55, slatsegB[1]);
 724     segmentation[0]->Add(id0+57, slatsegB[1]);
 725     segmentation[0]->Add(id0+58, slatsegB[1]);
 726
 727     // type 333000
 728     segmentation[0]->Add(id0   , slatsegB[2]);
 729     segmentation[0]->Add(id0+12, slatsegB[2]);
 730     segmentation[0]->Add(id0+50, slatsegB[2]);
 731     segmentation[0]->Add(id0+62, slatsegB[2]);
 732
 733     // type 222333
 734     segmentation[0]->Add(id0+ 3, slatsegB[3]);
 735     segmentation[0]->Add(id0+ 9, slatsegB[3]);
 736     segmentation[0]->Add(id0+53, slatsegB[3]);
 737     segmentation[0]->Add(id0+59, slatsegB[3]);
 738
 739     // type 223330
 740     segmentation[0]->Add(id0+ 2, slatsegB[4]);
 741     segmentation[0]->Add(id0+10, slatsegB[4]);
 742     segmentation[0]->Add(id0+52, slatsegB[4]);
 743     segmentation[0]->Add(id0+60, slatsegB[4]);
 744
 745     // type 333300
 746     segmentation[0]->Add(id0+ 1, slatsegB[5]);
 747     segmentation[0]->Add(id0+11, slatsegB[5]);
 748     segmentation[0]->Add(id0+51, slatsegB[5]);
 749     segmentation[0]->Add(id0+61, slatsegB[5]);
 750     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
 751
 752     // cathode 1
 753     // type 122330
 754     segmentation[1]->Add(id0+ 6, slatsegNB[0]);
 755     segmentation[1]->Add(id0+56, slatsegNB[0]);
 756
 757     // type 112233
 758     segmentation[1]->Add(id0+ 4, slatsegNB[1]);
 759     segmentation[1]->Add(id0+ 5, slatsegNB[1]);
 760     segmentation[1]->Add(id0+ 7, slatsegNB[1]);
 761     segmentation[1]->Add(id0+ 8, slatsegNB[1]);
 762     segmentation[1]->Add(id0+54, slatsegNB[1]);
 763     segmentation[1]->Add(id0+55, slatsegNB[1]);
 764     segmentation[1]->Add(id0+57, slatsegNB[1]);
 765     segmentation[1]->Add(id0+58, slatsegNB[1]);
 766
 767     // type 333000
 768     segmentation[1]->Add(id0   , slatsegNB[2]);
 769     segmentation[1]->Add(id0+12, slatsegNB[2]);
 770     segmentation[1]->Add(id0+50, slatsegNB[2]);
 771     segmentation[1]->Add(id0+62, slatsegNB[2]);
 772
 773     // type 222333
 774     segmentation[1]->Add(id0+ 3, slatsegNB[3]);
 775     segmentation[1]->Add(id0+ 9, slatsegNB[3]);
 776     segmentation[1]->Add(id0+53, slatsegNB[3]);
 777     segmentation[1]->Add(id0+59, slatsegNB[3]);
 778
 779     // type 223330
 780     segmentation[1]->Add(id0+ 2, slatsegNB[4]);
 781     segmentation[1]->Add(id0+10, slatsegNB[4]);
 782     segmentation[1]->Add(id0+52, slatsegNB[4]);
 783     segmentation[1]->Add(id0+60, slatsegNB[4]);
 784
 785     // type 333300
 786     segmentation[1]->Add(id0+ 1, slatsegNB[5]);
 787     segmentation[1]->Add(id0+11, slatsegNB[5]);
 788     segmentation[1]->Add(id0+51, slatsegNB[5]);
 789     segmentation[1]->Add(id0+61, slatsegNB[5]);
 790     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
 791
 792     fMUON->SetResponseModel(chamber, fResponse0);
 793     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
 794   }
 795 }
 796
 797 //__________________________________________________________________________
 798 void AliMUONFactoryV2::BuildStation6()
 799 {
 800  // Create Trigger geometry segmentation for given chamber and cathod
 801
 802
 803     AliMUONGeometrySegmentation *chamberSeg[2];
 804 // Cluster-size off
 805         AliMUONResponseTrigger* responseTrigger0 =  new AliMUONResponseTrigger;
 806 // Cluster-size on
 807 //  AliMUONResponseTriggerV1* responseTrigger0 =  new AliMUONResponseTriggerV1;
 808
 809     for (Int_t chamber = 10; chamber < 14; chamber++) {
 810
 811       //Trigger Segmentation
 812       AliMUONTriggerSegmentation *trigSegX[9];
 813       AliMUONTriggerSegmentation *trigSegY[9];
 814       for(Int_t i=0; i<9; i++) {
 815         trigSegX[i] = new AliMUONTriggerSegmentation(1);
 816         trigSegY[i] = new AliMUONTriggerSegmentation(0);
 817         fDESegmentations->Add(trigSegX[i]);
 818         fDESegmentations->Add(trigSegY[i]);
 819         trigSegX[i]->SetLineNumber(9-i);
 820         trigSegY[i]->SetLineNumber(9-i);
 821       }
 822
 823       AliMUONChamber *iChamber, *iChamber1;
 824       iChamber1 = &fMUON->Chamber(10);
 825       iChamber  = &fMUON->Chamber(chamber);
 826       Float_t zpos1= - iChamber1->Z();
 827       Float_t zpos = - iChamber->Z();
 828       Float_t zRatio = zpos / zpos1;
 829
 830       // init
 831       Float_t stripWidth[3]={0.,0.,0.};     // 1.0625 2.125 4.25
 832       Float_t stripLength[4]={0.,0.,0.,0.}; // 17. 34. 51. 68.
 833       for (Int_t i=0; i<3; i++)
 834         stripWidth[i]=AliMUONTriggerConstants::StripWidth(i)*zRatio;
 835       for (Int_t i=0; i<4; i++)
 836         stripLength[i]=AliMUONTriggerConstants::StripLength(i)*zRatio;
 837       Int_t nStrip[7]={0,0,0,0,0,0,0};
 838       Float_t stripYsize[7]={0.,0.,0.,0.,0.,0.,0.};
 839       Float_t stripXsize[7]={0.,0.,0.,0.,0.,0.,0.};
 840
 841       // chamber 8 0 cathode 0
 842       for (Int_t i=0; i<7; i++) nStrip[i]=16;
 843       for (Int_t i=0; i<7; i++) stripYsize[i]=stripWidth[2];
 844       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
 845       stripXsize[6]=stripLength[2];
 846       trigSegX[8]->Init(0,nStrip,stripYsize,stripXsize,0.);
 847       trigSegX[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
 848
 849       // chamber 8 7 1 0 cathode 1
 850       for (Int_t i=0; i<6; i++) nStrip[i]=8;
 851       nStrip[6]=16;
 852       for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
 853       for (Int_t i=0; i<7; i++) stripXsize[i]=stripWidth[2];
 854       trigSegY[8]->Init(0,nStrip,stripYsize,stripXsize,0.);
 855       trigSegY[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
 856       trigSegY[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
 857       trigSegY[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
 858
 859       // chamber 7 6 2 1 cathode 0
 860       for (Int_t i=0; i<6; i++) nStrip[i]=32;
 861       nStrip[6]=16;
 862       for (Int_t i=0; i<6; i++) stripYsize[i]=stripWidth[1];
 863       stripYsize[6]=stripWidth[2];
 864       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
 865       stripXsize[6]=stripLength[2];
 866       trigSegX[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
 867       trigSegX[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
 868       trigSegX[2]->Init(0,nStrip,stripYsize,stripXsize,0.);
 869       trigSegX[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
 870
 871       // chamber 6 2 cathode 1
 872       for (Int_t i=0; i<5; i++) nStrip[i]=16;
 873       for (Int_t i=5; i<6; i++) nStrip[i]=8;
 874       nStrip[6]=16;
 875       for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
 876       for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
 877       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
 878       trigSegY[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
 879       trigSegY[2]->Init(0,nStrip,stripYsize,stripXsize,0.);
 880
 881       // chamber 5 3 cathode 0
 882       nStrip[0]=48;
 883       for (Int_t i=1; i<3; i++) nStrip[i]=64;
 884       for (Int_t i=3; i<6; i++) nStrip[i]=32;
 885       nStrip[6]=16;
 886       for (Int_t i=0; i<3; i++) stripYsize[i]=stripWidth[0];
 887       for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
 888       stripYsize[6]=stripWidth[2];
 889       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
 890       stripXsize[6]=stripLength[2];
 891       trigSegX[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);
 892       trigSegX[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
 893
 894       // chamber 5 3 cathode 1
 895       for (Int_t i=0; i<5; i++) nStrip[i]=16;
 896       for (Int_t i=5; i<6; i++) nStrip[5]=8;
 897       nStrip[6]=16;
 898       stripYsize[0]=stripLength[2];
 899       for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
 900       for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
 901       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
 902       trigSegY[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);
 903       trigSegY[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
 904
 905       // chamber 4 cathode 0
 906       nStrip[0]=0;
 907       for (Int_t i=1; i<3; i++) nStrip[i]=64;
 908       for (Int_t i=3; i<6; i++) nStrip[i]=32;
 909       nStrip[6]=16;
 910       stripYsize[0]=0.;
 911       for (Int_t i=1; i<3; i++) stripYsize[i]=stripWidth[0];
 912       for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
 913       stripYsize[6]=stripWidth[2];
 914       stripXsize[0]=0;
 915       stripXsize[1]=stripLength[0];
 916       for (Int_t i=2; i<6; i++) stripXsize[i]=stripLength[1];
 917       stripXsize[6]=stripLength[2];
 918       trigSegX[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
 919
 920       // chamber 4 cathode 1
 921       nStrip[0]=0;
 922       nStrip[1]=8;
 923       for (Int_t i=2; i<5; i++) nStrip[i]=16;
 924       for (Int_t i=5; i<6; i++) nStrip[i]=8;
 925       nStrip[6]=16;
 926       stripYsize[0]=0.;
 927       for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
 928       stripXsize[0]=0.;
 929       for (Int_t i=1; i<5; i++) stripXsize[i]=stripWidth[1];
 930       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
 931       trigSegY[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
 932
 933       chamberSeg[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 934       chamberSeg[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
 935
 936       fMUON->SetNsec(chamber,2);
 937       Int_t icount=chamber-10;  // chamber counter (0 1 2 3)
 938       Int_t id0=(10+icount+1)*100;
 939
 940
 941       //  printf("in CreateTriggerSegmentation here 0 id0=%i \n",id0);
 942
 943       for (Int_t i = 0; i < 9; i++) {
 944
 945         // cathode 0
 946         chamberSeg[0]->Add(id0+i,     trigSegX[i]);
 947         chamberSeg[0]->Add(id0+50+i,  trigSegX[i]);
 948         fMUON->SetSegmentationModel(chamber, 1, chamberSeg[0]);
 949
 950         // cathode 1
 951         chamberSeg[1]->Add(id0+i,     trigSegY[i]);
 952         chamberSeg[1]->Add(id0+50+i,  trigSegY[i]);
 953         fMUON->SetSegmentationModel(chamber, 2, chamberSeg[1]);
 954
 955       }
 956
 957       fMUON->SetResponseModel(chamber, responseTrigger0);
 958       fMUON->Chamber(chamber).SetChargeCorrel(0); // same charge on cathodes
 959
 960
 961       //  printf("in CreateTriggerSegmentation here 1\n");
 962
 963       if (!id0) {
 964         AliWarning(Form("Segmentation for chamber %d is not yet defined",chamber));
 965         return ;
 966       }
 967     }
 968 }
 969 //__________________________________________________________________________
 970 void AliMUONFactoryV2::Build(AliMUON* where, const char* what)
 971 {
 972   //
 973   // Construct MUON from chambers, segmentation and responses
 974   //
 975
 976   fMUON = where;
 977   char tmp[20];
 978   strcpy(tmp, what);
 979
 980   if (strcmp(tmp, "default")==0) {
 981     // Set default parameters
 982     fMUON->SetIshunt(0);
 983     fMUON->SetMaxStepGas(0.1);
 984     fMUON->SetMaxStepAlu(0.1);
 985
 986     // Build stations
 987     BuildCommon();
 988     if (IsGeometryDefined(0))  BuildStation1();
 989     if (IsGeometryDefined(2))  BuildStation2();
 990     if (IsGeometryDefined(4))  BuildStation3();
 991     if (IsGeometryDefined(6))  BuildStation4();
 992     if (IsGeometryDefined(8))  BuildStation5();
 993     if (IsGeometryDefined(10)) BuildStation6();
 994   }
 995   else
 996     AliDebug(0,"Non default version of MUON selected. You have to construct yourself the MUON elements !!");
 997 }
 998
 999 //__________________________________________________________________________
1000 void AliMUONFactoryV2::BuildStation(AliMUON* where, Int_t stationNumber)
1001 {
1002   //
1003   // Construct MUON from chambers, segmentation and responses
1004   //
1005   // Version 0
1006   //
1007   // First define the number of planes that are segmented (1 or 2) by a call
1008   // to SetNsec.
1009   // Then chose for each chamber (chamber plane) the segmentation
1010   // and response model.
1011   // They should be equal for the two chambers of each station. In a future
1012   // version this will be enforced.
1013   //
1014
1015   fMUON = where;
1016   if (!fResponse0) BuildCommon();
1017
1018   switch (stationNumber) {
1019   case 1:  BuildStation1(); break;
1020   case 2:  BuildStation2(); break;
1021   case 3:  BuildStation3(); break;
1022   case 4:  BuildStation4(); break;
1023   case 5:  BuildStation5(); break;
1024   case 6:  BuildStation6(); break;
1025
1026   default: AliFatal("Wrong station number");
1027   }
1028 }