Amplifier Class Reference

#include <Amplifier.h>

Collaboration diagram for Amplifier:

Public Member Functions
	Amplifier ()
	~Amplifier ()
double	TThreshold () const
double	TFirst () const
double	AdcResponse () const
double	Charge () const
bool	PassedThreshold () const
void	Reset ()
void	SetTriggerElectron (double el)
void	SetIntegrationWindow (double win)
void	InitResponse (unsigned int bins, double binsize)
bool	AddClusters (const cluster_vec &clusters)

Private Types
typedef std::vector< double >	cluster_vec
typedef std::vector< double >::iterator	cluster_vec_it
typedef std::vector< double >::const_iterator	cluster_vec_const_it

Private Member Functions
void	InitAmplifierParameters ()
void	DoAdcResponse ()
double	Integral ()
void	AddCluster (int bin, double charge)

Static Private Member Functions
static double	ResponseFunction (double time)

Private Attributes
double	m_binsize = 0.0
double	m_triggerElectron = 0.0
double	m_threshold = 0.0
int	m_integrationWindow
double	m_responseMax = 0.0
cluster_vec	m_response
cluster_vec	m_signal
double	m_integral
double	m_adcResponse
cluster_vec_it	m_signal_th
cluster_vec_it	m_signal_stop
cluster_vec_const_it	m_cluster_stop
cluster_vec_const_it	m_cluster_begin
double	m_adcOffset = 0.0
double	m_adcFactor = 0.0
double	m_adcFraction = 0.0

Detailed Description

Definition at line 11 of file Amplifier.h.

Member Typedef Documentation

cluster_vec

typedef std::vector<double> Amplifier::cluster_vec

private

Definition at line 14 of file Amplifier.h.

cluster_vec_const_it

typedef std::vector<double>::const_iterator Amplifier::cluster_vec_const_it

private

Definition at line 16 of file Amplifier.h.

cluster_vec_it

typedef std::vector<double>::iterator Amplifier::cluster_vec_it

private

Definition at line 15 of file Amplifier.h.

Constructor & Destructor Documentation

Amplifier()

Amplifier::Amplifier

(

)

Definition at line 9 of file Amplifier.cxx.

                    : 
 m_integrationWindow(0),
 m_integral(0.),
 m_adcResponse(0.) 
{ 
  InitAmplifierParameters();
}

~Amplifier()

Amplifier::~Amplifier ( )

default

Member Function Documentation

AdcResponse()

double Amplifier::AdcResponse ( ) const

inline

Definition at line 90 of file Amplifier.h.

                                           { 
  return m_adcResponse; 
}

AddCluster()

void Amplifier::AddCluster ( int  bin, double  charge  )

private

Definition at line 154 of file Amplifier.cxx.

{
  // set start values for signal and response iterator
  cluster_vec_it response_it = m_response.begin();
  cluster_vec_it signal_it = m_signal.begin() + startbin;
 
  // fill signal vector until stop
  while(signal_it != m_signal_stop) {
    
    *signal_it += *response_it*charge;        // add response to signal
    if( *signal_it > m_threshold ) {   // check if signal above threshold
      
      // if new threshold found
      if( signal_it - m_signal_th < 0 || m_signal_th == m_signal.begin()){
    m_signal_th = signal_it;            // store new threshold
 
    /* if the bin (threshold + integration window) lies 
       within the signal vector, point the stop iterator to it.
       nessicary to make sure m_signal_stop never points beyond the m_signal.end() */
 
    if( m_signal.end()-signal_it > m_integrationWindow ){
      m_signal_stop = m_signal_th + m_integrationWindow;
      m_cluster_stop = m_cluster_begin + (m_signal_stop-m_signal.begin());
    }
      }
      ++signal_it;++response_it;             // go to next bin
      while(signal_it != m_signal_stop) {         // loop over the remaining bins
    *signal_it += *response_it*charge;
    ++signal_it;++response_it;
      }
      break;                     // exit loop
    }
    ++signal_it;++response_it;    
  }
}

AddClusters()

bool Amplifier::AddClusters

(

const cluster_vec &

clusters

)

Definition at line 132 of file Amplifier.cxx.

