[u/mrichter/AliRoot.git] / ITS / AliITSresponseSDD.h

#ifndef ALIITSRESPONSESDD_H
#define ALIITSRESPONSESDD_H
 
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */
 
/*
  $Id$
*/

#include "TArrayF.h"
#include "AliITSresponse.h"

// response for SDD

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(Float_t p1) {// Adc-count saturation value
	fMaxAdc=p1;}
    virtual Float_t MaxAdc() const {// Get maximum Adc-count value
	return fMaxAdc;}
    virtual void    SetChargeLoss(Float_t p1) {
	// Set Linear Charge Loss Steepness  // 0.01 for 20%
	fChargeLoss=p1;}
    Float_t ChargeLoss() const {// Get Charge Loss Coefficient
	return fChargeLoss;}
    virtual void    SetDynamicRange(Float_t p1) {// Set Dynamic Range
	fDynamicRange = p1;}
    virtual Float_t DynamicRange() const {// Get Dynamic Range
	return fDynamicRange;}
    // Diffusion coefficients
    virtual void    SetDiffCoeff(Float_t p1, Float_t p2) 
	{fDiffCoeff=p1;fDiffCoeff1=p2;}
    // Get diffusion coefficients
    virtual void DiffCoeff(Float_t&diff,Float_t&diff1) const 
	{diff = fDiffCoeff;diff1 = fDiffCoeff1;}
    virtual void    SetDriftSpeed(Float_t p1) {// Drift velocity
	fDriftSpeed=p1;}
    virtual Float_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(Float_t n, Float_t b){
	// Noise and baseline  // 10 for ALICE with beam test measurements 2001
	fNoise=n; fBaseline=b;}
    virtual void  GetNoiseParam(Float_t &n, Float_t &b) const {// get noise param
	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(Float_t  mv, Float_t /* b */){
       // Min value used in 2D - could be used as a threshold setting
	fMinVal = mv;}
    virtual void   Thresholds(Float_t &  mv, Float_t & /* b */) const 
      {mv = fMinVal;}
    virtual void  GiveCompressParam(Int_t *x) const;
    // Detector type response methods
    virtual void    SetNSigmaIntegration(Float_t p1) {
	// Set number of sigmas over which cluster disintegration is performed
	fNsigmas=p1;}
    virtual Float_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 Float_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
    void  SetNoiseAfterElectronics(Float_t n=2.38){
	// Noise after electronics (ADC units)
	// 2.36 for ALICE from beam test measurements 2001
	fNoiseAfterEl=n;}
    Float_t  GetNoiseAfterElectronics() const {
	// Noise after electronics (ADC units)
	return fNoiseAfterEl;}
    void  SetJitterError(Float_t jitter=20) {
	// set Jitter error (20 um for ALICE from beam test measurements 2001)
	fJitterError=jitter;}
    Float_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; }
    Float_t Gain(Int_t mod,Int_t chip,Int_t ch)const {return fGain[mod][chip][ch]; }
    // these functions should be move to AliITSsegmentationSDD
    const Int_t Modules() const{return fgkModules;}// Total number of SDD modules
    const Int_t Chips() const{return fgkChips;} // Number of chips/module
    const 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(Float_t *)
      {NotImplemented("SetDetParam");}
    virtual void   GetDetParam(Float_t *) const 
      {NotImplemented("GetDetParam");}
    virtual  void   SetNDetParam(Int_t /* n */)
      {NotImplemented("SetNDetParam");}
    virtual Int_t  NDetParam() const
      {NotImplemented("NDetParam"); return 0;}
    virtual void    SetSigmaSpread(Float_t, Float_t) 
      {NotImplemented("SetSigmaSpread");}
    virtual void    SigmaSpread(Float_t & /* p1 */,Float_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 Float_t fgkDynamicRangeDefault; // default for fDynamicRange
    static const Float_t fgkfChargeLossDefault; // default for fChargeLoss
    static const Float_t fgkDiffCoeffDefault; // default for fDiffCoeff
    static const Float_t fgkDiffCoeff1Default; // default for fDiffCoeff1 
    static const Float_t fgkTemperatureDefault; // default for fT (Kelvin)
    static const TString fgkParam1Default; // default for fParam1
    static const TString fgkParam2Default; // default for fParam2
    static const Float_t fgkNoiseDefault; // default for fNoise
    static const Float_t fgkBaselineDefault; // default for fBaseline
    static const TString fgkOptionDefault; // default for fOption
    static const Float_t fgkMinValDefault; // default for fMinVal
    static const Float_t fgkDriftSpeedDefault; // default for fDriftSpeed
    static const Float_t fgkNsigmasDefault; //default for fNsigmas
    static const Int_t fgkNcompsDefault; //default for fNcomps
    //

    Int_t fDeadModules;                   // Total number of dead SDD modules
    Int_t fDeadChips;                     // Number of dead chips
    Int_t fDeadChannels;                  // Number of dead channels
    Float_t   fGain[fgkModules][fgkChips][fgkChannels];   // Array for channel gains
    Int_t     fCPar[8];        // Hardware compression parameters
    Float_t   fNoise;          // Noise
    Float_t   fBaseline;       // Baseline
    Float_t   fNoiseAfterEl;   // Noise after electronics
    Float_t   fJitterError;    // jitter error
    Float_t   fDynamicRange;   // Set Dynamic Range 
    Float_t   fChargeLoss;     // Set Linear Coefficient for Charge Loss 
    Float_t   fDriftSpeed;     // Drift velocity
    Int_t     fElectronics;    // Electronics
    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
    TArrayF   *fGaus;          // Gaussian lookup table for signal generation
    Int_t      fNcomps;        // Number of samplings along the gaussian
    Float_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

 private:
    AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
    AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.


    ClassDef(AliITSresponseSDD,3) // SDD response 
    
    };
