[u/mrichter/AliRoot.git] / ITS / UPGRADE / AliITSUGeomTGeo.h

#ifndef ALIITSUGEOMTGEO_H
#define ALIITSUGEOMTGEO_H

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/////////////////////////////////////////////////////////////////////////
//  AliITSUGeomTGeo is a simple interface class to TGeoManager       //
//  It is used in the simulation and reconstruction in order to        //
//  query the TGeo ITS geometry                                        //
//                                                                     //
//  author - cvetan.cheshkov@cern.ch                                   //
//  15/02/2007                                                         //
//  adapted to ITSupg 18/07/2012 - ruben.shahoyan@cern.ch              //
//  RS: in order to preserve the static character of the class but     //
//  make it dynamically access geometry, we need to check in every     //
//  method if the structures are initialized. To be converted to       //
//  singleton at later stage.                                          //
//                                                                     //
//  Note on the upgrade detector types:                                //
//  The coarse type defines detectors served by different classes,     //
//  like Pix. Each such a detector type can have kMaxSegmPerDetType    //
//  segmentations (pitch etc.) whose parameteres are stored in the     //
//  AliITSsegmentation derived class (like AliITSUSegmentationPix)     //
//  This allows to have in the setup modules served by the same        //
//  classes but with different segmentations.                          //
//  The full detector type is composed as:                             //
//  CoarseType*kMaxSegmPerDetType + segmentationType                   //
//  The only requirement on the segmentationType that should be        //
//  < kMaxSegmPerDetType.                                              //
//  The methods like GetLayerDetTypeID return the full detector type   //
//                                                                     //
//                                                                     //
/////////////////////////////////////////////////////////////////////////

#include <TObject.h>
#include <TGeoMatrix.h>
#include <TString.h>
#include <TObjArray.h>

class TGeoPNEntry;
class TDatime;

class AliITSUGeomTGeo : public TObject {

 public:
  enum {kITSVNA, kITSVUpg}; // ITS version
  enum {kDetTypePix=0, kNDetTypes, kMaxSegmPerDetType=10}; // defined detector types (each one can have different segmentations)
  //
  AliITSUGeomTGeo(Bool_t build = kFALSE);
  virtual ~AliITSUGeomTGeo(); 
  AliITSUGeomTGeo(const AliITSUGeomTGeo &src);
  AliITSUGeomTGeo& operator=(const AliITSUGeomTGeo &geom);
  //
  Int_t  GetNModules()                                                    const {return fNModules;}
  Int_t  GetNDetectors(Int_t lay)                                         const {return fNDetectors[lay];}
  Int_t  GetNLadders(Int_t lay)                                           const {return fNLadders[lay];}
  Int_t  GetNLayers()                                                     const {return fNLayers;}
  
