TMatacq Class Reference

#include <TMatacq.h>

Inheritance diagram for TMatacq:

Public Member Functions
int	compute_trise ()
int	countBadPulses (Int_t)
int	doFit ()
void	enterdata (Int_t)
int	findPeak ()
double	getAmpl ()
double	getBaseLine ()
double	getFwhm ()
int	getPt1 ()
int	getPt2 ()
int	getPtm ()
double	getsigBaseLine ()
double	getsigTimpeak ()
double	getSlide ()
double	getTimax ()
double	getTimpeak ()
double	getTrise ()
double	getWidth20 ()
double	getWidth80 ()
void	init ()
void	printitermatacqData (Int_t, Int_t, Int_t)
void	printmatacqData (Int_t, Int_t, Int_t)
int	rawPulseAnalysis (Int_t, Double_t *)
	TMatacq (int, int, int, int, int, int, int, int, int, int, int)
	~TMatacq () override

Static Public Attributes
static constexpr int	NMAXSAMP = 100
static constexpr int	NSPARAB = 16

Private Member Functions
double	interpolate (double)

Private Attributes
double	ampl
int	bing [101]
double	bl
double	bong [NMAXSAMP]
double	comp_peak [1200]
double	comp_trise [1200]
int	endsample
int	firstsample
int	fNsamples
int	fNum_samp_aft_max
int	fNum_samp_bef_max
double	fv1 [NSPARAB]
double	fv2 [NSPARAB]
double	fv3 [NSPARAB]
int	lastsample
double	level1
double	level2
double	level3
double	meantrise
int	nevlasers
int	nevmtq0
int	nevmtq1
double	nsigcut
int	nslide
double	pkval
int	presample
int	samplemax
double	sigbl
double	sigpkval
double	sigtrise
double	slidingmean
int	status [1200]
double	t [NSPARAB]
double	timeatmax
double	trise
double	val [NSPARAB]
double	val_max
double	width20
double	width50
double	width80

Detailed Description

Definition at line 6 of file TMatacq.h.

Constructor & Destructor Documentation

TMatacq()

TMatacq::TMatacq	(	int	Ntot,
		int	Nsamp1,
		int	Nsamp2,
		int	cut,
		int	Nbef,
		int	Naft,
		int	niv1,
		int	niv2,
		int	niv3,
		int	nevl,
		int	NSlide
	)

Definition at line 31 of file TMatacq.cc.

References DMR_cfg::cut, dqmdumpme::k, getInfo::level1, and mps_update::status.

                                                                                                                       {
  fNsamples = Ntot;
  presample = Nsamp1;
  endsample = Nsamp2;
  nsigcut = (double)cut;
  fNum_samp_bef_max = Nbef;
  fNum_samp_aft_max = Naft;
  level1 = ((double)niv1) / 100.;
  level2 = ((double)niv2) / 100.;
  level3 = ((double)niv3) / 100.;
  nevlasers = nevl;
  slidingmean = 0.0;
  nslide = NSlide;
  for (int k = 0; k < nevlasers; k++)
    status[k] = 0;
  for (int k = 0; k < nevlasers; k++)
    status[k + nevlasers] = 0;

  nevmtq0 = 0;
  nevmtq1 = 0;
}

~TMatacq()

TMatacq::~TMatacq ( )

override

Definition at line 55 of file TMatacq.cc.

{}

Member Function Documentation

compute_trise()

int TMatacq::compute_trise

(

)

Definition at line 347 of file TMatacq.cc.

References relativeConstraints::error, createfilelist::int, and HiCaloJetParameters_cff::interpolate.

Referenced by EcalMatacqAnalyzer::endJob().

                           {
  int error;
  trise = 0.;

  double t20 = interpolate(ampl * level2);
  if (t20 < 0.) {
    error = (int)-t20;
    return error;
  }
  double t80 = interpolate(ampl * level3);
  if (t80 < 0.) {
    error = (int)-t80;
    return error;
  }
  trise = t80 - t20;

  return 0;
}

countBadPulses()

int TMatacq::countBadPulses

(

Int_t

)

Definition at line 484 of file TMatacq.cc.