#endif
Commit	Line	Data
b0f5e3fc	1	#ifndef ALIITSRESPONSESDD_H
b0f5e3fc	2	#define ALIITSRESPONSESDD_H
48058160	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	/*
	8	$Id$
	9	*/
88cb7938	10
7551c5b2	11	#include "TArrayF.h"
b0f5e3fc	12	#include "AliITSresponse.h"
	13
	14	// response for SDD
	15
48058160	16	class AliITSresponseSDD : public AliITSresponse {
48058160	17	public:
50d05d7b	18	//
	19	// Configuration methods
	20	//
b0f5e3fc	21
48058160	22	AliITSresponseSDD();
	23	AliITSresponseSDD(const char *dataType);
	24	virtual ~AliITSresponseSDD();
7551c5b2	25
703a4e51	26	// Implementation of virtual member functions declared in AliITSresponse
	27	virtual void SetElectronics(Int_t p1=1)
	28	{fElectronics=p1; /* Electronics: Pascal (1) or OLA (2) */ }
	29	virtual Int_t Electronics() const {// Electronics: 1 = Pascal; 2 = OLA
48058160	30	return fElectronics;}
703a4e51	31	virtual void SetMaxAdc(Float_t p1) {// Adc-count saturation value
48058160	32	fMaxAdc=p1;}
703a4e51	33	virtual Float_t MaxAdc() const {// Get maximum Adc-count value
48058160	34	return fMaxAdc;}
703a4e51	35	virtual void SetChargeLoss(Float_t p1) {
48058160	36	// Set Linear Charge Loss Steepness // 0.01 for 20%
48058160	37	fChargeLoss=p1;}
bff7ba71	38	Float_t ChargeLoss() const {// Get Charge Loss Coefficient
48058160	39	return fChargeLoss;}
703a4e51	40	virtual void SetDynamicRange(Float_t p1) {// Set Dynamic Range
	41	fDynamicRange = p1;}
	42	virtual Float_t DynamicRange() const {// Get Dynamic Range
48058160	43	return fDynamicRange;}
703a4e51	44	// Diffusion coefficients
	45	virtual void SetDiffCoeff(Float_t p1, Float_t p2)
	46	{fDiffCoeff=p1;fDiffCoeff1=p2;}
	47	// Get diffusion coefficients
	48	virtual void DiffCoeff(Float_t&diff,Float_t&diff1) const
	49	{diff = fDiffCoeff;diff1 = fDiffCoeff1;}
	50	virtual void SetDriftSpeed(Float_t p1) {// Drift velocity
48058160	51	fDriftSpeed=p1;}
703a4e51	52	virtual Float_t DriftSpeed() const {// drift speed
48058160	53	return fDriftSpeed;}
703a4e51	54	virtual void SetParamOptions(const char opt1,const char opt2){
48058160	55	// Parameters: "same" or read from "file"
48058160	56	fParam1=opt1; fParam2=opt2;}
703a4e51	57	virtual void ParamOptions(char opt1,char opt2) const {// options
48058160	58	strcpy(opt1,fParam1.Data()); strcpy(opt2,fParam2.Data());}
703a4e51	59	virtual void SetNoiseParam(Float_t n, Float_t b){
48058160	60	// Noise and baseline // 10 for ALICE with beam test measurements 2001
48058160	61	fNoise=n; fBaseline=b;}
703a4e51	62	virtual void GetNoiseParam(Float_t &n, Float_t &b) const {// get noise param
	63	n=fNoise; b=fBaseline;}
	64	virtual Bool_t Do10to8() const {// get 10 to 8 compression option
	65	return fBitComp;}
	66	void SetZeroSupp (const char *opt) {
	67	// Zero-suppression option - could be 1D, 2D or non-ZS
	68	fOption=opt;}
	69	const char *ZeroSuppOption() const {// Get zero-suppression option
	70	return fOption.Data();}
	71	virtual void SetThresholds(Float_t mv, Float_t /* b */){
	72	// Min value used in 2D - could be used as a threshold setting
	73	fMinVal = mv;}
	74	virtual void Thresholds(Float_t & mv, Float_t & /* b */) const
	75	{mv = fMinVal;}
	76	virtual void GiveCompressParam(Int_t *x) const;
	77	// Detector type response methods
	78	virtual void SetNSigmaIntegration(Float_t p1) {
	79	// Set number of sigmas over which cluster disintegration is performed
	80	fNsigmas=p1;}
	81	virtual Float_t NSigmaIntegration() const {
	82	// Get number of sigmas over which cluster disintegration is performed
	83	return fNsigmas;}
	84	virtual void SetNLookUp(Int_t p1);
	85	// Get number of intervals in which the gaussian lookup table is divided
	86	virtual Int_t GausNLookUp() const {return fNcomps;}
	87	virtual Float_t GausLookUp(Int_t i) const {
	88	if(i<0 \|\| i>=fNcomps) return 0.;return fGaus->At(i);}
	89
	90	//Declaration of member functions peculiar to this class
48058160	91	void SetNoiseAfterElectronics(Float_t n=2.38){
	92	// Noise after electronics (ADC units)
	93	// 2.36 for ALICE from beam test measurements 2001
	94	fNoiseAfterEl=n;}
703a4e51	95	Float_t GetNoiseAfterElectronics() const {
48058160	96	// Noise after electronics (ADC units)
	97	return fNoiseAfterEl;}
	98	void SetJitterError(Float_t jitter=20) {
	99	// set Jitter error (20 um for ALICE from beam test measurements 2001)
	100	fJitterError=jitter;}
703a4e51	101	Float_t JitterError() const {// set Jitter error
48058160	102	return fJitterError;}
	103	void SetDo10to8(Bool_t bitcomp=kTRUE) {
	104	// set the option for 10 to 8 bit compression
	105	fBitComp = bitcomp;}
48058160	106	// Compression parameters
48058160	107	void SetCompressParam(Int_t cp[8]);
48058160	108	void SetDeadChannels(Int_t nmodules=0, Int_t nchips=0, Int_t nchannels=0);
703a4e51	109	Int_t GetDeadModules() const { return fDeadModules; }
	110	Int_t GetDeadChips() const { return fDeadChips; }
	111	Int_t GetDeadChannels() const { return fDeadChannels; }
	112	Float_t Gain(Int_t mod,Int_t chip,Int_t ch)const {return fGain[mod][chip][ch]; }
48058160	113	// these functions should be move to AliITSsegmentationSDD
703a4e51	114	const Int_t Modules() const{return fgkModules;}// Total number of SDD modules
	115	const Int_t Chips() const{return fgkChips;} // Number of chips/module
	116	const Int_t Channels() const { return fgkChannels;}// Number of channels/chip
48058160	117	//********
	118	void PrintGains();
	119	void Print();
03898a57	120
703a4e51	121	// not implemented virtual methods (devlared in the mother class
	122	virtual void SetDetParam(Float_t *)
	123	{NotImplemented("SetDetParam");}
	124	virtual void GetDetParam(Float_t *) const
	125	{NotImplemented("GetDetParam");}
	126	virtual void SetNDetParam(Int_t /* n */)
	127	{NotImplemented("SetNDetParam");}
	128	virtual Int_t NDetParam() const
	129	{NotImplemented("NDetParam"); return 0;}
	130	virtual void SetSigmaSpread(Float_t, Float_t)
	131	{NotImplemented("SetSigmaSpread");}
	132	virtual void SigmaSpread(Float_t & /* p1 /,Float_t & / p2 */) const
	133	{NotImplemented("SigmaSpread");}
2aea926d	134
48058160	135	protected:
48058160	136	// these statis const should be move to AliITSsegmentationSDD
703a4e51	137	static const Int_t fgkModules = 520; // Total number of SDD modules
	138	static const Int_t fgkChips = 4; // Number of chips/module
	139	static const Int_t fgkChannels = 64; // Number of channels/chip
48058160	140	//*******
703a4e51	141	static const Int_t fgkMaxAdcDefault; // default for fMaxAdc
	142	static const Float_t fgkDynamicRangeDefault; // default for fDynamicRange
	143	static const Float_t fgkfChargeLossDefault; // default for fChargeLoss
	144	static const Float_t fgkDiffCoeffDefault; // default for fDiffCoeff
	145	static const Float_t fgkDiffCoeff1Default; // default for fDiffCoeff1
	146	static const Float_t fgkTemperatureDefault; // default for fT (Kelvin)
	147	static const TString fgkParam1Default; // default for fParam1
	148	static const TString fgkParam2Default; // default for fParam2
	149	static const Float_t fgkNoiseDefault; // default for fNoise
	150	static const Float_t fgkBaselineDefault; // default for fBaseline
	151	static const TString fgkOptionDefault; // default for fOption
	152	static const Float_t fgkMinValDefault; // default for fMinVal
	153	static const Float_t fgkDriftSpeedDefault; // default for fDriftSpeed
	154	static const Float_t fgkNsigmasDefault; //default for fNsigmas
	155	static const Int_t fgkNcompsDefault; //default for fNcomps
	156	//
	157
48058160	158	Int_t fDeadModules; // Total number of dead SDD modules
	159	Int_t fDeadChips; // Number of dead chips
	160	Int_t fDeadChannels; // Number of dead channels
703a4e51	161	Float_t fGain[fgkModules][fgkChips][fgkChannels]; // Array for channel gains
48058160	162	Int_t fCPar[8]; // Hardware compression parameters
	163	Float_t fNoise; // Noise
	164	Float_t fBaseline; // Baseline
	165	Float_t fNoiseAfterEl; // Noise after electronics
	166	Float_t fJitterError; // jitter error
	167	Float_t fDynamicRange; // Set Dynamic Range
	168	Float_t fChargeLoss; // Set Linear Coefficient for Charge Loss
48058160	169	Float_t fDriftSpeed; // Drift velocity
48058160	170	Int_t fElectronics; // Electronics
703a4e51	171	Float_t fMaxAdc; // Adc saturation value
	172	Float_t fDiffCoeff; // Diffusion Coefficient (scaling the time)
	173	Float_t fDiffCoeff1; // Diffusion Coefficient (constant term)
	174	Float_t fNsigmas; // Number of sigmas over which charge disintegration
48058160	175	// is performed
	176	TArrayF *fGaus; // Gaussian lookup table for signal generation
	177	Int_t fNcomps; // Number of samplings along the gaussian
703a4e51	178	Float_t fMinVal; // Min value used in 2D zero-suppression algo
48058160	179	Bool_t fBitComp; // 10 to 8 bit compression option
	180	TString fOption; // Zero-suppresion option (1D, 2D or none)
	181	TString fParam1; // Read baselines from file option
	182	TString fParam2; // Read compression algo thresholds from file
703a4e51	183
	184	private:
	185	AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
	186	AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.
	187
fa1750f9	188
800f9ca4	189	ClassDef(AliITSresponseSDD,3) // SDD response
50d05d7b	190
50d05d7b	191	};
b0f5e3fc	192	#endif