#include <Riostream.h>
+#include <TGeoManager.h>
+#include <TGeoGlobalMagField.h>
#include <TSystem.h>
#include <TVirtualMC.h>
+#include <TArrayI.h>
#include "AliConst.h"
#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();
}
*/
//End_Html
- if(fDebug) printf("%s: Create PIPEv0 geometry \n",ClassName());
+ AliDebugClass(1,"Create PIPEv0 geometry");
Int_t *idtmed = fIdtmed->GetArray();
//
// 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;
ppcon[61] = 0.;
ppcon[62] = 4.1;
- gMC->Gsvolu("QBPM", "PCON", idtmed[kAir], ppcon,63);
+ TVirtualMC::GetMC()->Gsvolu("QBPM", "PCON", idtmed[kAir], ppcon,63);
//
//
// The Vacuum
- gMC->Gsvolu("QBVA","TUBE", idtmed[kVac], ptube, 0);
+ TVirtualMC::GetMC()->Gsvolu("QBVA","TUBE", idtmed[kVac], ptube, 0);
ptube[0] = 0.0;
ptube[1] = kRinSt;
ptube[2] = (90.-hlenQbbe2)/2.;
dz = -90. + ptube[2];
- gMC->Gsposp ("QBVA", 1, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp ("QBVA", 1, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
dz = dz + ptube[2];
ptube[1] = kRinBe;
ptube[2] = hlenQbbe+hlenQbab;
dz = dz + ptube[2];
- gMC->Gsposp ("QBVA", 2, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp ("QBVA", 2, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
dz = dz + ptube[2];
ptube[1] = kRinSt;
ptube[2] = (800.-hlenQbbe1-2.*hlenQbab)/2.;
dz = dz + ptube[2];
- gMC->Gsposp ("QBVA", 3, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp ("QBVA", 3, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
//
// Be Pipe in central Alice
ptube[1] = kRoutBe;
ptube[2] = hlenQbbe;
- gMC->Gsvolu("QBBE","TUBE", idtmed[kBe], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QBBE","TUBE", idtmed[kBe], ptube, 3);
//
// Support Ring
ptube[0] = kRoutSB;
ptube[1] = 4.0;
ptube[2] = 0.6;
- gMC->Gsvolu("QBSR", "TUBE", idtmed[kAlu], ptube,3);
+ TVirtualMC::GetMC()->Gsvolu("QBSR", "TUBE", idtmed[kAlu], ptube,3);
// Inner support
ptube[0] = kRoutSB;
ptube[1] = 3.5;
- gMC->Gsvolu("QBSS", "TUBE", idtmed[kPA], ptube,3);
- gMC->Gspos("QBSS", 1, "QBSR", 0.0, 0.0, 0.0, 0, "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QBSS", "TUBE", idtmed[kPA], ptube,3);
+ TVirtualMC::GetMC()->Gspos("QBSS", 1, "QBSR", 0.0, 0.0, 0.0, 0, "ONLY");
- gMC->Gspos("QBSR", 1, "QBPM", 0.0, 0.0, 40., 0, "ONLY");
- gMC->Gspos("QBSR", 2, "QBPM", 0.0, 0.0, 150., 0, "ONLY");
- gMC->Gspos("QBSR", 3, "QBPM", 0.0, 0.0, 260., 0, "ONLY");
- gMC->Gspos("QBSR", 4, "QBPM", 0.0, 0.0,- 46., 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBSR", 1, "QBPM", 0.0, 0.0, 40., 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBSR", 2, "QBPM", 0.0, 0.0, 150., 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBSR", 3, "QBPM", 0.0, 0.0, 260., 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBSR", 4, "QBPM", 0.0, 0.0,- 46., 0, "ONLY");
//
// Flange and Fixed Point: non absorber side
//
ppcon[13] = kRinSt;
ppcon[14] = 3.6;
- gMC->Gsvolu("QB29", "PCON", idtmed[kAir], ppcon,15);
+ TVirtualMC::GetMC()->Gsvolu("QB29", "PCON", idtmed[kAir], ppcon,15);
// Flange
ptube[0] = kRinSt;
ptube[1] = 5.7;
ptube[2] = 1.75;
- gMC->Gsvolu("QF29","TUBE", idtmed[kInox], ptube, 3);
- gMC->Gspos("QF29", 1, "QB29", 0.0, 0.0, -hlenQb29+1.75, 0, "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QF29","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gspos("QF29", 1, "QB29", 0.0, 0.0, -hlenQb29+1.75, 0, "ONLY");
// Pipe
ptube[0] = kRinSt;
ptube[1] = 3.06;
ptube[2] = hlenQb29;
- gMC->Gsvolu("QS29","TUBE", idtmed[kInox], ptube, 3);
- gMC->Gspos("QS29", 1, "QB29", 0.0, 0.0, 0., 0, "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QS29","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gspos("QS29", 1, "QB29", 0.0, 0.0, 0., 0, "ONLY");
// Fixed point
ptube[0] = kRinSt;
ptube[1] = 3.5;
ptube[2] = 0.3;
- gMC->Gsvolu("QP29","TUBE", idtmed[kInox], ptube, 3);
- gMC->Gspos("QP29", 1, "QB29", 0.0, 0.0, -hlenQb29+9.75+3., 0, "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QP29","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gspos("QP29", 1, "QB29", 0.0, 0.0, -hlenQb29+9.75+3., 0, "ONLY");
//
//
ptube[1] = kRoutSt;
ptube[2] = hlenQb28;
- gMC->Gsvolu("QB28","TUBE", idtmed[kInox], ptube, 3);
-
+ TVirtualMC::GetMC()->Gsvolu("QB28","TUBE", idtmed[kInox], ptube, 3);
+//
+//
+// Air with high transport cuts outside QB28
+ ptube[0] = 25.;
+ ptube[1] = 100.;
+ ptube[2] = hlenQb28;
+ TVirtualMC::GetMC()->Gsvolu("QA28","TUBE", idtmed[kAirHigh], ptube, 3);
+ gGeoManager->SetVolumeAttribute("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)
}
ptube[2] = hlenQbab;
- gMC->Gsvolu("QBAB","TUBE", idtmed[fPipeMaterial], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QBAB","TUBE", idtmed[fPipeMaterial], ptube, 3);
// 2.5 mm thick SS tube for hanging pump
/*
ptube[1] = 3.15;
ptube[2] = hlenQb26;
- gMC->Gsvolu("QB26","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QB26","TUBE", idtmed[kInox], ptube, 3);
*/
//
// Bellows
pconQBE0[31] = kRinSB;
pconQBE0[32] = kRoutSB;
- gMC->Gsvolu("QBE0", "PCON", idtmed[kAir], pconQBE0, 33);
+ TVirtualMC::GetMC()->Gsvolu("QBE0", "PCON", idtmed[kAir], pconQBE0, 33);
//
// Undulated piece mother
ptube[0] = kRinSB;
ptube[1] = 3.60;
ptube[2] = kdzubA;
- gMC->Gsvolu("QBEM","TUBE", idtmed[kAir], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QBEM","TUBE", idtmed[kAir], ptube, 3);
dz = -kdzbA+kdzubA+2.5;
- gMC->Gspos("QBEM", 2 ,"QBE0", 0.0, 0.0, dz, 0 , "ONLY");
- gMC->Gspos("QBEM", 1 ,"QBE0", 0.0, 0.0, -dz, idrotm[2012], "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBEM", 2 ,"QBE0", 0.0, 0.0, dz, 0 , "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBEM", 1 ,"QBE0", 0.0, 0.0, -dz, idrotm[2012], "ONLY");
//
Float_t pund[30];
Float_t uw = 0.02;
pund[25] = pund[22];
pund[26] = pund[23];
- gMC->Gsvolu("QBEU", "PCON", idtmed[kInox], pund, 27);
+ TVirtualMC::GetMC()->Gsvolu("QBEU", "PCON", idtmed[kInox], pund, 27);
for (i = 0; i < 18; i++)
{
dz = -kdzubA+(1+2*i)*ut;
- gMC->Gspos("QBEU", i+1 ,"QBEM", 0.0, 0.0, dz, 0 , "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBEU", i+1 ,"QBEM", 0.0, 0.0, dz, 0 , "ONLY");
}
ptube[0] = kRinSB;
ptube[1] = kRinSB+uw;
ptube[2] = uz;
- gMC->Gsvolu("QBEW","TUBE", idtmed[kInox], ptube, 3);
- gMC->Gspos("QBEW", 1 ,"QBEM", 0.0, 0.0, kdzubA-uz, 0 , "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QBEW","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gspos("QBEW", 1 ,"QBEM", 0.0, 0.0, kdzubA-uz, 0 , "ONLY");
//
// BeamPipe
- gMC->Gsvolu("QBEP","TUBE", idtmed[kInox], ptube, 0);
+ TVirtualMC::GetMC()->Gsvolu("QBEP","TUBE", idtmed[kInox], ptube, 0);
ptube[0] = kRinSB;
ptube[1] = kRoutSB;
ptube[2] = 1.25;
- gMC->Gsposp("QBEP", 1 ,"QBE0", 0.0, 0.0, -kdzbA+1.25, 0 , "ONLY", ptube, 3);
- gMC->Gsposp("QBEP", 2 ,"QBE0", 0.0, 0.0, kdzbA-1.25, 0 , "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp("QBEP", 1 ,"QBE0", 0.0, 0.0, -kdzbA+1.25, 0 , "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp("QBEP", 2 ,"QBE0", 0.0, 0.0, kdzbA-1.25, 0 , "ONLY", ptube, 3);
ptube[2] = kdzbbA;
- gMC->Gsposp("QBEP", 3 ,"QBE0", 0.0, 0.0, 0., 0 , "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp("QBEP", 3 ,"QBE0", 0.0, 0.0, 0., 0 , "ONLY", ptube, 3);
//
//
// ----> End Bellow
//
// first the beryllium section
Float_t zpos = -(hlenQbbe2-hlenQbbe1)/2;
- gMC->Gspos("QBBE", 1, "QBPM", 0., 0., zpos, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBBE", 1, "QBPM", 0., 0., zpos, 0, "ONLY");
// next meta-metal transition QBT1 on on-absorber side
// zpos = zpos + hlenQbbe + hlenQbt1;
-// gMC->Gspos("QBT1", 1, "QBPM", 0., 0., zpos, 0, "ONLY");
+// TVirtualMC::GetMC()->Gspos("QBT1", 1, "QBPM", 0., 0., zpos, 0, "ONLY");
// Aluminium OR Al-be alloy section
zpos = hlenQbbe1+hlenQbab;
- gMC->Gspos("QBAB", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBAB", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
//
// inox flange at the start of bellow
zpos = zpos + hlenQbab + hlenQb29;
- gMC->Gspos("QB29", 1, "QBPM", 0.0, 0.0, zpos, idrotm[2012], "ONLY");
+ TVirtualMC::GetMC()->Gspos("QB29", 1, "QBPM", 0.0, 0.0, zpos, idrotm[2012], "ONLY");
//
// bellow section
zpos = zpos + hlenQb29 + hlenQbe0;
- gMC->Gspos("QBE0", 2 ,"QBPM", 0.0, 0.0, zpos, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBE0", 2 ,"QBPM", 0.0, 0.0, zpos, 0, "ONLY");
//
// inox flange at the end of bellow and start of thick inox for pump
zpos = zpos + hlenQbe0 + hlenQb29;
- gMC->Gspos("QB29", 2, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QB29", 2, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
//
//last inox section till 800 cm
zpos = zpos + hlenQb29 + hlenQb28;
- gMC->Gspos("QB28", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QB28", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QA28", 1, "ALIC", 0.0, 0.0, zpos, 0, "ONLY");
//******** end of placement on non-absorber side *********
//
//
// Beam pipes between elements
//
- gMC->Gsvolu("QB24","TUBE", idtmed[kInox], ptube, 0);
+ TVirtualMC::GetMC()->Gsvolu("QB24","TUBE", idtmed[kInox], ptube, 0);
ptube[0] = kRinSt;
ptube[1] = kRoutSt;
ptube[2] = hlenQb24[0];
bpbe[31] = kRinSt;
bpbe[32] = kRoutSt;
- gMC->Gsvolu("QA24","PCON", idtmed[kInox], bpbe, 33);
+ TVirtualMC::GetMC()->Gsvolu("QA24","PCON", idtmed[kInox], bpbe, 33);
dz = hlenQbbe2 + hlenQb24[0];
- gMC->Gspos("QA24", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QA24", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY");
//
// Bellow on absorber side
dz = dz+hlenQb24[0] + kdzbA;
- gMC->Gspos("QBE0", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QBE0", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY");
//
ptube[2] = hlenQb24[1];
dz = dz + kdzb + ptube[2];
- gMC->Gsposp("QB24", 2 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp("QB24", 2 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3);
dz = dz + ptube[2];
//
ptube[1] = 4.300;
ptube[2] = 1.4;
- gMC->Gsvolu("QFA0","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QFA0","TUBE", idtmed[kInox], ptube, 3);
dz = dz + ptube[2];
- gMC->Gspos("QFA0", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QFA0", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY");
dz = dz + ptube[2];
//
//
ptube[1] = kRoutSB;
ptube[2] = hlenQb24[2];
dz = dz + ptube[2];
- gMC->Gsposp("QB24", 3 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3);
+ TVirtualMC::GetMC()->Gsposp("QB24", 3 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3);
// --- Place the PIPE ghost volume (QBPM) in its mother volume (ALIC)
// by rotating it to 180 deg. and make it invisible
//
- gMC->Gspos("QBPM",1,"ALIC", 0, 0, 0, 0, "ONLY");
- gMC->Gsbool("QBPM", "L3DX");
- gMC->Gsbool("QBPM", "L3O3");
- gMC->Gsbool("QBPM", "L3O4");
+ TVirtualMC::GetMC()->Gspos("QBPM",1,"ALIC", 0, 0, 0, 0, "ONLY");
+ TVirtualMC::GetMC()->Gsbool("QBPM", "L3DX");
+ TVirtualMC::GetMC()->Gsbool("QBPM", "L3O3");
+ TVirtualMC::GetMC()->Gsbool("QBPM", "L3O4");
//
pbox[0] = 6.50;
pbox[1] = 6.75;
pbox[2] = 15.60;
- gMC->Gsvolu("QI32","BOX", idtmed[kInox], pbox, 3);
+ TVirtualMC::GetMC()->Gsvolu("QI32","BOX", idtmed[kInox], pbox, 3);
pbox[0] = 5.90;
pbox[1] = 6.15;
pbox[2] = 15.00;
- gMC->Gsvolu("QI42","BOX", idtmed[kGetter], pbox, 3);
- gMC->Gspos("QI42", 1, "QI32", 0.0, 0.0, 0.0, 0, "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QI42","BOX", idtmed[kGetter], pbox, 3);
+ TVirtualMC::GetMC()->Gspos("QI42", 1, "QI32", 0.0, 0.0, 0.0, 0, "ONLY");
// <-
ptube[0] = 0.0;
ptube[1] = 19.0;
ptube[2] = 2.5;
- gMC->Gsvolu("QI33","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QI33","TUBE", idtmed[kInox], ptube, 3);
ptube[0] = 0.0;
ptube[1] = 15.0;
ptube[2] = 2.5;
- gMC->Gsvolu("QI43","TUBE", idtmed[kAir], ptube, 3);
- gMC->Gspos("QI43", 1, "QI33", 0.0, 0.0, 0.0, 0, "ONLY");
+ TVirtualMC::GetMC()->Gsvolu("QI43","TUBE", idtmed[kAir], ptube, 3);
+ TVirtualMC::GetMC()->Gspos("QI43", 1, "QI33", 0.0, 0.0, 0.0, 0, "ONLY");
//
// Connecting tube ->
ptube[0] = 0.0;
ptube[1] = 4.5;
ptube[2] = 14.6;
- gMC->Gsvolu("QI34","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QI34","TUBE", idtmed[kInox], ptube, 3);
ptube[0] = 0.0;
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");
+ TVirtualMC::GetMC()->Gsvolu("QI44","TUBE", idtmed[kAir], ptube, 3);
+ TVirtualMC::GetMC()->Gspos("QI44", 1, "QI34", 0.0, 0.0, 0.0, 0, "ONLY");
// <-
//
ptube[0] = 4.6;
ptube[1] = 7.30;
ptube[2] = 2.15;
- gMC->Gsvolu("QI35","TUBE", idtmed[kInox], ptube, 3);
+ TVirtualMC::GetMC()->Gsvolu("QI35","TUBE", idtmed[kInox], ptube, 3);
// <-
- gMC->Gspos("QI32", 1, "QBPM", 0.0, -44.25, zPump, 0, "ONLY");
- gMC->Gspos("QI33", 1, "QBPM", 0.0, -35.00, zPump,idrotm[2002], "ONLY");
- gMC->Gspos("QI34", 1, "QBPM", 0.0, -17.90, zPump,idrotm[2002], "ONLY");
- gMC->Gspos("QI35", 1, "QBPM", 0.0, -24.35, zPump,idrotm[2002], "ONLY");
+ TVirtualMC::GetMC()->Gspos("QI32", 1, "QBPM", 0.0, -44.25, zPump, 0, "ONLY");
+ TVirtualMC::GetMC()->Gspos("QI33", 1, "QBPM", 0.0, -35.00, zPump,idrotm[2002], "ONLY");
+ TVirtualMC::GetMC()->Gspos("QI34", 1, "QBPM", 0.0, -17.90, zPump,idrotm[2002], "ONLY");
+ TVirtualMC::GetMC()->Gspos("QI35", 1, "QBPM", 0.0, -24.35, zPump,idrotm[2002], "ONLY");
- gMC->Gsatt("QBPM", "SEEN", 1);
- gMC->Gsatt("QBEM", "SEEN", 1);
+ gGeoManager->SetVolumeAttribute("QBPM", "SEEN", 1);
+ gGeoManager->SetVolumeAttribute("QBEM", "SEEN", 1);
}
// 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. };
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);
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
AliMixture(16, "VACUUM$ ", aAir, zAir, dAir1, 4, wAir);
//
AliMixture(22, "PA$", aPA, zPA, 1.14, -4, wPA);
//
-
+ // Kapton
+ AliMixture(23, "KAPTON", aKapton, zKapton, dKapton, 4, wKapton);
+ //
// ****************
// Defines tracking media parameters.
//
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);
//
// 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);
}