AtTPCIonDecay.h

Go to the documentation of this file.

#ifndef AtTPCIonDecay_H
#define AtTPCIonDecay_H
 
#include <FairGenerator.h>
#include <FairIon.h>
#include <FairParticle.h>
 
#include <Rtypes.h>
#include <TString.h>
 
#include <memory>
#include <vector>
 
class FairPrimaryGenerator;
class TBuffer;
class TClass;
class TMemberInspector;
 
class AtTPCIonDecay : public FairGenerator {
 
public:
   AtTPCIonDecay();
 
   // AtTPCIonDecay(const char* name,std::vector<Int_t> *z,std::vector<Int_t> *a,std::vector<Int_t> *q, Int_t mult,
   // std::vector<Double_t> *px,
   //   std::vector<Double_t>* py,std::vector<Double_t> *pz, std::vector<Double_t> *mass , Double_t ResEner, Int_t ZB,
   //   Int_t AB, Double_t PxB, Double_t PyB, Double_t PzB, Double_t BMass, Double_t TMass);
 
   AtTPCIonDecay(std::vector<std::vector<Int_t>> *z, std::vector<std::vector<Int_t>> *a,
                 std::vector<std::vector<Int_t>> *q, std::vector<std::vector<Double_t>> *mass, Int_t ZB, Int_t AB,
                 Double_t BMass, Double_t TMass, Double_t ExEnergy, std::vector<Double_t> *SepEne);
 
   AtTPCIonDecay(const AtTPCIonDecay &);
 
   AtTPCIonDecay &operator=(const AtTPCIonDecay &) { return *this; }
 
   virtual Bool_t ReadEvent(FairPrimaryGenerator *primGen);
   void SetSequentialDecay(Bool_t var) { fIsSequentialDecay = var; }
 
   virtual ~AtTPCIonDecay() = default;
 
private:
   static Int_t fgNIon;                                     
   std::vector<Int_t> fMult;                                // Multiplicity per decay channel
   Int_t fNbCases;                                          // Number of decay channel
   std::vector<Double_t> fPx, fPy, fPz;                     // Momentum components [GeV] per nucleon
   std::vector<std::vector<Double_t>> fMasses;              // Masses of the N products
   Double_t fVx, fVy, fVz;                                  // Vertex coordinates [cm]
   std::vector<std::vector<std::unique_ptr<FairIon>>> fIon; // Pointer to the FairIon to be generated
   std::vector<std::vector<std::unique_ptr<FairParticle>>> fParticle;
   std::vector<std::vector<TString>> fPType;
   std::vector<Int_t> fQ; // Electric charge [e]
   // std::vector<Int_t> fA;
   // std::vector<Int_t> fZ;
   Double_t fBeamEnergy{};      // Residual beam energy for phase calculation
   Double_t fBeamEnergy_buff{}; // Residual beam energy for phase calculation
   Int_t fZBeam{};
   Int_t fABeam{};
   Double_t fPxBeam;
   Double_t fPyBeam;
   Double_t fPzBeam;
   Bool_t fIsDecay{};
   Double_t fBeamMass{};
   Double_t fTargetMass{};
   Double_t fExEnergy{};
   std::vector<Double_t> fSepEne;
   Bool_t fIsSequentialDecay{}; //<! True if the decay generator is to be used after a reaction generator.
 
   ClassDef(AtTPCIonDecay, 3)
};
 
#endif