  Int_t  GetModuleIndex(Int_t lay,int detInLay)                           const {return GetFirstModIndex(lay)+detInLay;}
  Int_t  GetModuleIndex(Int_t lay,Int_t lad,Int_t det)                    const;
  Bool_t GetModuleId(Int_t index,Int_t &lay,Int_t &lad,Int_t &det)        const;
  Int_t  GetLayer(Int_t index)                                            const;
  Int_t  GetLadder(Int_t index)                                           const;
  Int_t  GetModIdInLayer(Int_t index)                                     const;
  Int_t  GetModIdInLadder(Int_t index)                                    const;
  //
  Int_t  GetLastModIndex(Int_t lay)                                       const {return fLastModIndex[lay];}
  Int_t  GetFirstModIndex(Int_t lay)                                      const {return (lay==0) ? 0:fLastModIndex[lay-1]+1;}
  //  
  const char *GetSymName(Int_t index)                                     const;
  const char *GetSymName(Int_t lay,Int_t lad,Int_t det)                   const;
  //
  // Attention: these are the matrices for the alignable volumes of the modules, i.e. not necessarily the sensors
  TGeoHMatrix* GetMatrix(Int_t index)                                     const;
  TGeoHMatrix* GetMatrix(Int_t lay,Int_t lad,Int_t det)                   const;
  Bool_t GetTranslation(Int_t index, Double_t t[3])                       const;
  Bool_t GetTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3])     const;
  Bool_t GetRotation(Int_t index, Double_t r[9])                          const;
  Bool_t GetRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9])        const;
  Bool_t GetOrigMatrix(Int_t index, TGeoHMatrix &m)                       const;
  Bool_t GetOrigMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m)     const;
  Bool_t GetOrigTranslation(Int_t index, Double_t t[3])                   const;
  Bool_t GetOrigTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const;
  Bool_t GetOrigRotation(Int_t index, Double_t r[9])                      const;
  Bool_t GetOrigRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9])    const;
  //
  const TGeoHMatrix* GetMatrixT2L(Int_t index);
  const TGeoHMatrix* GetMatrixT2L(Int_t lay,Int_t lad,Int_t det)  {return GetMatrixT2L( GetModuleIndex(lay,lad,det) );}
  const TGeoHMatrix* GetMatrixSens(Int_t index);
  const TGeoHMatrix* GetMatrixSens(Int_t lay,Int_t lad,Int_t det) {return GetMatrixSens( GetModuleIndex(lay,lad,det) );}
  //
  Bool_t GetTrackingMatrix(Int_t index, TGeoHMatrix &m);
  Bool_t GetTrackingMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m);
  //
  // Attention: these are transformations wrt sensitive volume!
  void   LocalToGlobal(Int_t index, const Double_t *loc, Double_t *glob);
  void   LocalToGlobal(Int_t lay, Int_t lad, Int_t det,const Double_t *loc, Double_t *glob);
  //
  void   GlobalToLocal(Int_t index, const Double_t *glob, Double_t *loc);
  void   GlobalToLocal(Int_t lay, Int_t lad, Int_t det,const Double_t *glob, Double_t *loc);
  //
  void   LocalToGlobalVect(Int_t index, const Double_t *loc, Double_t *glob);
  void   GlobalToLocalVect(Int_t index, const Double_t *glob, Double_t *loc);
  Int_t  GetLayerDetTypeID(Int_t lr)                                         const;
  Int_t  GetModuleDetTypeID(Int_t id)                                        const;
  //
  virtual void Print(Option_t *opt="")  const;

  static const char* GetITSVolPattern()                                 {return fgkITSVolName;}
  static const char* GetITSLayerPattern()                               {return fgkITSLrName;}
  static const char* GetITSLadderPattern()                              {return fgkITSLadName;}
  static const char* GetITSModulePattern()                              {return fgkITSModName;}
  static const char* GetITSSensorPattern()                              {return fgkITSSensName;}
  static const char* GetITSsegmentationFileName()                       {return fgITSsegmFileName.Data();}
  static const char* GetDetTypeName(Int_t i);
  static void        SetITSsegmentationFileName(const char* nm)         {fgITSsegmFileName = nm;}
  static UInt_t      ComposeDetTypeID(UInt_t segmId);
  //
  static const char *ComposeSymNameITS();
  static const char *ComposeSymNameLayer(Int_t lr);
  static const char *ComposeSymNameLadder(Int_t lr, Int_t lad);
  static const char *ComposeSymNameModule(Int_t lr, Int_t lad, int det);
  //
  // hack to avoid using AliGeomManager
  Int_t              LayerToVolUID(Int_t lay,int detInLay)        const {return GetModuleIndex(lay,detInLay);}
  static Int_t       ModuleVolUID(Int_t mod)                            {return mod;}
  //
 protected:
  void         FetchMatrices();
  TGeoHMatrix* ExtractMatrixT2L(Int_t index)                      const;
  TGeoHMatrix* ExtractMatrixSens(Int_t index)                     const;
  Bool_t       GetLayer(Int_t index,Int_t &lay,Int_t &index2)     const;
  TGeoPNEntry* GetPNEntry(Int_t index)                            const;
  Int_t        ExtractNumberOfDetectors(Int_t lay)                const;
  Int_t        ExtractNumberOfLadders(Int_t lay)                  const;
  Int_t        ExtractLayerDetType(Int_t lay)                     const;
  Int_t        ExtractNumberOfLayers()                            const;
  void         BuildITS();
  //
 private:
//
  //
  Int_t  fVersion;             // ITS Version 
  Int_t  fNLayers;             // number of layers
  Int_t  fNModules;            //[fNLayers] The total number of modules
  Int_t *fNLadders;            //[fNLayers] Array of the number of ladders/layer(layer)
  Int_t *fLrDetType;           //[fNLayers] Array of layer detector types
  Int_t *fNDetectors;          //[fNLayers] Array of the number of detector/ladder(layer)
  Int_t *fLastModIndex;        //[fNLayers] max ID of the detctor in the layer
  //
  TObjArray* fMatSens;         // Sensor's matrices pointers in the geometry
  TObjArray* fMatT2L;          // Tracking to Local matrices pointers in the geometry
  //
  static const char*  fgkITSVolName;             // ITS mother volume name
  static const char*  fgkITSLrName;              // ITS Layer name
  static const char*  fgkITSLadName;             // ITS Ladder name 
  static const char*  fgkITSModName;             // ITS Module name 
  static const char*  fgkITSSensName;            // ITS Sensor name 
  static const char*  fgkITSDetTypeName[kNDetTypes]; // ITS upg detType Names
  //
  static TString      fgITSsegmFileName;         // file name for segmentations
  //
  ClassDef(AliITSUGeomTGeo, 1) // ITS geometry based on TGeo
};

//_____________________________________________________________________________________________
inline const char *AliITSUGeomTGeo::GetSymName(Int_t lay,Int_t lad,Int_t det) const    
{
  // sym name
  return GetSymName(GetModuleIndex(lay,lad,det));
}

//_____________________________________________________________________________________________
inline TGeoHMatrix* AliITSUGeomTGeo::GetMatrix(Int_t lay,Int_t lad,Int_t det) const    
{
  // module current matrix
  return GetMatrix(GetModuleIndex(lay,lad,det));
}

//_____________________________________________________________________________________________
inline Bool_t AliITSUGeomTGeo::GetTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const    
{
  // translation
  return GetTranslation(GetModuleIndex(lay,lad,det),t); 
}

//_____________________________________________________________________________________________
inline Bool_t AliITSUGeomTGeo::GetRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const    
{
  // rot
  return GetRotation(GetModuleIndex(lay,lad,det),r); 
}

//_____________________________________________________________________________________________
inline Bool_t AliITSUGeomTGeo::GetOrigMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m) const    
{
  // orig matrix
  return GetOrigMatrix(GetModuleIndex(lay,lad,det),m); 
}

//_____________________________________________________________________________________________
inline Bool_t AliITSUGeomTGeo::GetOrigTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const    
{
  // orig trans
  return GetOrigTranslation(GetModuleIndex(lay,lad,det),t); 
}

//_____________________________________________________________________________________________
inline Bool_t AliITSUGeomTGeo::GetOrigRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const    
{
  // orig rot
  return GetOrigRotation(GetModuleIndex(lay,lad,det),r); 
}

//_____________________________________________________________________________________________
inline Bool_t AliITSUGeomTGeo::GetTrackingMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m)
{
  // tracking mat
  return GetTrackingMatrix(GetModuleIndex(lay,lad,det),m); 
}

//_____________________________________________________________________________________________
inline Int_t  AliITSUGeomTGeo::GetLayerDetTypeID(Int_t lr) const  
{
  // detector type ID of layer
  return fLrDetType[lr];
}

//_____________________________________________________________________________________________
inline Int_t  AliITSUGeomTGeo::GetModuleDetTypeID(Int_t id) const  
{
  // detector type ID of module
  return GetLayerDetTypeID(GetLayer(id));
} 

//_____________________________________________________________________________________________
inline const TGeoHMatrix* AliITSUGeomTGeo::GetMatrixSens(Int_t index)
{
  // access global to sensor matrix
  if (!fMatSens) FetchMatrices();
  return (TGeoHMatrix*)fMatSens->At(index);
}

//_____________________________________________________________________________________________
inline const TGeoHMatrix* AliITSUGeomTGeo::GetMatrixT2L(Int_t index)
{
  // access tracking to local matrix
  if (!fMatT2L) FetchMatrices();
  return (TGeoHMatrix*)fMatT2L->At(index);
}

//______________________________________________________________________
inline void AliITSUGeomTGeo::LocalToGlobal(Int_t index,const Double_t *loc, Double_t *glob)
{
  // sensor local to global 
  GetMatrixSens(index)->LocalToMaster(loc,glob);
}

//______________________________________________________________________
inline void AliITSUGeomTGeo::GlobalToLocal(Int_t index, const Double_t *glob, Double_t *loc)
{
  // global to sensor local 
  GetMatrixSens(index)->MasterToLocal(glob,loc);
}

//______________________________________________________________________
inline void AliITSUGeomTGeo::LocalToGlobalVect(Int_t index, const Double_t *loc, Double_t *glob)
{
  // sensor local to global 
  GetMatrixSens(index)->LocalToMasterVect(loc,glob);
}