References corrVsCorr::filename, mps_fire::i, and mps_update::status.

                                          {
  FILE *fmatacq;
  char filename[80];
  sprintf(filename, "badevtsMatacq%d.dat", gRunNumber);
  fmatacq = fopen(filename, "w");
  if (fmatacq == nullptr)
    printf("Error while opening file : %s\n", filename);

  int nevbad = 0;
  for (Int_t i = 0; i < nevmtq0 + nevmtq1; i++) {
    if (comp_trise[i] < meantrise - 3. * sigtrise) {
      nevbad++;
      fprintf(fmatacq, "%d \n", status[i]);
      status[i] = -1;
    }
    if (comp_trise[i] > meantrise + 3. * sigtrise) {
      nevbad++;
      fprintf(fmatacq, "%d \n", status[i]);
      status[i] = -1;
    }
  }

  int iret = fclose(fmatacq);
  printf(" Closing file : %d\n", iret);
  return nevbad;
}

doFit()

int TMatacq::doFit

(

)

Definition at line 213 of file TMatacq.cc.

References a, b, c, DeadROC_duringRun::f1, DeadROC_duringRun::f2, mps_fire::i, dqmiolumiharvest::j, heppy_batch::val, and ApeEstimator_cff::width.

Referenced by EcalMatacqAnalyzer::endJob().

                   {
  int testneg = 0;
  ampl = 0.;
  timeatmax = 0.;
  int heresamplemax = samplemax - firstsample;

  int beginfit = heresamplemax - fNum_samp_bef_max;
  int endfit = heresamplemax + fNum_samp_aft_max;

  int nval = endfit - beginfit + 1;
  if (nval != NSPARAB)
    return 201;
  for (int kn = beginfit; kn <= endfit; kn++) {
    if (bong[kn] <= 0)
      testneg = 1;
  }
  if (testneg == 1)
    return 202;

  for (int i = 0; i < nval; i++) {
    val[i] = bong[beginfit + i];
    fv1[i] = 1.;
    fv2[i] = (double)(i);
    fv3[i] = ((double)(i)) * ((double)(i));
  }

  TVectorD y(nval, val);
  TVectorD f1(nval, fv1);
  TVectorD f2(nval, fv2);
  TVectorD f3(nval, fv3);

  double bj[3];
  bj[0] = f1 * y;
  bj[1] = f2 * y;
  bj[2] = f3 * y;
  TVectorD b(3, bj);

  double aij[9];
  aij[0] = f1 * f1;
  aij[1] = f1 * f2;
  aij[2] = f1 * f3;
  aij[3] = f2 * f1;
  aij[4] = f2 * f2;
  aij[5] = f2 * f3;
  aij[6] = f3 * f1;
  aij[7] = f3 * f2;
  aij[8] = f3 * f3;
  TMatrixD a(3, 3, aij);

  double det1;
  a.InvertFast(&det1);

  TVectorD c = a * b;

  double *par = c.GetMatrixArray();
  if (par[2] != 0.) {
    timeatmax = -par[1] / (2. * par[2]);
    ampl = par[0] - par[2] * (timeatmax * timeatmax);
  }

  if (ampl <= 0.) {
    ampl = bong[heresamplemax];
    return 203;
  }

  if ((int)timeatmax > NSPARAB)
    return 204;
  if ((int)timeatmax < 0)
    return 205;

  timeatmax += (beginfit + firstsample);

  int tplus[3], tminus[3];
  double xampl[3];
  xampl[0] = 0.2 * ampl;
  xampl[1] = 0.5 * ampl;
  xampl[2] = 0.8 * ampl;

  int hitplus[3];
  int hitminus[3];
  double width[3];

  for (int i = 0; i < 3; i++) {
    hitplus[i] = 0;
    hitminus[i] = 0;
    width[i] = 0.0;
    tplus[i] = firstsample + lastsample;
    tminus[i] = 0;
  }

  // calculate first estimate of half width
  int im = heresamplemax;
  int iplusbound = firstsample + lastsample - im;

  for (int j = 0; j < 3; j++) {
    for (int i = 1; i < im; i++) {
      if (bong[im - i] < xampl[j] && bong[im - i + 1] > xampl[j]) {
        tminus[j] = im - i;
        hitminus[j]++;
        i = im;
      }
    }

    for (int i = 0; i < iplusbound; i++) {
      if (bong[im + i] > xampl[j] && bong[im + i + 1] < xampl[j]) {
        tplus[j] = im + i;
        hitplus[j]++;
        i = iplusbound;
      }
    }

    // interpolate to get a better estimate

    double slopeplus = double(bong[tplus[j] + 1] - bong[tplus[j]]);
    double slopeminus = double(bong[tminus[j] + 1] - bong[tminus[j]]);

    double timeminus = double(tminus[j]) + 0.5;
    if (slopeminus != 0)
      timeminus = tminus[j] + (xampl[j] - double(bong[tminus[j]])) / slopeminus;

    double timeplus = double(tplus[j]) + 0.5;
    if (slopeplus != 0)
      timeplus = tplus[j] + (xampl[j] - double(bong[tplus[j]])) / slopeplus;

    width[j] = timeplus - timeminus;
  }

  width20 = width[0];
  width50 = width[1];
  width80 = width[2];

  return 0;
}