{
  if(clusters.empty()){   // no clusters, no signal!
    return false;
  }
 
  // reset cluster iterators 
  m_cluster_begin = clusters.begin(); 
  m_cluster_stop = clusters.end();         // stop iterator can be changed in the AddCluster-routine
  cluster_vec_const_it cit = m_cluster_begin;  
  int bin(0);
  for( ; cit!=m_cluster_stop;++cit){       // iterate over all clusters
    if(*cit == 0.) continue;         // check if the cluster size in the current bin is not zero
    bin = cit - m_cluster_begin;           // calculate the corresponding bin number
    AddCluster(bin,*cit);            // add response of cluster to total detector response
  }
  Integral();                        // calculate the Q-integral
  DoAdcResponse();                   // get the ADC response
 
  return PassedThreshold();
}

Charge()

double Amplifier::Charge ( ) const

inline

Definition at line 93 of file Amplifier.h.

                                      {
  return m_integral; 
}

DoAdcResponse()

void Amplifier::DoAdcResponse ( )

private

Definition at line 201 of file Amplifier.cxx.

{
  if(m_integral > 0){
    m_adcResponse = m_adcOffset + m_adcFactor*(log(m_integral)/log(m_adcFraction)-1);
  }else{
    m_adcResponse = -1000;
  }
}

InitAmplifierParameters()

void Amplifier::InitAmplifierParameters ( )

private

Definition at line 20 of file Amplifier.cxx.

{
  m_threshold = 0.;
  m_responseMax = 0.;
  m_triggerElectron = 20.;
  m_adcOffset = 30.;
  m_adcFactor = 90.;
  m_adcFraction = 10.;
  m_binsize = 1.;
  //std::cout << "Amplifier Threshold " << m_triggerElectron << std::endl;
} 

InitResponse()

void Amplifier::InitResponse	(	unsigned int	bins,
		double	binsize
	)

Definition at line 97 of file Amplifier.cxx.

{
  m_binsize = binsize;  
  m_integrationWindow =  (int)(20./binsize)-1 ;
 
  m_response.resize(bins);
  m_signal.resize(bins);
  //std::cout << "Amplifier Initialize: new response vector with " << m_response.size() << " bins" << std::endl;
 
  cluster_vec_it it = m_response.begin();
  double i(0.);
  double resp(0),max(0.);//,integral(0.);
  while( it != m_response.end() ){
    resp = m_binsize*ResponseFunction(m_binsize*i);
    if(resp>max) max = resp;
    //integral+=resp;
    *it = resp;
    ++it;++i;
  }  
  m_responseMax = max;
  m_threshold = m_triggerElectron*max;
  //std::cout << "max " << max << " integral " << integral << " threshold " << m_threshold <<std::endl;
  Reset();
}

Integral()

double Amplifier::Integral ( )

private

Definition at line 190 of file Amplifier.cxx.

{
  cluster_vec_it s = m_signal_th;
  m_integral = 0.;
  while(s != m_signal_stop){         // sum charge within integration window
    m_integral += *s;
    ++s;
  }
  return m_integral;
} 

PassedThreshold()

bool Amplifier::PassedThreshold ( ) const

inline

Definition at line 105 of file Amplifier.h.

                                              { 
  return m_signal.begin() != m_signal_th; 
}

Reset()

void Amplifier::Reset

(

)

Definition at line 122 of file Amplifier.cxx.

{ 
  m_integral = -10000.;     
  std::vector<double>::iterator dit = m_signal.begin();
  for( ; dit != m_signal.end(); ++dit) 
    *dit = 0.;
  m_signal_th = m_signal.begin();
  m_signal_stop = m_signal.end();
}

ResponseFunction()

double Amplifier::ResponseFunction ( double  time )

staticprivate

Definition at line 32 of file Amplifier.cxx.