//______________________________________________________________________
inline void AliITSUGeomTGeo::GlobalToLocalVect(Int_t index, const Double_t *glob, Double_t *loc)
{
  // global to sensor local
  GetMatrixSens(index)->MasterToLocalVect(glob,loc);
}

//_____________________________________________________________________________________________
inline void AliITSUGeomTGeo::LocalToGlobal(Int_t lay, Int_t lad, Int_t det,const Double_t *loc, Double_t *glob)
{
  // Local2Master (sensor)
  LocalToGlobal(GetModuleIndex(lay,lad,det), loc, glob);
}

//_____________________________________________________________________________________________
inline void AliITSUGeomTGeo::GlobalToLocal(Int_t lay, Int_t lad, Int_t det,const Double_t *glob, Double_t *loc)
{
  // master2local (sensor)
  GlobalToLocal(GetModuleIndex(lay,lad,det), glob, loc);
}

//_____________________________________________________________________________________________
inline const char* AliITSUGeomTGeo::GetDetTypeName(Int_t i)
{
  if (i>=kNDetTypes) i/=kMaxSegmPerDetType; // full type is provided
  return fgkITSDetTypeName[i];
}

#endif
Commit	Line	Data
451f5018	1	#ifndef ALIITSUGEOMTGEO_H
451f5018	2	#define ALIITSUGEOMTGEO_H
a11ef2e4	3
451f5018	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	/////////////////////////////////////////////////////////////////////////
	8	// AliITSUGeomTGeo is a simple interface class to TGeoManager //
	9	// It is used in the simulation and reconstruction in order to //
	10	// query the TGeo ITS geometry //
	11	// //
	12	// author - cvetan.cheshkov@cern.ch //
	13	// 15/02/2007 //
	14	// adapted to ITSupg 18/07/2012 - ruben.shahoyan@cern.ch //
	15	// RS: in order to preserve the static character of the class but //
	16	// make it dynamically access geometry, we need to check in every //
	17	// method if the structures are initialized. To be converted to //
	18	// singleton at later stage. //
	19	// //
	20	// Note on the upgrade detector types: //
	21	// The coarse type defines detectors served by different classes, //
75875328	22	// like Pix. Each such a detector type can have kMaxSegmPerDetType //
451f5018	23	// segmentations (pitch etc.) whose parameteres are stored in the //
75875328	24	// AliITSsegmentation derived class (like AliITSUSegmentationPix) //
451f5018	25	// This allows to have in the setup modules served by the same //
	26	// classes but with different segmentations. //
	27	// The full detector type is composed as: //
	28	// CoarseType*kMaxSegmPerDetType + segmentationType //
	29	// The only requirement on the segmentationType that should be //
	30	// < kMaxSegmPerDetType. //
	31	// The methods like GetLayerDetTypeID return the full detector type //
	32	// //
	33	// //
	34	/////////////////////////////////////////////////////////////////////////
	35
	36	#include <TObject.h>
	37	#include <TGeoMatrix.h>
	38	#include <TString.h>
