Transition to NewIO
[u/mrichter/AliRoot.git] / TPC / AliTPCv0.cxx
CommitLineData
4c039060 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
88cb7938 16/* $Id$ */
73042f01 17//
4c039060 18
fe4da5cc 19///////////////////////////////////////////////////////////////////////////////
20// //
21// Time Projection Chamber version 0 -- "coarse" TPC //
22// //
23//Begin_Html
24/*
1439f98e 25<img src="picts/AliTPCv0Class.gif">
fe4da5cc 26*/
27//End_Html
28// //
29// //
30///////////////////////////////////////////////////////////////////////////////
31
88cb7938 32#include <Riostream.h>
33
fe4da5cc 34#include <TGeometry.h>
88cb7938 35#include <TMath.h>
fe4da5cc 36#include <TNode.h>
37#include <TTUBE.h>
88cb7938 38#include <TVirtualMC.h>
fe4da5cc 39
88cb7938 40#include "AliConst.h"
41#include "AliRun.h"
cc80f89e 42#include "AliTPCDigitsArray.h"
88cb7938 43#include "AliTPCParamSR.h"
44#include "AliTPCv0.h"
1283eee5 45
fe4da5cc 46ClassImp(AliTPCv0)
47
48//_____________________________________________________________________________
49AliTPCv0::AliTPCv0(const char *name, const char *title)
50 :AliTPC(name, title)
51{
52 //
53 // Standard creator for TPC version 0
54 //
73042f01 55
56 SetGasMixt(2,20,10,-1,0.9,0.1,0.);
57
58 if (fTPCParam)
59 fTPCParam->Write(fTPCParam->GetTitle());
fe4da5cc 60}
61
62//_____________________________________________________________________________
63void AliTPCv0::CreateGeometry()
64{
79575d04 65
fe4da5cc 66 //
37831078 67 // Creation of the TPC version 0, i.e. no sensitive volumes,
68 // only the material distribution
fe4da5cc 69 //
70 //Begin_Html
71 /*
37831078 72 <img src="picts/AliTPC.gif">
fe4da5cc 73 */
74 //End_Html
75 //Begin_Html
76 /*
1439f98e 77 <img src="picts/AliTPCv0Tree.gif">
fe4da5cc 78 */
79 //End_Html
80
37831078 81 //---------------------------------------------------------------
82 // TPC geometry
83 //
84 // Origin: Marek Kowalski, INP Cracow
85 // Marek.Kowalski@ifj.edu.pl
86 //---------------------------------------------------------------
37831078 87 Float_t dm[50];
1283eee5 88
37831078 89 Int_t *idtmed = fIdtmed->GetArray(); // TPC media
1283eee5 90
37831078 91 Int_t idrotm[120]; // rotation matrices
1283eee5 92
37831078 93 Int_t nRotMat = 0; // actual number of rotation matrices
1283eee5 94
1283eee5 95
37831078 96 //
79575d04 97 // Mother volume TPC (Air) - all volumes will be positioned in it
37831078 98 //
79575d04 99
37831078 100 dm[0]=0.;
101 dm[1]=360.;
79575d04 102 dm[2]=8.;
1283eee5 103
37831078 104 //
1283eee5 105
79575d04 106 dm[3]=-283.7;
107 dm[4]=65.6;
108 dm[5]=278.;
1283eee5 109
37831078 110 //
1283eee5 111
79575d04 112 dm[6]=-253.6;
113 dm[7]=65.6;
114 dm[8]=278.;
1283eee5 115
37831078 116 //
fe4da5cc 117
79575d04 118 dm[9]=-73.3;
119 dm[10]=60.9;
120 dm[11]=278.;
1283eee5 121
37831078 122 //
1283eee5 123
79575d04 124 dm[12]=-73.3;
125 dm[13]=56.9;
126 dm[14]=278.;
1283eee5 127
37831078 128 //
1283eee5 129
79575d04 130 dm[15]=73.3;
131 dm[16]=56.9;
132 dm[17]=278.;
1283eee5 133
37831078 134 //
1283eee5 135
79575d04 136 dm[18]=73.3;
137 dm[19]=60.9;
138 dm[20]=278.;
1283eee5 139
37831078 140 //
1283eee5 141
79575d04 142 dm[21]=253.6;
143 dm[22]=65.6;
144 dm[23]=278.;
1283eee5 145
37831078 146 //
1283eee5 147
79575d04 148 dm[24]=283.7;
149 dm[25]=65.6;
150 dm[26]=278.;
151
152 gMC->Gsvolu("TPC ","PCON",idtmed[0],dm,27);
1283eee5 153
79575d04 154 // outer part
1283eee5 155
37831078 156 //-------------------------------------------------------------------
79575d04 157 // Tpc Outer INsulator (CO2) - contains cont. vessel and field cage
37831078 158 //-------------------------------------------------------------------
1283eee5 159
37831078 160 dm[0]= 0.;
161 dm[1]= 360.;
162 dm[2]= 6.;
1283eee5 163
37831078 164 //
1283eee5 165
79575d04 166 dm[3]=-253.6;
167 dm[4]=258.;
168 dm[5]=275.5;
1283eee5 169
37831078 170 //
cc80f89e 171
79575d04 172 dm[6]=-250.6;
173 dm[7]=258.;
174 dm[8]=275.5;
1283eee5 175
37831078 176 //
1283eee5 177
79575d04 178 dm[9]=-250.6;
179 dm[10]=258.;
180 dm[11]=278.;
1283eee5 181
79575d04 182 //
1283eee5 183
79575d04 184 dm[12]=253.6;
185 dm[13]=258.;
186 dm[14]=278.;
1283eee5 187
37831078 188 //
1283eee5 189
79575d04 190 dm[15]=253.6;
191 dm[16]=264.8;
192 dm[17]=278.;
1283eee5 193
37831078 194 //
1283eee5 195
79575d04 196 dm[18]=256.6;
197 dm[19]=264.8;
198 dm[20]=278.;
1283eee5 199
79575d04 200 gMC->Gsvolu("TOIN","PCON",idtmed[3],dm,21);
1283eee5 201
37831078 202 //----------------------------------------------------------------
203 // Tpc Outer Contaiment Vessel
204 // mother volume - Al, daughters - composite (sandwich)
205 //----------------------------------------------------------------
79575d04 206
207 dm[0]=0.;
208 dm[1]=360.;
720b23b8 209 dm[2]=6.;
1283eee5 210
37831078 211 //
1283eee5 212
79575d04 213 dm[3]=-250.6;
214 dm[4]=270.4;
215 dm[5]=278.;
37831078 216
217 //
218
79575d04 219 dm[6]=-247.6;
220 dm[7]=270.4;
221 dm[8]=278.;
37831078 222
223 //
79575d04 224
225 dm[9]=-247.6;
cea9b4f7 226 dm[10]=274.8124;
79575d04 227 dm[11]=278.;
37831078 228
229 //
230
79575d04 231 dm[12]=253.6;
cea9b4f7 232 dm[13]=274.8124;
79575d04 233 dm[14]=278.;
37831078 234
720b23b8 235 //
37831078 236
79575d04 237 dm[15]=253.6;
238 dm[16]=264.8;
239 dm[17]=278.;
37831078 240
79575d04 241 //
720b23b8 242
79575d04 243 dm[18]=256.6;
244 dm[19]=264.8;
245 dm[20]=278.;
720b23b8 246
247 gMC->Gsvolu("TOCV","PCON",idtmed[4],dm,21);
37831078 248
79575d04 249 // Daughter volumes - sandwich
37831078 250
251 // Tpc SAndwich 1 - Al
79575d04 252
cea9b4f7 253 dm[0]=274.8124;
79575d04 254 dm[1]=278.;
255 dm[2]=252.1;
1283eee5 256
37831078 257 gMC->Gsvolu("TSA1","TUBE",idtmed[4],dm,3);
1283eee5 258
79575d04 259 // Tpc SAndwich 2 - epoxy glue (I use Lexan)
1283eee5 260
37831078 261 dm[0] += 5.e-3;
262 dm[1] -= 5.e-3;
79575d04 263
264 gMC->Gsvolu("TSA2","TUBE",idtmed[14],dm,3);
265
266 // Tpc SAndwich 3 - Tedlar
267
268 dm[0] += 0.01;
269 dm[1] -= 0.01;
37831078 270
79575d04 271 gMC->Gsvolu("TSA3","TUBE",idtmed[9],dm,3);
1283eee5 272
1283eee5 273
79575d04 274 // Tpc SAndwich 4 - fiber glass (G10)
1283eee5 275
cea9b4f7 276 dm[0] += 3.8e-3;
277 dm[1] -= 3.8e-3;
1283eee5 278
79575d04 279 gMC->Gsvolu("TSA4","TUBE",idtmed[12],dm,3);
37831078 280
79575d04 281 // Tpc SAndwich 5 - NOMEX honeycomb
37831078 282
79575d04 283 dm[0] += 0.075;
284 dm[1] -= 0.075;
37831078 285
79575d04 286 gMC->Gsvolu("TSA5","TUBE",idtmed[6],dm,3);
287
288 // 5->4->3->2->1->TCOV
37831078 289
79575d04 290
291 gMC->Gspos("TSA5",1,"TSA4",0.,0.,0.,0,"ONLY");
37831078 292 gMC->Gspos("TSA4",1,"TSA3",0.,0.,0.,0,"ONLY");
293 gMC->Gspos("TSA3",1,"TSA2",0.,0.,0.,0,"ONLY");
79575d04 294 gMC->Gspos("TSA2",1,"TSA1",0.,0.,0.,0,"ONLY");
1283eee5 295
79575d04 296 gMC->Gspos("TSA1",1,"TOCV",0.,0.,3.,0,"ONLY");
1283eee5 297
37831078 298 // TCOV-> TOIN
299
300 gMC->Gspos("TOCV",1,"TOIN",0.,0.,0.,0,"ONLY");
301
302 //-------------------------------------------------------
303 // Tpc Outer Field Cage
304 // mother volume - Al, daughters - composite (sandwich)
305 //-------------------------------------------------------
306
307 dm[0]=0.;
308 dm[1]=360.;
309 dm[2]=6.;
310
79575d04 311 //
312
313 dm[3]=-253.6;
314 dm[4]=258.;
315 dm[5]=275.5;
37831078 316
317 //
318
79575d04 319 dm[6]=-250.6;
320 dm[7]=258.;
321 dm[8]=275.5;
37831078 322
323 //
79575d04 324
325 dm[9]=-250.6;
326 dm[10]=258.;
cea9b4f7 327 dm[11]=260.0476;
1283eee5 328
37831078 329 //
1283eee5 330
79575d04 331 dm[12]=250.6;
332 dm[13]=258.;
cea9b4f7 333 dm[14]=260.0476;
1283eee5 334
37831078 335 //
336
79575d04 337 dm[15]=250.6;
338 dm[16]=258.;
dcf38560 339 dm[17]=269.6;
37831078 340
341 //
342
79575d04 343 dm[18]=253.6;
344 dm[19]=258.;
dcf38560 345 dm[20]=269.6;
1283eee5 346
79575d04 347 gMC->Gsvolu("TOFC","PCON",idtmed[4],dm,21);
1283eee5 348
37831078 349 // Daughter volumes
1283eee5 350
79575d04 351 // Tpc SAndwich 6 - Tedlar
37831078 352
353 dm[0]= 258.;
cea9b4f7 354 dm[1]= 260.0476;
79575d04 355 dm[2]= 252.1;
37831078 356
79575d04 357 gMC->Gsvolu("TSA6","TUBE",idtmed[9],dm,3);
37831078 358
79575d04 359 // Tpc SAndwich 7 - fiber glass
37831078 360
cea9b4f7 361 dm[0] += 3.8e-3;
362 dm[1] -= 3.8e-3;
37831078 363
79575d04 364 gMC->Gsvolu("TSA7","TUBE",idtmed[12],dm,3);
37831078 365
37831078 366
79575d04 367 // Tpc SAndwich 8 - NOMEX
37831078 368
369 dm[0] += 0.02;
370 dm[1] -= 0.02;
371
79575d04 372 gMC->Gsvolu("TSA8","TUBE",idtmed[6],dm,3);
37831078 373
79575d04 374 // 8->7->6->TOFC
37831078 375
79575d04 376 gMC->Gspos("TSA8",1,"TSA7",0.,0.,0.,0,"ONLY");
377 gMC->Gspos("TSA7",1,"TSA6",0.,0.,0.,0,"ONLY");
378 gMC->Gspos("TSA6",1,"TOFC",0.,0.,0.,0,"ONLY");
37831078 379
380 // TOFC->TOIN
381
79575d04 382 gMC->Gspos("TOFC",1,"TOIN",0.,0.,0.,0,"ONLY");
37831078 383
384 // TOIN->TPC
385
386 gMC->Gspos("TOIN",1,"TPC ",0.,0.,0.,0,"ONLY");
387
79575d04 388 // inner part
37831078 389
79575d04 390 //--------------------------------------------------------------------
391 // Tpc Inner INsulator (CO2) - inner f.c. will be placed there
392 // Inner containment vessel will be placed directly in the TPC
393 //--------------------------------------------------------------------
37831078 394
395 dm[0]=0.;
79575d04 396 dm[1]=360.;
397 dm[2]=4.;
1283eee5 398
79575d04 399 //
1283eee5 400
79575d04 401 dm[3]=-253.6;
402 dm[4]=65.9;
403 dm[5]=79.2;
1283eee5 404
37831078 405 //
1283eee5 406
79575d04 407 dm[6]=-73.3;
408 dm[7]=61.2;
409 dm[8]=79.2;
1283eee5 410
37831078 411 //
1283eee5 412
79575d04 413 dm[9]=73.3;
414 dm[10]=61.2;
415 dm[11]=79.2;
1283eee5 416
37831078 417 //
1283eee5 418
79575d04 419 dm[12]=253.6;
420 dm[13]=65.9;
421 dm[14]=79.2;
1283eee5 422
79575d04 423 gMC->Gsvolu("TIIN","PCON",idtmed[3],dm,15);
37831078 424
79575d04 425 // the middle part of the F.C. is thinner - carve out the strip - Ne-CO2
1283eee5 426
79575d04 427 dm[0]=79.16;
428 dm[1]=79.2;
429 dm[2]=88.;
1283eee5 430
79575d04 431 gMC->Gsvolu("TII1","TUBE",idtmed[1],dm,3);
1283eee5 432
79575d04 433 gMC->Gspos("TII1",1,"TIIN",0.,0.,0.,0,"ONLY");
1283eee5 434
79575d04 435 //-----------------------------------------------------
436 // Tpc Inner Field Cage
437 // mother volume - Al, daughters - composite (sandwich)
438 //------------------------------------------------------
1283eee5 439
79575d04 440 dm[0]=0.;
441 dm[1]=360.;
442 dm[2]=10.;
1283eee5 443
37831078 444 //
1283eee5 445
79575d04 446 dm[3]=-253.6;
447 dm[4]=70.3;
448 dm[5]=79.2;
1283eee5 449
37831078 450 //
1283eee5 451
79575d04 452 dm[6]=-250.6;
453 dm[7]=70.3;
454 dm[8]=79.2;
1283eee5 455
37831078 456 //
1283eee5 457
79575d04 458 dm[9]=-250.6;
cea9b4f7 459 dm[10]=77.0524;
79575d04 460 dm[11]=79.2;
1283eee5 461
37831078 462 //
1283eee5 463
79575d04 464 dm[12]=-88.;
cea9b4f7 465 dm[13]=77.0524;
79575d04 466 dm[14]=79.2;
1283eee5 467
37831078 468 //
1283eee5 469
79575d04 470 dm[15]=-88.;
cea9b4f7 471 dm[16]=77.0924;
79575d04 472 dm[17]=79.16;
1283eee5 473
37831078 474 //
1283eee5 475
79575d04 476 dm[18]=88.;
cea9b4f7 477 dm[19]=77.0924;
79575d04 478 dm[20]=79.16;
1283eee5 479
37831078 480 //
1283eee5 481
79575d04 482 dm[21]=88.;
cea9b4f7 483 dm[22]=77.0524;
79575d04 484 dm[23]=79.2;
1283eee5 485
37831078 486 //
1283eee5 487
79575d04 488 dm[24]=250.6;
cea9b4f7 489 dm[25]=77.0524;
79575d04 490 dm[26]=79.2;
1283eee5 491
37831078 492 //
1283eee5 493
79575d04 494 dm[27]=250.6;
495 dm[28]=70.3;
496 dm[29]=79.2;
37831078 497
498 //
499
79575d04 500 dm[30]=253.6;
501 dm[31]=70.3;
502 dm[32]=79.2;
37831078 503
79575d04 504 gMC->Gsvolu("TIFC","PCON",idtmed[4],dm,33);
37831078 505
79575d04 506 // daughter volumes - central part
37831078 507
79575d04 508 // Tpc Sandwich 9 -Tedlar
37831078 509
cea9b4f7 510 dm[0]=77.0924;
79575d04 511 dm[1]=79.16;
512 dm[2]=88.;
37831078 513
79575d04 514 gMC->Gsvolu("TSA9","TUBE",idtmed[9],dm,3);
1283eee5 515
79575d04 516 // Tpc Sandwich 10 - fiber glass (G10)
fe4da5cc 517
cea9b4f7 518 dm[0] += 3.8e-3;
519 dm[1] -= 3.8e-3;
1283eee5 520
79575d04 521 gMC->Gsvolu("TS10","TUBE",idtmed[12],dm,3);
1283eee5 522
79575d04 523 // Tpc Sandwich 11 - NOMEX
1283eee5 524
79575d04 525 dm[0] += 0.03;
526 dm[1] -= 0.03;
1283eee5 527
79575d04 528 gMC->Gsvolu("TS11","TUBE",idtmed[6],dm,3);
37831078 529
79575d04 530 // 11->10->9->TIFC
1283eee5 531
79575d04 532 gMC->Gspos("TS11",1,"TS10",0.,0.,0.,0,"ONLY");
533 gMC->Gspos("TS10",1,"TSA9",0.,0.,0.,0,"ONLY");
37831078 534
79575d04 535 gMC->Gspos("TSA9",1,"TIFC",0.,0.,0.,0,"ONLY");
37831078 536
79575d04 537 // daughter volumes - outer parts (reinforced)
37831078 538
79575d04 539 // Tpc Sandwich 12 -Tedlar
37831078 540
cea9b4f7 541 dm[0]=77.0524;
79575d04 542 dm[1]=79.2;
543 dm[2]=82.05;
37831078 544
79575d04 545 gMC->Gsvolu("TS12","TUBE",idtmed[9],dm,3);
37831078 546
79575d04 547 // Tpc Sandwich 13 - fiber glass (G10)
37831078 548
cea9b4f7 549 dm[0] += 3.8e-3;
550 dm[1] -= 3.8e-3;
37831078 551
79575d04 552 gMC->Gsvolu("TS13","TUBE",idtmed[12],dm,3);
37831078 553
79575d04 554 // Tpc Sandwich 14 - NOMEX
37831078 555
79575d04 556 dm[0] += 0.07;
557 dm[1] -= 0.07;
37831078 558
79575d04 559 gMC->Gsvolu("TS14","TUBE",idtmed[6],dm,3);
37831078 560
79575d04 561 // 14->13->12->TIFC
37831078 562
79575d04 563 gMC->Gspos("TS14",1,"TS13",0.,0.,0.,0,"ONLY");
564 gMC->Gspos("TS13",1,"TS12",0.,0.,0.,0,"ONLY");
37831078 565
79575d04 566 gMC->Gspos("TS12",1,"TIFC",0.,0.,170.05,0,"ONLY");
567 gMC->Gspos("TS12",2,"TIFC",0.,0.,-170.05,0,"ONLY");
37831078 568
79575d04 569 // place this inside the inner insulator
37831078 570
79575d04 571 gMC->Gspos("TIFC",1,"TIIN",0.,0.,0.,0,"ONLY");
37831078 572
79575d04 573 // and now in the TPC...
37831078 574
79575d04 575 gMC->Gspos("TIIN",1,"TPC ",0.,0.,0.,0,"ONLY");
37831078 576
79575d04 577 //---------------------------------------------------------
578 // Tpc Inner Containment vessel - Cones
579 //---------------------------------------------------------
37831078 580
79575d04 581 dm[0]=0.;
582 dm[1]=360.;
37831078 583 dm[2]=8.;
584
585 //
586
79575d04 587 dm[3]=71.8;
588 dm[4]=56.9;
589 dm[5]=59.4;
37831078 590
591 //
592
79575d04 593 dm[6]=73.;
594 dm[7]=56.9;
595 dm[8]=59.4;
37831078 596
597 //
598
79575d04 599 dm[9]=73.;
600 dm[10]=56.9;
601 dm[11]=61.2;
37831078 602
603 //
604
79575d04 605 dm[12]=73.3;
606 dm[13]=56.9;
607 dm[14]=61.2;
37831078 608
609 //
79575d04 610
611 dm[15]=73.3;
612 dm[16]=60.9;
613 dm[17]=61.2;
37831078 614
79575d04 615 //
37831078 616
79575d04 617 dm[18]=253.6;
618 dm[19]=65.6;
619 dm[20]=65.9;
37831078 620
621 //
622
79575d04 623 dm[21]=253.6;
624 dm[22]=65.6;
f4900c61 625 dm[23]=74.6;
37831078 626
627 //
628
f4900c61 629 dm[24]=256.6;
79575d04 630 dm[25]=65.6;
f4900c61 631 dm[26]=74.6;
37831078 632
79575d04 633 gMC->Gsvolu("TICC","PCON",idtmed[4],dm,27);
37831078 634
79575d04 635 Float_t phi1,phi2,phi3,theta1,theta2,theta3; // rotation angles
37831078 636
79575d04 637 // reflection matrix
638
639 theta1 = 90.;
640 phi1 = 0.;
641 theta2 = 90.;
642 phi2 = 270.;
643 theta3 = 180.;
644 phi3 = 0.;
37831078 645
79575d04 646 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
37831078 647
79575d04 648 gMC->Gspos("TICC",1,"TPC ",0.,0.,0.,0,"ONLY");
649 gMC->Gspos("TICC",2,"TPC ",0.,0.,0.,idrotm[nRotMat],"ONLY");
37831078 650
37831078 651
79575d04 652 //---------------------------------------------------------
653 // Tpc Inner Containment vessel - Middle part -Al
654 //---------------------------------------------------------
37831078 655
79575d04 656 dm[0]=0.;
657 dm[1]=360.;
658 dm[2]=6.;
37831078 659
660 //
661
79575d04 662 dm[3]=-71.6;
663 dm[4]=60.2;
664 dm[5]=61.2;
37831078 665
666 //
667
79575d04 668 dm[6]=-69.1;
669 dm[7]=60.2;
670 dm[8]=61.2;
37831078 671
672 //
673
79575d04 674 dm[9]=-69.1;
cea9b4f7 675 dm[10]=60.6224;
79575d04 676 dm[11]=61.2;
37831078 677
678 //
679
79575d04 680 dm[12]=69.1;
cea9b4f7 681 dm[13]=60.6224;
79575d04 682 dm[14]=61.2;
37831078 683
684 //
685
79575d04 686 dm[15]=69.1;
687 dm[16]=60.2;
688 dm[17]=61.2;
37831078 689
690 //
691
79575d04 692 dm[18]=71.6;
693 dm[19]=60.2;
694 dm[20]=61.2;
37831078 695
79575d04 696 gMC->Gsvolu("TICM","PCON",idtmed[4],dm,21);
37831078 697
79575d04 698 // Tpc Sandwich 15 - Al
37831078 699
cea9b4f7 700 dm[0]=60.6224;
79575d04 701 dm[1]=61.2;
702 dm[2]=70.1;
37831078 703
79575d04 704 gMC->Gsvolu("TS15","TUBE",idtmed[4],dm,3);
37831078 705
79575d04 706 // Tpc Sandwich 16 - epoxy glue
37831078 707
79575d04 708 dm[0] += 5.e-3;
709 dm[1] -= 5.e-3;
37831078 710
79575d04 711 gMC->Gsvolu("TS16","TUBE",idtmed[14],dm,3);
37831078 712
6a905708 713 // Tpc Sandwich 17 - Tedlar
37831078 714
79575d04 715 dm[0] += 0.01;
716 dm[1] -= 0.01;
37831078 717
6a905708 718 gMC->Gsvolu("TS17","TUBE",idtmed[9],dm,3);
37831078 719
cea9b4f7 720 // Tpc Sandwich 18 - carbon fiber
37831078 721
cea9b4f7 722 dm[0] += 3.8e-3;
723 dm[1] -= 3.8e-3;
37831078 724
cea9b4f7 725 gMC->Gsvolu("TS18","TUBE",idtmed[15],dm,3);
37831078 726
79575d04 727 // Tpc Sandwich 19 - Nomex
37831078 728
cea9b4f7 729 dm[0] += 0.02;
730 dm[1] -= 0.02;
37831078 731
90222263 732 gMC->Gsvolu("TS19","TUBE",idtmed[6],dm,3);
37831078 733
79575d04 734 // 19->18->17->16->15-> TICM
37831078 735
79575d04 736 gMC->Gspos("TS19",1,"TS18",0.,0.,0.,0,"ONLY");
37831078 737 gMC->Gspos("TS18",1,"TS17",0.,0.,0.,0,"ONLY");
79575d04 738 gMC->Gspos("TS17",1,"TS16",0.,0.,0.,0,"ONLY");
739 gMC->Gspos("TS16",1,"TS15",0.,0.,0.,0,"ONLY");
37831078 740
79575d04 741 gMC->Gspos("TS15",1,"TICM ",0.,0.,0.,0,"ONLY");
742
37831078 743
79575d04 744 // TPc inner cont. vessel Joints
37831078 745
79575d04 746 dm[0]=60.2;
747 dm[1]=61.2;
748 dm[2]=0.5;
37831078 749
79575d04 750 gMC->Gsvolu("TPJ1","TUBE",idtmed[4],dm,3);
37831078 751
79575d04 752 gMC->Gspos("TPJ1",1,"TPC ",0.,0.,72.3,0,"ONLY");
753 gMC->Gspos("TPJ1",2,"TPC ",0.,0.,-72.3,0,"ONLY");
37831078 754
79575d04 755 //
37831078 756
79575d04 757 dm[0]=0.;
758 dm[1]=360.;
759 dm[2]=4.;
37831078 760
79575d04 761 //
37831078 762
79575d04 763 dm[3]=70.8;
764 dm[4]=58.4;
765 dm[5]=60.1;
37831078 766
79575d04 767 //
37831078 768
79575d04 769 dm[6]=71.2;
770 dm[7]=58.4;
771 dm[8]=60.1;
37831078 772
79575d04 773 //
37831078 774
79575d04 775 dm[9]=71.2;
776 dm[10]=58.4;
777 dm[11]=59.4;
37831078 778
79575d04 779 //
37831078 780
79575d04 781 dm[12]=71.6;
782 dm[13]=58.4;
783 dm[14]=59.4;
37831078 784
79575d04 785 gMC->Gsvolu("TPJ2","PCON",idtmed[4],dm,15);
37831078 786
79575d04 787 gMC->Gspos("TPJ2",1,"TPC ",0.,0.,0.,0,"ONLY");
788 gMC->Gspos("TPJ2",2,"TPC ",0.,0.,0.,idrotm[nRotMat],"ONLY");
37831078 789
37831078 790
37831078 791
79575d04 792 // Tpc Inner Containment vessel Seal (Viton, I use Lexan for a time being)
37831078 793
79575d04 794 dm[0]=58.4;
795 dm[1]=61.2;
796 dm[2]=0.1;
37831078 797
79575d04 798 gMC->Gsvolu("TICS","TUBE",idtmed[14],dm,3);
37831078 799
79575d04 800 gMC->Gspos("TICS",1,"TPC ",0.,0.,71.7,0,"ONLY");
801 gMC->Gspos("TICS",2,"TPC ",0.,0.,-71.7,0,"ONLY");
37831078 802
79575d04 803 // TICM -> TPC
37831078 804
79575d04 805 gMC->Gspos("TICM",1,"TPC ",0.,0.,0.,0,"ONLY");
37831078 806
79575d04 807 //
37831078 808
79575d04 809 nRotMat++; // prepare for the next rotation matrix
37831078 810
1283eee5 811 //---------------------------------------------------------
37831078 812 // Tpc Dift Gas volume Nonsensitive (Ne-CO2 90/10)
813 // and its daughters (HV membrane, rods, readout chambers)
1283eee5 814 //---------------------------------------------------------
815
37831078 816 dm[0]= 79.2;
817 dm[1]= 258.0;
818 dm[2]= 253.6;
1283eee5 819
79575d04 820 gMC->Gsvolu("TDGN","TUBE",idtmed[1],dm,3);
1283eee5 821
37831078 822 // sector opening angles
1283eee5 823
37831078 824 Float_t innerOpenAngle = fTPCParam->GetInnerAngle();
825
826 // sector angle shift
827
828 Float_t innerAngleShift = fTPCParam->GetInnerAngleShift();
1283eee5 829
37831078 830 // number of sectors
831
832 Int_t nInnerSector = fTPCParam->GetNInnerSector()/2;
833 Int_t nOuterSector = fTPCParam->GetNOuterSector()/2;
834
835 // All above parameters are identical for inner and outer
836 // sectors. The distinction is kept for the historical reasons
837 // and eventually will disappear.
838
839 Float_t tanAlpha = TMath::Tan(0.5*innerOpenAngle);
840 Float_t cosAlpha = TMath::Sqrt(1.+tanAlpha*tanAlpha);
841 Float_t space;
842
843 //-------------------------------------------------------------------------
844 // Tpc Inner Readout Chambers
845 //-------------------------------------------------------------------------
846
847 dm[0]= 14.483;
848 dm[1]= 23.3345;
849 dm[2]= 1.6; // thickness
850 dm[3]= 25.1;
851
852 gMC->Gsvolu("TIRC","TRD1",idtmed[4],dm,4);
853
854 // this volume will be positioned in the empty space
855 // of the end-cap to avoid overlaps
856
857 dm[0]= 13.7305;
858 dm[1]= 21.1895;
859 dm[2]= 2.25;
860 dm[3]= 21.15;
861
862 gMC->Gsvolu("TIC1","TRD1",idtmed[4],dm,4);
863
864
865 //------------------------------------------------
866 // Tpc Inner readout chamber Pad Plane
867 //------------------------------------------------
868
869 dm[0]= 14.483;
870 dm[1]= 23.3345;
871 dm[2]= 0.5;
872 dm[3]= 25.1;
873
874 gMC->Gsvolu("TIPP","TRD1",idtmed[12],dm,4);
875
876 //
1283eee5 877
37831078 878 dm[0] -= 1.218511934;
879 dm[1] -= 1.218511934;
880 dm[2] = 0.35;
1283eee5 881
37831078 882 gMC->Gsvolu("TIC3","TRD1",idtmed[1],dm,4);
1283eee5 883
37831078 884 gMC->Gspos("TIC3",1,"TIPP",0.,0.15,0.,0,"ONLY");
1283eee5 885
37831078 886 gMC->Gspos("TIPP",1,"TIRC",0.,1.1,0.,0,"ONLY");
1283eee5 887
1283eee5 888
37831078 889 //----------------------------------------------
890 // Tpc Readout Chambers Empty spaces - for both
891 // inner and outer sectors
892 //----------------------------------------------
1283eee5 893
37831078 894 gMC->Gsvolu("TRCE","TRD1",idtmed[0],dm,0);
1283eee5 895
37831078 896 // Inner sector - 4 spaces
897
898
899 dm[3] = 4.7625;
900 dm[0] = 12.472;
901
902 Float_t rr = 90.52;
903 Float_t zz;
904
905 zz= -12.7875;
906
907 space = rr*tanAlpha-dm[0];
908
909 for(Int_t nsLow=0;nsLow<4;nsLow++){
910
911 rr += 9.525;
912 dm[1]= rr*tanAlpha - space;
913
914 dm[2]=0.8;
915
916 gMC->Gsposp("TRCE",nsLow+1,"TIRC",0.,-0.8,zz,0,"ONLY",dm,4);
917
918 //
919
920 dm[2] = 1.2;
921
922 gMC->Gsposp("TRCE",nsLow+5,"TIC1",0.,1.05,zz-2.1,0,"ONLY",dm,4);
923
924 rr += 0.4;
925 dm[0] = rr*tanAlpha - space;
926 zz += (0.4+9.525);
1283eee5 927
fe4da5cc 928 }
1283eee5 929
37831078 930 dm[0]= 12.472;
931 // dm[1] - this is the dm[1] from the previous TRCE
932 dm[2]= 1.05;
933 dm[3]= 19.65;
934
935 gMC->Gsposp("TRCE",9,"TIC1",0.,-1.2,0.,0,"ONLY",dm,4);
936
937 //
938 // TPc Space for Connectors
939 //
940
941 dm[0]= .3;
942 dm[1]= .3;
943 dm[2]= 4.5;
944
945 gMC->Gsvolu("TPSC","BOX ",idtmed[0],dm,3);
946
947 // TPC Connectors
948
949 dm[0]= .25;
950 dm[1]= .15;
951 dm[2]= 3.75;
952
953 gMC->Gsvolu("TPCC","BOX ",idtmed[13],dm,3);
954
955 gMC->Gspos("TPCC",1,"TPSC",0.,0.15,0.,0,"ONLY");
956
957 zz = -12.7875;
958
959
960 Float_t alpha;
961 Float_t astep;
962
37831078 963 // inner part of the inner sector - 2 x 20 holes
964
965 astep = 20.00096874/19.;
966
967 alpha = 10.00048437-astep;
968
969 Float_t x1,x2;
970
971 x1 = 13.31175725;
972 x1 -= 0.996357832;
973
974 x2 = 15.06180253;
975 x2 -= 1.163028812;
976
977 Int_t ncon;
1283eee5 978
37831078 979 for(ncon=0;ncon<20;ncon++){
1283eee5 980
37831078 981 phi1 = 0.;
982 theta1 = 90.+alpha;
983 phi2=90.;
984 theta2 = 90.;
985 phi3 = (alpha>0) ? 0. : 180.;
986 theta3 = TMath::Abs(alpha);
1283eee5 987
37831078 988 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1283eee5 989
37831078 990
991
992 gMC->Gspos("TPSC",ncon+1,"TIRC",x1,0.3,-12.7875,idrotm[nRotMat],"ONLY");
993 gMC->Gspos("TPSC",ncon+21,"TIRC",x2,0.3,-2.8625,idrotm[nRotMat],"ONLY");
994
995
996 x1 -= 1.296357833;
997 x2 -= 1.463028812;
998
999 alpha -= astep;
1000 nRotMat++;
1001
1002 }
1003
1004 // outer part of the inner sector - 2 x 25 holes
1005
1006 astep = 20.00096874/24.;
1007 alpha = 10.00048437-astep;
1283eee5 1008
37831078 1009 x1 = 16.81184781;
1010 x1 -= 1.016295986;
1283eee5 1011
37831078 1012 x2 = 18.5618931;
1013 x2 -= 1.150914854;
1283eee5 1014
37831078 1015 for(ncon=0;ncon<25;ncon++){
1283eee5 1016
37831078 1017 phi1 = 0.;
1018 theta1 = 90.+alpha;
1019 phi2=90.;
1020 theta2 = 90.;
1021 phi3 = (alpha>0) ? 0. : 180.;
1022 theta3 = TMath::Abs(alpha);
1283eee5 1023
37831078 1024 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1283eee5 1025
37831078 1026
1027
1028 gMC->Gspos("TPSC",ncon+41,"TIRC",x1,0.3,7.0625,idrotm[nRotMat],"ONLY");
1029 gMC->Gspos("TPSC",ncon+66,"TIRC",x2,0.3,16.9875,idrotm[nRotMat],"ONLY");
1283eee5 1030
1283eee5 1031
37831078 1032 x1 -= 1.316295986;
1033 x2 -= 1.450914854;
1283eee5 1034
37831078 1035 alpha -= astep;
1036 nRotMat++;
1283eee5 1037
1038 }
1283eee5 1039
37831078 1040 //--------------------------------------------------------------------------
1041 // TPC Outer Readout Chambers
1042 // this is NOT a final design
1043 //--------------------------------------------------------------------------
1044
1045 dm[0]= 23.3875;
1046 dm[1]= 43.524;
1047 dm[2]= 1.5; //thickness
1048 dm[3]= 57.1;
1049
1050 gMC->Gsvolu("TORC","TRD1",idtmed[4],dm,4);
1051
1052 //------------------------------------------------
1053 // Tpc Outer readout chamber Pad Plane
1054 //------------------------------------------------
1055
1056 dm[2]= 0.5;
1057
1058 gMC->Gsvolu("TOPP","TRD1",idtmed[12],dm,4);
1059
1060 dm[0] -= 1.218511934;
1061 dm[1] -= 1.218511934;
1062 dm[2] = 0.35;
1063
1064 gMC->Gsvolu("TOC3","TRD1",idtmed[1],dm,4);
1065
1066 gMC->Gspos("TOC3",1,"TOPP",0.,0.15,0.,0,"ONLY");
1067
1068 gMC->Gspos("TOPP",1,"TORC",0.,1.0,0.,0,"ONLY");
1069
1070 // empty space
1071
1072 dm[0]= 21.035;
1073 dm[1]= 38.7205;
1074 dm[2]= 0.7;
1075 dm[3]= 50.15;
1076
1077 gMC->Gsposp("TRCE",10,"TORC",0.,-0.8,-2.15,0,"ONLY",dm,4);
1078
1079 dm[0]= 22.2935;
1080 dm[1]= 40.5085;
1081 dm[2]= 2.25;
1082 dm[3]= 51.65;
1083
1084 gMC->Gsvolu("TOC1","TRD1",idtmed[4],dm,4);
1085
1086 dm[0]= 21.35;
1087 dm[1]= 38.7205;
1088 dm[2]= 2.25;
1089 dm[3]= 50.15;
1090
79575d04 1091 gMC->Gsposp("TRCE",11,"TOC1",0.,0.,0.,0,"ONLY",dm,4);
37831078 1092
1093 //-----------------------------------------------
1094 // Tpc Services Support Wheel
1095 //-----------------------------------------------
1096
1097 dm[0]=0.;
1098 dm[1]=360.;
1099 dm[2]=18.;
1100 dm[3]=2.;
1101
1102 dm[4]= -5.;
1103 dm[5]= 77.017;
1104 dm[6]= 255.267;
1105
1106 dm[7]= 5.;
1107 dm[8]= dm[5];
1108 dm[9]= dm[6];
1109
1110 gMC->Gsvolu("TSSW","PGON",idtmed[4],dm,10);
1111
1112 // Tpc Services Wheel Cover
1113
1114 dm[4]= -0.5;
1115 dm[7]= 0.5;
1116
1117 gMC->Gsvolu("TSWC","PGON",idtmed[4],dm,10);
1118
1119 // Tpc Service wheel Cover Empty space
1120
1121 dm[0]= 10.99;
1122 dm[1]= 39.599;
1123 dm[2]= .5;
1124 dm[3]= 81.125;
1125
1126 gMC->Gsvolu("TSCE","TRD1",idtmed[0],dm,4);
1127
1128 // Tpc services Wheel Empty Spaces
1129
1130 dm[0]= 13.18017507;
1131 dm[1]= 44.61045938;
1132 dm[2]= 4.;
1133 dm[3]= 89.125;
1134
1135 gMC->Gsvolu("TWES","TRD1",idtmed[0],dm,4);
1136
1137 // Tpc Services Wheel Bars
1138
1139 gMC->Gsvolu("TSWB","TRD1",idtmed[4],dm,0);
1140
1141 // bars-> TWES
1142
1143 dm[2]= 4.;
1144 dm[3]= .4;
1145
1146 dm[0]= 13.8149522;
1147 dm[1]= 13.95601379;
fe4da5cc 1148
37831078 1149 gMC->Gsposp("TSWB",1,"TWES",0.,0.,-85.125,0,"ONLY",dm,4);
1150
1151 dm[0]= 43.83462067;
1152 dm[1]= 43.97568225;
1153
1154 gMC->Gsposp("TSWB",2,"TWES",0.,0.,85.125,0,"ONLY",dm,4);
1155
1156 // TPc ELectronics - right now 30% X0 Si
1157
1158 dm[0]= 14.03813696;
1159 dm[1]= 43.3524075;
1160 dm[2]= 1.404;
1161 dm[3]= 83.125;
1162
1163 gMC->Gsvolu("TPEL","TRD1",idtmed[11],dm,4);
1164 gMC->Gspos("TPEL",1,"TWES",0.,0.,0.,0,"ONLY");
1165
1166 //--------------------------------------------------------------------------
1167 // End caps
1168 //--------------------------------------------------------------------------
1169
1170 // TPc Main Wheel - Al
1171
79575d04 1172 dm[0]= 74.9;
1173 dm[1]= 264.4;
37831078 1174 dm[2]= 3.0;
1175
1176 gMC->Gsvolu("TPMW","TUBE",idtmed[4],dm,3);
1177
37831078 1178 //--------------------------------------------------------------------------
1179 // Tpc Empty Space for the Readout chambers
1180 //--------------------------------------------------------------------------
1181
1182 Float_t rLow= 86.2;
1183 Float_t rUp= 243.5;
1184 Float_t dR = 0.5*(rUp-rLow);
1185
79575d04 1186 space= 1.5/cosAlpha; // wheel ribs are 3.0 cm wide
37831078 1187
1188 dm[0]= rLow*tanAlpha-space;
1189 dm[1]= rUp*tanAlpha-space;
1190 dm[2] = 3.0;
1191 dm[3]= dR;
1192
1193 gMC->Gsvolu("TESR","TRD1",idtmed[0],dm,4);
1194
1195 // TIC1->TESR
1196
1197 gMC->Gspos("TIC1",1,"TESR",0.,0.75,-dR+23.97,0,"ONLY");
1198
1199 // TOC1->TESR
1200
1201 gMC->Gspos("TOC1",1,"TESR",0.,0.75,dR-55.02,0,"ONLY");
1202
1203 // Tpc Empty Space Bars - Al (daughters of TESR)
1204
1205 Float_t zBar;
1206
1207 gMC->Gsvolu("TESB","TRD1",idtmed[4],dm,0);
1208
1209 // lower bar
1210
1211 dm[0]= rLow*tanAlpha-space;
1212 dm[1]= 88.7*tanAlpha-space;
79575d04 1213 dm[2]= 2.25;
37831078 1214 dm[3]= 1.275;
1215
1216 zBar = -dR+dm[3];
1217
79575d04 1218 gMC->Gsposp("TESB",1,"TESR",0.,0.75,zBar,0,"ONLY",dm,4);
37831078 1219
1220 // middle bar
1221
1222 dm[0]= 131.65*tanAlpha-space;
1223 dm[1]= 136.5*tanAlpha-space;
1224 dm[3]= 2.425;
1225
1226 zBar = -dR +131.65+dm[3]-rLow;
1227
79575d04 1228 gMC->Gsposp("TESB",2,"TESR",0.,0.75,zBar,0,"ONLY",dm,4);
37831078 1229
1230 // upper bar
1231
1232 dm[0]= 240.4*tanAlpha-space;
1233 dm[1]= rUp*tanAlpha-space;
1234 dm[3]= 1.55;
1235
1236 zBar = dR-dm[3];
1237
79575d04 1238 gMC->Gsposp("TESB",3,"TESR",0.,0.75,zBar,0,"ONLY",dm,4);
37831078 1239
1240 // positioning of the empty spaces into the main wheel
1241
37831078 1242 Float_t rCenter,xc,yc;
1243 Float_t rInner,rOuter; // center of the inner and outer chamber
1244
1245 rCenter = rLow+dR;
1246
1247 rInner = 108.07;
1248 rOuter = 190.68;
1249
1250
1251 for(Int_t ns=0; ns<nInnerSector;ns++){
1252
1253 phi1 = ns * innerOpenAngle + innerAngleShift;
1254 phi1 *= kRaddeg; // in degrees
1255
1256 phi1 = (Float_t)TMath::Nint(phi1) + 270.;
1257
1258 if (phi1 > 360.) phi1 -= 360.;
1259
1260 theta1 = 90.;
1261 phi2 = 90.;
1262 theta2 = 180.;
1263 phi3 = ns * innerOpenAngle + innerAngleShift;
1264 phi3 *= kRaddeg; // in degrees
1265
1266 phi3 = (Float_t)TMath::Nint(phi3);
1267
1268 if(phi3 > 360.) phi3 -= 360.;
1269
1270 theta3 = 90.;
1271
1272 // "holes"->End plate
1273
1274 xc = rCenter*TMath::Cos(phi3*kDegrad);
1275 yc = rCenter*TMath::Sin(phi3*kDegrad);
1276
1277 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1278
1279 gMC->Gspos("TESR",ns+1,"TPMW",xc,yc,0.,idrotm[nRotMat],"ONLY");
1280
1281 // TSCE->TSWC (services wheel volumes)
1282
1283 xc = 166.142*TMath::Cos(phi3*kDegrad);
1284 yc = 166.142*TMath::Sin(phi3*kDegrad);
1285
1286 gMC->Gspos("TSCE",ns+1,"TSWC",xc,yc,0.,idrotm[nRotMat],"ONLY");
1287 gMC->Gspos("TWES",ns+1,"TSSW",xc,yc,0.,idrotm[nRotMat],"ONLY");
1288
1289
1290 // readout chambers->TDGN (drift gas)
1291
1292 xc = rInner*TMath::Cos(phi3*kDegrad);
1293 yc = rInner*TMath::Sin(phi3*kDegrad);
1294
1295 gMC->Gspos("TIRC",ns+1,"TDGN",xc,yc,252.,idrotm[nRotMat],"ONLY");
1296
1297 xc = rOuter*TMath::Cos(phi3*kDegrad);
1298 yc = rOuter*TMath::Sin(phi3*kDegrad);
1299
1300 gMC->Gspos("TORC",ns+1,"TDGN",xc,yc,252.1,idrotm[nRotMat],"ONLY");
1301
1302 nRotMat++;
1303
1304 theta2 = 0.; // reflection
1305
1306 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1307
1308 xc = rInner*TMath::Cos(phi3*kDegrad);
1309 yc = rInner*TMath::Sin(phi3*kDegrad);
1310
1311 gMC->Gspos("TIRC",ns+nInnerSector+1,"TDGN",xc,yc,-252.,idrotm[nRotMat],"ONLY");
1312
1313 xc = rOuter*TMath::Cos(phi3*kDegrad);
1314 yc = rOuter*TMath::Sin(phi3*kDegrad);
1315
1316 gMC->Gspos("TORC",ns+nOuterSector+1,"TDGN",xc,yc,-252.1,idrotm[nRotMat],"ONLY");
1317
1318 nRotMat++;
1319
1320 }
37831078 1321 // TPMW->TPC
1322
1323 gMC->Gspos("TPMW",1,"TPC ",0.,0.,256.6,0,"ONLY");
79575d04 1324 gMC->Gspos("TPMW",2,"TPC ",0.,0.,-256.6,idrotm[0],"ONLY");
37831078 1325
79575d04 1326 //---------------------------------------------------------
1327 // Tpc High Voltage Membrane - 100 microns of mylar
1328 //---------------------------------------------------------
37831078 1329
79575d04 1330 dm[0]=82.8;
1331 dm[1]=252.;
1332 dm[2]=0.005;
37831078 1333
79575d04 1334 gMC->Gsvolu("THVM","TUBE",idtmed[8],dm,3);
1283eee5 1335
79575d04 1336 gMC->Gspos("THVM",1,"TDGN",0.,0.,0.,0,"ONLY");
1283eee5 1337
79575d04 1338 // Tpc High Voltage membrane Holders
37831078 1339
79575d04 1340 gMC->Gsvolu("THVH","TUBE",idtmed[4],dm,0);
37831078 1341
79575d04 1342
37831078 1343
79575d04 1344 // inner
37831078 1345
79575d04 1346 dm[0]=79.3;
1347 dm[1]=82.8;
1348 dm[2]=0.2;
37831078 1349
79575d04 1350 gMC->Gsposp("THVH",1,"TDGN",0.,0.,0.,0,"ONLY",dm,3);
1351
1352 // outer
1353
1354 dm[0]= 252.;
1355 dm[1]= 257.9;
1356 dm[2]= 0.4;
1357
1358 gMC->Gsposp("THVH",2,"TDGN",0.,0.,0.,0,"ONLY",dm,3);
37831078 1359
1360 //----------------------------------------------------------
1361 // TPc Support Rods - MAKROLON
1362 //----------------------------------------------------------
1363
1364 dm[0]= 0.9;
1365 dm[1]= 1.2;
37831078 1366
79575d04 1367 gMC->Gsvolu("TPSR","TUBE",idtmed[7],dm,0); // inner and outer rods differ
1368
37831078 1369
79575d04 1370 for(Int_t nrod=0;nrod<18;nrod++){
37831078 1371 Float_t angle=innerOpenAngle*(Float_t)nrod;
1372
79575d04 1373 xc=81.5*TMath::Cos(angle);
1374 yc=81.5*TMath::Sin(angle);
1283eee5 1375
79575d04 1376 dm[2]=126.7;
1283eee5 1377
79575d04 1378 gMC->Gsposp("TPSR",nrod+1,"TDGN",xc,yc,126.9,0,"ONLY",dm,3);
1379 gMC->Gsposp("TPSR",nrod+19,"TDGN",xc,yc,-126.9,0,"ONLY",dm,3);
1283eee5 1380
79575d04 1381 dm[2]=126.6;
1382
1383 xc=254.25*TMath::Cos(angle);
1384 yc=254.25*TMath::Sin(angle);
1385
1386 // rod number 54 contans the HV cable
1283eee5 1387
79575d04 1388 if(nrod<17) {
1389 gMC->Gsposp("TPSR",nrod+37,"TDGN",xc,yc,127.,0,"ONLY",dm,3);
1390 gMC->Gsposp("TPSR",nrod+54,"TDGN",xc,yc,-127.,0,"ONLY",dm,3);
1391 }
1392
37831078 1393 }
1394
1395 //----------------------------------------------------------
79575d04 1396 // Tpc High Voltage Rod - MAKROLON + Copper cable
37831078 1397 //----------------------------------------------------------
1283eee5 1398
37831078 1399 // rod with cable (Left)
1400
1401 dm[0]=0.;
1402 dm[1]=2.25;
79575d04 1403 dm[2]=126.6;
37831078 1404
1405 gMC->Gsvolu("THVL","TUBE",idtmed[7],dm,3);
1406
1407 // HV cable
1283eee5 1408
37831078 1409 dm[0]=0.;
1410 dm[1]=0.3;
79575d04 1411 dm[2]=126.6;
1283eee5 1412
79575d04 1413 gMC->Gsvolu("THVC","TUBE",idtmed[10],dm,3);
1283eee5 1414
37831078 1415 // empty space
1283eee5 1416
37831078 1417 dm[0]=0.3;
1418 dm[1]=1.;
79575d04 1419 dm[2]=126.6;
1283eee5 1420
37831078 1421 gMC->Gsvolu("THVE","TUBE",idtmed[1],dm,3);
1283eee5 1422
37831078 1423 gMC->Gspos("THVC",1,"THVL",0.,0.,0.,0,"ONLY");
1424 gMC->Gspos("THVE",1,"THVL",0.,0.,0.,0,"ONLY");
1283eee5 1425
37831078 1426 // rod without cable
1283eee5 1427
37831078 1428 dm[0]=1.8;
1429 dm[1]=2.25;
79575d04 1430 dm[2]=126.6;
1283eee5 1431
37831078 1432 gMC->Gsvolu("THVR","TUBE",idtmed[7],dm,3);
1283eee5 1433
79575d04 1434 gMC->Gspos("THVL",1,"TDGN",xc,yc,-127.,0,"ONLY");
1435 gMC->Gspos("THVR",1,"TDGN",xc,yc,127.,0,"ONLY");
37831078 1436
1437 gMC->Gspos("TDGN",1,"TPC ",0.,0.,0.,0,"ONLY");
79575d04 1438
37831078 1439 // services wheel cover -> wheel
1440
1441
1442 gMC->Gspos("TSWC",1,"TSSW",0.,0.,4.5,0,"ONLY");
1443 gMC->Gspos("TSWC",2,"TSSW",0.,0.,-4.5,0,"ONLY");
1444
1445
1446 // put the wheel into the TPC
1447
1448 gMC->Gspos("TSSW",1,"TPC ",0.,0.,278.7,0,"ONLY");
1449 gMC->Gspos("TSSW",2,"TPC ",0.,0.,-278.7,0,"ONLY");
1450
1451 //
1452
1453 gMC->Gsord("TPMW",6);
1454 gMC->Gsord("TSSW",6);
1455 gMC->Gsord("TSWC",6);
1456
1457 // put the TPC into ALIC (main mother volume)
1458
1459 gMC->Gspos("TPC ",1,"ALIC",0.,0.,0.,0,"ONLY");
79575d04 1460
1283eee5 1461} // end of function
1462
fe4da5cc 1463
1464//_____________________________________________________________________________
1465void AliTPCv0::CreateMaterials()
1466{
1467 //
1468 // Define materials for the TPC
1469 //
1470 AliTPC::CreateMaterials();
1471}
1472
1473//_____________________________________________________________________________
8c555625 1474void AliTPCv0::DrawDetector()
fe4da5cc 1475{
1476 //
1477 // Draw a shaded view of the Time Projection Chamber version 0
1478 //
1479
fe4da5cc 1480 // Set everything unseen
cfce8870 1481 gMC->Gsatt("*", "seen", -1);
fe4da5cc 1482 //
1483 // Set ALIC mother transparent
cfce8870 1484 gMC->Gsatt("ALIC","SEEN",0);
fe4da5cc 1485 //
1486 // Set the volumes visible
37831078 1487 gMC->Gsatt("TPC ","SEEN",0);
1488 gMC->Gsatt("TOIN","SEEN",1);
1489 gMC->Gsatt("TOIN","COLO",7);
37831078 1490 gMC->Gsatt("TOCV","SEEN",1);
1491 gMC->Gsatt("TOCV","COLO",4);
1492 gMC->Gsatt("TSA1","SEEN",0);
1493 gMC->Gsatt("TSA2","SEEN",0);
1494 gMC->Gsatt("TSA3","SEEN",0);
79575d04 1495 gMC->Gsatt("TSA4","SEEN",0);
1496 gMC->Gsatt("TSA5","SEEN",0);
37831078 1497 gMC->Gsatt("TOFC","SEEN",1);
1498 gMC->Gsatt("TOFC","COLO",4);
37831078 1499 gMC->Gsatt("TSA6","SEEN",0);
1500 gMC->Gsatt("TSA7","SEEN",0);
79575d04 1501 gMC->Gsatt("TSA8","SEEN",0);
37831078 1502 gMC->Gsatt("TIIN","SEEN",1);
79575d04 1503 gMC->Gsatt("TIIN","COLO",7);
1504 gMC->Gsatt("TII1","SEEN",0);
1505 gMC->Gsatt("TIFC","SEEN",1);
1506 gMC->Gsatt("TIFC","COLO",4);
1507 gMC->Gsatt("TSA9","SEEN",0);
37831078 1508 gMC->Gsatt("TS10","SEEN",0);
1509 gMC->Gsatt("TS11","SEEN",0);
1510 gMC->Gsatt("TS12","SEEN",0);
37831078 1511 gMC->Gsatt("TS13","SEEN",0);
1512 gMC->Gsatt("TS14","SEEN",0);
79575d04 1513 gMC->Gsatt("TICC","SEEN",0);
1514 gMC->Gsatt("TICM","SEEN",0);
37831078 1515 gMC->Gsatt("TS15","SEEN",0);
1516 gMC->Gsatt("TS16","SEEN",0);
37831078 1517 gMC->Gsatt("TS17","SEEN",0);
79575d04 1518 gMC->Gsatt("TS18","SEEN",0);
1519 gMC->Gsatt("TS19","SEEN",0);
1520 gMC->Gsatt("TPJ1","SEEN",0);
1521 gMC->Gsatt("TPJ2","SEEN",0);
1522 gMC->Gsatt("TICS","SEEN",0);
1523 gMC->Gsatt("TDGN","SEEN",0);
37831078 1524 gMC->Gsatt("TIRC","SEEN",0);
1525 gMC->Gsatt("TIC1","SEEN",1);
1526 gMC->Gsatt("TIPP","SEEN",0);
1527 gMC->Gsatt("TIC3","SEEN",0);
1528 gMC->Gsatt("TRCE","SEEN",0);
1529 gMC->Gsatt("TPSC","SEEN",0);
79575d04 1530 gMC->Gsatt("TPCC","SEEN",0);
37831078 1531 gMC->Gsatt("TORC","SEEN",0);
1532 gMC->Gsatt("TOPP","SEEN",0);
1533 gMC->Gsatt("TOC3","SEEN",0);
1534 gMC->Gsatt("TOC1","SEEN",1);
1535 gMC->Gsatt("TSSW","SEEN",1);
1536 gMC->Gsatt("TSWC","SEEN",1);
37831078 1537 gMC->Gsatt("TSSW","COLO",3);
1538 gMC->Gsatt("TSWC","COLO",3);
1539 gMC->Gsatt("TSCE","COLO",6);
79575d04 1540 gMC->Gsatt("TSCE","SEEN",1);
37831078 1541 gMC->Gsatt("TWES","SEEN",0);
1542 gMC->Gsatt("TSWB","SEEN",0);
1543 gMC->Gsatt("TPEL","SEEN",0);
1544 gMC->Gsatt("TPMW","SEEN",1);
37831078 1545 gMC->Gsatt("TESR","SEEN",1);
1546 gMC->Gsatt("TPMW","COLO",12);
37831078 1547 gMC->Gsatt("TIC1","COLO",5);
79575d04 1548 gMC->Gsatt("TOC1","COLO",5);
37831078 1549 gMC->Gsatt("TESB","SEEN",0);
1550 gMC->Gsatt("THVM","SEEN",1);
1551 gMC->Gsatt("THVM","COLO",11);
79575d04 1552 gMC->Gsatt("THVH","SEEN",0);
1553 gMC->Gsatt("TPSR","SEEN",0);
37831078 1554 gMC->Gsatt("THVL","SEEN",0);
79575d04 1555 gMC->Gsatt("THVC","SEEN",0);
37831078 1556 gMC->Gsatt("THVE","SEEN",0);
1557 gMC->Gsatt("THVR","SEEN",0);
79575d04 1558
fe4da5cc 1559 //
cfce8870 1560 gMC->Gdopt("hide", "on");
1561 gMC->Gdopt("shad", "on");
1562 gMC->Gsatt("*", "fill", 7);
1563 gMC->SetClipBox(".");
37831078 1564 gMC->SetClipBox("TPMW",-300,300,-300,300,254.,270.);
1565 gMC->SetClipBox("TESR",-300,300,-300,300,254.,270.);
1566 gMC->SetClipBox("TSSW",-300,300,-300,300,283.,284.);
1567 gMC->SetClipBox("TSWC",-300,300,-300,300,283.,284.);
1568 gMC->SetClipBox("*", 0, 300, -300, 300, -290, 290);
cfce8870 1569 gMC->DefaultRange();
1570 gMC->Gdraw("alic", 40, 30, 0, 12, 9.5, .025, .025);
1571 gMC->Gdhead(1111, "Time Projection Chamber");
1572 gMC->Gdman(18, 4, "MAN");
1573 gMC->Gdopt("hide","off");
fe4da5cc 1574}
1575
1576//_____________________________________________________________________________
1577void AliTPCv0::Init()
1578{
1579 //
1580 // Initialise Time Projection Chamber version 0
1581 //
13a8c8d7 1582
1583 printf("%s: *** TPC version 0 initialized***\n",ClassName());
1584
1585 // printf("TPC version 0 initialized\n");
fe4da5cc 1586}
1587
1588//_____________________________________________________________________________
1589void AliTPCv0::StepManager()
1590{
1591 //
1592 // Procedure called at each step in the TPC
1593 //
1594}