filterresponse()

Tseries filterresponse	(	const int &	size,
		const double &	dt,
		const Options &	opt
	)

Definition at line 42 of file filterresponse.cc.

References Options::integshift, Options::nointeg, Options::tfac, Options::tlim, Options::tshift, and Options::verbose.

Referenced by tdtapertransformation().

{
  if (opt.verbose)
  {
    cout << "  construct 1/sqrt(t); ";
    if (opt.nointeg)
    {
      cout << "tshift: " << opt.tshift << "; "
        << "tlim: " << opt.tlim << "; "
        << "tfac: " << opt.tfac << endl;
    }
    else
    {
      cout << "integshift: " << opt.integshift << " means " 
        << opt.integshift*dt << "s time shift" << endl;
    }
  }
  Tseries retval=Tseries(0,size-1);
  aff::Iterator<Tseries> I(retval);
  for (unsigned int i=0; i<retval.size(); ++i)
  {
    if (opt.nointeg)
    {
      double t=dt*(double(i)+opt.tshift);
      if (t < (opt.tlim*dt)) { t=opt.tfac*dt; }
      *I=sqrt(1./t);
    }
    else
    {
      if (i == 0)
      {
        double tref=1.-opt.integshift;
        *I=2.*sqrt(tref/dt);
      }
      else
      {
        double tref=(double(i)-opt.integshift);
        *I=2.*(sqrt(tref+1)-sqrt(tref))/sqrt(dt);
      }
    }
    ++I;
  }
  return(retval);
}

Here is the caller graph for this function: