Still fixing AliITSsimulationSSD.h for missing AliITSTableSSD's.
[u/mrichter/AliRoot.git] / ITS / AliITSsimulationSSD.h
1 #ifndef ALIITSSIMULATIONSSD_H
2 #define ALIITSSIMULATIONSSD_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice     */
5 /* $Id$ */
6
7 #include "AliITSsimulation.h"
8 #include "AliITSsegmentationSSD.h" // function used in inline functions
9
10 class AliITSMapA2;
11 class AliITSpList;
12 class AliITSTableSSD;
13 class AliITSdcsSSD;
14 class AliITSsegmentationSSD;
15 class AliITSresponseSSD;
16
17 class AliITSsimulationSSD: public AliITSsimulation {
18
19  public:
20     AliITSsimulationSSD(); // Default constructor
21     AliITSsimulationSSD(const AliITSsimulationSSD &source); // copy constructor
22     // operator =
23     AliITSsimulationSSD& operator=(const AliITSsimulationSSD &source);
24     //Standard Constructor
25     AliITSsimulationSSD(AliITSsegmentation *seg,AliITSresponse *resp);
26     //Destructor
27     virtual ~AliITSsimulationSSD();
28     // Initilize variables for this simulation
29     void Init(AliITSsegmentationSSD *seg,AliITSresponseSSD *resp);
30     //Digitizes all of the hits in a module
31     void DigitiseModule(AliITSmodule *mod,Int_t dummy0,Int_t dummy1);
32     // Computes the Summable Digits
33     void SDigitiseModule(AliITSmodule *mod,Int_t module,Int_t dummy);
34     // Computes the Charge on each Strip/ Analog/summable digits
35     void HitsToAnalogDigits(AliITSmodule *mod,AliITSpList *pList);
36     //Computes the signal from one hit
37     void HitToDigit(Int_t module,Double_t x0,Double_t y0,Double_t z0, 
38                     Double_t x,Double_t y,Double_t z,Double_t de,
39                     AliITSTableSSD *tav);
40     //returns a pointer to the SSD segmentation.
41     AliITSsegmentationSSD *GetSegmentation() {
42         return (AliITSsegmentationSSD*) fSegmentation;}
43     //Returns the ionization energy for Si in GeV.
44     Double_t GetIonizeE() const {return fIonE;}
45     //Sets the ionization energy for Si in GeV.
46     void SetIonizeE(Double_t e=3.62E-09){fIonE = e;}
47     //Returns the Diffusion constant h in cm**2/sec
48     Double_t GetDiffConst(Int_t i) const {return fDifConst[i];}
49     //Sets the Diffusion constant h in cm**2/sec
50     void SetDiffConst(Double_t h=11.0,Double_t e=30.0)
51         {fDifConst[0] = h;fDifConst[1]=e;}
52     //Returns the Drift velocity for the side i
53     Double_t GetDriftVelocity(Int_t i) const {return fDriftVel[i];}
54     //Sets the Drift velocity for the P and N sides
55     void SetDriftVelocity(Double_t v0=0.86E+06,Double_t v1=2.28E+06)
56         {fDriftVel[0] = v0;fDriftVel[1] = v1;}
57     // Standard ascii class print function
58     void Print(ostream *os);
59     // Standard ascii class read function
60     void Read(istream *is);
61
62  private:
63     // returns the number of steps needed to proplerly distribute the charge
64     // in a step
65     Int_t NumOfSteps(Double_t x,Double_t y,Double_t z,
66                      Double_t  &dex,Double_t &dey,Double_t &dez);
67     // Keepts track and orders tracks for a give strip.
68     void GetList(Int_t trk,Int_t ht,Int_t mod,AliITSpList *pLt,
69                  AliITSTableSSD *tav);
70     // sets thresholds and fills digits
71     void ChargeToSignal(AliITSpList *pList);
72     // Writes Summable Digits to a root file for later use.
73     void WriteSDigits(AliITSpList *pList);
74     // ReadSDigits and create Digits
75     void SDigitToDigit(Int_t module,AliITSpList *pList);
76     // Fills fMapA2 from pList AliITSpList
77     void FillMapFrompList(AliITSpList *pList);
78     // Diffuses the charge onto neighboring strips.
79     void    IntegrateGaussian(Int_t k,Double_t par,Double_t av,Double_t sigma, 
80                               Double_t inf, Double_t sup,
81                               AliITSTableSSD *tav);
82      // Applies noise to strips randomly
83     void    ApplyNoise(AliITSpList *pList,Int_t mod);
84      // Applies posible signal coupling between strips
85     void    ApplyCoupling(AliITSpList *pList,Int_t mod);
86     // Computes the integral of a gaussian using Error Function
87     Float_t F(Float_t av, Float_t x, Float_t s);
88     // returns, from the segmentation, the number of stips
89     Int_t GetNStrips() {return GetSegmentation()->Npx();}
90     // returns, from the segmentation, the strip pitch
91     Float_t GetStripPitch() {return GetSegmentation()->Dpx(0);}
92
93     // Data members
94  protected:
95     AliITSdcsSSD *fDCS;   // Class containing detector controle paramters
96
97  private:
98     AliITSMapA2 *fMapA2;      //! Map of ionization, used localy only
99     Double_t    fIonE;        // ionization energy of Si in GeV
100     Double_t    fDifConst[2]; // Diffusion constants [h,e] in cm**2/sec
101     Double_t    fDriftVel[2]; // Drift velocities [P,N sides] cm/sec
102
103     ClassDef(AliITSsimulationSSD,2) // SSD signal simulation class
104
105 };
106 // Input and output functions for standard C++ input/output.
107 ostream &operator<<(ostream &os,AliITSsimulationSSD &source);
108 istream &operator>>(istream &is,AliITSsimulationSSD &source);
109 #endif