enterdata()

void TMatacq::enterdata

(

Int_t

anevt

)

Definition at line 397 of file TMatacq.cc.

References mps_update::status.

                                   {
  if (anevt < 2 * nevlasers) {
    if (anevt < nevlasers) {
      nevmtq0++;
      status[nevmtq0 - 1] = anevt;
      comp_trise[nevmtq0 - 1] = trise;
      comp_peak[nevmtq0 - 1] = pkval;
    } else {
      nevmtq1++;
      status[nevmtq0 + nevmtq1 - 1] = anevt;
      comp_trise[nevmtq0 + nevmtq1 - 1] = trise;
      comp_peak[nevmtq0 + nevmtq1 - 1] = pkval;
    }
  } else {
    if (anevt < 3 * nevlasers) {
      nevmtq0++;
      status[nevmtq0 - 1] = anevt;
      comp_trise[nevmtq0 - 1] = trise;
      comp_peak[nevmtq0 - 1] = pkval;
    } else {
      nevmtq1++;
      status[nevmtq0 + nevmtq1 - 1] = anevt;
      comp_trise[nevmtq0 + nevmtq1 - 1] = trise;
      comp_peak[nevmtq0 + nevmtq1 - 1] = pkval;
    }
  }
}

findPeak()

int TMatacq::findPeak

(

)

Definition at line 136 of file TMatacq.cc.

References TMarkov::getPeakValue(), createfilelist::int, dqmdumpme::k, getInfo::level1, and TMarkov::peakFinder().

Referenced by EcalMatacqAnalyzer::endJob().

                      {
  int k;
  int nbinf = 0;
  int jfind = 0;
  int nbsup = 0;
  double thres = val_max * level1;

  for (k = 0, jfind = 0; k < NMAXSAMP; k++) {
    if (jfind == 0) {
      if (bong[k] > thres) {
        nbinf = k;
        jfind = 1;
      }
    }
  }
  if (jfind == 0)
    nbinf = 0;

  for (k = NMAXSAMP, jfind = 0; k > nbinf; k--) {
    if (jfind == 0) {
      if (bong[k] > thres) {
        nbsup = k;
        jfind = 1;
      }
    }
  }
  if (nbsup == 0)
    nbsup = nbinf;

  pkval = 0.;
  sigpkval = 0.5;
  if (nbsup == nbinf) {
    return 301;
  } else {
    if (nbinf > 4)
      nbinf -= 3;
    else
      nbinf = 1;
    if (nbsup < NMAXSAMP - 4)
      nbsup += 3;
    else
      nbsup = NMAXSAMP;

    for (k = 0; k < nbinf; k++)
      bong[k] = 0.;
    for (k = nbsup; k < NMAXSAMP; k++)
      bong[k] = 0.;

    for (k = 0, jfind = 0; k < NMAXSAMP; k++) {
      if (bong[k] > 0.)
        jfind++;
    }
    if (jfind == 0) {
      return 302;
    } else if (jfind == 1) {
      return 303;
    } else {
      for (k = 0; k < NMAXSAMP; k++) {
        if (k < 100)
          bing[k + 1] = (int)bong[k];
      }

      TMarkov *peak = new TMarkov();

      peak->peakFinder(&bing[0]);
      pkval = peak->getPeakValue(0);
      sigpkval = peak->getPeakValue(1);

      delete peak;

      pkval += (firstsample - 1);
    }
  }

  return 0;
}

getAmpl()

double TMatacq::getAmpl ( )

inline

Definition at line 68 of file TMatacq.h.

