ROOT logo

ALICE Silicon Drift Detector example (2D)

{
  KDetector det;
  det.SStep=3;
  det.nx=3500;
  det.ny=100;
  det.nz=1;
  Float_t dimX=35000;
  Float_t dimY=300;
  //init geometry
  det.EG=new TH3I("EG","EG",det.nx,0,dimX,det.ny,0,dimY,det.nz,0,1);
  det.EG->GetXaxis()->SetTitle("x [#mum]");
  det.EG->GetYaxis()->SetTitle("y [#mum]");
  det.EG->GetZaxis()->SetTitle("z [#mum]");
  //init material
  det.DM=new TH3I("DM","DM",det.nx,0,dimX,det.ny,0,dimY,det.nz,0,1);
  det.DM->GetXaxis()->SetTitle("x [#mum]");
  det.DM->GetYaxis()->SetTitle("y [#mum]");
  det.DM->GetZaxis()->SetTitle("z [#mum]");
  //init space charge histo
  det.NeffH=new TH3F("Neff","Neff",det.nx,0,dimX,det.ny,0,dimY,det.nz,0,1);
  det.NeffH->GetXaxis()->SetTitle("x [#mum]");
  det.NeffH->GetYaxis()->SetTitle("y [#mum]");
  det.NeffH->GetZaxis()->SetTitle("z [#mum]");

  // Collection electrode
  Float_t BackPos[3]={34930,299,0.5};
  Float_t BackSiz[3]={50,0.5,0.5};
  det.ElRectangle(BackPos,BackSiz,16385,0);  
 
  // Strip definition
  Float_t StripPos[3]={50,300,0};
  Float_t StripSiz[3]={30,0.5,0};

  // definition of votlages to be applied to the field strips

  det->Voltages->Set(291);
  for(int i=0; i<=290; i++)
    {
      StripPos[0]=i*120; 
     det->Voltages[i]=-1800+i*6;
     StripPos[1]=299;
     det.ElRectangle(StripPos,StripSiz,det->SetElecVolt(i),0);  
    
     StripPos[1]=1;
     det.ElRectangle(StripPos,StripSiz,det->SetElecVolt(i),0);  
    }

  //  SetUpMaterial
  for(int j=0;j<=det.ny;j++)
    for(int i=0;i<=det.nx;i++) 
      {
          det->DM->SetBinContent(i,j,1,0); 
          det->NeffH->SetBinContent(i,j,1,0.5);  // very high resistivity Neff=1e11 cm-3
      }

  det->SetBoundaryConditions();

    det->CalField(0);
    det->CalField(1);

    det->SetEntryPoint(27000,300,0.5);
    det->SetExitPoint(27000,1,0.5);

    det.SetDriftHisto(2000e-9);
    det.SetPrecision(1e-9);
    det->diff=1;
    det->ShowMipIR(300);
}
 SiDriftExample_2.C:1
 SiDriftExample_2.C:2
 SiDriftExample_2.C:3
 SiDriftExample_2.C:4
 SiDriftExample_2.C:5
 SiDriftExample_2.C:6
 SiDriftExample_2.C:7
 SiDriftExample_2.C:8
 SiDriftExample_2.C:9
 SiDriftExample_2.C:10
 SiDriftExample_2.C:11
 SiDriftExample_2.C:12
 SiDriftExample_2.C:13
 SiDriftExample_2.C:14
 SiDriftExample_2.C:15
 SiDriftExample_2.C:16
 SiDriftExample_2.C:17
 SiDriftExample_2.C:18
 SiDriftExample_2.C:19
 SiDriftExample_2.C:20
 SiDriftExample_2.C:21
 SiDriftExample_2.C:22
 SiDriftExample_2.C:23
 SiDriftExample_2.C:24
 SiDriftExample_2.C:25
 SiDriftExample_2.C:26
 SiDriftExample_2.C:27
 SiDriftExample_2.C:28
 SiDriftExample_2.C:29
 SiDriftExample_2.C:30
 SiDriftExample_2.C:31
 SiDriftExample_2.C:32
 SiDriftExample_2.C:33
 SiDriftExample_2.C:34
 SiDriftExample_2.C:35
 SiDriftExample_2.C:36
 SiDriftExample_2.C:37
 SiDriftExample_2.C:38
 SiDriftExample_2.C:39
 SiDriftExample_2.C:40
 SiDriftExample_2.C:41
 SiDriftExample_2.C:42
 SiDriftExample_2.C:43
 SiDriftExample_2.C:44
 SiDriftExample_2.C:45
 SiDriftExample_2.C:46
 SiDriftExample_2.C:47
 SiDriftExample_2.C:48
 SiDriftExample_2.C:49
 SiDriftExample_2.C:50
 SiDriftExample_2.C:51
 SiDriftExample_2.C:52
 SiDriftExample_2.C:53
 SiDriftExample_2.C:54
 SiDriftExample_2.C:55
 SiDriftExample_2.C:56
 SiDriftExample_2.C:57
 SiDriftExample_2.C:58
 SiDriftExample_2.C:59
 SiDriftExample_2.C:60
 SiDriftExample_2.C:61
 SiDriftExample_2.C:62
 SiDriftExample_2.C:63
 SiDriftExample_2.C:64
 SiDriftExample_2.C:65
 SiDriftExample_2.C:66
 SiDriftExample_2.C:67
 SiDriftExample_2.C:68
 SiDriftExample_2.C:69