#include <InductionGar.h>

Public Member Functions
	CEF ()

	~CEF ()

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

Public Attributes
ConvolveWithConst *	conv

std::vector< double >	EBranch

Detailed Description

Definition at line 27 of file InductionGar.h.

Constructor & Destructor Documentation

◆ CEF()

CEF::CEF ( )

Definition at line 415 of file InductionGar.cxx.

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

◆ ~CEF()

CEF::~CEF ( )

Definition at line 430 of file InductionGar.cxx.

         {
    delete conv;
}

Member Function Documentation

◆ Convolution_Ebranch_fft()

TH1D CEF::Convolution_Ebranch_fft ( const double * Input_Curr_plot_001 )

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

Parameters

Input_Curr_plot_001

Returns: TH1D

Definition at line 440 of file InductionGar.cxx.

{
    double xmin(0), xmax(1000);
    const int Npx(2000);
 
    TH1D h_signal("h_signalE", "", 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);
 
#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_Ebranch_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_Ebranch_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* CEF::conv

Definition at line 30 of file InductionGar.h.

Referenced by CEF(), Convolution_Ebranch_fft(), and ~CEF().

◆ EBranch

std::vector<double> CEF::EBranch

Definition at line 31 of file InductionGar.h.

Referenced by CEF().

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

Public Member Functions

Public Attributes