[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);
  //
  // 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	//
02d6eccc	120	// hack to avoid using AliGeomManager
cf457606	121	Int_t LayerToVolUID(Int_t lay,int detInLay) const {return GetModuleIndex(lay,detInLay);}
	122	static Int_t ModuleVolUID(Int_t mod) {return mod;}
	123	//
	124	protected:
	125	void FetchMatrices();
	126	TGeoHMatrix* ExtractMatrixT2L(Int_t index) const;
	127	TGeoHMatrix* ExtractMatrixSens(Int_t index) const;
	128	Bool_t GetLayer(Int_t index,Int_t &lay,Int_t &index2) const;
	129	TGeoPNEntry* GetPNEntry(Int_t index) const;
	130	Int_t ExtractNumberOfDetectors(Int_t lay) const;
	131	Int_t ExtractNumberOfLadders(Int_t lay) const;
	132	Int_t ExtractLayerDetType(Int_t lay) const;
	133	Int_t ExtractNumberOfLayers() const;
	134	void BuildITS();
02d6eccc	135	//
451f5018	136	private:
451f5018	137	//
451f5018	138	//
	139	Int_t fVersion; // ITS Version
	140	Int_t fNLayers; // number of layers
	141	Int_t fNModules; //[fNLayers] The total number of modules
	142	Int_t *fNLadders; //[fNLayers] Array of the number of ladders/layer(layer)
	143	Int_t *fLrDetType; //[fNLayers] Array of layer detector types
	144	Int_t *fNDetectors; //[fNLayers] Array of the number of detector/ladder(layer)
	145	Int_t *fLastModIndex; //[fNLayers] max ID of the detctor in the layer
	146	//
cf457606	147	TObjArray* fMatSens; // Sensor's matrices pointers in the geometry
	148	TObjArray* fMatT2L; // Tracking to Local matrices pointers in the geometry
	149	//
451f5018	150	static const char* fgkITSVolName; // ITS mother volume name
	151	static const char* fgkITSLrName; // ITS Layer name
	152	static const char* fgkITSLadName; // ITS Ladder name
	153	static const char* fgkITSModName; // ITS Module name
	154	static const char* fgkITSSensName; // ITS Sensor name
	155	static const char* fgkITSDetTypeName[kNDetTypes]; // ITS upg detType Names
	156	//
	157	static TString fgITSsegmFileName; // file name for segmentations
	158	//
	159	ClassDef(AliITSUGeomTGeo, 1) // ITS geometry based on TGeo
	160	};
	161
	162	//_____________________________________________________________________________________________
	163	inline const char *AliITSUGeomTGeo::GetSymName(Int_t lay,Int_t lad,Int_t det) const
	164	{
	165	// sym name
	166	return GetSymName(GetModuleIndex(lay,lad,det));
	167	}
	168
	169	//_____________________________________________________________________________________________
	170	inline TGeoHMatrix* AliITSUGeomTGeo::GetMatrix(Int_t lay,Int_t lad,Int_t det) const
	171	{
	172	// module current matrix
	173	return GetMatrix(GetModuleIndex(lay,lad,det));
	174	}
	175
	176	//_____________________________________________________________________________________________
	177	inline Bool_t AliITSUGeomTGeo::GetTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const
	178	{
	179	// translation
	180	return GetTranslation(GetModuleIndex(lay,lad,det),t);
	181	}
	182
	183	//_____________________________________________________________________________________________
	184	inline Bool_t AliITSUGeomTGeo::GetRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const
	185	{
	186	// rot
	187	return GetRotation(GetModuleIndex(lay,lad,det),r);
	188	}
	189
	190	//_____________________________________________________________________________________________
	191	inline Bool_t AliITSUGeomTGeo::GetOrigMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m) const
	192	{
	193	// orig matrix
	194	return GetOrigMatrix(GetModuleIndex(lay,lad,det),m);
	195	}
	196
	197	//_____________________________________________________________________________________________
	198	inline Bool_t AliITSUGeomTGeo::GetOrigTranslation(Int_t lay,Int_t lad,Int_t det, Double_t t[3]) const
	199	{
	200	// orig trans
	201	return GetOrigTranslation(GetModuleIndex(lay,lad,det),t);
	202	}
	203
	204	//_____________________________________________________________________________________________
	205	inline Bool_t AliITSUGeomTGeo::GetOrigRotation(Int_t lay,Int_t lad,Int_t det, Double_t r[9]) const
	206	{
	207	// orig rot
	208	return GetOrigRotation(GetModuleIndex(lay,lad,det),r);
	209	}
	210
	211	//_____________________________________________________________________________________________
