#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:

Public Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ cluster_vec

◆ cluster_vec_const_it

◆ cluster_vec_it

Constructor & Destructor Documentation

◆ Amplifier()

◆ ~Amplifier()

Member Function Documentation

◆ AdcResponse()

◆ AddCluster()

◆ AddClusters()

◆ Charge()

◆ DoAdcResponse()

◆ InitAmplifierParameters()

◆ InitResponse()

◆ Integral()

◆ PassedThreshold()

◆ Reset()

◆ ResponseFunction()

◆ SetIntegrationWindow()

◆ SetTriggerElectron()

◆ TFirst()

◆ TThreshold()

Member Data Documentation

◆ m_adcFactor

◆ m_adcFraction

◆ m_adcOffset

◆ m_adcResponse

◆ m_binsize

◆ m_cluster_begin

◆ m_cluster_stop

◆ m_integral

◆ m_integrationWindow

◆ m_response

◆ m_responseMax

◆ m_signal

◆ m_signal_stop

◆ m_signal_th

◆ m_threshold

◆ m_triggerElectron