THcHallCSpectrometer.h

#ifndef ROOT_THcHallCSpectrometer
#define ROOT_THcHallCSpectrometer

//////////////////////////////////////////////////////////////////////////
//
// THcHallCSpectrometer
//
//////////////////////////////////////////////////////////////////////////

#include "THaSpectrometer.h"

#include <vector>

#include "TClonesArray.h"
#include "THaNonTrackingDetector.h"
#include "THcHitList.h"
#include "THcRawHodoHit.h"
#include "THcScintillatorPlane.h"
#include "THcDC.h"

//#include "THaTrackingDetector.h"
//#include "THcHitList.h"
#include "THcSpacePoint.h"
#include "THcDriftChamberPlane.h"
#include "THcDriftChamber.h"
#include "TMath.h"

#include "THaSubDetector.h"
#include "TClonesArray.h"
#include <iostream>
#include <fstream>

//#include "spdlog/spdlog.h"
//#include "spdlog/sinks/stdout_color_sinks.h" //support for stdout logging
//#include "spdlog/sinks/basic_file_sink.h" // support for basic file logging

#include "hcana/Logger.h"

//class THaScintillator;

class THcHallCSpectrometer : public hcana::ConfigLogging<THaSpectrometer> {
private:
  //std::shared_ptr<spdlog::logger> _logger;

public:
  THcHallCSpectrometer( const char* name, const char* description );
  virtual ~THcHallCSpectrometer();

  virtual Int_t   ReadDatabase( const TDatime& date );
  virtual void    EnforcePruneLimits();
  virtual void    CalculateTargetQuantities(THaTrack* track,Double_t& gbeam_y,Double_t&  xptar,Double_t& ytar,Double_t& yptar,Double_t& delta);
  virtual Int_t   FindVertices( TClonesArray& tracks );
  virtual Int_t   TrackCalc();
  virtual Int_t   BestTrackSimple();
  virtual Int_t   BestTrackUsingScin();
  virtual Int_t   BestTrackUsingPrune();
  virtual Int_t   TrackTimes( TClonesArray* tracks );

  virtual Int_t   Decode( const THaEvData& );

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

  Bool_t SetTrSorting( Bool_t set = kFALSE );
  Bool_t GetTrSorting() const;

  // Mass of nominal detected particle type
  Double_t GetParticleMass() const {return fPartMass; }
  Double_t GetBetaAtPcentral() const { return
      fPcentral/TMath::Sqrt(fPcentral*fPcentral+fPartMass*fPartMass);}

  virtual void AddEvtType(int evtype);
  virtual void SetEvtType(int evtype);
  virtual Bool_t IsMyEvent(Int_t evtype) const;
  virtual Int_t GetNumTypes() { return eventtypes.size(); };
  virtual Bool_t IsPresent() {return fPresent;};

  Bool_t InsideDipoleExitWindow(Double_t x_fp, Double_t xp_fp, Double_t y_fp, Double_t yp_fp);

protected:
  void InitializeReconstruction();

  Bool_t SHMSDipoleExitWindow(Double_t x_dip, Double_t y_dip);
  Bool_t HMSDipoleExitWindow(Double_t x_dip, Double_t y_dip);
  Bool_t fUseSHMSDipoleExitWindow;
  Bool_t fUseHMSDipoleExitWindow;
 //  Bool_t*      fKeep;
  //  Int_t*       fReject;

  Double_t     fPruneDipoleExit;
  Double_t     fPartMass;
  Double_t     fPruneXp;
  Double_t     fPruneYp;
  Double_t     fPruneYtar;
  Double_t     fPruneDelta;
  Double_t     fPruneBeta;
  Double_t     fPruneDf;
  Double_t     fPruneChiBeta;
  Double_t     fPruneFpTime;
  Double_t     fPruneNPMT;
  Double_t     fSatCorr;
  Double_t     fPruneSelect;

  Int_t        fGoodTrack;
  Int_t        fSelUsingScin;
  Int_t        fSelUsingPrune;
  Int_t        fNPlanes;
  Int_t        fNtracks;

  Double_t     fSelNDegreesMin;
  Double_t     fSeldEdX1Min;
  Double_t     fSeldEdX1Max;
  Double_t     fSelBetaMin;
  Double_t     fSelBetaMax;
  Double_t     fSelEtMin;
  Double_t     fSelEtMax;
  Double_t     fScin2XZpos;
  Double_t     fScin2XdZpos;
  Double_t     fScin2YZpos;
  Double_t     fScin2YdZpos;

  Double_t     fHodoCenter4, fHodoCenter3;
  Double_t     fScin2YSpacing, fScin2XSpacing;

  //  Int_t**   fHodScinHit;                // [4] Array

  THcHodoscope* fHodo;
  THcDC* fDC;

  Int_t fNReconTerms;
  struct reconTerm {
    Double_t Coeff[4];
    Int_t Exp[5];
    reconTerm() {
      for(Int_t i=0;i<4;i++) {
	Coeff[i] = 0.0;
      }
      for(Int_t i=0;i<5;i++) {
	Exp[i] = 0;
      }
    }
  };
  std::vector<reconTerm> fReconTerms;
  //  Double_t fReconCoeff[fMaxReconElements][4];
  //  Int_t fReconExponents[fMaxReconElements][5];
  Double_t fAngSlope_x;
  Double_t fAngSlope_y;
  Double_t fAngOffset_x;
  Double_t fAngOffset_y;
  Double_t fDetOffset_x;
  Double_t fDetOffset_y;
  Double_t fZTrueFocus;
  Double_t fThetaOffset; // Zero order term in yptar optics  matrix (rad)
  Double_t fPhiOffset; // Zero order term in xptar optics matrix (rad)
  Double_t fDeltaOffset; // Zero order term in delta optics matrix (%)
  Double_t fThetaCentralOffset; // Offset of Central spectrometer angle (rad)
  Double_t fOopCentralOffset; //Offset of central out-of-plane angle (rad)
  Double_t fPCentralOffset; // Offset Central spectrometer momentum (%)
  Double_t fTheta_lab; // Central spectrometer angle (deg)
  Double_t fPhi_lab; // Central spectrometer angle (deg)
  Double_t fMispointing_x; // Spectrometer Verticcal Mispointing
  Double_t fMispointing_y; // Spectrometer Horizontal Mispointing
  // For spectrometer central momentum use fPcentral in THaSpectrometer.h
  //  THaScintillator *sc_ref;  // calculate time track hits this plane

  // Flag for fProperties indicating that tracks are to be sorted by chi2
  static const UInt_t kSortTracks = BIT(16);

  std::vector<Int_t> eventtypes;
  Bool_t fPresent;

  ClassDef(THcHallCSpectrometer,0) //A Hall C Spectrometer
};

#endif