3D detector example

{
 gStyle->SetCanvasPreferGL(kTRUE);

  // define a 3D detector with 5 electrodes
  // x=100 , y is 50 and thickness 120
  K3D *det=new K3D(5,100,50,120);
  // define the voltage
  det->Voltage=100;
  // define the drift mesh size and simulation mesh size in microns
  det->SetUpVolume(1,1);
  // define  columns #, postions, weigthing factor 2=0 , material Al=1
  det->SetUpColumn(0,0,0,5,75,2,1);
  det->SetUpColumn(1,100,0,5,75,2,1);
  det->SetUpColumn(2,0,50,5,75,2,1);
  det->SetUpColumn(3,100,50,5,75,2,1);
  det->SetUpColumn(4,50,25,5,-75,16385,1);
 Float_t Pos[3]={100,50,1};
  Float_t Size[3]={100,50,2};
  det->ElRectangle(Pos,Size,0,20);

  det->SetUpElectrodes();
  det->SetBoundaryConditions();


  //define the space charge
  TF3 *f2=new TF3("f2","x[0]*x[1]*x[2]*0+[0]",0,3000,0,3000,0,3000);
  f2->SetParameter(0,-2);
  det->NeffF=f2;

  // calculate weigting field
  // calculate electric field
   det->CalField(0);
   det->CalField(1);

   // set entry points of the track
   det->enp[0]=30;
   det->enp[1]=30;
   det->enp[2]=50;
   det->exp[0]=30;
   det->exp[1]=30;
   det->exp[2]=10;

   // switch on the diffusion
   det->diff=1;
  
   // Show mip track
   TCanvas c1; c1.cd();
   det.ShowMipIR(30);
   // Show electric potential
   TCanvas c2; c2.cd();
   det.Draw("EPxy",99).Draw("COLZ");
   // calcualte induced current
   TCanvas c3; c3.cd();
   det.MipIR(100);
   det->sum.Draw();
   det->neg.Draw("SAME");
   det->pos.Draw("SAME");


}