HypStatDecay.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                 HypStatDecay header file                  -----
// -----               Created 30/10/08  by A.Sanchez                 -----
// -------------------------------------------------------------------------

#ifndef HYPSTATDECAY_HH
#define HYPSTATDECAY_HH


//#include "THParticle.h"
#include <vector>
#include "TRandom.h"

class TClonesArray;
class TFile;
class TTree;

class HypStatDecay 
{

 public: 

  HypStatDecay();


  HypStatDecay(const Char_t* fileName);


  virtual ~HypStatDecay();


  void GetFragment(Int_t primGen); 
  void GetData(Int_t tr,TVector3 &p,Int_t &pid,Double_t &mass);
  void GetAZH(Int_t ion,Int_t &AI,Int_t &ZI,Int_t &L);
  
  
  Int_t  GetNtr() const { return fNtr;  } 
  //TClonesArray* GetFrag() const { return fFrag; }
  
  void SetNtr(Int_t t)          { fNtr    = t; };
  //void SetFrag(TClonesArray* frag)          { fFrag    = frag; };
  
  private:
   
  Int_t iEvent;   
  Int_t fev;
 const Char_t* fFileName;   

 TRandom r;
  double she;
  int an, zn, hn;
  
  
  Int_t fNtr;
  std::vector<int>fPid;
  TVector3 fMom;
  TLorentzVector v,vN;
  TVector3 tvL,tNvL;

  std::vector<double> fPx;
  std::vector<double> fPy;
  std::vector<double> fPz;
  //Double fPhx[1000],fPhy[1000],fPhz[1000];
  std::vector<double> fMass;
  std::vector<double>fEx;
  std::vector<int> A;
  std::vector<int>Z;
  std::vector<int>H;

  TLorentzVector GetPgCMSLab(float mass,float& Delta,TVector3& PL,TRandom& rd);
  TLorentzVector GetPNuCMSLab(float mass,float& Delta,TVector3& PL,TRandom& rd);
  Double_t GetEtot(Float_t mass,TVector3 P);

  void CloseInput();
  
  ClassDef(HypStatDecay,2);

};

#endif