#ifndef ALIITSRESPONSESDD_H
#define ALIITSRESPONSESDD_H
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+/* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
-/*
- $Id$
-*/
-#include "TArrayF.h"
-#include "AliITSresponse.h"
+#include <TObject.h>
+#include <AliLog.h>
-// response for SDD
+/* $Id$ */
+
+/////////////////////////////////////////////////////////////
+// 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:
+ enum {kVDCorr2Side = BIT(14),kVDCorrMult = BIT(15)}; // if bit set, the object contains separate corrections for 2 sides
+ //
+ AliITSresponseSDD();
+ virtual ~AliITSresponseSDD(){};
+
+ virtual void SetSideATimeZero(Float_t tzero){
+ SetLayer3ATimeZero(tzero);
+ SetLayer4ATimeZero(tzero);
+ }
+ virtual void SetSideCTimeZero(Float_t tzero){
+ SetLayer3CTimeZero(tzero);
+ SetLayer4CTimeZero(tzero);
+ }
+ virtual void SetLayer3ATimeZero(Float_t tzero){
+ for(Int_t iLad=1; iLad<=kNLaddersLay3; iLad++) SetHalfLadderATimeZero(3,iLad,tzero);
+ }
+ virtual void SetLayer3CTimeZero(Float_t tzero){
+ for(Int_t iLad=1; iLad<=kNLaddersLay3; iLad++) SetHalfLadderCTimeZero(3,iLad,tzero);
+ }
+ virtual void SetLayer4ATimeZero(Float_t tzero){
+ for(Int_t iLad=1; iLad<=kNLaddersLay4; iLad++) SetHalfLadderATimeZero(4,iLad,tzero);
+ }
+ virtual void SetLayer4CTimeZero(Float_t tzero){
+ for(Int_t iLad=1; iLad<=kNLaddersLay4; iLad++) SetHalfLadderCTimeZero(4,iLad,tzero);
+ }
+ 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(CheckModuleIndex(modIndex)) fTimeZero[modIndex-kNSPDmods]=tzero;
+ }
+
+ virtual void SetDeltaVDrift(Int_t modIndex, Float_t dv, Bool_t rightSide=kFALSE) {
+ int ind = GetVDIndex(modIndex,rightSide);
+ if (ind>=0) fDeltaVDrift[ind] = dv;
+ }
+
+ virtual Float_t GetDeltaVDrift(Int_t modIndex,Bool_t rightSide=kFALSE) const {
+ int ind = GetVDIndex(modIndex,rightSide);
+ return ind<0 ? 0.:fDeltaVDrift[ind];
+ }
+ //
+ Bool_t IsVDCorr2Side() const {return TestBit(kVDCorr2Side);}
+ Bool_t IsVDCorrMult() const {return TestBit(kVDCorrMult);}
+ void SetVDCorr2Side(Bool_t v=kTRUE) {SetBit(kVDCorr2Side,v);}
+ void SetVDCorrMult(Bool_t v=kTRUE) {SetBit(kVDCorrMult,v);}
+ //
+ static Float_t DefaultTimeOffset() {return fgkTimeOffsetDefault;}
+ virtual void SetTimeOffset(Float_t to){fTimeOffset = to;}
+ virtual Float_t GetTimeOffset()const {return fTimeOffset;}
+ virtual Float_t GetTimeZero(Int_t modIndex) const {
+ if(CheckModuleIndex(modIndex)) return fTimeZero[modIndex-kNSPDmods];
+ else return 0.;
+ }
+
+ virtual void SetADC2keV(Float_t conv){fADC2keV=conv;}
+ virtual Float_t GetADC2keV()const {return fADC2keV;}
+ virtual void SetADCtokeV(Int_t modIndex, Float_t conv){
+ if(CheckModuleIndex(modIndex)) fADCtokeV[modIndex-kNSPDmods]=conv;
+ }
+ virtual Float_t GetADCtokeV(Int_t modIndex) const {
+ if(CheckModuleIndex(modIndex)) return fADCtokeV[modIndex-kNSPDmods];
+ else return 0.;
+ }
+
+ virtual void SetChargevsTime(Float_t slope){fChargevsTime=slope;}
+ virtual Float_t GetChargevsTime()const {return fChargevsTime;}
+
+ virtual void SetADCvsDriftTime(Int_t modIndex, Float_t slope){
+ if(CheckModuleIndex(modIndex)) fADCvsDriftTime[modIndex-kNSPDmods]=slope;
+ }
+ virtual Float_t GetADCvsDriftTime(Int_t modIndex) const {
+ if(CheckModuleIndex(modIndex)) return fADCvsDriftTime[modIndex-kNSPDmods];
+ else return 0.;
+ }
+
+ static Float_t DefaultADC2keV() {return fgkADC2keVDefault;}
+ static Float_t DefaultChargevsTime() {return fgkChargevsTimeDefault;}
+ static Float_t DefaultADCvsDriftTime() {return fgkADCvsDrTimeDefault;}
+
+ static Float_t GetCarlosRXClockPeriod() {return fgkCarlosRXClockPeriod;}
+ void PrintChargeCalibrationParams() const;
+ void PrintTimeZeroes() const;
+ void PrintVdriftCorerctions() const;
-class AliITSresponseSDD : public AliITSresponse {
- public:
- //
- // Configuration methods
- //
- AliITSresponseSDD();
- AliITSresponseSDD(const char *dataType);
- virtual ~AliITSresponseSDD();
-
- // Implementation of virtual member functions declared in AliITSresponse
- virtual void SetElectronics(Int_t p1=1)
- {fElectronics=p1; /* Electronics: Pascal (1) or OLA (2) */ }
- virtual Int_t Electronics() const {// Electronics: 1 = Pascal; 2 = OLA
- return fElectronics;}
- virtual void SetMaxAdc(Double_t p1) {// Adc-count saturation value
- fMaxAdc=p1;}
- virtual Double_t MaxAdc() const {// Get maximum Adc-count value
- return fMaxAdc;}
- virtual void SetChargeLoss(Double_t p1) {
- // Set Linear Charge Loss Steepness // 0.01 for 20%
- fChargeLoss=p1;}
- Double_t ChargeLoss() const {// Get Charge Loss Coefficient
- return fChargeLoss;}
- virtual void SetDynamicRange(Double_t p1) {// Set Dynamic Range
- fDynamicRange = p1;}
- virtual Double_t DynamicRange() const {// Get Dynamic Range
- return fDynamicRange;}
- // Diffusion coefficients
- virtual void SetDiffCoeff(Double_t p1, Double_t p2)
- {fDiffCoeff=p1;fDiffCoeff1=p2;}
- // Get diffusion coefficients
- virtual void DiffCoeff(Double_t&diff,Double_t&diff1) const
- {diff = fDiffCoeff;diff1 = fDiffCoeff1;}
- virtual void SetDriftSpeed(Double_t p1) {// Drift velocity
- fDriftSpeed=p1;}
- virtual Double_t DriftSpeed() const {// drift speed
- return fDriftSpeed;}
- virtual void SetParamOptions(const char *opt1,const char *opt2){
- // Parameters: "same" or read from "file"
- fParam1=opt1; fParam2=opt2;}
- virtual void ParamOptions(char *opt1,char *opt2) const {// options
- strcpy(opt1,fParam1.Data()); strcpy(opt2,fParam2.Data());}
- virtual void SetNoiseParam(Double_t n, Double_t b){
- // Noise and baseline // 10 for ALICE with beam test measurements 2001
- fNoise=n; fBaseline=b;}
- // get noise param
- virtual void GetNoiseParam(Double_t &n, Double_t &b) const
- {n=fNoise; b=fBaseline;}
- virtual Bool_t Do10to8() const {// get 10 to 8 compression option
- return fBitComp;}
- void SetZeroSupp (const char *opt) {
- // Zero-suppression option - could be 1D, 2D or non-ZS
- fOption=opt;}
- const char *ZeroSuppOption() const {// Get zero-suppression option
- return fOption.Data();}
- virtual void SetThresholds(Double_t mv, Double_t /* b */){
- // Min value used in 2D - could be used as a threshold setting
- fMinVal = mv;}
- virtual void Thresholds(Double_t & mv, Double_t & /* b */) const
- {mv = fMinVal;}
- virtual void GiveCompressParam(Int_t *x) const;
- // Detector type response methods
- virtual void SetNSigmaIntegration(Double_t p1) {
- // Set number of sigmas over which cluster disintegration is performed
- fNsigmas=p1;}
- virtual Double_t NSigmaIntegration() const {
- // Get number of sigmas over which cluster disintegration is performed
- return fNsigmas;}
- virtual void SetNLookUp(Int_t p1);
- // Get number of intervals in which the gaussian lookup table is divided
- virtual Int_t GausNLookUp() const {return fNcomps;}
- virtual Double_t GausLookUp(Int_t i) const {
- if(i<0 || i>=fNcomps) return 0.;return fGaus->At(i);}
-
-//Declaration of member functions peculiar to this class
- Int_t Convert8to10(Int_t signal) const; //undo 10 to 8 bit SDD compresion
- void SetNoiseAfterElectronics(Double_t n=2.38){
- // Noise after electronics (ADC units)
- // 2.36 for ALICE from beam test measurements 2001
- fNoiseAfterEl=n;}
- Double_t GetNoiseAfterElectronics() const {
- // Noise after electronics (ADC units)
- return fNoiseAfterEl;}
- void SetJitterError(Double_t jitter=20) {
- // set Jitter error (20 um for ALICE from beam test measurements 2001)
- fJitterError=jitter;}
- Double_t JitterError() const {// set Jitter error
- return fJitterError;}
- void SetDo10to8(Bool_t bitcomp=kTRUE) {
- // set the option for 10 to 8 bit compression
- fBitComp = bitcomp;}
- // Compression parameters
- void SetCompressParam(Int_t cp[8]);
- void SetDeadChannels(Int_t nmodules=0, Int_t nchips=0, Int_t nchannels=0);
- Int_t GetDeadModules() const { return fDeadModules; }
- Int_t GetDeadChips() const { return fDeadChips; }
- Int_t GetDeadChannels() const { return fDeadChannels; }
- Double_t Gain(Int_t mod,Int_t chip,Int_t ch)const
- {return fGain[mod][chip][ch]; }
- // these functions should be move to AliITSsegmentationSDD
- Int_t Modules()const{return fgkModules;}//Total number of SDD modules
- Int_t Chips() const{return fgkChips;} // Number of chips/module
- Int_t Channels() const{ return fgkChannels;}//Number of channels/chip
- //********
- void PrintGains();
- void Print();
-
- // not implemented virtual methods (devlared in the mother class
- virtual void SetDetParam(Double_t *)
- {NotImplemented("SetDetParam");}
- virtual void GetDetParam(Double_t *) const
- {NotImplemented("GetDetParam");}
- virtual void SetNDetParam(Int_t /* n */)
- {NotImplemented("SetNDetParam");}
- virtual Int_t NDetParam() const
- {NotImplemented("NDetParam"); return 0;}
- virtual void SetSigmaSpread(Double_t, Double_t)
- {NotImplemented("SetSigmaSpread");}
- virtual void SigmaSpread(Double_t & /* p1 */,Double_t & /* p2 */) const
- {NotImplemented("SigmaSpread");}
protected:
- // these statis const should be move to AliITSsegmentationSDD
- static const Int_t fgkModules = 520; // Total number of SDD modules
- static const Int_t fgkChips = 4; // Number of chips/module
- static const Int_t fgkChannels = 64; // Number of channels/chip
- //*******
- static const Int_t fgkMaxAdcDefault; // default for fMaxAdc
- static const Double_t fgkDynamicRangeDefault; // default for fDynamicRange
- static const Double_t fgkfChargeLossDefault; // default for fChargeLoss
- static const Double_t fgkDiffCoeffDefault; // default for fDiffCoeff
- static const Double_t fgkDiffCoeff1Default; // default for fDiffCoeff1
- static const Double_t fgkTemperatureDefault; // default for fT (Kelvin)
- static const TString fgkParam1Default; // default for fParam1
- static const TString fgkParam2Default; // default for fParam2
- static const Double_t fgkNoiseDefault; // default for fNoise
- static const Double_t fgkBaselineDefault; // default for fBaseline
- static const TString fgkOptionDefault; // default for fOption
- static const Double_t fgkMinValDefault; // default for fMinVal
- static const Double_t fgkDriftSpeedDefault; // default for fDriftSpeed
- static const Double_t fgkNsigmasDefault; //default for fNsigmas
- static const Int_t fgkNcompsDefault; //default for fNcomps
//
+ virtual Int_t GetVDIndex(Int_t modIndex, Bool_t rightSide=kFALSE) const {
+ int ind = modIndex - kNSPDmods;
+ if(ind<0 || ind>=kNSDDmods) {AliError(Form("SDD module number %d out of range",modIndex)); return -1;}
+ return (rightSide && IsVDCorr2Side()) ? ind + kNSDDmods : ind;
+ }
+
+ virtual Bool_t CheckModuleIndex(Int_t modIndex) const {
+ if(modIndex<kNSPDmods || modIndex>=kNSPDmods+kNSDDmods){ AliError(Form("SDD module number %d out of range",modIndex)); return kFALSE;}
+ return kTRUE;
+ }
+
+ 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 fgkChargevsTimeDefault; // default for fChargevsTime
+ static const Float_t fgkADCvsDrTimeDefault; // default for fADCvsDriftTime
+ 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[2*kNSDDmods]; // Vdrift correction (um/ns) for each module left (<kNSDDmods) and right (>=kNSDDmods) sides
+ Float_t fADC2keV; // Conversion factor from ADC to keV
+ // --> obsolete, kept for backw. comp.
+ Float_t fChargevsTime; // --> obsolete, kept for backw. comp.
- Int_t fDeadModules; // Total number of dead SDD modules
- Int_t fDeadChips; // Number of dead chips
- Int_t fDeadChannels; // Number of dead channels
- Double_t fGain[fgkModules][fgkChips][fgkChannels];//Array for channel gains
- Int_t fCPar[8]; // Hardware compression parameters
- Double_t fNoise; // Noise
- Double_t fBaseline; // Baseline
- Double_t fNoiseAfterEl; // Noise after electronics
- Double_t fJitterError; // jitter error
- Double_t fDynamicRange; // Set Dynamic Range
- Double_t fChargeLoss; // Set Linear Coefficient for Charge Loss
- Double_t fDriftSpeed; // Drift velocity
- Int_t fElectronics; // Electronics
- Double_t fMaxAdc; // Adc saturation value
- Double_t fDiffCoeff; // Diffusion Coefficient (scaling the time)
- Double_t fDiffCoeff1; // Diffusion Coefficient (constant term)
- Double_t fNsigmas; // Number of sigmas over which charge disintegration
- // is performed
- TArrayF *fGaus; // Gaussian lookup table for signal generation
- Int_t fNcomps; // Number of samplings along the gaussian
- Double_t fMinVal; // Min value used in 2D zero-suppression algo
- 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
+ Float_t fADCvsDriftTime[kNSDDmods]; // Correction for zero suppression effect
+ Float_t fADCtokeV[kNSDDmods]; // ADC to keV conversion for each module
private:
- AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
- AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.
+ AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
+ AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.
- ClassDef(AliITSresponseSDD,3) // SDD response
-
+ ClassDef(AliITSresponseSDD,21)
+
};
#endif
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff