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

Go to the documentation of this file.

#ifndef MAGNETICFIELDSVC_H
#define MAGNETICFIELDSVC_H 

// Include files
#include <vector>
#include <string>
#include "CLHEP/Geometry/Vector3D.h"
#include "CLHEP/Geometry/Point3D.h"
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
//  backwards compatibility will be enabled ONLY in CLHEP 1.9
typedef HepGeom::Point3D<double> HepPoint3D;
typedef HepGeom::Vector3D<double> HepVector3D;
#endif
#include "MagneticField/ConnectionDB.h"

#ifndef BEAN
#include "GaudiKernel/Service.h"
//#include "GaudiKernel/IMagneticFieldSvc.h"
#include "MagneticField/IMagneticFieldSvc.h"
#include "GaudiKernel/IIncidentListener.h"
using namespace CLHEP;

// Forward declarations
template <class TYPE> class SvcFactory;

class IMagneticFieldSvc;
class MucMagneticField;
class IDataProviderSvc;
#else
class MucMagneticField;
typedef bool StatusCode;
#endif

#ifndef BEAN
class MagneticFieldSvc : public Service,
                         virtual public IMagneticFieldSvc,
                         virtual public IIncidentListener {

protected:
  
  MagneticFieldSvc( const std::string& name, ISvcLocator* svc );
  
  virtual ~MagneticFieldSvc();

public:

  virtual StatusCode initialize();
  
  void init_params();

  virtual StatusCode finalize();
  
  virtual StatusCode queryInterface( const InterfaceID& riid, 
                                     void** ppvInterface      );
  
  virtual const InterfaceID& type() const { return IMagneticFieldSvc::interfaceID(); };


  // Incident handler
  void handle(const Incident&);

#else
// -------------------------- BEAN ------------------------------------
class MagneticFieldSvc {

public:
  static MagneticFieldSvc* instance() {
    return (m_field) ? m_field : (m_field=new(std::nothrow)MagneticFieldSvc());
  }

private:
  MagneticFieldSvc();
  virtual ~MagneticFieldSvc();

public:
  void SetPath(std::string new_path)            {path = new_path;}
  std::string GetPath()                         {return path;}

  void GridDistance( int gridDistance = 5)      {m_gridDistance = gridDistance;} 
  void RunMode( int runmode = 2)                {m_runmode = runmode;}
  void IfRealField( bool ifRealField = true)    {m_ifRealField = ifRealField;}
  void OutLevel( int outlevel = 1)              {m_outlevel = outlevel;}

  void Cur_SCQ1_55( double Cur_SCQ1_55 = 349.4) {m_Cur_SCQ1_55 = Cur_SCQ1_55;}
  void Cur_SCQ1_89( double Cur_SCQ1_89 = 426.2) {m_Cur_SCQ1_89 = Cur_SCQ1_89;}
  void Cur_SCQ2_10( double Cur_SCQ2_10 = 474.2) {m_Cur_SCQ2_10 = Cur_SCQ2_10;}

  void UseDBFlag( bool useDB = true)            {m_useDB = useDB;}

  bool init_mucMagneticField();
  virtual StatusCode initialize();
  void init_params(int run);

  // handler for new run
  void handle(int new_run);
#endif

  virtual StatusCode fieldVector( const HepPoint3D& xyz, 
                                  HepVector3D& fvec ) const;

  virtual StatusCode uniFieldVector( const HepPoint3D& xyz,
                                  HepVector3D& fvec ) const;

  virtual double getReferField();

  virtual bool ifRealField() const;


private:
#ifndef BEAN
  friend class SvcFactory<MagneticFieldSvc>;
#else
  static MagneticFieldSvc* m_field;
#endif

  void fieldGrid( const HepPoint3D& xyz, HepVector3D& fvec ) const;
  void fieldGrid_TE( const HepPoint3D& xyz, HepVector3D& fvec ) const;

  StatusCode parseFile( );       
  StatusCode parseFile_TE( );       
  
  std::string path;             // path to data files

  std::string m_filename;        
  std::string m_filename_TE;        
  int m_runmode;                 
  int m_gridDistance;            
  int m_outlevel;
  double m_scale;

  double m_Cur_SCQ1_55;
  double m_Cur_SCQ1_89;
  double m_Cur_SCQ2_10;

  std::vector<double> m_Q;       
  std::vector<double> m_P;       
  std::vector<double> m_Q_TE;       
  std::vector<double> m_P_TE;       
  std::vector<double> m_Q_1;       
  std::vector<double> m_P_1;       
  std::vector<double> m_Q_2;       
  std::vector<double> m_P_2;       
  double m_Dxyz[3];              
  int    m_Nxyz[3];              
  double m_max_FL[3];    
  double m_min_FL[3];
  //for tof and emc
  double m_Dxyz_TE[3];              
  int    m_Nxyz_TE[3];              
  double m_max_FL_TE[3];   
  double m_min_FL_TE[3];

  double m_zOffSet;              
  double m_zOffSet_TE;              
  MucMagneticField* m_Mucfield;
  double m_zfield;
  bool m_ifRealField;

  //property
  bool m_useDB;
 
  bool m_turnOffField;
  bool m_uniField;

#ifndef BEAN
  IDataProviderSvc* m_eventSvc;
#endif

  //database
  FieldDBUtil::ConnectionDB* m_connect_run;
};

#endif  // MAGNETICFIELDSVC_H

---