def METHOD_NAME(self):
    self.Msun = 4.92549095e-6
    self.sample_rate = 4096
    self.segment_length = 256
    self.low_frequency_cutoff = 30.0
    # chirp params
    self.m1 = 2.0
    self.m2 = 2.5
    self.del_t = 1.0/self.sample_rate
    self.Dl = 40.0
    self.iota = 1.0
    self.phi_c = 2.0
    self.tc_indx = 86*self.sample_rate ## offset from the beginnig of a segment
    self.fmax = 1.0/(6.**1.5 *pi *(self.m1+self.m2)*self.Msun)
    self.zeta = 1.0
    self.thetaS = 0.5
    self.phiS = 2.781
    self.Fp = 0.5*cos(2.0*self.zeta)*(1.0 + cos(self.thetaS)*cos(self.thetaS))*cos(2.0*self.phiS) - \
            sin(2.*self.zeta)*cos(self.thetaS)*sin(2.*self.phiS)
    self.Fc = 0.5*sin(2.0*self.zeta)*(1.0 + cos(self.thetaS)*cos(self.thetaS))*cos(2.0*self.phiS) + \
             cos(2.*self.zeta)*sin(self.thetaS)*sin(2.*self.phiS)
    # params of sin-gaussian
    self.Q = 1.e-1
    self.om = 200.0*pi*2.0
    # use flat psd
    self.seg_len_idx = self.segment_length * self.sample_rate
    self.psd_len = int(self.seg_len_idx/2+1)
    self.Psd = np.ones(self.psd_len)*2.0e-46
    #  generate waveform and chirp signal
    hp, hc = get_td_waveform(approximant="SpinTaylorT5", mass1=self.m1, mass2=self.m2, \
                    delta_t=self.del_t, f_lower=self.low_frequency_cutoff, distance=self.Dl, \
                    inclination=self.iota, coa_phase=self.phi_c)
    # signal which is a noiseless data
    thp = np.zeros(self.seg_len_idx)
    thp[self.tc_indx:len(hp)+self.tc_indx] = hp
    thc = np.zeros(self.seg_len_idx)
    thc[self.tc_indx:len(hc)+self.tc_indx] = hc
    fct = 10.0/15.21377
    self.sig1 = fct*(self.Fp*thp + self.Fc*thc)
    #### template
    h = np.zeros(self.seg_len_idx)
    h[0:len(hp)] = hp
    hpt = TimeSeries(h, self.del_t)
    self.htilde = make_frequency_series(hpt)
    # generate sin-gaussian signal
    time = np.arange(0, len(hp))*self.del_t
    Nby2 = int(len(hp)/2)
    sngt = np.zeros(len(hp))
    for i in range(len(hp)):
        sngt[i] = 9.0e-21*exp(-(time[i]-time[Nby2])**2/self.Q)*sin(self.om*time[i])
    self.sig2 = np.zeros(self.seg_len_idx)
    self.sig2[self.tc_indx:len(sngt)+self.tc_indx] = sngt