#include <InductionGar2.h>

Public Member Functions
	CTF2 ()

	~CTF2 ()

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 20 of file InductionGar2.h.

Constructor & Destructor Documentation

◆ CTF2()

CTF2::CTF2 ( )

Definition at line 430 of file InductionGar2.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
}

◆ ~CTF2()

CTF2::~CTF2 ( )

Definition at line 445 of file InductionGar2.cxx.

            {
    delete conv;
}

Member Function Documentation

◆ Convolution_Tbranch_fft()

TH1D CTF2::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 455 of file InductionGar2.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::InductionGar2::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 << '\n';
        }
    }
#endif
 
    return h_signal;
}

Referenced by InductionGar2::setMultiElectrons().

Member Data Documentation

◆ conv

ConvolveWithConst* CTF2::conv

Definition at line 22 of file InductionGar2.h.

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

◆ TBranch

std::vector<double> CTF2::TBranch

Definition at line 23 of file InductionGar2.h.

Referenced by CTF2().

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

Public Member Functions

Public Attributes