#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>
-// response for SDD
+/* $Id$ */
-class AliITSresponseSDD :
- public AliITSresponse {
-public:
- //
- // Configuration methods
+/////////////////////////////////////////////////////////////
+// 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() {
- // destructor
- }
- AliITSresponseSDD(const AliITSresponseSDD &source); // copy constructor
- AliITSresponseSDD& operator=(const AliITSresponseSDD &source); // ass. op.
-
- virtual void SetMaxAdc(Float_t p1=1023) {
- // Adc-count saturation value
- fMaxAdc=p1;
- }
- virtual Float_t MaxAdc() {
- // Get maximum Adc-count value
- return fMaxAdc;
- }
-
- virtual void SetMagicValue(Float_t p1=450.) {
- // Set maximum Adc-magic value
- fTopValue=p1;
- }
- virtual Float_t MagicValue() {
- // Get maximum Adc-magic value
- return fTopValue;
- }
-
- virtual void SetDiffCoeff(Float_t p1=5.) {
- // Diffusion coefficient
- fDiffCoeff=p1;
- }
- virtual Float_t DiffCoeff() {
- // Get diffusion coefficient
- return fDiffCoeff;
- }
-
- virtual void SetQref(Float_t p1=120.) {
- // Coulomb repulsion
- fQref=p1;
- }
- virtual Float_t Qref() {
- // qref
- return fQref;
- }
-
- 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(char *data="simulated") {
- // Type of data - real or simulated
- fDataType=data;
- }
- virtual const char *DataType() {
- // Get data type
- return fDataType.Data();
- }
-
- virtual void SetParamOptions(Option_t *opt1="same",Option_t *opt2="same"){
- // Parameters: "same" or read from "file"
- fParam1=opt1; fParam2=opt2;
- }
- virtual void ParamOptions(Option_t *&opt1,Option_t *&opt2) {
- // options
- opt1=fParam1; opt2=fParam2;
- }
-
- virtual void SetNoiseParam(Float_t n=3., Float_t b=20.){
- // 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 SetZeroSupp(Option_t *opt="2D") {
- // Zero-suppression option - could be 1D, 2D or non-ZS
- fOption=opt;
- }
- virtual Option_t *ZeroSuppOption() {
- // Get zero-suppression option
- return fOption;
- }
- virtual void SetMinVal(Int_t mv=4) {
- // Min value used in 2D - could be used as a threshold setting
- fMinVal = mv;
- }
- virtual Int_t MinVal() {
- // min val
- return fMinVal;
- }
-
- virtual void SetFilenames(char *f1=0,char *f2=0, char *f3=0) {
- // Set filenames - input, output, parameters ....
- fFileName1=f1; fFileName2=f2; fFileName3=f3;
- }
- virtual void Filenames(const char*&input,const char*&baseline,const char*¶m) {
- // Filenames
- input=fFileName1; baseline=fFileName2; param=fFileName3;
- }
-
-
- virtual void SetOutputOption(Bool_t write=kFALSE) {
- // set output option
- fWrite = write;
- }
- Bool_t OutputOption() {
- // output option
- return fWrite;
+ 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];
}
//
- // Compression parameters
- virtual void SetCompressParam(Int_t cp[8]);
- void GiveCompressParam(Int_t *x);
-
- //
- // Detector type response methods
- virtual void SetNSigmaIntegration(Float_t p1) {
- // Set number of sigmas over which cluster disintegration is performed
- }
- virtual Float_t NSigmaIntegration() {
- // Get number of sigmas over which cluster disintegration is performed
- return 0.;
- }
- 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.;
- }
-
-
-protected:
-
- Int_t fCPar[8]; // Hardware compression parameters
- //Int_t fNDetPar; // Number of detector param
- //Float_t fDetPar[fNDetPar];
-
- 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
- Float_t fQref; // Coulomb repulsion
-
- 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
- Option_t *fOption; // Zero-suppresion option (1D, 2D or none)
- Option_t *fParam1; // Read baselines from file option
- Option_t *fParam2; // Read compression algo thresholds from file
-
- TString fDataType; // input keys : run, module #
- 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
+ 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;
+ protected:
+ //
+ 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.
+
+ 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.
+
+ ClassDef(AliITSresponseSDD,21)
+
+ };
+#endif