References ampl.

Referenced by EcalMatacqAnalyzer::endJob().

{ return ampl; }

getBaseLine()

double TMatacq::getBaseLine ( )

inline

Definition at line 62 of file TMatacq.h.

References bl.

Referenced by EcalMatacqAnalyzer::endJob().

{ return bl; }

getFwhm()

double TMatacq::getFwhm ( )

inline

Definition at line 72 of file TMatacq.h.

References width50.

Referenced by EcalMatacqAnalyzer::endJob().

{ return width50; }

getPt1()

int TMatacq::getPt1 ( )

inline

Definition at line 58 of file TMatacq.h.

References firstsample.

{ return firstsample; }

getPt2()

int TMatacq::getPt2 ( )

inline

Definition at line 59 of file TMatacq.h.

References lastsample.

{ return lastsample; }

getPtm()

int TMatacq::getPtm ( )

inline

Definition at line 60 of file TMatacq.h.

References samplemax.

{ return samplemax; }

getsigBaseLine()

double TMatacq::getsigBaseLine ( )

inline

Definition at line 63 of file TMatacq.h.

References sigbl.

Referenced by EcalMatacqAnalyzer::endJob().

{ return sigbl; }

getsigTimpeak()

double TMatacq::getsigTimpeak ( )

inline

Definition at line 66 of file TMatacq.h.

References sigpkval.

Referenced by EcalMatacqAnalyzer::endJob().

{ return sigpkval; }

getSlide()

double TMatacq::getSlide ( )

inline

Definition at line 75 of file TMatacq.h.

References slidingmean.

Referenced by EcalMatacqAnalyzer::endJob().

{ return slidingmean; }

getTimax()

double TMatacq::getTimax ( )

inline

Definition at line 69 of file TMatacq.h.

References timeatmax.

Referenced by EcalMatacqAnalyzer::endJob().

{ return timeatmax; }

getTimpeak()

double TMatacq::getTimpeak ( )

inline

Definition at line 65 of file TMatacq.h.

References pkval.

Referenced by EcalMatacqAnalyzer::endJob().

{ return pkval; }

getTrise()

double TMatacq::getTrise ( )

inline

Definition at line 71 of file TMatacq.h.

References trise.

Referenced by EcalMatacqAnalyzer::endJob().

{ return trise; }

getWidth20()

double TMatacq::getWidth20 ( )

inline

Definition at line 73 of file TMatacq.h.

References width20.

Referenced by EcalMatacqAnalyzer::endJob().

{ return width20; }

getWidth80()

double TMatacq::getWidth80 ( )

inline

Definition at line 74 of file TMatacq.h.

References width80.

Referenced by EcalMatacqAnalyzer::endJob().

{ return width80; }

init()

void TMatacq::init

(

)

Definition at line 17 of file TMatacq.cc.

References dqmdumpme::k, and submitPVValidationJobs::t.

                   {
  for (int k = 0; k < NMAXSAMP; k++)
    bong[k] = 0.;

  for (int k = 0; k <= 100; k++)
    bing[k] = 0;

  for (int k = 0; k < NSPARAB; k++)
    t[k] = (double)k;

  return;
}

interpolate()

double TMatacq::interpolate ( double  amplx )

private

Definition at line 366 of file TMatacq.cc.

References createfilelist::int, dqmdumpme::k, and slope.

                                        {
  double T;
  int kmax = (int)pkval - firstsample;

  int bin_low = 0;
  for (Int_t k = 0; k < kmax; k++)
    if (0. < bong[k] && bong[k] < amplx) {
      bin_low = k;
    }
  if (bin_low == 0)
    return -301.;
  int bin_high = 0;
  for (Int_t k = kmax; k >= 0; k--)
    if (bong[k] > amplx) {
      bin_high = k;
    }
  if (bin_high == 0)
    return -302.;
  if (bin_high < bin_low)
    return -303.;

  if (bin_low == bin_high) {
    T = (double)bin_high;
  } else {
    double slope = (bong[bin_high] - bong[bin_low]) / ((double)(bin_high - bin_low));
    T = (double)bin_high - (bong[bin_high] - amplx) / slope;
  }

  return T;
}

printitermatacqData()

void TMatacq::printitermatacqData	(	Int_t	,
		Int_t	,
		Int_t
	)

