ROOT logo

Use transfer function to get induced current from measurements

{
  ////////////////////////////////////////////////////////////////////////////
  // This is a demonstrator for getting the induced current  waveforms from //    
  // measured waveforms taking into account the transfer function           //
  //////////////////////////////////////////////////////////////////////////// 

  PSTCT stct("NTypeDiode-300um_R.rtct", 11.5,2);
  stct.CorrectBaseLine();   // Baseline correction 
  stct.PrintInfo();         // Information about the read data 
  // group a set of votlages into a MeasureWF structure
  // basically this can be done for any set x,y,z,U1 and U2 parameters
  MeasureWF *wf= stct->Projection(0,3,0,0,0,0,0, stct->NU1); 

  TGraph *gr[10];

  TCTConv *jj=new TCTConv(2000,0.1);      // initialize class
  jj->Load("TestTransferFunc-1500.dat");  // load transfer function

  for(Int_t i=0;i<5;i++)                  //
   {
     hism= wf->GetHA(12+i*4);             // get measured WF
     jj->GetIndCur(hism,0,1500);          // get induced current from WF
     gr[i]=jj->Draw(1);                   // get induced curent 
     gr[i]->SetLineColor(colori[i]);      // set color    
     if(i==0) gr[i]->Draw("AL"); else  gr[i]->Draw("L"); //draw

   }
}
 Get_Im_to_I.C:1
 Get_Im_to_I.C:2
 Get_Im_to_I.C:3
 Get_Im_to_I.C:4
 Get_Im_to_I.C:5
 Get_Im_to_I.C:6
 Get_Im_to_I.C:7
 Get_Im_to_I.C:8
 Get_Im_to_I.C:9
 Get_Im_to_I.C:10
 Get_Im_to_I.C:11
 Get_Im_to_I.C:12
 Get_Im_to_I.C:13
 Get_Im_to_I.C:14
 Get_Im_to_I.C:15
 Get_Im_to_I.C:16
 Get_Im_to_I.C:17
 Get_Im_to_I.C:18
 Get_Im_to_I.C:19
 Get_Im_to_I.C:20
 Get_Im_to_I.C:21
 Get_Im_to_I.C:22
 Get_Im_to_I.C:23
 Get_Im_to_I.C:24
 Get_Im_to_I.C:25
 Get_Im_to_I.C:26
 Get_Im_to_I.C:27
 Get_Im_to_I.C:28