AtVertexPropagator.h

Go to the documentation of this file.

#ifndef AtVertexPropagator_H
#define AtVertexPropagator_H
 
#include <Rtypes.h>
#include <TObject.h>
#include <TVector3.h>
 
#include <map>
#include <memory>
 
class TBuffer;
class TClass;
class TMemberInspector;
 
class AtVertexPropagator : public TObject {
 
private:
   static std::unique_ptr<AtVertexPropagator> fInstance;
 
   Int_t fGlobalEvtCnt;
   Int_t fBeamEvtCnt;
   Int_t fDecayEvtCnt;
 
   std::map<int, double> fTrackEn;
   std::map<int, double> fTrackAngle;
 
   Double_t fVx;
   Double_t fVy;
   Double_t fVz;
   Double_t fPx;
   Double_t fPy;
   Double_t fPz;
   Double_t fE;
   Double_t fBeamMass;
   Double_t fRndELoss;
   Double_t fBeamNomE;
   Double_t fInVx;
   Double_t fInVy;
   Double_t fInVz;
 
   Double_t fRecoilE;
   Double_t fRecoilA;
   Double_t fScatterE;
   Double_t fScatterA;
   Double_t fBURes1E;
   Double_t fBURes1A;
   Double_t fBURes2E;
   Double_t fBURes2A;
 
   Bool_t fIsValidKine;
   Bool_t fIsd2HeEvt;
   Int_t fAiso;
   Int_t fZiso;
   TVector3 fScatP;
   TVector3 fd2HeVtx;
   Double_t fExEjectile;
 
protected:
   AtVertexPropagator();
 
public:
   virtual ~AtVertexPropagator() = default;
 
   static AtVertexPropagator *Instance();
 
   void SetVertex(Double_t vx, Double_t vy, Double_t vz, Double_t invx, Double_t invy, Double_t invz, Double_t px,
                  Double_t py, Double_t pz, Double_t E);
   void SetBeamMass(Double_t m);
 
   void SetTrackEnergy(int trackID, double energy);
   void SetTrackAngle(int trackID, double angle);
 
   /*   void SetRecoilE(Double_t val);
      void SetRecoilA(Double_t val);
      void SetScatterE(Double_t val);
      void SetScatterA(Double_t val);
      void SetBURes1E(Double_t val); // Recoil(Scatt) breaks up. Residual 1
      void SetBURes1A(Double_t val);
      void SetBURes2E(Double_t val); // Recoil(Scatt) breaks up. Residual 2
      void SetBURes2A(Double_t val);
   */
 
   void SetRndELoss(Double_t eloss);
   void SetBeamNomE(Double_t ener);
   void ResetVertex();
   void SetMassNum(Int_t mnum);
   void SetAtomicNum(Int_t anum);
   void SetScatterP(TVector3 val);
   void SetScatterEx(Double_t val);
   void Setd2HeVtx(TVector3 val);
   void Setd2HeVtx(Double_t x0, Double_t y0, Double_t theta, Double_t phi);
 
   Int_t GetGlobalEvtCnt();
   Int_t GetBeamEvtCnt();
   Int_t GetDecayEvtCnt();
   Double_t GetBeamMass();
   Double_t GetVx();
   Double_t GetVy();
   Double_t GetVz();
   Double_t GetInVx();
   Double_t GetInVy();
   Double_t GetInVz();
   Double_t GetPx();
   Double_t GetPy();
   Double_t GetPz();
   Double_t GetEnergy();
   Double_t GetRndELoss();
   Double_t GetBeamNomE();
 
   Double_t GetTrackAngle(int trackID);
   Double_t GetTrackEnergy(int trackID);
 
   /*
   Double_t GetRecoilE();
   Double_t GetRecoilA();
   Double_t GetScatterE();
   Double_t GetScatterA();
   Double_t GetBURes1E(); // Recoil(Scatt) breaks up. Residual 1
   Double_t GetBURes1A();
   Double_t GetBURes2E(); // Recoil(Scatt) breaks up. Residual 2
   Double_t GetBURes2A();
   */
 
   Int_t GetMassNum();
   Int_t GetAtomicNum();
   TVector3 GetScatterP();
   Double_t GetScatterEx();
   TVector3 Getd2HeVtx();
 
   void IncGlobalEvtCnt();
   void IncBeamEvtCnt();
   void IncDecayEvtCnt();
 
   void SetValidKine(Bool_t val);
   Bool_t GetValidKine();
   Bool_t Getd2HeEvt();
 
   ClassDef(AtVertexPropagator, 2)
};
 
#endif