StRoot: StTrsMaker/include/StTrsWireHistogram.hh Source File

00001 /***************************************************************************
00002  *
00003  * $Id: StTrsWireHistogram.hh,v 1.9 2001/02/15 21:34:48 perev Exp $
00004  *
00005  * Author: brian, May 1998 
00006  ***************************************************************************
00007  *
00008  * Description: Collection of all StTrsMiniChargeSegment transported to
00009  *              the pad-plane
00010  *
00011  ***************************************************************************
00012  *
00013  * $Log: StTrsWireHistogram.hh,v $
00014  * Revision 1.9  2001/02/15 21:34:48  perev
00015  * clear improved
00016  *
00017  * Revision 1.8  2000/07/30 02:43:13  long
00018  * add erf function look up table,table builder,etc.
00019  *
00020  * Revision 1.7  1999/12/08 02:10:25  calderon
00021  * Modified to eliminate warnings on Linux.
00022  *
00023  * Revision 1.6  1999/07/19 21:39:31  lasiuk
00024  * - addEntry() distributes charge on a (user) settable range of wires
00025  * - setRangeOfWiresForChargeDistribution(int) added (default is 0)
00026  *
00027  * Revision 1.5  1999/07/09 03:46:49  lasiuk
00028  * add switch for singleElectron multiplication, gaussian random
00029  * number generator
00030  *
00031  * Revision 1.4  1999/06/16 14:26:51  fisyak
00032  * Add flags for egcs on Solaris
00033  *
00034  * Revision 1.3  1999/02/10 18:03:27  lasiuk
00035  * gas gain manual setting
00036  *
00037  * Revision 1.2  1998/11/16 14:50:22  lasiuk
00038  * remove diagnostics (mLastEntry, mLastWire)
00039  *
00040  * Revision 1.1  1998/11/10 17:12:13  fisyak
00041  * Put Brian trs versin into StRoot
00042  *
00043  * Revision 1.5  1998/11/08 17:07:53  lasiuk
00044  * change from boolean macros
00045  * use resize() for LINUX/ allocators for SUN
00046  * add typedefs for vector<> types
00047  *
00048  * Revision 1.4  1998/11/03 17:31:47  lasiuk
00049  * incorporate gas gain/fluctuations
00050  * add time delay of collection depending on charge position wrt wire position
00051  * rename wire()
00052  *
00053  * Revision 1.3  1998/10/22 00:23:26  lasiuk
00054  * Oct 22
00055  *
00056  * Revision 1.2  1998/06/30 22:56:59  lasiuk
00057  * added clear(), wire() and wireHistogram()
00058  *
00059  * Revision 1.1  1998/06/04 23:32:21  lasiuk
00060  * Initial Revision
00061  *
00062  **************************************************************************/
00063 #ifndef ST_TRS_WIRE_HISTOGRAM_HH
00064 #define ST_TRS_WIRE_HISTOGRAM_HH
00065 
00066 #include <vector>
00067 
00068 #if  defined(__sun) && ! defined(__GNUG__)
00069 #include <stdcomp.h>
00070 #endif
00071 
00072 #include "StThreeVector.hh"
00073 #include "Randomize.h"
00074 #include "TRandom.h"
00075 
00076 #include "StTpcGeometry.hh"
00077 #include "StTpcSlowControl.hh"
00078 #include "StTrsDeDx.hh"
00079 #include "StMagneticField.hh"
00080 #include "StTrsWireBinEntry.hh"
00081 
00082 #ifndef ST_NO_TEMPLATE_DEF_ARGS
00083 typedef vector<StTrsWireBinEntry> aTpcWire;
00084 typedef vector<aTpcWire>          aTpcWirePlane;
00085 #else
00086 typedef vector<StTrsWireBinEntry, allocator<StTrsWireBinEntry> > aTpcWire;
00087 typedef vector<aTpcWire, allocator<aTpcWire> >                   aTpcWirePlane;
00088 #endif
00089 
00090     
00091 class StTrsWireHistogram {
00092 public:
00093     static StTrsWireHistogram* instance(StTpcGeometry*, StTpcSlowControl*, StTrsDeDx*,StMagneticField*);
00094     static void dropit(){delete mInstance;mInstance=0;}
00095     ~StTrsWireHistogram();
00096 
00097     //StTrsWireHistogram(const StTrsWireHistogramy&);
00098     //StTrsWireHistogram& operator=(cont StTrsWireHistogram&);
00099 
00100     // access functions
00101     int minWire() const;
00102     int maxWire() const;
00103     
00104     // Book-keeping
00105     void                addEntry(StTrsWireBinEntry&,int );
00106     void                clear();
00107     
00108     // Wire Operations (and access)
00109     double              wireCoordinate(int) const;
00110     aTpcWire&           getWire(int)             ;
00111     aTpcWirePlane&      getWireHistogram()       ;
00112     void    FreqFunctionTableBuilder();
00113     double  table_fast(double) const;
00114  
00115 #ifndef ST_NO_TEMPLATE_DEF_ARGS 
00116   vector<double> mFreqFunctionTable; 
00117   
00118 #else
00119 
00120   vector<double, allocator<double> > mFreqFunctionTable; 
00121   
00122   
00123 #endif
00124   double         mNumberOfEntriesInTable;
00125   double         mRangeOfTable ;
00126     //Maybe needed later?
00127     //void sortTimeBins();
00128     //void putWire(vector<StTrsWireBinEntry>&);
00129 
00130   
00131     
00132     // Gas Gain
00133     void   setDoGasGain(bool)            ;
00134     void   setGasGainInnerSector(double) ;
00135     void   setGasGainOuterSector(double) ;
00136     void   setDoGasGainFluctuations(bool);
00137     void   setDoSingleElectronMultiplication(bool);
00138     double exponentialAvalanche(int, double)      ;
00139     double polyaAvalanche(int, double)      ; 
00140   double polya()      ;//b=0.4
00141     double gaussianAvalanche(int, double)         ;
00142   
00143     // Time Delay For Charge Collection
00144     void   setDoTimeDelay(bool);
00145     
00146 private:
00147     void   gasGainCalculation()                ;
00148     
00149     double noFluctuations(int)            const;
00150     TRandom *random;  
00151 private:
00152     StTrsWireHistogram(StTpcGeometry*, StTpcSlowControl*, StTrsDeDx*,StMagneticField*);
00153     
00154 private:
00155     int             mMin;                      // minimum bin filled
00156     int             mMax;                      // maximum bin filled
00157 
00158     StTpcGeometry*    mGeomDb;
00159     StTpcSlowControl* mSCDb;
00160     StTrsDeDx*        mGasDb;
00161      
00162   //  void   setRangeOfWiresForChargeDistribution(int);
00163 
00164     // Gas Gain
00165     bool            mDoGasGain;
00166     bool            mDoGasGainFluctuations;
00167     bool            mGasGainCalculationDone;
00168     bool            mDoSingleElectronMultiplication;
00169 
00170     // Time Delay
00171     bool            mDoTimeDelay;
00172     
00173     // Charge Distribution
00174     
00175     aTpcWirePlane      mSectorWires;
00176     
00177     int             mNumberOfInnerSectorAnodeWires; // from dataBase
00178     int             mNumberOfOuterSectorAnodeWires; // from dataBase
00179     int             mTotalNumberOfAnodeWires;
00180 
00181     double          mInnerSectorGasGain;
00182     double          mOuterSectorGasGain;
00183    
00184     static StTrsWireHistogram* mInstance;
00185 
00186     static HepJamesRandom  mEngine;
00187     static RandGauss       mGaussianDistribution;
00188     static RandExponential mExponentialDistribution;
00189      float           mOmegaTau;
00190      float dx[4]; 
00191      float dz[4];
00192     
00193 };
00194 int inline StTrsWireHistogram::minWire() const {return mMin;}
00195 int inline StTrsWireHistogram::maxWire() const {return mMax;}
00196 void inline StTrsWireHistogram::setDoGasGain(bool gg) {mDoGasGain = gg;}
00197 void inline StTrsWireHistogram::setDoTimeDelay(bool t) {mDoTimeDelay = t;}
00198 void inline StTrsWireHistogram::setDoSingleElectronMultiplication(bool t) {mDoSingleElectronMultiplication = t;}
00199 
00200 
00201 
00202 
00203 
00204 ;
00205 
00206 #endif