ITS/UPGRADE/AliITSUv0Layer.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 // This class Defines the Geometry for the ITS Upgrade using TGeo
  18 // This is a work class used to study different configurations
  19 // during the development of the new ITS structure.
  20 //
  21 //  Mario Sitta <sitta@to.infn.it>
  22 //  Chinorat Kobdaj (kobdaj@g.sut.ac.th)
  23 //*************************************************************************
  24
  25
  26 /* $Id: AliITSUv0Layer.cxx  */
  27 // General Root includes
  28 #include <TMath.h>
  29 // Root Geometry includes
  30 //#include <AliLog.h>
  31 #include <TGeoManager.h>
  32 #include <TGeoVolume.h>
  33 #include <TGeoPcon.h>
  34 #include <TGeoCone.h>
  35 #include <TGeoTube.h> // contaings TGeoTubeSeg
  36 #include <TGeoArb8.h>
  37 #include <TGeoXtru.h>
  38 #include <TGeoCompositeShape.h>
  39 #include <TGeoMatrix.h>
  40 #include "AliITSUv0Layer.h"
  41 #include "AliITSUGeomTGeo.h"
  42 #include <TGeoBBox.h>
  43 #include <TGeoShape.h>
  44 #include <TGeoTrd1.h>
  45 using namespace TMath;
  46
  47 const Double_t AliITSUv0Layer::fgkDefaultSensorThick = 300*fgkmicron;
  48 const Double_t AliITSUv0Layer::fgkDefaultLadderThick =   1*fgkcm;
  49
  50 ClassImp(AliITSUv0Layer)
  51
  52 #define SQ(A) (A)*(A)
  53
  54 //________________________________________________________________________
  55 AliITSUv0Layer::AliITSUv0Layer():
  56   AliITSv11Geometry(),
  57   fLayerNumber(0),
  58   fPhi0(0),
  59   fLayRadius(0),
  60   fZLength(0),
  61   fSensorThick(0),
  62   fLadderThick(0),
  63   fLadderWidth(0),
  64   fLadderTilt(0),
  65   fNLadders(0),
  66   fNModules(0),
  67   fDetTypeID(0),
  68   fIsTurbo(0),
  69   fBuildLevel(0),
  70   fStaveModel(AliITSUv0::kModelDummy)
  71 {
  72   //
  73   // Standard constructor
  74   //
  75 }
  76
  77 //________________________________________________________________________
  78 AliITSUv0Layer::AliITSUv0Layer(Int_t debug):
  79   AliITSv11Geometry(debug),
  80   fLayerNumber(0),
  81   fPhi0(0),
  82   fLayRadius(0),
  83   fZLength(0),
  84   fSensorThick(0),
  85   fLadderThick(0),
  86   fLadderWidth(0),
  87   fLadderTilt(0),
  88   fNLadders(0),
  89   fNModules(0),
  90   fDetTypeID(0),
  91   fIsTurbo(0),
  92   fBuildLevel(0),
  93   fStaveModel(AliITSUv0::kModelDummy)
  94 {
  95   //
  96   // Constructor setting debugging level
  97   //
  98 }
  99
 100 //________________________________________________________________________
 101 AliITSUv0Layer::AliITSUv0Layer(Int_t lay, Int_t debug):
 102   AliITSv11Geometry(debug),
 103   fLayerNumber(lay),
 104   fPhi0(0),
 105   fLayRadius(0),
 106   fZLength(0),
 107   fSensorThick(0),
 108   fLadderThick(0),
 109   fLadderWidth(0),
 110   fLadderTilt(0),
 111   fNLadders(0),
 112   fNModules(0),
 113   fDetTypeID(0),
 114   fIsTurbo(0),
 115   fBuildLevel(0),
 116   fStaveModel(AliITSUv0::kModelDummy)
 117 {
 118   //
 119   // Constructor setting layer number and debugging level
 120   //
 121 }
 122
 123 //________________________________________________________________________
 124 AliITSUv0Layer::AliITSUv0Layer(Int_t lay, Bool_t turbo, Int_t debug):
 125   AliITSv11Geometry(debug),
 126   fLayerNumber(lay),
 127   fPhi0(0),
 128   fLayRadius(0),
 129   fZLength(0),
 130   fSensorThick(0),
 131   fLadderThick(0),
 132   fLadderWidth(0),
 133   fLadderTilt(0),
 134   fNLadders(0),
 135   fNModules(0),
 136   fDetTypeID(0),
 137   fIsTurbo(turbo),
 138   fBuildLevel(0),
 139   fStaveModel(AliITSUv0::kModelDummy)
 140 {
 141   //
 142   // Constructor setting layer number and debugging level
 143   // for a "turbo" layer (i.e. where ladders overlap in phi)
 144   //
 145 }
 146
 147 //________________________________________________________________________
 148 AliITSUv0Layer::AliITSUv0Layer(const AliITSUv0Layer &s):
 149   AliITSv11Geometry(s.GetDebug()),
 150   fLayerNumber(s.fLayerNumber),
 151   fPhi0(s.fPhi0),
 152   fLayRadius(s.fLayRadius),
 153   fZLength(s.fZLength),
 154   fSensorThick(s.fSensorThick),
 155   fLadderThick(s.fLadderThick),
 156   fLadderWidth(s.fLadderWidth),
 157   fLadderTilt(s.fLadderTilt),
 158   fNLadders(s.fNLadders),
 159   fNModules(s.fNModules),
 160   fDetTypeID(s.fDetTypeID),
 161   fIsTurbo(s.fIsTurbo),
 162   fBuildLevel(s.fBuildLevel),
 163   fStaveModel(s.fStaveModel)
 164 {
 165   //
 166   // Copy constructor
 167   //
 168 }
 169
 170 //________________________________________________________________________
 171 AliITSUv0Layer& AliITSUv0Layer::operator=(const AliITSUv0Layer &s)
 172 {
 173   //
 174   // Assignment operator
 175   //
 176   if(&s == this) return *this;
 177
 178   fLayerNumber = s.fLayerNumber;
 179   fPhi0        = s.fPhi0;
 180   fLayRadius   = s.fLayRadius;
 181   fZLength     = s.fZLength;
 182   fSensorThick = s.fSensorThick;
 183   fLadderThick = s.fLadderThick;
 184   fLadderWidth = s.fLadderWidth;
 185   fLadderTilt  = s.fLadderTilt;
 186   fNLadders    = s.fNLadders;
 187   fNModules    = s.fNModules;
 188   fIsTurbo     = s.fIsTurbo;
 189   fDetTypeID   = s.fDetTypeID;
 190   fBuildLevel  = s.fBuildLevel;
 191   fStaveModel  = s.fStaveModel;
 192
 193   return *this;
 194 }
 195
 196 //________________________________________________________________________
 197 AliITSUv0Layer::~AliITSUv0Layer() {
 198   //
 199   // Destructor
 200   //
 201 }
 202
 203 //________________________________________________________________________
 204 void AliITSUv0Layer::CreateLayer(TGeoVolume *moth){
 205 //
 206 // Creates the actual Layer and places inside its mother volume
 207 //
 208 // Input:
 209 //         moth : the TGeoVolume owing the volume structure
 210 //
 211 // Output:
 212 //
 213 // Return:
 214 //
 215 // Created:      17 Jun 2011  Mario Sitta
 216 // Updated:      08 Jul 2011  Mario Sitta
 217 // Updated:      20 May 2013  Mario Sitta  Layer is Assembly instead of Tube
 218 //
 219   // Local variables
 220   char volname[30];
 221   Double_t xpos, ypos, zpos;
 222   Double_t alpha;
 223
 224
 225   // Check if the user set the proper parameters
 226   if (fLayRadius <= 0) AliFatal(Form("Wrong layer radius (%f)",fLayRadius));
 227   if (fZLength   <= 0) AliFatal(Form("Wrong layer length (%f)",fZLength));
 228   if (fNLadders  <= 0) AliFatal(Form("Wrong number of ladders (%d)",fNLadders));
 229   if (fNModules  <= 0) AliFatal(Form("Wrong number of modules (%d)",fNModules));
 230
 231   if (fLadderThick <= 0) {
 232     AliInfo(Form("Ladder thickness wrong or not set (%f), using default (%f)",
 233                  fLadderThick,fgkDefaultLadderThick));
 234     fLadderThick = fgkDefaultLadderThick;
 235   }
 236
 237   if (fSensorThick <= 0) {
 238     AliInfo(Form("Sensor thickness wrong or not set (%f), using default (%f)",
 239                  fSensorThick,fgkDefaultSensorThick));
 240     fSensorThick = fgkDefaultSensorThick;
 241   }
 242
 243   if (fSensorThick > fLadderThick) {
 244     AliWarning(Form("Sensor thickness (%f) is greater than ladder thickness (%f), fixing",
 245                  fSensorThick,fLadderThick));
 246     fSensorThick = fLadderThick;
 247   }
 248
 249
 250   // If a Turbo layer is requested, do it and exit
 251   if (fIsTurbo) {
 252     CreateLayerTurbo(moth);
 253     return;
 254   }
 255
 256
 257   // First create the ladder container
 258   alpha = (360./(2*fNLadders))*DegToRad();
 259
 260   //  fLadderWidth = fLayRadius*Tan(alpha);
 261
 262   snprintf(volname, 30, "%s%d", AliITSUGeomTGeo::GetITSLayerPattern(),fLayerNumber);
 263   TGeoVolume *layVol = new TGeoVolumeAssembly(volname);
 264   layVol->SetUniqueID(fDetTypeID);
 265
 266 //  layVol->SetVisibility(kFALSE);
 267   layVol->SetVisibility(kTRUE);
 268   layVol->SetLineColor(1);
 269
 270   TGeoVolume *laddVol = CreateLadder();
 271
 272
 273   // Now build up the layer
 274   alpha = 360./fNLadders;
 275   Double_t r = fLayRadius + ((TGeoBBox*)laddVol->GetShape())->GetDY();
 276   for (Int_t j=0; j<fNLadders; j++) {
 277     Double_t phi = j*alpha + fPhi0;
 278     xpos = r*CosD(phi);// r*SinD(-phi);
 279     ypos = r*SinD(phi);// r*CosD(-phi);
 280     zpos = 0.;
 281     phi += 90;
 282     layVol->AddNode(laddVol, j, new TGeoCombiTrans( xpos, ypos, zpos,
 283                                                     new TGeoRotation("",phi,0,0)));
 284   }
 285
 286
 287   // Finally put everything in the mother volume
 288   moth->AddNode(layVol, 1, 0);
 289
 290
 291   // Upgrade geometry is served
 292   return;
 293 }
 294
 295 //________________________________________________________________________
 296 void AliITSUv0Layer::CreateLayerTurbo(TGeoVolume *moth){
 297 //
 298 // Creates the actual Layer and places inside its mother volume
 299 // A so-called "turbo" layer is a layer where ladders overlap in phi
 300 // User can set width and tilt angle, no check is performed here
 301 // to avoid volume overlaps
 302 //
 303 // Input:
 304 //         moth : the TGeoVolume owing the volume structure
 305 //
 306 // Output:
 307 //
 308 // Return:
 309 //
 310 // Created:      08 Jul 2011  Mario Sitta
 311 // Updated:      08 Mar 2012  Mario Sitta  Correct way to compute container R
 312 // Updated:      20 May 2013  Mario Sitta  Layer is Assemgbly instead of Tube
 313 //
 314
 315
 316   // Local variables
 317   char volname[30];
 318   Double_t xpos, ypos, zpos;
 319   Double_t alpha;
 320
 321
 322   // Check if the user set the proper (remaining) parameters
 323   if (fLadderWidth <= 0)
 324     AliFatal(Form("Wrong ladder width (%f)",fLadderWidth));
 325   if (Abs(fLadderTilt) > 45)
 326     AliWarning(Form("Ladder tilt angle (%f) greater than 45deg",fLadderTilt));
 327
 328
 329   snprintf(volname, 30, "%s%d", AliITSUGeomTGeo::GetITSLayerPattern(), fLayerNumber);
 330   TGeoVolume *layVol = new TGeoVolumeAssembly(volname);
 331   layVol->SetUniqueID(fDetTypeID);
 332   layVol->SetVisibility(kTRUE);
 333   layVol->SetLineColor(1);
 334   TGeoVolume *laddVol = CreateLadder();
 335
 336
 337   // Now build up the layer
 338   alpha = 360./fNLadders;
 339   Double_t r = fLayRadius /* +module thick ?! */;
 340   for (Int_t j=0; j<fNLadders; j++) {
 341     Double_t phi = j*alpha + fPhi0;
 342     xpos = r*CosD(phi);// r*SinD(-phi);
 343     ypos = r*SinD(phi);// r*CosD(-phi);
 344     zpos = 0.;
 345     phi += 90;
 346     layVol->AddNode(laddVol, j, new TGeoCombiTrans( xpos, ypos, zpos,
 347                                                     new TGeoRotation("", phi-fLadderTilt,0,0)));
 348   }
 349
 350
 351   // Finally put everything in the mother volume
 352   moth->AddNode(layVol, 1, 0);
 353
 354   return;
 355 }
 356
 357 //________________________________________________________________________
 358 TGeoVolume* AliITSUv0Layer::CreateLadder(const TGeoManager * /*mgr*/){
 359 //
 360 // Creates the actual Ladder
 361 //
 362 // Input:
 363 //         mgr  : the GeoManager (used only to get the proper material)
 364 //
 365 // Output:
 366 //
 367 // Return:
 368 //
 369 // Created:      22 Jun 2011  Mario Sitta
 370 //
 371
 372   char volname[30];
 373   Double_t xlen, ylen, zlen;
 374   Double_t xpos, ypos, zpos, zmod;
 375   Double_t alpha;
 376
 377
 378   // First create all needed shapes
 379   alpha = (360./(2*fNLadders))*DegToRad();
 380
 381   // The ladder
 382   xlen = fLayRadius*Tan(alpha);
 383   if (fIsTurbo) xlen = 0.5*fLadderWidth;
 384   ylen = 0.5*fLadderThick;
 385   zlen = 0.5*fZLength;
 386
 387   Double_t yplus = 0.46;
 388   //  Double_t origin[]={0,-yplus/2,0};
 389   //  TGeoBBox *ladder = new TGeoBBox(xlen, ylen+yplus/2, zlen, origin);
 390   TGeoXtru *ladder = new TGeoXtru(2); //z sections
 391   Double_t xv[5] = {xlen,xlen,0,-xlen,-xlen};
 392   Double_t yv[5] = {ylen+0.09,-0.15,-yplus-fSensorThick,-0.15,ylen+0.09};
 393   ladder->DefinePolygon(5,xv,yv);
 394   ladder->DefineSection(0,-zlen,0,0,1.);
 395   ladder->DefineSection(1,+zlen,0,0,1.);
 396
 397
 398   // We have all shapes: now create the real volumes
 399 //  TGeoMedium *medAir = mgr->GetMedium("ITS_AIR$");
 400
 401   snprintf(volname, 30, "%s%d", AliITSUGeomTGeo::GetITSLadderPattern(), fLayerNumber);
 402 //  TGeoVolume *laddVol = new TGeoVolume(volname, ladder, medAir);
 403   TGeoVolume *laddVol = new TGeoVolumeAssembly(volname);
 404
 405   //  laddVol->SetVisibility(kFALSE);
 406   laddVol->SetVisibility(kTRUE);
 407   laddVol->SetLineColor(2);
 408   TGeoVolume *modVol = CreateModule(xlen,ylen, zlen);
 409
 410
 411   // Now build up the ladder
 412   zmod = ((TGeoBBox*)modVol->GetShape())->GetDZ();
 413   for (Int_t j=0; j<fNModules; j++) {
 414     xpos = 0.;
 415     ypos = 0.021;  // Remove small overlap - M.S: 21may13
 416     zpos = -ladder->GetDZ() + j*2*zmod + zmod;
 417     laddVol->AddNode(modVol, j, new TGeoTranslation(xpos, ypos, zpos));
 418   }
 419
 420
 421   // put mechanical stave structure, only inner barrel up to now
 422   if (fLayerNumber<3) {
 423     TGeoVolume *mechLaddVol = CreateStaveStruct(xlen,zlen);
 424     if (mechLaddVol)
 425       laddVol->AddNode(mechLaddVol, fNModules, new TGeoCombiTrans(0, -0.15-ylen, 0, new TGeoRotation("",0, 0, 180)));
 426   }
 427
 428
 429   // Done, return the ladder
 430   return laddVol;
 431 }
 432
 433 //________________________________________________________________________
 434 TGeoVolume* AliITSUv0Layer::CreateStaveStruct(const Double_t xlad,
 435                                               const Double_t zlad,
 436                                               const TGeoManager *mgr){
 437 //
 438 // Create the mechanical stave structure
 439 //
 440 // Input:
 441 //         xlad : X length
 442 //         zlad : Z length
 443 //         mgr  : the GeoManager (used only to get the proper material)
 444 //
 445 // Output:
 446 //
 447 // Return:
 448 //
 449 // Created:      22 Mar 2013  Chinorat Kobdaj
 450 // Updated:      26 Apr 2013  Mario Sitta
 451 //
 452
 453   TGeoVolume *mechLaddVol = 0;
 454
 455   switch (fStaveModel) {
 456     case AliITSUv0::kModelDummy:
 457       mechLaddVol = CreateStaveModelDummy(xlad,zlad,mgr);
 458       break;
 459     case AliITSUv0::kModel0:
 460       mechLaddVol = CreateStaveModel0(xlad,zlad,mgr);
 461       break;
 462     case AliITSUv0::kModel1:
 463       mechLaddVol = CreateStaveModel1(xlad,zlad,mgr);
 464       break;
 465     case AliITSUv0::kModel21:
 466       mechLaddVol = CreateStaveModel21(xlad,zlad,mgr);
 467       break;
 468     case AliITSUv0::kModel22:
 469       mechLaddVol = CreateStaveModel22(xlad,zlad,mgr);
 470       break;
 471     case AliITSUv0::kModel3:
 472       mechLaddVol = CreateStaveModel3(xlad,zlad,mgr);
 473       break;
 474     default:
 475       AliFatal(Form("Unknown stave model %d",fStaveModel));
 476       break;
 477   }
 478
 479   return mechLaddVol;
 480 }
 481
 482
 483 //________________________________________________________________________
 484 TGeoVolume* AliITSUv0Layer::CreateStaveModelDummy(const Double_t ,
 485                                               const Double_t ,
 486                                               const TGeoManager *) const {
 487 //
 488 // Create dummy stave
 489 //
 490 // Input:
 491 //         xlad : X length
 492 //         zlad : Z length
 493 //         mgr  : the GeoManager (used only to get the proper material)
 494 //
 495 // Output:
 496 //
 497 // Return:
 498 //
 499 // Created:      22 Mar 2013  Chinorat Kobdaj
 500 // Updated:      26 Apr 2013  Mario Sitta
 501 //
 502
 503   // Done, return the stave structur
 504   return 0;
 505 }
 506
 507 //________________________________________________________________________
 508 TGeoVolume* AliITSUv0Layer::CreateStaveModel0(const Double_t xlad,
 509                                               const Double_t zlad,
 510                                               const TGeoManager *mgr){
 511 //
 512 // Create the mechanical stave structure for Model 0 of TDR
 513 //
 514 // Input:
 515 //         xlad : X length
 516 //         zlad : Z length
 517 //         mgr  : the GeoManager (used only to get the proper material)
 518 //
 519 // Output:
 520 //
 521 // Return:
 522 //
 523 // Created:      22 Mar 2013  Chinorat Kobdaj
 524 // Updated:      26 Apr 2013  Mario Sitta
 525 //
 526
 527   // Materials defined in AliITSUv0
 528   TGeoMedium *medAir    = mgr->GetMedium("ITS_AIR$");
 529   TGeoMedium *medWater  = mgr->GetMedium("ITS_WATER$");
 530
 531   TGeoMedium *medM60J3K    = mgr->GetMedium("ITS_M60J3K$");
 532   TGeoMedium *medKapton    = mgr->GetMedium("ITS_KAPTON(POLYCH2)$");
 533   TGeoMedium *medGlue      = mgr->GetMedium("ITS_GLUE$");
 534   TGeoMedium *medFlexCable = mgr->GetMedium("ITS_FLEXCABLE$");
 535
 536   // Local parameters
 537   Double_t kConeOutRadius = 0.15/2;
 538   Double_t kConeInRadius = 0.1430/2;
 539   Double_t kStaveLength = zlad*2;
 540   Double_t kStaveWidth = xlad*2-kConeOutRadius*2;
 541   Double_t kWidth = kStaveWidth/4;//1/2 of kWidth
 542   Double_t kStaveHeight = 0.3;
 543   Double_t kHeight = kStaveHeight/2;
 544   Double_t kAlpha = 90-67;//90-33.69;
 545   Double_t kTheta = kAlpha*TMath::DegToRad();
 546   Double_t kS1 = kWidth/TMath::Sin(kTheta);
 547   Double_t kL1 = kWidth/TMath::Tan(kTheta);
 548   Double_t kS2 = TMath::Sqrt(kHeight*kHeight + kS1*kS1);//TMath::Sin(the2);
 549   Double_t kThe2 = TMath::ATan(kHeight/kS1);
 550   Double_t kBeta = kThe2*TMath::RadToDeg();
 551   // Int_t  loop = kStaveLength/(kL1);
 552   // Double_t s3 = kWidth/(2*TMath::Sin(kTheta));
 553   // Double_t s4 = 3*kWidth/(2*TMath::Sin(kTheta));
 554
 555   AliDebug(1, Form("BuildLevel %d\n",fBuildLevel));
 556
 557   char volname[30];
 558   snprintf(volname, 30, "%s%d_StaveStruct", AliITSUGeomTGeo::GetITSLadderPattern(), fLayerNumber);
 559
 560   Double_t z=0, y=-0.011+0.0150, x=0;
 561
 562    TGeoVolume *mechLaddVol = 0;
 563
 564   if (fBuildLevel < 5) {
 565
 566     // world (trapezoid)
 567     TGeoXtru *mechStruct = new TGeoXtru(2); //z sections
 568     Double_t xv[5] = {kStaveWidth/2+0.1,kStaveWidth/2+0.1,0,-kStaveWidth/2-0.1,-kStaveWidth/2-0.1};
 569     Double_t yv[5] = {-kConeOutRadius*2-0.07,0,kStaveHeight,0,-kConeOutRadius*2-0.07};
 570     mechStruct->DefinePolygon(5,xv,yv);
 571     mechStruct->DefineSection(0,-kStaveLength-0.1,0,0,1.);
 572     mechStruct->DefineSection(1,kStaveLength+0.1,0,0,1.);
 573
 574     mechLaddVol = new TGeoVolume(volname, mechStruct, medAir);
 575     mechLaddVol->SetLineColor(12);
 576     mechLaddVol->SetFillColor(12);
 577     mechLaddVol->SetVisibility(kTRUE);
 578
 579     // detailed structure ++++++++++++++
 580     //Pipe Kapton grey-35
 581     TGeoTube *coolTube = new TGeoTube(kConeInRadius,kConeOutRadius,kStaveLength/2);
 582     TGeoVolume *volCoolTube= new TGeoVolume("pipe", coolTube, medKapton);
 583     volCoolTube->SetFillColor(35);
 584     volCoolTube->SetLineColor(35);
 585     mechLaddVol->AddNode(volCoolTube,0,new TGeoTranslation(x+(kStaveWidth/2),y-(kHeight-kConeOutRadius),0));
 586     mechLaddVol->AddNode(volCoolTube,1,new TGeoTranslation(x-(kStaveWidth/2),y-(kHeight-kConeOutRadius),0));
 587   }
 588
 589   if (fBuildLevel < 4) {
 590     TGeoTube *coolTubeW = new TGeoTube(0.,kConeInRadius,kStaveLength/2);
 591     TGeoVolume *volCoolTubeW= new TGeoVolume("pipeWater", coolTubeW, medWater);
 592     volCoolTubeW->SetFillColor(4);
 593     volCoolTubeW->SetLineColor(4);
 594     mechLaddVol->AddNode(volCoolTubeW,0,new TGeoTranslation(x+(kStaveWidth/2),y-(kHeight-kConeOutRadius),0));
 595     mechLaddVol->AddNode(volCoolTubeW,1,new TGeoTranslation(x-(kStaveWidth/2),y-(kHeight-kConeOutRadius),0));
 596   }
 597
 598   //frequency of filament
 599   //n = 4 means very dense(4 filaments per interval)
 600   //n = 2 means dense(2 filaments per interval)
 601   Int_t n =4;
 602   Int_t loop = (Int_t)(kStaveLength/(4*kL1/n) + 2/n)-1;
 603   if (fBuildLevel < 3) {
 604     //Top CFRP Filament black-12 Carbon structure TGeoBBox (length,thickness,width)
 605     TGeoBBox *t2=new TGeoBBox(kS2,0.007/2,0.15/2);//(kS2,0.002,0.02);
 606     TGeoVolume *volT2=new TGeoVolume("TopFilament", t2, medM60J3K);
 607     volT2->SetLineColor(12);
 608     volT2->SetFillColor(12);
 609
 610     for(int i=1;i<loop;i++){  //i<60;i++){
 611       mechLaddVol->AddNode(volT2,4*i+0,
 612                                   new TGeoCombiTrans(x+kWidth,y+(2*kConeOutRadius),z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2,//z-14.25+(i*2*kL1),
 613                                                      new TGeoRotation("volT2",90,90-kAlpha,90-kBeta)));
 614       mechLaddVol->AddNode(volT2,4*i+1,
 615                                   new TGeoCombiTrans(x-kWidth,y+(2*kConeOutRadius),z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2,//z-14.25+(i*2*kL1),
 616                                                      new TGeoRotation("volT2",90,-90+kAlpha,-90+kBeta)));
 617       mechLaddVol->AddNode(volT2,4*i+2,
 618                                   new TGeoCombiTrans(x+kWidth,y+(2*kConeOutRadius),z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2,//z-14.25+(i*2*kL1),
 619                                                      new TGeoRotation("volT2",90,-90+kAlpha,90-kBeta)));
 620       mechLaddVol->AddNode(volT2,4*i+3,
 621                                   new TGeoCombiTrans(x-kWidth,y+(2*kConeOutRadius),z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2,//z-14.25+(i*2*kL1),
 622                                                      new TGeoRotation("volT2",90,90-kAlpha,-90+kBeta)));
 623     }
 624
 625
 626     //Bottom CFRP Filament black-12 Carbon structure  TGeoBBox (thickness,width,length)
 627     TGeoBBox *t1=new TGeoBBox(0.007/2,0.15/2,kS1);//(0.002,0.02,kS1);
 628     TGeoVolume *volT1=new TGeoVolume("CFRPBottom", t1, medM60J3K);
 629     volT1->SetLineColor(12);
 630     volT1->SetFillColor(12);
 631
 632     for(int i=1;i<loop;i++){
 633       mechLaddVol->AddNode(volT1,4*i+0,
 634                                   new TGeoCombiTrans(x+kWidth,y-kHeight,z-kStaveLength/2+((4/n)*kL1*i)+kS1/2, //z-14.25+(i*2*kL1),
 635                                                      new TGeoRotation("volT1",-90,kAlpha,0)));
 636       mechLaddVol->AddNode(volT1,4*i+1,
 637                                   new TGeoCombiTrans(x-kWidth,y-kHeight,z-kStaveLength/2+((4/n)*kL1*i)+kS1/2,  //z-14.25+(i*2*kL1),
 638                                                      new TGeoRotation("volT1",90,kAlpha,0)));
 639       mechLaddVol->AddNode(volT1,4*i+2,
 640                                   new TGeoCombiTrans(x+kWidth,y-kHeight,z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2, //z-14.25+(i*2*kL1),
 641                                                      new TGeoRotation("volT1",-90,-kAlpha,0)));
 642       mechLaddVol->AddNode(volT1,4*i+3,
 643                                   new TGeoCombiTrans(x-kWidth,y-kHeight,z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2, //z-14.25+(i*2*kL1),
 644                                                      new TGeoRotation("volT1",-90,+kAlpha,0)));
 645     }
 646   }
 647
 648   if (fBuildLevel < 2) {
 649     // Glue CFRP-Silicon layers TGeoBBox(thickness,width,kS1);
 650     TGeoBBox *tG=new TGeoBBox(0.0075/2,0.18/2,kS1);
 651     TGeoVolume *volTG=new TGeoVolume("Glue1", tG, medGlue);
 652     volTG->SetLineColor(5);
 653     volTG->SetFillColor(5);
 654
 655     for(int i=1;i<loop;i++){ //i<60;i++){
 656       mechLaddVol->AddNode(volTG,4*i+0,
 657                                   new TGeoCombiTrans(x+kWidth,y-0.16,z-kStaveLength/2+((4/n)*kL1*i)+kS1/2, //z-14.25+(2*kL1*i),
 658                                                      new TGeoRotation("volTG",-90,kAlpha,0)));
 659       mechLaddVol->AddNode(volTG,4*i+1,
 660                                   new TGeoCombiTrans(x-kWidth,y-0.16,z-kStaveLength/2+((4/n)*kL1*i)+kS1/2, //z-14.25+(2*kL1*i),
 661                                                      new TGeoRotation("volTG",90,kAlpha,0)));
 662       mechLaddVol->AddNode(volTG,4*i+2,
 663                                   new TGeoCombiTrans(x+kWidth,y-0.16,z-kStaveLength/2+((4/n)*i*kL1)+kS1/2, //z-14.25+(i*2*kL1),
 664                                                      new TGeoRotation("volTG",-90,-kAlpha,0)));
 665       mechLaddVol->AddNode(volTG,4*i+3,
 666                                   new TGeoCombiTrans(x-kWidth,y-0.16,z-kStaveLength/2+(i*(4/n)*kL1)+kS1/2, //z-14.25+(i*2*kL1),
 667                                                      new TGeoRotation("volTG",-90,+kAlpha,0)));
 668     }
 669
 670     TGeoBBox *glue = new TGeoBBox(xlad, 0.005/2, zlad);
 671     TGeoVolume *volGlue=new TGeoVolume("Glue2", glue, medGlue);
 672     volGlue->SetLineColor(5);
 673     volGlue->SetFillColor(5);
 674     //mechLaddVol->AddNode(volGlue, 0, new TGeoCombiTrans(x, y-0.16, z, new TGeoRotation("",0, 0, 0)));
 675     mechLaddVol->AddNode(volGlue, 1, new TGeoCombiTrans(x, y-0.165-fSensorThick-0.005, z, new TGeoRotation("",0, 0, 0)));
 676   }
 677
 678   if (fBuildLevel < 1) {
 679     //Flex cable brown-28 TGeoBBox(width,thickness,length);
 680     TGeoBBox *kapCable = new TGeoBBox(xlad, 0.01/2, zlad);
 681     TGeoVolume *volCable=new TGeoVolume("FlexCable", kapCable, medFlexCable);
 682     volCable->SetLineColor(28);
 683     volCable->SetFillColor(28);
 684     mechLaddVol->AddNode(volCable, 0, new TGeoCombiTrans(x, y-0.165-fSensorThick-0.005-0.01, z, new TGeoRotation("",0, 0, 0)));
 685  }
 686
 687   // Done, return the stave structur
 688   return mechLaddVol;
 689 }
 690
 691
 692 //________________________________________________________________________
 693 TGeoVolume* AliITSUv0Layer::CreateStaveModel1(const Double_t xlad,
 694                                               const Double_t zlad,
 695                                               const TGeoManager *mgr){
 696 //
 697 // Create the mechanical stave structure for Model 1 of TDR
 698 //
 699 // Input:
 700 //         xlad : X length
 701 //         zlad : Z length
 702 //         mgr  : the GeoManager (used only to get the proper material)
 703 //
 704 // Output:
 705 //
 706 // Return:
 707 //
 708 // Created:      22 Mar 2013  Chinorat Kobdaj
 709 // Updated:      26 Apr 2013  Mario Sitta
 710 //
 711
 712   // Materials defined in AliITSUv0
 713   TGeoMedium *medAir    = mgr->GetMedium("ITS_AIR$");
 714   TGeoMedium *medWater  = mgr->GetMedium("ITS_WATER$");
 715
 716   TGeoMedium *medM60J3K    = mgr->GetMedium("ITS_M60J3K$");
 717   TGeoMedium *medKapton    = mgr->GetMedium("ITS_KAPTON(POLYCH2)$");
 718   TGeoMedium *medGlue      = mgr->GetMedium("ITS_GLUE$");
 719   TGeoMedium *medFlexCable = mgr->GetMedium("ITS_FLEXCABLE$");
 720
 721   // Local parameters
 722   Double_t kConeOutRadius = 0.15/2;
 723   //    Double_t kConeInRadius = 0.1430/2;
 724   Double_t kStaveLength = zlad*2;
 725   //    Double_t kStaveWidth = xlad*2-kConeOutRadius*2;
 726   Double_t kStaveWidth = xlad*2;
 727   Double_t kWidth = kStaveWidth/4;//1/2 of kWidth
 728   Double_t kStaveHeight = 0.3;
 729   Double_t kHeight = kStaveHeight/2;
 730   Double_t kAlpha = 90-33.;//90-30;
 731   Double_t kTheta = kAlpha*TMath::DegToRad();
 732   Double_t kS1 = kWidth/TMath::Sin(kTheta);
 733   Double_t kL1 = kWidth/TMath::Tan(kTheta);
 734   Double_t kS2 = TMath::Sqrt(kHeight*kHeight + kS1*kS1);//TMath::Sin(the2);
 735   Double_t kThe2 = TMath::ATan(kHeight/kS1);
 736   Double_t kBeta = kThe2*TMath::RadToDeg();
 737   Int_t  loop = (Int_t)((kStaveLength/(2*kL1))/2);
 738
 739
 740   TGeoVolume *mechLaddVol = 0;
 741
 742   char volname[30];
 743   snprintf(volname, 30, "%s%d_StaveStruct", AliITSUGeomTGeo::GetITSLadderPattern(), fLayerNumber);
 744
 745
 746   // detailed structure ++++++++++++++
 747   Double_t z=0, y=-0.011+0.0150, x=0;
 748
 749   // Polimide micro channels numbers
 750   Double_t yMC = y-kHeight+0.01;
 751   Int_t nb = (Int_t)(kStaveWidth/0.1)+1;
 752   Double_t xladMC = (nb*0.1-0.08)/2;
 753
 754
 755   if (fBuildLevel < 5) {
 756     // world (trapezoid)
 757     TGeoXtru *mechStruct = new TGeoXtru(2); //z sections
 758     Double_t xv[5] = {kStaveWidth/2+0.1,kStaveWidth/2+0.1,0,-kStaveWidth/2-0.1,-kStaveWidth/2-0.1};
 759     Double_t yv[5] = {-kConeOutRadius*2-0.07,0,kStaveHeight,0,-kConeOutRadius*2-0.07};
 760     mechStruct->DefinePolygon(5,xv,yv);
 761     mechStruct->DefineSection(0,-kStaveLength-0.1,0,0,1.);
 762     mechStruct->DefineSection(1,kStaveLength+0.1,0,0,1.);
 763
 764     mechLaddVol = new TGeoVolume(volname, mechStruct, medAir);
 765     mechLaddVol->SetLineColor(12);
 766     mechLaddVol->SetFillColor(12);
 767     mechLaddVol->SetVisibility(kTRUE);
 768
 769     // Polimide micro channels numbers
 770     TGeoBBox *tM0=new TGeoBBox(xladMC, 0.005/2, zlad);
 771     TGeoVolume *volTM0=new TGeoVolume("MicroChanCover", tM0, medKapton);
 772     volTM0->SetLineColor(35);
 773     volTM0->SetFillColor(35);
 774     mechLaddVol->AddNode(volTM0, 0, new TGeoCombiTrans(x,-0.0125+yMC, z, new TGeoRotation("",0, 0, 0)));
 775     mechLaddVol->AddNode(volTM0, 1, new TGeoCombiTrans(x,+0.0125+yMC, z, new TGeoRotation("",0, 0, 0)));
 776
 777     TGeoBBox *tM0b=new TGeoBBox(0.02/2, 0.02/2, zlad);
 778     TGeoVolume *volTM0b=new TGeoVolume("MicroChanWalls", tM0b, medKapton);
 779     volTM0b->SetLineColor(35);
 780     volTM0b->SetFillColor(35);
 781     for (Int_t ib=0;ib<nb;ib++) {
 782       mechLaddVol->AddNode(volTM0b, ib, new TGeoCombiTrans(x+ib*0.1-xladMC+0.01,yMC, z, new TGeoRotation("",0, 0, 0)));
 783     }
 784
 785   }
 786
 787   if (fBuildLevel < 4) {
 788     // Water in Polimide micro channels
 789     TGeoBBox *water=new TGeoBBox(0.08/2, 0.02/2, zlad+0.1);
 790     TGeoVolume *volWater=new TGeoVolume("Water", water, medWater);
 791     volWater->SetLineColor(4);
 792     volWater->SetFillColor(4);
 793     for (Int_t ib=0;ib<(nb-1);ib++) {
 794       mechLaddVol->AddNode(volWater, ib, new TGeoCombiTrans(x+ib*0.1-xladMC+0.06,yMC, z, new TGeoRotation("",0, 0, 0)));
 795     }
 796   }
 797
 798   if (fBuildLevel < 3) {
 799     //Bottom filament CFRP black-12 Carbon structure TGeoBBox (thickness,width,length)
 800     Double_t filWidth = 0.04;
 801     Double_t filHeight= 0.02;
 802     TGeoBBox *t1=new TGeoBBox(filHeight/2,filWidth/2,kS1);
 803     TGeoVolume *volT1=new TGeoVolume("CFRPBottom", t1, medM60J3K);
 804     volT1->SetLineColor(12);
 805     volT1->SetFillColor(12);
 806     for(int i=0;i<loop;i++){//i<30;i++){
 807       mechLaddVol->AddNode(volT1,4*i+0,
 808                                   new TGeoCombiTrans(x+kWidth,y-kHeight+0.04+filHeight/2,z-kStaveLength/2+(4*kL1)+kS1/2,
 809                                                      new TGeoRotation("volT1",-90,kAlpha,0)));
 810       mechLaddVol->AddNode(volT1,4*i+1,
 811                                   new TGeoCombiTrans(x-kWidth,y-kHeight+0.04+filHeight/2,z-kStaveLength/2+(4*kL1*i)+kS1/2,
 812                                                      new TGeoRotation("volT1",90,kAlpha,0)));
 813       mechLaddVol->AddNode(volT1,4*i+2,
 814                                   new TGeoCombiTrans(x+kWidth,y-kHeight+0.04+filHeight/2,z-kStaveLength/2+2*kL1+(i*4*kL1)+kS1/2,
 815                                                      new TGeoRotation("volT1",-90,-kAlpha,0)));
 816       mechLaddVol->AddNode(volT1,4*i+3,
 817                                   new TGeoCombiTrans(x-kWidth,y-kHeight+0.04+filHeight/2,z-kStaveLength/2+2*kL1+(i*4*kL1)+kS1/2,
 818                                                      new TGeoRotation("volT1",-90,+kAlpha,0)));
 819     }
 820
 821       // Top filament CFRP black-12 Carbon structure TGeoBBox (length,thickness,width)
 822     TGeoBBox *t2=new TGeoBBox(kS2,filHeight/2,filWidth/2);
 823     TGeoVolume *volT2=new TGeoVolume("CFRPTop", t2, medM60J3K);
 824     volT2->SetLineColor(12);
 825     volT2->SetFillColor(12);
 826     for(int i=0;i<loop;i++){ //i<30;i++){
 827       mechLaddVol->AddNode(volT2,4*i+0,
 828                                   new TGeoCombiTrans(x+kWidth,y+0.04+filHeight/2,z-kStaveLength/2+(i*4*kL1)+kS1/2,
 829                                                      new TGeoRotation("volT2",90,90-kAlpha,90-kBeta)));
 830       mechLaddVol->AddNode(volT2,4*i+1,
 831                                   new TGeoCombiTrans(x-kWidth,y+0.04+filHeight/2,z-kStaveLength/2+(i*4*kL1)+kS1/2,
 832                                                      new TGeoRotation("volT2",90,-90+kAlpha,-90+kBeta)));
 833       mechLaddVol->AddNode(volT2,4*i+2,
 834                                   new TGeoCombiTrans(x+kWidth,y+0.04+filHeight/2,z-kStaveLength/2+2*kL1+(i*4*kL1)+kS1/2,
 835                                                      new TGeoRotation("volT2",90,-90+kAlpha,90-kBeta)));
 836       mechLaddVol->AddNode(volT2,4*i+3,
 837                                   new TGeoCombiTrans(x-kWidth,y+0.04+filHeight/2,z-kStaveLength/2+2*kL1+(i*4*kL1)+kS1/2,
 838                                                      new TGeoRotation("volT2",90,90-kAlpha,-90+kBeta)));
 839     }
 840   }
 841
 842   if (fBuildLevel < 2) {
 843     // Glue between filament and polimide micro channel
 844     TGeoBBox *t3=new TGeoBBox(0.01/2,0.04,kS1);
 845     TGeoVolume *volT3=new TGeoVolume("FilamentGlue", t3, medGlue);
 846     volT3->SetLineColor(5);
 847     volT3->SetFillColor(5);
 848     for(int i=0;i<loop;i++){//i<30;i++){
 849       mechLaddVol->AddNode(volT3,4*i+0,
 850                                   new TGeoCombiTrans(x+kWidth,y-kHeight+0.0325,z-kStaveLength/2+(4*kL1*i)+kS1/2,
 851                                                      new TGeoRotation("volT1",-90,kAlpha,0)));
 852       mechLaddVol->AddNode(volT3,4*i+1,
 853                                   new TGeoCombiTrans(x-kWidth,y-kHeight+0.0325,z-kStaveLength/2+(4*kL1*i)+kS1/2,
 854                                                      new TGeoRotation("volT1",90,kAlpha,0)));
 855       mechLaddVol->AddNode(volT3,4*i+2,
 856                                   new TGeoCombiTrans(x+kWidth,y-kHeight+0.0325,z-kStaveLength/2+2*kL1+(i*4*kL1)+kS1/2,
 857                                                      new TGeoRotation("volT1",-90,-kAlpha,0)));
 858       mechLaddVol->AddNode(volT3,4*i+3,
 859                                   new TGeoCombiTrans(x-kWidth,y-kHeight+0.0325,z-kStaveLength/2+2*kL1+(i*4*kL1)+kS1/2,
 860                                                      new TGeoRotation("volT1",-90,+kAlpha,0)));
 861     }
 862
 863     // Glue microchannel and sensor
 864     TGeoBBox *glueM = new TGeoBBox(xlad, 0.01/2, zlad);
 865     TGeoVolume *volGlueM=new TGeoVolume("MicroChanGlue", glueM, medGlue);
 866     volGlueM->SetLineColor(5);
 867     volGlueM->SetFillColor(5);
 868     mechLaddVol->AddNode(volGlueM, 0, new TGeoCombiTrans(x, y-0.16, z, new TGeoRotation("",0, 0, 0)));
 869
 870     // Glue sensor and kapton
 871     TGeoBBox *glue = new TGeoBBox(xlad, 0.005/2, zlad);
 872     TGeoVolume *volGlue=new TGeoVolume("SensorGlue", glue, medGlue);
 873     volGlue->SetLineColor(5);
 874     volGlue->SetFillColor(5);
 875     mechLaddVol->AddNode(volGlue, 1, new TGeoCombiTrans(x, y-0.165-fSensorThick-0.005, z, new TGeoRotation("",0, 0, 0)));
 876   }
 877
 878   if (fBuildLevel < 1) {
 879     TGeoBBox *kapCable = new TGeoBBox(xlad, 0.01/2, zlad);
 880     TGeoVolume *volCable=new TGeoVolume("FlexCable", kapCable, medFlexCable);
 881     volCable->SetLineColor(28);
 882     volCable->SetFillColor(28);
 883     mechLaddVol->AddNode(volCable, 0, new TGeoCombiTrans(x, y-0.165-fSensorThick-0.005-0.01, z, new TGeoRotation("",0, 0, 0)));
 884   }
 885
 886   // Done, return the stave structur
 887   return mechLaddVol;
 888
 889 }
 890
 891 //________________________________________________________________________
 892 TGeoVolume* AliITSUv0Layer::CreateStaveModel21(const Double_t xlad,
 893                                                const Double_t zlad,
 894                                                const TGeoManager *mgr){
 895 //
 896 // Create the mechanical stave structure for Model 2.1 of TDR
 897 //
 898 // Input:
 899 //         xlad : X length
 900 //         zlad : Z length
 901 //         mgr  : the GeoManager (used only to get the proper material)
 902 //
 903 // Output:
 904 //
 905 // Return:
 906 //
 907 // Created:      22 Mar 2013  Chinorat Kobdaj
 908 // Updated:      26 Apr 2013  Mario Sitta
 909 //
 910
 911   // Materials defined in AliITSUv0
 912   TGeoMedium *medAir    = mgr->GetMedium("ITS_AIR$");
 913   TGeoMedium *medWater  = mgr->GetMedium("ITS_WATER$");
 914
 915   TGeoMedium *medM60J3K    = mgr->GetMedium("ITS_M60J3K$");
 916   TGeoMedium *medKapton    = mgr->GetMedium("ITS_KAPTON(POLYCH2)$");
 917   TGeoMedium *medGlue      = mgr->GetMedium("ITS_GLUE$");
 918   TGeoMedium *medFlexCable = mgr->GetMedium("ITS_FLEXCABLE$");
 919   TGeoMedium *medK13D2U2k  = mgr->GetMedium("ITS_K13D2U2k$");
 920   TGeoMedium *medFGS003    = mgr->GetMedium("ITS_FGS003$");
 921   TGeoMedium *medCarbonFleece = mgr->GetMedium("ITS_CarbonFleece$");
 922
 923   // Local parameters
 924   Double_t kConeOutRadius =0.151384/2;
 925   Double_t kConeInRadius = 0.145034/2;
 926   Double_t kStaveLength = zlad;
 927   Double_t kStaveWidth = xlad*2;
 928   Double_t kWidth = (kStaveWidth+0.005)/4;
 929   Double_t kStaveHeigth = 0.33;//0.33;
 930   Double_t kHeight = (kStaveHeigth+0.025)/2;
 931   Double_t kAlpha = 57; //56.31;
 932   Double_t kTheta = kAlpha*TMath::DegToRad();
 933   Double_t kS1 = (kStaveWidth/4)/TMath::Sin(kTheta);
 934   Double_t kL1 = (kStaveWidth/4)/TMath::Tan(kTheta);
 935   Double_t kS2 = sqrt(kHeight*kHeight + kS1*kS1);//TMath::Sin(the2);
 936   Double_t kThe2 = TMath::ATan(kHeight/kS1);
 937   Double_t kBeta = kThe2*TMath::RadToDeg();
 938   // Double_t lay1 = 0.003157;
 939   Double_t kLay1 = 0.003;//Amec carbon
 940   // Double_t lay2 = 0.0043215;//C Fleece carbon
 941   Double_t kLay2 = 0.002;//C Fleece carbon
 942   Double_t kLay3 = 0.007;//K13D2U carbon
 943   Int_t  loop = (Int_t)(kStaveLength/(2*kL1));
 944
 945
 946   char volname[30];
 947   snprintf(volname, 30, "%s%d_StaveStruct", AliITSUGeomTGeo::GetITSLadderPattern(), fLayerNumber);
 948
 949   Double_t z=0, y=-(kConeOutRadius+0.03)+0.0385, x=0;
 950
 951   TGeoVolume *mechLaddVol = 0;
 952
 953   if (fBuildLevel < 5) {
 954     // world (trapezoid)
 955     TGeoXtru *mechStruct = new TGeoXtru(2); //z sections
 956     Double_t xv[5] = {kStaveWidth/2+0.1,kStaveWidth/2+0.1,0,-kStaveWidth/2-0.1,-kStaveWidth/2-0.1};
 957     Double_t yv[5] = {-kConeOutRadius*2-0.07,0,kStaveHeigth,0,-kConeOutRadius*2-0.07};
 958     mechStruct->DefinePolygon(5,xv,yv);
 959     mechStruct->DefineSection(0,-kStaveLength-0.1,0,0,1.);
 960     mechStruct->DefineSection(1,kStaveLength+0.1,0,0,1.);
 961
 962     mechLaddVol = new TGeoVolume(volname, mechStruct, medAir);
 963     mechLaddVol->SetLineColor(12);
 964     mechLaddVol->SetFillColor(12);
 965     mechLaddVol->SetVisibility(kTRUE);
 966
 967     //Pipe Kapton grey-35
 968     TGeoCone *cone1 = new TGeoCone(kStaveLength,kConeInRadius,kConeOutRadius,kConeInRadius,kConeOutRadius);
 969     TGeoVolume *volCone1= new TGeoVolume("PolyimidePipe", cone1, medKapton);
 970     volCone1->SetFillColor(35);
 971     volCone1->SetLineColor(35);
 972     mechLaddVol->AddNode(volCone1,1,new TGeoTranslation(x+0.25,y,z));
 973     mechLaddVol->AddNode(volCone1,2,new TGeoTranslation(x-0.25,y,z));
 974   }
 975
 976   if (fBuildLevel < 4) {
 977
 978     TGeoTube *coolTubeW = new TGeoTube(0.,kConeInRadius,kStaveLength);
 979     TGeoVolume *volCoolTubeW= new TGeoVolume("Water", coolTubeW, medWater);
 980     volCoolTubeW->SetFillColor(4);
 981     volCoolTubeW->SetLineColor(4);
 982     mechLaddVol->AddNode(volCoolTubeW,0,new TGeoTranslation(x-0.25,y,z));
 983     mechLaddVol->AddNode(volCoolTubeW,1,new TGeoTranslation(x+0.25,y,z));
 984   }
 985
 986   if (fBuildLevel < 3) {
 987     //top fillament
 988     // Top filament M60J black-12 Carbon structure TGeoBBox (length,thickness,width)
 989     TGeoBBox *t2=new TGeoBBox(kS2,0.02/2,0.04/2); //TGeoBBox *t2=new TGeoBBox(kS2,0.01,0.02);
 990     TGeoVolume *volT2=new TGeoVolume("TopFilament", t2, medM60J3K);
 991     volT2->SetLineColor(12);
 992     volT2->SetFillColor(12);
 993     for(int i=0;i<loop;i++){// i<28;i++){
 994       mechLaddVol->AddNode(volT2,i*4+1,new TGeoCombiTrans(x+kWidth,y+kHeight+(0.12/2)-0.014+0.007,z-kStaveLength+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,90-kAlpha,90-kBeta)));
 995       mechLaddVol->AddNode(volT2,i*4+2,new TGeoCombiTrans(x-kWidth,y+kHeight+(0.12/2)-0.014+0.007,z-kStaveLength+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,-90+kAlpha,-90+kBeta)));
 996       mechLaddVol->AddNode(volT2,i*4+3,new TGeoCombiTrans(x+kWidth,y+kHeight+(0.12/2)-0.014+0.007,z-kStaveLength+2*kL1+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,-90+kAlpha,90-kBeta)));
 997       mechLaddVol->AddNode(volT2,i*4+4,new TGeoCombiTrans(x-kWidth,y+kHeight+(0.12/2)-0.014+0.007,z-kStaveLength+2*kL1+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,90-kAlpha,-90+kBeta)));
 998 //    mechLaddVol->AddNode(volT2,i*4+1,new TGeoCombiTrans(x+kWidth+0.0036,y+kHeight-(0.12/2)+0.072,z+kStaveLength+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,90-kAlpha,90-kBeta)));
 999
