Material of muon filter ist Cast Iron
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 24 Sep 2012 09:45:30 +0000 (09:45 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 24 Sep 2012 09:45:30 +0000 (09:45 +0000)
https://savannah.cern.ch/bugs/?97388
  Bogdan VULPESCU

STRUCT/AliSHIL.cxx
STRUCT/AliSHILv3.cxx

index 350999f..a7f8efa 100644 (file)
@@ -109,7 +109,12 @@ void AliSHIL::CreateMaterials()
   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
   Float_t dAir = 1.20479E-3;
   Float_t dAir1 = 1.20479E-10;
-
+//
+// Cast iron
+//
+  Float_t acasti[4] = {55.847,12.011,28.085,54.938};
+  Float_t zcasti[4] = {26.,6.,14.,25.};
+  Float_t wcasti[4] = {0.929,0.035,0.031,0.005};
   
   Float_t epsil, stmin, tmaxfd, deemax, stemax;
   
@@ -181,6 +186,11 @@ void AliSHIL::CreateMaterials()
   AliMixture(14, "INSULATION1", ains, zins, 0.41, 4, wins);
   AliMixture(34, "INSULATION2", ains, zins, 0.41, 4, wins);
   AliMixture(54, "INSULATION3", ains, zins, 0.41, 4, wins);
+
+  //     Cast iron
+  AliMixture(18, "CAST IRON1", acasti, zcasti, 7.2, 4, wcasti);
+  AliMixture(38, "CAST IRON2", acasti, zcasti, 7.2, 4, wcasti);
+  AliMixture(58, "CAST IRON3", acasti, zcasti, 7.2, 4, wcasti);
   
   // **************** 
   //     Defines tracking media parameters. 
@@ -254,6 +264,11 @@ void AliSHIL::CreateMaterials()
   AliMedium(37, "CC_C1           ", 37, 0, 0,      sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(57, "CC_C2           ", 57, 0, 0,      sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
+  //    Cast iron 
+  AliMedium(18, "CAST_IRON0      ", 18, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(38, "CAST_IRON1      ", 38, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(58, "CAST_IRON2      ", 58, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  
 }
 
 //_____________________________________________________________________________
index 74d7b44..fb0fc6e 100644 (file)
@@ -92,6 +92,9 @@ void AliSHILv3::CreateGeometry()
 //
     TGeoMedium* kMedConcSh  = gGeoManager->GetMedium("SHIL_CC_C2");
 //
+    TGeoMedium* kMedCastiron   = gGeoManager->GetMedium("SHIL_CAST_IRON0");
+    TGeoMedium* kMedCastironSh = gGeoManager->GetMedium("SHIL_CAST_IRON2");
+//
     const Float_t kDegRad = TMath::Pi() / 180.;
     const Float_t kAngle02   = TMath::Tan( 2.00   * kDegRad);   
     const Float_t kAngle0071 = TMath::Tan( 0.71   * kDegRad);   
@@ -1184,7 +1187,7 @@ void AliSHILv3::CreateGeometry()
       TGeoCompositeShape* shMuonFilter = new TGeoCompositeShape("MuonFilter", "FilterO-FilterI:trFilter");
       //
       // !!!!! Needs to be inclined
-      TGeoVolume* voMuonFilter = new TGeoVolume("YMuonFilter", shMuonFilter, kMedSteel);
+      TGeoVolume* voMuonFilter = new TGeoVolume("YMuonFilter", shMuonFilter, kMedCastiron);
 
       // Inner part with higher transport cuts
       Float_t dzMuonFilterH = 50.;
@@ -1193,7 +1196,7 @@ void AliSHILv3::CreateGeometry()
       TGeoTube*   shMuonFilterIH  = new TGeoTube(0., 50., dzMuonFilterH + 5.);
       shMuonFilterIH->SetName("FilterIH");
       TGeoCompositeShape* shMuonFilterH = new TGeoCompositeShape("MuonFilterH", "FilterOH-FilterIH:trFilter");
-      TGeoVolume* voMuonFilterH = new TGeoVolume("YMuonFilterH", shMuonFilterH, kMedSteelSh);
+      TGeoVolume* voMuonFilterH = new TGeoVolume("YMuonFilterH", shMuonFilterH, kMedCastironSh);
       voMuonFilter->AddNode(voMuonFilterH, 1, gGeoIdentity);
       
 //