cePulser.h File Reference

#include "TF1.h"
#include <string>
#include <iostream>

Go to the source code of this file.

Functions
double	cePulserVoltage (int iaIn, const double rfLines[6], double r6, double vfpga)
double	cePulser (int ia, const double rfLines[6], double r6, double qvscale)
double	cePulserVoltage (double *x, double *pars)
double	cePulser (double *x, double *pars)
TF1 *	cePulserVoltageTF1 (std::string fname="cepulserVoltage", double tol=0.01)
int	initCePulserTF1 (TF1 *pf, unsigned int ipoff, double tol=0.01)
TF1 *	cePulserTF1 (std::string fname="cepulser", double tol=0.01)

Function Documentation

double cePulser ( int  ia, const double  rfLines[6], double  r6, double  qvscale  )

inline

Definition at line 61 of file cePulser.h.

                                                                            {
  if ( ia == 0 ) return 0.0;
  if ( ia > 63 ) return 0.0;
  if ( ia < 0 ) return cePulser(-ia, rfLines, r6, qvscale);
  return cePulserVoltage(ia, rfLines, r6, qvscale) -
         cePulserVoltage( 0, rfLines, r6, qvscale);
}

double cePulser ( double *  x, double *  pars  )

inline

Definition at line 84 of file cePulser.h.

                                         {
  int ia = x[0];
  //double* rfLines = pars;
  double r6 = pars[6];
  double qvscale = pars[7];
  return cePulser(ia, pars, r6, qvscale);
}

TF1* cePulserTF1 ( std::string  fname = "cepulser", double  tol = 0.01  )

inline

Definition at line 181 of file cePulser.h.

                                                               {
  using Index = unsigned int;
  const Index npar = 8;
  TF1* pf = new TF1(fname.c_str(), cePulser, -64, 64, npar);
  initCePulserTF1(pf, 0, tol);
  return pf;
}

double cePulserVoltage ( int  iaIn, const double  rfLines[6], double  r6, double  vfpga  )

inline

Definition at line 26 of file cePulser.h.

                                                                                   {
  using Index = unsigned int;
  if ( r6 < 0.0 ) return 0.0;
  if ( iaIn < 0 ) return 0.0;
  if ( iaIn > 63 ) return 0.0;
  Index ia = iaIn;
  double con1 = 0.0;                      // Conductance to V=0
  double con2 = r6 > 0.0 ? 1.0/r6 : 0.0;  // Conductance to V=Vfpga
  for ( Index ilin=0; ilin<6; ++ilin ) {
    double rfLine = rfLines[ilin];
    if ( rfLine <= 0.0 ) continue;    // Means this line is disconnected.
    double conLine = 1.0/rfLine;
    // If this bit is set then the path is to Vfpga.
    if ( (ia >> ilin) & 0x1) con2 += conLine;
    // Otherwise, the path is to ground.
    else con1 += conLine;
  }
  double fac = 1.0;
  if ( con1 > 0.0 ) {
    double rtot = 1.0/con1 + 1.0/con2;
    fac = 1.0/con1/rtot;
  }
  //const double vrail = 1.8;
  return fac*vfpga;
}

double cePulserVoltage ( double *  x, double *  pars  )

inline

Definition at line 73 of file cePulser.h.

                                                {
  int ia = x[0];
  //double* rfLines = pars;
  double r6 = pars[6];
  double qvscale = pars[7];
  return cePulserVoltage(ia, pars, r6, qvscale);
}

TF1* cePulserVoltageTF1 ( std::string  fname = "cepulserVoltage", double  tol = 0.01  )

inline

Definition at line 96 of file cePulser.h.

                                                                             {
  using Index = unsigned int;
  const Index npar = 8;
  TF1* pf = new TF1(fname.c_str(), cePulserVoltage, 0, 63, npar);
  pf->SetParName(0, "R0");
  pf->SetParName(1, "R1");
  pf->SetParName(2, "R2");
  pf->SetParName(3, "R3");
  pf->SetParName(4, "R4");
  pf->SetParName(5, "R5");
  pf->SetParName(6, "R6");
  pf->SetParName(7, "Vfpga");
  pf->SetParameter(0, 32.40);
  pf->SetParameter(1, 16.02);
  pf->SetParameter(2,  8.06);
  pf->SetParameter(3,  4.02);
  pf->SetParameter(4,  2.00);
  pf->SetParameter(5,  1.00);
  pf->SetParameter(6,  1.00);
  pf->SetParameter(7, 1.8);
  if ( tol > 0.0 && tol < 1.0 ) {
    double fmin = 1 - tol;
    double fmax = 1 + tol;
    for ( Index ipar=0; ipar<npar; ++ipar ) {
      double val = pf->GetParameter(ipar);
      pf->SetParLimits(ipar, fmin*val, fmax*val);
    }
  }
  return pf;
}

int initCePulserTF1 ( TF1 *  pf, unsigned int  ipoff, double  tol = 0.01  )

inline

Definition at line 135 of file cePulser.h.

                                                                   {
  using std::cout;
  using std::endl;
  const std::string myname = "initCePulserTF1: ";
  using Index = unsigned int;
  const Index npar = 8;
  if ( pf->GetNpar() < int(npar) ) {
    cout << myname << "Function has fewer than " << npar << " parameters." << endl;
    return 1;
  }
  if ( pf->GetNpar() < int(ipoff + npar) ) {
    cout << myname << "Function has fewer than " << npar << " parameters after offset "
         << ipoff << "." << endl;
    return 2;
  }
  pf->SetParName(ipoff + 0, "R0");
  pf->SetParName(ipoff + 1, "R1");
  pf->SetParName(ipoff + 2, "R2");
  pf->SetParName(ipoff + 3, "R3");
  pf->SetParName(ipoff + 4, "R4");
  pf->SetParName(ipoff + 5, "R5");
  pf->SetParName(ipoff + 6, "R6");
  pf->SetParName(ipoff + 7, "QVscale");
  pf->SetParameter(ipoff + 0, 32.40);
  pf->SetParameter(ipoff + 1, 16.02);
  pf->SetParameter(ipoff + 2,  8.06);
  pf->SetParameter(ipoff + 3,  4.02);
  pf->SetParameter(ipoff + 4,  2.00);
  pf->SetParameter(ipoff + 5,  1.00);
  pf->SetParameter(ipoff + 6,  1.00);
  pf->SetParameter(ipoff + 7, 1.8*183.0);
  if ( tol > 0.0 && tol < 1.0 ) {
    double fmin = 1 - tol;
    double fmax = 1 + tol;
    for ( Index ipar=ipoff; ipar<ipoff+npar; ++ipar ) {
      double val = pf->GetParameter(ipar);
      pf->SetParLimits(ipar, fmin*val, fmax*val);
    }
  }
  return 0;
}