TVirtualSpectrometer.h

Go to the documentation of this file.

#ifndef __TVIRTUALSPECTROMETER_h
#define __TVIRTUALSPECTROMETER_h

#include <iostream>

#include <TROOT.h>
#include <TObject.h>
#include <TTree.h>

#define alignMode_SVD         0  
#define alignMode_Linear      1  
#define alignMode_Parabolic   2  
class TVirtualSpectrometer : public TObject {
public:
  TVirtualSpectrometer( Int_t nOrbitBPMs, Int_t nChicaneBPMs );
  virtual ~TVirtualSpectrometer();

  Int_t            alignSystem( Long64_t iEv1, Long64_t iEv2, Int_t mode );
  
  void             setMagnetLocations( Double_t z1, Double_t z2 ) { _zB1 = z1; _zB2 = z2; }

  void             setMagnetFields( Double_t Bdl1, Double_t Bdl2 ) { _Bdl1 = Bdl1; _Bdl2 = Bdl2; }

  Int_t            getChicaneResiduals( void );

  Double_t         getEnergy( Long64_t iev );


  virtual Int_t    setTree( TTree *t ) = 0;


  virtual Int_t    getEvent( Long64_t iev ) = 0; // should fill the local vars from the tree

  void             showAlignment( void );

protected:
  TTree           *_dataTree;    
  Int_t            _nOrbitBPMs;  
  Double_t        *_xOrbitPos;   
  Double_t        *_zOrbitPos;   
  Int_t            _nChicaneBPMs; 
  Double_t        *_xChicanePos;  
  Double_t        *_zChicanePos;  
private:
  Int_t            _alignMode;    
  Int_t            _getSVDCoefficients( Long64_t iEv1, Long64_t iEv2 );
  Int_t            _getAlignmentConstants( Long64_t iEv1, Long64_t iEv2 );

  Double_t        *_xOrbitAlignOffset;   
  Double_t        *_xOrbitCorrPos;       
  Double_t        *_xChicaneAlignOffset; 
  Double_t        *_xChicaneCorrPos;     
  Double_t        *_xChicaneResidual;    
  Double_t       **_cSVD;           // the SVD constants

  Double_t         _zB1;
  Double_t         _zB2;
  Double_t         _Bdl1;
  Double_t         _Bdl2;

  ClassDef( TVirtualSpectrometer, 1 )
};

#endif /* #ifndef __TVIRTUALSPECTROMETER_h */

---