cf457606	212	inline Bool_t AliITSUGeomTGeo::GetTrackingMatrix(Int_t lay,Int_t lad,Int_t det, TGeoHMatrix &m)
451f5018	213	{
cf457606	214	// tracking mat
cf457606	215	return GetTrackingMatrix(GetModuleIndex(lay,lad,det),m);
451f5018	216	}
	217
	218	//_____________________________________________________________________________________________
cf457606	219	inline Int_t AliITSUGeomTGeo::GetLayerDetTypeID(Int_t lr) const
451f5018	220	{
cf457606	221	// detector type ID of layer
cf457606	222	return fLrDetType[lr];
451f5018	223	}
	224
	225	//_____________________________________________________________________________________________
cf457606	226	inline Int_t AliITSUGeomTGeo::GetModuleDetTypeID(Int_t id) const
451f5018	227	{
cf457606	228	// detector type ID of module
	229	return GetLayerDetTypeID(GetLayer(id));
	230	}
451f5018	231
451f5018	232	//_____________________________________________________________________________________________
cf457606	233	inline const TGeoHMatrix* AliITSUGeomTGeo::GetMatrixSens(Int_t index)
451f5018	234	{
cf457606	235	// access global to sensor matrix
	236	if (!fMatSens) FetchMatrices();
	237	return (TGeoHMatrix*)fMatSens->At(index);
451f5018	238	}
	239
	240	//_____________________________________________________________________________________________
cf457606	241	inline const TGeoHMatrix* AliITSUGeomTGeo::GetMatrixT2L(Int_t index)
451f5018	242	{
cf457606	243	// access tracking to local matrix
	244	if (!fMatT2L) FetchMatrices();
	245	return (TGeoHMatrix*)fMatT2L->At(index);
	246	}
	247
	248	//______________________________________________________________________
	249	inline void AliITSUGeomTGeo::LocalToGlobal(Int_t index,const Double_t loc, Double_t glob)
	250	{
	251	// sensor local to global
	252	GetMatrixSens(index)->LocalToMaster(loc,glob);
	253	}
	254
	255	//______________________________________________________________________
	256	inline void AliITSUGeomTGeo::GlobalToLocal(Int_t index, const Double_t glob, Double_t loc)
	257	{
	258	// global to sensor local
	259	GetMatrixSens(index)->MasterToLocal(glob,loc);
	260	}
	261
	262	//______________________________________________________________________
	263	inline void AliITSUGeomTGeo::LocalToGlobalVect(Int_t index, const Double_t loc, Double_t glob)
	264	{
	265	// sensor local to global
	266	GetMatrixSens(index)->LocalToMasterVect(loc,glob);
	267	}
	268
	269	//______________________________________________________________________
	270	inline void AliITSUGeomTGeo::GlobalToLocalVect(Int_t index, const Double_t glob, Double_t loc)
	271	{
	272	// global to sensor local
	273	GetMatrixSens(index)->MasterToLocalVect(glob,loc);
451f5018	274	}
	275
	276	//_____________________________________________________________________________________________
cf457606	277	inline void AliITSUGeomTGeo::LocalToGlobal(Int_t lay, Int_t lad, Int_t det,const Double_t loc, Double_t glob)
451f5018	278	{
cf457606	279	// Local2Master (sensor)
cf457606	280	LocalToGlobal(GetModuleIndex(lay,lad,det), loc, glob);
451f5018	281	}
	282
	283	//_____________________________________________________________________________________________
cf457606	284	inline void AliITSUGeomTGeo::GlobalToLocal(Int_t lay, Int_t lad, Int_t det,const Double_t glob, Double_t loc)
451f5018	285	{
cf457606	286	// master2local (sensor)
	287	GlobalToLocal(GetModuleIndex(lay,lad,det), glob, loc);
	288	}
451f5018	289
75875328	290	//_____________________________________________________________________________________________
	291	inline const char* AliITSUGeomTGeo::GetDetTypeName(Int_t i)
	292	{
	293	if (i>=kNDetTypes) i/=kMaxSegmPerDetType; // full type is provided
	294	return fgkITSDetTypeName[i];
	295	}
451f5018	296
451f5018	297	#endif