{
  /*
    amplifier f(t) = (t/tau)**3 exp(-t/tau) u(t)
    
    signal    I(t) = I0/(t+t0) u(t)
    
    u(t) step function 
    tau 15 ns, t0=12 nsec, I0=1
 
  */               
  //double TP=20.;
  double TP=15.;
  double T0=12./2.;
  /*
    FI1=1./3.
    FI2=2./3.
    signal = ( FI1 exp(t/t0)/t0 + FI2 exp(-t/(9*t0))/(9*t0))
    ion/electron fraction  electron delta function ion signal
  */ 
  double FI1=1./3.;
  double FI2=2./3.;
  double face=2.0; 
 
  /* third order amplifier
 
  f(t)= u(t) * t/tau**3* exp(-t/tau)/3!/tau
  I(t)= u(t) * exp(-t/t0)/ t0 
 
  */
  double TAU=TP/3.;
  double a=1./TAU;
  double b =1./T0;
  double a4 = a*a*a*a;
  double ab4 = a-b;
  ab4 = ab4*ab4*ab4*ab4;
  double dt = (a-b)*time;
  double dt2 = dt*dt;
  double dt3 = dt2*dt;
  // double dt4 = dt2*dt2; // not used
 
  double F3 = a4*time*time*time*exp(-a*time)/6.;
  double xnor = b*a4/ab4;
  double res31= xnor*(exp(-b*time) - exp(-a*time)  
              -dt*exp(-a*time)
              -dt2/2.*exp(-a*time)
              -dt3/6.*exp(-a*time));
 
  b =1./T0/9.;
  ab4 = a-b;
  ab4 = ab4*ab4*ab4*ab4;
  dt = (a-b)*time;
  dt2 = dt*dt;
  dt3 = dt2*dt;
  // dt4 = dt2*dt2; // not used
  xnor = b*a4/ ab4;
  double res32= xnor*(exp(-b*time) - exp(-a*time)  
              -dt*exp(-a*time)
              -dt2/2.*exp(-a*time)
              -dt3/6.*exp(-a*time));
  double resolution=(face*(res31*FI1+res32*FI2)+F3)/(1.+face);
 
  return resolution;
}

SetIntegrationWindow()

void Amplifier::SetIntegrationWindow ( double  win )

inline

Definition at line 102 of file Amplifier.h.

                                                       { 
  m_integrationWindow = (int)(win/m_binsize)-1; 
}

SetTriggerElectron()

void Amplifier::SetTriggerElectron ( double  el )

inline

Definition at line 97 of file Amplifier.h.

                                                    {
  m_triggerElectron = el;
  m_threshold = el*m_responseMax; 
  //std::cout << "Amplifier: Changing threshold: now triggering on " << el << "th electron "        << "new threshold " << m_threshold << std::endl;
}

TFirst()

double Amplifier::TFirst ( ) const

inline

Definition at line 84 of file Amplifier.h.

                                          { 
  cluster_vec_const_it cit = m_signal.begin();
  while( *cit == 0 )
    ++cit;
  return m_binsize*(cit-m_signal.begin()); ;
}

TThreshold()

double Amplifier::TThreshold ( ) const

inline

Definition at line 81 of file Amplifier.h.

                                          { 
  return m_binsize*(m_signal_th-m_signal.begin()); 
}

Member Data Documentation

m_adcFactor

double Amplifier::m_adcFactor = 0.0

private

Definition at line 76 of file Amplifier.h.

m_adcFraction

double Amplifier::m_adcFraction = 0.0

private

Definition at line 77 of file Amplifier.h.

m_adcOffset

double Amplifier::m_adcOffset = 0.0

private

Definition at line 75 of file Amplifier.h.

m_adcResponse

double Amplifier::m_adcResponse

private

Definition at line 65 of file Amplifier.h.

m_binsize

double Amplifier::m_binsize = 0.0

private

Definition at line 54 of file Amplifier.h.

m_cluster_begin

cluster_vec_const_it Amplifier::m_cluster_begin

private

Definition at line 71 of file Amplifier.h.

m_cluster_stop

cluster_vec_const_it Amplifier::m_cluster_stop

private

Definition at line 70 of file Amplifier.h.

m_integral

double Amplifier::m_integral

private

Definition at line 64 of file Amplifier.h.

m_integrationWindow

int Amplifier::m_integrationWindow

private

Definition at line 57 of file Amplifier.h.

m_response

cluster_vec Amplifier::m_response

private

Definition at line 62 of file Amplifier.h.

m_responseMax

double Amplifier::m_responseMax = 0.0

private

Definition at line 60 of file Amplifier.h.

m_signal

cluster_vec Amplifier::m_signal

private

Definition at line 63 of file Amplifier.h.

m_signal_stop

cluster_vec_it Amplifier::m_signal_stop

private

Definition at line 69 of file Amplifier.h.

m_signal_th

cluster_vec_it Amplifier::m_signal_th

private

Definition at line 68 of file Amplifier.h.

m_threshold

double Amplifier::m_threshold = 0.0

private

Definition at line 56 of file Amplifier.h.

m_triggerElectron

double Amplifier::m_triggerElectron = 0.0

private

Definition at line 55 of file Amplifier.h.

---

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

• Amplifier.h
• Amplifier.cxx