StRoot: StTpcDb/StTpcDb.h Source File

00001 /***************************************************************************
00002  *
00003  * $Id: StTpcDb.h,v 1.38 2012/05/03 23:56:48 fisyak Exp $
00004  *
00005  * Author:  David Hardtke
00006  ***************************************************************************
00007  *
00008  * Description: This is the interface between the database and the offline
00009  *              TPC software.  This classes takes care of the annoying
00010  *              calls to the root infrastucture, packages and manipulates
00011  *              the data, and returns the data to the user via simple
00012  *              interface classes.    
00013  *
00014  ***************************************************************************
00015  *
00016  * $Log: StTpcDb.h,v $
00017  * Revision 1.38  2012/05/03 23:56:48  fisyak
00018  * Set interpolation for one week only, fix sign of interpolation (thanks Gene), add TriggerId
00019  *
00020  * Revision 1.37  2011/07/21 16:48:53  fisyak
00021  * New schema for Sub Sector Alginement: SuperSectror position (defined by inner sub sector) and Outer sector position wrt SuperSectror position
00022  *
00023  * Revision 1.36  2011/01/18 14:39:43  fisyak
00024  * Clean up TpcDb interfaces and Tpc coordinate transformation
00025  *
00026  * Revision 1.35  2010/05/27 19:14:26  fisyak
00027  * Take out flavoring by 'sim' for tpcGlobalPosition,tpcSectorPosition and starClockOnl tables. remove usage tpcISTimeOffsets and tpcOSTimeOffsets tables
00028  *
00029  * Revision 1.34  2009/12/07 23:44:58  fisyak
00030  * Drop coordinate transformation for fortran, remove TpcHitErr
00031  *
00032  * Revision 1.33  2009/11/02 17:31:41  fisyak
00033  * use directly field from StarMagField, replace St_tpcGainC and St_tpcT0C by St_tpcPadGainT0C, add remove defaults in coordinate transformations
00034  *
00035  * Revision 1.32  2009/03/16 14:13:31  fisyak
00036  * Use StDetectorDb chairs for TpcGlobalPosition and TpcSectorPosition
00037  *
00038  * Revision 1.31  2008/09/10 15:46:36  fisyak
00039  * Recalculate Tpc drift velocity once per event, avoid expensive conversion to unix time
00040  *
00041  * Revision 1.30  2008/08/01 14:28:25  fisyak
00042  * Add new getT0, clean up
00043  *
00044  * Revision 1.29  2007/08/12 15:06:30  fisyak
00045  * Use separated East/West drift velocities only >= 2007, for back compartibility
00046  *
00047  * Revision 1.28  2007/07/12 20:21:09  fisyak
00048  * Drift velocity depends on TPC half, use online RHIC clock
00049  *
00050  * Revision 1.27  2007/03/21 17:27:01  fisyak
00051  * use TGeoHMatrix, change mode for switching drift velocities
00052  *
00053  * Revision 1.26  2004/10/27 21:44:28  fisyak
00054  * Add debug print for tables Validities, add access to ExB correction
00055  *
00056  * Revision 1.25  2004/02/23 00:35:00  fisyak
00057  * Add access to tpcPadResponse
00058  *
00059  * Revision 1.24  2004/01/14 22:54:30  fisyak
00060  * Add hooks for Pedestal and tpcGain
00061  *
00062  * Revision 1.23  2002/04/02 00:16:31  hardtke
00063  * New class that gets hit errors from database
00064  *
00065  * Revision 1.22  2002/02/06 18:39:13  hardtke
00066  * Add tpc Field Cage structure
00067  *
00068  * Revision 1.21  2001/08/14 18:18:03  hardtke
00069  * Add sector position structures
00070  *
00071  * Revision 1.20  2001/05/21 23:25:34  hardtke
00072  * Add tpcGlobalPosition to StTpcDb.  This includes the global position offset and the rotation w.r.t. the magnet
00073  *
00074  * Revision 1.19  2000/08/10 18:41:34  hardtke
00075  * only look for L0_trigger table once per event -- improves timing
00076  *
00077  * Revision 1.18  2000/08/09 14:54:54  hardtke
00078  * Add Clear option, set trigger table pointer to 0 after each event
00079  *
00080  * Revision 1.17  2000/08/08 19:15:23  hardtke
00081  * use correct trigger time offset in case of laser
00082  *
00083  * Revision 1.16  2000/05/11 17:17:27  hardtke
00084  * make trigger time offset available -- currently NOT different for beam and laser events
00085  *
00086  * Revision 1.15  2000/03/27 21:21:22  fine
00087  * Adjusted to ROOT 2.24
00088  *
00089  * Revision 1.14  2000/02/10 00:29:09  hardtke
00090  * Add tpg functions to StTpcDbMaker, fix a few bugs
00091  *
00092  * Revision 1.13  2000/01/25 16:01:10  fisyak
00093  * Devorce with StAF
00094  *
00095  * Revision 1.12  2000/01/11 15:49:52  hardtke
00096  * get Electronics table from Calibrations database, Fix error messages
00097  *
00098  * Revision 1.11  1999/12/16 22:00:53  hardtke
00099  * add CVS tags
00100  *
00101  **************************************************************************/
00102 #ifndef ClassStTpcDb
00103 #define ClassStTpcDb
00104 
00106 //                                                                      //
00107 // StTpcDb                                                              //
00108 //                                                                      //
00109 // This class implements to offline interface to the STAR database      //
00110 //                                                                      //
00112 
00113 #include "StMessMgr.h"
00114 #include "StDetectorDbMaker/St_tpcPadPlanesC.h"
00115 #include "StDetectorDbMaker/St_tpcWirePlanesC.h"
00116 #include "StDetectorDbMaker/St_tpcDimensionsC.h"
00117 #include "StDetectorDbMaker/St_tpcElectronicsC.h"
00118 #include "StDetectorDbMaker/St_tpcSlowControlSimC.h"
00119 #include "StDetectorDbMaker/St_tpcGlobalPositionC.h"
00120 #include "StDetectorDbMaker/St_tpcSectorPositionC.h"
00121 #include "StDetectorDbMaker/St_tpcFieldCageC.h"
00122 #include "StDetectorDbMaker/St_tpcPedestalC.h"
00123 #include "StDetectorDbMaker/St_tpcPadResponseC.h"
00124 #include "StDetectorDbMaker/St_tpcPadGainT0C.h"
00125 #include "StDbUtilities/StMagUtilities.h"
00126 #include "StDetectorDbMaker/St_trgTimeOffsetC.h"
00127 #include "TGeoMatrix.h"
00128 #include "TString.h"
00129 #include "tables/St_tpcDriftVelocity_Table.h"
00130 
00131 class StTpcDb;
00132 // Global pointers:
00133 R__EXTERN StTpcDb* gStTpcDb;
00134 class StTpcDb {
00135  public:    
00136   static StTpcDb* instance() {if (! gStTpcDb) new StTpcDb(); return gStTpcDb;}
00137   // Glob     = Global coordinate 
00138   // Tpc      = Tpc    -"-
00139   // SupS     = super sector misalignment(?)
00140   // SubS[io] = SubSector[io] misalignment
00141   // SecL     = sector -"- coordinate (y_p, x_p, DriftDistance - z_p);
00142   // Pad      = Pad -"- (x_p,y_p,z_p) (Sector12 coordinate system)
00143   // Tpc => Global is mTpc2GlobMatrix
00144   // Pad => SecL   is internal Flip matrix
00145   enum ETpcSectorRotationType {kUndefSector     =-2,
00146                                kFlip            =-1, // Flip * Subs[io] => SupS
00147                                kSupS2Tpc        = 0, // SupS => Tpc
00148                                kSupS2Glob       = 1, // SupS => Tpc => Glob; 
00149                                kSubSInner2SupS  = 2, // Subs[io] => SupS
00150                                kSubSOuter2SupS  = 3, // -"-
00151                                kSubSInner2Tpc   = 4, // (Subs[io] => SupS) => Tpc
00152                                kSubSOuter2Tpc   = 5, // -"-
00153                                kSubSInner2Glob  = 6, // (Subs[io] => SupS => Tpc) => Glob
00154                                kSubSOuter2Glob  = 7, // -"-
00155                                kPadInner2SupS   = 8, // (Pad => SecL) => (SubS[io] => SupS)
00156                                kPadOuter2SupS   = 9, // -"- 
00157                                kPadInner2Tpc    =10, // (Pad => SecL) => (SubS[io] => SupS => Tpc)
00158                                kPadOuter2Tpc    =11, // -"- 
00159                                kPadInner2Glob   =12, // (Pad => SecL) => (SubS[io] => SupS => Tpc => Glob)
00160                                kPadOuter2Glob   =13, // -"- 
00161                                
00162                                kTotalTpcSectorRotaions =14}; 
00163  private:
00164   Char_t                mBeg[1];         
00165   St_tpcDriftVelocity*  dvel;           
00166   StMagUtilities*       mExB;           
00167   Int_t                 m_Debug;        
00168   TGeoHMatrix          *mFlip;          
00169   TGeoHMatrix          *mTpc2GlobMatrix;
00170   TGeoHMatrix          *mTpcSectorRotations[24][kTotalTpcSectorRotaions]; 
00171   Float_t               mDriftVel[2];   
00172   UInt_t                mUc;            
00173   Int_t                 mTriggerId;     
00174   Char_t                mEnd[1];        
00175  private:
00176   StTpcDb();
00177  public:
00178   virtual ~StTpcDb();
00179   St_tpcPadPlanesC      *PadPlaneGeometry() {return St_tpcPadPlanesC::instance();}
00180   St_tpcWirePlanesC     *WirePlaneGeometry() {return St_tpcWirePlanesC::instance();}
00181   St_tpcDimensionsC     *Dimensions() {return St_tpcDimensionsC::instance();}
00182   St_tpcSlowControlSimC *SlowControlSim() {return St_tpcSlowControlSimC::instance();}
00183   St_tpcElectronicsC    *Electronics() {return St_tpcElectronicsC::instance();}
00184   St_tpcGlobalPositionC *GlobalPosition() {return St_tpcGlobalPositionC::instance();}
00185   St_tpcFieldCageC      *FieldCage() {return St_tpcFieldCageC::instance();}
00186   St_tpcSectorPositionC *SectorPosition() {return St_tpcSectorPositionC::instance();}
00187   St_tpcPedestalC       *Pedestal() {return St_tpcPedestalC::instance();}
00188   St_tpcPadGainT0C      *tpcGain() {return St_tpcPadGainT0C::instance();}
00189   St_tpcPadGainT0C      *tpcT0()   {return St_tpcPadGainT0C::instance();}
00190   St_tpcPadResponseC    *PadResponse() {return St_tpcPadResponseC::instance();}
00191   Float_t                triggerTimeOffset() {return St_trgTimeOffsetC::instance()->triggerTimeOffset();}
00192 
00193   //small pieces of data:
00194   void    SetDriftVelocity();
00195   Float_t DriftVelocity(Int_t sector=24);
00196   StMagUtilities* ExB() {return mExB;}
00197   void SetExB(StMagUtilities *m) {mExB = m;}
00198   void SetTpcRotations();
00199   void SetTpc2GlobalMatrix(TGeoHMatrix *m) {SetTpcRotationMatrix(m);}
00200   void SetTpcRotationMatrix(TGeoHMatrix *m, Int_t sector = 0, Int_t k = kSupS2Tpc) {
00201     if (sector == 0)  {if (m) *mTpc2GlobMatrix = *m;}
00202     else              {if (m) *mTpcSectorRotations[sector-1][k] = *m;}
00203   }
00204   void SetDebug(Int_t m) {m_Debug = m;}
00205   Int_t Debug() {return m_Debug;}
00206   void  SetTriggerId(Int_t m) {mTriggerId = m;}
00207   Int_t TriggerId() {return mTriggerId;}
00208   const TGeoHMatrix &Flip()                           const {return *mFlip;}
00209   const TGeoHMatrix &Tpc2GlobalMatrix()               const {return *mTpc2GlobMatrix;}
00210   const TGeoHMatrix &TpcRot(Int_t sector, Int_t k)    const {return *mTpcSectorRotations[sector-1][k];}
00211   const TGeoHMatrix &SupS2Tpc(Int_t sector = 1)       const {return TpcRot(sector,kSupS2Tpc);}
00212   const TGeoHMatrix &SupS2Glob(Int_t sector = 1)      const {return TpcRot(sector,kSupS2Glob);}
00213   const TGeoHMatrix &SubSInner2SupS(Int_t sector = 1) const {return TpcRot(sector,kSubSInner2SupS);}
00214   const TGeoHMatrix &SubSOuter2SupS(Int_t sector = 1) const {return TpcRot(sector,kSubSOuter2SupS);}
00215   const TGeoHMatrix &SubSInner2Tpc(Int_t sector = 1)  const {return TpcRot(sector,kSubSInner2Tpc);}
00216   const TGeoHMatrix &SubSOuter2Tpc(Int_t sector = 1)  const {return TpcRot(sector,kSubSOuter2Tpc);}
00217   const TGeoHMatrix &SubSInner2Glob(Int_t sector = 1) const {return TpcRot(sector,kSubSInner2Glob);}
00218   const TGeoHMatrix &SubSOuter2Glob(Int_t sector = 1) const {return TpcRot(sector,kSubSOuter2Glob);}
00219 
00220   const TGeoHMatrix &PadInner2SupS(Int_t sector = 1)  const {return TpcRot(sector,kPadInner2SupS);}
00221   const TGeoHMatrix &PadOuter2SupS(Int_t sector = 1)  const {return TpcRot(sector,kPadOuter2SupS);}
00222   const TGeoHMatrix &PadInner2Tpc(Int_t sector = 1)   const {return TpcRot(sector,kPadInner2Tpc);}
00223   const TGeoHMatrix &PadOuter2Tpc(Int_t sector = 1)   const {return TpcRot(sector,kPadOuter2Tpc);}
00224   const TGeoHMatrix &PadInner2Glob(Int_t sector = 1)  const {return TpcRot(sector,kPadInner2Glob);}
00225   const TGeoHMatrix &PadOuter2Glob(Int_t sector = 1)  const {return TpcRot(sector,kPadOuter2Glob);}
00226 
00227   const TGeoHMatrix &SubS2SupS(Int_t sector = 1, Int_t row = 1) const {Int_t k = (row <= 13) ? kSubSInner2SupS : kSubSOuter2SupS; return TpcRot(sector,k);}
00228   const TGeoHMatrix &SubS2Tpc(Int_t sector = 1, Int_t row = 1)  const {Int_t k = (row <= 13) ? kSubSInner2Tpc : kSubSOuter2Tpc; return TpcRot(sector,k);}
00229   const TGeoHMatrix &SubS2Glob(Int_t sector = 1, Int_t row = 1) const {Int_t k = (row <= 13) ? kSubSInner2Glob: kSubSOuter2Glob; return TpcRot(sector,k);}
00230 
00231   const TGeoHMatrix &Pad2SupS(Int_t sector = 1, Int_t row = 1)  const {Int_t k = (row <= 13) ? kPadInner2SupS: kPadOuter2SupS; return TpcRot(sector,k);}
00232   const TGeoHMatrix &Pad2Tpc(Int_t sector = 1, Int_t row = 1)   const {Int_t k = (row <= 13) ? kPadInner2Tpc: kPadOuter2Tpc; return TpcRot(sector,k);}
00233   const TGeoHMatrix &Pad2Glob(Int_t sector = 1, Int_t row = 1)  const {Int_t k = (row <= 13) ? kPadInner2Glob: kPadOuter2Glob; return TpcRot(sector,k);}
00234   ClassDef(StTpcDb,0)
00235 };
00236 #endif