GeoSwimApplication.h

Go to the documentation of this file.

//
// GeoSwimApplication
//
// January 20, 2006  M.Ishitsuka First version for GeoSwimmer
//                   ref. PTSimApplication.h by S. Kasahara 05/04
//

#ifndef GEOSWIMAPPLICATION_H
#define GEOSWIMAPPLICATION_H

#include <TLorentzVector.h>

#include "MCApplication/MCAppUser.h" // base class
#include "Validity/VldRange.h"
#include "GeoSwimmer/GeoSwimStack.h"

class TClonesArray;
class GeoSwimZCondition;

//----------GCTRAK
#define MAXMEC 30
//      PARAMETER (MAXMEC=30)
//      COMMON/GCTRAK/VECT(7),GETOT,GEKIN,VOUT(7),NMEC,LMEC(MAXMEC)
//     + ,NAMEC(MAXMEC),NSTEP ,MAXNST,DESTEP,DESTEL,SAFETY,SLENG
//     + ,STEP  ,SNEXT ,SFIELD,TOFG  ,GEKRAT,UPWGHT,IGNEXT,INWVOL
//     + ,ISTOP ,IGAUTO,IEKBIN, ILOSL, IMULL,INGOTO,NLDOWN,NLEVIN
//     + ,NLVSAV,ISTORY
typedef struct {
  Float_t  vect[7];
  Float_t  getot;
  Float_t  gekin;
  Float_t  vout[7];
  Int_t    nmec;
  Int_t    lmec[MAXMEC];
  Int_t    namec[MAXMEC];
  Int_t    nstep;
  Int_t    maxnst;
  Float_t  destep;
  Float_t  destel;
  Float_t  safety;
  Float_t  sleng;
  Float_t  step;
  Float_t  snext;
  Float_t  sfield;
  Float_t  tofg;
  Float_t  gekrat;
  Float_t  upwght;
  Int_t    ignext;
  Int_t    inwvol;
  Int_t    istop;
  Int_t    igauto;
  Int_t    iekbin;
  Int_t    ilosl;
  Int_t    imull;
  Int_t    ingoto;
  Int_t    nldown;
  Int_t    nlevin;
  Int_t    nlsav;
  Int_t    istory;
} Gctrak_t;

//----------GCFLAG
//      COMMON/GCFLAG/IDEBUG,IDEMIN,IDEMAX,ITEST,IDRUN,IDEVT,IEORUN
//     +        ,IEOTRI,IEVENT,ISWIT(10),IFINIT(20),NEVENT,NRNDM(2)
typedef struct {
  Int_t    idebug;
  Int_t    idemin;
  Int_t    idemax;
  Int_t    itest;
  Int_t    idrun;
  Int_t    idevt;
  Int_t    ieorun;
  Int_t    ieotri;
  Int_t    ievent;
  Int_t    iswit[10];
  Int_t    ifinit[20];
  Int_t    nevent;
  Int_t    nrndm[2];
} Gcflag_t;

class GeoSwimApplication : public MCAppUser {

 public:

  virtual Gctrak_t* Gctrak() const {return fGctrak;}
  // Constructors/Destructors
  GeoSwimApplication(const char* name, const char* title);
  virtual ~GeoSwimApplication();

  // Must supply the following, declared pure virtual in base class
  virtual void GeneratePrimaries() {}
  virtual void BeginEvent();
  virtual void BeginPrimary();
  virtual void PreTrack();
  virtual void Stepping();
  virtual void PostTrack();
  virtual void FinishPrimary();
  virtual void FinishEvent();

  GeoSwimStack* GetStack() const { return fGeoSwimStack; }

  void SetZCondition(GeoSwimZCondition& c);
  void SetDirection(Int_t idir);
  void SetHallSize(TVector3 min, TVector3 max);

  bool GetSatisfied();
  Double_t GetS();
  Double_t GetRange();
  
 private:
  // private methods
  bool IsExiting();

  // data members
  GeoSwimStack*        fGeoSwimStack;
  VldRange             fVldRange;
  Double_t             fZCondition;
  TLorentzVector       fInitMomentum;
  TLorentzVector       fInitPosition;
  Int_t                fDir;
  bool                 fSatisfied;
  Bool_t               fExitLiner;  // flag to determine when to stop tracking
  Int_t                fStep;
  Double_t             fS;            // Path
  Double_t             fRange;        // Range
  Float_t              fPreVect[7];

  TVector3             fHallMin;
  TVector3             fHallMax;

  Gctrak_t *fGctrak;          
  Gcflag_t *fGcflag;          
  ClassDef(GeoSwimApplication,0) // Interface to Monte Carlo application
};

#endif // GEOSWIMAPPLICATION_H

---