StRoot: St_geant_Maker/St_geant

00001 // $Id: St_geant_Maker.h,v 1.51 2011/09/11 20:57:14 fisyak Exp $
00002 // $Log: St_geant_Maker.h,v $
00003 // Revision 1.51  2011/09/11 20:57:14  fisyak
00004 // Add kinematics definition via MuDst, Clean up
00005 //
00006 // Revision 1.50  2010/08/10 16:35:33  fisyak
00007 // Add initialization of starsim parameter tables after opening zebra-file
00008 //
00009 // Revision 1.49  2010/05/27 13:36:14  fisyak
00010 // 3rd attemp to synchronize mag.field. Now take care that the maker can be not Active and do InitRun in Work
00011 //
00012 // Revision 1.47  2010/05/24 15:38:40  fisyak
00013 // Move geometry and mag.field initialization from Init into InitRun in order to allow mag. field settings from StMagFMaker::InitRun
00014 //
00015 // Revision 1.45  2010/05/10 14:19:52  fisyak
00016 // move geometry load from Init to InitRun in order to allow MagF maker to set mag.field
00017 //
00018 // Revision 1.44  2009/12/31 00:02:59  perev
00019 // Add the material name to the volume name
00020 //
00021 // Revision 1.43  2008/07/30 15:04:36  fisyak
00022 // Remove custom SetDebug, fix bug #1252
00023 //
00024 // Revision 1.42  2007/07/12 20:36:03  fisyak
00025 // Add includes for ROOT 5.16
00026 //
00027 // Revision 1.41  2007/07/12 20:35:30  fisyak
00028 // Add includes for ROOT 5.16
00029 //
00030 // Revision 1.40  2007/04/26 15:51:31  fisyak
00031 // Move creation of TGiant3 in ctor (fix byg 942)
00032 //
00033 // Revision 1.39  2007/03/03 00:35:43  fine
00034 // Fix the leak of the ROOT objects and introduce the method to return the source code filename  for the arbitrary geometry node
00035 //
00036 // Revision 1.38  2005/11/22 23:13:24  fisyak
00037 // Add default kinematics if there is no input fiels and if maker is active
00038 //
00039 // Revision 1.37  2005/10/06 19:23:07  fisyak
00040 // Add set date/time from fz-file
00041 //
00042 // Revision 1.36  2005/08/29 21:47:09  fisyak
00043 // Changes for VMC
00044 //
00045 // Revision 1.35  2005/04/13 22:27:11  fisyak
00046 // Add Hit description extractor (AgstHits)
00047 //
00048 // Revision 1.34  2005/02/07 21:09:20  fisyak
00049 // rename antique TGeant3 to TGiant3
00050 //
00051 // Revision 1.33  2004/02/10 23:16:34  potekhin
00052 // First version of Ag2Geom
00053 //
00054 // Revision 1.32  2003/10/01 23:54:08  potekhin
00055 // Added a declaration a a pointer to the structure geom_gdat,
00056 // needed for the propagation the GEANT run data --
00057 // geometry tag and field scale.
00058 //
00059 // Revision 1.31  2003/09/10 19:47:47  perev
00060 // ansi corrs
00061 //
00062 // Revision 1.30  2002/11/01 03:17:42  fine
00063 // the previous version has been restored. No need of the special flag
00064 //
00065 // Revision 1.28  2002/03/12 21:22:39  fisyak
00066 // Set only one StEvtHddr as default option (due to Embedding)
00067 //
00068 // Revision 1.27  2001/06/01 03:03:57  perev
00069 // overloaded GetDataSet -> FindDataSet
00070 //
00071 // Revision 1.26  2001/05/31 16:06:12  perev
00072 // hiding
00073 //
00074 // Revision 1.25  2000/09/23 03:07:13  fine
00075 // class TShape; statement introduced to fix a side aeffect of TDataSet clean up
00076 //
00077 // Revision 1.24  2000/03/26 02:43:22  fine
00078 //  adjusted to ROOT 2.24
00079 //
00080 // Revision 1.23  2000/02/07 18:58:38  fisyak
00081 // Set default NwGeant = 20 000 000 words
00082 //
00083 // Revision 1.22  2000/02/03 19:34:42  fisyak
00084 // Clean up St_geant_Maker::Init, move its parameters to ctor
00085 //
00086 // Revision 1.21  2000/02/03 16:15:54  fisyak
00087 // Add Kathy's histograms
00088 //
00089 
00090 #ifndef STAR_St_geant_Maker
00091 #define STAR_St_geant_Maker
00092 
00094 //                                                                      //
00095 // St_geant_Maker virtual base class for Maker                          //
00096 //                                                                      //
00098 #include "StMaker.h"
00099 //#define DetectorIndex
00100 #ifdef DetectorIndex
00101 #include "TArrayI.h"
00102 #endif
00103 class TVolume;
00104 class TGiant3;
00105 class TRotMatrix;
00106 class TH1F;
00107 class TShape;
00108 class TGeoVolume;
00109 class St_geom_gdat;
00110 class TFileSet;
00111 
00112 class St_geant_Maker : public StMaker {
00113 protected:
00114   Int_t  fNwGeant;     // No. of words in GCBANK common block
00115   Int_t  fNwPaw;       // No. of words in PAWC  common block
00116   Int_t  fIwType;      // HIGZ interface (=0 no HIGZ)
00117   TDataSet*   fVolume; 
00118   TGeoVolume* fTopGeoVolume; 
00119   TString fInputFile; // 
00120   TFileSet  *fGeoDirectory; // the pointer the STAR geometry source code
00121   StEvtHddr *fEvtHddr;
00122   Bool_t   fRemakeGeometry;// Flag whether we need to remake the geometry
00123   virtual TShape  *MakeShape(TString *name, Int_t ivo);
00124   virtual TVolume *MakeVolume(TString *name, Int_t ivo, Int_t Nlevel, Int_t *Names, Int_t *Numbers);
00125   virtual void ClearRootGeoms();
00126  private:
00127 
00128   St_geom_gdat *m_geom_gdat;
00129  public: 
00130                   St_geant_Maker(const char *name="geant",
00131                                  Int_t nwgeant=20,Int_t nwpaw=0, Int_t iwtype=0);
00132    virtual       ~St_geant_Maker(){};
00133    virtual Int_t  Finish(){SafeDelete(m_DataSet); return kStOK;}
00134    virtual Int_t  Init();
00135    virtual Int_t  InitRun(Int_t run);
00136    virtual void   SetDateTime(int idat=0,int itim=0);//
00137            void   SetFieldOpt(const char *opt) {mFieldOpt = opt;}
00138    virtual void   Do(const Char_t *option = "dcut cave x 0.1 10 10 0.03 0.03"); // *MENU 
00139    virtual void   Draw(const char* opt="IN");
00140    virtual Int_t  Make();
00141    virtual void   LoadGeometry (const Char_t *option = "detp geometry field_only");  // *MENU
00142    virtual void   SetNwGEANT (Int_t n=2) {fNwGeant = n;}
00143    virtual void   SetNwPAW   (Int_t n=0) {fNwPaw   = n;}
00144    virtual void   SetIwtype  (Int_t n=0) {fIwType  = n;}
00145    virtual Int_t  GetNwGEANT () {return fNwGeant;}
00146    virtual Int_t  GetNwPAW   () {return fNwPaw  ;}
00147    virtual Int_t  GetIwtype  () {return fIwType ;}
00148    Bool_t  Remake() const { return fRemakeGeometry;}
00149    void    SetRemake(Bool_t remake=kTRUE) {fRemakeGeometry=remake;}
00150    virtual Int_t  Skip(Int_t Nskip=1);                        // *MENU*
00151    virtual TDataSet *Work();
00152    virtual void   Mark(TVolume *topvol);
00153    virtual void   Call(const Char_t *name); // *MENU 
00154    virtual TRotMatrix *GetMatrix(float theta1, float phi1,
00155                                  float theta2, float phi2,
00156                                  float theta3, float phi3);
00157 
00158 
00159            Int_t SetInputFile(const char* file);
00160 
00161    TDataSet* GetVolume() { return fVolume; }
00162    const TFileSet *GetGeoDirectory() const { return fGeoDirectory;}
00163    TGeoVolume* GetTopGeoVolume() {return fTopGeoVolume;}
00164    static void RootMapTable(Char_t *Cdest,Char_t *Table, Char_t* Spec, 
00165                             Int_t &k, Char_t *iq);
00166    virtual void     Geometry();
00167    virtual Int_t    Agstroot();
00168    virtual Int_t    AgstHits();
00169    virtual Int_t    G2t_volume_id(const Char_t *name, Int_t *numbv);
00170    virtual Int_t    Agvolume(TVolume *&node,Float_t *&par,Float_t *&pos,Float_t *&mot,
00171                              Int_t &who, Int_t &copy,Float_t *&par1,Int_t &npar, char mat[21]);
00172    virtual void     Agnzgete (Int_t &ILK, Int_t &IDE,
00173                               Int_t &NPART, Int_t &IRUN,
00174                               Int_t &IEVT, const Char_t *CGNAM,
00175                               Float_t *VERT,Int_t &IWTFL,Float_t &WEIGH);
00176    
00177    virtual void     Gfxzrm(Int_t & Nlevel, 
00178                      Float_t &x, Float_t &y, Float_t &z,
00179                      Float_t &Theta1, Float_t & Phi1,
00180                      Float_t &Theta2, Float_t & Phi2,
00181                      Float_t &Theta3, Float_t & Phi3,
00182                      Float_t &Type);  
00183    virtual void     Dzddiv(Int_t& idiv ,Int_t &Ldummy,
00184                            const Char_t* path,const Char_t* opt,
00185                            Int_t& one,Int_t &two,Int_t &three,Int_t& iw);
00186 #ifdef DetectorIndex
00187   void        DetSetIndex();
00188   void        DumpIndex(const Char_t *name, const Char_t *vers, const Char_t *fmt, TArrayI &NVmax, TArrayI &Ids);
00189 #endif
00190   TString GetVolumeSrcFile(const char *volumeName) const;
00191   Int_t   KinematicsFromMuDst(Int_t flag=0);
00192   Int_t   SetDatimeFromMuDst();
00193   static St_geant_Maker *instance() {return fgGeantMk;}
00194  protected:
00195    virtual TDataSet  *FindDataSet (const char* logInput,
00196                                     const StMaker *uppMk=0,
00197                                     const StMaker *dowMk=0) const ;
00198    static TDataSet   *fgGeom; 
00199    static TGiant3    *geant3; 
00200    static St_geant_Maker *fgGeantMk; 
00201    TString           mInitialization; // !
00202    TString           mFieldOpt; // !
00203    virtual const char *GetCVS() const
00204    {static const char cvs[]="Tag $Name:  $ $Id: St_geant_Maker.h,v 1.51 2011/09/11 20:57:14 fisyak Exp $ built "__DATE__" "__TIME__ ; return cvs;}
00205 ClassDef(St_geant_Maker,0)   //StAF chain virtual base class for Makers
00206 };
00207 
00208 #endif
00209
St_geant_Maker/St_geant_Maker.h
