/**************************************************************************
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* *
* Author: The ALICE Off-line Project. *
* Contributors are mentioned in the code where appropriate. *
* *
* Permission to use, copy, modify and distribute this software and its *
* documentation strictly for non-commercial purposes is hereby granted *
* without fee, provided that the above copyright notice appears in all *
* copies and that both the copyright notice and this permission notice *
* appear in the supporting documentation. The authors make no claims *
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
/* $Id$ */
//_________________________________________________________________________
// Implementation version v0 of PHOS Manager class
// Layout EMC + PPSD has name GPS2
// Layout EMC + CPV has name IHEP
// An object of this class does not produce hits nor digits
// It is the one to use if you do not want to produce outputs in TREEH or TREED
//
//*-- Author: Yves Schutz (SUBATECH)
// --- ROOT system ---
#include "TBRIK.h"
#include "TNode.h"
#include "TRandom.h"
#include "TGeometry.h"
// --- Standard library ---
#include
PHOS in ALICE displayed by root
CPV perspective view | CPV front view |
*/ //END_HTML // Get pointer to the array containing media indexes Int_t *idtmed = fIdtmed->GetArray() - 699 ; // --- // --- Define PHOS box volume, fPUFPill with thermo insulating foam --- // --- Foam Thermo Insulating outer cover dimensions --- // --- Put it in bigbox = PHOS Float_t dphos[3] ; dphos[0] = fGeom->GetOuterBoxSize(0) / 2.0 ; dphos[1] = fGeom->GetOuterBoxSize(1) / 2.0 ; dphos[2] = fGeom->GetOuterBoxSize(2) / 2.0 ; gMC->Gsvolu("EMCA", "BOX ", idtmed[706], dphos, 3) ; Float_t yO = - fGeom->GetCPVBoxSize(1) / 2.0 ; gMC->Gspos("EMCA", 1, "PHOS", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Textolit Wall box, position inside EMCA --- // --- Textolit Wall box dimentions --- Float_t dptxw[3]; dptxw[0] = fGeom->GetTextolitBoxSize(0) / 2.0 ; dptxw[1] = fGeom->GetTextolitBoxSize(1) / 2.0 ; dptxw[2] = fGeom->GetTextolitBoxSize(2) / 2.0 ; gMC->Gsvolu("PTXW", "BOX ", idtmed[707], dptxw, 3); yO = ( fGeom->GetOuterBoxThickness(1) - fGeom->GetUpperPlateThickness() ) / 2. ; gMC->Gspos("PTXW", 1, "EMCA", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Upper Polystyrene Foam Plate, place inside PTXW --- // --- immediately below Foam Thermo Insulation Upper plate --- // --- Upper Polystyrene Foam plate thickness --- Float_t dpufp[3] ; dpufp[0] = fGeom->GetTextolitBoxSize(0) / 2.0 ; dpufp[1] = fGeom->GetSecondUpperPlateThickness() / 2. ; dpufp[2] = fGeom->GetTextolitBoxSize(2) /2.0 ; gMC->Gsvolu("PUFP", "BOX ", idtmed[703], dpufp, 3) ; yO = ( fGeom->GetTextolitBoxSize(1) - fGeom->GetSecondUpperPlateThickness() ) / 2.0 ; gMC->Gspos("PUFP", 1, "PTXW", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define air-filled box, place inside PTXW --- // --- Inner AIR volume dimensions --- Float_t dpair[3] ; dpair[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ; dpair[1] = fGeom->GetAirFilledBoxSize(1) / 2.0 ; dpair[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ; gMC->Gsvolu("PAIR", "BOX ", idtmed[798], dpair, 3) ; yO = ( fGeom->GetTextolitBoxSize(1) - fGeom->GetAirFilledBoxSize(1) ) / 2.0 - fGeom->GetSecondUpperPlateThickness() ; gMC->Gspos("PAIR", 1, "PTXW", 0.0, yO, 0.0, 0, "ONLY") ; // --- Dimensions of PbWO4 crystal --- Float_t xtlX = fGeom->GetCrystalSize(0) ; Float_t xtlY = fGeom->GetCrystalSize(1) ; Float_t xtlZ = fGeom->GetCrystalSize(2) ; Float_t dptcb[3] ; dptcb[0] = fGeom->GetNPhi() * ( xtlX + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 + fGeom->GetModuleBoxThickness() ; dptcb[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 + fGeom->GetModuleBoxThickness() / 2.0 ; dptcb[2] = fGeom->GetNZ() * ( xtlZ + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 + fGeom->GetModuleBoxThickness() ; gMC->Gsvolu("PTCB", "BOX ", idtmed[706], dptcb, 3) ; yO = fGeom->GetAirFilledBoxSize(1) / 2.0 - dptcb[1] - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetModuleBoxThickness() - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() ) ; gMC->Gspos("PTCB", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Crystal BLock filled with air, position it inside PTCB --- Float_t dpcbl[3] ; dpcbl[0] = fGeom->GetNPhi() * ( xtlX + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 ; dpcbl[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 ; dpcbl[2] = fGeom->GetNZ() * ( xtlZ + 2 * fGeom->GetGapBetweenCrystals() ) / 2.0 ; gMC->Gsvolu("PCBL", "BOX ", idtmed[798], dpcbl, 3) ; // --- Divide PCBL in X (phi) and Z directions -- gMC->Gsdvn("PROW", "PCBL", Int_t (fGeom->GetNPhi()), 1) ; gMC->Gsdvn("PCEL", "PROW", Int_t (fGeom->GetNZ()), 3) ; yO = -fGeom->GetModuleBoxThickness() / 2.0 ; gMC->Gspos("PCBL", 1, "PTCB", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define STeel (actually, it's titanium) Cover volume, place inside PCEL Float_t dpstc[3] ; dpstc[0] = ( xtlX + 2 * fGeom->GetCrystalWrapThickness() ) / 2.0 ; dpstc[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 ; dpstc[2] = ( xtlZ + 2 * fGeom->GetCrystalWrapThickness() + 2 * fGeom->GetCrystalHolderThickness() ) / 2.0 ; gMC->Gsvolu("PSTC", "BOX ", idtmed[704], dpstc, 3) ; gMC->Gspos("PSTC", 1, "PCEL", 0.0, 0.0, 0.0, 0, "ONLY") ; // --- // --- Define Tyvek volume, place inside PSTC --- Float_t dppap[3] ; dppap[0] = xtlX / 2.0 + fGeom->GetCrystalWrapThickness() ; dppap[1] = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 ; dppap[2] = xtlZ / 2.0 + fGeom->GetCrystalWrapThickness() ; gMC->Gsvolu("PPAP", "BOX ", idtmed[702], dppap, 3) ; yO = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 - ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() + fGeom->GetCrystalHolderThickness() ) / 2.0 ; gMC->Gspos("PPAP", 1, "PSTC", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define PbWO4 crystal volume, place inside PPAP --- Float_t dpxtl[3] ; dpxtl[0] = xtlX / 2.0 ; dpxtl[1] = xtlY / 2.0 ; dpxtl[2] = xtlZ / 2.0 ; gMC->Gsvolu("PXTL", "BOX ", idtmed[699], dpxtl, 3) ; yO = ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 - xtlY / 2.0 - fGeom->GetCrystalWrapThickness() ; gMC->Gspos("PXTL", 1, "PPAP", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define crystal support volume, place inside PPAP --- Float_t dpsup[3] ; dpsup[0] = xtlX / 2.0 + fGeom->GetCrystalWrapThickness() ; dpsup[1] = fGeom->GetCrystalSupportHeight() / 2.0 ; dpsup[2] = xtlZ / 2.0 + fGeom->GetCrystalWrapThickness() ; gMC->Gsvolu("PSUP", "BOX ", idtmed[798], dpsup, 3) ; yO = fGeom->GetCrystalSupportHeight() / 2.0 - ( xtlY + fGeom->GetCrystalSupportHeight() + fGeom->GetCrystalWrapThickness() ) / 2.0 ; gMC->Gspos("PSUP", 1, "PPAP", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define PIN-diode volume and position it inside crystal support --- // --- right behind PbWO4 crystal // --- PIN-diode dimensions --- Float_t dppin[3] ; dppin[0] = fGeom->GetPinDiodeSize(0) / 2.0 ; dppin[1] = fGeom->GetPinDiodeSize(1) / 2.0 ; dppin[2] = fGeom->GetPinDiodeSize(2) / 2.0 ; gMC->Gsvolu("PPIN", "BOX ", idtmed[705], dppin, 3) ; yO = fGeom->GetCrystalSupportHeight() / 2.0 - fGeom->GetPinDiodeSize(1) / 2.0 ; gMC->Gspos("PPIN", 1, "PSUP", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Upper Cooling Panel, place it on top of PTCB --- Float_t dpucp[3] ; // --- Upper Cooling Plate thickness --- dpucp[0] = dptcb[0] ; dpucp[1] = fGeom->GetUpperCoolingPlateThickness() ; dpucp[2] = dptcb[2] ; gMC->Gsvolu("PUCP", "BOX ", idtmed[701], dpucp,3) ; yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetUpperCoolingPlateThickness() ) / 2. - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetModuleBoxThickness() - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() - fGeom->GetUpperCoolingPlateThickness() ) ; gMC->Gspos("PUCP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Al Support Plate, position it inside PAIR --- // --- right beneath PTCB --- // --- Al Support Plate thickness --- Float_t dpasp[3] ; dpasp[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ; dpasp[1] = fGeom->GetSupportPlateThickness() / 2.0 ; dpasp[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ; gMC->Gsvolu("PASP", "BOX ", idtmed[701], dpasp, 3) ; yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetSupportPlateThickness() ) / 2. - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() + dpcbl[1] * 2 ) ; gMC->Gspos("PASP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Thermo Insulating Plate, position it inside PAIR --- // --- right beneath PASP --- // --- Lower Thermo Insulating Plate thickness --- Float_t dptip[3] ; dptip[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ; dptip[1] = fGeom->GetLowerThermoPlateThickness() / 2.0 ; dptip[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ; gMC->Gsvolu("PTIP", "BOX ", idtmed[706], dptip, 3) ; yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetLowerThermoPlateThickness() ) / 2. - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() + dpcbl[1] * 2 + fGeom->GetSupportPlateThickness() ) ; gMC->Gspos("PTIP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ; // --- // --- Define Textolit Plate, position it inside PAIR --- // --- right beneath PTIP --- // --- Lower Textolit Plate thickness --- Float_t dptxp[3] ; dptxp[0] = fGeom->GetAirFilledBoxSize(0) / 2.0 ; dptxp[1] = fGeom->GetLowerTextolitPlateThickness() / 2.0 ; dptxp[2] = fGeom->GetAirFilledBoxSize(2) / 2.0 ; gMC->Gsvolu("PTXP", "BOX ", idtmed[707], dptxp, 3) ; yO = ( fGeom->GetAirFilledBoxSize(1) - fGeom->GetLowerTextolitPlateThickness() ) / 2. - ( fGeom->GetIPtoCrystalSurface() - fGeom->GetIPtoOuterCoverDistance() - fGeom->GetUpperPlateThickness() - fGeom->GetSecondUpperPlateThickness() + dpcbl[1] * 2 + fGeom->GetSupportPlateThickness() + fGeom->GetLowerThermoPlateThickness() ) ; gMC->Gspos("PTXP", 1, "PAIR", 0.0, yO, 0.0, 0, "ONLY") ; } //____________________________________________________________________________ void AliPHOSv0::CreateGeometryforPPSD() { // Create the PHOS-PPSD geometry for GEANT //BEGIN_HTML /*
*/ //END_HTML // Get pointer to the array containing media indexes Int_t *idtmed = fIdtmed->GetArray() - 699 ; // The box containing all ppsd's for one PHOS module filled with air Float_t ppsd[3] ; ppsd[0] = fGeom->GetCPVBoxSize(0) / 2.0 ; ppsd[1] = fGeom->GetCPVBoxSize(1) / 2.0 ; ppsd[2] = fGeom->GetCPVBoxSize(2) / 2.0 ; gMC->Gsvolu("PPSD", "BOX ", idtmed[798], ppsd, 3) ; Float_t yO = fGeom->GetOuterBoxSize(1) / 2.0 ; gMC->Gspos("PPSD", 1, "PHOS", 0.0, yO, 0.0, 0, "ONLY") ; // Now we build a micromegas module // The box containing the whole module filled with epoxy (FR4) Float_t mppsd[3] ; mppsd[0] = fGeom->GetPPSDModuleSize(0) / 2.0 ; mppsd[1] = fGeom->GetPPSDModuleSize(1) / 2.0 ; mppsd[2] = fGeom->GetPPSDModuleSize(2) / 2.0 ; gMC->Gsvolu("MPPS", "BOX ", idtmed[708], mppsd, 3) ; // Inside mppsd : // 1. The Top Lid made of epoxy (FR4) Float_t tlppsd[3] ; tlppsd[0] = fGeom->GetPPSDModuleSize(0) / 2.0 ; tlppsd[1] = fGeom->GetLidThickness() / 2.0 ; tlppsd[2] = fGeom->GetPPSDModuleSize(2) / 2.0 ; gMC->Gsvolu("TLPS", "BOX ", idtmed[708], tlppsd, 3) ; Float_t y0 = ( fGeom->GetMicromegas1Thickness() - fGeom->GetLidThickness() ) / 2. ; gMC->Gspos("TLPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // 2. the upper panel made of composite material Float_t upppsd[3] ; upppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; upppsd[1] = fGeom->GetCompositeThickness() / 2.0 ; upppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; gMC->Gsvolu("UPPS", "BOX ", idtmed[709], upppsd, 3) ; y0 = y0 - fGeom->GetLidThickness() / 2. - fGeom->GetCompositeThickness() / 2. ; gMC->Gspos("UPPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // 3. the anode made of Copper Float_t anppsd[3] ; anppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; anppsd[1] = fGeom->GetAnodeThickness() / 2.0 ; anppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; gMC->Gsvolu("ANPS", "BOX ", idtmed[710], anppsd, 3) ; y0 = y0 - fGeom->GetCompositeThickness() / 2. - fGeom->GetAnodeThickness() / 2. ; gMC->Gspos("ANPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // 4. the conversion gap + avalanche gap filled with gas Float_t ggppsd[3] ; ggppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; ggppsd[1] = ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2.0 ; ggppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; gMC->Gsvolu("GGPS", "BOX ", idtmed[715], ggppsd, 3) ; // --- Divide GGPP in X (phi) and Z directions -- gMC->Gsdvn("GROW", "GGPS", fGeom->GetNumberOfPadsPhi(), 1) ; gMC->Gsdvn("GCEL", "GROW", fGeom->GetNumberOfPadsZ() , 3) ; y0 = y0 - fGeom->GetAnodeThickness() / 2. - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. ; gMC->Gspos("GGPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // 6. the cathode made of Copper Float_t cappsd[3] ; cappsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; cappsd[1] = fGeom->GetCathodeThickness() / 2.0 ; cappsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; gMC->Gsvolu("CAPS", "BOX ", idtmed[710], cappsd, 3) ; y0 = y0 - ( fGeom->GetConversionGap() + fGeom->GetAvalancheGap() ) / 2. - fGeom->GetCathodeThickness() / 2. ; gMC->Gspos("CAPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // 7. the printed circuit made of G10 Float_t pcppsd[3] ; pcppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2,.0 ; pcppsd[1] = fGeom->GetPCThickness() / 2.0 ; pcppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; gMC->Gsvolu("PCPS", "BOX ", idtmed[711], cappsd, 3) ; y0 = y0 - fGeom->GetCathodeThickness() / 2. - fGeom->GetPCThickness() / 2. ; gMC->Gspos("PCPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // 8. the lower panel made of composite material Float_t lpppsd[3] ; lpppsd[0] = ( fGeom->GetPPSDModuleSize(0) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; lpppsd[1] = fGeom->GetCompositeThickness() / 2.0 ; lpppsd[2] = ( fGeom->GetPPSDModuleSize(2) - fGeom->GetMicromegasWallThickness() ) / 2.0 ; gMC->Gsvolu("LPPS", "BOX ", idtmed[709], lpppsd, 3) ; y0 = y0 - fGeom->GetPCThickness() / 2. - fGeom->GetCompositeThickness() / 2. ; gMC->Gspos("LPPS", 1, "MPPS", 0.0, y0, 0.0, 0, "ONLY") ; // Position the fNumberOfModulesPhi x fNumberOfModulesZ modules (mppsd) inside PPSD to cover a PHOS module // the top and bottom one's (which are assumed identical) : Float_t yt = ( fGeom->GetCPVBoxSize(1) - fGeom->GetMicromegas1Thickness() ) / 2. ; Float_t yb = - ( fGeom->GetCPVBoxSize(1) - fGeom->GetMicromegas2Thickness() ) / 2. ; Int_t copyNumbertop = 0 ; Int_t copyNumberbot = fGeom->GetNumberOfModulesPhi() * fGeom->GetNumberOfModulesZ() ; Float_t x = ( fGeom->GetCPVBoxSize(0) - fGeom->GetPPSDModuleSize(0) ) / 2. ; for ( Int_t iphi = 1; iphi <= fGeom->GetNumberOfModulesPhi(); iphi++ ) { // the number of micromegas modules in phi per PHOS module Float_t z = ( fGeom->GetCPVBoxSize(2) - fGeom->GetPPSDModuleSize(2) ) / 2. ; for ( Int_t iz = 1; iz <= fGeom->GetNumberOfModulesZ(); iz++ ) { // the number of micromegas modules in z per PHOS module gMC->Gspos("MPPS", ++copyNumbertop, "PPSD", x, yt, z, 0, "ONLY") ; gMC->Gspos("MPPS", ++copyNumberbot, "PPSD", x, yb, z, 0, "ONLY") ; z = z - fGeom->GetPPSDModuleSize(2) ; } // end of Z module loop x = x - fGeom->GetPPSDModuleSize(0) ; } // end of phi module loop // The Lead converter between two air gaps // 1. Upper air gap Float_t uappsd[3] ; uappsd[0] = fGeom->GetCPVBoxSize(0) / 2.0 ; uappsd[1] = fGeom->GetMicro1ToLeadGap() / 2.0 ; uappsd[2] = fGeom->GetCPVBoxSize(2) / 2.0 ; gMC->Gsvolu("UAPPSD", "BOX ", idtmed[798], uappsd, 3) ; y0 = ( fGeom->GetCPVBoxSize(1) - 2 * fGeom->GetMicromegas1Thickness() - fGeom->GetMicro1ToLeadGap() ) / 2. ; gMC->Gspos("UAPPSD", 1, "PPSD", 0.0, y0, 0.0, 0, "ONLY") ; // 2. Lead converter Float_t lcppsd[3] ; lcppsd[0] = fGeom->GetCPVBoxSize(0) / 2.0 ; lcppsd[1] = fGeom->GetLeadConverterThickness() / 2.0 ; lcppsd[2] = fGeom->GetCPVBoxSize(2) / 2.0 ; gMC->Gsvolu("LCPPSD", "BOX ", idtmed[712], lcppsd, 3) ; y0 = y0 - fGeom->GetMicro1ToLeadGap() / 2. - fGeom->GetLeadConverterThickness() / 2. ; gMC->Gspos("LCPPSD", 1, "PPSD", 0.0, y0, 0.0, 0, "ONLY") ; // 3. Lower air gap Float_t lappsd[3] ; lappsd[0] = fGeom->GetCPVBoxSize(0) / 2.0 ; lappsd[1] = fGeom->GetLeadToMicro2Gap() / 2.0 ; lappsd[2] = fGeom->GetCPVBoxSize(2) / 2.0 ; gMC->Gsvolu("LAPPSD", "BOX ", idtmed[798], lappsd, 3) ; y0 = y0 - fGeom->GetLeadConverterThickness() / 2. - fGeom->GetLeadToMicro2Gap() / 2. ; gMC->Gspos("LAPPSD", 1, "PPSD", 0.0, y0, 0.0, 0, "ONLY") ; } //____________________________________________________________________________ void AliPHOSv0::CreateGeometryforCPV() { // Create the PHOS-CPV geometry for GEANT // Author: Yuri Kharlov 11 September 2000 //BEGIN_HTML /*
CPV perspective view | CPV front view |
One CPV module, perspective view | One CPV module, front view (extended in vertical direction) |
*/ //END_HTML Float_t par[5], x0,y0,z0 ; Int_t i,j,copy; // Get pointer to the array containing media indexes Int_t *idtmed = fIdtmed->GetArray() - 699 ; // --- Dummy box containing two rails on which PHOS support moves // --- Put these rails to the bottom of the L3 magnet par[0] = fGeom->GetRailRoadSize(0) / 2.0 ; par[1] = fGeom->GetRailRoadSize(1) / 2.0 ; par[2] = fGeom->GetRailRoadSize(2) / 2.0 ; gMC->Gsvolu("PRRD", "BOX ", idtmed[798], par, 3) ; y0 = -(fGeom->GetRailsDistanceFromIP() - fGeom->GetRailRoadSize(1) / 2.0) ; gMC->Gspos("PRRD", 1, "ALIC", 0.0, y0, 0.0, 0, "ONLY") ; // --- Dummy box containing one rail par[0] = fGeom->GetRailOuterSize(0) / 2.0 ; par[1] = fGeom->GetRailOuterSize(1) / 2.0 ; par[2] = fGeom->GetRailOuterSize(2) / 2.0 ; gMC->Gsvolu("PRAI", "BOX ", idtmed[798], par, 3) ; for (i=0; i<2; i++) { x0 = (2*i-1) * fGeom->GetDistanceBetwRails() / 2.0 ; gMC->Gspos("PRAI", i, "PRRD", x0, 0.0, 0.0, 0, "ONLY") ; } // --- Upper and bottom steel parts of the rail par[0] = fGeom->GetRailPart1(0) / 2.0 ; par[1] = fGeom->GetRailPart1(1) / 2.0 ; par[2] = fGeom->GetRailPart1(2) / 2.0 ; gMC->Gsvolu("PRP1", "BOX ", idtmed[716], par, 3) ; y0 = - (fGeom->GetRailOuterSize(1) - fGeom->GetRailPart1(1)) / 2.0 ; gMC->Gspos("PRP1", 1, "PRAI", 0.0, y0, 0.0, 0, "ONLY") ; y0 = (fGeom->GetRailOuterSize(1) - fGeom->GetRailPart1(1)) / 2.0 - fGeom->GetRailPart3(1); gMC->Gspos("PRP1", 2, "PRAI", 0.0, y0, 0.0, 0, "ONLY") ; // --- The middle vertical steel parts of the rail par[0] = fGeom->GetRailPart2(0) / 2.0 ; par[1] = fGeom->GetRailPart2(1) / 2.0 ; par[2] = fGeom->GetRailPart2(2) / 2.0 ; gMC->Gsvolu("PRP2", "BOX ", idtmed[716], par, 3) ; y0 = - fGeom->GetRailPart3(1) / 2.0 ; gMC->Gspos("PRP2", 1, "PRAI", 0.0, y0, 0.0, 0, "ONLY") ; // --- The most upper steel parts of the rail par[0] = fGeom->GetRailPart3(0) / 2.0 ; par[1] = fGeom->GetRailPart3(1) / 2.0 ; par[2] = fGeom->GetRailPart3(2) / 2.0 ; gMC->Gsvolu("PRP3", "BOX ", idtmed[716], par, 3) ; y0 = (fGeom->GetRailOuterSize(1) - fGeom->GetRailPart3(1)) / 2.0 ; gMC->Gspos("PRP3", 1, "PRAI", 0.0, y0, 0.0, 0, "ONLY") ; // --- The wall of the cradle // --- The wall is empty: steel thin walls and air inside par[1] = TMath::Sqrt( TMath::Power((fGeom->GetIPtoOuterCoverDistance() + fGeom->GetOuterBoxSize(1)),2) + TMath::Power((fGeom->GetOuterBoxSize(0)/2),2)) + 10.; par[0] = par[1] - fGeom->GetCradleWall(1) ; par[2] = fGeom->GetCradleWall(2) / 2.0 ; par[3] = fGeom->GetCradleWall(3) ; par[4] = fGeom->GetCradleWall(4) ; gMC->Gsvolu("PCRA", "TUBS", idtmed[716], par, 5) ; par[0] -= fGeom->GetCradleWallThickness() ; par[1] -= fGeom->GetCradleWallThickness() ; par[2] -= fGeom->GetCradleWallThickness() ; gMC->Gsvolu("PCRE", "TUBS", idtmed[798], par, 5) ; gMC->Gspos ("PCRE", 1, "PCRA", 0.0, 0.0, 0.0, 0, "ONLY") ; for (i=0; i<2; i++) { z0 = (2*i-1) * (fGeom->GetOuterBoxSize(2) + fGeom->GetCradleWall(2)) / 2.0 ; gMC->Gspos("PCRA", i, "ALIC", 0.0, 0.0, z0, 0, "ONLY") ; } // --- The "wheels" of the cradle par[0] = fGeom->GetCradleWheel(0) / 2; par[1] = fGeom->GetCradleWheel(1) / 2; par[2] = fGeom->GetCradleWheel(2) / 2; gMC->Gsvolu("PWHE", "BOX ", idtmed[716], par, 3) ; y0 = -(fGeom->GetRailsDistanceFromIP() - fGeom->GetRailRoadSize(1) - fGeom->GetCradleWheel(1)/2) ; for (i=0; i<2; i++) { z0 = (2*i-1) * ((fGeom->GetOuterBoxSize(2) + fGeom->GetCradleWheel(2)) / 2.0 + fGeom->GetCradleWall(2)); for (j=0; j<2; j++) { copy = 2*i + j; x0 = (2*j-1) * fGeom->GetDistanceBetwRails() / 2.0 ; gMC->Gspos("PWHE", copy, "ALIC", x0, y0, z0, 0, "ONLY") ; } } } //____________________________________________________________________________ void AliPHOSv0::Init(void) { // Just prints an information message Int_t i; printf("\n"); for(i=0;i<35;i++) printf("*"); printf(" PHOS_INIT "); for(i=0;i<35;i++) printf("*"); printf("\n"); // Here the PHOS initialisation code (if any!) if (fGeom!=0) cout << "AliPHOS" << Version() << " : PHOS geometry intialized for " << fGeom->GetName() << endl ; else cout << "AliPHOS" << Version() << " : PHOS geometry initialization failed !" << endl ; for(i=0;i<80;i++) printf("*"); printf("\n"); }