]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSsimulationSDD.h
Adding the track fit residuals as a consequence of the ExB distortions (Marian)
[u/mrichter/AliRoot.git] / ITS / AliITSsimulationSDD.h
index bc0e770d2dd1703a26d3da9e9d0418f10ae517ce..661ff2b0adcff705161acf28cd1b28b162f3e3be 100644 (file)
@@ -3,16 +3,18 @@
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice     */
 
+
 /* $Id$ */
 
-#include <TNtuple.h>
-#include <TArrayF.h>
+////////////////////////////////////////////////////////////
+// Simulation class for SDD                               //
+////////////////////////////////////////////////////////////
 
+#include <TNtuple.h>
 #include "AliITSsimulation.h"
 
 class TH1F;
 class TFile;
-class TVector;
 class TArrayI;
 class TArrayF;
 class AliITS;
@@ -23,25 +25,27 @@ class AliITSMapA2;
 class AliITSetfSDD;
 class AliITSsegmentationSDD;
 class AliITSInStream;
-class AliITSresponse;
-class AliITSresponseSDD;
+class AliITSCalibration;
+class AliITSCalibrationSDD;
 
 class AliITSsimulationSDD : public AliITSsimulation {
-
- public:
+  public:
     AliITSsimulationSDD(); // default constructor
     //Standard Constructor
-    AliITSsimulationSDD(AliITSsegmentation *seg, AliITSresponse *res);
-    // Copy opporator
-    AliITSsimulationSDD(AliITSsimulationSDD &source);
+    AliITSsimulationSDD(AliITSDetTypeSim* dettyp);
+    // Copy operator
+    AliITSsimulationSDD(const AliITSsimulationSDD &source);
     virtual ~AliITSsimulationSDD(); // Destructor
-    // = opporator
-    AliITSsimulationSDD& operator=(AliITSsimulationSDD &source);
+    // = operator
+    AliITSsimulationSDD& operator=(const AliITSsimulationSDD &source);
+    virtual AliITSsimulation& operator=(const AliITSsimulation &source);
     // Initilize variables for this simulation
-    void Init(AliITSsegmentationSDD *seg,AliITSresponseSDD *resp);
+    void Init();
 
-    // get the address of the array mapping the signal or pointers to arrays
-//    virtual AliITSMap*  HitMap(Int_t i);
+    // Get a pointer to the segmentation object
+    virtual AliITSsegmentation* GetSegmentationModel(Int_t /*dt*/){return fDetType->GetSegmentationModel(1);}
+    // set pointer to segmentation object
+    virtual void SetSegmentationModel(Int_t /*dt*/, AliITSsegmentation *seg){fDetType->SetSegmentationModel(1,seg);}
 
     // set the scale size factor for the smples in FFT
     virtual void SetScaleFourier(Int_t scale=4) {fScaleSize=scale;}
@@ -54,43 +58,37 @@ class AliITSsimulationSDD : public AliITSsimulation {
     virtual void SetCrosstalkFlag(Bool_t flag=kFALSE) {fCrosstalkFlag=flag;}
     // return crosstalk flag
     Bool_t CrosstalkFlag() const {return fCrosstalkFlag;}
-    // set compression parameters for 2D or 1D via response functions
-    void SetCompressParam();
     // retrieve compression parameters for 2D or 1D
-    void CompressionParam(Int_t i, Int_t &db, Int_t &tl, Int_t &th);
+    void CompressionParam(Int_t i, Int_t &db, Int_t &tl, Int_t &th) const;
     // retrieve compression parameters for 2D or 1D
-    void CompressionParam(Int_t i, Int_t &db, Int_t &tl);
+    void CompressionParam(Int_t i, Int_t &db, Int_t &tl) const;
 
-    virtual Int_t Convert10to8(Int_t signal) const;// 10 to 8 bit SDD compresion
-    virtual Int_t Convert8to10(Int_t signal) const;//undo 10 to 8 bit SDD compresion
+    virtual Int_t Convert10to8(Int_t signal) const;//10 to 8 bit SDD compresion
     virtual void ZeroSuppression(const char *opt); // Apply zero suppresion
     virtual void Init2D();   // initiilzes 2D compresion algorithm
     virtual void Compress2D(); // Applies 2D compresion algorithm
     virtual void Init1D();   // initilizes 1D compresion algorithm
     virtual void Compress1D(); // Applies 1D compresion algorithm
     virtual void StoreAllDigits(); // if No compresion run this.
-    virtual void ReadBaseline();  // read baseline values from a file
     // returns baseline and noise for a given anode i.
-    virtual void GetAnodeBaseline(Int_t i, Float_t &baseline, Float_t &noise);
+    //virtual void GetAnodeBaseline(Int_t i,Double_t &baseline,Double_t &noise) const;
     // local implementation of ITS->AddDigit. Specific for SDD
     virtual void AddDigit(Int_t i, Int_t j, Int_t signal);
     // Finds clulsters of signals. Use with regards to Compresion algorithms
     virtual void  FindCluster(Int_t i, Int_t j,Int_t signal,
-                             Int_t minval,Bool_t &cond);
+                              Int_t minval,Bool_t &cond);
 
     // get parameters for 1D - this could be changed when we get more
     // input from Torino after they have a look at the code 
-    virtual Int_t Tolerance(Int_t i) {return fTol[i];}//returns tolerance
-    virtual Int_t Disable(Int_t i)  {return fT2[i];}//high threshold  2D
+    virtual Int_t Tolerance(Int_t i) const {return fTol[i];}//returns tolerance
+    virtual Int_t Disable(Int_t i)  const {return fT2[i];}//high threshold  2D
     // Set the output file name - for 1D encoding 
     virtual void SetFileName(const char *filnam) {fFileName=filnam;}
 
-    // add baseline, noise, electronics and ADC saturation effects
-    void ChargeToSignal(Bool_t bAddNoise=kFALSE);
-    // add dead channels
-    void ApplyDeadChannels();
+    // add baseline, noise, gain, electronics and ADC saturation effects
+    void ChargeToSignal(Int_t mod,Bool_t bAddNoise=kFALSE, Bool_t bAddGain=kTRUE);
     // add crosstalk effect
-    void ApplyCrosstalk();
+    void ApplyCrosstalk(Int_t mod);
     
     // create maps to build the lists of tracks for each summable digit
     void InitSimulationModule( Int_t module, Int_t event );
@@ -111,10 +109,11 @@ class AliITSsimulationSDD : public AliITSsimulation {
     // Spread charge in a SDD module
     void HitsToAnalogDigits(AliITSmodule *mod);
     // Sorts tracks for the 3 most highly contributed one to be added to digit.
-    //void SortTracks(Int_t *tracks,Float_t *charges,Int_t *hits,Int_t ntracks);
+    //void SortTracks(Int_t *tracks,Float_t *charges,Int_t *hits
+    //                Int_t ntracks);
     // collects and returns the fired SDD cells (uses AliITSMapA2...).
     //void ListOfFiredCells(Int_t *arg,Double_t timeAmplitude,TObjArray *list,
-       //                TClonesArray *padr);
+    //           TClonesArray *padr);
 
     // Creates histograms of maps for debugging
     void CreateHistograms(Int_t scale);
@@ -126,7 +125,8 @@ class AliITSsimulationSDD : public AliITSsimulation {
     TObjArray*  GetHistArray() {return fHis;}
     // create a separate tree for background monitoring (2D) 
     virtual  void  MakeTreeB(Option_t *option="B") 
-       { fTreeB = new TNtuple("ntuple","2D backgr","nz:nl:nh:low:anode");}
+        { if(strstr(option,"B"))
+            fTreeB = new TNtuple("ntuple", "2D backgr","nz:nl:nh:low:anode");}
     // presently a dummy routine use TreeB() instead
     void GetTreeB(Int_t) { }
     // Return pointer to TreeB
@@ -143,18 +143,17 @@ class AliITSsimulationSDD : public AliITSsimulation {
     void SetDoFFT(Int_t doFFT=1) {fDoFFT=doFFT;}
 
     // Print SSD simulation Parameters
-    virtual void Print();
+    virtual void PrintStatus() const;
 
- private:
+  private:
+    // virtual void GetBaseline(Int_t mod);  // read baseline values from a file
+   // set compression parameters for 2D or 1D via response functions
+    void SetCompressParam();
     // Variables and pointers for local use only. Not Streamed out.
     AliITS         *fITS;          //! local pointer to ITS
-//    AliITSMapA1    *fHitMap1;      //! local pointer to map of digits
     AliITSMapA2    *fHitMap2;      //! local pointer to map of signals
     AliITSMapA2    *fHitSigMap2;   //! local pointer to map of signals
     AliITSMapA2    *fHitNoiMap2;   //! local pointer to map of signals
-//    AliITSpList    *fpList;        //! 
-//    TObjArray      *falist;        //
-//    TClonesArray   *fpadr;         //
     AliITSInStream *fStream;       //! input file stream
     AliITSetfSDD   *fElectronics;  //! local pointer to electronics simulation
     Double_t       *fInZR;         //! [fScaleSize*fMaxNofSamples] input of the
@@ -172,8 +171,8 @@ class AliITSsimulationSDD : public AliITSsimulation {
     TArrayI    fT1;           // low thresholds
     TArrayI    fT2;           // high thresholds(2D) or disable (1D) 
     TArrayI    fTol;          // tolerance
-    TArrayF    fBaseline;     // Baseline
-    TArrayF    fNoise;        // Noise value
+    // TArrayF    fBaseline;     // Baseline
+    //TArrayF    fNoise;        // Noise value
     TNtuple   *fTreeB;        // Background info tree for 2D
     TString    fParam;        // Compresion algorithm options
     TString    fFileName;     // File name for possible options above
@@ -184,10 +183,8 @@ class AliITSsimulationSDD : public AliITSsimulation {
     Int_t      fNofMaps;      // Number of anodes used ( 1-2*nanodes per wing )
     Int_t      fMaxNofSamples;// Number of time samples
     Int_t      fScaleSize;    // scale size factor for the samples in FFT
-//    Int_t      fModule;       //! in case bgr,noise,param,change module-by-mod.
-//    Int_t      fEvent;        //! solely for output from bgr monitoring of 2D
 
-  ClassDef(AliITSsimulationSDD,1)  // Simulation of SDD clusters
+    ClassDef(AliITSsimulationSDD,2)  // Simulation of SDD clusters
 
 };
 #endif