calib_barrel_atten.h

Go to the documentation of this file.

#ifndef CALIB_BARREL_ATTEN_H
#define CALIB_BARREL_ATTEN_H
 
#include "TofCalib.h"
 
const int nBarrelAtten = 2;
 
class calib_barrel_atten : public TofCalib {
public:
  calib_barrel_atten() : TofCalib( nBarrelAtten ) { m_name = string( "calib_barrel_atten" ); }
 
  void calculate_funcs( const Record* r ) {
    double z  = r->zrhit();
    double q1 = r->qleft();
    double q2 = r->qright();
 
    if ( q1 > 0.0 && q2 > 0.0 )
    {
      funcs[0] = 2.0 * z;
      funcs[1] = 1.0;
    }
    else
    {
      funcs[0] = 0.0;
      funcs[1] = 0.0;
    }
 
    return;
  }
 
  void calculate_y( const Record* r ) {
    double q1 = r->qleft();
    double q2 = r->qright();
    if ( q1 > 0.0 && q2 > 0.0 ) { y = log( q1 / q2 ); }
    else { y = 0.0; }
    return;
  }
 
  void updateData( RecordSet*& data ) {
    if ( data->size() > 0 )
    {
      std::vector<Record*>::iterator iter = data->begin();
      for ( ; iter != data->end(); iter++ )
      {
        double z        = ( *iter )->zrhit();
        double sintheta = ( *iter )->theta();
        double q1       = ( *iter )->qleft();
        double q2       = ( *iter )->qright();
        double latten   = 1.0 / X[0];
        double a1overa2 = exp( X[1] );
        double q0 =
            sintheta *
            ( q1 * exp( ( 115.0 - z ) / latten ) + q2 * exp( ( 115.0 + z ) / latten ) ) /
            ( 1.0 + a1overa2 );
        ( *iter )->setQ0( q0 );
      }
    }
    return;
  }
};
 
#endif