a11ef2e4	39	#include <TObjArray.h>
451f5018	40
	41	class TGeoPNEntry;
	42	class TDatime;
	43
	44	class AliITSUGeomTGeo : public TObject {
	45
	46	public:
	47	enum {kITSVNA, kITSVUpg}; // ITS version
02d6eccc	48	enum {kDetTypePix=0, kNDetTypes, kMaxSegmPerDetType=10}; // defined detector types (each one can have different segmentations)
451f5018	49	//
	50	AliITSUGeomTGeo(Bool_t build = kFALSE);
	51	virtual ~AliITSUGeomTGeo();
	52	AliITSUGeomTGeo(const AliITSUGeomTGeo &src);
	53	AliITSUGeomTGeo& operator=(const AliITSUGeomTGeo &geom);
	54	//
02d6eccc	55	Int_t GetNModules() const {return fNModules;}
	56	Int_t GetNDetectors(Int_t lay) const {return fNDetectors[lay];}
	57	Int_t GetNLadders(Int_t lay) const {return fNLadders[lay];}
	58	Int_t GetNLayers() const {return fNLayers;}
451f5018	59
02d6eccc	60	Int_t GetModuleIndex(Int_t lay,int detInLay) const {return GetFirstModIndex(lay)+detInLay;}
451f5018	61	Int_t GetModuleIndex(Int_t lay,Int_t lad,Int_t det) const;
	62	Bool_t GetModuleId(Int_t index,Int_t &lay,Int_t &lad,Int_t &det) const;
	63	Int_t GetLayer(Int_t index) const;
	64	Int_t GetLadder(Int_t index) const;
	65	Int_t GetModIdInLayer(Int_t index) const;
	66	Int_t GetModIdInLadder(Int_t index) const;
	67	//
	68	Int_t GetLastModIndex(Int_t lay) const {return fLastModIndex[lay];}
	69	Int_t GetFirstModIndex(Int_t lay) const {return (lay==0) ? 0:fLastModIndex[lay-1]+1;}
	70	//
	71	const char *GetSymName(Int_t index) const;
	72	const char *GetSymName(Int_t lay,Int_t lad,Int_t det) const;
	73	//
	74	// Attention: these are the matrices for the alignable volumes of the modules, i.e. not necessarily the sensors
	75	TGeoHMatrix* GetMatrix(Int_t index) const;
	76	TGeoHMatrix* GetMatrix(Int_t lay,Int_t lad,Int_t det) const;
	77	Bool_t GetTranslation(Int_t index, Double_t t[3]) const;
	78	Bool_t GetTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const;
	79	Bool_t GetRotation(Int_t index, Double_t r[9]) const;
	80	Bool_t GetRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const;
	81	Bool_t GetOrigMatrix(Int_t index, TGeoHMatrix &m) const;
	82	Bool_t GetOrigMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m) const;
	83	Bool_t GetOrigTranslation(Int_t index, Double_t t[3]) const;
	84	Bool_t GetOrigTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const;
	85	Bool_t GetOrigRotation(Int_t index, Double_t r[9]) const;
	86	Bool_t GetOrigRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const;
	87	//
cf457606	88	const TGeoHMatrix* GetMatrixT2L(Int_t index);
	89	const TGeoHMatrix* GetMatrixT2L(Int_t lay,Int_t lad,Int_t det) {return GetMatrixT2L( GetModuleIndex(lay,lad,det) );}
	90	const TGeoHMatrix* GetMatrixSens(Int_t index);
	91	const TGeoHMatrix* GetMatrixSens(Int_t lay,Int_t lad,Int_t det) {return GetMatrixSens( GetModuleIndex(lay,lad,det) );}
451f5018	92	//
cf457606	93	Bool_t GetTrackingMatrix(Int_t index, TGeoHMatrix &m);
cf457606	94	Bool_t GetTrackingMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m);
451f5018	95	//
451f5018	96	// Attention: these are transformations wrt sensitive volume!
cf457606	97	void LocalToGlobal(Int_t index, const Double_t loc, Double_t glob);
cf457606	98	void LocalToGlobal(Int_t lay, Int_t lad, Int_t det,const Double_t loc, Double_t glob);
451f5018	99	//
cf457606	100	void GlobalToLocal(Int_t index, const Double_t glob, Double_t loc);
cf457606	101	void GlobalToLocal(Int_t lay, Int_t lad, Int_t det,const Double_t glob, Double_t loc);
451f5018	102	//
cf457606	103	void LocalToGlobalVect(Int_t index, const Double_t loc, Double_t glob);
cf457606	104	void GlobalToLocalVect(Int_t index, const Double_t glob, Double_t loc);
451f5018	105	Int_t GetLayerDetTypeID(Int_t lr) const;
	106	Int_t GetModuleDetTypeID(Int_t id) const;
	107	//
02d6eccc	108	virtual void Print(Option_t *opt="") const;
02d6eccc	109
cf457606	110	static const char* GetITSVolPattern() {return fgkITSVolName;}
	111	static const char* GetITSLayerPattern() {return fgkITSLrName;}
	112	static const char* GetITSLadderPattern() {return fgkITSLadName;}
	113	static const char* GetITSModulePattern() {return fgkITSModName;}
	114	static const char* GetITSSensorPattern() {return fgkITSSensName;}
cf457606	115	static const char* GetITSsegmentationFileName() {return fgITSsegmFileName.Data();}
75875328	116	static const char* GetDetTypeName(Int_t i);
cf457606	117	static void SetITSsegmentationFileName(const char* nm) {fgITSsegmFileName = nm;}
451f5018	118	static UInt_t ComposeDetTypeID(UInt_t segmId);
451f5018	119	//
392efe73	120	static const char *ComposeSymNameITS();
	121	static const char *ComposeSymNameLayer(Int_t lr);
	122	static const char *ComposeSymNameLadder(Int_t lr, Int_t lad);
	123	static const char *ComposeSymNameModule(Int_t lr, Int_t lad, int det);
	124	//
