Changed AddSDigits to Module to return a Bool_t. True if there was a non noise
[u/mrichter/AliRoot.git] / ITS / AliITSsimulationSPD.h
CommitLineData
b0f5e3fc 1#ifndef ALIITSSIMULATIONSPD_H
2#define ALIITSSIMULATIONSPD_H
c7a4dac0 3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
b0f5e3fc 5
c7a4dac0 6/*
7$Id$
8*/
9#include "AliITSresponseSPD.h"
10#include "AliITSsegmentationSPD.h"
b0f5e3fc 11#include "AliITSsimulation.h"
12
13class AliITSMapA2;
c7a4dac0 14class AliITSpList;
b0f5e3fc 15class AliITSmodule;
16
17//-------------------------------------------------------------------
18
19class AliITSsimulationSPD : public AliITSsimulation {
20
c7a4dac0 21 public:
22 AliITSsimulationSPD(); // Default constructor
23 // Standard constructor
24 AliITSsimulationSPD(AliITSsegmentation *seg, AliITSresponse *res);
25 ~AliITSsimulationSPD();// destructor
26 AliITSsimulationSPD(const AliITSsimulationSPD &source); // copy constructo
27 // assignment operator
28 AliITSsimulationSPD& operator=(const AliITSsimulationSPD &source);
29 // Initilizes the variables
30 void Init(AliITSsegmentationSPD *seg, AliITSresponseSPD *resp);
b0f5e3fc 31
c7a4dac0 32 // Sum digitize module
3a97c582 33 // Create maps to build the lists of tracks for each summable digit
34 void InitSimulationModule(Int_t module,Int_t events);
3a97c582 35 // Digitize module from the sum of summable digits.
36 void FinishSDigitiseModule();
c7a4dac0 37 void SDigitiseModule(AliITSmodule *mod, Int_t dummy0,Int_t dummy1);
38 // digitize module. Also need to digitize modules with only noise.
39 void DigitiseModule(AliITSmodule *mod,Int_t dummy0, Int_t dummy1);
40 // sum digits to Digits.
41 void SDigitsToDigits(Int_t module,AliITSpList *pList);
42 // updates the Map of signal, adding the energy (ene) released by
43 // the current track
44 void UpdateMapSignal(Int_t row,Int_t col,Int_t trk,Int_t hit,Int_t mod,
45 Double_t ene,AliITSpList *pList);
46 // updates the Map of noise, adding the energy (ene) give my noise
47 void UpdateMapNoise(Int_t row,Int_t col,Int_t mod,Double_t ene,
48 AliITSpList *pList);
49 // Loops over all hits to produce Analog/floting point digits. This
50 // is also the first task in producing standard digits.
51 void HitsToAnalogDigits(AliITSmodule *mod,Int_t *frowpixel,
52 Int_t *fcolpixel,Double_t *fenepixel,
53 AliITSpList *pList);
54 // Steering function to determine the digits associated to a given
55 // hit (hitpos)
56 // The digits are created by charge sharing (ChargeSharing) and by
57 // capacitive coupling (SetCoupling). At all the created digits is
58 // associated the track number of the hit (ntrack)
59 void HitToDigit(AliITSmodule *mod, Int_t hitpos,Int_t *frowpixel,
60 Int_t *fcolpixel, Double_t *fenepixel,AliITSpList *pList);
61 // Take into account the geometrical charge sharing when the track
62 // crosses more than one pixel.
63 void ChargeSharing(Float_t x1l,Float_t z1l,Float_t x2l,Float_t z2l,
64 Int_t c1,Int_t r1,Int_t c2,Int_t r2,Float_t etot,
65 Int_t &npixel,Int_t *frowpixel,Int_t *fcolpixel,
66 Double_t *fenepixel);
67 // Take into account the coupling between adiacent pixels.
68 // The parameters probcol and probrow are the fractions of the
69 // signal in one pixel shared in the two adjacent pixels along
70 // the column and row direction, respectively.
71 void SetCoupling(Int_t row,Int_t col,Int_t ntrack,Int_t idhit,Int_t module,
72 AliITSpList *pList);
73 // The pixels are fired if the energy deposited inside them is above
74 // the threshold parameter ethr. Fired pixed are interpreted as digits
75 // and stored in the file digitfilename. One also needs to write out
76 // cases when there is only noise (nhits==0).
77 void CreateDigit(Int_t module,AliITSpList *pList);
78 // Set the electronic noise and threshold non-uniformities to all the
79 // pixels in a detector.
80 // The parameter fSigma is the squared sum of the sigma due to noise
81 // and the sigma of the threshold distribution among pixels.
82 void SetFluctuations(AliITSpList *pList,Int_t module);
83 // Apply a mask to the SPD module. 1% of the pixel channels are
84 // masked. When the database will be ready, the masked pixels
85 // should be read from it.
86 void SetMask();
87 // Create Histograms
88 void CreateHistograms();
89 // Reset histograms for this detector
90 void ResetHistograms();
91 // Fills the Summable digits Tree
92 void WriteSDigits(AliITSpList *pList);
93 // Fills fMap2A from the pList of Summable digits
94 void FillMapFrompList(AliITSpList *pList);
b0f5e3fc 95 // get hist array
c7a4dac0 96 TObjArray* GetHistArray() {return fHis;}
b0f5e3fc 97
c7a4dac0 98 private:
99 // Getters for data kept in fSegmentation and fResponse.
100 // Returns the Threshold in electrons
101 Double_t GetThreshold(){Float_t a=0.0,b=0.0;
102 ((AliITSresponseSPD*)fResponse)->Thresholds(a,b); return a;}
103 // Returns the threshold and rms noise.
104 void GetThresholds(Float_t &t,Float_t &s){
105 ((AliITSresponseSPD*)fResponse)->Thresholds(t,s);}
106 // Returns the couplings Columb and Row.
107 void GetCouplings(Float_t &cc,Float_t cr){
108 ((AliITSresponseSPD*)fResponse)->GetNoiseParam(cc,cr);}
109 // Returns the number of pixels in x
110 Int_t GetNPixelsX(){return ((AliITSsegmentationSPD*)fSegmentation)->Npx();}
111 // Returns the number of pixels in z
112 Int_t GetNPixelsZ(){return ((AliITSsegmentationSPD*)fSegmentation)->Npz();}
b0f5e3fc 113
c7a4dac0 114 private:
115/* Float_t fThresh; //! Threshold from fResponse
116 Float_t fSigma; //! Noise from fResponse
117 Float_t fCouplCol; //! Coupling along columns from fResponse
118 Float_t fCouplRow; //! Coupling along rows from fResponse
119 Int_t fNPixelsX; //! NPixelsX from fSegmentation
120 Int_t fNPixelsZ; //! NPixelsZ from fSegmentation
121*/
3a97c582 122// AliITSpList *fList; //! Internal map and sdigits
123// Int_t fModule; //! Module number being simulated
124// Int_t fEvent; //! Event number being simulated
c7a4dac0 125 AliITSMapA2 *fMapA2; //! MapA2 for Local internal use only
126 TObjArray *fHis; //! just in case for histogramming for Local
127 // internal use only
b0f5e3fc 128
c7a4dac0 129 ClassDef(AliITSsimulationSPD,1) // Simulation of SPD clusters
b0f5e3fc 130
131};
132
c7a4dac0 133#endif