SBSBPM.h

#ifndef SBSBPM_H 
#define SBSBPM_H

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// SBSBPM                                                                   //
//  Class for a BPM measuring signals from four antennas                     //
//  and  the position at the BPM (without phase correction)                  //
//  works with FADCs                                                         //
//  is good for unrastered beam, or to get average positions                 //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include <cstdlib>
#include <iostream> 
#include <vector>

#include "THaEvData.h"
#include "THaDetMap.h"
#include "THaBeamDet.h"

#include "VarDef.h"
#include "VarType.h"

#include "TVectorT.h"
#include "TMath.h"

using namespace std;

class SBSBPM : public THaBeamDet {

public:
	explicit SBSBPM( const char* name, const char* description = "",
                         THaApparatus* a = nullptr );
	virtual ~SBSBPM();

	virtual void       Clear( Option_t* ="" );
	virtual Int_t      Decode( const THaEvData& );
	virtual Int_t      Process();

	virtual TVector3 GetPosition()  const { return fPosition; }
	virtual TVector3 GetDirection()  const { return fDirection; }

	Double_t GetRawSignal0() {return fRawSignal(0);}
	Double_t GetRawSignal1() {return fRawSignal(1);}
	Double_t GetRawSignal2() {return fRawSignal(2);}
	Double_t GetRawSignal3() {return fRawSignal(3);}

	Double_t GetRotPosX() {return fRotPos(0); }
	Double_t GetRotPosY() {return fRotPos(1); }


protected:

	virtual Int_t  ReadDatabase( const TDatime& date );
	virtual Int_t  DefineVariables( EMode mode = kDefine );

	//  SBSBPM() {}
	//  SBSBPM( const SBSBPM& ) {}
	//  SBSBPM& operator=( const SBSBPM& ) { return *this; }

	TVectorD  fRawSignal;           // induced signal of the antennas
	TVectorD  fPedestals;
	TVectorD  fCorSignal;           // pedestal subtracted signal

	TVectorD  fRotPos;              // position in the BPM system, arbitrary units

	TMatrix fRot2HCSPos;     

	TVector3  fOffset;
	TVector3  fPosition;            // Beam position at the BPM (meters)
	TVector3  fDirection;           // Beam direction at the BPM (always points along z axis) 

	std::ofstream *fDebugFile;

	Int_t fNfired;
	Double_t fCalibRot;

	ClassDef(SBSBPM,0)          // Generic BPM class 
};

#endif