02d6eccc	125	// hack to avoid using AliGeomManager
cf457606	126	Int_t LayerToVolUID(Int_t lay,int detInLay) const {return GetModuleIndex(lay,detInLay);}
	127	static Int_t ModuleVolUID(Int_t mod) {return mod;}
	128	//
	129	protected:
	130	void FetchMatrices();
	131	TGeoHMatrix* ExtractMatrixT2L(Int_t index) const;
	132	TGeoHMatrix* ExtractMatrixSens(Int_t index) const;
	133	Bool_t GetLayer(Int_t index,Int_t &lay,Int_t &index2) const;
	134	TGeoPNEntry* GetPNEntry(Int_t index) const;
	135	Int_t ExtractNumberOfDetectors(Int_t lay) const;
	136	Int_t ExtractNumberOfLadders(Int_t lay) const;
	137	Int_t ExtractLayerDetType(Int_t lay) const;
	138	Int_t ExtractNumberOfLayers() const;
	139	void BuildITS();
02d6eccc	140	//
451f5018	141	private:
451f5018	142	//
451f5018	143	//
	144	Int_t fVersion; // ITS Version
	145	Int_t fNLayers; // number of layers
	146	Int_t fNModules; //[fNLayers] The total number of modules
	147	Int_t *fNLadders; //[fNLayers] Array of the number of ladders/layer(layer)
	148	Int_t *fLrDetType; //[fNLayers] Array of layer detector types
	149	Int_t *fNDetectors; //[fNLayers] Array of the number of detector/ladder(layer)
	150	Int_t *fLastModIndex; //[fNLayers] max ID of the detctor in the layer
	151	//
cf457606	152	TObjArray* fMatSens; // Sensor's matrices pointers in the geometry
	153	TObjArray* fMatT2L; // Tracking to Local matrices pointers in the geometry
	154	//
451f5018	155	static const char* fgkITSVolName; // ITS mother volume name
	156	static const char* fgkITSLrName; // ITS Layer name
	157	static const char* fgkITSLadName; // ITS Ladder name
	158	static const char* fgkITSModName; // ITS Module name
	159	static const char* fgkITSSensName; // ITS Sensor name
	160	static const char* fgkITSDetTypeName[kNDetTypes]; // ITS upg detType Names
	161	//
	162	static TString fgITSsegmFileName; // file name for segmentations
	163	//
	164	ClassDef(AliITSUGeomTGeo, 1) // ITS geometry based on TGeo
	165	};
	166
	167	//_____________________________________________________________________________________________
	168	inline const char *AliITSUGeomTGeo::GetSymName(Int_t lay,Int_t lad,Int_t det) const
	169	{
	170	// sym name
	171	return GetSymName(GetModuleIndex(lay,lad,det));
	172	}
	173
	174	//_____________________________________________________________________________________________
	175	inline TGeoHMatrix* AliITSUGeomTGeo::GetMatrix(Int_t lay,Int_t lad,Int_t det) const
	176	{
	177	// module current matrix
	178	return GetMatrix(GetModuleIndex(lay,lad,det));
	179	}
	180
	181	//_____________________________________________________________________________________________
	182	inline Bool_t AliITSUGeomTGeo::GetTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const
	183	{
	184	// translation
	185	return GetTranslation(GetModuleIndex(lay,lad,det),t);
	186	}
	187
	188	//_____________________________________________________________________________________________
	189	inline Bool_t AliITSUGeomTGeo::GetRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const
	190	{
	191	// rot
	192	return GetRotation(GetModuleIndex(lay,lad,det),r);
	193	}
	194
	195	//_____________________________________________________________________________________________
	196	inline Bool_t AliITSUGeomTGeo::GetOrigMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m) const
	197	{
	198	// orig matrix
	199	return GetOrigMatrix(GetModuleIndex(lay,lad,det),m);
	200	}
	201
	202	//_____________________________________________________________________________________________
	203	inline Bool_t AliITSUGeomTGeo::GetOrigTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const
	204	{
	205	// orig trans
	206	return GetOrigTranslation(GetModuleIndex(lay,lad,det),t);
	207	}
	208
	209	//_____________________________________________________________________________________________
	210	inline Bool_t AliITSUGeomTGeo::GetOrigRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const
	211	{
	212	// orig rot
	213	return GetOrigRotation(GetModuleIndex(lay,lad,det),r);
	214	}
	215
	216	//_____________________________________________________________________________________________