1000  }
1001
1002     //wall side structure out
1003     TGeoBBox *box4 = new TGeoBBox(0.03/2,0.12/2,kStaveLength-0.50);
1004     TGeoVolume *plate4 = new TGeoVolume("WallOut",box4,medM60J3K);
1005     plate4->SetFillColor(35);
1006     plate4->SetLineColor(35);
1007     mechLaddVol->AddNode(plate4,1,new TGeoCombiTrans(x+(2*kStaveWidth/4)-(0.03/2),y-0.0022-kConeOutRadius+0.12/2+0.007,z,new TGeoRotation("plate4",0,0,0)));
1008     mechLaddVol->AddNode(plate4,2,new TGeoCombiTrans(x-(2*kStaveWidth/4)+(0.03/2),y-0.0022-kConeOutRadius+0.12/2+0.007,z,new TGeoRotation("plate4",0,0,0)));
1009     //wall side in
1010     TGeoBBox *box5 = new TGeoBBox(0.015/2,0.12/2,kStaveLength-0.50);
1011     TGeoVolume *plate5 = new TGeoVolume("WallIn",box5,medM60J3K);
1012     plate5->SetFillColor(12);
1013     plate5->SetLineColor(12);
1014     mechLaddVol->AddNode(plate5,1,new TGeoCombiTrans(x+(2*kStaveWidth/4)-0.03-0.015/2,y-0.0022-kConeOutRadius+0.12/2+0.007,z,new TGeoRotation("plate5",0,0,0)));
1015     mechLaddVol->AddNode(plate5,2,new TGeoCombiTrans(x-(2*kStaveWidth/4)+0.03+0.015/2,y-0.0022-kConeOutRadius+0.12/2+0.007,z,new TGeoRotation("plate5",0,0,0)));
1016
1017      //Amec Thermasol red-2 cover tube FGS300
1018     TGeoConeSeg *cons1 = new TGeoConeSeg(kStaveLength-0.50,kConeOutRadius,kConeOutRadius+kLay1,kConeOutRadius,kConeOutRadius+kLay1,0,180);
1019     TGeoVolume *cone11 = new TGeoVolume("ThermasolPipeCover",cons1,medFGS003);
1020     cone11->SetFillColor(2);
1021     cone11->SetLineColor(2);
1022     mechLaddVol->AddNode(cone11,1,new TGeoCombiTrans(x+0.25,y,z,new TGeoRotation("Cone11",0,0,0)));
1023     mechLaddVol->AddNode(cone11,2,new TGeoCombiTrans(x-0.25,y,z,new TGeoRotation("Cone11",0,0,0)));
1024
1025     TGeoBBox *box2 = new TGeoBBox((0.50-(2*kConeOutRadius))/2,kLay1/2,kStaveLength-0.50);
1026     TGeoVolume *plate2 = new TGeoVolume("ThermasolMiddle",box2,medFGS003);
1027     plate2->SetFillColor(2);
1028     plate2->SetLineColor(2);
1029     mechLaddVol->AddNode(plate2,1,new TGeoCombiTrans(x,y-kConeOutRadius+(kLay1/2),z,new TGeoRotation("plate2",0,0,0)));
1030
1031     TGeoBBox *box21 = new TGeoBBox((0.75-0.25-kConeOutRadius-kLay1)/2,kLay1/2,kStaveLength-0.50);
1032     TGeoVolume *plate21 = new TGeoVolume("ThermasolLeftRight",box21,medFGS003);
1033     plate21->SetFillColor(2);
1034     plate21->SetLineColor(2);
1035     mechLaddVol->AddNode(plate21,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+(0.75-0.25-kConeOutRadius)/2-(kLay1/2),y-kConeOutRadius+(kLay1/2),z,new TGeoRotation("plate21",0,0,0)));
1036     mechLaddVol->AddNode(plate21,2,new TGeoCombiTrans(x-0.25-kConeOutRadius-(0.75-0.25-kConeOutRadius)/2+(kLay1/2),y-kConeOutRadius+(kLay1/2),z,new TGeoRotation("plate21",0,0,0)));
1037
1038     TGeoBBox *box22 = new TGeoBBox((kLay1/2),kConeOutRadius/2,kStaveLength-0.50);
1039     TGeoVolume *plate22 = new TGeoVolume("ThermasolVertical",box22,medFGS003);
1040     plate22->SetFillColor(2);
1041     plate22->SetLineColor(2);
1042     mechLaddVol->AddNode(plate22,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+(kLay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1043     mechLaddVol->AddNode(plate22,2,new TGeoCombiTrans(x+0.25-kConeOutRadius-(kLay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1044     mechLaddVol->AddNode(plate22,3,new TGeoCombiTrans(x-0.25+kConeOutRadius+(kLay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1045     mechLaddVol->AddNode(plate22,4,new TGeoCombiTrans(x-0.25-kConeOutRadius-(kLay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1046
1047     //C Fleece
1048     TGeoConeSeg *cons2 = new TGeoConeSeg(kStaveLength-0.50,kConeOutRadius+kLay1,kConeOutRadius+kLay1+kLay2,kConeOutRadius+kLay1,kConeOutRadius+kLay1+kLay2,0,180);
1049     TGeoVolume *cone12 = new TGeoVolume("CFleecePipeCover",cons2,medCarbonFleece);
1050     cone12->SetFillColor(28);
1051     cone12->SetLineColor(28);
1052     mechLaddVol->AddNode(cone12,1,new TGeoCombiTrans(x+0.25,y,z,new TGeoRotation("Cone12",0,0,0)));
1053     mechLaddVol->AddNode(cone12,2,new TGeoCombiTrans(x-0.25,y,z,new TGeoRotation("Cone12",0,0,0)));
1054
1055     TGeoBBox *box3 = new TGeoBBox((0.50-(2*(kConeOutRadius+kLay1)))/2,kLay2/2,kStaveLength-0.50);
1056     TGeoVolume *plate3 = new TGeoVolume("CFleeceMiddle",box3,medCarbonFleece);
1057     plate3->SetFillColor(28);
1058     plate3->SetLineColor(28);
1059     mechLaddVol->AddNode(plate3,1,new TGeoCombiTrans(x,y-kConeOutRadius+kLay1+(kLay2/2),z,new TGeoRotation("plate3",0,0,0)));
1060
1061     TGeoBBox *box31 = new TGeoBBox((0.75-0.25-kConeOutRadius-kLay1)/2,kLay2/2,kStaveLength-0.50);
1062     TGeoVolume *plate31 = new TGeoVolume("CFleeceLeftRight",box31,medCarbonFleece);
1063     plate31->SetFillColor(28);
1064     plate31->SetLineColor(28);
1065     mechLaddVol->AddNode(plate31,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+kLay1+(0.75-0.25-kConeOutRadius-kLay1)/2,y-kConeOutRadius+kLay1+(kLay2/2),z,new TGeoRotation("plate31",0,0,0)));
1066     mechLaddVol->AddNode(plate31,2,new TGeoCombiTrans(x-0.25-kConeOutRadius-kLay1-(0.75-0.25-kConeOutRadius-kLay1)/2,y-kConeOutRadius+kLay1+(kLay2/2),z,new TGeoRotation("plate31",0,0,0)));
1067
1068     TGeoBBox *box32 = new TGeoBBox((kLay2/2),(kConeOutRadius-kLay1)/2,kStaveLength-0.50);
1069     TGeoVolume *plate32 = new TGeoVolume("CFleeceVertical",box32,medCarbonFleece);
1070     plate32->SetFillColor(28);
1071     plate32->SetLineColor(28);
1072     mechLaddVol->AddNode(plate32,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+kLay1+(kLay2/2),y+(kLay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1073     mechLaddVol->AddNode(plate32,2,new TGeoCombiTrans(x+0.25-kConeOutRadius-kLay1-(kLay2/2),y+(kLay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1074     mechLaddVol->AddNode(plate32,3,new TGeoCombiTrans(x-0.25+kConeOutRadius+kLay1+(kLay2/2),y+(kLay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1075     mechLaddVol->AddNode(plate32,4,new TGeoCombiTrans(x-0.25-kConeOutRadius-kLay1-(kLay2/2),y+(kLay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1076
1077
1078     //K13D2U carbon plate
1079     TGeoBBox *box1 = new TGeoBBox(2*kWidth,kLay3/2,kStaveLength-0.50);
1080     TGeoVolume *plate1 = new TGeoVolume("CarbonPlate",box1,medK13D2U2k);
1081     plate1->SetFillColor(5);
1082     plate1->SetLineColor(5);
1083     mechLaddVol->AddNode(plate1,1,new TGeoCombiTrans(x,y-(kConeOutRadius+(kLay3/2)),z,new TGeoRotation("plate1",0,0,0)));
1084
1085     //C Fleece bottom plate
1086     TGeoBBox *box6 = new TGeoBBox(2*kWidth,kLay2/2,kStaveLength-0.50);
1087     TGeoVolume *plate6 = new TGeoVolume("CFleeceBottom",box6,medCarbonFleece);
1088     plate6->SetFillColor(2);
1089     plate6->SetLineColor(2);
1090     mechLaddVol->AddNode(plate6,1,new TGeoCombiTrans(x,y-(kConeOutRadius+kLay3+(kLay2/2)),z,new TGeoRotation("plate1",0,0,0)));
1091
1092
1093   }
1094
1095   if (fBuildLevel < 2) {
1096     //Glue layers and kapton
1097     TGeoBBox *glue = new TGeoBBox(kStaveWidth/2, 0.005/2, zlad);
1098     TGeoVolume *volGlue=new TGeoVolume("Glue", glue, medGlue);
1099     volGlue->SetLineColor(5);
1100     volGlue->SetFillColor(5);
1101     mechLaddVol->AddNode(volGlue, 0, new TGeoCombiTrans(x,y-(kConeOutRadius+kLay3+(kLay2/2)+(0.01/2)), z, new TGeoRotation("",0, 0, 0)));
1102     mechLaddVol->AddNode(volGlue, 1, new TGeoCombiTrans(x,y-(kConeOutRadius+kLay3+(kLay2/2)+0.01+fSensorThick+(0.01/2)), z, new TGeoRotation("",0, 0, 0)));
1103   }
1104
1105   if (fBuildLevel < 1) {
1106     TGeoBBox *kapCable = new TGeoBBox(kStaveWidth/2, 0.01/2, zlad);
1107     TGeoVolume *volCable=new TGeoVolume("FlexCable", kapCable, medFlexCable);
1108     volCable->SetLineColor(28);
1109     volCable->SetFillColor(28);
1110     mechLaddVol->AddNode(volCable, 0, new TGeoCombiTrans(x, y-(kConeOutRadius+kLay3+(kLay2/2)+0.01+fSensorThick+0.01+(0.01/2)), z, new TGeoRotation("",0, 0, 0)));
1111   }
1112
1113
1114   // Done, return the stave structure
1115   return mechLaddVol;
1116
1117 }
1118 // new model22
1119 //________________________________________________________________________
1120 TGeoVolume* AliITSUv0Layer::CreateStaveModel22(const Double_t xlad,
1121                                                const Double_t zlad,
1122                                                const TGeoManager *mgr){
1123 //
1124 // Create the mechanical stave structure for Model 2.2 of TDR
1125 //
1126 // Input:
1127 //         xlad : X length
1128 //         zlad : Z length
1129 //         mgr  : the GeoManager (used only to get the proper material)
1130 //
1131 // Output:
1132 //
1133 // Return:
1134 //
1135 // Created:      22 Mar 2013  Chinorat Kobdaj
1136 // Updated:      26 Apr 2013  Mario Sitta
1137 // Updated:      30 Apr 2013  Wanchaloem Poonsawat
1138 //
1139
1140   // Materials defined in AliITSUv0
1141   TGeoMedium *medAir    = mgr->GetMedium("ITS_AIR$");
1142   TGeoMedium *medWater  = mgr->GetMedium("ITS_WATER$");
1143
1144   TGeoMedium *medM60J3K    = mgr->GetMedium("ITS_M60J3K$");
1145   TGeoMedium *medKapton    = mgr->GetMedium("ITS_KAPTON(POLYCH2)$");
1146   TGeoMedium *medGlue      = mgr->GetMedium("ITS_GLUE$");
1147   TGeoMedium *medFlexCable = mgr->GetMedium("ITS_FLEXCABLE$");
1148   TGeoMedium *medK13D2U2k  = mgr->GetMedium("ITS_K13D2U2k$");
1149   TGeoMedium *medFGS003    = mgr->GetMedium("ITS_FGS003$");
1150   TGeoMedium *medCarbonFleece = mgr->GetMedium("ITS_CarbonFleece$");
1151
1152   // Local parameters
1153   Double_t kConeOutRadius =0.107/2;//0.107/2;
1154   Double_t kConeInRadius = 0.1015/2;//0.10105/2
1155   Double_t kStaveLength = zlad;
1156   Double_t kStaveWidth = xlad*2;
1157   Double_t kWidth = (kStaveWidth)/4;
1158   Double_t kStaveHeight = 0.283;//0.33;
1159   Double_t kHeight = (kStaveHeight)/2;
1160   Double_t kAlpha = 57;//56.31;
1161   Double_t kTheta = kAlpha*TMath::DegToRad();
1162   Double_t kS1 = ((kStaveWidth)/4)/TMath::Sin(kTheta);
1163   Double_t kL1 = (kStaveWidth/4)/TMath::Tan(kTheta);
1164   Double_t kS2 = sqrt(kHeight*kHeight + kS1*kS1);//TMath::Sin(kThe2);
1165   Double_t kThe2 = TMath::ATan(kHeight/(0.375-0.036));
1166   Double_t kBeta = kThe2*TMath::RadToDeg();
1167   Double_t klay1 = 0.003;//Amec carbon
1168   Double_t klay2 = 0.002;//C Fleece carbon
1169   Double_t klay3 = 0.007;//CFplate K13D2U carbon
1170   Double_t klay4 = 0.007;//GluekStaveLength/2
1171   Double_t klay5 = 0.01;//Flex cable
1172   Double_t kTopVertexMaxWidth = 0.072;
1173   Double_t kTopVertexHeight = 0.04;
1174   Double_t kSideVertexMWidth = 0.052;
1175   Double_t kSideVertexHeight = 0.11;
1176
1177
1178   Int_t  loop = (Int_t)(kStaveLength/(2*kL1));
1179
1180   char volname[30];
1181   snprintf(volname, 30, "%s%d_StaveStruct", AliITSUGeomTGeo::GetITSLadderPattern(), fLayerNumber);
1182
1183   Double_t z=0, y=-(2*kConeOutRadius)+klay1+klay2+fSensorThick/2-0.0004, x=0;
1184
1185   TGeoVolume *mechLaddVol = 0;
1186
1187   if (fBuildLevel < 5) {
1188     // world (trapezoid)
1189     TGeoXtru *mechStruct = new TGeoXtru(2); //z sections
1190     Double_t xv[6] = {kStaveWidth/2,kStaveWidth/2,0.012,-0.012,-kStaveWidth/2,-kStaveWidth/2};
1191     /* Double_t yv[6] = {-2*(kConeOutRadius+klay1+1.5*klay2+klay3+klay4+fSensorThick+klay5),
1192                         0-0.02,kStaveHeight+0.01,kStaveHeight+0.01,0-0.02,
1193                         -2*(kConeOutRadius+klay1+1.5*klay2+klay3+klay4+fSensorThick+klay5)};  // (kConeOutRadius*2)-0.0635 */
1194     Double_t yv[6] = {-(kConeOutRadius*2)-0.06395,0-0.02,kStaveHeight+0.01,kStaveHeight+0.01,0-0.02,-(kConeOutRadius*2)-0.06395};  // (kConeOutRadius*2)-0.064
1195     mechStruct->DefinePolygon(6,xv,yv);
1196     mechStruct->DefineSection(0,-kStaveLength,0,0,1.);
1197     mechStruct->DefineSection(1,kStaveLength,0,0,1.);
1198
1199     mechLaddVol = new TGeoVolume(volname, mechStruct, medAir);
1200     mechLaddVol->SetLineColor(12);
1201     mechLaddVol->SetFillColor(12);
1202     mechLaddVol->SetVisibility(kTRUE);
1203
1204     //Polyimide Pipe Kapton grey-35
1205     TGeoCone *cone1 = new TGeoCone(kStaveLength,kConeInRadius,kConeOutRadius-0.0001,kConeInRadius,kConeOutRadius-0.0001);
1206     TGeoVolume *volCone1= new TGeoVolume("PolyimidePipe", cone1, medKapton);
1207     volCone1->SetFillColor(35);
1208     volCone1->SetLineColor(35);
1209     mechLaddVol->AddNode(volCone1,1,new TGeoTranslation(x+0.25,y,z));
1210     mechLaddVol->AddNode(volCone1,2,new TGeoTranslation(x-0.25,y,z));
1211     }
1212
1213   if (fBuildLevel < 4) {
1214     TGeoTube *coolTubeW = new TGeoTube(0.,kConeInRadius-0.0001,kStaveLength);
1215     TGeoVolume *volCoolTubeW= new TGeoVolume("Water", coolTubeW, medWater);
1216     volCoolTubeW->SetFillColor(4);
1217     volCoolTubeW->SetLineColor(4);
1218     mechLaddVol->AddNode(volCoolTubeW,0,new TGeoTranslation(x-0.25,y,z));
1219     mechLaddVol->AddNode(volCoolTubeW,1,new TGeoTranslation(x+0.25,y,z));
1220   }
1221
1222   if (fBuildLevel < 3) {
1223     //top fillament
1224     // Top filament M60J black-12 Carbon structure TGeoBBox (length,thickness,width)
1225     TGeoBBox *t2=new TGeoBBox(kS2-0.028,0.02/2,0.02/2); //0.04/2//TGeoBBox *t2=new TGeoBBox(kS2,0.01,0.02);//kS2-0.03 old Config.C
1226     TGeoVolume *volT2=new TGeoVolume("TopFilament", t2, medM60J3K);
1227     volT2->SetLineColor(12);
1228     volT2->SetFillColor(12);
1229     for(int i=0;i<loop;i++){// i<28;i++){
1230       // 1) Front Left Top Filament
1231        mechLaddVol->AddNode(volT2,i*4+1,new TGeoCombiTrans(x+kWidth+0.0036,y+kHeight+0.01,z-kStaveLength+0.1+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,90-kAlpha,90-kBeta)));
1232       // 2) Front Right Top Filament
1233       mechLaddVol->AddNode(volT2,i*4+2,new TGeoCombiTrans(x-kWidth-0.0036,y+kHeight+0.01,z-kStaveLength+0.1+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,-90+kAlpha,-90+kBeta)));
1234       // 3) Back Left  Top Filament
1235       mechLaddVol->AddNode(volT2,i*4+3,new TGeoCombiTrans(x+kWidth+0.0036,y+kHeight+0.01,z-kStaveLength+0.1+2*kL1+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,-90+kAlpha,90-kBeta)));
1236       // 4) Back Right Top Filament
1237       mechLaddVol->AddNode(volT2,i*4+4,new TGeoCombiTrans(x-kWidth-0.0036,y+kHeight+0.01,z-kStaveLength+0.1+2*kL1+(i*4*kL1)+kS1/2, new TGeoRotation("volT2",90,90-kAlpha,-90+kBeta)));
1238    }
1239
1240      //Vertex  structure
1241
1242       //top ver trd1
1243       TGeoTrd1 *trd1 = new TGeoTrd1(0,kTopVertexMaxWidth/2,kStaveLength,kTopVertexHeight/2);
1244       TGeoVolume *ibdv = new TGeoVolume("TopVertex",trd1,medM60J3K);
1245       ibdv->SetFillColor(12);
1246       ibdv->SetLineColor(12);
1247       mechLaddVol->AddNode(ibdv,1,new TGeoCombiTrans(x,y+kStaveHeight+0.03,z,new TGeoRotation("ibdv",0.,-90,0)));//y+kStaveHeight+0.056
1248
1249       //left trd2
1250       TGeoTrd1 *trd2 = new TGeoTrd1(0,kSideVertexMWidth/2,kStaveLength, kSideVertexHeight/2);
1251       TGeoVolume *ibdv2 = new TGeoVolume("LeftVertex",trd2,medM60J3K);
1252       ibdv2->SetFillColor(12);
1253       ibdv2->SetLineColor(12);
1254       mechLaddVol->AddNode(ibdv2,1,new TGeoCombiTrans(x+kStaveWidth/2-0.06,y-0.0355,z,new TGeoRotation("ibdv2",-103.3,90,0))); //x-kStaveWidth/2-0.09 old Config.C y-0.0355,
1255
1256       //right trd3
1257       TGeoTrd1 *trd3 = new TGeoTrd1(0,kSideVertexMWidth/2,kStaveLength, kSideVertexHeight/2);
1258       TGeoVolume *ibdv3 = new TGeoVolume("RightVertex",trd3,medM60J3K);
1259       ibdv3->SetFillColor(12);
1260       ibdv3->SetLineColor(12);
1261       mechLaddVol->AddNode(ibdv3,1,new TGeoCombiTrans(x-kStaveWidth/2+0.06,y-0.0355,z,new TGeoRotation("ibdv3",103.3,90,0))); //x-kStaveWidth/2+0.09 old Config.C
1262
1263      //Carbon Fleece
1264       TGeoConeSeg *cons2 = new TGeoConeSeg(zlad,kConeOutRadius+klay1,kConeOutRadius+klay1+klay2,kConeOutRadius+klay1,kConeOutRadius+klay1+klay2,0,180);
1265       TGeoVolume *cone12 = new TGeoVolume("CarbonFleecePipeCover",cons2,medCarbonFleece);
1266       cone12->SetFillColor(28);
1267       cone12->SetLineColor(28);
1268       mechLaddVol->AddNode(cone12,1,new TGeoCombiTrans(x+0.25,y,z,new TGeoRotation("cone12",0,0,0)));
1269       mechLaddVol->AddNode(cone12,2,new TGeoCombiTrans(x-0.25,y,z,new TGeoRotation("cone12",0,0,0)));
1270
1271       TGeoBBox *box3 = new TGeoBBox((0.50-(2*(kConeOutRadius+klay1)))/2,klay2/2,zlad);//kStaveLength-0.50);
1272       TGeoVolume *plate3 = new TGeoVolume("CarbonFleeceMiddle",box3,medCarbonFleece);
1273       plate3->SetFillColor(28);
1274       plate3->SetLineColor(28);
1275       mechLaddVol->AddNode(plate3,1,new TGeoCombiTrans(x,y-kConeOutRadius+klay1+(klay2/2),z,new TGeoRotation("plate3",0,0,0)));
1276
1277       TGeoBBox *box31 = new TGeoBBox((0.75-0.25-kConeOutRadius-klay1)/2+0.0025,klay2/2,zlad);
1278       TGeoVolume *plate31 = new TGeoVolume("CarbonFleeceLeftRight",box31,medCarbonFleece);
1279       plate31->SetFillColor(28);
1280       plate31->SetLineColor(28);
1281       mechLaddVol->AddNode(plate31,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+klay1+(0.75-0.25-kConeOutRadius-klay1)/2,y-kConeOutRadius+klay1+(klay2/2),z,new TGeoRotation("plate31",0,0,0)));
1282       mechLaddVol->AddNode(plate31,2,new TGeoCombiTrans(x-0.25-kConeOutRadius-klay1-(0.75-0.25-kConeOutRadius-klay1)/2,y-kConeOutRadius+klay1+(klay2/2),z,new TGeoRotation("plate31",0,0,0)));
1283
1284       TGeoBBox *box32 = new TGeoBBox((klay2/2),(kConeOutRadius-klay1)/2,zlad);
1285       TGeoVolume *plate32 = new TGeoVolume("CarbonFleeceVertical",box32,medCarbonFleece);
1286       plate32->SetFillColor(28);
1287       plate32->SetLineColor(28);
1288       mechLaddVol->AddNode(plate32,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+klay1+(klay2/2),y+(klay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1289       mechLaddVol->AddNode(plate32,2,new TGeoCombiTrans(x+0.25-kConeOutRadius-klay1-(klay2/2),y+(klay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1290       mechLaddVol->AddNode(plate32,3,new TGeoCombiTrans(x-0.25+kConeOutRadius+klay1+(klay2/2),y+(klay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1291       mechLaddVol->AddNode(plate32,4,new TGeoCombiTrans(x-0.25-kConeOutRadius-klay1-(klay2/2),y+(klay1-kConeOutRadius)/2,z,new TGeoRotation("plate32",0,0,0)));
1292
1293      //Amec Thermasol red-2 cover tube FGS300 or Carbon Paper
1294       TGeoConeSeg *cons1 = new TGeoConeSeg(zlad,kConeOutRadius,kConeOutRadius+klay1-0.0001,kConeOutRadius,kConeOutRadius+klay1-0.0001,0,180);//kConeOutRadius+klay1-0.0001
1295       TGeoVolume *cone11 = new TGeoVolume("ThermasolPipeCover",cons1,medFGS003);
1296       cone11->SetFillColor(2);
1297       cone11->SetLineColor(2);
1298       mechLaddVol->AddNode(cone11,1,new TGeoCombiTrans(x+0.25,y,z,new TGeoRotation("cone11",0,0,0)));
1299       mechLaddVol->AddNode(cone11,2,new TGeoCombiTrans(x-0.25,y,z,new TGeoRotation("cone11",0,0,0)));
1300
1301       TGeoBBox *box2 = new TGeoBBox((0.50-(2*kConeOutRadius))/2,(klay1/2),zlad);//kStaveLength-0.50);
1302       TGeoVolume *plate2 = new TGeoVolume("ThermasolMiddle",box2,medFGS003);
1303       plate2->SetFillColor(2);
1304       plate2->SetLineColor(2);
1305       mechLaddVol->AddNode(plate2,1,new TGeoCombiTrans(x,y-kConeOutRadius+(klay1/2),z,new TGeoRotation("plate2",0,0,0)));
1306
1307       TGeoBBox *box21 = new TGeoBBox((0.75-0.25-kConeOutRadius-klay1)/2+0.0025,(klay1/2),zlad);
1308       TGeoVolume *plate21 = new TGeoVolume("ThermasolLeftRight",box21,medFGS003);
1309       plate21->SetFillColor(2);
1310       plate21->SetLineColor(2);
1311       mechLaddVol->AddNode(plate21,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+(0.75-0.25-kConeOutRadius)/2-(klay1/2)+0.0025,y-kConeOutRadius+(klay1/2),z,new TGeoRotation("plate21",0,0,0)));
1312       mechLaddVol->AddNode(plate21,2,new TGeoCombiTrans(x-0.25-kConeOutRadius-(0.75-0.25-kConeOutRadius)/2+(klay1/2)-0.0025,y-kConeOutRadius+(klay1/2),z,new TGeoRotation("plate21",0,0,0)));
1313
1314       TGeoBBox *box22 = new TGeoBBox((klay1/2),kConeOutRadius/2,zlad);
1315       TGeoVolume *plate22 = new TGeoVolume("ThermasolVertical",box22,medFGS003);
1316       plate22->SetFillColor(2);
1317       plate22->SetLineColor(2);
1318       mechLaddVol->AddNode(plate22,1,new TGeoCombiTrans(x+0.25+kConeOutRadius+(klay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1319       mechLaddVol->AddNode(plate22,2,new TGeoCombiTrans(x+0.25-kConeOutRadius-(klay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1320       mechLaddVol->AddNode(plate22,3,new TGeoCombiTrans(x-0.25+kConeOutRadius+(klay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1321       mechLaddVol->AddNode(plate22,4,new TGeoCombiTrans(x-0.25-kConeOutRadius-(klay1/2),y-kConeOutRadius/2,z,new TGeoRotation("plate22",0,0,0)));
1322
1323      //K13D2U CF plate
1324       TGeoBBox *box1 = new TGeoBBox(2*kWidth,(klay3)/2,zlad);
1325       TGeoVolume *plate1 = new TGeoVolume("CFPlate",box1,medK13D2U2k);
1326       plate1->SetFillColor(5);
1327       plate1->SetLineColor(5);
1328       mechLaddVol->AddNode(plate1,1,new TGeoCombiTrans(x,y-(kConeOutRadius+(klay3/2)),z,new TGeoRotation("plate1",0,0,0)));
1329
1330      //C Fleece bottom plate
1331       TGeoBBox *box6 = new TGeoBBox(2*kWidth,(klay2)/2,zlad);
1332       TGeoVolume *plate6 = new TGeoVolume("CarbonFleeceBottom",box6,medCarbonFleece);
1333       plate6->SetFillColor(2);
1334       plate6->SetLineColor(2);
1335       mechLaddVol->AddNode(plate6,1,new TGeoCombiTrans(x,y-(kConeOutRadius+klay3+(klay2/2)),z,new TGeoRotation("plate6",0,0,0)));
1336
1337     }
1338    if (fBuildLevel < 2) {
1339       //Glue klayers and kapton
1340      TGeoBBox *glue = new TGeoBBox(kStaveWidth/2, (klay4)/2, zlad);
1341       TGeoVolume *volGlue=new TGeoVolume("Glue", glue, medGlue);
1342       volGlue->SetLineColor(5);
1343       volGlue->SetFillColor(5);
1344       // mechLaddVol->AddNode(volGlue, 0, new TGeoCombiTrans(x,y-(kConeOutRadius+klay3+klay2+(klay4/2)), z, new TGeoRotation("",0, 0, 0)));
1345       mechLaddVol->AddNode(volGlue, 0, new TGeoCombiTrans(x,y-(kConeOutRadius+klay3+klay2+(klay4)/2)+0.00005, z, new TGeoRotation("",0, 0, 0)));
1346     }
1347
1348      if (fBuildLevel < 1) {
1349      //Flex Cable or Bus
1350       TGeoBBox *kapCable = new TGeoBBox(kStaveWidth/2, klay5/2, zlad);//klay5/2
1351       TGeoVolume *volCable=new TGeoVolume("FlexCable", kapCable, medFlexCable);
1352       volCable->SetLineColor(28);
1353       volCable->SetFillColor(28);
1354       //      mechLaddVol->AddNode(volCable, 0, new TGeoCombiTrans(x, y-(kConeOutRadius+klay3+klay2+klay4+fSensorThick+(klay5)/2)+0.0002, z, new TGeoRotation("",0, 0, 0)));
1355       mechLaddVol->AddNode(volCable, 0, new TGeoCombiTrans(x, y-(kConeOutRadius+klay3+klay2+klay4+fSensorThick+(klay5)/2)+0.01185, z, new TGeoRotation("",0, 0, 0)));
1356       }
1357     // Done, return the stave structe
1358     return mechLaddVol;
1359 }
1360
1361 // model3
1362 //________________________________________________________________________
1363 TGeoVolume* AliITSUv0Layer::CreateStaveModel3(const Double_t xlad,
1364                                                const Double_t zlad,
1365                                                const TGeoManager *mgr){
1366 //
1367 // Create the mechanical stave structure for Model 3 of TDR
1368 //
1369 // Input:
1370 //         xlad : X length
1371 //         zlad : Z length
1372 //         mgr  : the GeoManager (used only to get the proper material)
1373 //
1374 // Output:
1375 //
1376 // Return:
1377 //
1378 // Created:      28 May 2013  Chinorat Kobdaj
1379 // Updated:                   Mario Sitta
1380 // Updated:                   Wanchaloem Poonsawat
1381 //
1382
1383   // Materials defined in AliITSUv0
1384   TGeoMedium *medAir    = mgr->GetMedium("ITS_AIR$");
1385   TGeoMedium *medWater  = mgr->GetMedium("ITS_WATER$");
1386
1387   TGeoMedium *medM60J3K    = mgr->GetMedium("ITS_M60J3K$");
1388   TGeoMedium *medKapton    = mgr->GetMedium("ITS_KAPTON(POLYCH2)$");
1389   TGeoMedium *medGlue      = mgr->GetMedium("ITS_GLUE$");
1390   TGeoMedium *medFlexCable = mgr->GetMedium("ITS_FLEXCABLE$");
1391   //TGeoMedium *medK13D2U2k  = mgr->GetMedium("ITS_K13D2U2k$");
1392   //TGeoMedium *medFGS003    = mgr->GetMedium("ITS_FGS003$");
1393   //TGeoMedium *medCarbonFleece = mgr->GetMedium("ITS_CarbonFleece$");
1394
1395   // Local parameters
1396     Double_t kConeOutRadius = 0.15/2;
1397     Double_t kStaveLength = zlad*2;
1398     Double_t kStaveWidth = xlad*2;
1399     Double_t w = kStaveWidth/4;//1/2 of W
1400     Double_t staveHeight = 0.3;
1401     Double_t h = staveHeight/2;
1402     Double_t alpha = 90-33.;//90-30;
1403     Double_t the1 = alpha*TMath::DegToRad();
1404     Double_t s1 = w/TMath::Sin(the1);
1405     Double_t l = w/TMath::Tan(the1);
1406     Double_t s2 = TMath::Sqrt(h*h + s1*s1);//TMath::Sin(the2);
1407     Double_t the2 = TMath::ATan(h/s1);
1408     Double_t beta = the2*TMath::RadToDeg();
1409     Double_t klay4 = 0.007; //Glue
1410     Double_t klay5 = 0.01; //Flexcable
1411     Int_t  loop = (kStaveLength/(2*l))/2;
1412     Double_t hh = 0.01;
1413        Double_t ang1 = 0*TMath::DegToRad();
1414        Double_t ang2 = 0*TMath::DegToRad();
1415        Double_t ang3 = 0*TMath::DegToRad();
1416        Int_t modules = 4;
1417        Double_t headWidth=0.25;
1418        Double_t smcLength=kStaveLength/modules-2*headWidth;//6.25;
1419        Double_t smcWidth=kStaveWidth;
1420        Double_t smcSide1Thick=0.03;
1421        Double_t vaporThick=0.032;
1422        Double_t liquidThick=0.028;
1423        Double_t smcSide2Thick=0.01;
1424        Double_t smcSide3Thick=0.0055;
1425        Double_t smcSide4Thick=0.0095;
1426        Double_t smcSide5Thick=0.0075;
1427        Double_t smcSpace=0.01;
1428
1429
1430     char volname[30];
1431     snprintf(volname, 30, "%s%d_StaveStruct", AliITSUGeomTGeo::GetITSLadderPattern(), fLayerNumber);
1432
1433     // detailed structure ++++++++++++++
1434     Double_t z=0, y=0-0.007, x=0;
1435
1436     // Polimide micro channels numbers
1437     Double_t yMC = y-h+0.01;
1438     Int_t nb = (Int_t)(kStaveWidth/0.1)+1;
1439     Double_t xladMC = (nb*0.1-0.08)/2;
1440
1441
1442     TGeoVolume *mechLaddVol = 0;
1443     if (fBuildLevel < 5) {
1444       // world (trapezoid)
1445       TGeoXtru *mechStruct = new TGeoXtru(2); //z sections
1446       Double_t xv[5] = {kStaveWidth/2+0.1,kStaveWidth/2+0.1,0,-kStaveWidth/2-0.1,-kStaveWidth/2-0.1};
1447       Double_t yv[5] = {-kConeOutRadius*2-0.07,0,staveHeight,0,-kConeOutRadius*2-0.07};
1448       mechStruct->DefinePolygon(5,xv,yv);
1449       mechStruct->DefineSection(0,-kStaveLength-0.1,0,0,1.);
1450       mechStruct->DefineSection(1,kStaveLength+0.1,0,0,1.);
1451       mechLaddVol = new TGeoVolume(volname, mechStruct, medAir);
1452       mechLaddVol->SetLineColor(12);
1453       mechLaddVol->SetFillColor(12);
1454       mechLaddVol->SetVisibility(kTRUE);
1455
1456        // Silicon micro channels numbers
1457
1458       TGeoBBox *tM0a=new TGeoBBox(smcWidth/2, 0.003/2, headWidth/2);
1459       TGeoVolume *volTM0a=new TGeoVolume("microChanTop1", tM0a, medKapton);
1460       volTM0a->SetLineColor(35);
1461       volTM0a->SetFillColor(35);
1462
1463       for(Int_t  mo=1; mo<=modules; mo++) {
1464       mechLaddVol->AddNode(volTM0a, 0, new TGeoCombiTrans(x,yMC+0.03, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth+smcLength/2+(headWidth/2), new TGeoRotation("",ang1, ang2, ang3)));
1465       mechLaddVol->AddNode(volTM0a, 1, new TGeoCombiTrans(x,yMC+0.03, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth-smcLength/2-(headWidth/2), new TGeoRotation("",ang1, ang2, ang3)));
1466       }
1467       TGeoBBox *tM0c=new TGeoBBox(0.3/2, 0.003/2,smcLength/2);
1468       TGeoVolume *volTM0c=new TGeoVolume("microChanTop2", tM0c, medKapton);
1469       volTM0c->SetLineColor(35);
1470       volTM0c->SetFillColor(35);
1471       for(Int_t  mo=1; mo<=modules; mo++) {
1472       mechLaddVol->AddNode(volTM0c, 0, new TGeoCombiTrans(x+(smcWidth/2)-(0.3/2),yMC+0.03, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));
1473       mechLaddVol->AddNode(volTM0c, 1, new TGeoCombiTrans(x-(smcWidth/2)+(0.3/2),yMC+0.03, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1474       }
1475       TGeoBBox *tM0c1=new TGeoBBox(0.2225/2, 0.003/2,smcLength/2);
1476       TGeoVolume *volTM0c1=new TGeoVolume("microChanBot1", tM0c1, medKapton);
1477       volTM0c1->SetLineColor(6);
1478       volTM0c1->SetFillColor(6);
1479       for(Int_t  mo=1; mo<=modules; mo++) {
1480       mechLaddVol->AddNode(volTM0c1, 0, new TGeoCombiTrans(x+smcWidth/2-(smcSide1Thick)-(vaporThick)-(smcSide2Thick)-(smcSide3Thick)-(0.2225/2),yMC+0.03-hh-(0.003), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1481       mechLaddVol->AddNode(volTM0c1, 1, new TGeoCombiTrans(x-smcWidth/2+(smcSide1Thick)+(liquidThick)+(smcSide2Thick)+(smcSide4Thick)+(0.2225/2),yMC+0.03-hh-(0.003), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1482       }
1483       TGeoBBox *tM0c2=new TGeoBBox(0.072/2, 0.003/2,smcLength/2);
1484       TGeoVolume *volTM0c2=new TGeoVolume("microChanBot2", tM0c2, medKapton);
1485       volTM0c2->SetLineColor(35);
1486       volTM0c2->SetFillColor(35);
1487       for(Int_t  mo=1; mo<=modules; mo++) {
1488       mechLaddVol->AddNode(volTM0c2, 0, new TGeoCombiTrans(x+smcWidth/2-(0.072/2),yMC+0.03-(0.035+0.0015)-(0.003)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1489       }
1490       TGeoBBox *tM0c2r=new TGeoBBox(0.068/2, 0.003/2,smcLength/2);
1491       TGeoVolume *volTM0c2r=new TGeoVolume("microChanBot3", tM0c2r, medKapton);
1492       volTM0c2r->SetLineColor(35);
1493       volTM0c2r->SetFillColor(35);
1494       for(Int_t  mo=1; mo<=modules; mo++) {
1495       mechLaddVol->AddNode(volTM0c2r, 0, new TGeoCombiTrans(x-smcWidth/2+(0.068/2),yMC+0.03-(0.035+0.0015)-(0.003)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1496       }
1497       TGeoBBox *tM0d=new TGeoBBox(smcSide1Thick/2, 0.035/2,smcLength/2);
1498       TGeoVolume *volTM0d=new TGeoVolume("microChanSide1", tM0d, medKapton);
1499       volTM0d->SetLineColor(12);
1500       volTM0d->SetFillColor(12);
1501       for(Int_t  mo=1; mo<=modules; mo++) {
1502       mechLaddVol->AddNode(volTM0d, 0, new TGeoCombiTrans(x+smcWidth/2-(smcSide1Thick/2),yMC+0.03-0.0015-(0.035)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1503       mechLaddVol->AddNode(volTM0d, 1, new TGeoCombiTrans(x-smcWidth/2+(smcSide1Thick/2),yMC+0.03-0.0015-(0.035)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1504       }
1505
1506       TGeoBBox *tM0d1=new TGeoBBox(smcSide2Thick/2, 0.035/2,smcLength/2);
1507       TGeoVolume *volTM0d1=new TGeoVolume("microChanSide2", tM0d1, medKapton);
1508       volTM0d1->SetLineColor(12);
1509       volTM0d1->SetFillColor(12);
1510       for(Int_t  mo=1; mo<=modules; mo++) {
1511       mechLaddVol->AddNode(volTM0d1, 0, new TGeoCombiTrans(x+smcWidth/2-(smcSide1Thick)-(vaporThick)-(smcSide2Thick/2),yMC+0.03-(0.003+0.035)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1512       mechLaddVol->AddNode(volTM0d1, 1, new TGeoCombiTrans(x-smcWidth/2+(smcSide1Thick)+(liquidThick)+(smcSide2Thick/2),yMC+0.03-(0.003+0.035)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1513       }
1514       TGeoBBox *tM0d2=new TGeoBBox(smcSide3Thick/2, (hh+0.003)/2, smcLength/2);
1515       TGeoVolume *volTM0d2=new TGeoVolume("microChanSide3", tM0d2, medKapton);
1516       volTM0d2->SetLineColor(12);
1517       volTM0d2->SetFillColor(12);
1518       for(Int_t  mo=1; mo<=modules; mo++) {
1519       mechLaddVol->AddNode(volTM0d2, 0, new TGeoCombiTrans(x+smcWidth/2-(smcSide1Thick)-(vaporThick)-(smcSide2Thick)-(smcSide3Thick/2),yMC+0.03-(0.003+hh+0.003)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1520       }
1521       TGeoBBox *tM0d2r=new TGeoBBox(smcSide4Thick/2, (hh+0.003)/2, smcLength/2);
1522       TGeoVolume *volTM0d2r=new TGeoVolume("microChanSide4", tM0d2r, medKapton);
1523       volTM0d2r->SetLineColor(12);
1524       volTM0d2r->SetFillColor(12);
1525       for(Int_t  mo=1; mo<=modules; mo++) {
1526       mechLaddVol->AddNode(volTM0d2r, 0, new TGeoCombiTrans(x-smcWidth/2+(smcSide1Thick)+(liquidThick)+(smcSide2Thick)+(smcSide4Thick/2),yMC+0.03-(0.003+hh+0.003)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1527       }
1528       TGeoBBox *tM0e=new TGeoBBox(smcSide5Thick/2, hh/2,smcLength/2);
1529       TGeoVolume *volTM0e=new TGeoVolume("microChanSide5", tM0e, medKapton);
1530       volTM0e->SetLineColor(12);
1531       volTM0e->SetFillColor(12);
1532       for(Int_t  mo=1; mo<=modules; mo++) {
1533       for (Int_t ie=0;ie<11;ie++) {
1534         mechLaddVol->AddNode(volTM0e, 0, new TGeoCombiTrans(x-(ie*(smcSpace+smcSide5Thick))+smcWidth/2-(smcSide1Thick)-(vaporThick)-(smcSide2Thick)-(smcSide3Thick)-smcSpace-(smcSide5Thick/2),yMC+0.03-(0.003+hh)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1535         mechLaddVol->AddNode(volTM0e, 1, new TGeoCombiTrans(x+(ie*(smcSpace+smcSide5Thick))-smcWidth/2+(smcSide1Thick)+(liquidThick)+(smcSide2Thick)+(smcSide4Thick)+smcSpace+(smcSide5Thick/2),yMC+0.03-(0.003+hh)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1536          }
1537       }
1538
1539       TGeoBBox *tM0f=new TGeoBBox(0.02/2, hh/2, smcLength/2);
1540       TGeoVolume *volTM0f=new TGeoVolume("microChanTop3", tM0f, medKapton);
1541       //Double_t smcChannels=12;
1542       Double_t smcCloseWallvapor=smcWidth/2-smcSide1Thick-vaporThick-smcSide2Thick-smcSide3Thick-12*smcSpace-11*smcSide5Thick;
1543       Double_t smcCloseWallliquid=smcWidth/2-smcSide1Thick-liquidThick-smcSide2Thick-smcSide4Thick-12*smcSpace-11*smcSide5Thick;
1544       volTM0f->SetLineColor(12);
1545       volTM0f->SetFillColor(12);
1546       for(Int_t  mo=1; mo<=modules; mo++) {
1547        mechLaddVol->AddNode(volTM0f, 0, new TGeoCombiTrans(x+smcCloseWallvapor-(0.02)/2,yMC+0.03-(0.003+hh)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1548        mechLaddVol->AddNode(volTM0f, 1, new TGeoCombiTrans(x-smcCloseWallliquid+(0.02)/2,yMC+0.03-(0.003+hh)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1549       }
1550       //Head(back) microchannel
1551
1552       TGeoBBox *tM0hb=new TGeoBBox(smcWidth/2, 0.025/2, headWidth/2);
1553       TGeoVolume *volTM0hb=new TGeoVolume("microChanHeadBackBottom1", tM0hb, medKapton);
1554       volTM0hb->SetLineColor(4);
1555       volTM0hb->SetFillColor(4);
1556       for(Int_t  mo=1; mo<=modules; mo++) {
1557       mechLaddVol->AddNode(volTM0hb, 0, new TGeoCombiTrans(x,yMC+0.03-0.0145-(0.025/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth+smcLength/2+(headWidth/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1558       mechLaddVol->AddNode(volTM0hb, 1, new TGeoCombiTrans(x,yMC+0.03-0.0145-(0.025)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth-smcLength/2-(headWidth/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1559       }
1560       TGeoBBox *tM0h1=new TGeoBBox(smcWidth/2, 0.013/2, 0.05/2);
1561       TGeoVolume *volTM0h1=new TGeoVolume("microChanHeadBackBottom2", tM0h1, medKapton);
1562       volTM0h1->SetLineColor(5);
1563       volTM0h1->SetFillColor(5);
1564       for(Int_t  mo=1; mo<=modules; mo++) {
1565       mechLaddVol->AddNode(volTM0h1, 0, new TGeoCombiTrans(x,yMC+0.03-0.0015-(0.013/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth-smcLength/2-headWidth+(0.05/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1566       }
1567       TGeoBBox *tM0h2=new TGeoBBox(smcWidth/2, 0.003/2, 0.18/2);
1568       TGeoVolume *volTM0h2=new TGeoVolume("microChanHeadBackBottom7", tM0h2, medKapton);
1569       volTM0h2->SetLineColor(6);
1570       volTM0h2->SetFillColor(6);
1571       for(Int_t  mo=1; mo<=modules; mo++) {
1572       mechLaddVol->AddNode(volTM0h2, 0, new TGeoCombiTrans(x,yMC+0.03-0.0015-0.01-(0.003/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth-smcLength/2-0.02-(0.18/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1573       }
1574       TGeoBBox *tM0h3=new TGeoBBox(smcWidth/2, 0.013/2, 0.02/2);
1575       TGeoVolume *volTM0h3=new TGeoVolume("microChanHeadBackBottom3", tM0h3, medKapton);
1576       volTM0h3->SetLineColor(5);
1577       volTM0h3->SetFillColor(5);
1578       for(Int_t  mo=1; mo<=modules; mo++) {
1579       mechLaddVol->AddNode(volTM0h3, 0, new TGeoCombiTrans(x,yMC+0.03-0.0015-(0.013/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth-smcLength/2-(0.02/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1580       }
1581       TGeoBBox *tM0b1=new TGeoBBox(smcWidth/2, 0.013/2, 0.03/2);
1582       TGeoVolume *volTM0b1=new TGeoVolume("microChanHeadBackBottom4", tM0b1, medKapton);
1583       volTM0b1->SetLineColor(5);
1584       volTM0b1->SetFillColor(5);
1585       for(Int_t  mo=1; mo<=modules; mo++) {
1586       mechLaddVol->AddNode(volTM0b1, 0, new TGeoCombiTrans(x,yMC+0.03-0.0015-(0.013/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth+smcLength/2+headWidth-(0.03/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1587       }
1588       TGeoBBox *tM0b2=new TGeoBBox(smcWidth/2, 0.003/2, 0.2/2);
1589       TGeoVolume *volTM0b2=new TGeoVolume("microChanHeadBackBottom5", tM0b2, medKapton);
1590       volTM0b2->SetLineColor(6);
1591       volTM0b2->SetFillColor(6);
1592       for(Int_t  mo=1; mo<=modules; mo++) {
1593       mechLaddVol->AddNode(volTM0b2, 0, new TGeoCombiTrans(x,yMC+0.03-0.0015-0.01-(0.003/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth+smcLength/2+0.02+(0.2/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1594       }
1595       TGeoBBox *tM0b3=new TGeoBBox(smcWidth/2, 0.013/2, 0.02/2);
1596       TGeoVolume *volTM0b3=new TGeoVolume("microChanHeadBackBottom6", tM0b3, medKapton);
1597       volTM0b3->SetLineColor(5);
1598       volTM0b3->SetFillColor(5);
1599       for(Int_t  mo=1; mo<=modules; mo++) {
1600       mechLaddVol->AddNode(volTM0b3, 0, new TGeoCombiTrans(x,yMC+0.03-0.0015-(0.013/2), z+(mo-3)*kStaveLength/4+smcLength/2+headWidth+smcLength/2+(0.02/2), new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1601       }
1602
1603       TGeoBBox *tM0b=new TGeoBBox(0.02/2, 0.02/2, zlad);
1604       TGeoVolume *volTM0b=new TGeoVolume("microChanWalls", tM0b, medKapton);
1605       volTM0b->SetLineColor(35);
1606       volTM0b->SetFillColor(35);
1607       for (Int_t ib=0;ib<nb;ib++) {
1608         //mechLaddVol->AddNode(volTM0b, ib, new TGeoCombiTrans(x+ib*0.1-xladMC+0.01,yMC, z, new TGeoRotation("",0, 0, 0)));
1609       }
1610
1611       }
1612
1613     if (fBuildLevel < 4) {
1614
1615       //**********cooling  inlet outlet
1616
1617       TGeoBBox *tM0dv=new TGeoBBox(vaporThick/2, 0.035/2,smcLength/2);
1618       TGeoVolume *volTM0dv=new TGeoVolume("microChanVapor", tM0dv, medWater);
1619       volTM0dv->SetLineColor(2);
1620       volTM0dv->SetFillColor(2);
1621       for(Int_t  mo=1; mo<=modules; mo++) {
1622       mechLaddVol->AddNode(volTM0dv, 0, new TGeoCombiTrans(x+smcWidth/2-(smcSide1Thick)-(vaporThick/2),yMC+0.03-0.0015-(0.035)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1623       }
1624       TGeoBBox *tM0dl=new TGeoBBox(liquidThick/2, 0.035/2,smcLength/2);
1625       TGeoVolume *volTM0dl=new TGeoVolume("microChanLiquid", tM0dl, medWater);
1626       volTM0dl->SetLineColor(3);
1627       volTM0dl->SetFillColor(3);
1628       for(Int_t  mo=1; mo<=modules; mo++) {
1629       mechLaddVol->AddNode(volTM0dl, 0, new TGeoCombiTrans(x-smcWidth/2+(smcSide1Thick)+(liquidThick/2),yMC+0.03-0.0015-(0.035)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1630       }
1631       // small cooling fluid now using water wait for freeon value
1632       TGeoBBox *tM0dlq=new TGeoBBox(smcSpace/2, hh/2,smcLength/2);
1633       TGeoVolume *volTM0dlq=new TGeoVolume("smallLiquid", tM0dlq, medWater);
1634       volTM0dlq->SetLineColor(3);
1635       volTM0dlq->SetFillColor(3);
1636       TGeoBBox *tM0dvp=new TGeoBBox(smcSpace/2, hh/2,smcLength/2);
1637       TGeoVolume *volTM0dvp=new TGeoVolume("microChanVapor", tM0dvp, medWater);
1638       volTM0dvp->SetLineColor(2);
1639       volTM0dvp->SetFillColor(2);
1640       for(Int_t  mo=1; mo<=modules; mo++) {
1641       for (Int_t is=0;is<12;is++) {
1642         mechLaddVol->AddNode(volTM0dlq, 0, new TGeoCombiTrans(x+(is*(smcSpace+smcSide5Thick))-smcWidth/2+(smcSide1Thick)+(vaporThick)+(smcSide2Thick)+(smcSide3Thick)+smcSpace/2,yMC+0.03-(0.003+hh)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1643         mechLaddVol->AddNode(volTM0dvp, 1, new TGeoCombiTrans(x-(is*(smcSpace+smcSide5Thick))+smcWidth/2-(smcSide1Thick)-(vaporThick)-(smcSide2Thick)-(smcSide3Thick)-smcSpace/2,yMC+0.03-(0.003+hh)/2, z+(mo-3)*kStaveLength/4+smcLength/2+headWidth, new TGeoRotation("",ang1, ang2, ang3)));//("",0, 0, 0)));
1644       }
1645       }
1646
1647       //*************
1648
1649     }
1650
1651     if (fBuildLevel < 3) {
1652
1653       //Bottom filament CFRP black-12 Carbon structure TGeoBBox (thickness,width,length)
1654
1655       Double_t filWidth = 0.04;
1656       Double_t filHeight= 0.02;
1657       TGeoBBox *t1=new TGeoBBox(filHeight/2,filWidth/2,s1);
1658       TGeoVolume *volT1=new TGeoVolume("bottomFilament", t1, medM60J3K);
1659       volT1->SetLineColor(12);
1660       volT1->SetFillColor(12);
1661       for(int i=0;i<loop;i++){//i<30;i++){
1662         mechLaddVol->AddNode(volT1,4*i+0,
1663                                     new TGeoCombiTrans(x+w,y-h+0.04+filHeight/2,z-kStaveLength/2+(4*l*i)+s1/2,
1664                                                        new TGeoRotation("volT1",-90,alpha,0)));
1665         mechLaddVol->AddNode(volT1,4*i+1,
1666                                     new TGeoCombiTrans(x-w,y-h+0.04+filHeight/2,z-kStaveLength/2+(4*l*i)+s1/2,
1667                                                        new TGeoRotation("volT1",90,alpha,0)));
1668         mechLaddVol->AddNode(volT1,4*i+2,
1669                                     new TGeoCombiTrans(x+w,y-h+0.04+filHeight/2,z-kStaveLength/2+2*l+(i*4*l)+s1/2,
1670                                                        new TGeoRotation("volT1",-90,-alpha,0)));
1671         mechLaddVol->AddNode(volT1,4*i+3,
1672                                     new TGeoCombiTrans(x-w,y-h+0.04+filHeight/2,z-kStaveLength/2+2*l+(i*4*l)+s1/2,
1673                                                        new TGeoRotation("volT1",-90,+alpha,0)));
1674         }
1675
1676      // Top filament CERP black-12 Carbon structure TGeoBBox (length,thickness,width)
1677
1678       TGeoBBox *t2=new TGeoBBox(s2,filHeight/2,filWidth/2);
1679       TGeoVolume *volT2=new TGeoVolume("topFilament", t2, medM60J3K);
1680       volT2->SetLineColor(12);
1681       volT2->SetFillColor(12);
1682       for(int i=0;i<loop;i++){ //i<30;i++){
1683         mechLaddVol->AddNode(volT2,4*i+0,
1684                                     new TGeoCombiTrans(x+w,y+0.04+filHeight/2,z-kStaveLength/2+(i*4*l)+s1/2,
1685                                                        new TGeoRotation("volT2",90,90-alpha,90-beta)));
1686         mechLaddVol->AddNode(volT2,4*i+1,
1687                                     new TGeoCombiTrans(x-w,y+0.04+filHeight/2,z-kStaveLength/2+(i*4*l)+s1/2,
1688                                                        new TGeoRotation("volT2",90,-90+alpha,-90+beta)));
1689         mechLaddVol->AddNode(volT2,4*i+2,
1690                                     new TGeoCombiTrans(x+w,y+0.04+filHeight/2,z-kStaveLength/2+2*l+(i*4*l)+s1/2,
1691                                                        new TGeoRotation("volT2",90,-90+alpha,90-beta)));
1692         mechLaddVol->AddNode(volT2,4*i+3,
1693                                     new TGeoCombiTrans(x-w,y+0.04+filHeight/2,z-kStaveLength/2+2*l+(i*4*l)+s1/2,
1694                                                        new TGeoRotation("volT2",90,90-alpha,-90+beta)));
1695         }
1696     }
1697
1698     if (fBuildLevel < 2) {
1699
1700       // Glue Filament and Silicon MicroChannel
1701       TGeoBBox *tM0=new TGeoBBox(xladMC/5, klay4/2, zlad);
1702       TGeoVolume *volTM0=new TGeoVolume("glueFM", tM0,medGlue );
1703       volTM0->SetLineColor(5);
1704       volTM0->SetFillColor(5);
1705       mechLaddVol->AddNode(volTM0, 0, new TGeoCombiTrans(x-xlad/2-0.25,0.03+yMC, z, new TGeoRotation("",0, 0, 0)));
1706       mechLaddVol->AddNode(volTM0, 1, new TGeoCombiTrans(x+xlad/2+0.25,0.03+yMC, z, new TGeoRotation("",0, 0, 0)));
1707
1708
1709       // Glue microchannel and sensor
1710       TGeoBBox *glueM = new TGeoBBox(xladMC/5, klay4/2, zlad);
1711       TGeoVolume *volGlueM=new TGeoVolume("glueMS", glueM, medGlue);
1712       volGlueM->SetLineColor(5);
1713       volGlueM->SetFillColor(5);
1714       mechLaddVol->AddNode(volGlueM, 0, new TGeoCombiTrans(x-xlad/2-0.25,yMC-0.01, z, new TGeoRotation("",0, 0, 0)));
1715       mechLaddVol->AddNode(volGlueM, 1, new TGeoCombiTrans(x+xlad/2+0.25,yMC-0.01, z, new TGeoRotation("",0, 0, 0)));
1716
1717        // Glue sensor and kapton
1718       TGeoBBox *glue = new TGeoBBox(xlad, klay4/2, zlad);
1719       TGeoVolume *volGlue=new TGeoVolume("glueSensorBus", glue, medGlue);
1720       volGlue->SetLineColor(5);
1721       volGlue->SetFillColor(5);
1722        mechLaddVol->AddNode(volGlue, 1, new TGeoCombiTrans(x, y-0.154-fSensorThick-klay4/2, z, new TGeoRotation("",0, 0, 0)));
1723     }
1724
1725     if (fBuildLevel < 1) {
1726       TGeoBBox *kapCable = new TGeoBBox(xlad, klay5/2, zlad);
1727       TGeoVolume *volCable=new TGeoVolume("Flexcable", kapCable, medFlexCable);
1728       volCable->SetLineColor(28);
1729       volCable->SetFillColor(28);
1730       mechLaddVol->AddNode(volCable, 0, new TGeoCombiTrans(x, y-0.154-fSensorThick-klay4-klay5/2, z, new TGeoRotation("",0, 0, 0)));
1731     }
1732
1733   // Done, return the stave structur
1734     return mechLaddVol;
1735  }
1736
1737 //________________________________________________________________________
1738 TGeoVolume* AliITSUv0Layer::CreateModule(const Double_t xlad,
1739                                          const Double_t ylad,
1740                                          const Double_t zlad,
1741                                          const TGeoManager *mgr){
1742 //
1743 // Creates the actual Module
1744 //
1745 // Input:
1746 //         xlad,zlad : the ladder dimensions
1747 //         mgr  : the GeoManager (used only to get the proper material)
1748 //
1749 // Output:
1750 //
1751 // Return:
1752 //
1753 // Created:      22 Jun 2011  Mario Sitta
1754 //
1755
1756   char volname[30];
1757   Double_t xlen, ylen, zlen;
1758   Double_t xpos, ypos, zpos;
1759
1760
1761   // First create all needed shapes
1762
1763   // The module
1764   TGeoBBox *module = new TGeoBBox(xlad,  ylad, zlad/fNModules);
1765
1766   // The sensor
1767   xlen = module->GetDX();
1768   ylen = 0.5*fSensorThick;
1769   zlen = module->GetDZ();
1770   TGeoBBox *sensor = new TGeoBBox(xlen, ylen, zlen);
1771
1772
1773   // We have all shapes: now create the real volumes
1774   //TGeoMedium *medAir = mgr->GetMedium("ITS_AIR$");
1775   TGeoMedium *medSi  = mgr->GetMedium("ITS_SI$");
1776
1777   snprintf(volname, 30, "%s%d", AliITSUGeomTGeo::GetITSModulePattern(), fLayerNumber);
1778   // TGeoVolume *modVol = new TGeoVolume(volname, module, medAir);
1779   TGeoVolume *modVol = new TGeoVolume(volname, module, medSi);
1780   modVol->SetVisibility(kTRUE);
1781   modVol->SetLineColor(1);
1782
1783   snprintf(volname, 30, "%s%d", AliITSUGeomTGeo::GetITSSensorPattern(), fLayerNumber);
1784   TGeoVolume *sensVol = new TGeoVolume(volname, sensor, medSi);
1785   sensVol->SetVisibility(kTRUE);
1786   sensVol->SetLineColor(8);
1787   sensVol->SetLineWidth(1);
1788   sensVol->SetFillColor(sensVol->GetLineColor());
1789   sensVol->SetFillStyle(4000); // 0% transparent
1790
1791
1792   // Now build up the module
1793   xpos = 0.;
1794   ypos = -module->GetDY() + sensor->GetDY();
1795   zpos = 0.;
1796
1797   modVol->AddNode(sensVol, 1, new TGeoTranslation(xpos, ypos, zpos));
1798
1799   // Done, return the module
1800   return modVol;
1801 }
1802
1803 //________________________________________________________________________
1804 Double_t AliITSUv0Layer::RadiusOfTurboContainer(){
1805 //
1806 // Computes the inner radius of the air container for the Turbo configuration
1807 // as the radius of either the circle tangent to the ladder or the circle
1808 // passing for the ladder's lower vertex
1809 //
1810 // Input:
1811 //         none (all needed parameters are class members)
1812 //
1813 // Output:
1814 //
1815 // Return:
1816 //        the radius of the container if >0, else flag to use the lower vertex
1817 //
1818 // Created:      08 Mar 2012  Mario Sitta
1819 //
1820
1821   Double_t rr, delta, z, lladd, rladd;
1822
1823   if (fLadderThick > 89.) // Very big angle: avoid overflows since surely
1824     return -1;            // the radius from lower vertex is the right value
1825
1826   rladd = fLayRadius + 0.5*fLadderThick;
1827   delta = (0.5*fLadderThick)/CosD(fLadderTilt);
1828   z     = (0.5*fLadderThick)*TanD(fLadderTilt);
1829
1830   rr = rladd - delta;
1831   lladd = (0.5*fLadderWidth) - z;
1832
1833   if ( (rr*SinD(fLadderTilt) < lladd) )
1834     return (rr*CosD(fLadderTilt));
1835   else
1836     return -1;
1837 }
1838
1839 //________________________________________________________________________
1840 void AliITSUv0Layer::SetLadderTilt(const Double_t t)
1841 {
1842 //
1843 // Sets the Ladder tilt angle (for turbo layers only)
1844 //
1845 // Input:
1846 //         t :  the ladder tilt angle
1847 //
1848 // Output:
1849 //
1850 // Return:
1851 //
1852 // Created:      08 Jul 2011  Mario Sitta
1853 //
1854
1855   if (fIsTurbo)
1856     fLadderTilt = t;
1857   else
1858     AliError("Not a Turbo layer");
1859
1860 }
1861
1862 //________________________________________________________________________
1863 void AliITSUv0Layer::SetLadderWidth(const Double_t w){
1864 //
1865 // Sets the Ladder width (for turbo layers only)
1866 //
1867 // Input:
1868 //         w :  the ladder width
1869 //
1870 // Output:
1871 //
1872 // Return:
1873 //
1874 // Created:      08 Jul 2011  Mario Sitta
1875 //
1876
1877   if (fIsTurbo)
1878     fLadderWidth = w;
1879   else
1880     AliError("Not a Turbo layer");
1881
1882 }