/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/MagneticField/MagneticField-00-01-38/MagneticField/MagFieldReader.h

Go to the documentation of this file.

#ifndef MAGFIELDREADER_H 
#define MAGFIELDREADER_H 1

// Include files
// from STL
#include <string>

// from Gaudi
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "BesTimerSvc/IBesTimerSvc.h"
#include "BesTimerSvc/BesTimerSvc.h"

// Forward references
class IMagneticFieldSvc;

class MagFieldReader : public Algorithm {
public:
  MagFieldReader( const std::string& name, ISvcLocator* pSvcLocator );

  virtual ~MagFieldReader( ){ }; 

  virtual StatusCode initialize();    
  virtual StatusCode execute   ();    
  virtual StatusCode finalize  ();    

protected:

private:

  StatusCode readField();
  
  // Pointer to the magnetic field service
  IMagneticFieldSvc* m_pIMF;

  // range and steps in z to get the magnetic field 
  double m_zMin, m_zMax, m_step;

  // range in x and y
  double m_xMin, m_xMax, m_yMin, m_yMax;

  NTuple::Tuple*        m_ntuple;
  NTuple::Tuple*        m_tuple1;
  NTuple::Tuple*        m_tuple2;
  NTuple::Tuple*        m_tuple3;
  NTuple::Tuple*        m_tuple4;
  NTuple::Item<float>   m_x, m_y, m_z, m_Bx, m_By, m_Bz, m_r;
  NTuple::Item<float>   m_x2, m_y2, m_z2, m_Bx2, m_By2, m_Bz2, m_r2;
  NTuple::Item<float>   m_x3, m_y3, m_z3, m_Bx3, m_By3, m_Bz3, m_r3, m_phi3;
  NTuple::Item<float> m_sigma_bx,m_sigma_by,m_sigma_bz;  
  NTuple::Item<float> m_time;

  std::string m_filename;

  IBesTimerSvc* m_timersvc;
  BesTimer* m_timer;
};
#endif // MAGFIELDREADER_H

---