]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TPC/AliTPCv1.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Add ResetDecayTable() and SsetDecayTable() methods.
[u/mrichter/AliRoot.git] / TPC / AliTPCv1.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
16/*
17$Log$
7ab50f40 18Revision 1.24 2000/11/30 11:48:50 kowal2
19TLorentzVector.h adde to satisfy the latest changes by Federico
20
a383b30d 21Revision 1.23 2000/11/14 10:48:57 kowal2
22Correct material used for TSA4. Thanks to J. Barbosa.
23
2c668420 24Revision 1.22 2000/11/06 17:24:10 kowal2
25Corrected bug in the outer containment vessel and
26the outer field cage geometry.
27Thanks to J. Barbosa.
28
720b23b8 29Revision 1.21 2000/11/02 16:55:24 kowal2
30Corrected bug in the inner containment vessel geometry.
31Thanks to J. Belikov
32
08db548a 33Revision 1.20 2000/11/02 07:24:11 kowal2
34Correction in the TPC geometry.
35Changes due to the new hit structure.
36
6a905708 37Revision 1.19 2000/07/10 20:57:39 hristov
38Update of TPC code and macros by M.Kowalski
39
37831078 40Revision 1.18 2000/06/30 12:07:50 kowal2
41Updated from the TPC-PreRelease branch
42
73042f01 43Revision 1.17.2.4 2000/06/26 07:39:42 kowal2
44Changes to obey the coding rules
45
46Revision 1.17.2.3 2000/06/25 08:38:41 kowal2
47Splitted from AliTPCtracking
48
49Revision 1.17.2.2 2000/06/16 12:58:13 kowal2
50Changed parameter settings
51
52Revision 1.17.2.1 2000/06/09 07:15:07 kowal2
53
54Defaults loaded automatically (hard-wired)
55Optional parameters can be set via macro called in the constructor
56
57Revision 1.17 2000/05/15 10:00:30 kowal2
58Corrected bug in the TPC geometry, thanks to Ivana Hrivnacova
59
e94dd9ca 60Revision 1.16 2000/04/17 09:37:33 kowal2
61removed obsolete AliTPCDigitsDisplay.C
62
cc80f89e 63Revision 1.15.8.2 2000/04/10 08:29:48 kowal2
64
65Different pad geometry for different sectors
66Updated readouch chambers geometry
67
68Revision 1.15.8.1 2000/04/10 07:56:53 kowal2
69Not used anymore - removed
70
71Revision 1.15 1999/11/04 17:28:07 fca
72Correct barrel part of HV Degrader
73
a5371c42 74Revision 1.14 1999/10/08 06:27:23 fca
75Corrected bug in the HV degrader geometry, thanks to G.Tabary
76
d183a600 77Revision 1.13 1999/10/04 13:39:54 fca
78Correct array index problem
79
cdde5b1e 80Revision 1.12 1999/09/29 09:24:34 fca
81Introduction of the Copyright and cvs Log
82
4c039060 83*/
84
73042f01 85//
fe4da5cc 86///////////////////////////////////////////////////////////////////////////////
87// //
88// Time Projection Chamber version 1 -- detailed TPC and fast simulation //
37831078 89// //
fe4da5cc 90//Begin_Html
91/*
1439f98e 92<img src="picts/AliTPCv1Class.gif">
fe4da5cc 93*/
94//End_Html
95// //
96// //
97///////////////////////////////////////////////////////////////////////////////
98
99#include <TMath.h>
100#include <TGeometry.h>
101#include <TNode.h>
102#include <TBRIK.h>
103#include <TTUBE.h>
104#include "AliTPCv1.h"
105#include "AliRun.h"
106#include <iostream.h>
107#include <fstream.h>
a383b30d 108#include "TLorentzVector.h"
fe4da5cc 109
110#include "AliMC.h"
111#include "AliConst.h"
112
8c555625 113#include "AliTPCParam.h"
73042f01 114#include "AliTPCParamSR.h"
cc80f89e 115#include "AliTPCDigitsArray.h"
73042f01 116#include <TInterpreter.h>
8c555625 117
fe4da5cc 118ClassImp(AliTPCv1)
6a905708 119
120//_____________________________________________________________________________
121
122 AliTPCv1::AliTPCv1(){
123
124 fHitType = 1;
125
126}
fe4da5cc 127
128//_____________________________________________________________________________
129AliTPCv1::AliTPCv1(const char *name, const char *title)
130 :AliTPC(name, title) {
131 //
132 // Standard constructor for Time Projection Chamber
133 //
37831078 134 fIdSens=fIdLSec=fIdUSec=0;
73042f01 135
136 SetGasMixt(2,20,10,-1,0.9,0.1,0.);
137
138 if (fTPCParam)
139 fTPCParam->Write(fTPCParam->GetTitle());
6a905708 140
141 fHitType=1;
fe4da5cc 142}
143
144//_____________________________________________________________________________
145void AliTPCv1::CreateGeometry()
146{
37831078 147 //-------------------------------------------------------------
fe4da5cc 148 // Creates geometry for Time Projection Chamber version 1
37831078 149 // Detailed geometry -- Fast simulation - space points
150 // at the pad-rows centers
151 // Origin M.Kowalski, INP Cracow, Marek.Kowalski@ifj.edu.pl
152 //-----------------------------------------------------------
153
fe4da5cc 154 //Begin_Html
155 /*
37831078 156 <img src="picts/AliTPC.gif">
fe4da5cc 157 */
158 //End_Html
159 //Begin_Html
160 /*
1439f98e 161 <img src="picts/AliTPCv1Tree.gif">
fe4da5cc 162 */
163 //End_Html
37831078 164 Float_t dm[50];
fe4da5cc 165
37831078 166 Int_t *idtmed = fIdtmed->GetArray(); // TPC media
8c555625 167
37831078 168 Int_t idrotm[120]; // rotation matrices
8c555625 169
37831078 170 Int_t nRotMat = 0; // actual number of rotation matrices
1283eee5 171
1283eee5 172
37831078 173 //
174 // Mother volume (Air) - all volumes will be positioned in it
175 //
176
177 dm[0]=0.;
178 dm[1]=360.;
179 dm[2]=12.;
1283eee5 180
37831078 181 //
182
183 dm[3]= -283.7;
184 dm[4]= 66.2;
185 dm[5]= 277.95;
1283eee5 186
37831078 187 //
1283eee5 188
37831078 189 dm[6]= -255.6;
190 dm[7]= 66.2;
191 dm[8]= 277.95;
1283eee5 192
37831078 193 //
1283eee5 194
37831078 195 dm[9]= -73.3;
196 dm[10]= 59.0;
197 dm[11]= 277.95;
1283eee5 198
37831078 199 //
1283eee5 200
37831078 201 dm[12]= -73.3;
202 dm[13]= 56.9;
203 dm[14]= 277.95;
1283eee5 204
37831078 205 //
fe4da5cc 206
37831078 207 dm[15]= -72.1;
208 dm[16]= 56.9;
209 dm[17]= 277.95;
1283eee5 210
37831078 211 //
1283eee5 212
37831078 213 dm[18]= -72.1;
214 dm[19]= 60.65;
215 dm[20]= 277.95;
1283eee5 216
37831078 217 //
1283eee5 218
37831078 219 dm[21]= 72.1;
220 dm[22]= 60.65;
221 dm[23]= 277.95;
1283eee5 222
37831078 223 //
1283eee5 224
37831078 225 dm[24]= 72.1;
226 dm[25]= 56.9;
227 dm[26]= 277.95;
1283eee5 228
37831078 229 //
1283eee5 230
37831078 231 dm[27]= 73.3;
232 dm[28]= 56.9;
233 dm[29]= 277.95;
1283eee5 234
37831078 235 //
1283eee5 236
37831078 237 dm[30]= 73.3;
7ab50f40 238 dm[31]= 59.;
37831078 239 dm[32]= 277.95;
1283eee5 240
37831078 241 //
1283eee5 242
37831078 243 dm[33]= 250.4;
244 dm[34]= 66.0;
245 dm[35]= 277.95;
1283eee5 246
37831078 247 //
1283eee5 248
37831078 249 dm[36]= 283.7;
250 dm[37]= 66.0;
251 dm[38]= 277.95;
1283eee5 252
1283eee5 253
37831078 254 gMC->Gsvolu("TPC ","PCON",idtmed[0],dm,39);
1283eee5 255
1283eee5 256
37831078 257 //-------------------------------------------------------------------
258 // Tpc Outer INsulator (CO2)
259 //-------------------------------------------------------------------
1283eee5 260
37831078 261 dm[0]= 0.;
262 dm[1]= 360.;
263 dm[2]= 6.;
1283eee5 264
37831078 265 //
266
267 dm[3]= -253.6;
268 dm[4]= 258.;
269 dm[5]= 266.65;
1283eee5 270
37831078 271 //
1283eee5 272
37831078 273 dm[6]= -253.;
274 dm[7]= 258.;
275 dm[8]= 266.65;
1283eee5 276
37831078 277 dm[9]= -253.;
278 dm[10]= 258.;
279 dm[11]= 277.97;
1283eee5 280
37831078 281 dm[12]= 253.6;
282 dm[13]= 258.;
283 dm[14]= 277.97;
1283eee5 284
37831078 285 //
1283eee5 286
37831078 287 dm[15]= 253.6;
288 dm[16]= 265.2;
289 dm[17]= 277.95;
cc80f89e 290
37831078 291 //
1283eee5 292
37831078 293 dm[18]= 255.6;
294 dm[19]= 265.2;
295 dm[20]= 277.95;
1283eee5 296
1283eee5 297
37831078 298 gMC->Gsvolu("TOIN","PCON",idtmed[3],dm,21);
1283eee5 299
37831078 300 //---------------------------------------------------------------
301 // shreds (G10) - TPC Rings
302 //---------------------------------------------------------------
cc80f89e 303
37831078 304 gMC->Gsvolu("TPCR","TUBE",idtmed[12],dm,0);
1283eee5 305
37831078 306 dm[0]= 258.;
307 dm[1]= 266.65;
308 dm[2]= 0.3;
1283eee5 309
37831078 310 gMC->Gsposp("TPCR",1,"TOIN",0.,0.,-253.3,0,"ONLY",dm,3); // left bottom
1283eee5 311
37831078 312 //
1283eee5 313
37831078 314 dm[0]= 258.;
315 dm[1]= 270.9;
316 dm[2]= 0.3;
1283eee5 317
37831078 318 gMC->Gsposp("TPCR",2,"TOIN",0.,0.,253.3,0,"ONLY",dm,3); // right
1283eee5 319
37831078 320 //
1283eee5 321
37831078 322 dm[0]= 272.2;
323 dm[1]= 277.95;
324 dm[2]= 0.3;
1283eee5 325
37831078 326 gMC->Gsposp("TPCR",3,"TOIN",0.,0.,-250.7,0,"ONLY",dm,3); // left top
1283eee5 327
37831078 328 //----------------------------------------------------------------
329 // Tpc Outer Contaiment Vessel
330 // mother volume - Al, daughters - composite (sandwich)
331 //----------------------------------------------------------------
332
333 dm[0]= 0.;
334 dm[1]= 360.;
720b23b8 335 dm[2]=6.;
1283eee5 336
37831078 337 //
1283eee5 338
37831078 339 dm[3]= -250.4;
340 dm[4]= 272.2;
341 dm[5]= 277.95;
1283eee5 342
37831078 343 //
8c555625 344
37831078 345 dm[6]= -248.4;
346 dm[7]= 272.2;
347 dm[8]= 277.95;
8c555625 348
37831078 349 //
720b23b8 350 dm[9]= -248.4;
351 dm[10]= 274.81;
352 dm[11]= 277.95;
1283eee5 353
37831078 354 //
1283eee5 355
37831078 356 dm[12]= 253.6;
720b23b8 357 dm[13]= 274.81;
37831078 358 dm[14]= 277.95;
1283eee5 359
720b23b8 360 //
8c555625 361
720b23b8 362 dm[15]= 253.6;
37831078 363 dm[16]= 265.2;
364 dm[17]= 277.95;
8c555625 365
720b23b8 366 //
367
368 dm[18]= 255.6;
369 dm[19]= 265.2;
370 dm[20]= 277.95;
371
372 gMC->Gsvolu("TOCV","PCON",idtmed[4],dm,21);
1283eee5 373
37831078 374 // Daughter volumes
1283eee5 375
37831078 376 // Tpc SAndwich 1 - Al
377
378 dm[0]= 274.81;
379 dm[1]= 277.95;
380 dm[2]= 251.7;
1283eee5 381
37831078 382 gMC->Gsvolu("TSA1","TUBE",idtmed[4],dm,3);
1283eee5 383
37831078 384 // Tpc SAndwich 2 - Tedlar
1283eee5 385
37831078 386 dm[0] += 5.e-3;
387 dm[1] -= 5.e-3;
388
389 gMC->Gsvolu("TSA2","TUBE",idtmed[9],dm,3);
1283eee5 390
37831078 391 // Tpc SAndwich 3 - Kevlar
1283eee5 392
37831078 393 dm[0] += 5e-3;
394 dm[1] -= 5.e-3;
1283eee5 395
37831078 396 gMC->Gsvolu("TSA3","TUBE",idtmed[5],dm,3);
1283eee5 397
37831078 398 // Tpc SAndwich 4 - NOMEX honeycomb
1283eee5 399
37831078 400 dm[0] += 0.06;
401 dm[1] -= 0.06;
1283eee5 402
2c668420 403 gMC->Gsvolu("TSA4","TUBE",idtmed[6],dm,3);
37831078 404
405 // 4->3->2->1->TOCV
1283eee5 406
37831078 407 gMC->Gspos("TSA4",1,"TSA3",0.,0.,0.,0,"ONLY");
408 gMC->Gspos("TSA3",1,"TSA2",0.,0.,0.,0,"ONLY");
409 gMC->Gspos("TSA2",1,"TSA1",0.,0.,0.,0,"ONLY");
1283eee5 410
37831078 411 gMC->Gspos("TSA1",1,"TOCV",0.,0.,2.6,0,"ONLY");
1283eee5 412
37831078 413 // TCOV-> TOIN
1283eee5 414
37831078 415 gMC->Gspos("TOCV",1,"TOIN",0.,0.,0.,0,"ONLY");
1283eee5 416
37831078 417 //-------------------------------------------------------
418 // Tpc Outer Field Cage
419 // mother volume - Al, daughters - composite (sandwich)
420 //-------------------------------------------------------
1283eee5 421
37831078 422 dm[0]=0.;
423 dm[1]=360.;
424 dm[2]=6.;
cc80f89e 425
37831078 426 dm[3]= -253.;
427 dm[4]= 258.;
428 dm[5]= 277.95;
1283eee5 429
37831078 430 //
1283eee5 431
37831078 432 dm[6]= -251.;
433 dm[7]= 258.;
434 dm[8]= 277.95;
8c555625 435
37831078 436 //
437
438 dm[9]= -251.;
439 dm[10]= 258.;
440 dm[11]= 260.05;
1283eee5 441
37831078 442 //
1283eee5 443
37831078 444 dm[12]= 251.;
445 dm[13]= 258.;
446 dm[14]= 260.05;
1283eee5 447
37831078 448 //
8c555625 449
37831078 450 dm[15]= 251.;
451 dm[16]= 258.;
720b23b8 452 dm[17]= 270.9;
1283eee5 453
37831078 454 //
1283eee5 455
37831078 456 dm[18]= 253.;
457 dm[19]= 258.;
458 dm[20]= 270.9;
1283eee5 459
37831078 460 gMC->Gsvolu("TOFC","PCON",idtmed[4],dm,21);
1283eee5 461
37831078 462 // Daughter volumes
1283eee5 463
6a905708 464 // Tpc SAndwich 5 - Tedlar
62a73ee5 465
37831078 466 dm[0]= 258.;
467 dm[1]= 260.05;
468 dm[2]= 251.7;
1283eee5 469
6a905708 470 gMC->Gsvolu("TSA5","TUBE",idtmed[9],dm,3);
1283eee5 471
6a905708 472 // Tpc SAndwich 6 - Kevlar
1283eee5 473
37831078 474 dm[0] += 5.e-3;
475 dm[1] -= 5.e-3;
1283eee5 476
6a905708 477 gMC->Gsvolu("TSA6","TUBE",idtmed[5],dm,3);
1283eee5 478
fe4da5cc 479
6a905708 480 // Tpc SAndwich 7 - NOMEX
1283eee5 481
37831078 482 dm[0] += 0.02;
483 dm[1] -= 0.02;
1283eee5 484
6a905708 485 gMC->Gsvolu("TSA7","TUBE",idtmed[6],dm,3);
1283eee5 486
6a905708 487 // 7->6->5->TOFC
1283eee5 488
37831078 489 gMC->Gspos("TSA7",1,"TSA6",0.,0.,0.,0,"ONLY");
490 gMC->Gspos("TSA6",1,"TSA5",0.,0.,0.,0,"ONLY");
1283eee5 491
37831078 492 gMC->Gspos("TSA5",1,"TOFC",0.,0.,0.,0,"ONLY");
1283eee5 493
37831078 494 // TOFC->TOIN
1283eee5 495
37831078 496 gMC->Gspos("TOFC",1,"TOIN",0.,0.,0.,0,"ONLY");
1283eee5 497
37831078 498 // TOIN->TPC
1283eee5 499
37831078 500 gMC->Gspos("TOIN",1,"TPC ",0.,0.,0.,0,"ONLY");
fe4da5cc 501
37831078 502 //--------------------------------------------------------------------
503 // Tpc Inner INsulator (CO2)
504 //--------------------------------------------------------------------
1283eee5 505
1283eee5 506
37831078 507 dm[0]=0.;
508 dm[1]= 360.;
509 dm[2]= 15.;
1283eee5 510
37831078 511 //
1283eee5 512
37831078 513 dm[3]= -255.6;
514 dm[4]= 66.2;
515 dm[5]= 74.8;
1283eee5 516
37831078 517 //
1283eee5 518
37831078 519 Float_t tanL = (66.2-59.0)/(255.6-73.3); // tangent of the left cone part
1283eee5 520
37831078 521 dm[6]= -253.6;
522 dm[7]= 59.0+ (253.6-73.3)*tanL;
523 dm[8]= 74.8;
1283eee5 524
37831078 525 //
1283eee5 526
37831078 527 dm[9]= -253.6;
528 dm[10]= dm[7];
529 dm[11]= 79.2;
1283eee5 530
37831078 531 //
1283eee5 532
37831078 533 dm[12]= -73.3;
534 dm[13]= 59.0;
535 dm[14]= 79.2;
1283eee5 536
37831078 537 //
1283eee5 538
37831078 539 dm[15]= -73.3;
540 dm[16]= 56.9;
541 dm[17]= 79.2;
1283eee5 542
37831078 543 //
1283eee5 544
37831078 545 dm[18]= -72.1;
546 dm[19]= 59.6;
547 dm[20]= 79.2;
1283eee5 548
549 //
37831078 550
551 dm[21]= -72.1;
552 dm[22]= 60.65;
553 dm[23]= 79.2;
554
1283eee5 555 //
556
37831078 557 dm[24]= 72.1;
558 dm[25]= 60.65;
559 dm[26]= 79.2;
560
561 //
1283eee5 562
37831078 563 dm[27]= 72.1;
564 dm[28]= 59.6;
565 dm[29]= 79.2;
1283eee5 566
37831078 567 //
568
569 dm[30]= 73.3;
570 dm[31]= 56.9;
571 dm[32]= 79.2;
572
573 //
1283eee5 574
37831078 575 dm[33]= 73.3;
576 dm[34]= 59.0;
577 dm[35]= 79.2;
1283eee5 578
37831078 579 //
580
581 dm[36]= 250.4;
582 dm[37]= 66.0;
583 dm[38]= 79.2;
584
585 //
fe4da5cc 586
37831078 587 dm[39]= 253.0;
588 dm[40]= 66.0;
589 dm[41]= 79.2;
1283eee5 590
37831078 591 //
592
593 dm[42]= 253.0;
594 dm[43]= 75.3;
595 dm[44]= 79.2;
596
597 //
1283eee5 598
37831078 599 dm[45]= 253.6;
600 dm[46]= 75.3;
601 dm[47]= 79.2;
1283eee5 602
37831078 603 gMC->Gsvolu("TIIN","PCON",idtmed[3],dm,48);
604
605
606 //--------------------------------------------------------------------
607 // Tpc Inner Containment vessel, Left part
608 // mother volume - Al, daughter - composite (sandwich)
609 //--------------------------------------------------------------------
610
611 dm[0]= 0.;
612 dm[1]= 360.;
613 dm[2]= 8.;
1283eee5 614
615 //
37831078 616
617 dm[3]= -255.6;
618 dm[4]= 66.2;
619 dm[5]= 74.8;
620
1283eee5 621 //
622
37831078 623 Float_t cosL = 1./TMath::Sqrt(1.+tanL*tanL);
624 Float_t sandThick = 2.14; // cone composite thickness
1283eee5 625
1283eee5 626
37831078 627 //
1283eee5 628
37831078 629 dm[6]= -253.6;
630 dm[7]= 59.0 + (253.6-73.3)*tanL;
631 dm[8]= 74.8;
1283eee5 632
37831078 633 //
1283eee5 634
37831078 635 dm[9]= -253.6;
636 dm[10]= dm[7];
637 dm[11]= dm[7]+sandThick/cosL;
1283eee5 638
37831078 639 //
1283eee5 640
37831078 641 dm[12]= -75.6;
642 dm[13]= 59.0+(75.6-73.3)*tanL;
643 dm[14]= dm[13]+sandThick/cosL;
1283eee5 644
37831078 645 //
1283eee5 646
37831078 647 dm[15]= -75.6;
648 dm[16]= dm[13];
649 dm[17]= 60.65;
1283eee5 650
37831078 651 //
1283eee5 652
37831078 653 dm[18]= -73.3;
654 dm[19]= 59.0;
655 dm[20]= 60.65;
1283eee5 656
37831078 657 //
1283eee5 658
37831078 659 dm[21]= -73.3;
660 dm[22]= 56.9;
661 dm[23]= 60.65;
1283eee5 662
37831078 663 //
1283eee5 664
37831078 665 dm[24]= -72.1;
666 dm[25]= 56.9;
667 dm[26]= 60.65;
1283eee5 668
37831078 669 gMC->Gsvolu("TICL","PCON",idtmed[4],dm,27);
1283eee5 670
37831078 671 // Daughter volumes
1283eee5 672
37831078 673 // Tpc SAndwich 9 - Al
1283eee5 674
37831078 675 dm[0]= 0.;
676 dm[1]= 360.;
677 dm[2]= 2.;
1283eee5 678
37831078 679 //
1283eee5 680
37831078 681 dm[3]= - 254.3;
682 dm[4]= 59.0+(254.3-73.3)*tanL;
683 dm[5]= dm[4]+sandThick/cosL;
1283eee5 684
37831078 685 //
1283eee5 686
37831078 687 dm[6]= -78.3;
688 dm[7]= 59.0+(78.3-73.3)*tanL;
689 dm[8]= dm[7]+sandThick/cosL;
1283eee5 690
37831078 691 //
1283eee5 692
37831078 693 gMC->Gsvolu("TSA9","PCON",idtmed[4],dm,9);
1283eee5 694
37831078 695 // Tpc SAndwich 10 - Tedlar
1283eee5 696
37831078 697 dm[4]+= 5.e-3/cosL;
698 dm[5]-= 5.e-3/cosL;
1283eee5 699
37831078 700 //
1283eee5 701
37831078 702 dm[7]+= 5.e-3/cosL;
703 dm[8]+= 5.e-3/cosL;
1283eee5 704
37831078 705 gMC->Gsvolu("TS10","PCON",idtmed[9],dm,9);
1283eee5 706
37831078 707 // Tpc SAndwich 11 - Kevlar
1283eee5 708
37831078 709 dm[4]+= 5.e-3/cosL;
710 dm[5]-= 5.e-3/cosL;
1283eee5 711
37831078 712 //
713
714 dm[7]+= 5.e-3/cosL;
715 dm[8]+= 5.e-3/cosL;
1283eee5 716
37831078 717 gMC->Gsvolu("TS11","PCON",idtmed[5],dm,9);
1283eee5 718
37831078 719 // Tpc SAndwich 12 - NOMEX
720
721 dm[4]+= 0.06/cosL;
722 dm[5]-= 0.06/cosL;
1283eee5 723
724 //
725
37831078 726 dm[7]+= 0.06/cosL;
727 dm[8]+= 0.06/cosL;
728
729 gMC->Gsvolu("TS12","PCON",idtmed[6],dm,9);
1283eee5 730
37831078 731 // 12->11->10->9
732
733 gMC->Gspos("TS12",1,"TS11",0.,0.,0.,0,"ONLY");
734 gMC->Gspos("TS11",1,"TS10",0.,0.,0.,0,"ONLY");
735 gMC->Gspos("TS10",1,"TSA9",0.,0.,0.,0,"ONLY");
736
737 // TSA9->TICL
738
739 gMC->Gspos("TSA9",1,"TICL",0.,0.,0.,0,"ONLY");
1283eee5 740
37831078 741 //--------------------------------------------------------------------
742 // Tpc Inner Containment vessel, Right part
743 // mother volume - Al, daughter - composite (sandwich)
744 //--------------------------------------------------------------------
745
746 dm[0]= 0.;
747 dm[1]= 360.;
748 dm[2]=8.;
1283eee5 749
37831078 750 //
1283eee5 751
37831078 752 dm[3]= 72.1;
753 dm[4]= 56.9;
754 dm[5]= 60.65;
1283eee5 755
37831078 756 //
1283eee5 757
37831078 758 dm[6]= 73.3;
759 dm[7]= 56.9;
760 dm[8]= 60.65;
1283eee5 761
37831078 762 //
1283eee5 763
37831078 764 dm[9]= 73.3;
765 dm[10]= 59.0;
766 dm[11]= 60.65;
1283eee5 767
37831078 768 //
1283eee5 769
37831078 770 Float_t tanR = (66.0-59.0)/(250.5-73.3); // to avoid accuracy problems
771 Float_t cosR = 1./TMath::Sqrt(1.+tanR*tanR); //as above
1283eee5 772
37831078 773 //
1283eee5 774
37831078 775 dm[12]= 75.6;
776 dm[13]= 59.0+(75.6-73.3)*tanR;
777 dm[14]= 60.65;
1283eee5 778
37831078 779 //
1283eee5 780
37831078 781 dm[15]= 75.6;
782 dm[16]= dm[13];
783 dm[17]= dm[16]+sandThick/cosR;
1283eee5 784
37831078 785 //
1283eee5 786
37831078 787 dm[18]= 248.4;
788 dm[19]= 59.0+(248.4-73.3)*tanR;
789 dm[20]= dm[19]+sandThick/cosR;
1283eee5 790
37831078 791 //
1283eee5 792
37831078 793 dm[21]= 248.4;
794 dm[22]= dm[19];
795 dm[23]= 70.2;
1283eee5 796
37831078 797 //
1283eee5 798
37831078 799 dm[24]= 250.4;
800 dm[25]= 66.0;
801 dm[26]= 70.2;
1283eee5 802
37831078 803 gMC->Gsvolu("TICR","PCON",idtmed[4],dm,27);
1283eee5 804
1283eee5 805
1283eee5 806
37831078 807 // Daughter volumes
1283eee5 808
37831078 809 // Tpc SAndwich 13 - Al
1283eee5 810
37831078 811 dm[0]= 0.;
812 dm[1]= 360.;
813 dm[2]= 2.;
1283eee5 814
37831078 815 //
1283eee5 816
37831078 817 dm[3]= 78.3;
818 dm[4]= 59.0+(78.3-73.3)*tanR;
819 dm[5]= dm[4]+sandThick/cosR;
1283eee5 820
37831078 821 //
1283eee5 822
37831078 823 dm[6]= 249.1;
824 dm[7]= 59.0+(249.1-73.3)*tanR;
825 dm[8]= dm[7]+sandThick/cosR;
1283eee5 826
37831078 827 //
1283eee5 828
37831078 829 gMC->Gsvolu("TS13","PCON",idtmed[4],dm,9);
1283eee5 830
37831078 831 // Tpc SAndwich 14 - Tedlar
fe4da5cc 832
37831078 833 dm[4]+= 5.e-3/cosR;
834 dm[5]-= 5.e-3/cosR;
fe4da5cc 835
37831078 836 //
1283eee5 837
37831078 838 dm[7]+= 5.e-3/cosR;
839 dm[8]+= 5.e-3/cosR;
1283eee5 840
37831078 841 gMC->Gsvolu("TS14","PCON",idtmed[9],dm,9);
1283eee5 842
37831078 843 // Tpc SAndwich 15 - Kevlar
1283eee5 844
37831078 845 dm[4]+= 5.e-3/cosR;
846 dm[5]-= 5.e-3/cosR;
1283eee5 847
37831078 848 //
fe4da5cc 849
37831078 850 dm[7]+= 5.e-3/cosR;
851 dm[8]+= 5.e-3/cosR;
1283eee5 852
37831078 853 gMC->Gsvolu("TS15","PCON",idtmed[5],dm,9);
1283eee5 854
37831078 855 // Tpc SAndwich 16 - NOMEX
1283eee5 856
37831078 857 dm[4]+= 0.06/cosR;
858 dm[5]-= 0.06/cosR;
1283eee5 859
37831078 860 //
1283eee5 861
37831078 862 dm[7]+= 0.06/cosR;
863 dm[8]+= 0.06/cosR;
1283eee5 864
37831078 865 gMC->Gsvolu("TS16","PCON",idtmed[6],dm,9);
1283eee5 866
37831078 867 // 16->15->14->13
1283eee5 868
37831078 869 gMC->Gspos("TS16",1,"TS15",0.,0.,0.,0,"ONLY");
870 gMC->Gspos("TS15",1,"TS14",0.,0.,0.,0,"ONLY");
871 gMC->Gspos("TS14",1,"TS13",0.,0.,0.,0,"ONLY");
1283eee5 872
37831078 873 // TS12->TICR
1283eee5 874
37831078 875 gMC->Gspos("TS13",1,"TICR",0.,0.,0.,0,"ONLY");
1283eee5 876
37831078 877 //------------------------------------------------------
878 // Tpc Inner Field Cage
879 // mother volume - Al, daughters - composite (sandwich)
880 //------------------------------------------------------
1283eee5 881
37831078 882 dm[0]= 0.;
883 dm[1]= 360.;
884 dm[2]=6.;
1283eee5 885
37831078 886 //
1283eee5 887
37831078 888 dm[3]= -253.0;
889 dm[4]= 70.7;
890 dm[5]= 79.2;
fe4da5cc 891
fe4da5cc 892 //
37831078 893
894 dm[6]= -251.0;
895 dm[7]= 70.7;
896 dm[8]= 79.2;
897
fe4da5cc 898 //
899
37831078 900 dm[9]= -251.0;
901 dm[10]= 77.15;
902 dm[11]= 79.2;
903
fe4da5cc 904 //
cc80f89e 905
37831078 906 dm[12]= 251.0;
907 dm[13]= 77.15;
908 dm[14]= 79.2;
fe4da5cc 909
fe4da5cc 910 //
37831078 911
912 dm[15]= 251.0;
913 dm[16]= 66.0;
914 dm[17]= 79.2;
915
fe4da5cc 916 //
fe4da5cc 917
37831078 918 dm[18]= 253.0;
919 dm[19]= 66.0;
920 dm[20]= 79.2;
921
922 gMC->Gsvolu("TIFC","PCON",idtmed[4],dm,21);
923
924 // Daughter volumes
925
6a905708 926 // Tpc Sandwich 17 - Tedlar
37831078 927
928 dm[0]= 77.15;
929 dm[1]= 79.2;
930 dm[2]= 251.7;
931
6a905708 932 gMC->Gsvolu("TS17","TUBE",idtmed[9],dm,3);
37831078 933
6a905708 934 // Tpc Sandwich 18 - Kevlar
37831078 935
936 dm[0]+= 5.e-3;
937 dm[1]-= 5.e-3;
938
6a905708 939 gMC->Gsvolu("TS18","TUBE",idtmed[5],dm,3);
37831078 940
37831078 941
6a905708 942 // Tpc Sandwich 19 - NOMEX
37831078 943
6a905708 944 dm[0]+= 0.02;
945 dm[1]-= 0.02;
37831078 946
6a905708 947 gMC->Gsvolu("TS19","TUBE",idtmed[6],dm,3);
37831078 948
6a905708 949 // 19->18->17
37831078 950
37831078 951 gMC->Gspos("TS19",1,"TS18",0.,0.,0.,0,"ONLY");
952 gMC->Gspos("TS18",1,"TS17",0.,0.,0.,0,"ONLY");
953
954 // TS17->TIFC
955
956 gMC->Gspos("TS17",1,"TIFC",0.,0.,0.,0,"ONLY");
957
958 // TPC Rings
959
960 dm[0]= 70.7;
961 dm[1]= 79.2;
962 dm[2]= 0.3;
963
964 gMC->Gsposp("TPCR",4,"TIIN",0.,0.,-253.3,0,"ONLY",dm,3);
965
966 dm[0]= 66.0;
967 dm[1]= 70.2;
968
969 gMC->Gsposp("TPCR",5,"TIIN",0.,0.,250.7,0,"ONLY",dm,3);
970
971 dm[0]= 75.3;
972 dm[1]= 79.2;
973
974 gMC->Gsposp("TPCR",6,"TIIN",0.,0.,253.3,0,"ONLY",dm,3);
975
976 // TICL->TIIN
977
978 gMC->Gspos("TICL",1,"TIIN",0.,0.,0.,0,"ONLY");
979
980 // TICR->TIIN
981
982 gMC->Gspos("TICR",1,"TIIN",0.,0.,0.,0,"ONLY");
983
984 // TIFC->TIIN
985
986 gMC->Gspos("TIFC",1,"TIIN",0.,0.,0.,0,"ONLY");
987
988 // Tpc Sandwich 21 - Al (central barrel)
989
08db548a 990 dm[0]= 60.65;
991 dm[1]= 61.21;
37831078 992 dm[2]= 75.2;
993
994 gMC->Gsvolu("TS21","TUBE",idtmed[4],dm,3);
995
996 // Tpc Sandwich 22 - Tedlar (central barrel)
997
998 dm[0]+= 5.e-3;
999 dm[1]-= 5.e-3;
1000
1001 gMC->Gsvolu("TS22","TUBE",idtmed[9],dm,3);
1002
1003 // Tpc Sandwich 23 - Kevlar (central barrel)
1004
1005 dm[0]+= 5.e-3;
1006 dm[1]-= 5.e-3;
1007
1008 gMC->Gsvolu("TS23","TUBE",idtmed[5],dm,3);
1009
1010 // Tpc Sandwich 24 - NOMEX (central barrel)
1011
1012 dm[0]+= 0.02;
1013 dm[1]-= 0.02;
1014
1015 gMC->Gsvolu("TS24","TUBE",idtmed[6],dm,3);
1016
1017 // 24->23->22->21
1018
1019 gMC->Gspos("TS24",1,"TS23",0.,0.,0.,0,"ONLY");
1020 gMC->Gspos("TS23",1,"TS22",0.,0.,0.,0,"ONLY");
1021 gMC->Gspos("TS22",1,"TS21",0.,0.,0.,0,"ONLY");
1022
1023 gMC->Gspos("TS21",1,"TIIN",0.,0.,0.,0,"ONLY");
1024
1025 // put everything into the TPC
1026
1027 gMC->Gspos("TIIN",1,"TPC ",0.,0.,0.,0,"ONLY");
1028
1029
1030 //---------------------------------------------------------
1031 // Tpc Dift Gas volume Nonsensitive (Ne-CO2 90/10)
1032 // and its daughters (HV membrane, rods, readout chambers)
1033 //---------------------------------------------------------
1034
1035 dm[0]= 79.2;
1036 dm[1]= 258.0;
1037 dm[2]= 253.6;
1038
1039 gMC->Gsvolu("TDGN","TUBE",idtmed[1],dm,3);
1040
1041 // sector opening angles
1042
1043 Float_t innerOpenAngle = fTPCParam->GetInnerAngle();
1044 Float_t outerOpenAngle = fTPCParam->GetOuterAngle();
1045
1046 // sector angle shift
1047
1048 Float_t innerAngleShift = fTPCParam->GetInnerAngleShift();
1049
1050 // number of sectors
1051
1052 Int_t nInnerSector = fTPCParam->GetNInnerSector()/2;
1053 Int_t nOuterSector = fTPCParam->GetNOuterSector()/2;
1054
1055 // All above parameters are identical for inner and outer
1056 // sectors. The distinction is kept for the historical reasons
1057 // and eventually will disappear.
1058
1059 Float_t tanAlpha = TMath::Tan(0.5*innerOpenAngle);
1060 Float_t cosAlpha = TMath::Sqrt(1.+tanAlpha*tanAlpha);
1061 Float_t space;
1062
1063 //-------------------------------------------------------------------------
1064 // Tpc Inner Readout Chambers
1065 //-------------------------------------------------------------------------
1066
1067 dm[0]= 14.483;
1068 dm[1]= 23.3345;
1069 dm[2]= 1.6; // thickness
1070 dm[3]= 25.1;
1071
1072 gMC->Gsvolu("TIRC","TRD1",idtmed[4],dm,4);
1073
1074 // this volume will be positioned in the empty space
1075 // of the end-cap to avoid overlaps
1076
1077 dm[0]= 13.7305;
1078 dm[1]= 21.1895;
1079 dm[2]= 2.25;
1080 dm[3]= 21.15;
1081
1082 gMC->Gsvolu("TIC1","TRD1",idtmed[4],dm,4);
1083
1084
1085 //------------------------------------------------
1086 // Tpc Inner readout chamber Pad Plane
1087 //------------------------------------------------
1088
1089 dm[0]= 14.483;
1090 dm[1]= 23.3345;
1091 dm[2]= 0.5;
1092 dm[3]= 25.1;
1093
1094 gMC->Gsvolu("TIPP","TRD1",idtmed[12],dm,4);
1095
1096 //
1097
1098 dm[0] -= 1.218511934;
1099 dm[1] -= 1.218511934;
1100 dm[2] = 0.35;
1101
1102 gMC->Gsvolu("TIC3","TRD1",idtmed[1],dm,4);
1103
1104 gMC->Gspos("TIC3",1,"TIPP",0.,0.15,0.,0,"ONLY");
1105
1106 gMC->Gspos("TIPP",1,"TIRC",0.,1.1,0.,0,"ONLY");
1107
1108
1109 //----------------------------------------------
1110 // Tpc Readout Chambers Empty spaces - for both
1111 // inner and outer sectors
1112 //----------------------------------------------
1113
1114 gMC->Gsvolu("TRCE","TRD1",idtmed[0],dm,0);
1115
1116 // Inner sector - 4 spaces
1117
1118
1119 dm[3] = 4.7625;
1120 dm[0] = 12.472;
1121
1122 Float_t rr = 90.52;
1123 Float_t zz;
1124
1125 zz= -12.7875;
1126
1127 space = rr*tanAlpha-dm[0];
1128
1129 for(Int_t nsLow=0;nsLow<4;nsLow++){
1130
1131 rr += 9.525;
1132 dm[1]= rr*tanAlpha - space;
1133
1134 dm[2]=0.8;
1135
1136 gMC->Gsposp("TRCE",nsLow+1,"TIRC",0.,-0.8,zz,0,"ONLY",dm,4);
1137
1138 //
1139
1140 dm[2]= 1.2;
1141
1142 gMC->Gsposp("TRCE",nsLow+5,"TIC1",0.,1.05,zz-2.1,0,"ONLY",dm,4);
1143
1144 rr += 0.4;
1145 dm[0] = rr*tanAlpha - space;
1146 zz += (0.4+9.525);
1147
1148 }
1149
1150 dm[0]= 12.472;
1151 // dm[1] - this is the dm[1] from the previous TRCE
1152 dm[2]= 1.05;
1153 dm[3]= 19.65;
1154
1155 gMC->Gsposp("TRCE",9,"TIC1",0.,-1.2,0.,0,"ONLY",dm,4);
1156
1157 //
1158 // TPc Space for Connectors
1159 //
1160
1161 dm[0]= .3;
1162 dm[1]= .3;
1163 dm[2]= 4.5;
1164
1165 gMC->Gsvolu("TPSC","BOX ",idtmed[0],dm,3);
1166
1167 // TPC Connectors
1168
1169 dm[0]= .25;
1170 dm[1]= .15;
1171 dm[2]= 3.75;
1172
1173 gMC->Gsvolu("TPCC","BOX ",idtmed[13],dm,3);
1174
1175 gMC->Gspos("TPCC",1,"TPSC",0.,0.15,0.,0,"ONLY");
1176
1177 zz = -12.7875;
1178
1179
1180 Float_t alpha;
1181 Float_t astep;
1182
1183 Float_t phi1,phi2,phi3,theta1,theta2,theta3; // rotation angles
1184
1185 // inner part of the inner sector - 2 x 20 holes
1186
1187 astep = 20.00096874/19.;
1188
1189 alpha = 10.00048437-astep;
1190
1191 Float_t x1,x2;
1192
1193 x1 = 13.31175725;
1194 x1 -= 0.996357832;
1195
1196 x2 = 15.06180253;
1197 x2 -= 1.163028812;
1198
1199 Int_t ncon;
1200
1201 for(ncon=0;ncon<20;ncon++){
1202
1203 phi1 = 0.;
1204 theta1 = 90.+alpha;
1205 phi2=90.;
1206 theta2 = 90.;
1207 phi3 = (alpha>0) ? 0. : 180.;
1208 theta3 = TMath::Abs(alpha);
1209
1210 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1211
1212
1213
1214 gMC->Gspos("TPSC",ncon+1,"TIRC",x1,0.3,-12.7875,idrotm[nRotMat],"ONLY");
1215 gMC->Gspos("TPSC",ncon+21,"TIRC",x2,0.3,-2.8625,idrotm[nRotMat],"ONLY");
1216
1217
1218 x1 -= 1.296357833;
1219 x2 -= 1.463028812;
1220
1221 alpha -= astep;
1222 nRotMat++;
1223
1224 }
1225
1226 // outer part of the inner sector - 2 x 25 holes
1227
1228 astep = 20.00096874/24.;
1229 alpha = 10.00048437-astep;
1230
1231 x1 = 16.81184781;
1232 x1 -= 1.016295986;
1233
1234 x2 = 18.5618931;
1235 x2 -= 1.150914854;
1236
1237 for(ncon=0;ncon<25;ncon++){
1238
1239 phi1 = 0.;
1240 theta1 = 90.+alpha;
1241 phi2=90.;
1242 theta2 = 90.;
1243 phi3 = (alpha>0) ? 0. : 180.;
1244 theta3 = TMath::Abs(alpha);
1245
1246 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1247
1248
1249
1250 gMC->Gspos("TPSC",ncon+41,"TIRC",x1,0.3,7.0625,idrotm[nRotMat],"ONLY");
1251 gMC->Gspos("TPSC",ncon+66,"TIRC",x2,0.3,16.9875,idrotm[nRotMat],"ONLY");
1252
1253
1254 x1 -= 1.316295986;
1255 x2 -= 1.450914854;
1256
1257 alpha -= astep;
1258 nRotMat++;
1259
1260 }
1261
1262 //--------------------------------------------------------------------------
1263 // TPC Outer Readout Chambers
1264 // this is NOT a final design
1265 //--------------------------------------------------------------------------
1266
1267 dm[0]= 23.3875;
1268 dm[1]= 43.524;
1269 dm[2]= 1.5; //thickness
1270 dm[3]= 57.1;
1271
1272 gMC->Gsvolu("TORC","TRD1",idtmed[4],dm,4);
1273
1274 //------------------------------------------------
1275 // Tpc Outer readout chamber Pad Plane
1276 //------------------------------------------------
1277
1278 dm[2]= 0.5;
1279
1280 gMC->Gsvolu("TOPP","TRD1",idtmed[12],dm,4);
1281
1282 dm[0] -= 1.218511934;
1283 dm[1] -= 1.218511934;
1284 dm[2] = 0.35;
1285
1286 gMC->Gsvolu("TOC3","TRD1",idtmed[1],dm,4);
1287
1288 gMC->Gspos("TOC3",1,"TOPP",0.,0.15,0.,0,"ONLY");
1289
1290 gMC->Gspos("TOPP",1,"TORC",0.,1.0,0.,0,"ONLY");
1291
1292 // empty space
1293
1294
1295 dm[0]= 21.035;
1296 dm[1]= 38.7205;
1297 dm[2]= 0.7;
1298 dm[3]= 50.15;
1299
1300 gMC->Gsposp("TRCE",11,"TORC",0.,-0.8,-2.15,0,"ONLY",dm,4);
1301
1302 dm[0]= 22.2935;
1303 dm[1]= 40.5085;
1304 dm[2]= 2.25;
1305 dm[3]= 51.65;
1306
1307 gMC->Gsvolu("TOC1","TRD1",idtmed[4],dm,4);
1308
1309
1310 dm[0]= 21.35;
1311 dm[1]= 38.7205;
1312 dm[2]= 2.25;
1313 dm[3]= 50.15;
1314
1315 gMC->Gsposp("TRCE",11,"TOC1",0.,0.,0.,0,"ONLY",dm,4);
1316
1317
1318 //-----------------------------------------------
1319 // Tpc Services Support Wheel
1320 //-----------------------------------------------
1321
1322 dm[0]=0.;
1323 dm[1]=360.;
1324 dm[2]=18.;
1325 dm[3]=2.;
1326
1327 dm[4]= -5.;
1328 dm[5]= 77.017;
1329 dm[6]= 255.267;
1330
1331 dm[7]= 5.;
1332 dm[8]= dm[5];
1333 dm[9]= dm[6];
1334
1335 gMC->Gsvolu("TSSW","PGON",idtmed[4],dm,10);
1336
1337 // Tpc Services Wheel Cover
1338
1339 dm[4]= -0.5;
1340 dm[7]= 0.5;
1341
1342 gMC->Gsvolu("TSWC","PGON",idtmed[4],dm,10);
1343
1344 // Tpc Service wheel Cover Empty space
1345
1346 dm[0]= 10.99;
1347 dm[1]= 39.599;
1348 dm[2]= .5;
1349 dm[3]= 81.125;
1350
1351 gMC->Gsvolu("TSCE","TRD1",idtmed[0],dm,4);
1352
1353 // Tpc services Wheel Empty Spaces
1354
1355 dm[0]= 13.18017507;
1356 dm[1]= 44.61045938;
1357 dm[2]= 4.;
1358 dm[3]= 89.125;
1359
1360 gMC->Gsvolu("TWES","TRD1",idtmed[0],dm,4);
1361
1362 // Tpc Services Wheel Bars
1363
1364 gMC->Gsvolu("TSWB","TRD1",idtmed[4],dm,0);
1365
1366 // bars-> TWES
1367
1368 dm[2]= 4.;
1369 dm[3]= .4;
1370
1371 dm[0]= 13.8149522;
1372 dm[1]= 13.95601379;
1373
1374 gMC->Gsposp("TSWB",1,"TWES",0.,0.,-85.125,0,"ONLY",dm,4);
1375
1376 dm[0]= 43.83462067;
1377 dm[1]= 43.97568225;
1378
1379 gMC->Gsposp("TSWB",2,"TWES",0.,0.,85.125,0,"ONLY",dm,4);
1380
1381 // TPc ELectronics - right now 30% X0 Si
1382
1383 dm[0]= 14.03813696;
1384 dm[1]= 43.3524075;
1385 dm[2]= 1.404;
1386 dm[3]= 83.125;
1387
1388 gMC->Gsvolu("TPEL","TRD1",idtmed[11],dm,4);
1389 gMC->Gspos("TPEL",1,"TWES",0.,0.,0.,0,"ONLY");
1390
1391 //--------------------------------------------------------------------------
1392 // End caps
1393 //--------------------------------------------------------------------------
1394
1395 // TPc Main Wheel - Al
1396
1397 dm[0]= 75.3;
1398 dm[1]= 264.8;
1399 dm[2]= 3.0;
1400
1401 gMC->Gsvolu("TPMW","TUBE",idtmed[4],dm,3);
1402
1403 // TPc Extra Wheel (to avoid overlapping) - Al
1404
1405 dm[0]= 264.8;
1406 dm[1]= 277.0;
1407 dm[2]= 1.95;
1408
1409 gMC->Gsvolu("TPEW","TUBE",idtmed[4],dm,3);
1410
1411 //--------------------------------------------------------------------------
1412 // Tpc Empty Space for the Readout chambers
1413 //--------------------------------------------------------------------------
1414
1415 Float_t rLow= 86.2;
1416 Float_t rUp= 243.5;
1417 Float_t dR = 0.5*(rUp-rLow);
1418
1419 space= 1.4/cosAlpha; // wheel ribs are 2.8 cm wide
1420
1421 dm[0]= rLow*tanAlpha-space;
1422 dm[1]= rUp*tanAlpha-space;
1423 dm[2]= 3.0;
1424 dm[3]= dR;
1425
1426 gMC->Gsvolu("TESR","TRD1",idtmed[0],dm,4);
1427
1428 // TIC1->TESR
1429
1430
1431 gMC->Gspos("TIC1",1,"TESR",0.,0.75,-dR+23.97,0,"ONLY");
1432
1433
1434 // TOC1->TESR
1435
1436 gMC->Gspos("TOC1",1,"TESR",0.,0.75,dR-55.02,0,"ONLY");
1437
1438 // Tpc Empty Space Bars - Al (daughters of TESR)
1439
1440 Float_t zBar;
1441
1442 gMC->Gsvolu("TESB","TRD1",idtmed[4],dm,0);
1443
1444 // lower bar
1445
1446 dm[0]= rLow*tanAlpha-space;
1447 dm[1]= 88.7*tanAlpha-space;
1448 dm[2]= 0.95;
1449 dm[3]= 1.275;
1450
1451 zBar = -dR+dm[3];
1452
1453 gMC->Gsposp("TESB",1,"TESR",0.,2.05,zBar,0,"ONLY",dm,4);
1454
1455 // middle bar
1456
1457 dm[0]= 131.65*tanAlpha-space;
1458 dm[1]= 136.5*tanAlpha-space;
1459 dm[3]= 2.425;
1460
1461 zBar = -dR +131.65+dm[3]-rLow;
1462
1463 gMC->Gsposp("TESB",2,"TESR",0.,2.05,zBar,0,"ONLY",dm,4);
1464
1465 // upper bar
1466
1467 dm[0]= 240.4*tanAlpha-space;
1468 dm[1]= rUp*tanAlpha-space;
1469 dm[3]= 1.55;
1470
1471 zBar = dR-dm[3];
1472
1473 gMC->Gsposp("TESB",3,"TESR",0.,2.05,zBar,0,"ONLY",dm,4);
1474
1475 //------------------------------------------------------
1476 // TPc Lower "S" Sectors
1477 //------------------------------------------------------
1478
1479 Float_t deadSpace;
1480
1481
1482 Float_t inSecLowEdge = fTPCParam->GetInnerRadiusLow();
1483 Float_t inSecUpEdge = fTPCParam->GetInnerRadiusUp();
1484
1485 dm[0] = inSecLowEdge*TMath::Tan(0.5*innerOpenAngle)-0.01;
1486 dm[1] = inSecUpEdge*TMath::Tan(0.5*innerOpenAngle)-0.01;
1487 dm[2] = 0.5*(250. - 0.3);
1488 dm[3] = 0.5*(inSecUpEdge-inSecLowEdge);
1489
1490 gMC->Gsvolu("TPLS", "TRD1", idtmed[1], dm, 4); // nonsensitive
1491
1492 //----------------------------------------------------------
1493 // TPc Upper Sectors
1494 //----------------------------------------------------------
1495
1496 Float_t ouSecLowEdge = fTPCParam->GetOuterRadiusLow();
1497 Float_t ouSecUpEdge = fTPCParam->GetOuterRadiusUp();
1498
1499 dm[0] = ouSecLowEdge*TMath::Tan(0.5*outerOpenAngle)-0.01;
1500 dm[1] = ouSecUpEdge*TMath::Tan(0.5*outerOpenAngle)-0.01;
1501 dm[2] = 0.5*(250. - 0.3);
1502 dm[3] = 0.5*(ouSecUpEdge-ouSecLowEdge);
1503
1504 gMC->Gsvolu("TPUS", "TRD1", idtmed[1], dm, 4); // nonsensitive
1505
1506 //--------------------------------------------------------
1507 // TPc Sensitive Strips, the same for S and L sectors
1508 //--------------------------------------------------------
1509
1510 gMC->Gsvolu("TPSS","TRD1",idtmed[2],dm,0); // sensitive
1511
1512 Int_t nofStrips,nstr;
1513 Float_t r1,r2,zs;
1514 Float_t stripThick = 0.01; // 100 microns
1515
1516 // inner sector
1517
1518 nofStrips = fTPCParam->GetNRowLow();
1519 deadSpace = fTPCParam->GetInnerWireMount();
1520
1521 dm[2] = 0.5*(250. - 0.3);
1522 dm[3] = 0.5 * stripThick;
1523
1524 for(nstr=0;nstr<nofStrips;nstr++){
1525
1526 r1 = fTPCParam->GetPadRowRadiiLow(nstr);
1527 r2 = r1 + stripThick;
1528 dm[0] = r1 * TMath::Tan(0.5*innerOpenAngle) - deadSpace;
1529 dm[1] = r2 * TMath::Tan(0.5*innerOpenAngle) - deadSpace;
1530 zs = -inSecLowEdge -0.5*(inSecUpEdge-inSecLowEdge);
1531 zs += r1;
1532 zs += dm[3];
1533
1534 gMC->Gsposp("TPSS", nstr+1, "TPLS", 0., 0., zs, 0, "ONLY", dm, 4);
1535
1536 }
1537
1538 Int_t nsSave = nofStrips;
1539
1540 // outer sector
1541
1542 nofStrips = fTPCParam->GetNRowUp();
1543 deadSpace = fTPCParam->GetOuterWireMount();
1544
1545 dm[2] = 0.5*(250. - 0.3);
1546 dm[3] = 0.5 * stripThick;
1547
1548 for(nstr=0;nstr<nofStrips;nstr++){
1549
1550 r1 = fTPCParam->GetPadRowRadiiUp(nstr);
1551 r2 = r1 + stripThick;
1552 dm[0] = r1 * TMath::Tan(0.5*outerOpenAngle) - deadSpace;
1553 dm[1] = r2 * TMath::Tan(0.5*outerOpenAngle) - deadSpace;
1554 zs = -ouSecLowEdge -0.5*(ouSecUpEdge-ouSecLowEdge);
1555 zs += r1;
1556 zs += dm[3];
1557
1558 gMC->Gsposp("TPSS", nstr+1+nsSave, "TPUS", 0., 0., zs, 0, "ONLY", dm, 4);
1559
1560 }
1561
1562
1563 //-------------------------------------------------------
1564 // positioning of the empty spaces into the main wheel
1565 // and readout chambers and sectors into the drift gas
1566 //-------------------------------------------------------
1567
1568 Float_t rCenter,xc,yc;
1569 Float_t rInner,rOuter; // center of the inner and outer chamber
1570
1571 rCenter = rLow+dR;
1572
1573 rInner = 108.07;
1574 rOuter = 190.68;
1575
1576 for(Int_t ns=0; ns<nInnerSector;ns++){
1577
1578 phi1 = ns * innerOpenAngle + innerAngleShift;
1579 phi1 *= kRaddeg; // in degrees
1580
1581 phi1 = (Float_t)TMath::Nint(phi1) + 270.;
1582
1583 if (phi1 > 360.) phi1 -= 360.;
1584
1585 theta1 = 90.;
1586 phi2 = 90.;
1587 theta2 = 180.;
1588 phi3 = ns * innerOpenAngle + innerAngleShift;
1589 phi3 *= kRaddeg; // in degrees
1590
1591 phi3 = (Float_t)TMath::Nint(phi3);
1592
1593 if(phi3 > 360.) phi3 -= 360.;
1594
1595 theta3 = 90.;
1596
1597 // "holes"->End plate
1598
1599 xc = rCenter*TMath::Cos(phi3*kDegrad);
1600 yc = rCenter*TMath::Sin(phi3*kDegrad);
1601
1602 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1603
1604 gMC->Gspos("TESR",ns+1,"TPMW",xc,yc,0.,idrotm[nRotMat],"ONLY");
1605
1606 // TSCE->TSWC (services wheel volumes)
1607
1608 xc = 166.142*TMath::Cos(phi3*kDegrad);
1609 yc = 166.142*TMath::Sin(phi3*kDegrad);
1610
1611 gMC->Gspos("TSCE",ns+1,"TSWC",xc,yc,0.,idrotm[nRotMat],"ONLY");
1612 gMC->Gspos("TWES",ns+1,"TSSW",xc,yc,0.,idrotm[nRotMat],"ONLY");
1613
1614
1615 // readout chambers->TDGN (drift gas)
1616
1617 xc = rInner*TMath::Cos(phi3*kDegrad);
1618 yc = rInner*TMath::Sin(phi3*kDegrad);
1619
1620 gMC->Gspos("TIRC",ns+1,"TDGN",xc,yc,252.,idrotm[nRotMat],"ONLY");
1621
1622 // here lower sectors
1623
1624 gMC->Gspos("TPLS",ns+1,"TDGN",xc,yc,125.15,idrotm[nRotMat],"ONLY");
1625 gMC->Gspos("TPLS",ns+nInnerSector+1,"TDGN",xc,yc,-125.15,idrotm[nRotMat],"ONLY");
1626
1627 xc = rOuter*TMath::Cos(phi3*kDegrad);
1628 yc = rOuter*TMath::Sin(phi3*kDegrad);
1629
1630 gMC->Gspos("TORC",ns+1,"TDGN",xc,yc,252.1,idrotm[nRotMat],"ONLY");
1631
1632 // here upper sectors
1633
1634 gMC->Gspos("TPUS",ns+1,"TDGN",xc,yc,125.15,idrotm[nRotMat],"ONLY");
1635 gMC->Gspos("TPUS",ns+nOuterSector+1,"TDGN",xc,yc,-125.15,idrotm[nRotMat],"ONLY");
1636
1637 nRotMat++;
1638
1639 theta2 = 0.; // reflection
1640
1641 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1642
1643 xc = rInner*TMath::Cos(phi3*kDegrad);
1644 yc = rInner*TMath::Sin(phi3*kDegrad);
1645
1646 gMC->Gspos("TIRC",ns+nInnerSector+1,"TDGN",xc,yc,-252.,idrotm[nRotMat],"ONLY");
1647
1648 xc = rOuter*TMath::Cos(phi3*kDegrad);
1649 yc = rOuter*TMath::Sin(phi3*kDegrad);
1650
1651 gMC->Gspos("TORC",ns+nOuterSector+1,"TDGN",xc,yc,-252.1,idrotm[nRotMat],"ONLY");
1652
1653 nRotMat++;
1654
1655 }
1656
1657 // reflection matrix
1658
1659 theta1 = 90.;
1660 phi1 = 0.;
1661 theta2 = 90.;
1662 phi2 = 270.;
1663 theta3 = 180.;
1664 phi3 = 0.;
1665
1666 AliMatrix(idrotm[nRotMat], theta1, phi1, theta2, phi2, theta3, phi3);
1667
1668
1669 // TPMW->TPC
1670
1671 gMC->Gspos("TPMW",1,"TPC ",0.,0.,256.6,0,"ONLY");
1672 gMC->Gspos("TPMW",2,"TPC ",0.,0.,-256.6,idrotm[nRotMat],"ONLY");
1673 gMC->Gspos("TPEW",1,"TPC ",0.,0.,257.65,0,"ONLY");
1674 gMC->Gspos("TPEW",2,"TPC ",0.,0.,-257.5,0,"ONLY");
1675
1676
1677
1678 //-------------------------------------------------------
1679 // Tpc High Voltage Membrane - NOMEX honeycomb
1680 //-------------------------------------------------------
1681
1682 dm[0]=0.,
1683 dm[1]=360.;
1684 dm[2]=18.;
1685 dm[3]=2.;
1686
1687 //
1688
1689 dm[4]= -0.3;
1690 dm[5]= 81.156;
1691 dm[6]= 253.386;
1692
1693 //
1694
1695 dm[7]= 0.3;
1696 dm[8]= dm[5];
1697 dm[9]= dm[6];
1698
1699 gMC->Gsvolu("THVM","PGON",idtmed[6],dm,10);
1700
1701 gMC->Gspos("THVM",1,"TDGN",0.,0.,0.,0,"ONLY");
1702
1703 //----------------------------------------------------------
1704 // TPc Support Rods - MAKROLON
1705 //----------------------------------------------------------
1706
1707 dm[0]= 0.9;
1708 dm[1]= 1.2;
1709 dm[2]= 126.65;
1710
1711 gMC->Gsvolu("TPSR","TUBE",idtmed[7],dm,3);
1712
1713 for(Int_t nrod=1;nrod<18;nrod++){
1714 Float_t angle=innerOpenAngle*(Float_t)nrod;
1715
1716 xc=82.4*TMath::Cos(angle);
1717 yc=82.4*TMath::Sin(angle);
1718
1719 gMC->Gspos("TPSR",nrod,"TDGN",xc,yc,126.95,0,"ONLY");
1720 gMC->Gspos("TPSR",nrod+17,"TDGN",xc,yc,-126.95,0,"ONLY");
1721
1722 xc=254.2*TMath::Cos(angle);
1723 yc=254.2*TMath::Sin(angle);
1724
1725 gMC->Gspos("TPSR",nrod+34,"TDGN",xc,yc,126.95,0,"ONLY");
1726 gMC->Gspos("TPSR",nrod+51,"TDGN",xc,yc,-126.95,0,"ONLY");
1727
1728 }
1729
1730 //----------------------------------------------------------
1731 // Tpc High Voltage rod - MAKROLON + Copper cable
1732 //----------------------------------------------------------
1733
1734 // rod with cable (Left)
1735
1736 dm[0]=0.;
1737 dm[1]=2.25;
1738 dm[2]=126.65;
1739
1740 gMC->Gsvolu("THVL","TUBE",idtmed[7],dm,3);
1741
1742 // HV cable
1743
1744 dm[0]=0.;
1745 dm[1]=0.3;
1746 dm[2]=126.65;
1747
1748 gMC->Gsvolu("THVC","TUBE",idtmed[10],dm,3);
1749
1750 // empty space
1751
1752 dm[0]=0.3;
1753 dm[1]=1.;
1754 dm[2]=126.65;
1755
1756 gMC->Gsvolu("THVE","TUBE",idtmed[1],dm,3);
1757
1758 gMC->Gspos("THVC",1,"THVL",0.,0.,0.,0,"ONLY");
1759 gMC->Gspos("THVE",1,"THVL",0.,0.,0.,0,"ONLY");
1760
1761 // rod without cable
1762
1763 dm[0]=1.8;
1764 dm[1]=2.25;
1765 dm[2]=126.65;
1766
1767 gMC->Gsvolu("THVR","TUBE",idtmed[7],dm,3);
1768
1769
1770
1771 gMC->Gspos("THVL",1,"TDGN",82.4,0.,-126.95,0,"ONLY");
1772 gMC->Gspos("THVL",2,"TDGN",254.2,0.,-126.95,0,"ONLY");
1773
1774 gMC->Gspos("THVR",1,"TDGN",82.4,0.,126.95,0,"ONLY");
1775 gMC->Gspos("THVR",2,"TDGN",254.2,0.,126.95,0,"ONLY");
1776
1777
1778
1779 gMC->Gspos("TDGN",1,"TPC ",0.,0.,0.,0,"ONLY");
1780
1781 // services wheel cover -> wheel
1782
1783
1784 gMC->Gspos("TSWC",1,"TSSW",0.,0.,4.5,0,"ONLY");
1785 gMC->Gspos("TSWC",2,"TSSW",0.,0.,-4.5,0,"ONLY");
1786
1787
1788 // put the wheel into the TPC
1789
1790 gMC->Gspos("TSSW",1,"TPC ",0.,0.,278.7,0,"ONLY");
1791 gMC->Gspos("TSSW",2,"TPC ",0.,0.,-278.7,0,"ONLY");
1792
1793 gMC->Gsord("TPMW",6);
1794 gMC->Gsord("TPLS",3);
1795 gMC->Gsord("TPUS",3);
1796 gMC->Gsord("TDGN",6);
1797 gMC->Gsord("TSSW",6);
1798 gMC->Gsord("TSWC",6);
1799
1800 // put the TPC into ALIC (main mother volume)
1801
1802 gMC->Gspos("TPC ", 1, "ALIC", 0, 0, 0, 0, "ONLY");
1803
1804
1805} // end of function
1806
1807//_____________________________________________________________________________
1808void AliTPCv1::DrawDetector()
1809{
1810 //
1811 // Draw a shaded view of the Time Projection Chamber version 1
1812 //
1813
1814 // Set everything unseen
1815 gMC->Gsatt("*", "seen", -1);
1816 //
1817 // Set ALIC mother transparent
1818 gMC->Gsatt("ALIC","SEEN",0);
1819 //
1820 // Set the volumes visible
1821 gMC->Gsatt("TPC ","SEEN",0);
1822 gMC->Gsatt("TOIN","SEEN",1);
1823 gMC->Gsatt("TOIN","COLO",7);
1824 gMC->Gsatt("TPCR","SEEN",0);
1825 gMC->Gsatt("TOCV","SEEN",1);
1826 gMC->Gsatt("TOCV","COLO",4);
1827 gMC->Gsatt("TSA1","SEEN",0);
1828 gMC->Gsatt("TSA2","SEEN",0);
1829 gMC->Gsatt("TSA3","SEEN",0);
1830 gMC->Gsatt("TSA4","SEEN",0);
1831 gMC->Gsatt("TOFC","SEEN",1);
1832 gMC->Gsatt("TOFC","COLO",4);
1833 gMC->Gsatt("TSA5","SEEN",0);
1834 gMC->Gsatt("TSA6","SEEN",0);
1835 gMC->Gsatt("TSA7","SEEN",0);
37831078 1836 gMC->Gsatt("TIIN","COLO",7);
1837 gMC->Gsatt("TIIN","SEEN",1);
1838 gMC->Gsatt("TICL","SEEN",0);
1839 gMC->Gsatt("TSA9","SEEN",0);
1840 gMC->Gsatt("TS10","SEEN",0);
1841 gMC->Gsatt("TS11","SEEN",0);
1842 gMC->Gsatt("TS12","SEEN",0);
1843 gMC->Gsatt("TICR","SEEN",0);
1844 gMC->Gsatt("TS13","SEEN",0);
1845 gMC->Gsatt("TS14","SEEN",0);
1846 gMC->Gsatt("TS15","SEEN",0);
1847 gMC->Gsatt("TS16","SEEN",0);
1848 gMC->Gsatt("TIFC","SEEN",1);
1849 gMC->Gsatt("TIFC","COLO",4);
1850 gMC->Gsatt("TS17","SEEN",0);
1851 gMC->Gsatt("TS18","SEEN",0);
1852 gMC->Gsatt("TS19","SEEN",0);
37831078 1853 gMC->Gsatt("TS21","SEEN",0);
1854 gMC->Gsatt("TS22","SEEN",0);
1855 gMC->Gsatt("TS23","SEEN",0);
1856 gMC->Gsatt("TS24","SEEN",0);
1857 gMC->Gsatt("TDGN","SEEN",0);
1858 gMC->Gsatt("TIRC","SEEN",0);
1859 gMC->Gsatt("TIC1","SEEN",1);
1860 gMC->Gsatt("TIPP","SEEN",0);
1861 gMC->Gsatt("TIC3","SEEN",0);
1862 gMC->Gsatt("TRCE","SEEN",0);
1863 gMC->Gsatt("TPSC","SEEN",0);
1864 gMC->Gsatt("TPCC","SEEN",0);
1865 gMC->Gsatt("TORC","SEEN",0);
1866 gMC->Gsatt("TOPP","SEEN",0);
1867 gMC->Gsatt("TOC3","SEEN",0);
1868 gMC->Gsatt("TOC1","SEEN",1);
1869 gMC->Gsatt("TSSW","SEEN",1);
1870 gMC->Gsatt("TSWC","SEEN",1);
1871 gMC->Gsatt("TSCE","SEEN",1);
1872 gMC->Gsatt("TSSW","COLO",3);
1873 gMC->Gsatt("TSWC","COLO",3);
1874 gMC->Gsatt("TSCE","COLO",6);
1875 gMC->Gsatt("TWES","SEEN",0);
1876 gMC->Gsatt("TSWB","SEEN",0);
1877 gMC->Gsatt("TPEL","SEEN",0);
1878 gMC->Gsatt("TPMW","SEEN",1);
1879 gMC->Gsatt("TPEW","SEEN",1);
1880 gMC->Gsatt("TESR","SEEN",1);
1881 gMC->Gsatt("TPMW","COLO",12);
1882 gMC->Gsatt("TPEW","COLO",12);
1883 gMC->Gsatt("TWES","COLO",5);
1884 gMC->Gsatt("TIC1","COLO",5);
1885 gMC->Gsatt("TOC1","COLO",5);
1886 gMC->Gsatt("TESB","SEEN",0);
1887 gMC->Gsatt("TPLS","SEEN",0);
1888 gMC->Gsatt("TPUS","SEEN",0);
1889 gMC->Gsatt("TPSS","SEEN",0);
1890 gMC->Gsatt("THVM","SEEN",1);
1891 gMC->Gsatt("THVM","COLO",11);
1892 gMC->Gsatt("TPSR","SEEN",0);
1893 gMC->Gsatt("THVL","SEEN",0);
1894 gMC->Gsatt("THVC","SEEN",0);
1895 gMC->Gsatt("THVE","SEEN",0);
1896 gMC->Gsatt("THVR","SEEN",0);
1897
1898 gMC->Gdopt("hide", "on");
1899 gMC->Gdopt("shad", "on");
1900 gMC->Gsatt("*", "fill", 7);
1901 gMC->SetClipBox(".");
1902 gMC->SetClipBox("TPMW",-300,300,-300,300,254.,270.);
1903 gMC->SetClipBox("TESR",-300,300,-300,300,254.,270.);
1904 gMC->SetClipBox("TSSW",-300,300,-300,300,283.,284.);
1905 gMC->SetClipBox("TSWC",-300,300,-300,300,283.,284.);
1906 gMC->SetClipBox("*", 0, 300, -300, 300, -290, 290);
1907 gMC->DefaultRange();
1908 gMC->Gdraw("alic", 40, 30, 0, 12, 9.5, .025, .025);
1909 gMC->Gdhead(1111, "Time Projection Chamber");
1910 gMC->Gdman(18, 4, "MAN");
1911 gMC->Gdopt("hide","off");
1912}
1913
1914//_____________________________________________________________________________
1915void AliTPCv1::CreateMaterials()
1916{
1917 //
1918 // Define materials for Time Projection Chamber
1919 //
1920 AliTPC::CreateMaterials();
1921}
1922
1923//_____________________________________________________________________________
1924void AliTPCv1::Init()
fe4da5cc 1925{
1926 //
1927 // Initialises TPC detector after it has been created
1928 //
cfce8870 1929
37831078 1930 fIdSens=gMC->VolId("TPSS"); // sensitive strip
1931 fIdLSec=gMC->VolId("TPLS"); // lower sector
1932 fIdUSec=gMC->VolId("TPUS"); // upper sector
fe4da5cc 1933
1934 printf("TPC version 1 initialized\n");
1935}
1936
1937//_____________________________________________________________________________
1938void AliTPCv1::StepManager()
1939{
1940 //
1941 // Called at every step in the Time Projection Chamber
1942 //
0a6d8768 1943 Int_t copy, id, i;
fe4da5cc 1944 Float_t hits[4];
1945 Int_t vol[2];
0a6d8768 1946 TLorentzVector p;
6a905708 1947 //
1283eee5 1948
0a6d8768 1949 if(gMC->TrackCharge() && gMC->IsTrackEntering()) {
37831078 1950
1951 //
1952 // Only charged tracks entering the sensitive volume
fe4da5cc 1953 //
37831078 1954
1955 id = gMC->CurrentVolID(copy);
1956
1957 if(id == fIdSens){
1958
1959 // sensitive volume (strip)
1960
1961 vol[1]= copy-1; // row number (absolute)
1962
1963 // sector type
1964
1965 id = gMC->CurrentVolOffID(1,copy);
1966
1967 if(id == fIdLSec){
1968
1969 // lower sector
1970
1971 vol[0] = copy-1; // sector number
1972
cc80f89e 1973 }
37831078 1974 else {
1975
1976 // upper sector
cc80f89e 1977
37831078 1978 vol[0] = copy-1+fTPCParam->GetNInnerSector(); // sector number
1979 vol[1] -= fTPCParam->GetNRowLow(); // row number
1283eee5 1980
37831078 1981 }
1982
1983 if(vol[1] == 0){
1984
1985 // because Jouri wants to have this
1986
1987 gMC->TrackMomentum(p);
1988 hits[0]=p[0];
1989 hits[1]=p[1];
1990 hits[2]=p[2];
1991 hits[3]=0.; // this hit has no energy loss
6a905708 1992 // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits);
1993
1994 AddHit(gAlice->CurrentTrack(), vol,hits); // M.I.
37831078 1995
1996 gMC->TrackPosition(p);
1997 hits[0]=p[0];
1998 hits[1]=p[1];
1999 hits[2]=p[2];
2000 hits[3]=0.; // this hit has no energy loss
6a905708 2001 // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits);
2002
2003 AddHit(gAlice->CurrentTrack(), vol,hits); // M.I.
37831078 2004
2005 }
2006
2007 gMC->TrackPosition(p);
2008 for(i=0;i<3;++i) hits[i]=p[i];
2009 hits[3]=1; //I'd like to have something positive here (I.Belikov)
6a905708 2010 // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits);
2011
2012 AddHit(gAlice->CurrentTrack(), vol,hits); // M.I.
37831078 2013
2014 }
2015
2016 }
2017
fe4da5cc 2018}
37831078 2019//_____________________________________________________________________________
1283eee5 2020
2021
2022