#include "CgemDigitizerSvc/InductionGar.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/DataSvc.h"
#include "CLHEP/Units/PhysicalConstants.h"
#include "G4DigiManager.hh"
#include "Randomize.hh"
#include "G4ios.hh"
#include <iomanip>
#include <iostream>
#include <fstream>
#include <cmath>
#include "TRandom.h"

Typedefs
typedef std::vector< int >	Vint

Functions
double	rectangle2 (double t, double x, double y, int Npx)

double	T_branch2 (double t)

TH1F	Convolution_Tbranch (double *Input_Curr_plot_001)

Typedef Documentation

◆ Vint

typedef std::vector<int> Vint

Definition at line 24 of file InductionGar.cxx.

Function Documentation

◆ Convolution_Tbranch()

TH1F Convolution_Tbranch ( double * Input_Curr_plot_001 )

Definition at line 52 of file InductionGar.cxx.

{
    double xmin(0), xmax(1000);
    const int Npx(2000);
 
    TH1F h_signal("h_signal","",Npx,xmin,xmax);
  double Input_Time_plot_001[Npx];
    for(int i=0; i<Npx; i++) Input_Time_plot_001[i]=h_signal.GetBinCenter(i+1);
 
    double xmin_f1(0.), xmax_f1(1000);
    int    nStep_f1(200);
    double step_f1=(xmax_f1-xmin_f1)/nStep_f1;
    double x_f1(0.0);
 
    double x[Npx];
    double y_001[Npx];
 
    for(int i=0; i<Npx; i++)
    {
        x[i]=xmin+(xmax-xmin)/Npx*(i+0.5);
        double fun_001=0.0;
 
        for(int j=0; j<nStep_f1; j++)
        {
            x_f1=xmin_f1+step_f1*(j+0.5);
            fun_001+=rectangle2(x_f1,Input_Time_plot_001,Input_Curr_plot_001,Npx)*T_branch2(x[i]-x_f1)*step_f1;
 
        }                                                                              
        y_001[i]=fun_001;
    }
 
    for(int init=0; init<Npx; init++){
        h_signal.SetBinContent(init+1,y_001[init]);
    }
 
    return h_signal;
}

Referenced by InductionGar::setMultiElectrons().

◆ rectangle2()

double rectangle2	(	double	t,
		double *	x,
		double *	y,
		int	Npx
	)

Definition at line 29 of file InductionGar.cxx.

{
 
    double I = y[0];
    int id = -1;
    for(int i = 0; i < Npx-1; i++){
        if(t>=x[i]&&t<x[i+1]){id = i; break;}
    }
          if(id != -1){I=y[id]+(y[id+1]-y[id])*(t-x[id])/(x[id+1]-x[id]);} 
  
    return I;
}

Referenced by Convolution_Tbranch().

◆ T_branch2()

double T_branch2 ( double t )

Definition at line 42 of file InductionGar.cxx.

{
    double result = 0.;
        t=t-10.;
    if(t>0) {
        result = 2000.*(0.00181928*exp(-t/3.)-0.0147059*exp(-t/20.)+0.0128866*exp(-t/100.));
    }
    return result;
}