#include <InductionGar.h>

Public Member Functions
	CTF ()

	~CTF ()

TH1D	Convolution_Tbranch_fft (const double *Input_Curr_plot_001)
	return convolve of Input_Curr_plot_001 and T_Branch2, a function; tries to mimic Convolution_Tbranch has a option to test correctness.

Public Attributes
ConvolveWithConst *	conv

std::vector< double >	TBranch

Detailed Description

Definition at line 16 of file InductionGar.h.

Constructor & Destructor Documentation

◆ CTF()

CTF::CTF ( )

Definition at line 360 of file InductionGar.cxx.

        {
    double xmin(0), xmax(1000);
    const int Npx(2000);
    TBranch.resize(Npx);// TBranch at x<0 is 0 
    double *p_TBranch=&(TBranch.front());
    for (int i = 0; i < Npx; i++) {
        double x=xmin+(xmax-xmin)*i/Npx;
        TBranch[i]=T_branch2(x);
    }// set initial TBranch values.
    // basically, a N and a (2N-1) needs (3N-2) bins to avoid overlap.
    conv=new ConvolveWithConst( p_TBranch,Npx,Npx);
    //TBranch_fft.resize(Npx*3-2);
 
    //conv.convolve
}

◆ ~CTF()

CTF::~CTF ( )

Definition at line 375 of file InductionGar.cxx.

         {
    delete conv;
}

Member Function Documentation

◆ Convolution_Tbranch_fft()

TH1D CTF::Convolution_Tbranch_fft ( const double * Input_Curr_plot_001 )

return convolve of Input_Curr_plot_001 and T_Branch2, a function; tries to mimic Convolution_Tbranch has a option to test correctness.

Parameters

Input_Curr_plot_001

Returns: TH1D

Definition at line 385 of file InductionGar.cxx.

{
    double xmin(0), xmax(1000);
    const int Npx(2000);
 
    TH1D h_signal("h_signal", "", Npx, xmin, xmax);
    double *p_signal=h_signal.GetArray()+1;//considering underflow.
    this->conv->convolve(p_signal,Input_Curr_plot_001,0,Npx,-0.5);// the '-' before 0.5: from original function.
 
#ifdef INDUCTIONGAR_DEBUG_CONVOLUTION_TBRANCH_FFT_TEST
    const double dx_abs_threhold=1e-6;// think about the accurancy of float; warn if  |x| > threhold
    const double dx_ref_threhold=1e-6;//think about the accurancy of float; warn if |dx/x| > threhold
    TH1F h_signal_ref=Convolution_Tbranch_noskip(Input_Curr_plot_001);
    bool found_mismatch=false;
    for (int i=1; i<=Npx; i++){
        if (abs((h_signal[i]-h_signal_ref[i])/h_signal_ref[i]) > dx_ref_threhold
                and abs(h_signal[i]-h_signal_ref[i]) > dx_abs_threhold ) {
            if (not found_mismatch){
                std::cout<<"CgemDigitizerSvc::InductionGar::Convolution_Tbranch_fft found results not match:  ";
                found_mismatch=true;
            }
            std::cout<<"ref: "<<h_signal_ref[i]<<"; this: "<<h_signal[i]<<"; at i="<<i<<";threhold abs="<<dx_abs_threhold<<"; ref="<<dx_ref_threhold<<std::endl;
        }
    }
#endif
 
    return h_signal;
}

Referenced by InductionGar::setMultiElectrons().

Member Data Documentation

◆ conv

ConvolveWithConst* CTF::conv

Definition at line 18 of file InductionGar.h.

Referenced by Convolution_Tbranch_fft(), CTF(), and ~CTF().

◆ TBranch

std::vector<double> CTF::TBranch

Definition at line 19 of file InductionGar.h.

Referenced by CTF().

The documentation for this class was generated from the following files:

Public Member Functions

Public Attributes