TrgEclDigitizer Class Reference

FAM module
More...

#include <TrgEclDigitizer.h>

Inheritance diagram for TrgEclDigitizer:

Collaboration diagram for TrgEclDigitizer:

Public Member Functions
	TrgEclDigitizer ()
	Constructor.
virtual	~TrgEclDigitizer ()
	Destructor.
void	setup (int)
	setup fam module
void	getTCHit (int)
	get TC Hits from Xtal hits
void	digitization01 (std::vector< std::vector< double >> &, std::vector< std::vector< double >> &)
	fit method, digi with 125ns interval
void	digitization02 (std::vector< std::vector< double >> &, std::vector< std::vector< double >> &)
	original no fit method, digi with 12ns interval
void	save (int)
	save fitting result into tables
void	setWaveform (int wave)
	Set flag of waveform table.
void	setFADC (int fadc)
	Set flag of waveform table.
double	FADC (int, double)
	FADC
double	SimplifiedFADC (int, double)
	FADC. More...
double	interFADC (double)
	Faster FADC using interpolation.
double	ShapeF (double, double, double, double, double, double, double)
	return shape using FADC function
double	ShapeF (double, double)
	Return shape using Simplified FADC.
double	u_max (double, double)
	Find max value between 2 vals;.

Private Attributes
double	TimeRange
	time range(default : -4000 ~ 4000 ns )
double	TCEnergy [576][80]
	TC Energy converted from Xtarl Energy [GeV].
double	TCTiming [576][80]
	TC Timing converted from Xtarl Timing [GeV].
double	TCEnergy_tot [576]
	TC Energy converted from Xtarl Energy [GeV].
double	TCTiming_tot [576]
	TC Timing converted from Xtarl Timing [GeV].
double	TCRawEnergy [576][60]
	Input TC energy[GeV].
double	TCRawTiming [576][60]
	Input TC timing[ns]
double	TCRawBkgTag [576][60]
	Input Beambackground tag
TrgEclMapping *	_TCMap
	Object of TC Mapping.
TrgEclDataBase *	_DataBase
	Object of DataBase.
std::vector< std::vector< double > >	MatrixParallel
	Noise Matrix of Parallel and Serial Noise. More...
std::vector< std::vector< double > >	MatrixSerial
	Noise Low triangle Matrix of Serial noise
double	TCBkgContribution [576][80]
	Beambackground contribution.
double	TCSigContribution [576][80]
	Signal contribution.
int	TCBeambkgTag [576][80]
	Beambackground tag.
int	_waveform
	Flag of waveform table.
double	WaveForm [576][64]
	TC Energy converted from Xtarl Energy [GeV].
int	_FADC
	Flag of choosing the method of waveform generation function 0: use simplifiedFADC, 1: use interFADC(interpolation)
int	_BeambkgTag
	Flag of saving beam background tag or not.

Detailed Description

FAM module

Definition at line 24 of file TrgEclDigitizer.h.

Member Function Documentation

SimplifiedFADC()

double SimplifiedFADC	(	int	flag_gen,
		double	timing
	)

FADC.

Faster FADC

Definition at line 838 of file TrgEclDigitizer.cc.

{
 
  //--------------------------------------
  //
  // o "timing" unit is [us]
  // o flag_gen = 0(=signal), 1(=parallel), 2(=serial)
  // o return value(PDF) is [GeV]
  // o Generate signal shape using a simplified function.
  //
  //
  //--------------------------------------
  double tsh, dd;
  static double tc, tc2, tsc, tris;
  static double amp, dft, as;
 
 
  //  int im, ij;
 
  static int ifir = 0;
 
  if (ifir == 0) {
 
    const double td =  0.10;  // diff time    (  0.10)
    const double t1 =  0.10;  // integ1 real  (  0.10)
    // b1 = 30.90;  // integ1 imag  ( 30.90)
    const double t2 =  0.01;  // integ2 real  (  0.01)
    // b2 = 30.01;  // integ2 imag  ( 30.01)
    double ts =  1.00;  // scint decay  (  1.00)
    dft = 0.600; // diff delay   ( 0.600)
    as = 0.548;  // diff frac    ( 0.548)
    //
    amp = 1.0;
    tris = 0.01;
    // ts0 = 0.0;
    tsc = ts;
    double    fm = 0;
    //
    int  im = 0;
    int  ij = 0;
    tc = 0;
    double tt = u_max(u_max(td, t1), u_max(t2, ts)) * 2;
    int flag_once = 0;
    while (flag_once == 0) {
      double dt = tt / 1000;
 
      double tm = 0;
      for (int j = 1; j <= 1000; j++) {
        tc2 = tc - dft;
        double fff =
          (ShapeF(tc, tsc) -
           ShapeF(tc, tris) * 0.01) -
          (ShapeF(tc2, tsc) -
           ShapeF(tc2, tris) * 0.01) * as;
        if (fff > fm) {
          fm = fff;
          tm = tc;
          im = j;
        }
        tc = tc + dt;
      }
      if (im >= 1000) {
        tt = 2 * tt;
        flag_once = 0;
        continue;
      }
      if (ij == 0) {
        ij = 1;
        tc = 0.99 * tm;
        dt = tm * 0.02 / 1000;
        flag_once = 0;
        continue;
      }
      flag_once = 1;
    }
    amp = 1.0 / fm;
    ifir = 1;
  }
  //
  //
  double pdf = 0;
  if (flag_gen == 0) {
    //-----<signal>
    tc2 = timing - dft;
    pdf = amp * (
            (ShapeF(timing, tsc) -
             ShapeF(timing, tris) * 0.01) -
            (ShapeF(tc2,    tsc) -
             ShapeF(tc2,      tris) * 0.01) * as);
  } else if (flag_gen == 1) {
    //-----<parallel>
    tc2 = timing - dft;
    tsh = 0.001;
    pdf = amp * (
            ShapeF(timing, tsh) -
            ShapeF(tc2,     tsh) * as);
    pdf = pdf * 0.001; // GeV
  } else {
    //-----<serial>
    tc2 = timing - dft;
    tsh = 0.001;
    pdf = amp * (
            ShapeF(timing, tsh) -
            ShapeF(tc2,    tsh) * as);
    //
    tc = timing - 0.01;
    if (tc < 0) { tc = 0; }
    dd = timing - tc;
    tc2 = tc - dft;
    pdf = (amp * (
             ShapeF(tc, tsh) -
             ShapeF(tc2, tsh) * as) - pdf) / dd;
    pdf = pdf * 0.001; // GeV
  }
 
  return pdf;
}

Member Data Documentation

MatrixParallel

std::vector<std::vector<double> > MatrixParallel

private

Noise Matrix of Parallel and Serial Noise.

Noise Low triangle Matrix of Parallel noise

Definition at line 96 of file TrgEclDigitizer.h.

---

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

• trg/ecl/include/TrgEclDigitizer.h
• trg/ecl/src/TrgEclDigitizer.cc