Definition at line 511 of file TMatacq.cc.

References corrVsCorr::filename, mps_fire::i, mathSSE::sqrt(), contentValuesCheck::ss, and mps_update::status.

                                                                          {
  FILE *fmatacq;
  char filename[80];
  int i;
  double ss;
  sprintf(filename, "runiterMatacq%d.pedestal", gRunNumber);
  fmatacq = fopen(filename, "w");
  if (fmatacq == nullptr)
    printf("Error while opening file : %s\n", filename);

  int nevmtqgood = 0;
  double sumtrise = 0.;
  double sumtrise2 = 0.;
  int timestop = timestart + 3;
  double mintrise = 10000.;
  double maxtrise = 0.;
  for (i = 0; i < nevmtq0; i++) {
    if (status[i] >= 0) {
      nevmtqgood++;
      if (comp_trise[i] > maxtrise) {
        maxtrise = comp_trise[i];
      }
      if (comp_trise[i] < mintrise) {
        mintrise = comp_trise[i];
      }
      sumtrise += comp_trise[i];
      sumtrise2 += comp_trise[i] * comp_trise[i];
    }
  }
  meantrise = sumtrise / ((double)nevmtqgood);
  ss = (sumtrise2 / ((double)nevmtqgood) - meantrise * meantrise);
  if (ss < 0.)
    ss = 0.;
  sigtrise = sqrt(ss);
  fprintf(fmatacq,
          "%d %d %d %d %f %f %f %f\n",
          nevmtqgood,
          color,
          timestart,
          timestop,
          meantrise,
          sigtrise,
          mintrise,
          maxtrise);

  nevmtqgood = 0;
  sumtrise = 0.;
  sumtrise2 = 0.;
  mintrise = 10000.;
  maxtrise = 0.;
  for (i = nevmtq0; i < nevmtq0 + nevmtq1; i++) {
    if (status[i] >= 0) {
      nevmtqgood++;
      if (comp_trise[i] > maxtrise) {
        maxtrise = comp_trise[i];
      }
      if (comp_trise[i] < mintrise) {
        mintrise = comp_trise[i];
      }
      sumtrise += comp_trise[i];
      sumtrise2 += comp_trise[i] * comp_trise[i];
    }
  }
  meantrise = sumtrise / ((double)nevmtqgood);
  ss = (sumtrise2 / ((double)nevmtqgood) - meantrise * meantrise);
  if (ss < 0.)
    ss = 0.;
  sigtrise = sqrt(ss);
  fprintf(fmatacq,
          "%d %d %d %d %f %f %f %f\n",
          nevmtqgood,
          color,
          timestart,
          timestop,
          meantrise,
          sigtrise,
          mintrise,
          maxtrise);

  int iret = fclose(fmatacq);
  printf(" Closing file : %d\n", iret);
}

printmatacqData()

void TMatacq::printmatacqData	(	Int_t	,
		Int_t	,
		Int_t
	)

Definition at line 425 of file TMatacq.cc.

References corrVsCorr::filename, mps_fire::i, mathSSE::sqrt(), and contentValuesCheck::ss.

                                                                      {
  FILE *fmatacq;
  char filename[80];
  int i;
  double ss;
  sprintf(filename, "runMatacq%d.pedestal", gRunNumber);
  fmatacq = fopen(filename, "w");
  if (fmatacq == nullptr)
    printf("Error while opening file : %s\n", filename);

  double sumtrise = 0.;
  double sumtrise2 = 0.;
  int timestop = timestart + 3;
  double mintrise = 10000.;
  double maxtrise = 0.;
  for (i = 0; i < nevmtq0; i++) {
    if (comp_trise[i] > maxtrise) {
      maxtrise = comp_trise[i];
    }
    if (comp_trise[i] < mintrise) {
      mintrise = comp_trise[i];
    }
    sumtrise += comp_trise[i];
    sumtrise2 += comp_trise[i] * comp_trise[i];
  }
  meantrise = sumtrise / ((double)nevmtq0);
  ss = (sumtrise2 / ((double)nevmtq0) - meantrise * meantrise);
  if (ss < 0.)
    ss = 0.;
  sigtrise = sqrt(ss);
  fprintf(
      fmatacq, "%d %d %d %d %f %f %f %f\n", nevmtq0, color, timestart, timestop, meantrise, sigtrise, mintrise, maxtrise);

  sumtrise = 0.;
  sumtrise2 = 0.;
  mintrise = 10000.;
  maxtrise = 0.;
  for (i = nevmtq0; i < nevmtq0 + nevmtq1; i++) {
    if (comp_trise[i] > maxtrise) {
      maxtrise = comp_trise[i];
    }
    if (comp_trise[i] < mintrise) {
      mintrise = comp_trise[i];
    }
    sumtrise += comp_trise[i];
    sumtrise2 += comp_trise[i] * comp_trise[i];
  }
  meantrise = sumtrise / ((double)nevmtq1);
  ss = (sumtrise2 / ((double)nevmtq1) - meantrise * meantrise);
  if (ss < 0.)
    ss = 0.;
  sigtrise = sqrt(ss);
  fprintf(
      fmatacq, "%d %d %d %d %f %f %f %f\n", nevmtq1, color, timestart, timestop, meantrise, sigtrise, mintrise, maxtrise);

  int iret = fclose(fmatacq);
  printf(" Closing file : %d\n", iret);
}

