]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STRUCT/AliPIPEv0.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Add SetOwner() for output lists
[u/mrichter/AliRoot.git] / STRUCT / AliPIPEv0.cxx
diff --git a/STRUCT/AliPIPEv0.cxx b/STRUCT/AliPIPEv0.cxx
index 1b8b7762e3c852aa9f291c361764e71287006f44..23318df0d5c24d734dd95e5ead01b1f45db7a01f 100644 (file)
--- a/STRUCT/AliPIPEv0.cxx
+++ b/STRUCT/AliPIPEv0.cxx
@@ -23,6 +23,7 @@
 
 #include <Riostream.h>
 
+#include <TGeoGlobalMagField.h>
 #include <TSystem.h>
 #include <TVirtualMC.h>
 
@@ -30,22 +31,24 @@
 #include "AliMagF.h"
 #include "AliPIPEv0.h"
 #include "AliRun.h"
+#include "AliLog.h"
  
 ClassImp(AliPIPEv0)
+
  
 //_____________________________________________________________________________
-AliPIPEv0::AliPIPEv0()
+AliPIPEv0::AliPIPEv0():
+    fPipeMaterial(kBe)   
 {
 // Constructor
-    SetPipeMaterial();
 }
 
 //_____________________________________________________________________________
 AliPIPEv0::AliPIPEv0(const char *name, const char *title)
-  : AliPIPE(name,title)
+    : AliPIPE(name,title),
+      fPipeMaterial(kBe)   
 {
 // Constructor
-    SetPipeMaterial();
 }
 
  
@@ -65,7 +68,7 @@ void AliPIPEv0::CreateGeometry()
 */
 //End_Html
 
-    if(fDebug) printf("%s: Create PIPEv0 geometry \n",ClassName());
+    AliDebugClass(1,"Create PIPEv0 geometry");
   
 
     Int_t *idtmed = fIdtmed->GetArray();
@@ -92,15 +95,15 @@ void AliPIPEv0::CreateGeometry()
 //
 //  distance between bellows
 //  total size of bellow section
-    const Float_t kdzb  = 15.0;
+    const Float_t kdzb  = 14.6;
 //  size of undulated region 
 //  
 //  Absorber side
 //   
 //  distance between bellows
-    const Float_t kdzbbA =  5.0;
+    const Float_t kdzbbA =  4.6;
 //  total size of bellow section
-    const Float_t kdzbA  = 15.0;
+    const Float_t kdzbA  = 14.6;
 //  size of undulated region 
     const Float_t kdzubA =  3.75;
 
@@ -113,8 +116,8 @@ void AliPIPEv0::CreateGeometry()
 //
 //    Float_t hlenQbt1 = 5.5/2.;
 //
-//  Pipe outside central region (non-absober side)
-    Float_t hlenQbab = 157.5;
+//  Pipe outside central region (non-absorber side)
+    Float_t hlenQbab = 157.5 + 23./2.;
 //
 //  Flange non-absorber side
     Float_t hlenQb29 = 11.5/2.+1.75 + 5.0;
@@ -123,7 +126,7 @@ void AliPIPEv0::CreateGeometry()
     Float_t hlenQbe0 = kdzbA;
 //
 //  Inox pipe between Be and Bellow (absorber side)
-    Float_t hlenQb24[3] = {10.5/2., 1.8, 3.3};
+    Float_t hlenQb24[3] = {11.3/2., 1.8, 3.3};
 //
 //
     Float_t hlenQb28 = (800.-hlenQbbe1-2.*hlenQbab-4.*hlenQb29-2.*hlenQbe0)/2.;
@@ -349,8 +352,15 @@ void AliPIPEv0::CreateGeometry()
     ptube[2] =   hlenQb28;    
 
     gMC->Gsvolu("QB28","TUBE", idtmed[kInox], ptube, 3);
+//
+//
+// Air with high transport cuts outside QB28
+    ptube[0] =   25.;
+    ptube[1] =   100.;
+    ptube[2] =   hlenQb28;    
 
-
+    gMC->Gsvolu("QA28","TUBE", idtmed[kAirHigh], ptube, 3);
+    gMC->Gsatt("QA28", "SEEN", 0);
 //  Al-Be (40-60 wgt%, rho=2.7 g/cm**3) beam pipe
 //
 //  This section is under study (A.M. 1/2/2002)
