#include <InductionGar2.h>

Public Member Functions
	CEF2 ()

	~CEF2 ()

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

Constructor & Destructor Documentation

◆ CEF2()

CEF2::CEF2 ( )

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

◆ ~CEF2()

CEF2::~CEF2 ( )

Definition at line 498 of file InductionGar2.cxx.

            {
    delete conv;
}

Member Function Documentation

◆ Convolution_Ebranch_fft()

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

Referenced by InductionGar2::setMultiElectrons().

Member Data Documentation

◆ conv

ConvolveWithConst* CEF2::conv

Definition at line 34 of file InductionGar2.h.

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

◆ EBranch

std::vector<double> CEF2::EBranch

Definition at line 35 of file InductionGar2.h.

Referenced by CEF2().

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

Public Member Functions

Public Attributes