calib_etf_left.h

Go to the documentation of this file.

#ifndef CALIB_ETF_LEFT_H
#define CALIB_ETF_LEFT_H
 
#include "TofCalib.h"
 
// const int nEtfLeft = 2;   // event start time
const int nEtfLeft = 13; // end cap tof calibration
 
class calib_etf_left : public TofCalib {
public:
  calib_etf_left() : TofCalib( nEtfLeft ) { m_name = string( "calib_etf_left" ); }
 
  void calculate_funcs( const Record* r ) {
    double q = r->qleft();
    double z = r->zrhit();
 
    if ( q < 0.0 )
    {
      for ( int i = 0; i < nEtfLeft; i++ ) { funcs[i] = 1.0; }
    }
    else
    {
      funcs[0] = 1.0;
      //      funcs[1] = z;   // event start time
      funcs[1]  = 1. / sqrt( q );
      funcs[2]  = q;
      funcs[3]  = q * q;
      funcs[4]  = q * q * q;
      funcs[5]  = q * q * q * q;
      funcs[6]  = z / sqrt( q );
      funcs[7]  = 1.0 / q;
      funcs[8]  = z / q;
      funcs[9]  = z * q;
      funcs[10] = z;
      funcs[11] = z * z;
      funcs[12] = z * z * z;
    }
 
    return;
  }
 
  void calculate_y( const Record* r ) {
    y = r->tleft() - r->texp();
    return;
  }
 
  void updateData( RecordSet*& data ) {
    if ( data->size() > 0 )
    {
      std::vector<Record*>::iterator iter = data->begin();
      for ( ; iter != data->end(); iter++ )
      {
        calculate_funcs( ( *iter ) );
        double tcorr = 0.0;
        for ( int i = 0; i < X.num_row(); i++ ) { tcorr += X[i] * funcs[i]; }
        ( *iter )->settleft( ( *iter )->tleft() - tcorr - ( *iter )->texp() );
      }
    }
    return;
  }
};
 
#endif