X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STRUCT%2FAliPIPEv1.cxx;h=75ec264151738b7f830a2bd3828513a5f4db4616;hb=b3d81bee32cc657834ce18d89b361d51b0d0fbbe;hp=20ca0fca4dee550e9b6a2dfb883faa0b1b454e26;hpb=ad51aeb0c02b3e46da92b7c5d59c38d3518352f2;p=u%2Fmrichter%2FAliRoot.git diff --git a/STRUCT/AliPIPEv1.cxx b/STRUCT/AliPIPEv1.cxx index 20ca0fca4de..75ec2641517 100644 --- a/STRUCT/AliPIPEv1.cxx +++ b/STRUCT/AliPIPEv1.cxx @@ -1,6 +1,23 @@ +/************************************************************************** + * 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$ */ + /////////////////////////////////////////////////////////////////////////////// // // -// Beam pipe class // +// Beam pipe class. Test version // // // //Begin_Html /* @@ -10,9 +27,11 @@ // // /////////////////////////////////////////////////////////////////////////////// +#include + +#include "AliMagF.h" #include "AliPIPEv1.h" #include "AliRun.h" -#include "AliMC.h" ClassImp(AliPIPEv1) @@ -50,24 +69,22 @@ void AliPIPEv1::CreateGeometry() */ //End_Html - AliMC* pMC = AliMC::GetMC(); - - Float_t tpar[3], dzmo, zpos, abs_d, abs_l; + Float_t tpar[3], dzmo, zpos, absorberDistance, absorberEnd; Float_t r2, dr; - const Double_t z_flange = 150; + const Double_t kZFlange = 150; Int_t *idtmed = fIdtmed->GetArray()-1999; - abs_d = 90.; // DEFINES DRIFT LENGTH + absorberDistance = 90.; // DEFINES DRIFT LENGTH //z_nose = 102.; //z_cone = 285.; //theta1 = 24.; // 1. angle defining the front absorber //theta2 = 5.; // 2. angle defining the front absorbe //acc_max = 9.; // ANGLE POLAIRE MAXIMUM //acc_min = 2.; // ANGLE POLAIRE MINIMUM DE DETECTION - abs_l = 503.; + absorberEnd = 503.; //d_steel = 1.; // THICKNESS OF STEEL SUPPORT //d_poly = 7.5; //d_pb = 2.5; @@ -109,75 +126,75 @@ void AliPIPEv1::CreateGeometry() tpar[0] = 0.; tpar[1] = 3.; - tpar[2] = (abs_d + 700.) / 2.; - dzmo = tpar[2] - abs_d; - pMC->Gsvolu("QQMO", "TUBE", idtmed[2015], tpar, 3); - pMC->Gspos("QQMO", 1, "ALIC", 0., 0., -dzmo, 0, "ONLY"); + tpar[2] = (absorberDistance + 700.) / 2.; + dzmo = tpar[2] - absorberDistance; + gMC->Gsvolu("QQMO", "TUBE", idtmed[2015], tpar, 3); + gMC->Gspos("QQMO", 1, "ALIC", 0., 0., -dzmo, 0, "ONLY"); // BEAM PIPE IN DRIFT SPACE - // -30-z_flange + // -30-kZFlange tpar[0] = 0.; tpar[1] = 3.; tpar[2] = 30; - pMC->Gsvolu("QDT1", "TUBE", idtmed[2015], tpar, 3); + gMC->Gsvolu("QDT1", "TUBE", idtmed[2015], tpar, 3); tpar[0] = 2.9; - pMC->Gsvolu("QTB1", "TUBE", idtmed[2004], tpar, 3); - pMC->Gspos("QTB1", 1, "QDT1", 0., 0., 0., 0, "ONLY"); - pMC->Gspos("QDT1", 1, "QQMO", 0., 0., dzmo, 0, "ONLY"); + gMC->Gsvolu("QTB1", "TUBE", idtmed[2004], tpar, 3); + gMC->Gspos("QTB1", 1, "QDT1", 0., 0., 0., 0, "ONLY"); + gMC->Gspos("QDT1", 1, "QQMO", 0., 0., dzmo, 0, "ONLY"); // 30-90 tpar[0] = 0.; tpar[1] = 3.; tpar[2] = 30.; - pMC->Gsvolu("QDT2", "TUBE", idtmed[2015], tpar, 3); + gMC->Gsvolu("QDT2", "TUBE", idtmed[2015], tpar, 3); tpar[0] = 2.9; - pMC->Gsvolu("QTB2", "TUBE", idtmed[2004], tpar, 3); - pMC->Gspos("QTB2", 1, "QDT2", 0., 0., 0., 0, "ONLY"); - pMC->Gspos("QDT2", 1, "QQMO", 0., 0., dzmo, 0, "ONLY"); + gMC->Gsvolu("QTB2", "TUBE", idtmed[2004], tpar, 3); + gMC->Gspos("QTB2", 1, "QDT2", 0., 0., 0., 0, "ONLY"); + gMC->Gspos("QDT2", 1, "QQMO", 0., 0., dzmo, 0, "ONLY"); // beam pipe outside absorber on the left side - // -30 - Z_FLANGE + // -30 - kZFlange tpar[0] = 0.; tpar[1] = 3.; - tpar[2] = (z_flange - 30)/2; - pMC->Gsvolu("QDT5", "TUBE", idtmed[2015], tpar, 3); + tpar[2] = (kZFlange - 30)/2; + gMC->Gsvolu("QDT5", "TUBE", idtmed[2015], tpar, 3); tpar[0] = 2.9; zpos = -30. - tpar[2] + dzmo; - pMC->Gsvolu("QTB5", "TUBE", idtmed[2004], tpar, 3); - pMC->Gspos("QTB5", 1, "QDT5", 0., 0., 0., 0, "ONLY"); - pMC->Gspos("QDT5", 1, "QQMO", 0., 0., zpos, 0, "ONLY"); + gMC->Gsvolu("QTB5", "TUBE", idtmed[2004], tpar, 3); + gMC->Gspos("QTB5", 1, "QDT5", 0., 0., 0., 0, "ONLY"); + gMC->Gspos("QDT5", 1, "QQMO", 0., 0., zpos, 0, "ONLY"); // STRAIGHT STEEL PIECE - zpos = -z_flange; + zpos = -kZFlange; r2 = 2.9; dr = .015; tpar[0] = 0.; tpar[1] = r2 + dr; tpar[2] = (zpos + 700.) / 2.; - pMC->Gsvolu("QDT7", "TUBE", idtmed[2015], tpar, 3); + gMC->Gsvolu("QDT7", "TUBE", idtmed[2015], tpar, 3); tpar[0] = r2; - pMC->Gsvolu("QTB7", "TUBE", idtmed[2018], tpar, 3); - pMC->Gspos("QTB7", 1, "QDT7", 0., 0., 0., 0, "ONLY"); + gMC->Gsvolu("QTB7", "TUBE", idtmed[2018], tpar, 3); + gMC->Gspos("QTB7", 1, "QDT7", 0., 0., 0., 0, "ONLY"); zpos = zpos - tpar[2] + dzmo; - pMC->Gspos("QDT7", 1, "QQMO", 0., 0., zpos, 0, "ONLY"); + gMC->Gspos("QDT7", 1, "QQMO", 0., 0., zpos, 0, "ONLY"); // flange dn 63 tpar[0] = 3.; tpar[1] = 5.7; tpar[2] = 2.; - pMC->Gsvolu("QN63", "TUBE", idtmed[2018], tpar, 3); - zpos = tpar[2] - z_flange; - pMC->Gspos("QN63", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); + gMC->Gsvolu("QN63", "TUBE", idtmed[2018], tpar, 3); + zpos = tpar[2] - kZFlange; + gMC->Gspos("QN63", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); // Replace Absorber or Shield by Beam-Pipe @@ -185,83 +202,44 @@ void AliPIPEv1::CreateGeometry() if (gAlice->GetModule("ABSO") == 0) { - pMC->Gspos("QN63", 2, "ALIC", 0., 0., z_flange, 0, "ONLY"); + gMC->Gspos("QN63", 2, "ALIC", 0., 0., kZFlange, 0, "ONLY"); r2 = 2.9; dr = .1; tpar[0] = 0.; tpar[1] = r2 + dr; - tpar[2] = (z_flange - abs_d) / 2.; - pMC->Gsvolu("QDT8", "TUBE", idtmed[2015], tpar, 3); + tpar[2] = (kZFlange - absorberDistance) / 2.; + gMC->Gsvolu("QDT8", "TUBE", idtmed[2015], tpar, 3); tpar[0] = r2; - pMC->Gsvolu("QTB8", "TUBE", idtmed[2004], tpar, 3); - pMC->Gspos("QTB8", 1, "QDT8", 0., 0., 0., 0, "ONLY"); - zpos = abs_d + tpar[2]; - pMC->Gspos("QDT8", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); + gMC->Gsvolu("QTB8", "TUBE", idtmed[2004], tpar, 3); + gMC->Gspos("QTB8", 1, "QDT8", 0., 0., 0., 0, "ONLY"); + zpos = absorberDistance + tpar[2]; + gMC->Gspos("QDT8", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); dr = .015; tpar[0] = 0.; tpar[1] = r2 + dr; - tpar[2] = (abs_l - z_flange) / 2.; - pMC->Gsvolu("QDTS", "TUBE", idtmed[2015], tpar, 3); + tpar[2] = (absorberEnd - kZFlange) / 2.; + gMC->Gsvolu("QDTS", "TUBE", idtmed[2015], tpar, 3); tpar[0] = r2; - pMC->Gsvolu("QTBS", "TUBE", idtmed[2018], tpar, 3); - pMC->Gspos("QTBS", 1, "QDTS", 0., 0., 0., 0, "ONLY"); - zpos = tpar[2] + z_flange; - pMC->Gspos("QDTS", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); + gMC->Gsvolu("QTBS", "TUBE", idtmed[2018], tpar, 3); + gMC->Gspos("QTBS", 1, "QDTS", 0., 0., 0., 0, "ONLY"); + zpos = tpar[2] + kZFlange; + gMC->Gspos("QDTS", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); } if (gAlice->GetModule("SHIL") == 0) { r2 = 2.9; dr = .015; tpar[0] = 0.; tpar[1] = r2 + dr; - tpar[2] = (700. - abs_l) / 2.; - pMC->Gsvolu("QDT9", "TUBE", idtmed[2015], tpar, 3); + tpar[2] = (700. - absorberEnd) / 2.; + gMC->Gsvolu("QDT9", "TUBE", idtmed[2015], tpar, 3); tpar[0] = r2; - pMC->Gsvolu("QTB9", "TUBE", idtmed[2018], tpar, 3); - pMC->Gspos("QTB9", 1, "QDT9", 0., 0., 0., 0, "ONLY"); - zpos = abs_l + tpar[2]; - pMC->Gspos("QDT9", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); + gMC->Gsvolu("QTB9", "TUBE", idtmed[2018], tpar, 3); + gMC->Gspos("QTB9", 1, "QDT9", 0., 0., 0., 0, "ONLY"); + zpos = absorberEnd + tpar[2]; + gMC->Gspos("QDT9", 1, "ALIC", 0., 0., zpos, 0, "ONLY"); } } -//_____________________________________________________________________________ -void AliPIPEv1::DrawModule() -{ - // - // Draw a shaded view of the Beam Pipe - // - - AliMC* pMC = AliMC::GetMC(); - - // Set everything unseen - pMC->Gsatt("*", "seen", -1); - // - // Set ALIC mother transparent - pMC->Gsatt("ALIC","SEEN",0); - // - // Set the volumes visible - pMC->Gsatt("QQMO","seen",1); - pMC->Gsatt("QDT1","seen",1); - pMC->Gsatt("QTB1","seen",1); - pMC->Gsatt("QDT2","seen",1); - pMC->Gsatt("QTB2","seen",1); - pMC->Gsatt("QDT5","seen",1); - pMC->Gsatt("QTB5","seen",1); - pMC->Gsatt("QDT7","seen",1); - pMC->Gsatt("QTB7","seen",1); - pMC->Gsatt("QN63","seen",1); - // - pMC->Gdopt("hide", "on"); - pMC->Gdopt("shad", "on"); - pMC->Gsatt("*", "fill", 7); - pMC->SetClipBox("."); - pMC->SetClipBox("*", 0, 3000, -3000, 3000, -6000, 6000); - pMC->DefaultRange(); - pMC->Gdraw("alic", 40, 30, 0, 3, 5, .04, .04); - pMC->Gdhead(1111, "Beam Pipe"); - pMC->Gdman(16, 6, "MAN"); - pMC->Gdopt("hide","off"); -} - //_____________________________________________________________________________ void AliPIPEv1::CreateMaterials() { @@ -269,8 +247,8 @@ void AliPIPEv1::CreateMaterials() // Create materials for beam pipe // - Int_t ISXFLD = gAlice->Field()->Integ(); - Float_t SXMGMX = gAlice->Field()->Max(); + Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ(); + Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max(); Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 }; Float_t zsteel[4] = { 26.,24.,28.,14. }; @@ -299,18 +277,18 @@ void AliPIPEv1::CreateMaterials() // Air - AliMedium(15, "AIR_L3_US", 15, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin); + AliMedium(15, "AIR_L3_US", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // Beryllium - AliMedium(5, "BE_L3_US", 5, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin); + AliMedium(5, "BE_L3_US", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // Vacuum - AliMedium(16, "VA_L3_US", 16, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin); + AliMedium(16, "VA_L3_US", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // Steel - AliMedium(19, "ST_L3_US", 19, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin); + AliMedium(19, "ST_L3_US", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); }