StRoot: StSvtSimulationMaker/StSvtSignal.hh Source File

00001 #ifndef STSVTSIGNAL_HH
00002 #define STSVTSIGNAL_HH
00003 
00004 
00005 #include <Stiostream.h>
00006 #include <string.h>
00007 #include <math.h>
00008 
00009 class StSvtElectronCloud;
00010 
00011 class StSvtSignal
00012 {
00013 public:
00014   StSvtSignal();
00015   ~StSvtSignal();
00016 
00017   void setOption(int option);
00018   void setAnodeTimeBinSizes(double timBinSize, double anodeSize);
00019   void setDriftVelocity(double driftVelocity);
00020   void setCloud(StSvtElectronCloud* elCloud);
00021   double chargeFraction(int an, double anHit);
00022   int timeCenterAndWidth(double anHit,double timeHit);
00023   void setTimeWidth(double timWidth);
00024   void calcConvSignal(double chargeOnAnode);
00025   void selemonSignal(int nMin,int nMax, double tStep, double charge);
00026   void rykovSignal(int nMin,int nMax, double tStep);
00027   void pasaRelatedStuff();
00028   void doPasaOnly(int option);
00029   void peakingTimeR();
00030   void halfWidthAtHalfMaxR();
00031   void peakingTimeS();
00032   void halfWidthAtHalfMaxS();
00033   void unNormPasaConst();
00034   void normPasaConst();
00035   void arrays();
00036   double signal(double t);
00037   double getShortSignal(double localTime);
00038   double getLongSignal(double localTime);
00039   double useArrays5And6(double ds1,double dsc);
00040   double useArrays3And4Or1And2(double ds1,double dsc);
00041   double numConvInt(int nMin , int n, int numOfIntPoints, double tStep, double t);
00042   double analConvInt(double tim, double sigmat, double tc);
00043   double simpsonInt(int numOfIntPoints,double lowlim, double step, double t); 
00044   double gausInput(double tim);
00045   double pasaRes(double tim);
00046   double freq(double num);
00047   double prob1(double anOrTimeDiff , double  sigma);
00048   double prob2(double num , double  sigma);
00049   int getLowTBin();
00050   int getHiTBin();
00051   int getTrackId() const                {return mTrackId;}
00052   double getTimeCenter();
00053   double getTimeWidth();
00054   double getPeak();
00055   double getMinUnderShoot();
00056   double getSignal(int n);
00057   void resetPeakAndUnderShoot();
00058   void resetSignal(int lBin, int hBin);
00059 
00060 private:
00061   double mTotalHitCharge;
00062   double mAnRightEdge;
00063   double mAnLeftEdge;
00064   double mFractionOfCharge;
00065   double mCollectedCharge;
00066 
00067   int mOption;
00068   int mLowTBin;
00069   int mHiTBin;
00070   int mTrackId;
00071 
00072   double mDriftVel;
00073   double mTimeBinSize;
00074   double mAnodeSize;
00075 
00076   double mSigmaMajor;
00077   double mSigmaMinor;
00078   double mTimeCenter;
00079   double mTimeWidth;
00080   double mPeakSignal;
00081   double mMinUnderShoot;
00082 
00083   double GAP_TWIDTH;
00084   double mPasaGain;
00085 
00086   //for selemon's piece of code
00087   double mPasaNorm;
00088   double mPeakTimeS;
00089   double mPasaMaxS;         
00090   double mFwhmS;
00091 
00092   //for Rykoves piece of code( fortran converted to c++)
00093   double mPasaMaxR;
00094   double mPeakTimeR;
00095   double mFwhmR;
00096     
00097   double mTau_s;
00098   double mTau_l;
00099 
00100   double mC1;
00101   double mC2;
00102   double mC3;
00103 
00104   double mSignal[128];    
00105   double mPasa[5];
00106   double mArray1[9];
00107   double mArray2[9];
00108   double mArray3[7];
00109   double mArray4[7];
00110   double mArray5[6];
00111   double mArray6[6];
00112 
00113   //ClassDef(StSvtSignal,1)
00114 
00115 };
00116 
00117 /*inline void StSvtSignal::arrays()
00118 {
00119  mArray1[0] = 2.46196981473530512524*::pow(0.1,10);
00120  mArray1[1] = 0.564189564831068821977;
00121  mArray1[2] = 7.46321056442269912687;
00122  mArray1[3] = 48.6371970985681366614;
00123  mArray1[4] = 196.520832956077098242;
00124  mArray1[5] = 526.445194995477358631;
00125  mArray1[6] = 934.528527171957607540;
00126  mArray1[7] = 1027.55188689515710272;
00127  mArray1[8] = 557.535335369399327526;
00128 
00129  mArray2[0] = 1.0;
00130  mArray2[1] = 13.2281951154744992508;
00131  mArray2[2] = 86.7072140885989742329;
00132  mArray2[3] = 354.937778887819891062;
00133  mArray2[4] = 975.708501743205489753;
00134  mArray2[5] = 1823.90916687909736289;
00135  mArray2[6] = 2246.33760818710981792;
00136  mArray2[7] = 1656.66309194161350182;
00137  mArray2[8] = 557.535340817727675546;
00138 
00139  mArray3[0] = 0.0;
00140  mArray3[1] = 0.564189583547755073984;
00141  mArray3[2] = 1.275366707599781044160;
00142  mArray3[3] = 5.019050422511804774140;
00143  mArray3[4] = 6.160210979930535851950;
00144  mArray3[5] = 7.409742699504489391600;
00145  mArray3[6] = 2.978866653721002406700;
00146 
00147  mArray4[0] = 1.0;
00148  mArray4[1] = 2.26052863220117276590;
00149  mArray4[2] = 9.39603524938001434673;
00150  mArray4[3] = 12.0489539808096656605;
00151  mArray4[4] = 17.0814450747565897222;
00152  mArray4[5] = 9.60896809063285878198;
00153  mArray4[6] = 3.36907645100081516050;
00154 
00155  mArray5[0] = 0.0;
00156  mArray5[1] = 9.60497373987051638749;
00157  mArray5[2] = 90.0260197203842689217;
00158  mArray5[3] = 2232.00534594684319226;
00159  mArray5[4] = 7003.32514112805075473;
00160  mArray5[5] = 55592.3013010394962768;
00161 
00162  mArray6[0] = 1.0;
00163  mArray6[1] = 33.5617141647503099647;
00164  mArray6[2] = 521.357949780152679795;
00165  mArray6[3] = 4594.32382970980127987;
00166  mArray6[4] = 22629.0000613890934246;
00167  mArray6[5] = 49267.3942608635921086;  
00168  
00169  }*/
00170 
00171 #endif