StRoot: StEEmcUtil/StEEmcSmd/EEmcSmdGeom.h Source File

00001 
00032 #ifndef EEMCSMDGEOM_H
00033 #define EEMCSMDGEOM_H
00034 
00035 #include <iostream>
00036 
00037 #include "TObject.h"
00038 #include "TVector3.h"
00039 #include "TString.h"
00040 
00041 #include "StEEmcUtil/EEmcGeom/EEmcGeomDefs.h"
00042 
00043 #include <vector>
00044 #ifndef ST_NO_NAMESPACES
00045 using std::vector;
00046 using std::cout;
00047 #endif
00048 
00049 struct StructEEmcSmdParam {
00050   float zPlane[kEEmcNumSmdPlanes];  // Z of planes
00051   float rOffset[kEEmcNumSmdPlanes]; // radius offsets
00052   float stripWidth;                 // width of strip
00053 };
00054 
00055 struct StructEEmcStripId{
00056   int  stripId;                     // strip Id  1-288 (283 for edge sector)
00057   int  UVId;                        // 1 for U and 2 for V
00058   int  sectorId;                    // sector Id 1-12 
00059   int  planeId;                     // plane Id (depth) 1-3
00060 };
00061 
00062 struct StructEEmcStrip{
00063   StructEEmcStripId stripStructId; // including 4 Ids
00064   TVector3    end1;          // one end of strip         
00065   TVector3    end2;          // the other end of strip
00066   float length;                    // length of strip 
00067 };
00068 
00069 ostream& operator<<(ostream &os, const StructEEmcStrip &strip);
00070 
00071 // Define vector of strip pointers and its Iterator
00072 #ifndef ST_NO_TEMPLATE_DEF_ARGS
00073 typedef vector<int> intVec;  
00074 typedef vector<StructEEmcStrip*> EEmcStripPtrVec;  
00075 #else
00076 typedef vector<int,allocator<int>> intVec; 
00077 typedef vector<StructEEmcStrip*,allocator<StructEEmcStrip*>> EEmcStripPtrVec; 
00078 #endif   
00079 typedef vector<StructEEmcStrip*>::iterator EEmcStripPtrVecIter; 
00080 
00081 
00082 struct StructEEmcSmdSector {
00083   int   sectorId;       
00084   int   planeId;        
00085   float phiMin;    // minimun phi
00086   float phiMax;    // maximum phi
00087   float rMin;      // minimum radius
00088   float rMax;      // maximum radius
00089   EEmcStripPtrVec stripPtrVec;  
00090 };
00091 
00092 class EEmcSmdGeom : public TObject {
00093 public:  
00094 
00095   EEmcSmdGeom();
00096   virtual ~EEmcSmdGeom(); 
00097 
00098 protected:
00099 
00100   StructEEmcSmdParam     mEEmcSmdParam;  
00101   StructEEmcSmdSector    mEEmcSector[kEEmcNumSmdUVs][kEEmcNumSectors]; 
00102 
00103   EEmcStripPtrVec       mStripPtrVector; 
00104   bool                   mIsSectorIn[kEEmcNumSectors];    
00105   int   kEEmcSmdMap_iPlane[kEEmcNumSmdUVs][kEEmcNumSectors];
00106 
00107   void buildSmdGeom();
00108 
00109   static EEmcSmdGeom* sInstance;
00110 
00111 public:
00114   const StructEEmcStrip* getDcaStripPtr(const Int_t iPlane, const TVector3& point, Float_t* dca) const;
00115   const StructEEmcStrip* getDcaStripPtr(const Int_t iPlane, const Int_t iSec, const TVector3& point, Float_t* dca) const;
00116   
00117   
00118   static EEmcSmdGeom* instance();   // handle the only instance
00119   static EEmcSmdGeom* instance(intVec sectorIdVec);   
00120   
00121   void init();         // init the dbase
00122   
00124   void buildStripPtrVector();
00125 
00127   void setSectors(const intVec sectorIdVec); 
00128   
00130   bool IsSectorIn(const Int_t iSec) const {return mIsSectorIn[iSec];}
00131   
00133   StructEEmcStrip initStrip() const;
00134 
00136   StructEEmcSmdParam &getEEmcSmdParam() {return mEEmcSmdParam;}
00137   const StructEEmcSmdParam &getEEmcSmdParam() const {return mEEmcSmdParam;}
00138 
00140   StructEEmcSmdSector &getEEmcSector(const Int_t iUV, const Int_t iSec) {return mEEmcSector[iUV][iSec];}
00141   const StructEEmcSmdSector &getEEmcSector(const Int_t iUV, const Int_t iSec) const {return mEEmcSector[iUV][iSec];}
00142 
00144   Int_t getEEmcISec(const Int_t iPlane, const TVector3& point) const;
00145   
00147   StructEEmcStrip* getStripPtr(const Int_t iStrip, const Int_t iUV, const Int_t iSec);
00148   const StructEEmcStrip* getStripPtr(const Int_t iStrip, const Int_t iUV, const Int_t iSec) const;
00149   
00158   const  StructEEmcStrip* getDca2Strip(const Int_t iUV, const TVector3& point, Float_t* dca) const;
00159 
00160   
00169   TVector3 getIntersection ( Int_t iSec, Int_t iUStrip, Int_t iVStrip, const TVector3 &vertex ) const;
00170 
00178   TVector3 getIntersection ( const StructEEmcStrip *u, const StructEEmcStrip *v, const TVector3 &vertex ) const;
00179 
00181   TVector3 getIntersection ( Int_t iSec, Int_t iUStrip, Int_t iVStrip ) const;
00183   TVector3 getIntersection ( Int_t iSec, Float_t iUStrip, Float_t iVStrip ) const
00184     { return getIntersection( iSec, (Int_t)iUStrip, (Int_t)iVStrip ); }
00186   TVector3 getIntersection ( const StructEEmcStrip *u, const StructEEmcStrip *v ) const;
00187 
00188 
00192   Int_t getNStrips ( Int_t iSec, Int_t iUV ) const { return getEEmcSector(iUV,iSec).stripPtrVec.size();  }
00193 
00195   bool matchStrips(const StructEEmcStripId &stripStructId1, const StructEEmcStripId &stripStructId2, Int_t nTolerance) const;
00196 
00197   // mehtod for C-scripts
00198 
00200   TVector3  getstripEnd(const StructEEmcStrip &strip, const Int_t endId) const;
00201 
00202   //
00203   // methods of printout 
00204   //
00205 
00206   void printGeom(ostream& os = cout) const;
00207   void printSector(const StructEEmcSmdSector Sector, ostream& os = cout) const;
00208   void printStrip(const StructEEmcStrip Strip, ostream& os = cout) const;
00209   void printStripId(const StructEEmcStripId StripId, ostream& os = cout) const;
00210   //void printSectorPhis(const Int_t iPlane, const Int_t iSec,ostream& os = cout);
00211 
00212   ClassDef(EEmcSmdGeom,1)  // STAR Endcap Electromagnetic Calorimeter SMD Geometry Class
00213 };
00214 
00215 #endif
00216  
00217 /*******************************************************************
00218  *
00219  *
00220  * $Log: EEmcSmdGeom.h,v $
00221  * Revision 1.10  2010/08/26 22:48:55  ogrebeny
00222  * Improved constness
00223  *
00224  * Revision 1.9  2007/07/12 19:30:15  fisyak
00225  * Add includes for ROOT 5.16
00226  *
00227  * Revision 1.8  2007/02/01 13:47:39  balewski
00228  * bug fix in getDca2Strip(), more methodhs are public
00229  *
00230  * Revision 1.7  2007/01/26 00:51:08  balewski
00231  * too strong protection
00232  *
00233  * Revision 1.6  2007/01/25 22:33:21  balewski
00234  * add:
00235  * - better writeup
00236  * - new simpler to use method calculating dca fo track to strip, it is just a wrapper, some approximations were used, may fail at the sector boundary
00237  *
00238  * Revision 1.5  2007/01/12 23:53:14  jwebb
00239  * Fix applied to eliminate parralax error in the EEmcSmdGeom::getIntersection()
00240  * method.
00241  *
00242  * Revision 1.3  2004/07/01 15:35:50  jwebb
00243  * Added placeholder method getIntersection(Int_t,Float_t,Float_t).  For now
00244  * it just calls the getIntersection with strips converted to integers.
00245  * Later we will write code to handle fractional strips.
00246  *
00247  * Revision 1.2  2004/02/03 22:57:54  jwebb
00248  * Added StEEmcSmdGeom::instance(), which is sort of needed...
00249  *
00250  * Revision 1.1  2004/01/29 15:26:10  jwebb
00251  * The StEEmcSmdGeom class was split into two classes.  All StRoot-independent
00252  * code has been moved to EEmcSmdGeom.  TVector3 replaces StThreeVectorD in
00253  * all function calls in EEmcSmdGeom.  StThreeVectorD wrappers are provided
00254  * in StEEmcSmdGeom, for integration into Star framework.
00255  *
00256  * Revisions:
00257  *
00258  * 01/28/04 Jason Webb -- Renamed to EEmcSmdGeom, StRoot dependent code moved 
00259  * to a derived class StEEmcSmdGeom.  The user interface for StEEmcSmdGeom
00260  * should remain unchanged.  Revision history for StEEmcSmdGeom moved  to end
00261  * of header file.  
00262  *
00263  * Log: StEEmcSmdGeom.h,v 
00264  * Revision 1.6  2003/12/05 00:06:11  jwebb
00265  * Member function added to return a vector pointing to the intersection of
00266  * two strips.
00267  *
00268  * Revision 1.5  2003/10/15 15:26:03  wzhang
00269  * improved and reorganized
00270  *
00271  * Revision 1.4  2003/08/22 15:14:03  wzhang
00272  * Added ClassDef and method stripEnd
00273  *
00274  * Revision 1.3  2003/06/11 18:58:15  wzhang
00275  * added geometry methods for StiEEmc
00276  *
00277  * Revision 1.2  2003/04/04 15:33:31  wzhang
00278  * included EEmcGeomDefs.h & improved codes
00279  *
00280  * Revision 1.1  2003/03/28 15:50:00  balewski
00281  * first
00282  *
00283  *
00284  ****************************************************************/