Masking of not calibrated chambers
[u/mrichter/AliRoot.git] / T0 / AliT0v1.cxx
index 2bef8bf..fe10a62 100644 (file)
@@ -41,6 +41,7 @@
 #include <TLorentzVector.h>
 #include <TMath.h>
 #include <TVirtualMC.h>
+#include <TString.h>
 
 #include "AliLog.h"
 #include "AliMagF.h"
@@ -261,12 +262,14 @@ void AliT0v1::CreateGeometry()
     new TGeoTube("0HOLE", 0, 1.51, 6.5);
     TGeoTranslation *tr [12];
     Double_t angle  = 2 * TMath::Pi() / 12;
-    Char_t nameTr[40];
+    // Char_t nameTr[40];
+    TString nameTr;
     for (Int_t itr=0; itr<12; itr++) {
-      sprintf (nameTr,"0TR%i",itr+1);
+      //      sprintf (nameTr,"0TR%i",itr+1);
+      nameTr = Form("0TR%i",itr+1);
       x = 6.5 * TMath::Sin(itr * angle);
       y = 6.5 * TMath::Cos(itr * angle);
-      tr[itr] = new TGeoTranslation(nameTr,x,y,6.5);
+      tr[itr] = new TGeoTranslation(nameTr.Data(),x,y,6.5);
       tr[itr]->RegisterYourself();
       stlin->AddNode(ins,itr+13,tr[itr]);
     }
@@ -477,7 +480,8 @@ void AliT0v1::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-11;
+                                                                               
 // Scintillator CH
    Float_t ascin[2]={1.01,12.01};
    Float_t zscin[2]={1,6};
@@ -519,7 +523,7 @@ void AliT0v1::CreateMaterials()
     
 //*** Definition Of avaible T0 materials ***
    AliMaterial(0, "Steel$", 55.850,26.,7.87,1.76,999);
-   AliMaterial(1, "Vacuum$", 1.e-16,1.e-16,1.e-16,1.e16,999);
+   AliMixture(1, "Vacuum$", aAir, zAir, dAir1,4,wAir);
    AliMixture(2, "Air$", aAir, zAir, dAir,4,wAir);
    AliMaterial(10, "CarbonPlastic$", 12.01, 6.0, 2.26, 18.8,999); 
    AliMaterial(11, "Aliminium$", 26.98, 13.0, 2.7, 8.9,999); 
@@ -569,32 +573,6 @@ void AliT0v1::CreateMaterials()
    
    
 }
-//---------------------------------------------------------------------
-void AliT0v1::DrawDetector()
-{
-//
-// Draw a shaded view of the Forward multiplicity detector version 0
-//
-  
-  //Set ALIC mother transparent
-  gMC->Gsatt("ALIC","SEEN",0);
-  //
-  //Set volumes visible
-  //  gMC->Gsatt("0STR","SEEN",0);
-  //  gMC->Gsatt("0INS","SEEN",0);
-  // gMC->Gsatt("0PMT","SEEN",1);
-  // gMC->Gsatt("0DIV","SEEN",1);
-  //
-  gMC->Gdopt("hide","off");
-  gMC->Gdopt("shad","on");
-  gMC->SetClipBox(".");
-  gMC->SetClipBox("*",0,1000,-1000,1000,-1000,1000);
-  gMC->DefaultRange();
-  gMC->Gdraw("alic",90,0,0,-35,9.5,.6,0.6);
-  //gMC->Gdraw("alic",0,0,0,10,9.5,.8,0.8); //other side view
-  gMC->Gdhead(1111,"T-Zero detector");
-  gMC->Gdopt("hide","off");
-}
 
 //-------------------------------------------------------------------
 void AliT0v1::DefineOpticalProperties()
@@ -606,7 +584,8 @@ void AliT0v1::DefineOpticalProperties()
    const Int_t kNbins=31;
    
    Float_t rindexSiO2[kNbins],  efficAll[kNbins], rindexAir[kNbins], absorAir[kNbins],rindexCathodeNext[kNbins], absorbCathodeNext[kNbins];
