Adding some further mother volumes to speed-up the overlap checking and particle...
[u/mrichter/AliRoot.git] / ITS / AliITSresponseSDD.h
index 0c67be9..b9d2f4e 100644 (file)
 #ifndef ALIITSRESPONSESDD_H
 #define ALIITSRESPONSESDD_H
+/* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
 
-#include "AliITSresponse.h"
+#include <TObject.h>
+#include <AliLog.h>
 
-class TString;
+/* $Id$ */
 
-// response for SDD
+/////////////////////////////////////////////////////////////
+//  Base settings for the ITS response classes.            //  
+//  The data member of this class are static and set once  //
+//  for all the modules.                                   //    
+///////////////////////////////////////////////////////////// 
+
+class AliITSresponseSDD : public TObject {
+ public:
 
-class AliITSresponseSDD :
-  public AliITSresponse {
-public:
-  //
-  // Configuration methods
-  //
-  
   AliITSresponseSDD();
-  virtual ~AliITSresponseSDD() { 
-    // destructor
-  }
-  AliITSresponseSDD(const AliITSresponseSDD &source); // copy constructor
-  AliITSresponseSDD& operator=(const AliITSresponseSDD &source); // ass. op.
-  
-  virtual void    SetMaxAdc(Float_t p1=1024) {
-    // Adc-count saturation value
-    fMaxAdc=p1;
-  }
-  virtual Float_t MaxAdc()  {
-    // Get maximum Adc-count value
-    return fMaxAdc;
-  }                       
-  
-  virtual void    SetMagicValue(Float_t p1=900.) {
-    // Set maximum Adc-top value
-    fTopValue=p1;
-    //it was 96.95
-  }
-  virtual Float_t MagicValue()  {
-    // Get maximum Adc-top value
-    return fTopValue;
-  }                       
-  
-  virtual void    SetDiffCoeff(Float_t p1=2.8,Float_t p2=28.) {
-    // Diffusion coefficients
-    fDiffCoeff=p1;
-    fDiffCoeff1=p2;
-  }
-  virtual void DiffCoeff(Float_t&diff,Float_t&diff1) {
-    // Get diffusion coefficients
-    diff = fDiffCoeff;
-    diff1 = fDiffCoeff1;
-  } 
-  
-  virtual void    SetDriftSpeed(Float_t p1=7.5) {
-    // Drift velocity
-    fDriftSpeed=p1;
-  }
-  virtual Float_t DriftSpeed() {
-    // drift speed
-    return fDriftSpeed;
-  } 
-  
-  virtual void    SetTemperature(Float_t p1=23.) {
-    // Temperature
-    fTemperature=p1;
-  }
-  virtual Float_t Temperature() {
-    // Get temperature
-    return fTemperature;
-  } 
-  
-  virtual void    SetDataType(const char *data="simulated") {
-    // Type of data - real or simulated
-    fDataType=data;
-  }
-  virtual const char  *DataType() const {
-    // Get data type
-    return fDataType.Data();
-  } 
-  
-  virtual void SetParamOptions(const char *opt1="same",const char *opt2="same"){
-    // Parameters: "same" or read from "file" 
-    fParam1=opt1; fParam2=opt2;
-  }
-  virtual void   ParamOptions(char *opt1,char *opt2) {
-    // options
-    strcpy(opt1,fParam1.Data()); strcpy(opt2,fParam2.Data());
-  }
-  
-  virtual  void  SetNoiseParam(Float_t n=1.75, Float_t b=10.){
-    // Noise and baseline
-    fNoise=n; fBaseline=b;
-  }   
-  virtual  void  GetNoiseParam(Float_t &n, Float_t &b) {
-    // get noise param
-    n=fNoise; b=fBaseline;
-  }  
-
-  virtual  void  SetDo10to8(Bool_t bitcomp=kTRUE) {
-    // set the option for 10 to 8 bit compression
-    fBitComp = bitcomp;
-  }
+  virtual ~AliITSresponseSDD(){};
 
-  Bool_t Do10to8() {
-    // get 10 to 8 compression option
-    return fBitComp;
-  }   
-  
-  virtual void    SetZeroSupp (const char *opt="2D") {
-    // Zero-suppression option - could be 1D, 2D or non-ZS 
-    fOption=opt;
-  }
-  virtual const char *ZeroSuppOption() const {
-    // Get zero-suppression option
-    return fOption.Data();
-  }
-  virtual  void  SetMinVal(Int_t mv=4) {
-    // Min value used in 2D - could be used as a threshold setting
-    fMinVal = mv;
+  virtual void SetSideATimeZero(Float_t tzero){
+    SetLayer3ATimeZero(tzero);
+    SetLayer4ATimeZero(tzero);
   }
-  virtual Int_t  MinVal() {
-    // min val
-    return fMinVal;
+  virtual void SetSideCTimeZero(Float_t tzero){
+    SetLayer3CTimeZero(tzero);
+    SetLayer4CTimeZero(tzero);
   }
-  
-  virtual void   SetFilenames(const char *f1="",const char *f2="",const char *f3="") {
-    // Set filenames - input, output, parameters ....
-    fFileName1=f1; fFileName2=f2; fFileName3=f3;
+  virtual void SetLayer3ATimeZero(Float_t tzero){
+    for(Int_t iLad=1; iLad<=kNLaddersLay3; iLad++) SetHalfLadderATimeZero(3,iLad,tzero);      
   }
-  virtual void   Filenames(char *input,char *baseline,char *param) {
-    // Filenames
-   strcpy(input,fFileName1.Data());  strcpy(baseline,fFileName2.Data());  
-   strcpy(param,fFileName3.Data());
-  }     
-  
-  
-  virtual  void  SetOutputOption(Bool_t write=kFALSE) {
-    // set output option
-    fWrite = write;
+  virtual void SetLayer3CTimeZero(Float_t tzero){
+    for(Int_t iLad=1; iLad<=kNLaddersLay3; iLad++) SetHalfLadderCTimeZero(3,iLad,tzero);
   }
-  Bool_t OutputOption()  {
-    // output option
-    return fWrite;
+  virtual void SetLayer4ATimeZero(Float_t tzero){
+    for(Int_t iLad=1; iLad<=kNLaddersLay4; iLad++) SetHalfLadderATimeZero(4,iLad,tzero);      
   }
-  // 
-  // Compression parameters
-  virtual  void  SetCompressParam(Int_t cp[8]); 
-  void  GiveCompressParam(Int_t *x);
-  
-  //  
-  // Detector type response methods
-  virtual void    SetNSigmaIntegration(Float_t p1=4.) {
-    // Set number of sigmas over which cluster disintegration is performed
-    fNsigmas=p1;
+  virtual void SetLayer4CTimeZero(Float_t tzero){
+    for(Int_t iLad=1; iLad<=kNLaddersLay4; iLad++) SetHalfLadderCTimeZero(4,iLad,tzero);
   }
-  virtual Float_t NSigmaIntegration() {
-    // Get number of sigmas over which cluster disintegration is performed
-    return fNsigmas;
+  virtual void SetHalfLadderATimeZero(Int_t lay, Int_t lad, Float_t tzero);
+  virtual void SetHalfLadderCTimeZero(Int_t lay, Int_t lad, Float_t tzero);
+  virtual void SetModuleTimeZero(Int_t modIndex, Float_t tzero){
+    if(modIndex<kNSPDmods || modIndex>kNSPDmods+kNSDDmods) AliError(Form("SDD module number %d out of range",modIndex));
+    fTimeZero[modIndex-kNSPDmods]=tzero;
   }
-  virtual void    SetSigmaSpread(Float_t p1, Float_t p2) {
-    // Set sigmas of the charge spread function
-  }
-  virtual void    SigmaSpread(Float_t &s1, Float_t &s2) {
-    // Get sigmas for the charge spread 
-  }
-  
-  virtual Float_t IntPH(Float_t eloss) {
-    // Pulse height from scored quantity (eloss)
-    return 0.;
-  }
-  virtual Float_t IntXZ(AliITSsegmentation *) {
-    // Charge disintegration 
-    return 0.;
+
+  virtual void SetDeltaVDrift(Int_t modIndex, Float_t dv){
+    if(modIndex<kNSPDmods || modIndex>kNSPDmods+kNSDDmods) AliError(Form("SDD module number %d out of range",modIndex));
+    fDeltaVDrift[modIndex-kNSPDmods]=dv;
   }
-  
-  
-protected:
-  
-  Int_t     fCPar[8];        // Hardware compression parameters
-  
-  Float_t   fNoise;          // Noise
-  Float_t   fBaseline;       // Baseline
-  Float_t   fTopValue;       // still unclear to me 
-  Float_t   fTemperature;    // Temperature 
-  Float_t   fDriftSpeed;     // Drift velocity
-  
-  Float_t    fMaxAdc;        // Adc saturation value
-  Float_t    fDiffCoeff;     // Diffusion Coefficient (scaling the time)
-  Float_t    fDiffCoeff1;    // Diffusion Coefficient (constant term)
-  Float_t    fNsigmas;       // Number of sigmas over which charge disintegration 
-                             // is performed 
-  
-  Int_t      fZeroSuppFlag;  // Zero-suppression flag
-  Int_t      fMinVal;        // Min value used in 2D zero-suppression algo
-  
-  Bool_t     fWrite;         // Write option for the compression algorithms
-  Bool_t     fBitComp;       // 10 to 8 bit compression option
-
-  TString    fOption;        // Zero-suppresion option (1D, 2D or none)
-  TString    fParam1;        // Read baselines from file option
-  TString    fParam2;        // Read compression algo thresholds from file 
-  
-  TString         fDataType;         // data type - real or simulated
-  TString         fFileName1;        // input keys : run, module #
-  TString         fFileName2;        // baseline & noise val or output coded                                                 // signal or monitored bgr.
-  TString         fFileName3;        // param values or output coded signal 
-  
-  ClassDef(AliITSresponseSDD,1) // SDD response 
-    
-    };
-#endif
 
+  virtual void SetTimeOffset(Float_t to){fTimeOffset = to;}
+  virtual Float_t GetTimeOffset()const {return fTimeOffset;}
+  virtual Float_t GetTimeZero(Int_t modIndex){
+    if(modIndex<kNSPDmods || modIndex>kNSPDmods+kNSDDmods){
+      AliError(Form("SDD module number %d out of range",modIndex));
+      return 0.;
+    }
+    return fTimeZero[modIndex-kNSPDmods];
+  }
+  virtual Float_t GetDeltaVDrift(Int_t modIndex){
+    if(modIndex<kNSPDmods || modIndex>kNSPDmods+kNSDDmods){
+      AliError(Form("SDD module number %d out of range",modIndex));
+      return 0.;
+    }
+    return fDeltaVDrift[modIndex-kNSPDmods];
+  }
+  static Float_t DefaultTimeOffset() {return fgkTimeOffsetDefault;}
 
+  virtual void SetADC2keV(Float_t conv){fADC2keV=conv;}
+  virtual Float_t GetADC2keV()const {return fADC2keV;}
+  static Float_t DefaulttADC2keV() {return fgkADC2keVDefault;}
 
+  static Float_t GetCarlosRXClockPeriod() {return fgkCarlosRXClockPeriod;}
 
+ protected:
 
+  enum {kNSPDmods = 240};
+  enum {kNSDDmods = 260};
+  enum {kNLaddersLay3 = 14};
+  enum {kNLaddersLay4 = 22};
 
 
+  static const Float_t fgkTimeOffsetDefault;   // default for fTimeOffset
+  static const Float_t fgkADC2keVDefault;      // default for fADC2keV
+  static const Float_t fgkCarlosRXClockPeriod; // clock period for CarlosRX
+
+  Float_t  fTimeOffset;             // Time offset due to electronic delays 
+                                    // --> obsolete, kept for backw. comp. 
+  Float_t  fTimeZero[kNSDDmods];    // Time Zero for each module
+  Float_t  fDeltaVDrift[kNSDDmods]; // Vdrift correction (um/ns) for each module
+  Float_t  fADC2keV;                // Conversion factor from ADC to keV
+
+ private:
+
+  AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
+  AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.
+
+  ClassDef(AliITSresponseSDD,17) 
+     
+    };
+#endif