rawPulseAnalysis()

int TMatacq::rawPulseAnalysis	(	Int_t	Nsamp,
		Double_t *	adc
	)

Definition at line 57 of file TMatacq.cc.

References gpuClustering::adc, mps_fire::i, init, dqmdumpme::k, mathSSE::sqrt(), and contentValuesCheck::ss.

Referenced by EcalMatacqAnalyzer::endJob().

{
  using namespace std;

  //  std::cout << "Entering rawPulseAnalysis" << std::endl;

  int k, ithr;
  double dsum = 0., dsum2 = 0.;

  //  std::cout << "calling init" << std::endl;
  init();
  //  std::cout << ".......done" << std::endl;

  if (Nsamp != fNsamples) {
    printf("found different number of samples fNsamples=%d Nsamp=%d\n", fNsamples, Nsamp);
    return 100;
  }

  for (k = 0; k < presample; k++) {
    dsum += adc[k];
    dsum2 += adc[k] * adc[k];
  }
  bl = dsum / ((double)presample);
  double ss = (dsum2 / ((double)presample) - bl * bl);
  if (ss < 0.)
    ss = 0.;
  sigbl = sqrt(ss);
  for (ithr = 0, k = presample; k < endsample; k++) {
    if (adc[k] > (bl + nsigcut * sigbl) && ithr == 0) {
      ithr = 1;
      firstsample = k;
    }
  }

  if (ithr == 0)
    return 101;

  for (ithr = 0, k = firstsample; k < Nsamp; k++) {
    if (adc[k] < (bl + nsigcut * sigbl) && ithr == 0) {
      ithr = 1;
      lastsample = k;
    }
  }
  if (ithr == 0)
    lastsample = Nsamp;

  if (lastsample > firstsample + NMAXSAMP)
    lastsample = firstsample + NMAXSAMP;

  val_max = 0.;
  samplemax = 0;
  for (Int_t is = firstsample; is < lastsample; is++) {
    bong[is - firstsample] = adc[is] - bl;
    if (bong[is - firstsample] > val_max) {
      val_max = bong[is - firstsample];
      samplemax = is;
    }
  }
  if (samplemax == 0)
    return 103;
  if (samplemax > lastsample)
    return 104;
  if (samplemax < firstsample)
    return 105;

  int endslide = samplemax - nslide;
  int beginslide = nslide;
  int islidingmean = 0;
  slidingmean = 0.0;

  for (int i = beginslide; i < endslide; i++) {
    slidingmean += adc[i];
    islidingmean += 1;
  }
  if (islidingmean != 0)
    slidingmean /= double(islidingmean);

  return 0;
}

Member Data Documentation

ampl

double TMatacq::ampl

private

Definition at line 24 of file TMatacq.h.

Referenced by getAmpl().

bing

int TMatacq::bing[101]

private

Definition at line 17 of file TMatacq.h.

bl

double TMatacq::bl

private

Definition at line 23 of file TMatacq.h.

Referenced by getBaseLine().

bong

double TMatacq::bong[NMAXSAMP]

private

Definition at line 20 of file TMatacq.h.

comp_peak

double TMatacq::comp_peak[1200]

private

Definition at line 32 of file TMatacq.h.

comp_trise

