StRoot: StTpcRSMaker/StTpcRSMaker.h Source File

00001 #ifndef STAR_ST_TRS_MAKER_HH
00002 #define STAR_ST_TRS_MAKER_HH
00003 
00005 //                                                                      //
00006 // StTpcRSMaker virtual base class for Maker                          //
00007 //                                                                      //
00009 #ifndef StMaker_H
00010 #include "StMaker.h"
00011 #endif
00012 #ifndef HEP_SYSTEM_OF_UNITS_H
00013 #include "SystemOfUnits.h"
00014 #endif
00015 #ifndef ST_NO_NAMESPACES
00016 using namespace units;
00017 #endif
00018 #include "StTpcRawData.h"
00019 #include "TF1F.h"
00020 #include "TH1.h"
00021 #include "TTree.h"
00022 #include "PAI.h"
00023 // g2t tables
00024 #include "tables/St_g2t_tpc_hit_Table.h"
00025 #include "tables/St_g2t_track_Table.h"
00026 #include "tables/St_g2t_vertex_Table.h" 
00027 #include "StTpcDb/StTpcDb.h"
00028 #include "StdEdxY2Maker/StTpcdEdxCorrection.h" 
00029 #include "StMagF.h"
00030 #include "TArrayF.h"
00031 class Altro;
00032 struct SignalSum_t {
00033   Float_t      Sum;
00034   Short_t      Adc;
00035   Short_t  TrackId;
00036 };
00037 class StTpcRSMaker : public StMaker {
00038  public:
00039   enum EMode {kPAI         = 0,// switch to PAI from GEANT
00040               kBICHSEL     = 1,// switch to Bichsel from GEANT 
00041               kGAINOAtALL  = 2,// do not use GAIN at all
00042               kdEdxCorr    = 3,// do use TpcdEdxCorrection
00043               kDistortion  = 4,// include distortions
00044               kNoToflight  = 5 // don't account for particle time of flight
00045   };
00046   StTpcRSMaker(const char *name="TpcRS");
00047   virtual              ~StTpcRSMaker();
00048   virtual Int_t         InitRun(int runnumber);
00049   virtual Int_t         Make();
00050   virtual Int_t         Finish();
00051   TF1F *GetShaperResponse(Int_t io = 0, Int_t sector = 1) {return (TF1F *) mShaperResponses[io][sector-1];}          
00052   TF1F *GetChargeFraction(Int_t io = 0)     {return (TF1F *) mChargeFraction[io];}     
00053   TF1F *GetPadResponseFunction(Int_t io = 0){return (TF1F *) mPadResponseFunction[io];}
00054   TF1F *GetPolya(Int_t io = 0)       {return (TF1F *) mPolya[io];}
00055   TF1F *GetTimeShape0(Int_t io = 0)  {return fgTimeShape0[io];}
00056   TF1F *GetTimeShape3(Int_t io = 0)  {return fgTimeShape3[io];}
00057   TF1  *GetHeed()                    {return mHeed;}
00058   Double_t GetNoPrimaryClusters(Double_t betaGamma, Int_t charge);
00059   virtual void Print(Option_t *option="") const;
00060   void DigitizeSector(Int_t sector);
00061   void SetLaserScale(Double_t m=1) {mLaserScale = m;}
00062   static Int_t    AsicThresholds(Short_t ADCs[__MaxNumberOfTimeBins__]);
00063   static Int_t    SearchT(const void *elem1, const void **elem2);
00064   static Int_t    CompareT(const void **elem1, const void **elem2);
00065   static Double_t shapeEI(Double_t *x, Double_t *par=0);
00066   static Double_t shapeEI_I(Double_t *x, Double_t *par=0);
00067   static Double_t shapeEI3(Double_t *x, Double_t *par=0);
00068   static Double_t shapeEI3_I(Double_t *x, Double_t *par=0);
00069   static Double_t fei(Double_t t, Double_t t0, Double_t T);
00070   static Double_t polya(Double_t *x, Double_t *par);
00071   SignalSum_t  *GetSignalSum();
00072   SignalSum_t  *ResetSignalSum();
00073   void SettauIntegrationX(Double_t p =      74.6e-9, Int_t io=0) {mtauIntegrationX[io] = p;}
00074   void SettauCX(Double_t           p =    1000.0e-9, Int_t io=0) {mtauCX[io] = p;}
00075   void SetCutEle(Double_t p = 1e-4)                  {mCutEle = p;}
00076   static Double_t Ec(Double_t *x, Double_t *p); // minimal energy to create an ion pair
00077   static TF1 *fEc(Double_t w = 26.2);           // HEED function to generate Ec
00078  private:
00079   static Double_t ShaperFunc(Double_t *x, Double_t *p);
00080   static Double_t PadResponseFunc(Double_t *x, Double_t *p);
00081   static Double_t Gatti(Double_t *x, Double_t *p);
00082   static Double_t InducedCharge(Double_t s, Double_t h, Double_t ra, Double_t Va, Double_t &t0);
00083   static TF1F     *fgTimeShape3[2];   
00084   static TF1F     *fgTimeShape0[2];   
00085   Char_t   beg[1];                    
00086   TTree   *fTree;                     
00087   SignalSum_t     *m_SignalSum;       
00088   TH1D*    mdNdx;                     
00089   TH1D*    mdNdE;                     
00090   TF1F  *mShaperResponses[2][24];     
00091   TF1F  *mChargeFraction[2];          
00092   TF1F  *mPadResponseFunction[2];     
00093   TF1F  *mPolya[2];                   
00094   TF1F  *mGG;                         
00095   TF1   *mHeed;                       
00096   StTpcdEdxCorrection *m_TpcdEdxCorrection; // !
00097   PAI  *mPAI;                         
00098   Double_t InnerAlphaVariation;       
00099   Double_t OuterAlphaVariation;       
00100   Altro *mAltro;                      
00101   Char_t end[1];                      
00102   Double_t             mLaserScale;   
00103   // local variables
00104   Int_t NumberOfInnerRows;            
00105   Int_t numberOfSectors;              
00106   Int_t NumberOfRows;                 
00107   Int_t NoPads;                       
00108   Int_t numberOfTimeBins;             
00109   Int_t    numberOfInnerSectorAnodeWires; 
00110   Double_t firstInnerSectorAnodeWire; 
00111   Double_t lastInnerSectorAnodeWire;  
00112   Int_t    numberOfOuterSectorAnodeWires; 
00113   Double_t firstOuterSectorAnodeWire; 
00114   Double_t lastOuterSectorAnodeWire;  
00115   Double_t anodeWirePitch;            
00116   Double_t numberOfElectronsPerADCcount; 
00117   Double_t anodeWireRadius;           
00118   const Double_t minSignal;           
00119   Double_t innerSectorAnodeVoltage;   
00120   Double_t outerSectorAnodeVoltage;   
00121   const Double_t ElectronRange;       
00122   const Double_t ElectronRangeEnergy; 
00123   const Double_t ElectronRangePower;  
00124   Double_t      mtauIntegrationX[2];  
00125   Double_t      mtauCX[2];            
00126   const Int_t NoOfSectors;            
00127   Int_t       NoOfRows;               
00128   Int_t       NoOfInnerRows;          
00129   const Int_t NoOfPads;               
00130   const Int_t NoOfTimeBins;           
00131   Double_t   mCutEle;                 
00132  public:    
00133   virtual const char *GetCVS() const {
00134     static const char cvs[]= 
00135       "Tag $Name:  $ $Id: StTpcRSMaker.h,v 1.25 2012/05/07 15:36:22 fisyak Exp $ built __DATE__ __TIME__"; 
00136       return cvs;
00137   }
00138   ClassDef(StTpcRSMaker,0)   //StAF chain virtual base class for Makers
00139 };
00140 #endif
00141 // $Id: StTpcRSMaker.h,v 1.25 2012/05/07 15:36:22 fisyak Exp $
00142 // $Log: StTpcRSMaker.h,v $
00143 // Revision 1.25  2012/05/07 15:36:22  fisyak
00144 // Remove hardcoded TPC parameters
00145 //
00146 // Revision 1.24  2012/04/03 14:05:18  fisyak
00147 // Speed up using  GetSaveL (__PAD_BLOCK__), sluggish shape histograms, Heed electron generation
00148 //
00149 // Revision 1.23  2011/12/13 17:23:22  fisyak
00150 // remove YXTProd, add WIREHISTOGRAM and WIREMAP, use particle definition from StarClassLibrary
00151 //
00152 // Revision 1.22  2011/10/14 23:27:51  fisyak
00153 // Back to standard version
00154 //
00155 // Revision 1.20  2011/09/18 22:39:48  fisyak
00156 // Extend dN/dx table (H.Bichsel 09/12/2011) to fix bug #2174 and #2181, clean-up
00157 //
00158 // Revision 1.19  2011/03/17 14:29:31  fisyak
00159 // Add extrapolation in region beta*gamma < 0.3
00160 //
00161 // Revision 1.18  2010/06/14 23:34:26  fisyak
00162 // Freeze at Version V
00163 //
00164 // Revision 1.17  2010/03/22 23:45:06  fisyak
00165 // Freeze version with new parameters table
00166 //
00167 // Revision 1.16  2010/03/17 15:53:16  fisyak
00168 // Move StTpcdEdxCorrection to StdEdxY2Maker to avoid dependence of StTpcDb on StDetectorDbMaker
00169 //
00170 // Revision 1.15  2010/03/16 19:41:46  fisyak
00171 // Move diffusion and sec/row correction in DB, clean up
00172 //
00173 // Revision 1.14  2010/02/26 18:53:33  fisyak
00174 // Take longitudinal Diffusion from Laser track fit, add Gating Grid
00175 //
00176 // Revision 1.13  2010/01/26 19:47:26  fisyak
00177 // Include dE/dx calibration and distortions in the simulation
00178 //
00179 // Revision 1.12  2009/10/30 21:12:00  fisyak
00180 // Freeze version rcf9108.F, Clean up
00181 //
00182 // Revision 1.11  2009/10/03 21:29:09  fisyak
00183 // Clean up, move all TpcT related macro into StTpcMcAnalysisMaker
00184 //
00185 // Revision 1.10  2009/09/21 13:20:39  fisyak
00186 // Variant O4, no mSigmaJitter, 100 keV
00187 //
00188 // Revision 1.9  2009/09/01 15:06:44  fisyak
00189 // Version N
00190 //
00191 // Revision 1.8  2009/08/24 20:16:41  fisyak
00192 // Freeze with new Altro parameters
00193 //
00194 // Revision 1.7  2008/12/29 15:24:55  fisyak
00195 // Freeze ~/WWW/star/Tpc/TpcRS/ComparisonMIP31
00196 //
00197 // Revision 1.6  2008/12/18 23:06:38  fisyak
00198 // Take care about references to TGiant
00199 //
00200 // Revision 1.5  2008/08/18 15:54:26  fisyak
00201 // Version 20
00202 //
00203 // Revision 1.4  2008/07/30 23:53:19  fisyak
00204 // Freeze
00205 //
00206 // Revision 1.3  2008/07/18 16:21:17  fisyak
00207 // Remove TF1F
00208 //
00209 // Revision 1.2  2008/06/19 22:45:43  fisyak
00210 // Freeze problem with TPX parameterization
00211 //
00212 // Revision 1.1.1.1  2008/04/28 14:39:47  fisyak
00213 // Start new Tpc Response Simulator
00214 //
00215 // Revision 1.14  2008/04/24 10:42:04  fisyak
00216 // Fix binning issues
00217 //
00218 // Revision 1.13  2008/04/04 15:00:11  fisyak
00219 // Freeze before shaper modifications
00220 //
00221 // Revision 1.12  2005/02/07 21:40:31  fisyak
00222 // rename antique TGeant3 to TGiant3
00223 //
00224 // Revision 1.11  2005/01/26 21:45:31  fisyak
00225 // Freeze correction made in June
00226 //
00227 // Revision 1.9  2004/05/29 21:16:27  fisyak
00228 // Fix pad direction, add sorting for ADC/cluster nonlinearity, replace product by sum of logs
00229 //
00230 // Revision 1.8  2004/05/04 13:39:06  fisyak
00231 // Add TF1
00232 //
00233 // Revision 1.7  2004/04/22 01:05:03  fisyak
00234 // Freeze the version before modification parametrization for K3
00235 //
00236 // Revision 1.6  2004/04/06 01:50:13  fisyak
00237 // Switch from Double_t to Float_t for sum
00238 //
00239 // Revision 1.5  2004/03/30 19:30:04  fisyak
00240 // Add Laser
00241 //
00242 // Revision 1.4  2004/03/21 19:00:43  fisyak
00243 // switch to GEANT step length
00244 //
00245 // Revision 1.3  2004/03/20 17:57:15  fisyak
00246 // Freeze the version of PAI model
00247 //
00248 // Revision 1.2  2004/03/17 20:47:43  fisyak
00249 // Add version with TrsCluster TTree
00250 //
00251 // Revision 1.1.1.1  2004/03/05 20:51:25  fisyak
00252 // replacement for Trs
00253 //