1 #ifndef ALIITSSIMULATIONSPD_H
2 #define ALIITSSIMULATIONSPD_H
4 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 ////////////////////////////////////////////////////////////
10 // Simulation class for SPD //
11 ////////////////////////////////////////////////////////////
13 #include "TObjArray.h"
14 #include "AliITSsimulation.h"
15 #include "AliITSCalibrationSPD.h"
16 #include "AliITSsegmentationSPD.h"
20 class AliITSresponseSPD;
22 //-------------------------------------------------------------------
24 class AliITSsimulationSPD : public AliITSsimulation {
26 AliITSsimulationSPD();
27 // AliITSsimulationSPD(AliITSsegmentation *seg,AliITSCalibration *res);
28 AliITSsimulationSPD(AliITSDetTypeSim *dettyp);
29 virtual ~AliITSsimulationSPD();
31 AliITSsimulationSPD(const AliITSsimulationSPD &source);
33 AliITSsimulationSPD& operator=(const AliITSsimulationSPD &s);
34 virtual AliITSsimulation& operator=(const AliITSsimulation &source);
35 // Initilizes the variables
38 // General User calling routines
39 // Initilize simulation for a specific event and module
40 void InitSimulationModule(Int_t module, Int_t event);
41 // Finish and write S Digitization
42 void FinishSDigitiseModule();
43 // From hits to Digits, without creating SDigits
44 void DigitiseModule(AliITSmodule *mod,Int_t,Int_t);
46 // More or less Internal Routines
47 // Create S Digits from specific module
48 void SDigitiseModule(AliITSmodule *mod, Int_t mask, Int_t event);
49 // Write S Digits to the tree of SDigits.
51 // fill pList from hits, charge sharing, diffusion, coupling
52 void HitToSDigit(AliITSmodule *mod);
53 // fill pList from hits, charge sharing, diffusion, coupling (faster method optimized by Bjorne)
54 void HitToSDigitFast(AliITSmodule *mod);
55 // Removes dead pixels from pList
56 void RemoveDeadPixels(AliITSmodule *mod);
57 // Take pList of signals and apply noise... create Digis
58 void FrompListToDigits();
60 void CreateHistograms();
61 void FillHistograms(Int_t ix,Int_t iz,Double_t v=1.0);
62 void ResetHistograms();
63 TH1F* GetHistogram(Int_t i){return (TH1F*)(fHis->At(i));}// get histogram
64 TObjArray* GetHistArray() {return fHis;}// get hist array
65 TString& GetHistName(){return fSPDname;}
66 void SetHistName(TString &n){fSPDname = n;}
68 // For backwards compatibility
69 void SDigitsToDigits(){ FinishSDigitiseModule();};
70 void HitToDigit(AliITSmodule *mod){
71 // Standard interface to DigitiseModule Inputs: AliITSmodule *mod Pointer to this module
72 DigitiseModule(mod,GetModuleNumber(),0);};
75 void SpreadCharge(Double_t x0,Double_t z0,Int_t ix0,Int_t iz0,
76 Double_t el,Double_t sig,Int_t t,Int_t hi);
77 void SpreadChargeAsym(Double_t x0,Double_t z0,Int_t ix0,Int_t iz0,
78 Double_t el,Double_t sigx,Double_t sigz,Int_t t,Int_t hi);
79 void UpdateMapSignal(Int_t ix,Int_t iz,Int_t trk,Int_t ht,Double_t signal){
80 // This function adds a signal to the pList from the pList class
81 // Inputs: iz column number ix row number trk track number ht hit number signal signal strength
82 GetMap()->AddSignal(iz,ix,trk,ht,GetModuleNumber(),signal);};
83 void UpdateMapNoise(Int_t ix,Int_t iz,Float_t noise){
84 // This function adds noise to data in the MapA2 as well as the pList
85 // Inputs: iz column number ix row number noise electronic noise generated by FrompListToDigits
86 GetMap()->AddNoise(iz,ix,GetModuleNumber(),noise);}
87 // Get a pointer to the segmentation object
88 virtual AliITSsegmentation* GetSegmentationModel(Int_t /*dt*/){return fDetType->GetSegmentationModel(0);}
89 // set pointer to segmentation objec
90 virtual void SetSegmentationModel(Int_t /*dt*/, AliITSsegmentation *seg){fDetType->SetSegmentationModel(0,seg);}
91 // Bari-Salerno Coupling parameters
92 // "New" coupling routine Tiziano Virgili
93 void SetCoupling(Int_t col,Int_t row,Int_t ntrack,Int_t idhit);
94 // "Old" coupling routine Rocco Caliandro
95 void SetCouplingOld(Int_t col, Int_t row,Int_t ntrack,Int_t idhit);
96 // Getters for data kept in fSegmentation and fResponse.
97 // Returns the Threshold in electrons
98 Double_t GetThreshold(){
99 Double_t th,sig;AliITSCalibrationSPD* res=(AliITSCalibrationSPD*)GetCalibrationModel(GetModuleNumber());
100 res->Thresholds(th,sig);return th;};
101 // Returns the couplings Columb and Row.
102 void GetCouplings(Double_t &cc,Double_t &cr){
103 AliITSCalibrationSPD* res = (AliITSCalibrationSPD*)GetCalibrationModel(GetModuleNumber());
104 res->GetCouplingParam(cc,cr);};
105 // Returns the number of pixels in x
106 Int_t GetNPixelsX(){return GetSegmentationModel(0)->Npx();};
107 // Returns the number of pixels in z
108 Int_t GetNPixelsZ(){return GetSegmentationModel(0)->Npz();};
110 TObjArray *fHis; //! just in case for histogramming
111 TString fSPDname; //! Histogram name
112 Int_t fCoupling; // Sets the coupling to be used.
113 // ==1 use SetCoupling, ==2 use SetCouplingOld
114 ClassDef(AliITSsimulationSPD,2) // Simulation of SPD clusters