-            
+   Double_t efficMet[kNbins], aReflMet[kNbins];
+           
    // quartz 20mm
    Float_t aAbsSiO2[kNbins]={29.0, 28.6, 28.3, 27.7, 27.3, 26.7, 26.4, 
                             25.9, 25.3, 24.9, 24.5, 23.7, 
@@ -621,6 +600,13 @@ void AliT0v1::DefineOpticalProperties()
                             5.8,  5.97, 6.16, 6.36, 6.57, 
                             6.8,  7.04, 7.3,  7.58, 7.89, 
                             8.22, 8.57, 8.97, 9.39 };  
+  Double_t dPckov[kNbins]  ={3.87, 3.94, 4.02, 4.11, 4.19, 4.29, 4.38,
+                             4.48, 4.58, 4.69, 4.81, 4.93,
+                             5.05, 5.19, 5.33, 5.48, 5.63,
+                             5.8,  5.97, 6.16, 6.36, 6.57,
+                             6.8,  7.04, 7.3,  7.58, 7.89,
+                             8.22, 8.57, 8.97, 9.39 };
+
    /*     
    Float_t effCathode[kNbins]={0.11, 0.13, 0.15, 0.16, 0.18, 0.19, 0.20,
                              0.21, 0.22, 0.23, 0.24, 0.26, 
@@ -630,26 +616,37 @@ void AliT0v1::DefineOpticalProperties()
                              0.17, 0.17, 0.2, 0.23};
    */     
    //  Float_t aAbsSiO2[kNbins]; //quartz 30mm
-  for(i=0;i<kNbins;i++)
+ for(i=0;i<kNbins;i++)
+
     {
-      aPckov[i]=aPckov[i]*1e-9;//Photons energy bins 4 eV - 8.5 eV step 0.1 eV   
-      rindexAir[i]=0;
+      aPckov[i]=aPckov[i]*1e-9;//Photons energy bins 4 eV - 8.5 eV step 0.1 eV
+      dPckov[i]=dPckov[i]*1e-9;//Photons energy bins 4 eV - 8.5 eV step 0.1 eV 
+      //      rindexAir[i]=0.0001;
+      rindexAir[i] = 1.;
       rindexSiO2[i]=1.458; //refractive index for qwarts
       rindexCathodeNext[i]=0;
       efficAll[i]=1.;
-
+      efficMet[i]=0.;
+      aReflMet[i]=1.;
       //      aAbsSiO2[i]=28.5; //quartz 30mm
       absorAir[i]=0.3;      
       absorbCathodeNext[i]=0;
 
     }
   
-  
-   gMC->SetCerenkov (idtmed[kOpGlass], kNbins, aPckov, aAbsSiO2, efficAll, rindexSiO2 );
+  gMC->SetCerenkov (idtmed[kOpGlass], kNbins, aPckov, aAbsSiO2, efficAll, rindexSiO2 );
    // gMC->SetCerenkov (idtmed[kOpGlassCathode], kNbins, aPckov, aAbsSiO2, effCathode, rindexSiO2 );
    gMC->SetCerenkov (idtmed[kOpGlassCathode], kNbins, aPckov, aAbsSiO2,efficAll , rindexSiO2 );
   gMC->SetCerenkov (idtmed[kOpAir], kNbins, aPckov,absorAir , efficAll,rindexAir );
    gMC->SetCerenkov (idtmed[kOpAirNext], kNbins, aPckov,absorbCathodeNext , efficAll, rindexCathodeNext);
+
+   //Define a boarder for radiator optical properties
+   gMC->DefineOpSurface("surfRd", kUnified /*kGlisur*/,kDielectric_metal,kPolished, 0.);
+   gMC->SetMaterialProperty("surfRd", "EFFICIENCY", kNbins, dPckov, efficMet);
+   gMC->SetMaterialProperty("surfRd", "REFLECTIVITY", kNbins, dPckov, aReflMet);
+   gMC->SetBorderSurface("0TOPborder", "0TOP",1,"0TOO",1, "surfRd");
+
+
 }
 
 //-------------------------------------------------------------------