#include <TLorentzVector.h>
#include <TMath.h>
#include <TVirtualMC.h>
+#include <TString.h>
#include "AliLog.h"
#include "AliMagF.h"
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]);
}
}
-//---------------------------------------------------------------------
-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()
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,
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,
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");
+
+
}
//-------------------------------------------------------------------