double TMatacq::comp_trise[1200]

private

Definition at line 32 of file TMatacq.h.

endsample

int TMatacq::endsample

private

Definition at line 16 of file TMatacq.h.

firstsample

int TMatacq::firstsample

private

Definition at line 16 of file TMatacq.h.

Referenced by getPt1().

fNsamples

int TMatacq::fNsamples

private

Definition at line 12 of file TMatacq.h.

fNum_samp_aft_max

int TMatacq::fNum_samp_aft_max

private

Definition at line 14 of file TMatacq.h.

fNum_samp_bef_max

int TMatacq::fNum_samp_bef_max

private

Definition at line 13 of file TMatacq.h.

fv1

double TMatacq::fv1[NSPARAB]

private

Definition at line 22 of file TMatacq.h.

fv2

double TMatacq::fv2[NSPARAB]

private

Definition at line 22 of file TMatacq.h.

fv3

double TMatacq::fv3[NSPARAB]

private

Definition at line 22 of file TMatacq.h.

lastsample

int TMatacq::lastsample

private

Definition at line 16 of file TMatacq.h.

Referenced by getPt2().

level1

double TMatacq::level1

private

Definition at line 19 of file TMatacq.h.

level2

double TMatacq::level2

private

Definition at line 19 of file TMatacq.h.

level3

double TMatacq::level3

private

Definition at line 19 of file TMatacq.h.

meantrise

double TMatacq::meantrise

private

Definition at line 28 of file TMatacq.h.

nevlasers

int TMatacq::nevlasers

private

Definition at line 30 of file TMatacq.h.

nevmtq0

int TMatacq::nevmtq0

private

Definition at line 30 of file TMatacq.h.

nevmtq1

int TMatacq::nevmtq1

private

Definition at line 30 of file TMatacq.h.

NMAXSAMP

constexpr int TMatacq::NMAXSAMP = 100

static

Definition at line 8 of file TMatacq.h.

nsigcut

double TMatacq::nsigcut

private

Definition at line 18 of file TMatacq.h.

nslide

int TMatacq::nslide

private

Definition at line 34 of file TMatacq.h.

NSPARAB

constexpr int TMatacq::NSPARAB = 16

static

Definition at line 9 of file TMatacq.h.

pkval

double TMatacq::pkval

private

Definition at line 25 of file TMatacq.h.

Referenced by getTimpeak().

presample

int TMatacq::presample

private

Definition at line 16 of file TMatacq.h.

samplemax

int TMatacq::samplemax

private

Definition at line 16 of file TMatacq.h.

Referenced by getPtm().

sigbl

double TMatacq::sigbl

private

Definition at line 23 of file TMatacq.h.

Referenced by getsigBaseLine().

sigpkval

double TMatacq::sigpkval

private

Definition at line 25 of file TMatacq.h.

Referenced by getsigTimpeak().

sigtrise

double TMatacq::sigtrise

private

Definition at line 28 of file TMatacq.h.

slidingmean

double TMatacq::slidingmean

private

Definition at line 33 of file TMatacq.h.

Referenced by getSlide().

status

int TMatacq::status[1200]

private

Definition at line 31 of file TMatacq.h.

Referenced by dirstructure.Comparison::__make_image(), and dirstructure.Comparison::__repr__().

t

double TMatacq::t[NSPARAB]

private

Definition at line 21 of file TMatacq.h.

timeatmax

double TMatacq::timeatmax

private

Definition at line 24 of file TMatacq.h.

Referenced by getTimax().

trise

double TMatacq::trise

private

Definition at line 26 of file TMatacq.h.

Referenced by getTrise().

val

double TMatacq::val[NSPARAB]

private

Definition at line 21 of file TMatacq.h.

Referenced by value.Value::__eq__(), value.Value::__iadd__(), value.Value::__idiv__(), value.Value::__isub__(), value.Value::__str__(), and value.Value::relerr().

val_max

double TMatacq::val_max

private

Definition at line 23 of file TMatacq.h.

width20

double TMatacq::width20

private

Definition at line 27 of file TMatacq.h.

Referenced by getWidth20().

width50

double TMatacq::width50

private

Definition at line 27 of file TMatacq.h.

Referenced by getFwhm().

width80

double TMatacq::width80

private

Definition at line 27 of file TMatacq.h.

Referenced by getWidth80().