cf457606	217	inline Bool_t AliITSUGeomTGeo::GetTrackingMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m)
451f5018	218	{
cf457606	219	// tracking mat
cf457606	220	return GetTrackingMatrix(GetModuleIndex(lay,lad,det),m);
451f5018	221	}
	222
	223	//_____________________________________________________________________________________________
cf457606	224	inline Int_t AliITSUGeomTGeo::GetLayerDetTypeID(Int_t lr) const
451f5018	225	{
cf457606	226	// detector type ID of layer
cf457606	227	return fLrDetType[lr];
451f5018	228	}
	229
	230	//_____________________________________________________________________________________________
cf457606	231	inline Int_t AliITSUGeomTGeo::GetModuleDetTypeID(Int_t id) const
451f5018	232	{
cf457606	233	// detector type ID of module
	234	return GetLayerDetTypeID(GetLayer(id));
	235	}
451f5018	236
451f5018	237	//_____________________________________________________________________________________________
cf457606	238	inline const TGeoHMatrix* AliITSUGeomTGeo::GetMatrixSens(Int_t index)
451f5018	239	{
cf457606	240	// access global to sensor matrix
	241	if (!fMatSens) FetchMatrices();
	242	return (TGeoHMatrix*)fMatSens->At(index);
451f5018	243	}
	244
	245	//_____________________________________________________________________________________________
cf457606	246	inline const TGeoHMatrix* AliITSUGeomTGeo::GetMatrixT2L(Int_t index)
451f5018	247	{
cf457606	248	// access tracking to local matrix
	249	if (!fMatT2L) FetchMatrices();
	250	return (TGeoHMatrix*)fMatT2L->At(index);
	251	}
	252
	253	//______________________________________________________________________
	254	inline void AliITSUGeomTGeo::LocalToGlobal(Int_t index,const Double_t loc, Double_t glob)
	255	{
	256	// sensor local to global
	257	GetMatrixSens(index)->LocalToMaster(loc,glob);
	258	}
	259
	260	//______________________________________________________________________
	261	inline void AliITSUGeomTGeo::GlobalToLocal(Int_t index, const Double_t glob, Double_t loc)
	262	{
	263	// global to sensor local
	264	GetMatrixSens(index)->MasterToLocal(glob,loc);
	265	}
	266
	267	//______________________________________________________________________
	268	inline void AliITSUGeomTGeo::LocalToGlobalVect(Int_t index, const Double_t loc, Double_t glob)
	269	{
	270	// sensor local to global
	271	GetMatrixSens(index)->LocalToMasterVect(loc,glob);
	272	}
	273
	274	//______________________________________________________________________
	275	inline void AliITSUGeomTGeo::GlobalToLocalVect(Int_t index, const Double_t glob, Double_t loc)
	276	{
	277	// global to sensor local
	278	GetMatrixSens(index)->MasterToLocalVect(glob,loc);
451f5018	279	}
	280
	281	//_____________________________________________________________________________________________
cf457606	282	inline void AliITSUGeomTGeo::LocalToGlobal(Int_t lay, Int_t lad, Int_t det,const Double_t loc, Double_t glob)
451f5018	283	{
cf457606	284	// Local2Master (sensor)
cf457606	285	LocalToGlobal(GetModuleIndex(lay,lad,det), loc, glob);
451f5018	286	}
	287
	288	//_____________________________________________________________________________________________
cf457606	289	inline void AliITSUGeomTGeo::GlobalToLocal(Int_t lay, Int_t lad, Int_t det,const Double_t glob, Double_t loc)
451f5018	290	{
cf457606	291	// master2local (sensor)
	292	GlobalToLocal(GetModuleIndex(lay,lad,det), glob, loc);
	293	}
451f5018	294
75875328	295	//_____________________________________________________________________________________________
	296	inline const char* AliITSUGeomTGeo::GetDetTypeName(Int_t i)
	297	{
	298	if (i>=kNDetTypes) i/=kMaxSegmPerDetType; // full type is provided
	299	return fgkITSDetTypeName[i];
	300	}
451f5018	301
451f5018	302	#endif