AgataGeometryTransformer.h

Go to the documentation of this file.

#ifndef Gw_AgataGeometryTransformer
#define Gw_AgataGeometryTransformer

#include "GwLogMessage.h"

#ifndef ROOT_TArrayF
#include "TArrayF.h"
#endif

#ifndef ROOT_TObject
#include "TObject.h"
#endif

#ifndef ROOT_TString
#include "TString.h"
#endif

class TGeoCombiTrans;
class TObjArray;
class TGeoVolume;
class TGeoMaterial;

namespace Gw {
   
class AgataGeometryTransformer : public TObject
{
          
public:
   
   AgataGeometryTransformer();
   virtual  ~AgataGeometryTransformer();
   
   
   Bool_t  ReadTransformations(const TString fileName);
   
   TGeoCombiTrans* FillTransMatrix(Double_t* trans, Double_t* rot);
   
   // Transformation methods 
   
   void Global2Local(Int_t detID,
                                         Double_t xg, Double_t yg, Double_t zg, 
                                         Double_t& xl, Double_t& yl, Double_t& zl) const;
   
   void Global2Local(Int_t detID,
                                         Float_t xg, Float_t yg, Float_t zg, 
                                         Float_t& xl, Float_t& yl, Float_t& zl) const;
   
   void Local2Global(Int_t detID,
                                         Double_t xl, Double_t yl, Double_t zl, 
                                         Double_t& xg, Double_t& yg, Double_t& zg) const;
   
   void Local2Global(Int_t detID,
                                         Float_t xl, Float_t yl, Float_t zl, 
                                         Float_t& xg, Float_t& yg, Float_t& zg) const;
   
   
   Int_t GetDepthId(Double_t z);
   
   Int_t GetSectorId(Double_t x, Double_t y);
   
public:
   

   static int BuildAgataCrystals(const char *filename = "asolid", const char *basecrysname = "gePoly");
   

   static int BuildAgataClusters(const char *filename = "aclust", const char *baseclusname = "agataCluster", const char *basecrysname = "gePoly");
   

   static TGeoVolume *BuildAgataDetector(const char *filename = "aeuler", const char *baseclusname = "agataCluster", const char *agataName = "Agata");
   

   static TGeoVolume *ImportAgata(const char *asolid = "asolid", const char *aclust = "aclust", const char *aeuler = "aeuler", 
                                                                  const char *baseclusname = "agataCluster", const char* agataName = "Agata");

   static TGeoVolume* AddDanteModule(const Float_t dz= 0., const Float_t rotz1 = 0., const Float_t roty = 0., const Float_t rotz2 = 0, 
                                                                         const Float_t rotx = 0., const char* basemoduleName = "Module", const char *danteName = "Dante");
   static TGeoVolume* AddBeamPipe(const Float_t dZ = 18, const Float_t minR = 2., const Float_t maxR = 2.3, const char *pipeName = "Pipe");
   
   static TGeoVolume* AddTarget(const Float_t dx = 0.8, const Float_t dy = 0.8, const Float_t dz = 0.05, 
                                                                const char *targetName = "Target", TGeoMaterial* mat = 0x0);
   static TGeoVolume* AddFloor(Float_t dx = 20, const Float_t dy = 0.1, const Float_t dz = 25, const char *floorName = "Floor");
   
   static TGeoVolume* AddPrismaMCP(const char *prismaName = "PrismaMCP", Float_t dx = 10., Float_t dy = 10., Float_t dz = 1., 
                                                                   Float_t phi = 135., TGeoMaterial* mat = 0x0);
   static TGeoVolume* AddPrismaQpole(const char *prismaName = "PrismaQpole", Float_t dz = 21., Float_t rmin = 15., Float_t rmax = 15.7);
   
   static TGeoVolume* AddPrismaDipoleOld(const char *prismaName = "PrismaDipole", Float_t dz = 10., Float_t rmin = 70., 
                                                                                 Float_t rmax = 170., Float_t phi1 = -20, Float_t phi2 = 65);
   static TGeoVolume* AddPrismaDipole(const char *prismaName = "PrismaDipole", Float_t dx = 50., Float_t dy = 10.,Float_t dz = 1.);
   
   static TGeoVolume* AddPrismaMWPC(const char *prismaName = "PrismaMWPC", Float_t dzFcIc = 72., Float_t dx = 75., Float_t dy = 10., 
                                                                        Float_t dz = 65., TGeoMaterial* mat = 0x0);
   
   static TGeoVolume* BuildAgata(const char* agataPath = "./", const char *asolid = "asolid", const char *aclust = "aclust", const char *aeuler = "aeuler", 
                                                                 const char *baseclusname = "agataCluster", const char* agataName = "Agata");
   
   static TGeoVolume* BuildDante(const char* danteConf = "Dante.conf", const char* basemoduleName = "Module", const char *danteName = "Dante");
   
   static TGeoVolume* BuildPrisma(const char* prismaConf = "solver.conf", Float_t dzMCP = 25., Float_t dzQpole = 54., Float_t dzDipole = 160., Float_t dzDipMWPC = 328.5, const char* nameMCP = "PrismaMCP", 
                                                                  const char* nameQPole = "PrismaQpole", const char* nameDipole = "PrismaDipole", const char* nameMWPC  = "PrismaMWPC",  const char* prismaName = "Prisma");
   static void SetDefaultTransparency(Char_t transparency) { fgDefaultTransp = transparency; }

//private:
   
   static Bool_t ReadPrismaConf(Float_t& dzMCP, Float_t& dzQpole, Float_t& dzDipole, Float_t& dzDipMWPC, const char* prismaConf = "solver.conf");
   
private:
   
   TObjArray* fMatrixList;   
   
   mutable LogMessage fLog;  
   
private: 
   
   static const Int_t   fgkMaxDetectorID; 
   static const Int_t   fgkMaxDepthID;    
   static const Float_t fgkDepth[];       
   static const Char_t* fgkGeomName;      
   static const Char_t* fgkGeomTitle;     
   static       Char_t  fgDefaultTransp; // default transparency value

};
} // end namespace

#endif