| 53d2a7cd |
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$ */ |
| 53d2a7cd |
17 | |
| 18 | /////////////////////////////////////////////////////////////////////////////// |
| 19 | // // |
| 20 | // ALICE Cosmic Ray Trigger // |
| 21 | // // |
| 22 | // This class contains the functions for version 0 of the ALICE Cosmic Ray // |
| 23 | // Trigger. This vesion is suposed to work as standalone module // |
| 24 | // // |
| 25 | // |
| 26 | // Authors: |
| 27 | // |
| 28 | // Arturo Fernandez <afernand@fcfm.buap.mx> |
| 29 | // Enrique Gamez <egamez@fcfm.buap.mx> |
| 30 | // |
| 31 | // Universidad Autonoma de Puebla |
| 32 | // |
| 33 | // |
| 34 | //Begin_Html |
| 35 | /* |
| 36 | <img src="picts/AliCRTv1Class.gif"> |
| 37 | </pre> |
| 38 | <br clear=left> |
| 39 | <p>The responsible person for this module is |
| 40 | <a href="mailto:egamez@fcfm.buap.mx">Enrique Gamez</a>. |
| 41 | </font> |
| 42 | <pre> |
| 43 | */ |
| 44 | //End_Html |
| 45 | // // |
| 46 | /////////////////////////////////////////////////////////////////////////////// |
| 47 | |
| 53d2a7cd |
48 | |
| 53d2a7cd |
49 | #include <TBRIK.h> |
| 116cbefd |
50 | #include <TGeometry.h> |
| 53d2a7cd |
51 | #include <TLorentzVector.h> |
| 116cbefd |
52 | #include <TNode.h> |
| 53 | #include <TPDGCode.h> |
| 53d2a7cd |
54 | |
| 53d2a7cd |
55 | #include "AliCRTConstants.h" |
| 116cbefd |
56 | #include "AliCRTv1.h" |
| 57 | #include "AliConst.h" |
| 58 | #include "AliMagF.h" |
| 59 | #include "AliRun.h" |
| 53d2a7cd |
60 | |
| 61 | ClassImp(AliCRTv1) |
| 62 | |
| 63 | //_____________________________________________________________________________ |
| 64 | AliCRTv1::AliCRTv1() : AliCRTv0() |
| 65 | { |
| 66 | // |
| 67 | // Default constructor for CRT |
| 68 | // |
| 69 | fCRTStatus = kTRUE; |
| 70 | fRICHStatus = kFALSE; |
| 71 | fTPCStatus = kFALSE; |
| 72 | fMagnetStatus = kTRUE; |
| 73 | |
| 74 | fCRTModule = kFALSE; |
| 75 | } |
| 76 | |
| 77 | //_____________________________________________________________________________ |
| 78 | AliCRTv1::AliCRTv1(const char *name, const char *title) |
| 79 | : AliCRTv0(name,title) |
| 80 | { |
| 81 | // |
| 82 | // Standard constructor for CRT |
| 83 | // |
| 84 | //Begin_Html |
| 85 | /* |
| 86 | <img src="picts/AliCRTv1.gif"> |
| 87 | */ |
| 88 | //End_Html |
| 89 | fCRTStatus = kTRUE; |
| 90 | fCRTModule = kFALSE; |
| 91 | |
| 92 | fRICHStatus = kFALSE; |
| 93 | fTPCStatus = kFALSE; |
| 94 | fMagnetStatus = kFALSE; |
| 95 | } |
| 96 | |
| 97 | //_____________________________________________________________________________ |
| 98 | AliCRTv1::AliCRTv1(const AliCRTv1& crt) |
| 99 | { |
| 100 | // |
| 101 | // Copy ctor. |
| 102 | // |
| 103 | crt.Copy(*this); |
| 104 | } |
| 105 | |
| 106 | //_____________________________________________________________________________ |
| 107 | AliCRTv1& AliCRTv1::operator= (const AliCRTv1& crt) |
| 108 | { |
| 109 | // |
| 110 | // Asingment operator |
| 111 | // |
| 112 | crt.Copy(*this); |
| 113 | return *this; |
| 114 | } |
| 115 | |
| 116 | //_____________________________________________________________________________ |
| 117 | void AliCRTv1::CreateGeometry() |
| 118 | { |
| 119 | // |
| 120 | // Create geometry for the CRT array |
| 121 | // |
| 122 | |
| 123 | Int_t idrotm[2499]; // The rotation matrix. |
| 124 | |
| 125 | Int_t * idtmed = fIdtmed->GetArray() - 1099 ; |
| 126 | |
| 127 | // |
| 128 | // Shafts. |
| 129 | this->CreateShafts(); |
| 130 | |
| 131 | // |
| 132 | // Molasse. |
| 133 | this->CreateMolasse(); |
| 134 | |
| 135 | |
| 136 | // |
| 137 | // Scintillators |
| 138 | |
| 139 | Float_t box[3]; |
| 140 | box[0] = AliCRTConstants::fgCageLenght/2.; // Half Length of the box along the X axis, cm. |
| 141 | box[1] = AliCRTConstants::fgCageHeight/2.; // Half Length of the box along the Y axis, cm. |
| 142 | box[2] = AliCRTConstants::fgCageWidth/2.; // Half Length of the box along the Z axis, cm. |
| 143 | |
| 144 | // |
| 145 | // Create a big voluem with air barrel above the magnet |
| 146 | Float_t barrel[10]; |
| 147 | Float_t magnetSides = 3.; |
| 148 | Float_t planesPerpendicularToZ = 2.; |
| 149 | Float_t rMin = 790.; |
| 150 | Float_t rMax = rMin + 20.; // 20 cm width |
| 151 | barrel[0] = 22.5; |
| 152 | barrel[1] = 45*magnetSides; |
| 153 | barrel[2] = magnetSides; |
| 154 | barrel[3] = planesPerpendicularToZ; |
| 155 | barrel[4] = -600.; |
| 156 | barrel[5] = rMin; |
| 157 | barrel[6] = rMax; |
| 158 | barrel[7] = 600.; |
| 159 | barrel[8] = rMin; |
| 160 | barrel[9] = rMax; |
| 161 | gMC->Gsvolu("CRT4", "PGON", idtmed[1114], barrel, 10); |
| 162 | gMC->Gspos("CRT4", 1 , "CRT", 0., -30., 0., 0, "ONLY"); |
| 163 | |
| 164 | |
| 165 | // Create the current sicuiitllator arry |
| 166 | // Define the Scintillators. as a big box. |
| 167 | Float_t scint[3]; |
| 168 | scint[0] = AliCRTConstants::fgActiveAreaLenght/2.; // Half Length in X |
| 169 | scint[1] = AliCRTConstants::fgActiveAreaHeight/2.; // Half Length in Y |
| 170 | scint[2] = AliCRTConstants::fgActiveAreaWidth/2.; // Half Length in Z |
| 171 | gMC->Gsvolu("CRT1", "BOX ", idtmed[1112], scint, 3); // Scintillators |
| 172 | // |
| 173 | // -- X axis. |
| 174 | // we'll start dawing from the center. |
| 175 | Float_t initX = 0.; |
| 176 | |
| 177 | // |
| 178 | // -- Y axis |
| 179 | Float_t gapY = 30.; // 30 cms. above the barrel. |
| 180 | // For the height we staimate the from the center of the ceiling, |
| 181 | // if were a cilinder, must be about 280cm. |
| 182 | Float_t barrelc = 790.; // Barrel radius. |
| 183 | Float_t height = barrelc + gapY - 30.; |
| 184 | Float_t initY = height; |
| 185 | |
| 186 | // |
| 187 | // -- Z axis. |
| 188 | // we'll start dawing from the center. |
| 189 | |
| 190 | // |
| 191 | // Put 4 modules on the top of the magnet |
| 192 | Int_t step = 4; |
| 193 | for ( Int_t i = 1 ; i <= 4 ; i++ ) { |
| 194 | gMC->Gspos("CRT1", i, "CRT", initX, initY, (i-step)*box[2], 0, "ONLY"); |
| 195 | step--; |
| 196 | } |
| 197 | |
| 198 | // Modules on the barrel sides. |
| 199 | // Because the openenig angle for each face is 22.5, and if we want to |
| 200 | // put the modules right in the middle |
| 201 | Float_t xtragap = 10.; |
| 202 | Float_t initXside = (height+xtragap)*TMath::Sin(2*22.5*kDegrad);//rigthside |
| 203 | Float_t initYside = (height+xtragap)*TMath::Cos(2*22.5*kDegrad); |
| 204 | |
| 205 | // Put 4 modules on the left side of the magnet |
| 206 | // The rotation matrix parameters, for the left side. |
| 207 | AliMatrix(idrotm[232], 90., 315., 90., 45., 0., 337.5); |
| 208 | Int_t stepl = 4; |
| 209 | for ( Int_t i = 1 ; i <= 4 ; i++ ) { |
| 210 | gMC->Gspos("CRT1", i+4, "CRT", initXside, initYside, (i-stepl)*box[2], |
| 211 | idrotm[232], "ONLY"); |
| 212 | stepl--; |
| 213 | } |
| 214 | |
| 215 | // Put 4 modules on the right side of the magnet |
| 216 | // The rotation matrix parameters for the right side. |
| 217 | AliMatrix(idrotm[231], 90., 45., 90., 315., 180., 202.5); |
| 218 | Int_t stepr = 4; |
| 219 | for ( Int_t i = 1 ; i <= 4 ; i++ ) { |
| 220 | gMC->Gspos("CRT1", i+8, "CRT", -initXside, initYside, (i-stepr)*box[2], |
| 221 | idrotm[231], "ONLY"); |
| 222 | stepr--; |
| 223 | } |
| 224 | |
| 225 | this->CreateMagnetGeometry(); |
| 226 | this->CreateRICHGeometry(); |
| 227 | this->CreateTPCGeometry(); |
| 228 | |
| 229 | } |
| 230 | |
| 231 | //_____________________________________________________________________________ |
| 232 | void AliCRTv1::CreateMagnetGeometry() |
| 233 | { |
| 234 | |
| 235 | cout<<"\n\n\tYou are requiring the CRT with the Magnet Activated!\n\n"; |
| 236 | |
| 237 | Int_t idrotm[2499]; // The rotation matrix. |
| 238 | |
| 239 | Int_t * idtmed = fIdtmed->GetArray() - 1099 ; |
| 240 | |
| 241 | // Disable the CRT StepManager method. |
| 242 | fCRTStatus = kFALSE; |
| 243 | |
| 244 | Float_t barrel[10]; |
| 245 | Float_t magnetSides = 3.; |
| 246 | Float_t planesPerpendicularToZ = 2.; |
| 247 | //Float_t rMin = 790.; |
| 248 | //Float_t rMax = rMin + 20.; // 20 cm width |
| 249 | |
| 250 | // MAgnet |
| 251 | // Create the upper faces of the magnet. |
| 252 | barrel[0] = 22.5; |
| 253 | barrel[1] = 360.; |
| 254 | barrel[2] = 8.; |
| 255 | barrel[3] = 2.; |
| 256 | barrel[4] = -600.; |
| 257 | barrel[5] = 580.; |
| 258 | barrel[6] = 790.; |
| 259 | barrel[7] = 600.; |
| 260 | barrel[8] = 580.; |
| 261 | barrel[9] = 790.; |
| 262 | gMC->Gsvolu("C3MO", "PGON", idtmed[1114], barrel, 10); |
| 263 | gMC->Gspos("C3MO", 1, "CRT", 0., -30., 0., 0, "ONLY"); |
| 264 | |
| 265 | // Define coils |
| 266 | |
| 267 | barrel[5] = 585.; |
| 268 | barrel[6] = 690.; |
| 269 | barrel[8] = 585.; |
| 270 | barrel[9] = 690.; |
| 271 | gMC->Gsvolu("C3CO", "PGON", idtmed[1108], barrel, 10); //Aluminium |
| 272 | gMC->Gspos("C3CO", 1, "C3MO", 0., 0., 0., 0, "ONLY"); |
| 273 | |
| 274 | barrel[5] = 580.; |
| 275 | barrel[6] = 585.; |
| 276 | barrel[8] = 580.; |
| 277 | barrel[9] = 585.; |
| 278 | gMC->Gsvolu("C3C1", "PGON", idtmed[1128], barrel, 10);// Aluminium |
| 279 | gMC->Gspos("C3C1", 1, "C3MO", 0., 0., 0., 0, "ONLY"); |
| 280 | |
| 281 | // Define yoke |
| 282 | |
| 283 | barrel[5] = 690.; |
| 284 | barrel[6] = 790.; |
| 285 | barrel[8] = 690.; |
| 286 | barrel[9] = 790.; |
| 287 | gMC->Gsvolu("C3YO", "PGON", idtmed[1109], barrel, 10); // Iron |
| 288 | gMC->Gspos("C3YO", 1, "C3MO", 0., 0., 0., 0, "ONLY"); |
| 289 | |
| 290 | |
| 291 | // Now create one inside the magnet as L3C1 |
| 292 | // voulme for tracking. |
| 293 | barrel[0] = 22.5; |
| 294 | barrel[1] = 45*magnetSides; |
| 295 | barrel[2] = magnetSides; |
| 296 | barrel[3] = planesPerpendicularToZ; |
| 297 | barrel[4] = -600.; |
| 298 | barrel[5] = 575.; |
| 299 | barrel[6] = 580.; |
| 300 | barrel[7] = 600.; |
| 301 | barrel[8] = 575.; |
| 302 | barrel[9] = 580.; |
| 303 | gMC->Gsvolu("C3CI", "PGON", idtmed[1134], barrel, 10); |
| 304 | gMC->Gspos("C3CI", 1 , "CRT", 0., -30., 0., 0, "ONLY"); |
| 305 | |
| 306 | // And a detector layer in the door 10 cm thick |
| 307 | // Volume for tracking. |
| 308 | barrel[0] = 22.5; |
| 309 | barrel[1] = 360.; |
| 310 | barrel[2] = 8.; |
| 311 | barrel[3] = 2.; |
| 312 | barrel[4] = 590.; |
| 313 | barrel[5] = 0.; |
| 314 | barrel[6] = 580.; |
| 315 | barrel[7] = 600.; |
| 316 | barrel[8] = barrel[5]; |
| 317 | barrel[9] = barrel[6]; |
| 318 | gMC->Gsvolu("C3C2", "PGON", idtmed[1154], barrel, 10); // Air |
| 319 | gMC->Gspos("C3C2", 1, "CRT", 0., -30., 0., 0, "ONLY"); |
| 320 | AliMatrix(idrotm[1010], 90., 0., 90., 90., 180., 0.); |
| 321 | gMC->Gspos("C3C2", 2, "CRT", 0., -30., 0., idrotm[1010], "ONLY"); |
| 322 | |
| 323 | |
| 324 | |
| 325 | barrel[4] = 600.; |
| 326 | barrel[5] = 0.; |
| 327 | barrel[6] = 790.; |
| 328 | barrel[7] = 700.; |
| 329 | barrel[8] = barrel[5]; |
| 330 | barrel[9] = barrel[6]; |
| 331 | gMC->Gsvolu("C3DO", "PGON", idtmed[1174], barrel, 10); // Air |
| 332 | gMC->Gspos("C3DO", 1, "CRT", 0., -30., 0., 0, "ONLY"); |
| 333 | AliMatrix(idrotm[1010], 90., 0., 90., 90., 180., 0.); |
| 334 | gMC->Gspos("C3DO", 2, "CRT", 0., -30., 0., idrotm[1010], "ONLY"); |
| 335 | |
| 336 | barrel[4] = 610.; |
| 337 | barrel[5] = 0.; |
| 338 | barrel[6] = 790.; |
| 339 | barrel[7] = 700.; |
| 340 | barrel[8] = barrel[5]; |
| 341 | barrel[9] = barrel[6]; |
| 342 | gMC->Gsvolu("C3FR", "PGON", idtmed[1149], barrel, 10); // Iron |
| 343 | gMC->Gspos("C3FR", 1, "C3DO", 0., 0., 0., 0, "ONLY"); |
| 344 | // INNER LAYER |
| 345 | |
| 346 | barrel[4] = 600.; |
| 347 | barrel[7] = 610.; |
| 348 | gMC->Gsvolu("C3IR", "PGON", idtmed[1149], barrel, 10); //Iron |
| 349 | gMC->Gspos("C3IR", 1, "C3DO", 0., 0., 0., 0, "ONLY"); |
| 350 | |
| 351 | } |
| 352 | |
| 353 | //_____________________________________________________________________________ |
| 354 | void AliCRTv1::CreateTPCGeometry() |
| 355 | { |
| 356 | cout<<"\n\n\tYou are requiring the CRT with the TPC Activated!\n\n"; |
| 357 | Int_t * idtmed = fIdtmed->GetArray() - 1099 ; |
| 358 | |
| 359 | // Disable the CRT StepManager method. |
| 360 | fCRTStatus = kFALSE; |
| 361 | // Disable the MAgnet |
| 362 | fMagnetStatus = kFALSE; |
| 363 | // Disable th RICH |
| 364 | fRICHStatus = kFALSE; |
| 365 | |
| 366 | // TPC |
| 367 | // Tpc SAndwich 1 - Al |
| 368 | // TSA1 |
| 369 | Float_t tube[5]; |
| 370 | tube[0]=274.8124; |
| 371 | tube[1]=278.; |
| 372 | tube[2]=252.1; |
| 373 | tube[3] = 0.; |
| 374 | tube[4] = 180.; |
| 375 | gMC->Gsvolu("CSA1","TUBS",idtmed[1154],tube,5); |
| 376 | // TSA1->TOCV (0.,0.,3.) ->TOIN (0.,0.,0.)->TPC (0.,0.,0.)->ALIC(0.,0.,0.) |
| 377 | gMC->Gspos("CSA1 ",1,"CRT",0.,0.,0.,0,"ONLY"); |
| 378 | |
| 379 | } |
| 380 | |
| 381 | //_____________________________________________________________________________ |
| 382 | void AliCRTv1::CreateRICHGeometry() |
| 383 | { |
| 384 | |
| 385 | cout<<"\n\n\tYou are requiring the CRT with the RICH Activated!\n\n"; |
| 386 | |
| 387 | Int_t idrotm[2499]; // The rotation matrix. |
| 388 | |
| 389 | Int_t * idtmed = fIdtmed->GetArray() - 1099 ; |
| 390 | |
| 391 | // Disable the CRT StepManager method. |
| 392 | fCRTStatus = kFALSE; |
| 393 | // Disable the MAgnet |
| 394 | fMagnetStatus = kFALSE; |
| 395 | |
| 396 | |
| 397 | // now create volume to simulate the HMPID volume. CSI |
| 398 | Float_t csi_length = 160*.8 + 2.6; |
| 399 | Float_t csi_width = 144*.84 + 2*2.6; |
| 400 | Float_t tbox[3]; |
| 401 | tbox[0] = csi_width/2; |
| 402 | tbox[1] = 11.5; |
| 403 | tbox[2] = csi_length/2; |
| 404 | gMC->Gsvolu("CRIC ", "BOX ", idtmed[1174], tbox, 3); |
| 405 | |
| 406 | Double_t dOffset = 490+1.267 - 8/2; // distance from center of mother volume ALIC to methane |
| 407 | |
| 408 | Double_t dAlpha = 19.5; // angle between centers of chambers - y-z plane |
| 409 | Double_t dAlphaRad = dAlpha*kDegrad; |
| 410 | |
| 411 | Double_t dBeta = 20.; // angle between center of chambers - y-x plane |
| 412 | Double_t dBetaRad = dBeta*kDegrad; |
| 413 | |
| 414 | Double_t dRotAngle = 60.; // the whole RICH is to be rotated in x-y plane + means clockwise rotation |
| 415 | Double_t dRotAngleRad = dRotAngle*kDegrad; |
| 416 | |
| 417 | |
| 418 | TRotMatrix *pRotMatrix; // tmp pointer |
| 419 | |
| 420 | TVector3 vector(0,dOffset,0); // Position of chamber 2 without rotation |
| 421 | |
| 422 | // Chamber 0 standalone (no other chambers in this row) |
| 423 | AliMatrix(idrotm[1000],90, -dRotAngle+360,90-dAlpha, 90-dRotAngle, dAlpha, -90+300); |
| 424 | pRotMatrix=new TRotMatrix("rot993","rot993",90,-dRotAngle, 90-dAlpha,90-dRotAngle,dAlpha, -90); |
| 425 | |
| 426 | vector.SetXYZ(0,dOffset,0); vector.RotateX(dAlphaRad); |
| 427 | vector.RotateZ(-dRotAngleRad); |
| 428 | |
| 429 | gMC->Gspos("CRIC",1,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1000], "ONLY"); |
| 430 | |
| 431 | // Chamber 1 |
| 432 | AliMatrix(idrotm[1001],90,-dBeta-dRotAngle,90,90-dBeta-dRotAngle, 0,0); |
| 433 | |
| 434 | pRotMatrix=new TRotMatrix("rot994","rot994",90,-dBeta-dRotAngle,90,90-dBeta-dRotAngle,0,0); |
| 435 | |
| 436 | vector.SetXYZ(0,dOffset,0); vector.RotateZ(-dBetaRad); |
| 437 | vector.RotateZ(-dRotAngleRad); |
| 438 | |
| 439 | gMC->Gspos("CRIC",2,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1001], "ONLY"); |
| 440 | |
| 441 | // Chamber 2 the top one with no Alpha-Beta rotation |
| 442 | AliMatrix(idrotm[1002],90,-dRotAngle,90,90-dRotAngle,0,0); |
| 443 | |
| 444 | pRotMatrix=new TRotMatrix("rot995","rot995",90,-dRotAngle,90,90-dRotAngle,0,0); |
| 445 | |
| 446 | vector.SetXYZ(0,dOffset,0); |
| 447 | vector.RotateZ(-dRotAngleRad); |
| 448 | |
| 449 | gMC->Gspos("CRIC",3,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1002], "ONLY"); |
| 450 | |
| 451 | // Chamber 3 |
| 452 | AliMatrix(idrotm[1003],90,dBeta-dRotAngle,90.,90+dBeta-dRotAngle,0,0); |
| 453 | pRotMatrix=new TRotMatrix("rot996","rot996", 90,dBeta-dRotAngle,90.,90+dBeta-dRotAngle,0,0); |
| 454 | |
| 455 | vector.SetXYZ(0,dOffset,0); vector.RotateZ(dBetaRad); |
| 456 | vector.RotateZ(-dRotAngleRad); |
| 457 | |
| 458 | gMC->Gspos("CRIC",4,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1003], "ONLY"); |
| 459 | |
| 460 | // Chamber 4 |
| 461 | AliMatrix(idrotm[1004],90,360-dBeta-dRotAngle,108.2,90-dBeta-dRotAngle,18.2,90-dBeta-60); |
| 462 | pRotMatrix=new TRotMatrix("rot997","rot997",90,360-dBeta-dRotAngle,108.2,90-dBeta-dRotAngle,18.2,90-dBeta); |
| 463 | |
| 464 | vector.SetXYZ(0,dOffset,0); vector.RotateZ(-dBetaRad); vector.RotateX(-dAlphaRad); |
| 465 | vector.RotateZ(-dRotAngleRad); |
| 466 | |
| 467 | gMC->Gspos("CRIC",5,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1004], "ONLY"); |
| 468 | |
| 469 | // Chamber 5 |
| 470 | AliMatrix(idrotm[1005],90,-dRotAngle+360,90+dAlpha,90-dRotAngle,dAlpha,90-60); |
| 471 | |
| 472 | pRotMatrix=new TRotMatrix("rot998","rot998",90,-dRotAngle,90+dAlpha,90-dRotAngle,dAlpha,90); |
| 473 | |
| 474 | vector.SetXYZ(0,dOffset,0); vector.RotateX(-dAlphaRad); |
| 475 | vector.RotateZ(-dRotAngleRad); |
| 476 | |
| 477 | gMC->Gspos("CRIC",6,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1005], "ONLY"); |
| 478 | |
| 479 | // Chamber 6 |
| 480 | AliMatrix(idrotm[1006],90,dBeta-dRotAngle+360,108.2,90+dBeta-dRotAngle,18.2,90+dBeta-60); |
| 481 | |
| 482 | pRotMatrix=new TRotMatrix("rot999","rot999",90,dBeta-dRotAngle,108.2,90+dBeta-dRotAngle,18.2,90+dBeta); |
| 483 | |
| 484 | vector.SetXYZ(0,dOffset,0); vector.RotateZ(dBetaRad); vector.RotateX(-dAlphaRad); |
| 485 | vector.RotateZ(-dRotAngleRad); |
| 486 | |
| 487 | gMC->Gspos("CRIC",7,"CRT",vector.X(),vector.Y(),vector.Z(),idrotm[1006], "ONLY"); |
| 488 | |
| 489 | } |
| 490 | |
| 491 | //_____________________________________________________________________________ |
| 492 | void AliCRTv1::CreateMolasse() |
| 493 | { |
| 494 | // |
| 495 | // |
| 496 | // |
| 497 | |
| 498 | Int_t idrotm[2499]; // The rotation matrix. |
| 499 | |
| 500 | Int_t * idtmed = fIdtmed->GetArray() - 1099 ; |
| 501 | |
| 502 | // |
| 503 | // Molasse |
| 504 | // |
| 505 | |
| 506 | // Exactly above the hall |
| 507 | Float_t tspar[5]; |
| 508 | tspar[0] = 1170.; |
| 509 | tspar[1] = 1170. + 375.; |
| 510 | tspar[2] = (1900.+1150.)/2.+100.; |
| 511 | tspar[3] = 0.; |
| 512 | tspar[4] = 180.; |
| 513 | gMC->Gsvolu("CMO1", "TUBS", idtmed[1123], tspar, 5); |
| 514 | gMC->Gspos("CMO1", 1, "CRT", 0., 500., 1900.-tspar[2]+400., 0, "MANY"); |
| 515 | |
| 516 | Float_t tbox[3]; |
| 517 | tbox[0] = 1250.; |
| 518 | tbox[1] = (4420. - 1670.)/2.; |
| 519 | tbox[2] = (1900.+1150.)/2. + 200.; |
| 520 | gMC->Gsvolu("CM12", "BOX", idtmed[1123], tbox, 3); |
| 521 | gMC->Gspos("CM12", 1, "CRT", 0., 4420. -tbox[1], 1900.-tbox[2]+400., 0, "MANY"); |
| 522 | |
| 523 | AliMatrix(idrotm[2003], 0., 0., 90., 0., 90., 90.); |
| 524 | // Along the PM25 |
| 525 | Float_t tube[3]; |
| 526 | tube[0] = 455. + 100.; |
| 527 | tube[1] = 555. + 375.; |
| 528 | tube[2] = (5150. - 1166.)/2.; |
| 529 | gMC->Gsvolu("CMO2", "TUBE", idtmed[1123], tube, 3); |
| 530 | gMC->Gspos("CMO2", 1, "CRT", -2100., 4420.-tube[2], 0., idrotm[2003], "MANY"); |
| 531 | |
| 532 | |
| 533 | // Along the PGC2 |
| 534 | tube[0] = 650.; |
| 535 | tube[1] = 2987.7; |
| 536 | tube[2] = (5150. - 690.)/2.; |
| 537 | gMC->Gsvolu("CMO3", "TUBE", idtmed[1123], tube, 3); |
| 538 | gMC->Gspos("CMO3", 1, "CRT", 375., 4420.-tube[2], 1900.+2987.7, idrotm[2003], "MANY"); |
| 539 | // Behind the PGC2 up to the end of the M. volume. |
| 540 | tbox[0] = 12073.; |
| 541 | tbox[1] = 2575. + 95.; |
| 542 | tbox[2] = (12073. - 1900.-2987.7-650.)/2.; |
| 543 | gMC->Gsvolu("CMO7", "BOX", idtmed[1123], tbox, 3); |
| 544 | gMC->Gspos("CMO7", 1, "CRT", 0., 4420.-tbox[1], 1900.+2987.7+650.+tbox[2], 0, "MANY"); |
| 545 | |
| 546 | // Along the PX24 , upper part. |
| 547 | tube[0] = 1250.; |
| 548 | tube[1] = 2300; |
| 549 | tube[2] = 2575. - 1300. + 95.; |
| 550 | gMC->Gsvolu("CMO4", "TUBE", idtmed[1123], tube, 3); |
| 551 | gMC->Gspos("CMO4", 1, "CRT", 0., 404.+1300.+tube[2], -2300., idrotm[2003], "MANY"); |
| 552 | |
| 553 | // Along the PX24 , lower part |
| 554 | tspar[0] = 1250.; |
| 555 | tspar[1] = 2300; |
| 556 | tspar[2] = 1300.; |
| 557 | tspar[3] = kRaddeg*TMath::ASin(1070./1150.); |
| 558 | tspar[4] = 360. - tspar[3]; |
| 559 | gMC->Gsvolu("CMO5", "TUBS", idtmed[1123], tspar, 5); |
| 560 | gMC->Gspos("CMO5", 1, "CRT", 0., 404., -2300., idrotm[2003], "MANY"); |
| 561 | // behind the PX24 |
| 562 | tbox[0] = 12073.; |
| 563 | tbox[1] = 2575. + 95.; |
| 564 | tbox[2] = 8523./2.; |
| 565 | gMC->Gsvolu("CMO6", "BOX", idtmed[1123], tbox, 3); |
| 566 | gMC->Gspos("CMO6", 1, "CRT", 0., 4420.-tbox[1], -3550.-tbox[2], 0, "MANY"); |
| 567 | |
| 568 | |
| 569 | // On the right side of th hall |
| 570 | tbox[0] = (12073. - 1250.)/2.; |
| 571 | tbox[1] = 2575. + 95.; |
| 572 | tbox[2] = (8437.7+650.)/2.; |
| 573 | gMC->Gsvolu("CMO8", "BOX", idtmed[1123], tbox, 3); |
| 574 | gMC->Gspos("CMO8", 1, "CRT", 1250.+tbox[0], 4420.-tbox[1], -3550.+tbox[2], 0, "MANY"); |
| 575 | |
| 576 | // on the left side of the hall, behind |
| 577 | tbox[0] = (12073. - 2755.)/2.; |
| 578 | tbox[1] = 2575. + 95.; |
| 579 | tbox[2] = (8437.7+650.)/2.; |
| 580 | gMC->Gsvolu("CMO9", "BOX", idtmed[1123], tbox, 3); |
| 581 | gMC->Gspos("CMO9", 1, "CRT", -2755.-tbox[0], 4420.-tbox[1], -3550.+tbox[2], 0, "MANY"); |
| 582 | |
| 583 | |
| 584 | // Molasse betwen the PX24 & PM25 on the left side. |
| 585 | tbox[0] = (2755. - 1250.)/2.; |
| 586 | tbox[1] = 2575. + 95.; |
| 587 | tbox[2] = (3550. - 555.)/2.; |
| 588 | gMC->Gsvolu("CM10", "BOX", idtmed[1123], tbox, 3); |
| 589 | gMC->Gspos("CM10", 1, "CRT", -1250.-tbox[0], 4420.-tbox[1], -tbox[2]-555., 0, "MANY"); |
| 590 | |
| 591 | |
| 592 | // Molasse betwen the PGC2 & PM25 on the left side. |
| 593 | tbox[0] = (2755. - 1250.)/2.; |
| 594 | tbox[1] = 2575. + 95.; |
| 595 | tbox[2] = (1900.+2987.7 - 555. + 650.)/2.; |
| 596 | gMC->Gsvolu("CM11", "BOX", idtmed[1123], tbox, 3); |
| 597 | gMC->Gspos("CM11", 1, "CRT", -1250.-tbox[0], 4420.-tbox[1], 555.+tbox[2], 0, "MANY"); |
| 598 | |
| 599 | |
| 600 | } |
| 601 | |
| 602 | //_____________________________________________________________________________ |
| 603 | void AliCRTv1::CreateShafts() |
| 604 | { |
| 605 | // |
| 606 | // |
| 607 | // |
| 608 | Int_t idrotm[2499]; // The rotation matrix. |
| 609 | |
| 610 | Int_t * idtmed = fIdtmed->GetArray() - 1099 ; |
| 611 | |
| 612 | // Create a mother volume. |
| 613 | Float_t pbox[3]; |
| 614 | //pbox[0] = AliCRTConstants::fgDepth*TMath::Tan(67.5*kDegrad); |
| 615 | pbox[0] = 12073.; |
| 616 | pbox[1] = AliCRTConstants::fgDepth; |
| 617 | pbox[2] = pbox[0]; |
| 618 | gMC->Gsvolu("CRT", "BOX", idtmed[1114], pbox, 3); |
| 619 | gMC->Gspos("CRT", 1, "ALIC", 0., 0., 0., 0, "ONLY"); |
| 620 | |
| 621 | // HAll ceiling |
| 622 | Float_t ptubs[5]; |
| 623 | ptubs[0] = 1070.; |
| 624 | ptubs[1] = 1170.; |
| 625 | ptubs[2] = 1900.; |
| 626 | ptubs[3] = 0.; |
| 627 | ptubs[4] = 180.; |
| 628 | gMC->Gsvolu("CHC1", "TUBS", idtmed[1116], ptubs, 5); |
| 629 | gMC->Gspos("CHC1", 1, "CRT", 0., 500., 0., 0, "ONLY"); |
| 630 | |
| 631 | |
| 632 | // |
| 633 | // Acces shafts |
| 634 | // |
| 635 | AliMatrix(idrotm[2001], 0., 0., 90., 0., 90., 90.); |
| 636 | |
| 637 | // PX24 |
| 638 | ptubs[0] = 1150.; |
| 639 | ptubs[1] = 1250.; |
| 640 | ptubs[2] = 1300.; |
| 641 | ptubs[3] = kRaddeg*TMath::ASin(1070./ptubs[0]); |
| 642 | ptubs[4] = 360 - ptubs[3]; |
| 643 | gMC->Gsvolu("CSF1", "TUBS", idtmed[1116], ptubs, 5); |
| 644 | gMC->Gspos("CSF1", 1, "CRT", 0., 404., -2300., idrotm[2001], "MANY"); |
| 645 | |
| 646 | Float_t ptube[3]; |
| 647 | ptube[0] = ptubs[0]; |
| 648 | ptube[1] = ptubs[1]; |
| 649 | ptube[2] = 2575. - ptubs[2] + 95.; |
| 650 | gMC->Gsvolu("CSF2", "TUBE", idtmed[1116], ptube, 3); |
| 651 | gMC->Gspos("CSF2", 1, "CRT", 0., 404.+ptubs[2]+ptube[2], -2300., idrotm[2001], "MANY"); |
| 652 | |
| 653 | // Concrete walls along the shaft |
| 654 | pbox[0] = 585./2.; |
| 655 | pbox[1] = 2575. + 95.; |
| 656 | pbox[2] = 20.; |
| 657 | gMC->Gsvolu("CSW1", "BOX", idtmed[1116], pbox, 3); |
| 658 | gMC->Gspos("CSW1", 1, "CRT", -290-pbox[0], 404.-1300.+pbox[1], -3450.+210.*2, 0, "MANY"); |
| 659 | |
| 660 | // |
| 661 | pbox[0] = 750./2.; |
| 662 | pbox[1] = 2575. + 95.; |
| 663 | pbox[2] = 20.; |
| 664 | gMC->Gsvolu("CSW3", "BOX", idtmed[1116], pbox, 3); |
| 665 | gMC->Gspos("CSW3", 1, "CRT", 420.-290.+pbox[0], 404.-1300.+pbox[1], -3450.+210.*2, 0, "MANY"); |
| 666 | |
| 667 | // |
| 668 | pbox[0] = 60.; |
| 669 | pbox[1] = 2575. + 95.; |
| 670 | pbox[2] = 210.; |
| 671 | gMC->Gsvolu("CSW2", "BOX", idtmed[1116], pbox, 3); |
| 672 | gMC->Gspos("CSW2", 1, "CRT", -290-pbox[0], 404.-1300.+pbox[1], -3450.+pbox[2], 0, "MANY"); |
| 673 | gMC->Gspos("CSW2", 2, "CRT", 420.-290.+pbox[0], 404.-1300.+pbox[1], -3450.+pbox[2], 0, "MANY"); |
| 674 | |
| 675 | |
| 676 | // |
| 677 | pbox[0] = 1000.; |
| 678 | pbox[1] = 80.; |
| 679 | pbox[2] = 200.; |
| 680 | gMC->Gsvolu("CSP1", "BOX", idtmed[1116], pbox, 3); |
| 681 | gMC->Gspos("CSP1", 1, "CRT", 0., 2600.-700., -1150-pbox[2], 0, "MANY"); |
| 682 | |
| 683 | // |
| 684 | pbox[0] = 340.8; |
| 685 | pbox[1] = 300./2.; |
| 686 | pbox[2] = 460./2.; |
| 687 | gMC->Gsvolu("CSP2", "BOX", idtmed[1116], pbox, 3); |
| 688 | gMC->Gspos("CSP2", 1, "CRT", 0., 2950.-700., -3450+pbox[2], 0, "MANY"); |
| 689 | |
| 690 | // |
| 691 | pbox[0] = 600.; |
| 692 | pbox[1] = 150.; |
| 693 | pbox[2] = 75.; |
| 694 | gMC->Gsvolu("CSP3", "BOX", idtmed[1116], pbox, 3); |
| 695 | gMC->Gspos("CSP3", 1, "CRT", 0., 2950.-700., -1150.-210.-pbox[2], 0, "MANY"); |
| 696 | |
| 697 | // |
| 698 | pbox[0] = 600.; |
| 699 | pbox[1] = 250.; |
| 700 | pbox[2] = 38.; |
| 701 | gMC->Gsvolu("CSP4", "BOX", idtmed[1116], pbox, 3); |
| 702 | gMC->Gspos("CSP4", 1, "CRT", 0., 2950.-700.+155.+pbox[1], -1150.-210.-pbox[2], 0, "MANY"); |
| 703 | |
| 704 | |
| 705 | // Shielding plug |
| 706 | pbox[0] = 850.; |
| 707 | pbox[1] = 90.; |
| 708 | pbox[2] = 720.; |
| 709 | gMC->Gsvolu("CSP5", "BOX", idtmed[1116], pbox, 3); |
| 710 | gMC->Gspos("CSP5", 1, "CRT", 0., 2950.-700., -3450.+460.+pbox[2], 0, "MANY"); |
| 711 | |
| 712 | // |
| 713 | pbox[0] = 80.; |
| 714 | pbox[1] = 150.; |
| 715 | pbox[2] = 720.; |
| 716 | gMC->Gsvolu("CSP6", "BOX", idtmed[1116], pbox, 3); |
| 717 | gMC->Gspos("CSP6", 1, "CRT", 1150.-600., 2950.-700., -3450.+460.+pbox[2], 0, "MANY"); |
| 718 | gMC->Gspos("CSP6", 2, "CRT", -1150.+600., 2950.-700., -3450.+460.+pbox[2], 0, "MANY"); |
| 719 | |
| 720 | |
| 721 | // |
| 722 | pbox[0] = 130.; |
| 723 | pbox[1] = 60.; |
| 724 | pbox[2] = 750.; |
| 725 | gMC->Gsvolu("CSP7", "BOX", idtmed[1116], pbox, 3); |
| 726 | gMC->Gspos("CSP7", 1, "CRT", 850.+pbox[0], 2950.-700.+100., -3450.+460.+pbox[2], 0, "MANY"); |
| 727 | gMC->Gspos("CSP7", 2, "CRT", -850.-pbox[0], 2950.-700.+100., -3450.+460.+pbox[2], 0, "MANY"); |
| 728 | |
| 729 | |
| 730 | // PM25 Acces Shaft |
| 731 | ptube[0] = 910./2.; |
| 732 | ptube[1] = ptube[0] + 100.; |
| 733 | ptube[2] = (5150. - 1166.)/2.; |
| 734 | gMC->Gsvolu("CSF3", "TUBE", idtmed[1116], ptube, 3); |
| 735 | gMC->Gspos("CSF3", 1, "CRT", -2100., AliCRTConstants::fgDepth-ptube[2], 0., idrotm[2001], "MANY"); |
| 736 | |
| 737 | // PGC2 Access Shaft |
| 738 | ptube[0] = 1100./2.; |
| 739 | ptube[1] = ptube[0] + 100.; |
| 740 | ptube[2] = (5150. - 690.)/2.; |
| 741 | gMC->Gsvolu("CSF4", "TUBE", idtmed[1116], ptube, 3); |
| 742 | gMC->Gspos("CSF4", 1, "CRT", 375., AliCRTConstants::fgDepth-ptube[2], 1900.+2987.7, idrotm[2001], "MANY"); |
| 743 | |
| 744 | } |
| 745 | |
| 746 | //_____________________________________________________________________________ |
| 747 | void AliCRTv1::DrawDetector() |
| 748 | { |
| 749 | // |
| 750 | // Draw a shaded view of the L3 magnet |
| 751 | // |
| 752 | cout << "AliCRTv1::DrawModule() : Drawing the module" << endl; |
| 753 | |
| 754 | |
| 755 | Int_t able = 1; |
| 756 | Int_t enable = 0; |
| 757 | gMC->Gsatt("*", "seen", -1); |
| 758 | gMC->Gsatt("alic", "seen", 0); |
| 759 | |
| 760 | gMC->Gsatt("ALIC","seen",enable); |
| 761 | gMC->Gsatt("CRT", "seen",enable); |
| 762 | gMC->Gsatt("L3MO","seen", able); // L3 Magnet |
| 763 | //gMC->Gsatt("CRT1","seen", able); // Scintillators |
| 764 | gMC->Gsatt("CRT4","seen", able); // Scintillators barrel |
| 765 | |
| 766 | // Draw the molasse volumes |
| 767 | gMC->Gsatt("CMO1","seen",enable); // Exactly above the HALL |
| 768 | gMC->Gsatt("CMO2","seen",enable); // Molasse, along the PM25 |
| 769 | gMC->Gsatt("CMO3","seen",enable); // molasse along the PGC2 |
| 770 | gMC->Gsatt("CMO4","seen",enable); // Molasse, behind the PX24 upper part |
| 771 | gMC->Gsatt("CMO5","seen",enable); // molasse behind px24, lower part |
| 772 | gMC->Gsatt("CMO6","seen",enable); // behind the PX24 |
| 773 | gMC->Gsatt("CMO7","seen",enable); // behind the PGC2 |
| 774 | gMC->Gsatt("CMO8","seen",enable); // on the right side. |
| 775 | gMC->Gsatt("CMO9","seen",enable); // on the left side. |
| 776 | gMC->Gsatt("CM10","seen",enable); // betwen PX24 & PM25. |
| 777 | gMC->Gsatt("CM11","seen",enable); // betwen PGC2 & PM25. |
| 778 | gMC->Gsatt("CM12","seen",enable); // box above the hall. |
| 779 | |
| 780 | gMC->Gdopt("hide", "on"); |
| 781 | gMC->Gdopt("edge","off"); |
| 782 | gMC->Gdopt("shad", "on"); |
| 783 | gMC->Gsatt("*", "fill", 7); |
| 784 | gMC->SetClipBox("ALIC", 0, 3000, -3000, 3000, -6000, 6000); |
| 785 | gMC->DefaultRange(); |
| 786 | gMC->Gdraw("alic", 40, 30, 0, 10, 9.5, .009, .009); |
| 787 | gMC->Gdhead(1111, "View of CRT(ACORDE)"); |
| 788 | gMC->Gdman(18, 4, "MAN"); |
| 789 | |
| 790 | |
| 791 | } |
| 792 | |
| 793 | //_____________________________________________________________________________ |
| 794 | void AliCRTv1::Init() |
| 795 | { |
| 796 | // |
| 797 | // Initialise L3 magnet after it has been built |
| 798 | Int_t i; |
| 799 | // |
| 800 | if(fDebug) { |
| 801 | printf("\n%s: ",ClassName()); |
| 802 | for(i=0;i<35;i++) printf("*"); |
| 803 | printf(" CRTv1_INIT "); |
| 804 | for(i=0;i<35;i++) printf("*"); |
| 805 | printf("\n%s: ",ClassName()); |
| 806 | // |
| 807 | // Here the CRTv1 initialisation code (if any!) |
| 808 | for(i=0;i<80;i++) printf("*"); |
| 809 | printf("\n"); |
| 810 | } |
| 811 | |
| 812 | } |
| 813 | |
| 814 | //____________________________________________________________________________ |
| 815 | void AliCRTv1::StepManager() |
| 816 | { |
| 817 | // |
| 818 | // Called for every step in the Cosmic Ray Trigger |
| 819 | // |
| 820 | static Int_t vol[5]; |
| 821 | Int_t ipart; |
| 822 | TLorentzVector pos; |
| 823 | TLorentzVector mom; |
| 824 | |
| 825 | static Float_t hits[14]; |
| 826 | static Float_t eloss; |
| 827 | static Float_t elossMag; |
| 828 | |
| 829 | if ( !gMC->IsTrackAlive() ) return; |
| 830 | |
| 831 | if (gMC->IsNewTrack()) { |
| 832 | // Reset the deposited energy |
| 833 | eloss = 0.; |
| 834 | elossMag = 0.; |
| 835 | } |
| 836 | |
| 837 | // Add th energy loss in each step. |
| 838 | eloss += gMC->Edep(); |
| 839 | |
| 840 | gMC->TrackPosition(pos); |
| 841 | |
| 842 | // |
| 843 | // CRT |
| 844 | // |
| 845 | |
| 846 | if ( gMC->IsTrackEntering() && (strcmp(gMC->CurrentVolName(),"CRT4") == 0) |
| 847 | &&(gMC->TrackPid() == kMuonMinus || gMC->TrackPid() == kMuonPlus) ) { |
| 848 | |
| 849 | // Get current particle id(ipart),track position (pos) and momentum (mom) |
| 850 | gMC->TrackPosition(pos); |
| 851 | gMC->TrackMomentum(mom); |
| 852 | ipart = gMC->TrackPid(); |
| 853 | |
| 854 | vol[0] = 1; |
| 855 | vol[1] = 0; |
| 856 | vol[2] = 0; |
| 857 | vol[3] = 0; |
| 858 | vol[4] = 0; |
| 859 | |
| 860 | ipart = gMC->TrackPid(); |
| 861 | hits[0] = (Float_t)ipart; // (fId) |
| 862 | |
| 863 | hits[1] = pos[0]; // X coordinate (fX) |
| 864 | hits[2] = pos[1]; // Y coordinate (fY) |
| 865 | hits[3] = pos[2]; // Z coordinate (fZ) |
| 866 | hits[4] = mom[0]; // Px (fpxug) |
| 867 | hits[5] = mom[1]; // Py (fpyug) |
| 868 | hits[6] = mom[2]; // Pz (fpzug) |
| 869 | |
| 870 | hits[7] = gMC->GetMedium(); //layer(flay) |
| 871 | hits[8] = eloss; // Energy loss |
| 872 | |
| 873 | hits[9] = 1; // CRT mother activated. |
| 874 | hits[10] = 0; |
| 875 | hits[11] = 0; |
| 876 | hits[12] = 0; |
| 877 | hits[13] = 0; |
| 878 | |
| 642f15cf |
879 | //hits[9] = gAlice->GetCurrentTrackNumber(); |
| 53d2a7cd |
880 | |
| 642f15cf |
881 | AddHit(gAlice->GetCurrentTrackNumber(),vol, hits); |
| 53d2a7cd |
882 | |
| 883 | eloss = 0.; |
| 884 | |
| 885 | } else if (gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"CRT1")==0) |
| 886 | &&(gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus)) { |
| 887 | |
| 888 | vol[0] = 0; |
| 889 | vol[1] = 1; |
| 890 | vol[2] = 0; |
| 891 | vol[3] = 0; |
| 892 | vol[4] = 0; |
| 893 | |
| 894 | hits[9] = 0; // CRT mother activated. |
| 895 | hits[10] = 1; |
| 896 | hits[11] = 0; |
| 897 | hits[12] = 0; |
| 898 | hits[13] = 0; |
| 899 | |
| 900 | //hits[10] = 1; |
| 901 | |
| 642f15cf |
902 | //AddHit(gAlice->GetCurrentTrackNumber(),vol, hits); |
| 53d2a7cd |
903 | |
| 904 | //eloss = 0.; |
| 905 | |
| 906 | |
| 907 | } else if (gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"C3CI")==0) |
| 908 | &&(gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus)) { |
| 909 | |
| 910 | // |
| 911 | // Inside the magnet, upper part. |
| 912 | // |
| 913 | |
| 914 | // Get current particle id(ipart),track position (pos) and momentum (mom) |
| 915 | |
| 916 | vol[0] = 0; |
| 917 | vol[1] = 0; |
| 918 | vol[2] = 1; |
| 919 | vol[3] = 0; |
| 920 | vol[4] = 0; |
| 921 | |
| 922 | hits[9] = 0; // CRT mother activated. |
| 923 | hits[10] = 0; |
| 924 | hits[11] = 1; |
| 925 | hits[12] = 0; |
| 926 | hits[13] = 0; |
| 927 | |
| 642f15cf |
928 | AddHit(gAlice->GetCurrentTrackNumber(),vol, hits); |
| 53d2a7cd |
929 | |
| 930 | //eloss = 0.; |
| 931 | |
| 932 | } else if ( gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"CRIC")==0) |
| 933 | && (gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus) ) { |
| 934 | |
| 935 | // |
| 936 | // HMPID |
| 937 | // |
| 938 | |
| 939 | // Get current particle id(ipart),track position (pos) and momentum (mom) |
| 940 | |
| 941 | vol[0] = 0; |
| 942 | vol[1] = 0; |
| 943 | vol[2] = 0; |
| 944 | vol[3] = 1; |
| 945 | vol[4] = 0; |
| 946 | |
| 947 | hits[9] = 0; |
| 948 | hits[10] = 0; |
| 949 | hits[11] = 0; |
| 950 | hits[12] = 1; |
| 951 | hits[13] = 0; |
| 952 | |
| 642f15cf |
953 | AddHit(gAlice->GetCurrentTrackNumber(),vol, hits); |
| 53d2a7cd |
954 | |
| 955 | //eloss = 0.; |
| 956 | |
| 957 | |
| 958 | } else if (gMC->IsTrackEntering()&&(strcmp(gMC->CurrentVolName(),"CSA1")==0) |
| 959 | &&(gMC->TrackPid()==kMuonMinus || gMC->TrackPid()==kMuonPlus)) { |
| 960 | |
| 961 | // |
| 962 | // TPC |
| 963 | // |
| 964 | |
| 965 | // Get current particle id(ipart),track position (pos) and momentum (mom) |
| 966 | |
| 967 | vol[0] = 0; |
| 968 | vol[1] = 0; |
| 969 | vol[2] = 0; |
| 970 | vol[3] = 0; |
| 971 | vol[4] = 1; |
| 972 | |
| 973 | hits[9] = 0; |
| 974 | hits[10] = 0; |
| 975 | hits[11] = 0; |
| 976 | hits[12] = 0; |
| 977 | hits[13] = 1; |
| 978 | |
| 979 | |
| 642f15cf |
980 | AddHit(gAlice->GetCurrentTrackNumber(),vol, hits); |
| 53d2a7cd |
981 | |
| 982 | //eloss = 0.; |
| 983 | |
| 984 | } else { |
| 985 | return; |
| 986 | } |
| 987 | |
| 988 | |
| 989 | } |
git://git.uio.no / u / mrichter / AliRoot.git / blame