PHOS/AliPHOSv2.cxx

   1 //-*-C++-*-
   2 //_________________________________________________________________________
   3 // Manager and hit classes for PHOS
   4 //*-- Author : Maxim Volkov, RRC KI
   5 // AliPHOSv2 derives directly from AliDetector, because too much functionality
   6 // has been put in AliPHOS for my liking.
   7 //////////////////////////////////////////////////////////////////////////////
   8
   9 // --- ROOT system ---
  10 #include "TH1.h"
  11 #include "TRandom.h"
  12 #include "TFile.h"
  13 #include "TTree.h"
  14 #include "TBRIK.h"
  15 #include "TNode.h"
  16
  17 // --- Standard library ---
  18 #include <stdio.h>
  19 #include <string.h>
  20 #include <stdlib.h>
  21
  22 // --- galice header files ---
  23 #include "AliPHOSv2.h"
  24 #include "AliRun.h"
  25 #include "AliConst.h"
  26 #include "AliMC.h"
  27
  28 ///////////////////////////////////////////////////////////////////////////////
  29
  30 ClassImp(AliPHOSv2)
  31
  32 ///////////////////////////////////////////////////////////////////////////////
  33
  34 AliPHOSv2::~AliPHOSv2(void)
  35 {
  36
  37   //  fNV       = 0;
  38   //  fNH       = 0;
  39   fIshunt   = 0;
  40
  41   delete fHits;
  42
  43 }
  44
  45 ///////////////////////////////////////////////////////////////////////////////
  46
  47 AliPHOSv2::AliPHOSv2()
  48 {
  49
  50   //  fNV       = 0;
  51   //  fNH       = 0;
  52   fIshunt   = 0;
  53   fHits     = 0;
  54   fDigits   = 0;
  55
  56 }
  57
  58 ///////////////////////////////////////////////////////////////////////////////
  59
  60 AliPHOSv2::AliPHOSv2(const char *name, const char *title):
  61   AliDetector(name,title)
  62 {
  63
  64   //Begin_Html
  65   /*
  66     <img src="gif/aliphos.gif">
  67   */
  68   //End_Html
  69
  70   fHits   = new TClonesArray("AliPHOShitv2",  405);
  71
  72   //  fNV         =  4;
  73   //  fNH         =  4;
  74   fIshunt     =  1; // All hits are associated with primary particles
  75
  76   DefPars(); // Set geometry parameters
  77
  78   //  SetMarkerColor(kGreen);
  79   //  SetMarkerStyle(2);
  80   //  SetMarkerSize(0.4);
  81
  82 }
  83
  84 //////////////////////////////////////////////////////////////////////////////
  85
  86 void AliPHOSv2::DefPars(void)
  87 {
  88
  89   // Initialization of GEANT3 geometry parameters
  90   fXtlSize[0]=2.2;
  91   fXtlSize[1]=18.0;
  92   fXtlSize[2]=2.2;
  93
  94   fWrapThickness=0.01;
  95
  96   fPINSize[0]=1.0;
  97   fPINSize[1]=0.1;
  98   fPINSize[2]=1.0;
  99
 100   fCPVThickness=1.0;
 101
 102   fPHOSFoam[0]=214.6;
 103   fPHOSFoam[1]=80.0;
 104   fPHOSFoam[2]=260.0;
 105
 106   fPHOStxwall[0]=209.0;
 107   fPHOStxwall[1]=71.0;
 108   fPHOStxwall[2]=250.0;
 109
 110   fPHOSAir[0]=206.0;
 111   fPHOSAir[1]=66.0;
 112   fPHOSAir[2]=244.0;
 113
 114   fRadius[0]=447.0;
 115   fRadius[1]=460.0;
 116
 117   fPHOSextra[0]=0.005; // Titanium cover thickness
 118   fPHOSextra[1]=6.95; // Crystal support height
 119   fPHOSextra[2]=4.0; // Thermo Insulating outer cover Upper plate thickness
 120   fPHOSextra[3]=5.0; // Upper Polystyrene Foam plate thickness
 121   fPHOSextra[4]=2.0; // Thermo insulating Crystal Block wall thickness
 122   fPHOSextra[5]=0.06; // Upper Cooling Plate thickness
 123   fPHOSextra[6]=10.0; // Al Support Plate thickness
 124   fPHOSextra[7]=3.0; // Lower Thermo Insulating Plate thickness
 125   fPHOSextra[8]=1.0; // Lower Textolit Plate thickness
 126   fPHOSextra[9]=0.03; // 1/2 total gap between adjacent crystals
 127
 128   fNphi=88;
 129   fNz=104;
 130
 131   fNModules=4;
 132
 133   fPHOSAngle[0]=0.0; // Module position angles are set in CreateGeometry()
 134   fPHOSAngle[1]=0.0;
 135   fPHOSAngle[2]=0.0;
 136   fPHOSAngle[3]=0.0;
 137
 138 }
 139
 140 //////////////////////////////////////////////////////////////////////////////
 141
 142 void AliPHOSv2::Init(void)
 143 {
 144
 145   Int_t i;
 146
 147   printf("\n");
 148   for(i=0;i<35;i++) printf("*");
 149   printf(" PHOS_INIT ");
 150   for(i=0;i<35;i++) printf("*");
 151   printf("\n");
 152
 153   // Here the PHOS initialisation code (if any!)
 154   for(i=0;i<80;i++) printf("*");
 155   printf("\n");
 156
 157 }
 158
 159 //////////////////////////////////////////////////////////////////////////////
 160
 161 void AliPHOSv2::BuildGeometry()
 162 {
 163
 164   // Stolen completely from A. Zvyagine
 165
 166   TNode *Node, *Top;
 167
 168   const int kColorPHOS = kRed;
 169
 170   Top=gAlice->GetGeometry()->GetNode("alice");
 171
 172   // PHOS
 173   Float_t pphi=12.9399462;
 174   new TRotMatrix("rot988","rot988",90,-3*pphi,90,90-3*pphi,0,0);
 175   new TRotMatrix("rot989","rot989",90,-  pphi,90,90-  pphi,0,0);
 176   new TRotMatrix("rot990","rot990",90,   pphi,90,90+  pphi,0,0);
 177   new TRotMatrix("rot991","rot991",90, 3*pphi,90,90+3*pphi,0,0);
 178   new TBRIK("S_PHOS","PHOS box","void",107.3,40,130);
 179   Top->cd();
 180   Node=new TNode("PHOS1","PHOS1","S_PHOS",-317.824921,-395.014343,0,"rot988");
 181   Node->SetLineColor(kColorPHOS);
 182   fNodes->Add(Node);
 183   Top->cd();
 184   Node=new TNode("PHOS2","PHOS2","S_PHOS",-113.532333,-494.124908,0,"rot989");
 185   fNodes->Add(Node);
 186   Node->SetLineColor(kColorPHOS);
 187   Top->cd();
 188   Node=new TNode("PHOS3","PHOS3","S_PHOS", 113.532333,-494.124908,0,"rot990");
 189   Node->SetLineColor(kColorPHOS);
 190   fNodes->Add(Node);
 191   Top->cd();
 192   Node=new TNode("PHOS4","PHOS4","S_PHOS", 317.824921,-395.014343,0,"rot991");
 193   Node->SetLineColor(kColorPHOS);
 194   fNodes->Add(Node);
 195
 196 }
 197
 198 //////////////////////////////////////////////////////////////////////////////
 199
 200 void AliPHOSv2::CreateMaterials()
 201 {
 202
 203   // DEFINITION OF AVAILABLE PHOS MATERIALS
 204
 205   AliMC* pMC=AliMC::GetMC();
 206
 207   Int_t   ISXFLD=gAlice->Field()->Integ();
 208   Float_t SXMGMX=gAlice->Field()->Max();
 209
 210   // --- The PbWO4 crystals ---
 211   Float_t AX[3]={207.19, 183.85, 16.0};
 212   Float_t ZX[3]={82.0, 74.0, 8.0};
 213   Float_t WX[3]={1.0, 1.0, 4.0};
 214   Float_t DX=8.28;
 215
 216   // --- Titanium ---
 217   Float_t ATIT[3]={47.88, 26.98, 54.94};
 218   Float_t ZTIT[3]={22.0, 13.0, 25.0};
 219   Float_t WTIT[3]={69.0, 6.0, 1.0};
 220   Float_t DTIT=4.5;
 221
 222   // --- The polysterene scintillator (CH) ---
 223   Float_t AP[2]={12.011, 1.00794};
 224   Float_t ZP[2]={6.0, 1.0};
 225   Float_t WP[2]={1.0, 1.0};
 226   Float_t DP=1.032;
 227
 228   // --- Tyvek (CnH2n) ---
 229   Float_t AT[2]={12.011, 1.00794};
 230   Float_t ZT[2]={6.0, 1.0};
 231   Float_t WT[2]={1.0, 2.0};
 232   Float_t DT=0.331;
 233
 234   // --- Polystyrene foam ---
 235   Float_t AF[2]={12.011, 1.00794};
 236   Float_t ZF[2]={6.0, 1.0};
 237   Float_t WF[2]={1.0, 1.0};
 238   Float_t DF=0.12;
 239
 240   // --- Foam thermo insulation ---
 241   Float_t ATI[2]={12.011, 1.00794};
 242   Float_t ZTI[2]={6.0, 1.0};
 243   Float_t WTI[2]={1.0, 1.0};
 244   Float_t DTI=0.1;
 245
 246   // --- Textolit ---
 247   Float_t ATX[4]={16.0, 28.09, 12.011, 1.00794};
 248   Float_t ZTX[4]={8.0, 14.0, 6.0, 1.0};
 249   Float_t WTX[4]={292.0, 68.0, 462.0, 736.0};
 250   Float_t DTX=1.75;
 251
 252   Int_t *idtmed = gAlice->Idtmed();
 253
 254   AliMixture(0, "PbWO4$", AX, ZX, DX, -3, WX);
 255   AliMixture(1, "Polystyrene$", AP, ZP, DP, -2, WP);
 256   AliMaterial(2, "Al$", 26.98, 13., 2.7, 8.9, 999., 0, 0);
 257   // ---                     Absorption length^ is ignored ---
 258   AliMixture(3, "Tyvek$", AT, ZT, DT, -2, WT);
 259   AliMixture(4, "Foam$", AF, ZF, DF, -2, WF);
 260   AliMixture(5, "Titanium$", ATIT, ZTIT, DTIT, -3, WTIT);
 261   AliMaterial(6, "Si$", 28.09, 14., 2.33, 9.36, 42.3, 0, 0);
 262   AliMixture(7, "Thermo Insul.$", ATI, ZTI, DTI, -2, WTI);
 263   AliMixture(8, "Textolit$", ATX, ZTX, DTX, -4, WTX);
 264   AliMaterial(99, "Air$", 14.61, 7.3, 0.001205, 30420., 67500., 0, 0);
 265
 266   AliMedium(700, "PHOS Xtal    $", 0, 1,
 267             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.1, 0.1, 0, 0);
 268   AliMedium(701, "CPV scint.   $", 1, 1,
 269             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.1, 0.1, 0, 0);
 270   AliMedium(702, "Al parts     $", 2, 0,
 271             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.001, 0.001, 0, 0);
 272   AliMedium(703, "Tyvek wrapper$", 3, 0,
 273             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.001, 0.001, 0, 0);
 274   AliMedium(704, "Polyst. foam $", 4, 0,
 275             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.1, 0.1, 0, 0);
 276   AliMedium(705, "Titan. cover $", 5, 0,
 277             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.0001, 0.0001, 0, 0);
 278   AliMedium(706, "Si PIN       $", 6, 0,
 279             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.01, 0.01, 0, 0);
 280   AliMedium(707, "Thermo Insul.$", 7, 0,
 281             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.1, 0.1, 0, 0);
 282   AliMedium(708, "Textolit     $", 8, 0,
 283             ISXFLD, SXMGMX, 10.0, 0.1, 0.1, 0.1, 0.1, 0, 0);
 284   AliMedium(799, "Air          $", 99, 0,
 285             ISXFLD, SXMGMX, 10.0, 1.0, 0.1, 0.1, 10.0, 0, 0);
 286
 287   // --- Set decent energy thresholds for gamma and electron tracking
 288   pMC->Gstpar(idtmed[700],"CUTGAM",0.5E-4);
 289   pMC->Gstpar(idtmed[700],"CUTELE",1.0E-4);
 290   // --- Generate explicitly delta rays in the titan cover ---
 291   pMC->Gstpar(idtmed[705],"LOSS",3.);
 292   pMC->Gstpar(idtmed[705],"DRAY",1.);
 293   // --- and in aluminium parts ---
 294   pMC->Gstpar(idtmed[702],"LOSS",3.);
 295   pMC->Gstpar(idtmed[702],"DRAY",1.);
 296
 297 }
 298
 299 //////////////////////////////////////////////////////////////////////////////
 300
 301 void AliPHOSv2::CreateGeometry()
 302 {
 303
 304   AliMC* pMC = AliMC::GetMC();
 305
 306   AliPHOSv2 *PHOS_tmp = (AliPHOSv2*)gAlice->GetModule("PHOS");
 307   if(PHOS_tmp==NULL){
 308
 309     fprintf(stderr,"PHOS detector not found!\n");
 310     return;
 311
 312   }
 313
 314   // --- Dimensions of volumes ---
 315   Float_t DPHOS[3], DPTXW[3], DPAIR[3];
 316   Float_t DPUFP[3], DPUCP[3], DPASP[3], DPTIP[3], DPTXP[3];
 317   Float_t DPTCB[3], DPCBL[3], DPSTC[3], DPPAP[3], DPXTL[3], DPSUP[3],
 318     DPPIN[3];
 319   Float_t DPCPV[3], DPCPA[3];
 320
 321   Float_t R, YO, XP1, YP1, PPHI, angle;
 322   Int_t IDROTM[99];
 323   Int_t i;
 324
 325   Double_t const RADDEG=180.0/kPI;
 326   //  Double_t const DEGRAD=kPI/180.0;
 327
 328   // --- Dimensions of PbWO4 crystal ---
 329   //      PARAMETER(XTL_X=2.2,XTL_Y=18.,XTL_Z=2.2)
 330   Float_t XTL_X=GetCrystalSize(0);
 331   Float_t XTL_Y=GetCrystalSize(1);
 332   Float_t XTL_Z=GetCrystalSize(2);
 333
 334   // --- Tyvek wrapper thickness
 335   //      PARAMETER(PAP_THICK=0.01)
 336   Float_t PAP_THICK=GetWrapThickness();
 337
 338   // --- Steel (titanium) cover thickness ---
 339   //      PARAMETER(STE_THICK=0.005)
 340   Float_t STE_THICK=GetPHOSextra(0);
 341
 342   // --- Crystal support height ---
 343   //      PARAMETER(SUP_Y=6.95)
 344   Float_t SUP_Y=GetPHOSextra(1);
 345
 346   // --- PIN-diode dimensions ---
 347   //      PARAMETER(PIN_X=1.4,PIN_Y=0.4,PIN_Z=1.4)
 348   Float_t PIN_X=GetPINSize(0);
 349   Float_t PIN_Y=GetPINSize(1);
 350   Float_t PIN_Z=GetPINSize(2);
 351
 352   // --- CPV thickness ---
 353   //      PARAMETER(CPV_Y=0.5)
 354   Float_t CPV_Y=GetCPVThickness();
 355
 356   // --- Foam Thermo Insulating outer cover dimensions ---
 357   //      PARAMETER(FTI_X=214.6,FTI_Y=80.,FTI_Z=260.)
 358   Float_t FTI_X=GetPHOSFoam(0);
 359   Float_t FTI_Y=GetPHOSFoam(1);
 360   Float_t FTI_Z=GetPHOSFoam(2);
 361
 362   // --- Thermo Insulating outer cover Upper plate thickness ---
 363   //      PARAMETER(FTIU_THICK=4.)
 364   Float_t FTIU_THICK=GetPHOSextra(2);
 365
 366   // --- Textolit Wall box dimentions ---
 367   //      PARAMETER(TXW_X=209.,TXW_Y=71.,TXW_Z=250.)
 368   Float_t TXW_X=GetPHOStxwall(0);
 369   Float_t TXW_Y=GetPHOStxwall(1);
 370   Float_t TXW_Z=GetPHOStxwall(2);
 371
 372   // --- Inner AIR volume dimensions ---
 373   //      PARAMETER(AIR_X=206.,AIR_Y=66.,AIR_Z=244.)
 374   Float_t AIR_X=GetPHOSAir(0);
 375   Float_t AIR_Y=GetPHOSAir(1);
 376   Float_t AIR_Z=GetPHOSAir(2);
 377
 378   // --- Upper Polystyrene Foam plate thickness ---
 379   //      PARAMETER(UFP_Y=5.)
 380   Float_t UFP_Y=GetPHOSextra(3);
 381
 382   // --- Thermo insulating Crystal Block wall thickness ---
 383   //      PARAMETER(TCB_THICK=2.)
 384   Float_t TCB_THICK=GetPHOSextra(4);
 385
 386   // --- Upper Cooling Plate thickness ---
 387   //      PARAMETER(UCP_Y=0.06)
 388   Float_t UCP_Y=GetPHOSextra(5);
 389
 390   // --- Al Support Plate thickness ---
 391   //      PARAMETER(ASP_Y=10.)
 392   Float_t ASP_Y=GetPHOSextra(6);
 393
 394   // --- Lower Thermo Insulating Plate thickness ---
 395   //      PARAMETER(TIP_Y=3.)
 396   Float_t TIP_Y=GetPHOSextra(7);
 397
 398   // --- Lower Textolit Plate thickness ---
 399   //      PARAMETER(TXP_Y=1.)
 400   Float_t TXP_Y=GetPHOSextra(8);
 401
 402   // --- 1/2 total gap between adjacent crystals
 403   Float_t TOTAL_GAP=GetPHOSextra(9);
 404
 405   // --- Distance from IP to Foam Thermo Insulating top plate ---
 406   //      PARAMETER(FTI_R=467.)
 407   Float_t FTI_R=GetRadius(0);
 408
 409   // --- Distance from IP to Crystal Block top Surface (needs to be 460.) ---
 410   //      PARAMETER(CBS_R=480.)
 411   Float_t CBS_R=GetRadius(1);
 412
 413   // Get pointer to the array containing media indeces
 414   Int_t *IDTMED = gAlice->Idtmed();
 415
 416   // --- Define PHOS box volume, fill with thermo insulating foam ---
 417   DPHOS[0]=FTI_X/2.0;
 418   DPHOS[1]=FTI_Y/2.0;
 419   DPHOS[2]=FTI_Z/2.0;
 420   pMC->Gsvolu("PHOS", "BOX ", IDTMED[706], DPHOS, 3);
 421
 422   // --- Define Textolit Wall box, position inside PHOS ---
 423   DPTXW[0]=TXW_X/2.0;
 424   DPTXW[1]=TXW_Y/2.0;
 425   DPTXW[2]=TXW_Z/2.0;
 426   pMC->Gsvolu("PTXW", "BOX ", IDTMED[707], DPTXW, 3);
 427   YO=(FTI_Y-TXW_Y)/2.0-FTIU_THICK;
 428   pMC->Gspos("PTXW", 1, "PHOS", 0.0, YO, 0.0, 0, "ONLY");
 429
 430   // --- Define Upper Polystyrene Foam Plate, place inside PTXW ---
 431   // --- immediately below Foam Thermo Insulation Upper plate ---
 432   DPUFP[0]=TXW_X/2.0;
 433   DPUFP[1]=UFP_Y/2.0;
 434   DPUFP[2]=TXW_Z/2.0;
 435   pMC->Gsvolu("PUFP", "BOX ", IDTMED[703], DPUFP, 3);
 436   YO=(TXW_Y-UFP_Y)/2.0;
 437   pMC->Gspos("PUFP", 1, "PTXW", 0.0, YO, 0.0, 0, "ONLY");
 438
 439   // --- Define air-filled box, place inside PTXW ---
 440   DPAIR[0]=AIR_X/2.0;
 441   DPAIR[1]=AIR_Y/2.0;
 442   DPAIR[2]=AIR_Z/2.0;
 443   pMC->Gsvolu("PAIR", "BOX ", IDTMED[798], DPAIR, 3);
 444   YO=(TXW_Y-AIR_Y)/2.0-UFP_Y;
 445   pMC->Gspos("PAIR", 1, "PTXW", 0.0, YO, 0.0, 0, "ONLY");
 446
 447   // --- Define Thermo insulating Crystal Box, position inside PAIR ---
 448   DPTCB[0]=GetNphi()*(XTL_X+2*TOTAL_GAP)/2.0+TCB_THICK;
 449   DPTCB[1]=(XTL_Y+SUP_Y+PAP_THICK+STE_THICK)/2.0+TCB_THICK/2.0;
 450   DPTCB[2]=GetNz()*(XTL_Z+2*TOTAL_GAP)/2.0+TCB_THICK;
 451   pMC->Gsvolu("PTCB", "BOX ", IDTMED[706], DPTCB, 3);
 452   YO=AIR_Y/2.0-DPTCB[1]-
 453     (CBS_R-FTI_R-TCB_THICK-FTIU_THICK-UFP_Y);
 454   pMC->Gspos("PTCB", 1, "PAIR", 0.0, YO, 0.0, 0, "ONLY");
 455
 456   // --- Define Crystal BLock filled with air, position it inside PTCB ---
 457   DPCBL[0]=GetNphi()*(XTL_X+2*TOTAL_GAP)/2.0;
 458   DPCBL[1]=(XTL_Y+SUP_Y+PAP_THICK+STE_THICK)/2.0;
 459   DPCBL[2]=GetNz()*(XTL_Z+2*TOTAL_GAP)/2.0;
 460   pMC->Gsvolu("PCBL", "BOX ", IDTMED[798], DPCBL, 3);
 461
 462   // --- Divide PCBL in X (phi) and Z directions --
 463   pMC->Gsdvn("PROW", "PCBL", Int_t (GetNphi()), 1);
 464   pMC->Gsdvn("PCEL", "PROW", Int_t (GetNz()), 3);
 465   YO=-TCB_THICK/2.0;
 466   pMC->Gspos("PCBL", 1, "PTCB", 0.0, YO, 0.0, 0, "ONLY");
 467
 468   // --- Define STeel (actually, it's titanium) Cover volume, place inside PCEL
 469   DPSTC[0]=(XTL_X+2*PAP_THICK)/2.0;
 470   DPSTC[1]=(XTL_Y+SUP_Y+PAP_THICK+STE_THICK)/2.0;
 471   DPSTC[2]=(XTL_Z+2*PAP_THICK+2*STE_THICK)/2.0;
 472   pMC->Gsvolu("PSTC", "BOX ", IDTMED[704], DPSTC, 3);
 473   pMC->Gspos("PSTC", 1, "PCEL", 0.0, 0.0, 0.0, 0, "ONLY");
 474
 475   // --- Define Tyvek volume, place inside PSTC ---
 476   DPPAP[0]=XTL_X/2.0+PAP_THICK;
 477   DPPAP[1]=(XTL_Y+SUP_Y+PAP_THICK)/2.0;
 478   DPPAP[2]=XTL_Z/2.0+PAP_THICK;
 479   pMC->Gsvolu("PPAP", "BOX ", IDTMED[702], DPPAP, 3);
 480   YO=(XTL_Y+SUP_Y+PAP_THICK)/2.0-(XTL_Y+SUP_Y+PAP_THICK+STE_THICK)/2.0;
 481   pMC->Gspos("PPAP", 1, "PSTC", 0.0, YO, 0.0, 0, "ONLY");
 482
 483   // --- Define PbWO4 crystal volume, place inside PPAP ---
 484   DPXTL[0]=XTL_X/2.0;
 485   DPXTL[1]=XTL_Y/2.0;
 486   DPXTL[2]=XTL_Z/2.0;
 487   pMC->Gsvolu("PXTL", "BOX ", IDTMED[699], DPXTL, 3);
 488   YO=(XTL_Y+SUP_Y+PAP_THICK)/2.0-XTL_Y/2.0-PAP_THICK;
 489   pMC->Gspos("PXTL", 1, "PPAP", 0.0, YO, 0.0, 0, "ONLY");
 490
 491   // --- Define crystal support volume, place inside PPAP ---
 492   DPSUP[0]=XTL_X/2.0+PAP_THICK;
 493   DPSUP[1]=SUP_Y/2.0;
 494   DPSUP[2]=XTL_Z/2.0+PAP_THICK;
 495   pMC->Gsvolu("PSUP", "BOX ", IDTMED[798], DPSUP, 3);
 496   YO=SUP_Y/2.0-(XTL_Y+SUP_Y+PAP_THICK)/2.0;
 497   pMC->Gspos("PSUP", 1, "PPAP", 0.0, YO, 0.0, 0, "ONLY");
 498
 499   // --- Define PIN-diode volume and position it inside crystal support ---
 500   // --- right behind PbWO4 crystal
 501   DPPIN[0]=PIN_X/2.0;
 502   DPPIN[1]=PIN_Y/2.0;
 503   DPPIN[2]=PIN_Z/2.0;
 504   pMC->Gsvolu("PPIN", "BOX ", IDTMED[705], DPPIN, 3);
 505   YO=SUP_Y/2.0-PIN_Y/2.0;
 506   pMC->Gspos("PPIN", 1, "PSUP", 0.0, YO, 0.0, 0, "ONLY");
 507
 508   // --- Define Upper Cooling Panel, place it on top of PTCB ---
 509   DPUCP[0]=DPTCB[0];
 510   DPUCP[1]=UCP_Y/2.0;
 511   DPUCP[2]=DPTCB[2];
 512   pMC->Gsvolu("PUCP", "BOX ", IDTMED[701], DPUCP,3);
 513   YO=(AIR_Y-UCP_Y)/2.0-(CBS_R-FTI_R-TCB_THICK-FTIU_THICK-UFP_Y-UCP_Y);
 514   pMC->Gspos("PUCP", 1, "PAIR", 0.0, YO, 0.0, 0, "ONLY");
 515
 516   // --- Define Al Support Plate, position it inside PAIR ---
 517   // --- right beneath PTCB ---
 518   DPASP[0]=AIR_X/2.0;
 519   DPASP[1]=ASP_Y/2.0;
 520   DPASP[2]=AIR_Z/2.0;
 521   pMC->Gsvolu("PASP", "BOX ", IDTMED[701], DPASP, 3);
 522   YO=(AIR_Y-ASP_Y)/2.0-(CBS_R-FTI_R-FTIU_THICK-UFP_Y+DPCBL[1]*2);
 523   pMC->Gspos("PASP", 1, "PAIR", 0.0, YO, 0.0, 0, "ONLY");
 524
 525   // --- Define Thermo Insulating Plate, position it inside PAIR ---
 526   // --- right beneath PASP ---
 527   DPTIP[0]=AIR_X/2.0;
 528   DPTIP[1]=TIP_Y/2.0;
 529   DPTIP[2]=AIR_Z/2.0;
 530   pMC->Gsvolu("PTIP", "BOX ", IDTMED[706], DPTIP, 3);
 531   YO=(AIR_Y-TIP_Y)/2.0-(CBS_R-FTI_R-FTIU_THICK-UFP_Y+DPCBL[1]*2+ASP_Y);
 532   pMC->Gspos("PTIP", 1, "PAIR", 0.0, YO, 0.0, 0, "ONLY");
 533
 534   // --- Define Textolit Plate, position it inside PAIR ---
 535   // --- right beneath PTIP ---
 536   DPTXP[0]=AIR_X/2.0;
 537   DPTXP[1]=TXP_Y/2.0;
 538   DPTXP[2]=AIR_Z/2.0;
 539   pMC->Gsvolu("PTXP", "BOX ", IDTMED[707], DPTXP, 3);
 540   YO=(AIR_Y-TXP_Y)/2.0-
 541     (CBS_R-FTI_R-FTIU_THICK-UFP_Y+DPCBL[1]*2+ASP_Y+TIP_Y);
 542   pMC->Gspos("PTXP", 1, "PAIR", 0.0, YO, 0.0, 0, "ONLY");
 543
 544   // --- Define CPV volume, DON'T PLACE IT YET ---
 545   // --- Divide in X and Z direction (same way as PCBL) ---
 546   DPCPV[0]=DPCBL[0];
 547   DPCPV[1]=CPV_Y/2.0;
 548   DPCPV[2]=DPCBL[2];
 549   //  pMC->Gsvolu("PCPV", "BOX ", IDTMED[700], DPCPV, 3);
 550   pMC->Gsvolu("PCPV", "BOX ", IDTMED[798], DPCPV, 3);
 551   pMC->Gsdvn("PCRO", "PCPV", Int_t (GetNphi()), 1);
 552   pMC->Gsdvn("PCCE", "PCRO", Int_t (GetNz()), 3);
 553
 554   // Define CPV sensitive pad. It has the same size as PCCE.
 555   DPCPA[0]=DPCBL[0]/GetNphi();
 556   DPCPA[1]=CPV_Y/2.0;
 557   DPCPA[2]=DPCBL[2]/GetNz();
 558   pMC->Gsvolu("PCPA", "BOX ", IDTMED[700], DPCPA, 3);
 559   pMC->Gspos("PCPA", 1, "PCCE", 0.0, 0.0, 0.0, 0, "ONLY");
 560
 561   // --- Position various PHOS units in ALICE setup ---
 562   // --- PHOS itself first ---
 563   PPHI=TMath::ATan(FTI_X/(2.0*FTI_R));
 564   PPHI*=RADDEG;
 565
 566   for(i=1; i<=GetNModules(); i++){
 567
 568     angle=PPHI*2*(i-GetNModules()/2.0-0.5);
 569     AliMatrix(IDROTM[i-1], 90.0, angle, 90.0, 90.0+angle, 0.0, 0.0);
 570
 571     // --- Position various PHOS units in ALICE setup ---
 572     // --- PHOS itself first ---
 573     R=FTI_R+FTI_Y/2.0;
 574     XP1=R*TMath::Sin(angle/RADDEG);
 575     YP1=-R*TMath::Cos(angle/RADDEG);
 576     pMC->Gspos("PHOS", i, "ALIC", XP1, YP1, 0.0, IDROTM[i-1], "ONLY");
 577
 578     // --- Now position PCPV so that its plates are right on top of ---
 579     // --- corresponding PHOS modules (previously called cradles) ---
 580     R=FTI_R-CPV_Y/2.0;
 581     XP1=R*TMath::Sin(angle/RADDEG);
 582     YP1=-R*TMath::Cos(angle/RADDEG);
 583     pMC->Gspos("PCPV", i, "ALIC", XP1, YP1, 0.0, IDROTM[i-1], "ONLY");
 584     GetModuleAngle(i-1)=angle-90.0;
 585
 586   }
 587
 588   // --- Set volumes seen without their descendants for drawing ---
 589   pMC->Gsatt("PCEL", "SEEN", -2);
 590   pMC->Gsatt("PCCE", "SEEN", -2);
 591
 592 }
 593
 594 //////////////////////////////////////////////////////////////////////////////
 595
 596 void AliPHOSv2::StepManager(void)
 597 {
 598
 599   AliMC *pMC=AliMC::GetMC();
 600   Int_t blrc[4]; // (box, layer, row, column) indices
 601   Float_t xyze[4]; // position wrt MRS and energy deposited
 602
 603   Int_t *IDTMED=gAlice->Idtmed();
 604
 605   if(pMC->GetMedium()==IDTMED[700]){ // We are inside a CPV sensitive pad
 606
 607     pMC->TrackPosition(xyze);
 608     xyze[3]=pMC->Edep();
 609
 610     pMC->CurrentVolOff(3, (Text_t*)NULL, blrc[0]);
 611     blrc[1]=1; // CPV corresponds to layer 1
 612     pMC->CurrentVolOff(2, (Text_t*)NULL, blrc[2]);
 613     pMC->CurrentVolOff(1, (Text_t*)NULL, blrc[3]);
 614
 615     AddHit(gAlice->CurrentTrack(), blrc, xyze);
 616
 617   }
 618
 619   if(pMC->GetMedium()==IDTMED[699]){ // We are inside a PWO crystal
 620
 621     pMC->TrackPosition(xyze);
 622     xyze[3]=pMC->Edep();
 623
 624     pMC->CurrentVolOff(9, (Text_t*)NULL, blrc[0]);
 625     blrc[1]=2; // PWO crystals correspond to layer 2
 626     pMC->CurrentVolOff(4, (Text_t*)NULL, blrc[2]);
 627     pMC->CurrentVolOff(3, (Text_t*)NULL, blrc[3]);
 628
 629     AddHit(gAlice->CurrentTrack(), blrc, xyze);
 630
 631   }
 632 }
 633
 634 ////////////////////////////////////////////////////////////////////////////
 635
 636 void AliPHOSv2::AddHit(Int_t track, Int_t *vol, Float_t *hits)
 637 {
 638
 639   Int_t hitCounter;
 640   TClonesArray &lhits = *fHits;
 641   AliPHOShitv2 *newHit,*curHit;
 642
 643   newHit=new AliPHOShitv2(fIshunt, track, vol, hits);
 644
 645   for(hitCounter=0;hitCounter<fNhits;hitCounter++){
 646     curHit=(AliPHOShitv2*)lhits[hitCounter];
 647     if(*curHit==*newHit){
 648       *curHit=*curHit+*newHit;
 649       delete newHit;
 650       return;
 651     }
 652   }
 653
 654   new(lhits[fNhits++]) AliPHOShitv2(*newHit);
 655   delete newHit;
 656
 657 }
 658
 659 ///////////////////////////////////////////////////////////////////////////////
 660
 661 ClassImp(AliPHOShitv2)
 662
 663 //////////////////////////////////////////////////////////////////////////////
 664
 665 AliPHOShitv2::AliPHOShitv2(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
 666 AliHit(shunt, track)
 667 {
 668
 669    Int_t i;
 670    for (i=0;i<4;i++)fVolume[i]=vol[i];
 671
 672    fX       = hits[0];
 673    fY       = hits[1];
 674    fZ       = hits[2];
 675    fELOS    = hits[3];
 676
 677 }
 678
 679 //////////////////////////////////////////////////////////////////////////////
 680
 681 Bool_t AliPHOShitv2::operator==(AliPHOShitv2 const &rValue) const
 682 {
 683
 684   Int_t volCounter;
 685
 686   //  if(fDet!=rValue.GetDet()) return kFALSE;
 687   if(fTrack!=rValue.GetTrack()) return kFALSE;
 688
 689   for(volCounter=0;volCounter<4;volCounter++)
 690     if(fVolume[volCounter]!=rValue.GetVolume(volCounter))return kFALSE;
 691
 692   return kTRUE;
 693
 694 }
 695
 696 /////////////////////////////////////////////////////////////////////////////
 697
 698 AliPHOShitv2 const AliPHOShitv2::operator+(AliPHOShitv2 const &rValue) const
 699 {
 700
 701   AliPHOShitv2 added(*this);
 702
 703   added.fELOS+=rValue.GetEnergy();
 704   return added;
 705
 706 }
 707
 708 //////////////////////////////////////////////////////////////////////////////