MUON/AliMUONSegFactory.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 // --------------------------
  23 // Class AliMUONSegFactory
  24 // --------------------------
  25 // New factory for building segmentations at all levels
  26 // Authors: Ivana Hrivnacova, IPN Orsay
  27
  28 #include "AliMUONSegFactory.h"
  29 #include "AliMUONConstants.h"
  30 #include "AliMUONGeometryStore.h"
  31 #include "AliMUONGeometryTransformer.h"
  32 #include "AliMUONGeometryModule.h"
  33 #include "AliMUONSegmentation.h"
  34 #include "AliMUONGeometrySegmentation.h"
  35 #include "AliMUONSt12QuadrantSegmentation.h"
  36 #include "AliMUONSt345SlatSegmentation.h"
  37 #include "AliMUONSt345SlatSegmentationV2.h"
  38 #include "AliMUONTriggerSegmentation.h"
  39 #include "AliMUONTriggerSegmentationV2.h"
  40 #include "AliMUONTriggerConstants.h"
  41
  42 #include "AliMpDEManager.h"
  43 #include "AliMpDEIterator.h"
  44
  45 #include "AliLog.h"
  46
  47 #include <Riostream.h>
  48 #include <TSystem.h>
  49 #include <TObjString.h>
  50 #include <TMap.h>
  51
  52 ClassImp(AliMUONSegFactory)
  53
  54 //______________________________________________________________________________
  55 AliMUONSegFactory::AliMUONSegFactory(const AliMUONGeometryTransformer* geometry)
  56     : TObject(),
  57       fMpSegFactory(),
  58       fDESegmentations(),
  59       fSegmentation(0),
  60       fkTransformer(geometry)
  61 {
  62 /// Standard constructor
  63
  64 }
  65
  66 //______________________________________________________________________________
  67 AliMUONSegFactory::AliMUONSegFactory(const TString& volPathsFileName,
  68                                      const TString& transformsFileName)
  69     : TObject(),
  70       fMpSegFactory(),
  71       fDESegmentations(),
  72       fSegmentation(0),
  73       fkTransformer(0)
  74 {
  75 /// Standard constructor
  76
  77   // Transformer
  78   AliMUONGeometryTransformer* transformer = new AliMUONGeometryTransformer(true);
  79   transformer->ReadGeometryData(volPathsFileName, transformsFileName);
  80   fkTransformer = transformer;
  81 }
  82
  83 //______________________________________________________________________________
  84   AliMUONSegFactory::AliMUONSegFactory()
  85     : TObject(),
  86       fMpSegFactory(),
  87       fDESegmentations(),
  88       fSegmentation(0),
  89       fkTransformer(0)
  90 {
  91 /// Default constructor
  92 }
  93
  94 //______________________________________________________________________________
  95 AliMUONSegFactory::AliMUONSegFactory(const AliMUONSegFactory& rhs)
  96  : TObject(rhs)
  97 {
  98 /// Protected copy constructor
  99
 100   AliFatal("Not implemented.");
 101 }
 102
 103 //______________________________________________________________________________
 104
 105 AliMUONSegFactory::~AliMUONSegFactory()
 106 {
 107 /// Destructor
 108
 109   //delete fSegmentation;
 110         // The segmentation is supposed to be deleted in the client code
 111 }
 112
 113 //______________________________________________________________________________
 114 AliMUONSegFactory&  AliMUONSegFactory::operator=(const AliMUONSegFactory& rhs)
 115 {
 116   // Protected assignement operator
 117
 118   if (this == &rhs) return *this;
 119
 120   AliFatal("Not implemented.");
 121
 122   return *this;
 123 }
 124
 125 //
 126 // Private methods
 127 //
 128
 129 //______________________________________________________________________________
 130 Bool_t AliMUONSegFactory::IsGeometryDefined(Int_t ichamber)
 131 {
 132 // Return true, if det elements for the chamber with the given ichamber Id
 133 // are defined in geometry (the geometry builder for this chamber was activated)
 134
 135   if ( ! fkTransformer ||
 136        ! fkTransformer->GetModuleTransformer(ichamber, false) )
 137
 138     return kFALSE;
 139
 140   return kTRUE;
 141 }
 142
 143 //__________________________________________________________________________
 144 AliMUONSegmentation* AliMUONSegFactory::Segmentation()
 145 {
 146 /// Return the segmentation container, create it if it does not yet exist
 147
 148   if ( ! fSegmentation )
 149     fSegmentation = new AliMUONSegmentation(AliMUONConstants::NCh());
 150
 151   return fSegmentation;
 152 }
 153
 154 //
 155 // public methods
 156 //
 157
 158 //______________________________________________________________________________
 159 AliMpVSegmentation*
 160 AliMUONSegFactory::CreateMpSegmentation(Int_t detElemId, Int_t cath)
 161 {
 162 /// Create mapping segmentation for given detElemId and cath
 163 /// using mapping manager
 164
 165   AliMpVSegmentation* mpSegmentation
 166     = fMpSegFactory.CreateMpSegmentation(detElemId, cath);
 167
 168   Segmentation()->AddMpSegmentation(mpSegmentation);
 169
 170   return mpSegmentation;
 171 }
 172
 173 //______________________________________________________________________________
 174 AliMUONVGeometryDESegmentation*
 175 AliMUONSegFactory::CreateDESegmentation(Int_t detElemId, Int_t cath)
 176 {
 177 // Create DE segmentation, operating in local DE reference frame
 178
 179   // Check detElemId & cath
 180   if ( ! AliMpDEManager::IsValid(detElemId, cath, true) ) return 0;
 181
 182   // Check if transformer is defined
 183   if ( ! fkTransformer) {
 184     AliErrorStream() << "Geometry transformer not defined" << endl;
 185     return 0;
 186   }
 187
 188   // Only return it, if DE segmentation for this detElemId and cath
 189   // was already defined
 190   //
 191   const AliMUONVGeometryDESegmentation* kdeSegmentation
 192     = Segmentation()->GetDESegmentation(detElemId, cath, false);
 193   if ( kdeSegmentation )
 194     return const_cast<AliMUONVGeometryDESegmentation*>(kdeSegmentation);
 195
 196   // Get module, create it if it does not exist
 197   //
 198   Int_t moduleId = AliMUONGeometryStore::GetModuleId(detElemId);
 199   AliMUONGeometrySegmentation* moduleSegmentation
 200     = Segmentation()->GetModuleSegmentation(moduleId, cath, false);
 201   if (! moduleSegmentation) {
 202     moduleSegmentation
 203       = new AliMUONGeometrySegmentation(
 204                fkTransformer->GetModuleTransformer(moduleId));
 205     Segmentation()->AddModuleSegmentation(moduleId, cath, moduleSegmentation);
 206   }
 207
 208   // Get DE segmentation for this DE type, create it if it does not exist
 209   //
 210   AliMUONVGeometryDESegmentation* deSegmentation = 0;
 211   TString deName = AliMpDEManager::GetDEName(detElemId, cath);
 212   TObject* objSegmentation = fDESegmentations.Get(deName);
 213   if ( objSegmentation )
 214     deSegmentation = (AliMUONVGeometryDESegmentation*)objSegmentation;
 215
 216   if ( !deSegmentation ) {
 217
 218     // Get/Create mapping segmentation via mapping manager
 219     AliMpVSegmentation* mpSegmentation
 220       = CreateMpSegmentation(detElemId, cath);
 221
 222     AliMpStationType stationType = AliMpDEManager::GetStationType(detElemId);
 223     AliMpPlaneType planeType = AliMpDEManager::GetPlaneType(detElemId, cath);
 224
 225     switch (stationType) {
 226
 227       case kStation1:
 228       case kStation2:
 229         deSegmentation = new AliMUONSt12QuadrantSegmentation(
 230                                  mpSegmentation, stationType, planeType);
 231         //cout << "   new AliMUONSt12QuadrantSegmentation "
 232         //     << StationTypeName(stationType) << "  "
 233         //     << PlaneTypeName(planeType) << "  "
 234         //     << deName << endl;
 235
 236         break;
 237
 238       case kStation345:
 239         deSegmentation = new AliMUONSt345SlatSegmentationV2(
 240                                  mpSegmentation, detElemId, planeType);
 241         //cout << "   new AliMUONSt345SlatSegmentationV2 "
 242         //     << StationTypeName(stationType) << "  "
 243         //     << PlaneTypeName(planeType) << "  "
 244         //     << deName << endl;
 245         break;
 246
 247       case kStationTrigger:
 248         deSegmentation = new AliMUONTriggerSegmentationV2(
 249                                  mpSegmentation, detElemId, planeType);
 250         //cout << "   new AliMUONTriggerSegmentationV2 "
 251         //     << StationTypeName(stationType) << "  "
 252         //     << PlaneTypeName(planeType) << "  "
 253         //     << deName << endl;
 254         break;
 255     }
 256
 257     // Map new DE segmentation
 258     fDESegmentations.Add(deName, deSegmentation);
 259     Segmentation()->AddDESegmentation(deSegmentation);
 260   }
 261
 262   // Add  DE segmentation to module
 263   //
 264   moduleSegmentation->Add(detElemId, deName, deSegmentation);
 265
 266   return deSegmentation;
 267 }
 268
 269
 270 //______________________________________________________________________________
 271 AliMUONGeometrySegmentation*
 272 AliMUONSegFactory::CreateModuleSegmentation(Int_t moduleId, Int_t cath)
 273 {
 274 // Create module segmentation, operating in global reference frame
 275 // Detection elements are defined via DE names map.
 276
 277   // Check cathod & module Id
 278   if ( ! AliMpDEManager::IsValidCathod(cath, true) ||
 279        ! AliMpDEManager::IsValidModuleId(moduleId, true) ) return 0;
 280
 281   AliMpDEIterator it;
 282   for ( it.First(moduleId); ! it.IsDone(); it.Next() )
 283     CreateDESegmentation(it.CurrentDE(), cath);
 284
 285   return Segmentation()->GetModuleSegmentation(moduleId, cath);
 286 }
 287
 288 //______________________________________________________________________________
 289 AliMUONSegmentation*
 290 AliMUONSegFactory::CreateSegmentation(const TString& option)
 291 {
 292 /// Create segmentations on all levels and return their container.
 293
 294   // Check options
 295   if ( option != "default"   &&
 296        option != "FactoryV2" &&
 297        option != "FactoryV3" &&
 298        option != "FactoryV4" &&
 299        option != "new") {
 300
 301     AliErrorStream() << "Option " << option << " not defined." << endl;
 302     return 0;
 303   }
 304
 305   if ( option == "FactoryV2" ) {
 306     // Default segmentation version
 307     for (Int_t moduleId = 0; moduleId<4; moduleId++)
 308       for (Int_t cath = 0; cath < 2; cath++) {
 309         if ( IsGeometryDefined(moduleId) )
 310           CreateModuleSegmentation( moduleId, cath);
 311     }
 312     if ( IsGeometryDefined(4) )  BuildStation3();
 313     if ( IsGeometryDefined(6) )  BuildStation4();
 314     if ( IsGeometryDefined(8) )  BuildStation5();
 315     if ( IsGeometryDefined(10))  BuildStation6();
 316   }
 317
 318   if ( option == "FactoryV3" ) {
 319     // New slat segmentation
 320     for (Int_t moduleId = 0; moduleId<10; moduleId++)
 321       for (Int_t cath = 0; cath < 2; cath++) {
 322         if ( IsGeometryDefined(moduleId) )
 323            CreateModuleSegmentation( moduleId, cath);
 324     }
 325     if ( IsGeometryDefined(10) ) BuildStation6();
 326   }
 327
 328   if (option == "default" || option == "new" || option == "FactoryV4" ) {
 329
 330     for (Int_t moduleId = 0; moduleId<AliMUONConstants::NCh(); moduleId++)
 331       for (Int_t cath = 0; cath < 2; cath++) {
 332         if ( IsGeometryDefined(moduleId) )
 333           CreateModuleSegmentation( moduleId, cath);
 334     }
 335   }
 336   return Segmentation();
 337 }
 338
 339 //
 340 //  Functions for building old segmentations (not based on mapping)
 341 //
 342
 343 //__________________________________________________________________________
 344 void AliMUONSegFactory::BuildStation3()
 345 {
 346   //--------------------------------------------------------
 347   // Configuration for Chamber TC5/6  (Station 3) ----------
 348   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 349
 350   AliMUONGeometrySegmentation* segmentation[2];
 351
 352   //Slats Segmentations
 353   AliMUONSt345SlatSegmentation *slatsegB[4]; // Types of segmentation for St3
 354   AliMUONSt345SlatSegmentation *slatsegNB[4];
 355   // Bending
 356
 357   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones
 358   for(Int_t i=0; i<4; i++) {
 359     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
 360     Segmentation()->AddDESegmentation(slatsegB[i]);
 361     slatsegB[i]->SetPadSize(10.,0.5);
 362     slatsegB[i]->SetPadDivision(ndiv);
 363     slatsegB[i]->SetId(1); // Id elt ????
 364     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
 365     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
 366     Segmentation()->AddDESegmentation(slatsegNB[i]);
 367     slatsegNB[i]->SetPadSize(1./1.4,10.); // Nbending
 368     slatsegNB[i]->SetPadDivision(ndiv);
 369     slatsegNB[i]->SetId(1);
 370     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
 371   }
 372
 373   // Type 112200 for 500, 501, 508, 509, 510, 517
 374   // in Ch5 (similar for Ch6) for the futur official numbering
 375   // Type 112200 for 503, 504, 505, 555, 554, 553
 376   // in Ch5 (similar for Ch6) actual numbering in the code to be changed in jan05
 377   Int_t n0[4] = { 0, 2, 2, 0 };
 378   slatsegB[0]->SetPcbBoards(n0);
 379   slatsegB[0]->Init(0);
 380   slatsegNB[0]->SetPcbBoards(n0);
 381   slatsegNB[0]->Init(0);
 382
 383   // Type 122200 for 502, 507, 511, 516 (similar in Ch6)
 384   // for future official numbering of ALICE
 385   // Type 122200 for 502, 506, 556, 552 (similiarin Ch6)
 386   // for actual numbering in muon code to be changed in jan05
 387   Int_t n1[4] = { 0, 1, 3, 0 };
 388   slatsegB[1]->SetPcbBoards(n1);
 389   slatsegB[1]->Init(0);
 390   slatsegNB[1]->SetPcbBoards(n1);
 391   slatsegNB[1]->Init(0);
 392
 393   // Type 222000 for 503, 506, 512, 515 (similar in Ch6)
 394   // for future official numbering of ALICE
 395   // Type 222000 for 501, 507, 557, 551 (similiarin Ch6)
 396   // for actual numbering in muon code to be changed in jan05
 397   Int_t n2[4] = { 0, 0, 3, 0 };
 398   slatsegB[2]->SetPcbBoards(n2);
 399   slatsegB[2]->Init(0);
 400   slatsegNB[2]->SetPcbBoards(n2);
 401   slatsegNB[2]->Init(0);
 402
 403   // Type 220000 for 504, 505, 513, 514 (similar in Ch6)
 404   // for future official numbering of ALICE
 405   // Type 220000 for 500, 508, 558, 550 (similiarin Ch6)
 406   // for actual numbering in muon code to be changed in jan05
 407   Int_t n3[4] = { 0, 0, 2, 0 };
 408   slatsegB[3]->SetPcbBoards(n3);
 409   slatsegB[3]->Init(0);
 410   slatsegNB[3]->SetPcbBoards(n3);
 411   slatsegNB[3]->Init(0);
 412
 413   for (Int_t chamber = 4; chamber < 6; chamber++) {
 414
 415     const AliMUONGeometryModuleTransformer* kModuleTransformer
 416       = fkTransformer->GetModuleTransformer(chamber);
 417
 418     segmentation[0] = new AliMUONGeometrySegmentation(kModuleTransformer);
 419     segmentation[1] = new AliMUONGeometrySegmentation(kModuleTransformer);
 420
 421     // id detection elt for chamber 1
 422     Int_t id0 = (chamber+1)*100;
 423
 424     // cathode 0
 425     // type 220000
 426     segmentation[0]->Add(id0+14, "Undefined", slatsegB[3]);
 427     segmentation[0]->Add(id0+ 4, "Undefined", slatsegB[3]);
 428     segmentation[0]->Add(id0+13, "Undefined", slatsegB[3]);
 429     segmentation[0]->Add(id0+ 5, "Undefined", slatsegB[3]);
 430     // type 222000
 431     segmentation[0]->Add(id0+15, "Undefined", slatsegB[2]);
 432     segmentation[0]->Add(id0+ 3, "Undefined", slatsegB[2]);
 433     segmentation[0]->Add(id0+12, "Undefined", slatsegB[2]);
 434     segmentation[0]->Add(id0+ 6, "Undefined", slatsegB[2]);
 435     // type 122200
 436     segmentation[0]->Add(id0+16, "Undefined", slatsegB[1]);
 437     segmentation[0]->Add(id0+ 2, "Undefined", slatsegB[1]);
 438     segmentation[0]->Add(id0+11, "Undefined", slatsegB[1]);
 439     segmentation[0]->Add(id0+ 7, "Undefined", slatsegB[1]);
 440     // type 112200
 441     segmentation[0]->Add(id0+17, "Undefined", slatsegB[0]);
 442     segmentation[0]->Add(id0,    "Undefined", slatsegB[0]);
 443     segmentation[0]->Add(id0+ 1, "Undefined", slatsegB[0]);
 444     segmentation[0]->Add(id0+10, "Undefined", slatsegB[0]);
 445     segmentation[0]->Add(id0+ 9, "Undefined", slatsegB[0]);
 446     segmentation[0]->Add(id0+ 8, "Undefined", slatsegB[0]);
 447     Segmentation()->AddModuleSegmentation(chamber, 0, segmentation[0]);
 448
 449     // cathode 1
 450     // type 220000
 451     segmentation[1]->Add(id0+14, "Undefined", slatsegNB[3]);
 452     segmentation[1]->Add(id0+ 4, "Undefined", slatsegNB[3]);
 453     segmentation[1]->Add(id0+13, "Undefined", slatsegNB[3]);
 454     segmentation[1]->Add(id0+ 5, "Undefined", slatsegNB[3]);
 455     // type 222000
 456     segmentation[1]->Add(id0+15, "Undefined", slatsegNB[2]);
 457     segmentation[1]->Add(id0+ 3, "Undefined", slatsegNB[2]);
 458     segmentation[1]->Add(id0+12, "Undefined", slatsegNB[2]);
 459     segmentation[1]->Add(id0+ 6, "Undefined", slatsegNB[2]);
 460     // type 122200
 461     segmentation[1]->Add(id0+16, "Undefined", slatsegNB[1]);
 462     segmentation[1]->Add(id0+ 2, "Undefined", slatsegNB[1]);
 463     segmentation[1]->Add(id0+11, "Undefined", slatsegNB[1]);
 464     segmentation[1]->Add(id0+ 7, "Undefined", slatsegNB[1]);
 465     // type 112200
 466     segmentation[1]->Add(id0+17, "Undefined", slatsegNB[0]);
 467     segmentation[1]->Add(id0,    "Undefined", slatsegNB[0]);
 468     segmentation[1]->Add(id0+ 1, "Undefined", slatsegNB[0]);
 469     segmentation[1]->Add(id0+10, "Undefined", slatsegNB[0]);
 470     segmentation[1]->Add(id0+ 9, "Undefined", slatsegNB[0]);
 471     segmentation[1]->Add(id0+ 8, "Undefined", slatsegNB[0]);
 472     Segmentation()->AddModuleSegmentation(chamber, 1, segmentation[1]);
 473   }
 474 }
 475
 476 //__________________________________________________________________________
 477 void AliMUONSegFactory::BuildStation4()
 478 {
 479   //--------------------------------------------------------
 480   // Configuration for Chamber TC7/8  (Station 4) ----------
 481   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 482
 483
 484   AliMUONGeometrySegmentation* segmentation[2];
 485
 486   //Slats Segmentations
 487   AliMUONSt345SlatSegmentation *slatsegB[7]; // Types of segmentation for St4
 488   AliMUONSt345SlatSegmentation *slatsegNB[7];
 489   // Bending
 490
 491   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones
 492   for(Int_t i = 0; i < 7; i++) {
 493     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
 494     Segmentation()->AddDESegmentation(slatsegB[i]);
 495     slatsegB[i]->SetPadSize(10.,0.5);
 496     slatsegB[i]->SetPadDivision(ndiv);
 497     slatsegB[i]->SetId(1);
 498     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
 499     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
 500     Segmentation()->AddDESegmentation(slatsegNB[i]);
 501     slatsegNB[i]->SetPadSize(1./1.4,10.);
 502     slatsegNB[i]->SetPadDivision(ndiv);
 503     slatsegNB[i]->SetId(1);
 504     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
 505   }
 506
 507   Int_t n4[4] = { 0, 1, 2, 2 };
 508   slatsegB[0]->SetPcbBoards(n4);
 509   slatsegB[0]->Init(0); // 0 detection element id
 510   slatsegNB[0]->SetPcbBoards(n4);
 511   slatsegNB[0]->Init(0); // 0 detection element id
 512
 513   // Type 112233 for 701, 712, 714, 725 in Ch7 (similar for Ch8)
 514   // for the futur official numbering
 515   // Type 112233 for 705, 707, 755, 757 in Ch7 (similar for Ch8)
 516   // actual numbering in the code to be changed in jan05
 517   // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9
 518   // (similar for Ch10) for the futur official numbering
 519   // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9
 520   // (similar for Ch10) actual numbering in the code to be changed in jan05
 521   Int_t n5[4] = { 0, 2, 2, 2 };
 522   slatsegB[1]->SetPcbBoards(n5);
 523   slatsegB[1]->Init(0); // 0 detection element id
 524   slatsegNB[1]->SetPcbBoards(n5);
 525   slatsegNB[1]->Init(0); // 0 detection element id
 526
 527   // Type 112230 for 702, 711, 715, 724 in Ch7 (similar for Ch8)
 528   // for the futur official numbering
 529   // Type 112230 for 704, 708, 754, 758 in Ch7 (similar for Ch8)
 530   // actual numbering in the code to be changed in jan05
 531   Int_t n6[4] = { 0, 2, 2, 1 };
 532   slatsegB[2]->SetPcbBoards(n6);
 533   slatsegB[2]->Init(0); // 0 detection element id
 534   slatsegNB[2]->SetPcbBoards(n6);
 535   slatsegNB[2]->Init(0); // 0 detection element id
 536
 537   // Type 222330 for 703, 710, 716, 723 in Ch7 (similar for Ch8)
 538   // for the futur official numbering
 539   // Type 222330 for 703, 709, 753, 759 in Ch7 (similar for Ch8)
 540   // actual numbering in the code to be changed in jan05
 541   Int_t n7[4] = { 0, 0, 3, 2 };
 542   slatsegB[3]->SetPcbBoards(n7);
 543   slatsegB[3]->Init(0); // 0 detection element id
 544   slatsegNB[3]->SetPcbBoards(n7);
 545   slatsegNB[3]->Init(0); // 0 detection element id
 546
 547   // Type 223300 for 704, 709, 717, 722 in Ch7 (similar for Ch8)
 548   // for the futur official numbering
 549   // Type 223300 for 702, 710, 752, 760 in Ch7 (similar for Ch8)
 550   // actual numbering in the code to be changed in jan05
 551   Int_t n8[4] = { 0, 0, 2, 2 };
 552   slatsegB[4]->SetPcbBoards(n8);
 553   slatsegB[4]->Init(0); // 0 detection element id
 554   slatsegNB[4]->SetPcbBoards(n8);
 555   slatsegNB[4]->Init(0); // 0 detection element id
 556
 557   // Type 333000 for 705, 708, 718, 721 in Ch7 (similar for Ch8)
 558   // for the futur official numbering
 559   // Type 333000 for 701, 711, 751, 761 in Ch7 (similar for Ch8)
 560   // actual numbering in the code to be changed in jan05
 561   // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10)
 562   // for the futur official numbering
 563   // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10)
 564   // actual numbering in the code to be changed in jan05
 565   Int_t n9[4] = { 0, 0, 0, 3 };
 566   slatsegB[5]->SetPcbBoards(n9);
 567   slatsegB[5]->Init(0); // 0 detection element id
 568   slatsegNB[5]->SetPcbBoards(n9);
 569   slatsegNB[5]->Init(0); // 0 detection element id
 570
 571   // Type 330000 for 706, 707, 719, 720 in Ch7 (similar for Ch8)
 572   // for the futur official numbering
 573   // Type 330000 for 700, 712, 750, 762 in Ch7 (similar for Ch8)
 574   // actual numbering in the code to be changed in jan05
 575   Int_t n10[4] = { 0, 0, 0, 2 };
 576   slatsegB[6]->SetPcbBoards(n10);
 577   slatsegB[6]->Init(0); // 0 detection element id
 578   slatsegNB[6]->SetPcbBoards(n10);
 579   slatsegNB[6]->Init(0); // 0 detection element id
 580
 581
 582   for (Int_t chamber = 6; chamber < 8; chamber++) {
 583
 584     const AliMUONGeometryModuleTransformer* kModuleTransformer
 585       = fkTransformer->GetModuleTransformer(chamber);
 586
 587     segmentation[0] = new AliMUONGeometrySegmentation(kModuleTransformer);
 588     segmentation[1] = new AliMUONGeometrySegmentation(kModuleTransformer);
 589
 590     // id detection elt for chamber 1
 591     Int_t id0 = (chamber+1)*100;
 592
 593     //--------------------------------------------------------
 594     // Configuration for Chamber TC6/7  (Station 4) ----------
 595
 596     // cathode 0
 597     // type 122330
 598     segmentation[0]->Add(id0+13, "Undefined", slatsegB[0]);
 599     segmentation[0]->Add(id0   , "Undefined", slatsegB[0]);
 600
 601     // type 112233
 602     segmentation[0]->Add(id0+14, "Undefined", slatsegB[1]);
 603     segmentation[0]->Add(id0+12, "Undefined", slatsegB[1]);
 604     segmentation[0]->Add(id0+25, "Undefined", slatsegB[1]);
 605     segmentation[0]->Add(id0+ 1, "Undefined", slatsegB[1]);
 606
 607     // type 112230
 608     segmentation[0]->Add(id0+15, "Undefined", slatsegB[2]);
 609     segmentation[0]->Add(id0+11, "Undefined", slatsegB[2]);
 610     segmentation[0]->Add(id0+24, "Undefined", slatsegB[2]);
 611     segmentation[0]->Add(id0+ 2, "Undefined", slatsegB[2]);
 612
 613     // type 222330
 614     segmentation[0]->Add(id0+16, "Undefined", slatsegB[3]);
 615     segmentation[0]->Add(id0+10, "Undefined", slatsegB[3]);
 616     segmentation[0]->Add(id0+23, "Undefined", slatsegB[3]);
 617     segmentation[0]->Add(id0+ 3, "Undefined", slatsegB[3]);
 618
 619     // type 223300
 620     segmentation[0]->Add(id0+17, "Undefined", slatsegB[4]);
 621     segmentation[0]->Add(id0+ 9, "Undefined", slatsegB[4]);
 622     segmentation[0]->Add(id0+22, "Undefined", slatsegB[4]);
 623     segmentation[0]->Add(id0+ 4, "Undefined", slatsegB[4]);
 624
 625     // type 333000
 626     segmentation[0]->Add(id0+18, "Undefined", slatsegB[5]);
 627     segmentation[0]->Add(id0+ 8, "Undefined", slatsegB[5]);
 628     segmentation[0]->Add(id0+21, "Undefined", slatsegB[5]);
 629     segmentation[0]->Add(id0+ 5, "Undefined", slatsegB[5]);
 630
 631     // type 330000
 632     segmentation[0]->Add(id0+19, "Undefined", slatsegB[6]);
 633     segmentation[0]->Add(id0+ 7, "Undefined", slatsegB[6]);
 634     segmentation[0]->Add(id0+20, "Undefined", slatsegB[6]);
 635     segmentation[0]->Add(id0+ 6, "Undefined", slatsegB[6]);
 636     Segmentation()->AddModuleSegmentation(chamber, 0, segmentation[0]);
 637
 638     // cathode 1
 639     // type 122330
 640     segmentation[1]->Add(id0+13, "Undefined", slatsegNB[0]);
 641     segmentation[1]->Add(id0   , "Undefined", slatsegNB[0]);
 642
 643     // type 112233
 644     segmentation[1]->Add(id0+14, "Undefined", slatsegNB[1]);
 645     segmentation[1]->Add(id0+12, "Undefined", slatsegNB[1]);
 646     segmentation[1]->Add(id0+25, "Undefined", slatsegNB[1]);
 647     segmentation[1]->Add(id0+ 1, "Undefined", slatsegNB[1]);
 648
 649     // type 112230
 650     segmentation[1]->Add(id0+15, "Undefined", slatsegNB[2]);
 651     segmentation[1]->Add(id0+11, "Undefined", slatsegNB[2]);
 652     segmentation[1]->Add(id0+24, "Undefined", slatsegNB[2]);
 653     segmentation[1]->Add(id0+ 2, "Undefined", slatsegNB[2]);
 654
 655     // type 222330
 656     segmentation[1]->Add(id0+16, "Undefined", slatsegNB[3]);
 657     segmentation[1]->Add(id0+10, "Undefined", slatsegNB[3]);
 658     segmentation[1]->Add(id0+23, "Undefined", slatsegNB[3]);
 659     segmentation[1]->Add(id0+ 3, "Undefined", slatsegNB[3]);
 660
 661     // type 223300
 662     segmentation[1]->Add(id0+17, "Undefined", slatsegNB[4]);
 663     segmentation[1]->Add(id0+ 9, "Undefined", slatsegNB[4]);
 664     segmentation[1]->Add(id0+22, "Undefined", slatsegNB[4]);
 665     segmentation[1]->Add(id0+ 4, "Undefined", slatsegNB[4]);
 666
 667     // type 333000
 668     segmentation[1]->Add(id0+18, "Undefined", slatsegNB[5]);
 669     segmentation[1]->Add(id0+ 8, "Undefined", slatsegNB[5]);
 670     segmentation[1]->Add(id0+21, "Undefined", slatsegNB[5]);
 671     segmentation[1]->Add(id0+ 5, "Undefined", slatsegNB[5]);
 672
 673     // type 330000
 674     segmentation[1]->Add(id0+19, "Undefined", slatsegNB[6]);
 675     segmentation[1]->Add(id0+ 7, "Undefined", slatsegNB[6]);
 676     segmentation[1]->Add(id0+20, "Undefined", slatsegNB[6]);
 677     segmentation[1]->Add(id0+ 6, "Undefined", slatsegNB[6]);
 678     Segmentation()->AddModuleSegmentation(chamber, 1, segmentation[1]);
 679   }
 680 }
 681
 682 //__________________________________________________________________________
 683 void AliMUONSegFactory::BuildStation5()
 684 {
 685   //--------------------------------------------------------
 686   // Configuration for Chamber TC9/10  (Station 5) ---------
 687   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 688
 689   AliMUONGeometrySegmentation* segmentation[2];
 690
 691   //Slats Segmentations
 692   AliMUONSt345SlatSegmentation *slatsegB[6]; // Types of segmentation for St5
 693   AliMUONSt345SlatSegmentation *slatsegNB[6];
 694   // Bending
 695
 696   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones
 697   for(Int_t i = 0; i < 6; i++) {
 698     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
 699     Segmentation()->AddDESegmentation(slatsegB[i]);
 700     slatsegB[i]->SetPadSize(10.,0.5);
 701     slatsegB[i]->SetPadDivision(ndiv);
 702     slatsegB[i]->SetId(1);
 703     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
 704     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
 705     Segmentation()->AddDESegmentation(slatsegNB[i]);
 706     slatsegNB[i]->SetPadSize(1./1.4,10.);
 707     slatsegNB[i]->SetPadDivision(ndiv);
 708     slatsegNB[i]->SetId(1);
 709     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
 710   }
 711
 712   // Type 122330 for 900, 913 in Ch9 (similar for Ch10)
 713   // for the futur official numbering
 714   // Type 122330 for 906, 956 in Ch9 (similar for Ch10)
 715   // actual numbering in the code to be changed in jan05
 716   Int_t n4[4] = { 0, 1, 2, 2 };
 717   slatsegB[0]->SetPcbBoards(n4);
 718   slatsegB[0]->Init(0); // 0 detection element id
 719   slatsegNB[0]->SetPcbBoards(n4);
 720   slatsegNB[0]->Init(0); // 0 detection element id
 721
 722   // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9
 723   // (similar for Ch10) for the futur official numbering
 724   // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9
 725   // (similar for Ch10) actual numbering in the code to be changed in jan05
 726   Int_t n5[4] = { 0, 2, 2, 2 };
 727   slatsegB[1]->SetPcbBoards(n5);
 728   slatsegB[1]->Init(0); // 0 detection element id
 729   slatsegNB[1]->SetPcbBoards(n5);
 730   slatsegNB[1]->Init(0); // 0 detection element id
 731
 732   // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10)
 733   // for the futur official numbering
 734   // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10)
 735   // actual numbering in the code to be changed in jan05
 736   Int_t n9[4] = { 0, 0, 0, 3 };
 737   slatsegB[2]->SetPcbBoards(n9);
 738   slatsegB[2]->Init(0); // 0 detection element id
 739   slatsegNB[2]->SetPcbBoards(n9);
 740   slatsegNB[2]->Init(0); // 0 detection element id
 741
 742   // Type 222333 for 903, 910, 916, 923 in Ch9 (similar for Ch10)
 743   // for the futur official numbering
 744   // Type 222333 for 903, 909, 953, 959 in Ch9 (similar for Ch10)
 745   // actual numbering in the code to be changed in jan05
 746   Int_t n11[4] = { 0, 0, 3, 3 };
 747   slatsegB[3]->SetPcbBoards(n11);
 748   slatsegB[3]->Init(0); // 0 detection element id
 749   slatsegNB[3]->SetPcbBoards(n11);
 750   slatsegNB[3]->Init(0); // 0 detection element id
 751
 752   // Type 223330 for 904, 909, 917, 922 in Ch9 (similar for Ch10)
 753   // for the futur official numbering
 754   // Type 223330 for 902, 910, 952, 960 in Ch9 (similar for Ch10)
 755   // actual numbering in the code to be changed in jan05
 756   Int_t n12[4] = { 0, 0, 2, 3 };
 757   slatsegB[4]->SetPcbBoards(n12);
 758   slatsegB[4]->Init(0); // 0 detection element id
 759   slatsegNB[4]->SetPcbBoards(n12);
 760   slatsegNB[4]->Init(0); // 0 detection element id
 761
 762   // Type 333300 for 905, 908, 918, 921 in Ch9 (similar for Ch10)
 763   // for the futur official numbering
 764   // Type 333300 for 901, 911, 951, 961 in Ch9 (similar for Ch10)
 765   // actual numbering in the code to be changed in jan05
 766   Int_t n13[4] = { 0, 0, 0, 4 };
 767   slatsegB[5]->SetPcbBoards(n13);
 768   slatsegB[5]->Init(0); // 0 detection element id
 769   slatsegNB[5]->SetPcbBoards(n13);
 770   slatsegNB[5]->Init(0); // 0 detection element id
 771
 772   for (Int_t chamber = 8; chamber < 10; chamber++) {
 773
 774     const AliMUONGeometryModuleTransformer* kModuleTransformer
 775       = fkTransformer->GetModuleTransformer(chamber);
 776
 777     segmentation[0] = new AliMUONGeometrySegmentation(kModuleTransformer);
 778     segmentation[1] = new AliMUONGeometrySegmentation(kModuleTransformer);
 779
 780     // id detection elt for chamber 1
 781     Int_t id0 = (chamber+1)*100;
 782
 783     //--------------------------------------------------------
 784     // Configuration for Chamber TC8/9  (Station 5) ----------
 785
 786     // cathode 0
 787     // type 122330
 788     segmentation[0]->Add(id0+13, "Undefined", slatsegB[0]);
 789     segmentation[0]->Add(id0   , "Undefined", slatsegB[0]);
 790
 791     // type 112233
 792     segmentation[0]->Add(id0+15, "Undefined", slatsegB[1]);
 793     segmentation[0]->Add(id0+14, "Undefined", slatsegB[1]);
 794     segmentation[0]->Add(id0+12, "Undefined", slatsegB[1]);
 795     segmentation[0]->Add(id0+11, "Undefined", slatsegB[1]);
 796     segmentation[0]->Add(id0+24, "Undefined", slatsegB[1]);
 797     segmentation[0]->Add(id0+25, "Undefined", slatsegB[1]);
 798     segmentation[0]->Add(id0+ 1, "Undefined", slatsegB[1]);
 799     segmentation[0]->Add(id0+ 2, "Undefined", slatsegB[1]);
 800
 801     // type 333000
 802     segmentation[0]->Add(id0+19, "Undefined", slatsegB[2]);
 803     segmentation[0]->Add(id0+ 7, "Undefined", slatsegB[2]);
 804     segmentation[0]->Add(id0+20, "Undefined", slatsegB[2]);
 805     segmentation[0]->Add(id0+ 6, "Undefined", slatsegB[2]);
 806
 807     // type 222333
 808     segmentation[0]->Add(id0+16, "Undefined", slatsegB[3]);
 809     segmentation[0]->Add(id0+10, "Undefined", slatsegB[3]);
 810     segmentation[0]->Add(id0+23, "Undefined", slatsegB[3]);
 811     segmentation[0]->Add(id0+ 3, "Undefined", slatsegB[3]);
 812
 813     // type 223330
 814     segmentation[0]->Add(id0+17, "Undefined", slatsegB[4]);
 815     segmentation[0]->Add(id0+ 9, "Undefined", slatsegB[4]);
 816     segmentation[0]->Add(id0+22, "Undefined", slatsegB[4]);
 817     segmentation[0]->Add(id0+ 4, "Undefined", slatsegB[4]);
 818
 819     // type 333300
 820     segmentation[0]->Add(id0+18, "Undefined", slatsegB[5]);
 821     segmentation[0]->Add(id0+ 8, "Undefined", slatsegB[5]);
 822     segmentation[0]->Add(id0+21, "Undefined", slatsegB[5]);
 823     segmentation[0]->Add(id0+ 5, "Undefined", slatsegB[5]);
 824     Segmentation()->AddModuleSegmentation(chamber, 0, segmentation[0]);
 825
 826     // cathode 1
 827     // type 122330
 828     segmentation[1]->Add(id0+13, "Undefined", slatsegNB[0]);
 829     segmentation[1]->Add(id0   , "Undefined", slatsegNB[0]);
 830
 831     // type 112233
 832     segmentation[1]->Add(id0+15, "Undefined", slatsegNB[1]);
 833     segmentation[1]->Add(id0+14, "Undefined", slatsegNB[1]);
 834     segmentation[1]->Add(id0+12, "Undefined", slatsegNB[1]);
 835     segmentation[1]->Add(id0+11, "Undefined", slatsegNB[1]);
 836     segmentation[1]->Add(id0+24, "Undefined", slatsegNB[1]);
 837     segmentation[1]->Add(id0+25, "Undefined", slatsegNB[1]);
 838     segmentation[1]->Add(id0+ 1, "Undefined", slatsegNB[1]);
 839     segmentation[1]->Add(id0+ 2, "Undefined", slatsegNB[1]);
 840
 841     // type 333000
 842     segmentation[1]->Add(id0+19 , "Undefined", slatsegNB[2]);
 843     segmentation[1]->Add(id0+ 7, "Undefined", slatsegNB[2]);
 844     segmentation[1]->Add(id0+20, "Undefined", slatsegNB[2]);
 845     segmentation[1]->Add(id0+ 6, "Undefined", slatsegNB[2]);
 846
 847     // type 222333
 848     segmentation[1]->Add(id0+16, "Undefined", slatsegNB[3]);
 849     segmentation[1]->Add(id0+10, "Undefined", slatsegNB[3]);
 850     segmentation[1]->Add(id0+23, "Undefined", slatsegNB[3]);
 851     segmentation[1]->Add(id0+ 3, "Undefined", slatsegNB[3]);
 852
 853     // type 223330
 854     segmentation[1]->Add(id0+17, "Undefined", slatsegNB[4]);
 855     segmentation[1]->Add(id0+ 9, "Undefined", slatsegNB[4]);
 856     segmentation[1]->Add(id0+22, "Undefined", slatsegNB[4]);
 857     segmentation[1]->Add(id0+ 4, "Undefined", slatsegNB[4]);
 858
 859     // type 333300
 860     segmentation[1]->Add(id0+18, "Undefined", slatsegNB[5]);
 861     segmentation[1]->Add(id0+ 8, "Undefined", slatsegNB[5]);
 862     segmentation[1]->Add(id0+21, "Undefined", slatsegNB[5]);
 863     segmentation[1]->Add(id0+ 5, "Undefined", slatsegNB[5]);
 864     Segmentation()->AddModuleSegmentation(chamber, 1, segmentation[1]);
 865   }
 866 }
 867
 868 //__________________________________________________________________________
 869 void AliMUONSegFactory::BuildStation6()
 870 {
 871   //--------------------------------------------------------
 872   // Configuration for Trigger stations
 873   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 874
 875
 876     AliMUONGeometrySegmentation *chamberSeg[2];
 877
 878     for (Int_t chamber = 10; chamber < 14; chamber++) {
 879
 880       //Trigger Segmentation
 881       AliMUONTriggerSegmentation *trigSegX[9];
 882       AliMUONTriggerSegmentation *trigSegY[9];
 883       for(Int_t i=0; i<9; i++) {
 884         trigSegX[i] = new AliMUONTriggerSegmentation(1);
 885         trigSegY[i] = new AliMUONTriggerSegmentation(0);
 886         Segmentation()->AddDESegmentation(trigSegX[i]);
 887         Segmentation()->AddDESegmentation(trigSegY[i]);
 888         trigSegX[i]->SetLineNumber(9-i);
 889         trigSegY[i]->SetLineNumber(9-i);
 890       }
 891
 892       //AliMUONChamber *iChamber, *iChamber1;
 893       //iChamber1 = &fMUON->Chamber(10);
 894       //iChamber  = &fMUON->Chamber(chamber);
 895       //Float_t zpos1= iChamber1->Z();
 896       //Float_t zpos = iChamber->Z();
 897       Float_t zpos1= AliMUONConstants::DefaultChamberZ(10);
 898       Float_t zpos = AliMUONConstants::DefaultChamberZ(chamber);
 899       Float_t zRatio = zpos / zpos1;
 900
 901       // init
 902       Float_t stripWidth[3]={0.,0.,0.};     // 1.0625 2.125 4.25
 903       Float_t stripLength[4]={0.,0.,0.,0.}; // 17. 34. 51. 68.
 904       for (Int_t i=0; i<3; i++)
 905         stripWidth[i]=AliMUONTriggerConstants::StripWidth(i)*zRatio;
 906       for (Int_t i=0; i<4; i++)
 907         stripLength[i]=AliMUONTriggerConstants::StripLength(i)*zRatio;
 908       Int_t nStrip[7]={0,0,0,0,0,0,0};
 909       Float_t stripYsize[7]={0.,0.,0.,0.,0.,0.,0.};
 910       Float_t stripXsize[7]={0.,0.,0.,0.,0.,0.,0.};
 911
 912       // chamber 8 0 cathode 0
 913       for (Int_t i=0; i<7; i++) nStrip[i]=16;
 914       for (Int_t i=0; i<7; i++) stripYsize[i]=stripWidth[2];
 915       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
 916       stripXsize[6]=stripLength[2];
 917       trigSegX[8]->Init(0,nStrip,stripYsize,stripXsize,0.);
 918       trigSegX[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
 919
 920       // chamber 8 7 1 0 cathode 1
 921       for (Int_t i=0; i<6; i++) nStrip[i]=8;
 922       nStrip[6]=16;
 923       for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
 924       for (Int_t i=0; i<7; i++) stripXsize[i]=stripWidth[2];
 925       trigSegY[8]->Init(0,nStrip,stripYsize,stripXsize,0.);
 926       trigSegY[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
 927       trigSegY[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
 928       trigSegY[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
 929
 930       // chamber 7 6 2 1 cathode 0
 931       for (Int_t i=0; i<6; i++) nStrip[i]=32;
 932       nStrip[6]=16;
 933       for (Int_t i=0; i<6; i++) stripYsize[i]=stripWidth[1];
 934       stripYsize[6]=stripWidth[2];
 935       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
 936       stripXsize[6]=stripLength[2];
 937       trigSegX[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
 938       trigSegX[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
 939       trigSegX[2]->Init(0,nStrip,stripYsize,stripXsize,0.);
 940       trigSegX[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
 941
 942       // chamber 6 2 cathode 1
 943       for (Int_t i=0; i<5; i++) nStrip[i]=16;
 944       for (Int_t i=5; i<6; i++) nStrip[i]=8;
 945       nStrip[6]=16;
 946       for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
 947       for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
 948       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
 949       trigSegY[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
 950       trigSegY[2]->Init(0,nStrip,stripYsize,stripXsize,0.);
 951
 952       // chamber 5 3 cathode 0
 953       nStrip[0]=48;
 954       for (Int_t i=1; i<3; i++) nStrip[i]=64;
 955       for (Int_t i=3; i<6; i++) nStrip[i]=32;
 956       nStrip[6]=16;
 957       for (Int_t i=0; i<3; i++) stripYsize[i]=stripWidth[0];
 958       for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
 959       stripYsize[6]=stripWidth[2];
 960       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
 961       stripXsize[6]=stripLength[2];
 962       trigSegX[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);
 963       trigSegX[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
 964
 965       // chamber 5 3 cathode 1
 966       for (Int_t i=0; i<5; i++) nStrip[i]=16;
 967       for (Int_t i=5; i<6; i++) nStrip[5]=8;
 968       nStrip[6]=16;
 969       stripYsize[0]=stripLength[2];
 970       for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
 971       for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
 972       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
 973       trigSegY[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);
 974       trigSegY[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
 975
 976       // chamber 4 cathode 0
 977       nStrip[0]=0;
 978       for (Int_t i=1; i<3; i++) nStrip[i]=64;
 979       for (Int_t i=3; i<6; i++) nStrip[i]=32;
 980       nStrip[6]=16;
 981       stripYsize[0]=0.;
 982       for (Int_t i=1; i<3; i++) stripYsize[i]=stripWidth[0];
 983       for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
 984       stripYsize[6]=stripWidth[2];
 985       stripXsize[0]=0;
 986       stripXsize[1]=stripLength[0];
 987       for (Int_t i=2; i<6; i++) stripXsize[i]=stripLength[1];
 988       stripXsize[6]=stripLength[2];
 989       trigSegX[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
 990
 991       // chamber 4 cathode 1
 992       nStrip[0]=0;
 993       nStrip[1]=8;
 994       for (Int_t i=2; i<5; i++) nStrip[i]=16;
 995       for (Int_t i=5; i<6; i++) nStrip[i]=8;
 996       nStrip[6]=16;
 997       stripYsize[0]=0.;
 998       for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
 999       stripXsize[0]=0.;
1000       for (Int_t i=1; i<5; i++) stripXsize[i]=stripWidth[1];
1001       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
1002       trigSegY[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
1003
1004       const AliMUONGeometryModuleTransformer* kModuleTransformer
1005         = fkTransformer->GetModuleTransformer(chamber);
1006
1007       chamberSeg[0] = new AliMUONGeometrySegmentation(kModuleTransformer);
1008       chamberSeg[1] = new AliMUONGeometrySegmentation(kModuleTransformer);
1009
1010       Int_t icount=chamber-10;  // chamber counter (0 1 2 3)
1011       Int_t id0=(10+icount+1)*100;
1012
1013
1014       //  printf("in CreateTriggerSegmentation here 0 id0=%i \n",id0);
1015       chamberSeg[0]->Add(id0+0,  "Undefined", trigSegX[4]);
1016       chamberSeg[0]->Add(id0+1,  "Undefined", trigSegX[5]);
1017       chamberSeg[0]->Add(id0+2,  "Undefined", trigSegX[6]);
1018       chamberSeg[0]->Add(id0+3,  "Undefined", trigSegX[7]);
1019       chamberSeg[0]->Add(id0+4,  "Undefined", trigSegX[8]);
1020       chamberSeg[0]->Add(id0+5,  "Undefined", trigSegX[8]);
1021       chamberSeg[0]->Add(id0+6,  "Undefined", trigSegX[7]);
1022       chamberSeg[0]->Add(id0+7,  "Undefined", trigSegX[6]);
1023       chamberSeg[0]->Add(id0+8,  "Undefined", trigSegX[5]);
1024       chamberSeg[0]->Add(id0+9,  "Undefined", trigSegX[4]);
1025       chamberSeg[0]->Add(id0+10, "Undefined", trigSegX[3]);
1026       chamberSeg[0]->Add(id0+11, "Undefined", trigSegX[2]);
1027       chamberSeg[0]->Add(id0+12, "Undefined", trigSegX[1]);
1028       chamberSeg[0]->Add(id0+13, "Undefined", trigSegX[0]);
1029       chamberSeg[0]->Add(id0+14, "Undefined", trigSegX[0]);
1030       chamberSeg[0]->Add(id0+15, "Undefined", trigSegX[1]);
1031       chamberSeg[0]->Add(id0+16, "Undefined", trigSegX[2]);
1032       chamberSeg[0]->Add(id0+17, "Undefined", trigSegX[3]);
1033
1034       chamberSeg[1]->Add(id0+0,  "Undefined", trigSegY[4]);
1035       chamberSeg[1]->Add(id0+1,  "Undefined", trigSegY[5]);
1036       chamberSeg[1]->Add(id0+2,  "Undefined", trigSegY[6]);
1037       chamberSeg[1]->Add(id0+3,  "Undefined", trigSegY[7]);
1038       chamberSeg[1]->Add(id0+4,  "Undefined", trigSegY[8]);
1039       chamberSeg[1]->Add(id0+5,  "Undefined", trigSegY[8]);
1040       chamberSeg[1]->Add(id0+6,  "Undefined", trigSegY[7]);
1041       chamberSeg[1]->Add(id0+7,  "Undefined", trigSegY[6]);
1042       chamberSeg[1]->Add(id0+8,  "Undefined", trigSegY[5]);
1043       chamberSeg[1]->Add(id0+9,  "Undefined", trigSegY[4]);
1044       chamberSeg[1]->Add(id0+10, "Undefined", trigSegY[3]);
1045       chamberSeg[1]->Add(id0+11, "Undefined", trigSegY[2]);
1046       chamberSeg[1]->Add(id0+12, "Undefined", trigSegY[1]);
1047       chamberSeg[1]->Add(id0+13, "Undefined", trigSegY[0]);
1048       chamberSeg[1]->Add(id0+14, "Undefined", trigSegY[0]);
1049       chamberSeg[1]->Add(id0+15, "Undefined", trigSegY[1]);
1050       chamberSeg[1]->Add(id0+16, "Undefined", trigSegY[2]);
1051       chamberSeg[1]->Add(id0+17, "Undefined", trigSegY[3]);
1052
1053       Segmentation()->AddModuleSegmentation(chamber, 0, chamberSeg[0]);
1054       Segmentation()->AddModuleSegmentation(chamber, 1, chamberSeg[1]);
1055
1056       //  printf("in CreateTriggerSegmentation here 1\n");
1057       if (!id0) {
1058         AliWarning(Form("Segmentation for chamber %d is not yet defined",chamber));
1059         return ;
1060       }
1061     }
1062 }
1063