@@ -536,6 +546,7 @@ void AliPIPEv0::CreateGeometry()
     //last inox section till 800 cm
     zpos = zpos + hlenQb29 + hlenQb28;
     gMC->Gspos("QB28", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY"); 
+    gMC->Gspos("QA28", 1, "ALIC", 0.0, 0.0, zpos, 0, "ONLY"); 
     
 //******** end of placement on non-absorber side *********
     //
@@ -670,12 +681,12 @@ void AliPIPEv0::CreateGeometry()
 // 
 // Connecting tube ->
     ptube[0] =  0.0;
-    ptube[1] =  4.9;
+    ptube[1] =  4.5;
     ptube[2] = 14.6;
     gMC->Gsvolu("QI34","TUBE", idtmed[kInox], ptube, 3);
     
     ptube[0] =  0.0;
-    ptube[1] =  4.3;
+    ptube[1] =  3.9;
     ptube[2] = 14.6;
     gMC->Gsvolu("QI44","TUBE", idtmed[kAir], ptube, 3);
     gMC->Gspos("QI44", 1, "QI34", 0.0, 0.0, 0.0, 0, "ONLY");
@@ -683,7 +694,7 @@ void AliPIPEv0::CreateGeometry()
 
   //
   // Flange ->
-    ptube[0] =  5.41;
+    ptube[0] =  4.6;
     ptube[1] =  7.30;
     ptube[2] =  2.15;
     gMC->Gsvolu("QI35","TUBE", idtmed[kInox], ptube, 3);
@@ -706,9 +717,9 @@ void AliPIPEv0::CreateMaterials()
   // Define materials for beam pipe
   //
 
-  if(fDebug) printf("%s: Create PIPEv0 materials \n",ClassName());
-  Int_t   isxfld = gAlice->Field()->Integ();
-  Float_t sxmgmx = gAlice->Field()->Max();
+  AliDebugClass(1,"Create PIPEv0 materials");
+  Int_t   isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
+  Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
   // Steel (Inox)  
   Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
   Float_t zsteel[4] = { 26.,24.,28.,14. };
@@ -729,7 +740,14 @@ void AliPIPEv0::CreateMaterials()
   Float_t zAir[4]={6.,7.,8.,18.};
   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
   Float_t dAir = 1.20479E-3;
-
+  Float_t dAir1 = 1.20479E-10;
+  //
+  // Kapton
+  //
+  Float_t aKapton[4]={1.00794,12.0107, 14.010,15.9994};
+  Float_t zKapton[4]={1.,6.,7.,8.};
+  Float_t wKapton[4]={0.026362,0.69113,0.07327,0.209235};
+  Float_t dKapton = 1.42;
   //
   //     Berillium 
   AliMaterial(5, "BERILLIUM$", 9.01, 4., 1.848, 35.3, 36.7);
@@ -741,10 +759,11 @@ void AliPIPEv0::CreateMaterials()
   AliMaterial(9,  "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
   //
   //     Air 
-  AliMixture(15, "AIR$      ", aAir, zAir, dAir, 4, wAir);
+  AliMixture(15, "AIR$",      aAir, zAir, dAir, 4, wAir);
+  AliMixture(35, "AIR_HIGH$", aAir, zAir, dAir, 4, wAir);
   //
   //     Vacuum 
-  AliMaterial(16, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
+  AliMixture(16, "VACUUM$ ", aAir, zAir, dAir1, 4, wAir);
   //
   //     stainless Steel 
   AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
@@ -758,7 +777,9 @@ void AliPIPEv0::CreateMaterials()
   //   
   AliMixture(22, "PA$", aPA, zPA, 1.14, -4, wPA);
   //
-
+  //     Kapton
+  AliMixture(23, "KAPTON", aKapton, zKapton, dKapton, 4, wKapton);
+  //
   // **************** 
   //     Defines tracking media parameters. 
   //
@@ -780,7 +801,8 @@ void AliPIPEv0::CreateMaterials()
   AliMedium(9, "ALU",      9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   //
   //    Air 
-  AliMedium(15, "AIR",    15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(15, "AIR",     15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(35, "AIR_HIFG",35, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   //
   //    Vacuum 
   AliMedium(16, "VACUUM", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
@@ -796,6 +818,9 @@ void AliPIPEv0::CreateMaterials()
   //
   //   Polyamid
   AliMedium(22, "PA"  ,   22, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  //
+  //   KAPTON
+  AliMedium(23, "KAPTON", 23, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
 }