1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //_________________________________________________________________________
19 // Implementation version v0 of PHOS Manager class
20 // Layout EMC + PPSD has name GPS2
21 // An object of this class does not produce hits nor digits
22 // It is the one to use if you do not want to produce outputs in TREEH or TREED
24 //*-- Author: Yves Schutz (SUBATECH)
27 // --- ROOT system ---
34 // --- Standard library ---
39 #include <strstream.h>
41 // --- AliRoot header files ---
43 #include "AliPHOSv0.h"
49 //____________________________________________________________________________
50 AliPHOSv0::AliPHOSv0(const char *name, const char *title):
53 // ctor : title is used to identify the layout
54 // GPS2 = 5 modules (EMC + PPSD)
56 // gets an instance of the geometry parameters class
58 if (strcmp(title,"") != 0 )
59 fGeom = AliPHOSGeometry::GetInstance(title, "") ;
61 cout << "AliPHOS" << Version() << " : PHOS geometry intialized for " << fGeom->GetName() << endl ;
63 cout << "AliPHOS" << Version() << " : PHOS geometry initialization failed !" << endl ;
66 //____________________________________________________________________________
67 void AliPHOSv0::BuildGeometry()
69 // Build the PHOS geometry for the ROOT display
73 PHOS in ALICE displayed by root
79 <IMG Align=BOTTOM ALT="All Views" SRC="../images/AliPHOSv0AllViews.gif">
84 <IMG Align=BOTTOM ALT="Front View" SRC="../images/AliPHOSv0FrontView.gif">
89 <IMG Align=BOTTOM ALT="3D View 1" SRC="../images/AliPHOSv03DView1.gif">
94 <IMG Align=BOTTOM ALT="3D View 2" SRC="../images/AliPHOSv03DView2.gif">
100 this->BuildGeometryforPHOS() ;
101 if ( ( strcmp(fGeom->GetName(), "GPS2" ) == 0 ) )
102 this->BuildGeometryforPPSD() ;
104 cout << "AliPHOSv0::BuildGeometry : no charged particle identification system installed" << endl;
108 //____________________________________________________________________________
109 void AliPHOSv0:: BuildGeometryforPHOS(void)
111 // Build the PHOS-EMC geometry for the ROOT display
113 const Int_t kColorPHOS = kRed ;
114 const Int_t kColorXTAL = kBlue ;
116 Double_t const kRADDEG = 180.0 / kPI ;
118 new TBRIK( "OuterBox", "PHOS box", "void", fGeom->GetOuterBoxSize(0)/2,
119 fGeom->GetOuterBoxSize(1)/2,
120 fGeom->GetOuterBoxSize(2)/2 );
122 // Textolit Wall box, position inside PHOS
124 new TBRIK( "TextolitBox", "PHOS Textolit box ", "void", fGeom->GetTextolitBoxSize(0)/2,
125 fGeom->GetTextolitBoxSize(1)/2,
126 fGeom->GetTextolitBoxSize(2)/2);
128 // Polystyrene Foam Plate
130 new TBRIK( "UpperFoamPlate", "PHOS Upper foam plate", "void", fGeom->GetTextolitBoxSize(0)/2,
131 fGeom->GetSecondUpperPlateThickness()/2,
132 fGeom->GetTextolitBoxSize(2)/2 ) ;
136 new TBRIK( "AirFilledBox", "PHOS air filled box", "void", fGeom->GetAirFilledBoxSize(0)/2,
137 fGeom->GetAirFilledBoxSize(1)/2,
138 fGeom->GetAirFilledBoxSize(2)/2 );
142 Float_t xtlX = fGeom->GetCrystalSize(0) ;
143 Float_t xtlY = fGeom->GetCrystalSize(1) ;
144 Float_t xtlZ = fGeom->GetCrystalSize(2) ;
146 Float_t xl = fGeom->GetNPhi() * ( xtlX + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 + fGeom->GetModuleBoxThickness() ;
147 Float_t yl = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0
148 + fGeom->GetModuleBoxThickness() / 2.0 ;
149 Float_t zl = fGeom->GetNZ() * ( xtlZ + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 + fGeom->GetModuleBoxThickness() ;
151 new TBRIK( "CrystalsBox", "PHOS crystals box", "void", xl, yl, zl ) ;
153 // position PHOS into ALICE
155 Float_t r = fGeom->GetIPtoOuterCoverDistance() + fGeom->GetOuterBoxSize(1) / 2.0 ;
157 Float_t pphi = TMath::ATan( fGeom->GetOuterBoxSize(0) / ( 2.0 * fGeom->GetIPtoOuterCoverDistance() ) ) ;
159 TNode * top = gAlice->GetGeometry()->GetNode("alice") ;
161 char * nodename = new char[20] ;
162 char * rotname = new char[20] ;
164 for( Int_t i = 1; i <= fGeom->GetNModules(); i++ ) {
165 Float_t angle = pphi * 2 * ( i - fGeom->GetNModules() / 2.0 - 0.5 ) ;
166 sprintf(rotname, "%s%d", "rot", number++) ;
167 new TRotMatrix(rotname, rotname, 90, angle, 90, 90 + angle, 0, 0);
169 sprintf(nodename,"%s%d", "Module", i) ;
170 Float_t x = r * TMath::Sin( angle / kRADDEG ) ;
171 Float_t y = -r * TMath::Cos( angle / kRADDEG ) ;
172 TNode * outerboxnode = new TNode(nodename, nodename, "OuterBox", x, y, 0, rotname ) ;
173 outerboxnode->SetLineColor(kColorPHOS) ;
174 fNodes->Add(outerboxnode) ;
176 // now inside the outer box the textolit box
177 y = ( fGeom->GetOuterBoxThickness(1) - fGeom->GetUpperPlateThickness() ) / 2. ;
178 sprintf(nodename,"%s%d", "TexBox", i) ;
179 TNode * textolitboxnode = new TNode(nodename, nodename, "TextolitBox", 0, y, 0) ;
180 textolitboxnode->SetLineColor(kColorPHOS) ;
181 fNodes->Add(textolitboxnode) ;
182 // upper foam plate inside outre box
184 sprintf(nodename, "%s%d", "UFPlate", i) ;
185 y = ( fGeom->GetTextolitBoxSize(1) - fGeom->GetSecondUpperPlateThickness() ) / 2.0 ;
186 TNode * upperfoamplatenode = new TNode(nodename, nodename, "UpperFoamPlate", 0, y, 0) ;
187 upperfoamplatenode->SetLineColor(kColorPHOS) ;
188 fNodes->Add(upperfoamplatenode) ;
189 // air filled box inside textolit box (not drawn)
190 textolitboxnode->cd();
191 y = ( fGeom->GetTextolitBoxSize(1) - fGeom->GetAirFilledBoxSize(1) ) / 2.0 - fGeom->GetSecondUpperPlateThickness() ;
192 sprintf(nodename, "%s%d", "AFBox", i) ;
193 TNode * airfilledboxnode = new TNode(nodename, nodename, "AirFilledBox", 0, y, 0) ;
194 fNodes->Add(airfilledboxnode) ;
195 // crystals box inside air filled box
196 airfilledboxnode->cd() ;
197 y = fGeom->GetAirFilledBoxSize(1) / 2.0 - yl
198 - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetModuleBoxThickness()
199 - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() ) ;
200 sprintf(nodename, "%s%d", "XTBox", i) ;
201 TNode * crystalsboxnode = new TNode(nodename, nodename, "CrystalsBox", 0, y, 0) ;
202 crystalsboxnode->SetLineColor(kColorXTAL) ;
203 fNodes->Add(crystalsboxnode) ;
210 //____________________________________________________________________________
211 void AliPHOSv0:: BuildGeometryforPPSD(void)
213 // Build the PHOS-PPSD geometry for the ROOT display
217 PPSD displayed by root
220 <LI> Zoom on PPSD: Front View
223 <IMG Align=BOTTOM ALT="PPSD Front View" SRC="../images/AliPHOSv0PPSDFrontView.gif">
225 <LI> Zoom on PPSD: Perspective View
228 <IMG Align=BOTTOM ALT="PPSD Prespective View" SRC="../images/AliPHOSv0PPSDPerspectiveView.gif">
233 Double_t const kRADDEG = 180.0 / kPI ;
235 const Int_t kColorPHOS = kRed ;
236 const Int_t kColorPPSD = kGreen ;
237 const Int_t kColorGas = kBlue ;
238 const Int_t kColorAir = kYellow ;
240 // Box for a full PHOS module
242 new TBRIK( "PPSDBox", "PPSD box", "void", fGeom->GetPPSDBoxSize(0)/2,
243 fGeom->GetPPSDBoxSize(1)/2,
244 fGeom->GetPPSDBoxSize(2)/2 );
246 // Box containing one micromegas module
248 new TBRIK( "PPSDModule", "PPSD module", "void", fGeom->GetPPSDModuleSize(0)/2,
249 fGeom->GetPPSDModuleSize(1)/2,
250 fGeom->GetPPSDModuleSize(2)/2 );
253 new TBRIK ( "TopLid", "Micromegas top lid", "void", fGeom->GetPPSDModuleSize(0)/2,
254 fGeom->GetLidThickness()/2,
255 fGeom->GetPPSDModuleSize(2)/2 ) ;
256 // composite panel (top and bottom)
258 new TBRIK ( "TopPanel", "Composite top panel", "void", ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() )/2,
259 fGeom->GetCompositeThickness()/2,
260 ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() )/2 ) ;
262 new TBRIK ( "BottomPanel", "Composite bottom panel", "void", ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() )/2,
263 fGeom->GetCompositeThickness()/2,
264 ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() )/2 ) ;
265 // gas gap (conversion and avalanche)
267 new TBRIK ( "GasGap", "gas gap", "void", ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() )/2,
268 ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() )/2,
269 ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() )/2 ) ;
273 new TBRIK ( "Anode", "Anode", "void", ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() )/2,
274 fGeom->GetAnodeThickness()/2,
275 ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() )/2 ) ;
277 new TBRIK ( "Cathode", "Cathode", "void", ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() )/2,
278 fGeom->GetCathodeThickness()/2,
279 ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() )/2 ) ;
282 new TBRIK ( "PCBoard", "Printed Circuit", "void", ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() )/2,
283 fGeom->GetPCThickness()/2,
284 ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() )/2 ) ;
285 // Gap between Lead and top micromegas
287 new TBRIK ( "LeadToM", "Air Gap top", "void", fGeom->GetPPSDBoxSize(0)/2,
288 fGeom->GetMicro1ToLeadGap()/2,
289 fGeom->GetPPSDBoxSize(2)/2 ) ;
291 // Gap between Lead and bottom micromegas
293 new TBRIK ( "MToLead", "Air Gap bottom", "void", fGeom->GetPPSDBoxSize(0)/2,
294 fGeom->GetLeadToMicro2Gap()/2,
295 fGeom->GetPPSDBoxSize(2)/2 ) ;
298 new TBRIK ( "Lead", "Lead converter", "void", fGeom->GetPPSDBoxSize(0)/2,
299 fGeom->GetLeadConverterThickness()/2,
300 fGeom->GetPPSDBoxSize(2)/2 ) ;
302 // position PPSD into ALICE
304 char * nodename = new char[20] ;
305 char * rotname = new char[20] ;
307 Float_t r = fGeom->GetIPtoTopLidDistance() + fGeom->GetPPSDBoxSize(1) / 2.0 ;
309 TNode * top = gAlice->GetGeometry()->GetNode("alice") ;
311 for( Int_t i = 1; i <= fGeom->GetNModules(); i++ ) { // the number of PHOS modules
312 Float_t angle = fGeom->GetPHOSAngle(i) ;
313 sprintf(rotname, "%s%d", "rotg", number++) ;
314 new TRotMatrix(rotname, rotname, 90, angle, 90, 90 + angle, 0, 0);
316 sprintf(nodename, "%s%d", "Moduleg", i) ;
317 Float_t x = r * TMath::Sin( angle / kRADDEG ) ;
318 Float_t y = -r * TMath::Cos( angle / kRADDEG ) ;
319 TNode * ppsdboxnode = new TNode(nodename , nodename ,"PPSDBox", x, y, 0, rotname ) ;
320 ppsdboxnode->SetLineColor(kColorPPSD) ;
321 fNodes->Add(ppsdboxnode) ;
323 // inside the PPSD box:
324 // 1. fNumberOfModulesPhi x fNumberOfModulesZ top micromegas
325 x = ( fGeom->GetPPSDBoxSize(0) - fGeom->GetPPSDModuleSize(0) ) / 2. ;
327 for ( Int_t iphi = 1; iphi <= fGeom->GetNumberOfModulesPhi(); iphi++ ) { // the number of micromegas modules in phi per PHOS module
328 Float_t z = ( fGeom->GetPPSDBoxSize(2) - fGeom->GetPPSDModuleSize(2) ) / 2. ;
330 for ( Int_t iz = 1; iz <= fGeom->GetNumberOfModulesZ(); iz++ ) { // the number of micromegas modules in z per PHOS module
331 y = ( fGeom->GetPPSDBoxSize(1) - fGeom->GetMicromegas1Thickness() ) / 2. ;
332 sprintf(nodename, "%s%d%d%d", "Mic1", i, iphi, iz) ;
333 micro1node = new TNode(nodename, nodename, "PPSDModule", x, y, z) ;
334 micro1node->SetLineColor(kColorPPSD) ;
335 fNodes->Add(micro1node) ;
336 // inside top micromegas
339 y = ( fGeom->GetMicromegas1Thickness() - fGeom->GetLidThickness() ) / 2. ;
340 sprintf(nodename, "%s%d%d%d", "Lid", i, iphi, iz) ;
341 TNode * toplidnode = new TNode(nodename, nodename, "TopLid", 0, y, 0) ;
342 toplidnode->SetLineColor(kColorPPSD) ;
343 fNodes->Add(toplidnode) ;
344 // b. composite panel
345 y = y - fGeom->GetLidThickness() / 2. - fGeom->GetCompositeThickness() / 2. ;
346 sprintf(nodename, "%s%d%d%d", "CompU", i, iphi, iz) ;
347 TNode * compupnode = new TNode(nodename, nodename, "TopPanel", 0, y, 0) ;
348 compupnode->SetLineColor(kColorPPSD) ;
349 fNodes->Add(compupnode) ;
351 y = y - fGeom->GetCompositeThickness() / 2. - fGeom->GetAnodeThickness() / 2. ;
352 sprintf(nodename, "%s%d%d%d", "Ano", i, iphi, iz) ;
353 TNode * anodenode = new TNode(nodename, nodename, "Anode", 0, y, 0) ;
354 anodenode->SetLineColor(kColorPHOS) ;
355 fNodes->Add(anodenode) ;
357 y = y - fGeom->GetAnodeThickness() / 2. - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. ;
358 sprintf(nodename, "%s%d%d%d", "GGap", i, iphi, iz) ;
359 TNode * ggapnode = new TNode(nodename, nodename, "GasGap", 0, y, 0) ;
360 ggapnode->SetLineColor(kColorGas) ;
361 fNodes->Add(ggapnode) ;
363 y = y - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. - fGeom->GetCathodeThickness() / 2. ;
364 sprintf(nodename, "%s%d%d%d", "Cathode", i, iphi, iz) ;
365 TNode * cathodenode = new TNode(nodename, nodename, "Cathode", 0, y, 0) ;
366 cathodenode->SetLineColor(kColorPHOS) ;
367 fNodes->Add(cathodenode) ;
368 // g. printed circuit
369 y = y - fGeom->GetCathodeThickness() / 2. - fGeom->GetPCThickness() / 2. ;
370 sprintf(nodename, "%s%d%d%d", "PC", i, iphi, iz) ;
371 TNode * pcnode = new TNode(nodename, nodename, "PCBoard", 0, y, 0) ;
372 pcnode->SetLineColor(kColorPPSD) ;
373 fNodes->Add(pcnode) ;
374 // h. composite panel
375 y = y - fGeom->GetPCThickness() / 2. - fGeom->GetCompositeThickness() / 2. ;
376 sprintf(nodename, "%s%d%d%d", "CompDown", i, iphi, iz) ;
377 TNode * compdownnode = new TNode(nodename, nodename, "BottomPanel", 0, y, 0) ;
378 compdownnode->SetLineColor(kColorPPSD) ;
379 fNodes->Add(compdownnode) ;
380 z = z - fGeom->GetPPSDModuleSize(2) ;
382 } // end of Z module loop
383 x = x - fGeom->GetPPSDModuleSize(0) ;
385 } // end of phi module loop
389 y = ( fGeom->GetPPSDBoxSize(1) - 2 * fGeom->GetMicromegas1Thickness() - fGeom->GetMicro1ToLeadGap() ) / 2. ;
390 sprintf(nodename, "%s%d", "GapUp", i) ;
391 TNode * gapupnode = new TNode(nodename, nodename, "LeadToM", 0, y, 0) ;
392 gapupnode->SetLineColor(kColorAir) ;
393 fNodes->Add(gapupnode) ;
395 y = y - fGeom->GetMicro1ToLeadGap() / 2. - fGeom->GetLeadConverterThickness() / 2. ;
396 sprintf(nodename, "%s%d", "LeadC", i) ;
397 TNode * leadcnode = new TNode(nodename, nodename, "Lead", 0, y, 0) ;
398 leadcnode->SetLineColor(kColorPPSD) ;
399 fNodes->Add(leadcnode) ;
401 y = y - fGeom->GetLeadConverterThickness() / 2. - fGeom->GetLeadToMicro2Gap() / 2. ;
402 sprintf(nodename, "%s%d", "GapDown", i) ;
403 TNode * gapdownnode = new TNode(nodename, nodename, "MToLead", 0, y, 0) ;
404 gapdownnode->SetLineColor(kColorAir) ;
405 fNodes->Add(gapdownnode) ;
406 // 5. fNumberOfModulesPhi x fNumberOfModulesZ bottom micromegas
407 x = ( fGeom->GetPPSDBoxSize(0) - fGeom->GetPPSDModuleSize(0) ) / 2. - fGeom->GetPhiDisplacement() ;
409 for ( Int_t iphi = 1; iphi <= fGeom->GetNumberOfModulesPhi(); iphi++ ) {
410 Float_t z = ( fGeom->GetPPSDBoxSize(2) - fGeom->GetPPSDModuleSize(2) ) / 2. - fGeom->GetZDisplacement() ;;
412 for ( Int_t iz = 1; iz <= fGeom->GetNumberOfModulesZ(); iz++ ) {
413 y = - ( fGeom->GetPPSDBoxSize(1) - fGeom->GetMicromegas2Thickness() ) / 2. ;
414 sprintf(nodename, "%s%d%d%d", "Mic2", i, iphi, iz) ;
415 micro2node = new TNode(nodename, nodename, "PPSDModule", x, y, z) ;
416 micro2node->SetLineColor(kColorPPSD) ;
417 fNodes->Add(micro2node) ;
418 // inside bottom micromegas
421 y = ( fGeom->GetMicromegas2Thickness() - fGeom->GetLidThickness() ) / 2. ;
422 sprintf(nodename, "%s%d", "Lidb", i) ;
423 TNode * toplidbnode = new TNode(nodename, nodename, "TopLid", 0, y, 0) ;
424 toplidbnode->SetLineColor(kColorPPSD) ;
425 fNodes->Add(toplidbnode) ;
426 // b. composite panel
427 y = y - fGeom->GetLidThickness() / 2. - fGeom->GetCompositeThickness() / 2. ;
428 sprintf(nodename, "%s%d", "CompUb", i) ;
429 TNode * compupbnode = new TNode(nodename, nodename, "TopPanel", 0, y, 0) ;
430 compupbnode->SetLineColor(kColorPPSD) ;
431 fNodes->Add(compupbnode) ;
433 y = y - fGeom->GetCompositeThickness() / 2. - fGeom->GetAnodeThickness() / 2. ;
434 sprintf(nodename, "%s%d", "Anob", i) ;
435 TNode * anodebnode = new TNode(nodename, nodename, "Anode", 0, y, 0) ;
436 anodebnode->SetLineColor(kColorPPSD) ;
437 fNodes->Add(anodebnode) ;
439 y = y - fGeom->GetAnodeThickness() / 2. - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. ;
440 sprintf(nodename, "%s%d", "GGapb", i) ;
441 TNode * ggapbnode = new TNode(nodename, nodename, "GasGap", 0, y, 0) ;
442 ggapbnode->SetLineColor(kColorGas) ;
443 fNodes->Add(ggapbnode) ;
445 y = y - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. - fGeom->GetCathodeThickness() / 2. ;
446 sprintf(nodename, "%s%d", "Cathodeb", i) ;
447 TNode * cathodebnode = new TNode(nodename, nodename, "Cathode", 0, y, 0) ;
448 cathodebnode->SetLineColor(kColorPPSD) ;
449 fNodes->Add(cathodebnode) ;
450 // g. printed circuit
451 y = y - fGeom->GetCathodeThickness() / 2. - fGeom->GetPCThickness() / 2. ;
452 sprintf(nodename, "%s%d", "PCb", i) ;
453 TNode * pcbnode = new TNode(nodename, nodename, "PCBoard", 0, y, 0) ;
454 pcbnode->SetLineColor(kColorPPSD) ;
455 fNodes->Add(pcbnode) ;
457 y = y - fGeom->GetPCThickness() / 2. - fGeom->GetCompositeThickness() / 2. ;
458 sprintf(nodename, "%s%d", "CompDownb", i) ;
459 TNode * compdownbnode = new TNode(nodename, nodename, "BottomPanel", 0, y, 0) ;
460 compdownbnode->SetLineColor(kColorPPSD) ;
461 fNodes->Add(compdownbnode) ;
462 z = z - fGeom->GetPPSDModuleSize(2) ;
464 } // end of Z module loop
465 x = x - fGeom->GetPPSDModuleSize(0) ;
467 } // end of phi module loop
476 //____________________________________________________________________________
477 void AliPHOSv0::CreateGeometry()
479 // Create the PHOS geometry for Geant
481 AliPHOSv0 *phostmp = (AliPHOSv0*)gAlice->GetModule("PHOS") ;
483 if ( phostmp == NULL ) {
485 fprintf(stderr, "PHOS detector not found!\n") ;
489 // Get pointer to the array containing media indeces
490 Int_t *idtmed = fIdtmed->GetArray() - 699 ;
493 bigbox[0] = fGeom->GetOuterBoxSize(0) / 2.0 ;
494 bigbox[1] = ( fGeom->GetOuterBoxSize(1) + fGeom->GetPPSDBoxSize(1) ) / 2.0 ;
495 bigbox[2] = fGeom->GetOuterBoxSize(2) / 2.0 ;
497 gMC->Gsvolu("PHOS", "BOX ", idtmed[798], bigbox, 3) ;
499 this->CreateGeometryforPHOS() ;
500 if ( strcmp( fGeom->GetName(), "GPS2") == 0 )
501 this->CreateGeometryforPPSD() ;
503 cout << "AliPHOSv0::CreateGeometry : no charged particle identification system installed" << endl;
505 // --- Position PHOS mdules in ALICE setup ---
508 Double_t const kRADDEG = 180.0 / kPI ;
510 for( Int_t i = 1; i <= fGeom->GetNModules(); i++ ) {
512 Float_t angle = fGeom->GetPHOSAngle(i) ;
513 AliMatrix(idrotm[i-1], 90.0, angle, 90.0, 90.0+angle, 0.0, 0.0) ;
515 Float_t r = fGeom->GetIPtoOuterCoverDistance() + ( fGeom->GetOuterBoxSize(1) + fGeom->GetPPSDBoxSize(1) ) / 2.0 ;
517 Float_t xP1 = r * TMath::Sin( angle / kRADDEG ) ;
518 Float_t yP1 = -r * TMath::Cos( angle / kRADDEG ) ;
520 gMC->Gspos("PHOS", i, "ALIC", xP1, yP1, 0.0, idrotm[i-1], "ONLY") ;
526 //____________________________________________________________________________
527 void AliPHOSv0::CreateGeometryforPHOS()
529 // Create the PHOS-EMC geometry for GEANT
533 Geant3 geometry tree of PHOS-EMC in ALICE
536 <IMG Align=BOTTOM ALT="EMC geant tree" SRC="../images/EMCinAlice.gif">
541 // Get pointer to the array containing media indexes
542 Int_t *idtmed = fIdtmed->GetArray() - 699 ;
545 // --- Define PHOS box volume, fPUFPill with thermo insulating foam ---
546 // --- Foam Thermo Insulating outer cover dimensions ---
547 // --- Put it in bigbox = PHOS
550 dphos[0] = fGeom->GetOuterBoxSize(0) / 2.0 ;
551 dphos[1] = fGeom->GetOuterBoxSize(1) / 2.0 ;
552 dphos[2] = fGeom->GetOuterBoxSize(2) / 2.0 ;
554 gMC->Gsvolu("EMCA", "BOX ", idtmed[706], dphos, 3) ;
556 Float_t yO = - fGeom->GetPPSDBoxSize(1) / 2.0 ;
558 gMC->Gspos("EMCA", 1, "PHOS", 0.0, yO, 0.0, 0, "ONLY") ;
561 // --- Define Textolit Wall box, position inside EMCA ---
562 // --- Textolit Wall box dimentions ---
566 dptxw[0] = fGeom->GetTextolitBoxSize(0) / 2.0 ;
567 dptxw[1] = fGeom->GetTextolitBoxSize(1) / 2.0 ;
568 dptxw[2] = fGeom->GetTextolitBoxSize(2) / 2.0 ;
570 gMC->Gsvolu("PTXW", "BOX ", idtmed[707], dptxw, 3);
572 yO = ( fGeom->GetOuterBoxThickness(1) - fGeom->GetUpperPlateThickness() ) / 2. ;
574 gMC->Gspos("PTXW", 1, "EMCA", 0.0, yO, 0.0, 0, "ONLY") ;
577 // --- Define Upper Polystyrene Foam Plate, place inside PTXW ---
578 // --- immediately below Foam Thermo Insulation Upper plate ---
580 // --- Upper Polystyrene Foam plate thickness ---
583 dpufp[0] = fGeom->GetTextolitBoxSize(0) / 2.0 ;
584 dpufp[1] = fGeom->GetSecondUpperPlateThickness() / 2. ;
585 dpufp[2] = fGeom->GetTextolitBoxSize(2) /2.0 ;
587 gMC->Gsvolu("PUFP", "BOX ", idtmed[703], dpufp, 3) ;
589 yO = ( fGeom->GetTextolitBoxSize(1) - fGeom->GetSecondUpperPlateThickness() ) / 2.0 ;
591 gMC->Gspos("PUFP", 1, "PTXW", 0.0, yO, 0.0, 0, "ONLY") ;
594 // --- Define air-filled box, place inside PTXW ---
595 // --- Inner AIR volume dimensions ---
599 dpair[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ;
600 dpair[1] = fGeom->GetAirFilledBoxSize(1) / 2.0 ;
601 dpair[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ;
603 gMC->Gsvolu("PAIR", "BOX ", idtmed[798], dpair, 3) ;
605 yO = ( fGeom->GetTextolitBoxSize(1) - fGeom->GetAirFilledBoxSize(1) ) / 2.0 - fGeom->GetSecondUpperPlateThickness() ;
607 gMC->Gspos("PAIR", 1, "PTXW", 0.0, yO, 0.0, 0, "ONLY") ;
609 // --- Dimensions of PbWO4 crystal ---
611 Float_t xtlX = fGeom->GetCrystalSize(0) ;
612 Float_t xtlY = fGeom->GetCrystalSize(1) ;
613 Float_t xtlZ = fGeom->GetCrystalSize(2) ;
616 dptcb[0] = fGeom->GetNPhi() * ( xtlX + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 + fGeom->GetModuleBoxThickness() ;
617 dptcb[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0
618 + fGeom->GetModuleBoxThickness() / 2.0 ;
619 dptcb[2] = fGeom->GetNZ() * ( xtlZ + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 + fGeom->GetModuleBoxThickness() ;
621 gMC->Gsvolu("PTCB", "BOX ", idtmed[706], dptcb, 3) ;
623 yO = fGeom->GetAirFilledBoxSize(1) / 2.0 - dptcb[1]
624 - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetModuleBoxThickness()
625 - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() ) ;
627 gMC->Gspos("PTCB", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ;
630 // --- Define Crystal BLock filled with air, position it inside PTCB ---
633 dpcbl[0] = fGeom->GetNPhi() * ( xtlX + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 ;
634 dpcbl[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 ;
635 dpcbl[2] = fGeom->GetNZ() * ( xtlZ + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 ;
637 gMC->Gsvolu("PCBL", "BOX ", idtmed[798], dpcbl, 3) ;
639 // --- Divide PCBL in X (phi) and Z directions --
640 gMC->Gsdvn("PROW", "PCBL", Int_t (fGeom->GetNPhi()), 1) ;
641 gMC->Gsdvn("PCEL", "PROW", Int_t (fGeom->GetNZ()), 3) ;
643 yO = -fGeom->GetModuleBoxThickness() / 2.0 ;
645 gMC->Gspos("PCBL", 1, "PTCB", 0.0, yO, 0.0, 0, "ONLY") ;
648 // --- Define STeel (actually, it's titanium) Cover volume, place inside PCEL
651 dpstc[0] = ( xtlX + 2 * fGeom->GetCrystalWrapThickness() ) / 2.0 ;
652 dpstc[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 ;
653 dpstc[2] = ( xtlZ + 2 * fGeom->GetCrystalWrapThickness() + 2 * fGeom->GetCrystalHolderThickness() ) / 2.0 ;
655 gMC->Gsvolu("PSTC", "BOX ", idtmed[704], dpstc, 3) ;
657 gMC->Gspos("PSTC", 1, "PCEL", 0.0, 0.0, 0.0, 0, "ONLY") ;
660 // --- Define Tyvek volume, place inside PSTC ---
663 dppap[0] = xtlX / 2.0 + fGeom->GetCrystalWrapThickness() ;
664 dppap[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 ;
665 dppap[2] = xtlZ / 2.0 + fGeom->GetCrystalWrapThickness() ;
667 gMC->Gsvolu("PPAP", "BOX ", idtmed[702], dppap, 3) ;
669 yO = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0
670 - ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 ;
672 gMC->Gspos("PPAP", 1, "PSTC", 0.0, yO, 0.0, 0, "ONLY") ;
675 // --- Define PbWO4 crystal volume, place inside PPAP ---
678 dpxtl[0] = xtlX / 2.0 ;
679 dpxtl[1] = xtlY / 2.0 ;
680 dpxtl[2] = xtlZ / 2.0 ;
682 gMC->Gsvolu("PXTL", "BOX ", idtmed[699], dpxtl, 3) ;
684 yO = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 - xtlY / 2.0 - fGeom->GetCrystalWrapThickness() ;
686 gMC->Gspos("PXTL", 1, "PPAP", 0.0, yO, 0.0, 0, "ONLY") ;
689 // --- Define crystal support volume, place inside PPAP ---
692 dpsup[0] = xtlX / 2.0 + fGeom->GetCrystalWrapThickness() ;
693 dpsup[1] = fGeom->GetCrystalSupportHeight() / 2.0 ;
694 dpsup[2] = xtlZ / 2.0 + fGeom->GetCrystalWrapThickness() ;
696 gMC->Gsvolu("PSUP", "BOX ", idtmed[798], dpsup, 3) ;
698 yO = fGeom->GetCrystalSupportHeight() / 2.0 - ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 ;
700 gMC->Gspos("PSUP", 1, "PPAP", 0.0, yO, 0.0, 0, "ONLY") ;
703 // --- Define PIN-diode volume and position it inside crystal support ---
704 // --- right behind PbWO4 crystal
706 // --- PIN-diode dimensions ---
710 dppin[0] = fGeom->GetPinDiodeSize(0) / 2.0 ;
711 dppin[1] = fGeom->GetPinDiodeSize(1) / 2.0 ;
712 dppin[2] = fGeom->GetPinDiodeSize(2) / 2.0 ;
714 gMC->Gsvolu("PPIN", "BOX ", idtmed[705], dppin, 3) ;
716 yO = fGeom->GetCrystalSupportHeight() / 2.0 - fGeom->GetPinDiodeSize(1) / 2.0 ;
718 gMC->Gspos("PPIN", 1, "PSUP", 0.0, yO, 0.0, 0, "ONLY") ;
721 // --- Define Upper Cooling Panel, place it on top of PTCB ---
723 // --- Upper Cooling Plate thickness ---
725 dpucp[0] = dptcb[0] ;
726 dpucp[1] = fGeom->GetUpperCoolingPlateThickness() ;
727 dpucp[2] = dptcb[2] ;
729 gMC->Gsvolu("PUCP", "BOX ", idtmed[701], dpucp,3) ;
731 yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetUpperCoolingPlateThickness() ) / 2.
732 - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetModuleBoxThickness()
733 - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() - fGeom->GetUpperCoolingPlateThickness() ) ;
735 gMC->Gspos("PUCP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ;
738 // --- Define Al Support Plate, position it inside PAIR ---
739 // --- right beneath PTCB ---
740 // --- Al Support Plate thickness ---
743 dpasp[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ;
744 dpasp[1] = fGeom->GetSupportPlateThickness() / 2.0 ;
745 dpasp[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ;
747 gMC->Gsvolu("PASP", "BOX ", idtmed[701], dpasp, 3) ;
749 yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetSupportPlateThickness() ) / 2.
750 - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance()
751 - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() + dpcbl[1] * 2 ) ;
753 gMC->Gspos("PASP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ;
756 // --- Define Thermo Insulating Plate, position it inside PAIR ---
757 // --- right beneath PASP ---
758 // --- Lower Thermo Insulating Plate thickness ---
761 dptip[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ;
762 dptip[1] = fGeom->GetLowerThermoPlateThickness() / 2.0 ;
763 dptip[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ;
765 gMC->Gsvolu("PTIP", "BOX ", idtmed[706], dptip, 3) ;
767 yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetLowerThermoPlateThickness() ) / 2.
768 - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetUpperPlateThickness()
769 - fGeom->GetSecondUpperPlateThickness() + dpcbl[1] * 2 + fGeom->GetSupportPlateThickness() ) ;
771 gMC->Gspos("PTIP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ;
774 // --- Define Textolit Plate, position it inside PAIR ---
775 // --- right beneath PTIP ---
776 // --- Lower Textolit Plate thickness ---
779 dptxp[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ;
780 dptxp[1] = fGeom->GetLowerTextolitPlateThickness() / 2.0 ;
781 dptxp[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ;
783 gMC->Gsvolu("PTXP", "BOX ", idtmed[707], dptxp, 3) ;
785 yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetLowerTextolitPlateThickness() ) / 2.
786 - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetUpperPlateThickness()
787 - fGeom->GetSecondUpperPlateThickness() + dpcbl[1] * 2 + fGeom->GetSupportPlateThickness()
788 + fGeom->GetLowerThermoPlateThickness() ) ;
790 gMC->Gspos("PTXP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ;
794 //____________________________________________________________________________
795 void AliPHOSv0::CreateGeometryforPPSD()
797 // Create the PHOS-PPSD geometry for GEANT
802 Geant3 geometry tree of PHOS-PPSD in ALICE
805 <IMG Align=BOTTOM ALT="PPSD geant tree" SRC="../images/PPSDinAlice.gif">
810 // Get pointer to the array containing media indexes
811 Int_t *idtmed = fIdtmed->GetArray() - 699 ;
813 // The box containing all ppsd's for one PHOS module filled with air
815 ppsd[0] = fGeom->GetPPSDBoxSize(0) / 2.0 ;
816 ppsd[1] = fGeom->GetPPSDBoxSize(1) / 2.0 ;
817 ppsd[2] = fGeom->GetPPSDBoxSize(2) / 2.0 ;
819 gMC->Gsvolu("PPSD", "BOX ", idtmed[798], ppsd, 3) ;
821 Float_t yO = fGeom->GetOuterBoxSize(1) / 2.0 ;
823 gMC->Gspos("PPSD", 1, "PHOS", 0.0, yO, 0.0, 0, "ONLY") ;
825 // Now we build a micromegas module
826 // The box containing the whole module filled with epoxy (FR4)
829 mppsd[0] = fGeom->GetPPSDModuleSize(0) / 2.0 ;
830 mppsd[1] = fGeom->GetPPSDModuleSize(1) / 2.0 ;
831 mppsd[2] = fGeom->GetPPSDModuleSize(2) / 2.0 ;
833 gMC->Gsvolu("MPPS", "BOX ", idtmed[708], mppsd, 3) ;
836 // 1. The Top Lid made of epoxy (FR4)
839 tlppsd[0] = fGeom->GetPPSDModuleSize(0) / 2.0 ;
840 tlppsd[1] = fGeom->GetLidThickness() / 2.0 ;
841 tlppsd[2] = fGeom->GetPPSDModuleSize(2) / 2.0 ;
843 gMC->Gsvolu("TLPS", "BOX ", idtmed[708], tlppsd, 3) ;
845 Float_t y0 = ( fGeom->GetMicromegas1Thickness() - fGeom->GetLidThickness() ) / 2. ;
847 gMC->Gspos("TLPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
849 // 2. the upper panel made of composite material
852 upppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
853 upppsd[1] = fGeom->GetCompositeThickness() / 2.0 ;
854 upppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
856 gMC->Gsvolu("UPPS", "BOX ", idtmed[709], upppsd, 3) ;
858 y0 = y0 - fGeom->GetLidThickness() / 2. - fGeom->GetCompositeThickness() / 2. ;
860 gMC->Gspos("UPPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
862 // 3. the anode made of Copper
865 anppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
866 anppsd[1] = fGeom->GetAnodeThickness() / 2.0 ;
867 anppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
869 gMC->Gsvolu("ANPS", "BOX ", idtmed[710], anppsd, 3) ;
871 y0 = y0 - fGeom->GetCompositeThickness() / 2. - fGeom->GetAnodeThickness() / 2. ;
873 gMC->Gspos("ANPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
875 // 4. the conversion gap + avalanche gap filled with gas
878 ggppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
879 ggppsd[1] = ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2.0 ;
880 ggppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
882 gMC->Gsvolu("GGPS", "BOX ", idtmed[715], ggppsd, 3) ;
884 // --- Divide GGPP in X (phi) and Z directions --
885 gMC->Gsdvn("GROW", "GGPS", fGeom->GetNumberOfPadsPhi(), 1) ;
886 gMC->Gsdvn("GCEL", "GROW", fGeom->GetNumberOfPadsZ() , 3) ;
888 y0 = y0 - fGeom->GetAnodeThickness() / 2. - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. ;
890 gMC->Gspos("GGPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
893 // 6. the cathode made of Copper
896 cappsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
897 cappsd[1] = fGeom->GetCathodeThickness() / 2.0 ;
898 cappsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
900 gMC->Gsvolu("CAPS", "BOX ", idtmed[710], cappsd, 3) ;
902 y0 = y0 - ( fGeom->GetAvalancheGap() + fGeom->GetAvalancheGap() ) / 2. - fGeom->GetCathodeThickness() / 2. ;
904 gMC->Gspos("CAPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
906 // 7. the printed circuit made of G10
909 pcppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2,.0 ;
910 pcppsd[1] = fGeom->GetPCThickness() / 2.0 ;
911 pcppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
913 gMC->Gsvolu("PCPS", "BOX ", idtmed[711], cappsd, 3) ;
915 y0 = y0 - fGeom->GetCathodeThickness() / 2. - fGeom->GetPCThickness() / 2. ;
917 gMC->Gspos("PCPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
919 // 8. the lower panel made of composite material
922 lpppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
923 lpppsd[1] = fGeom->GetCompositeThickness() / 2.0 ;
924 lpppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ;
926 gMC->Gsvolu("LPPS", "BOX ", idtmed[709], lpppsd, 3) ;
928 y0 = y0 - fGeom->GetPCThickness() / 2. - fGeom->GetCompositeThickness() / 2. ;
930 gMC->Gspos("LPPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ;
932 // Position the fNumberOfModulesPhi x fNumberOfModulesZ modules (mppsd) inside PPSD to cover a PHOS module
933 // the top and bottom one's (which are assumed identical) :
935 Float_t yt = ( fGeom->GetPPSDBoxSize(1) - fGeom->GetMicromegas1Thickness() ) / 2. ;
936 Float_t yb = - ( fGeom->GetPPSDBoxSize(1) - fGeom->GetMicromegas2Thickness() ) / 2. ;
938 Int_t copyNumbertop = 0 ;
939 Int_t copyNumberbot = fGeom->GetNumberOfModulesPhi() * fGeom->GetNumberOfModulesZ() ;
941 Float_t x = ( fGeom->GetPPSDBoxSize(0) - fGeom->GetPPSDModuleSize(0) ) / 2. ;
943 for ( Int_t iphi = 1; iphi <= fGeom->GetNumberOfModulesPhi(); iphi++ ) { // the number of micromegas modules in phi per PHOS module
944 Float_t z = ( fGeom->GetPPSDBoxSize(2) - fGeom->GetPPSDModuleSize(2) ) / 2. ;
946 for ( Int_t iz = 1; iz <= fGeom->GetNumberOfModulesZ(); iz++ ) { // the number of micromegas modules in z per PHOS module
947 gMC->Gspos("MPPS", ++copyNumbertop, "PPSD", x, yt, z, 0, "ONLY") ;
948 gMC->Gspos("MPPS", ++copyNumberbot, "PPSD", x, yb, z, 0, "ONLY") ;
949 z = z - fGeom->GetPPSDModuleSize(2) ;
950 } // end of Z module loop
951 x = x - fGeom->GetPPSDModuleSize(0) ;
952 } // end of phi module loop
954 // The Lead converter between two air gaps
958 uappsd[0] = fGeom->GetPPSDBoxSize(0) / 2.0 ;
959 uappsd[1] = fGeom->GetMicro1ToLeadGap() / 2.0 ;
960 uappsd[2] = fGeom->GetPPSDBoxSize(2) / 2.0 ;
962 gMC->Gsvolu("UAPPSD", "BOX ", idtmed[798], uappsd, 3) ;
964 y0 = ( fGeom->GetPPSDBoxSize(1) - 2 * fGeom->GetMicromegas1Thickness() - fGeom->GetMicro1ToLeadGap() ) / 2. ;
966 gMC->Gspos("UAPPSD", 1, "PPSD", 0.0, y0, 0.0, 0, "ONLY") ;
971 lcppsd[0] = fGeom->GetPPSDBoxSize(0) / 2.0 ;
972 lcppsd[1] = fGeom->GetLeadConverterThickness() / 2.0 ;
973 lcppsd[2] = fGeom->GetPPSDBoxSize(2) / 2.0 ;
975 gMC->Gsvolu("LCPPSD", "BOX ", idtmed[712], lcppsd, 3) ;
977 y0 = y0 - fGeom->GetMicro1ToLeadGap() / 2. - fGeom->GetLeadConverterThickness() / 2. ;
979 gMC->Gspos("LCPPSD", 1, "PPSD", 0.0, y0, 0.0, 0, "ONLY") ;
984 lappsd[0] = fGeom->GetPPSDBoxSize(0) / 2.0 ;
985 lappsd[1] = fGeom->GetLeadToMicro2Gap() / 2.0 ;
986 lappsd[2] = fGeom->GetPPSDBoxSize(2) / 2.0 ;
988 gMC->Gsvolu("LAPPSD", "BOX ", idtmed[798], lappsd, 3) ;
990 y0 = y0 - fGeom->GetLeadConverterThickness() / 2. - fGeom->GetLeadToMicro2Gap() / 2. ;
992 gMC->Gspos("LAPPSD", 1, "PPSD", 0.0, y0, 0.0, 0, "ONLY") ;
997 //____________________________________________________________________________
998 void AliPHOSv0::Init(void)
1000 // Just prints an information message
1005 for(i=0;i<35;i++) printf("*");
1006 printf(" PHOS_INIT ");
1007 for(i=0;i<35;i++) printf("*");
1010 // Here the PHOS initialisation code (if any!)
1012 for(i=0;i<80;i++) printf("*");