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

#ifndef ALIITSCALIBRATIONSSD_H
#define ALIITSCALIBRATIONSSD_H
/* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */ 
#include "AliITSCalibration.h"
#include "AliITSNoiseSSD.h"
#include "AliITSPedestalSSD.h"
#include "AliITSGainSSD.h"
#include "AliITSBadChannelsSSD.h"
#include "AliITSresponseSSD.h"
#include "TArrayF.h"
#include "TArrayI.h"

/* $Id$ */
//////////////////////////////////////////////
// Response class for SSD                   //
//                                          //
//////////////////////////////////////////////
class AliITSCalibrationSSD : public AliITSCalibration {

 public:
    AliITSCalibrationSSD();
    AliITSCalibrationSSD(const char *dataType);
    virtual ~AliITSCalibrationSSD();

    virtual  void   SetNoiseParam(Double_t np, Double_t nn) {
    	// set noise par
    	fNoiseP=np; fNoiseN=nn;
    }

    virtual void    GetNoiseParam(Double_t &np, Double_t &nn) const {
    	// get noise par
    	np=fNoiseP; nn=fNoiseN;
    }

    void AddNoiseP(Int_t c, Float_t n) { fNoise->AddNoiseP(c,n);}       
    TArrayF GetNoiseP() {return fNoise->GetNoiseP(); }
    Float_t GetNoiseP(Int_t n) {return fNoise->GetNoiseP(n); }
    void AddNoiseN(Int_t c, Float_t n) { fNoise->AddNoiseN(c,n);}
    TArrayF GetNoiseN() {return fNoise->GetNoiseN(); }
    Float_t GetNoiseN(Int_t n) {return fNoise->GetNoiseN(n); }
    void SetNoise( AliITSNoiseSSD* noise) {fNoise=noise;}

    void AddPedestalP(Int_t c, Float_t n) { fPedestal->AddPedestalP(c,n);}       
    TArrayF GetPedestalP() {return fPedestal->GetPedestalP(); }
    Float_t GetPedestalP(Int_t n) {return fPedestal->GetPedestalP(n); }
    void AddPedestalN(Int_t c, Float_t n) { fPedestal->AddPedestalN(c,n);}
    TArrayF GetPedestalN() {return fPedestal->GetPedestalN(); }
    Float_t GetPedestalN(Int_t n) {return fPedestal->GetPedestalN(n); }
    void SetPedestal( AliITSPedestalSSD* pedestal) {fPedestal=pedestal;}

    void AddGainP(Int_t c, Float_t n) { fGain->AddGainP(c,n);}       
    TArrayF GetGainP() {return fGain->GetGainP(); }
    Float_t GetGainP(Int_t n) {return fGain->GetGainP(n); }
    void AddGainN(Int_t c, Float_t n) { fGain->AddGainN(c,n);}
    TArrayF GetGainN() {return fGain->GetGainN(); }
    Float_t GetGainN(Int_t n) {return fGain->GetGainN(n); }
    void SetGain( AliITSGainSSD* gain) {fGain=gain;}

    void   SetBad() {
      fIsBad = kTRUE;
      for(Int_t i=0;i<fgkChipsPerModule;i++) fIsChipBad[i]=kTRUE;
    }
    virtual Bool_t IsBad() const { return fIsBad; }
    void   SetChipBad(Int_t nChip) {
      fIsChipBad[nChip] = kTRUE;
    }
    virtual Bool_t IsChipBad(Int_t nChip) const {
      return fIsChipBad[nChip];
    }
    Int_t ChipsPerModule() const{return fgkChipsPerModule;} // Number of chips/module
    Int_t ChannelsPerChip() const{ return fgkChannelsPerChip;}//Number of channels/chip

    TArrayI GetBadPChannelsList() { return fBadChannels->GetBadPChannelsList(); } const
    TArrayI GetBadNChannelsList() { return fBadChannels->GetBadNChannelsList(); } const
    void SetBadChannels( AliITSBadChannelsSSD* badchannels) {fBadChannels=badchannels;}

    void SetNoisePThreshold(Int_t threshold) { fNoisePThreshold = threshold;}
    void AddNoisyPChannel(Int_t c, Int_t n) { fNoisyPChannelsList.AddAt(n,c);}
    TArrayI GetNoisyPChannelsList() const {return fNoisyPChannelsList; }
    void SetNoiseNThreshold(Int_t threshold) { fNoiseNThreshold = threshold;}
    void AddNoisyNChannel(Int_t c, Int_t n) { fNoisyNChannelsList.AddAt(n,c);}
    TArrayI GetNoisyNChannelsList() const {return fNoisyNChannelsList; }

    void SetNDeadPChannelsList(Int_t n) { fDeadPChannelsList.Set(n); }
    void AddDeadPChannel(Int_t c, Int_t n) { fDeadPChannelsList.AddAt(n,c);}
    TArrayI GetDeadPChannelsList() const {return fDeadPChannelsList; }
    void SetNDeadNChannelsList(Int_t n) { fDeadNChannelsList.Set(n); }
    void AddDeadNChannel(Int_t c, Int_t n) { fDeadNChannelsList.AddAt(n,c);}
    TArrayI GetDeadNChannelsList() const {return fDeadNChannelsList; }
    //


    virtual  void   SetNDetParam(Int_t npar) {
	// set number of param
	fNPar=npar;
    }
    virtual  void   SetDetParam(Double_t *par);

    // Parameters options
    virtual Int_t   NDetParam() const {
	// number of param
	return fNPar;
    }
    virtual void    GetDetParam(Double_t *dpar) const; 

    virtual void    SetSigmaSpread(Double_t p1, Double_t p2) {
    	// Set sigmas of the charge spread function: Pside-Nside
    	// square of (microns)
    	fSigmaP=p1; fSigmaN=p2;
    }

    virtual void    SigmaSpread(Double_t &sP, Double_t &sN) const {
      // Get sigmas for the charge spread 
      sP=fSigmaP; sN=fSigmaN;
    }

    virtual void   SetThresholds(Double_t /* a */, Double_t /* b */)
      {NotImplemented("SetThresholds");}
    virtual void   Thresholds(Double_t & /* a */, Double_t & /* b */) const 
      {NotImplemented("Thresholds");}
  
    virtual void SetParamOptions(const char *opt1, const char *opt2) {((AliITSresponseSSD*)fResponse)->SetParamOptions(opt1,opt2);}
    virtual void GetParamOptions(char *opt1,char *opt2) const {((AliITSresponseSSD*)fResponse)->ParamOptions(opt1,opt2);}
    virtual void SetADCpereV(Double_t a=120./24888.9) {((AliITSresponseSSD*)fResponse)->SetADCpereV(a);}
    virtual Double_t GetDEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->DEvToADC(eV);}
    virtual Int_t IEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->IEvToADC(eV);} 

    virtual void SetKeVperADC(Double_t a=86.4/120.0) {((AliITSresponseSSD*)fResponse)->SetKeVperADC(a);}
    virtual Double_t ADCToKeV(Double_t adc) const {return ((AliITSresponseSSD*)fResponse)->ADCToKeV(adc);}

  virtual Double_t GetCouplingPR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPR();}
  virtual Double_t GetCouplingPL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPL();}
  virtual Double_t GetCouplingNR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNR();}
  virtual Double_t GetCouplingNL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNL();}
  virtual void SetCouplings(Double_t pr, Double_t pl, Double_t nr, Double_t nl) 
    { ((AliITSresponseSSD*)fResponse)->SetCouplings(pr,pl,nr,nl);}

  virtual Int_t GetZSThreshold() const {return ((AliITSresponseSSD*)fResponse)->GetZSThreshold();}
  virtual void SetZSThreshold(Int_t zsth) 
    { ((AliITSresponseSSD*)fResponse)->SetZSThreshold(zsth);}

protected:

    static const Int_t fgkChipsPerModule  = 12;    // Number of chips/module
    static const Int_t fgkChannelsPerChip = 128;   // Number of channels/chip

    static const Double_t fgkNoiseNDefault; // default for fNoiseN
    static const Double_t fgkNoisePDefault; // default for fNoiseP
    static const Int_t fgkNParDefault; // default for fNPar
    static const Double_t fgkSigmaPDefault; //default for fSigmaP
    static const Double_t fgkSigmaNDefault; //default for fSigmaP

    Int_t   fNPar;            // Number of detector param 
    Double_t *fDetPar;         //[fNPar] Array of parameters

    Double_t fNoiseP;          // Noise on Pside
    Double_t fNoiseN;          // Noise on Nside
    Double_t fSigmaP;          // Sigma charge spread on Pside
    Double_t fSigmaN;          // Sigma charge spread on Nside
    
    AliITSNoiseSSD *fNoise;
    AliITSPedestalSSD *fPedestal;
    AliITSGainSSD *fGain;
    AliITSBadChannelsSSD *fBadChannels;

    Bool_t   fIsBad;                         // module is dead or alive ?
    Bool_t   fIsChipBad[fgkChipsPerModule];  // chip is dead or alive ?

    TArrayF fGainP;           // Gain for P side channels
    TArrayF fGainN;           // Gain for N side channels

    TArrayF fNoisP;           // Noise for P side channels
    TArrayF fNoisN;           // Noise for N side channels

    TArrayF fPedP;
    TArrayF fPedN;

    Float_t fNoisePThreshold;     // need to decide if channel is noisy  
    TArrayI  fNoisyPChannelsList; // list of P side noisy channels
    Float_t fNoiseNThreshold;     // need to decide if channel is noisy  
    TArrayI  fNoisyNChannelsList; // list of N side noisy channels

    TArrayI  fDeadNChannelsList;  // list of P side dead channels
    TArrayI  fDeadPChannelsList;  // list of N side dead channels

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

    ClassDef(AliITSCalibrationSSD,3) //Response class for SSD
};
#endif
Commit	Line	Data
fcf95fc7	1	#ifndef ALIITSCALIBRATIONSSD_H
fcf95fc7	2	#define ALIITSCALIBRATIONSSD_H
b42cfa25	3	/* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
b42cfa25	4	* See cxx source for full Copyright notice */
fcf95fc7	5	#include "AliITSCalibration.h"
fb4dfab9	6	#include "AliITSNoiseSSD.h"
b42cfa25	7	#include "AliITSPedestalSSD.h"
fb4dfab9	8	#include "AliITSGainSSD.h"
fb4dfab9	9	#include "AliITSBadChannelsSSD.h"
fcf95fc7	10	#include "AliITSresponseSSD.h"
569a17d8	11	#include "TArrayF.h"
14a74335	12	#include "TArrayI.h"
fcf95fc7	13
b42cfa25	14	/* $Id$ */
fcf95fc7	15	//////////////////////////////////////////////
	16	// Response class for SSD //
	17	// //
	18	//////////////////////////////////////////////
	19	class AliITSCalibrationSSD : public AliITSCalibration {
	20
	21	public:
	22	AliITSCalibrationSSD();
	23	AliITSCalibrationSSD(const char *dataType);
	24	virtual ~AliITSCalibrationSSD();
14a74335	25
fcf95fc7	26	virtual void SetNoiseParam(Double_t np, Double_t nn) {
14a74335	27	// set noise par
14a74335	28	fNoiseP=np; fNoiseN=nn;
fcf95fc7	29	}
14a74335	30
fcf95fc7	31	virtual void GetNoiseParam(Double_t &np, Double_t &nn) const {
14a74335	32	// get noise par
14a74335	33	np=fNoiseP; nn=fNoiseN;
fcf95fc7	34	}
fcf95fc7	35
fb4dfab9	36	void AddNoiseP(Int_t c, Float_t n) { fNoise->AddNoiseP(c,n);}
	37	TArrayF GetNoiseP() {return fNoise->GetNoiseP(); }
	38	Float_t GetNoiseP(Int_t n) {return fNoise->GetNoiseP(n); }
	39	void AddNoiseN(Int_t c, Float_t n) { fNoise->AddNoiseN(c,n);}
	40	TArrayF GetNoiseN() {return fNoise->GetNoiseN(); }
	41	Float_t GetNoiseN(Int_t n) {return fNoise->GetNoiseN(n); }
	42	void SetNoise( AliITSNoiseSSD* noise) {fNoise=noise;}
	43
b42cfa25	44	void AddPedestalP(Int_t c, Float_t n) { fPedestal->AddPedestalP(c,n);}
	45	TArrayF GetPedestalP() {return fPedestal->GetPedestalP(); }
	46	Float_t GetPedestalP(Int_t n) {return fPedestal->GetPedestalP(n); }
	47	void AddPedestalN(Int_t c, Float_t n) { fPedestal->AddPedestalN(c,n);}
	48	TArrayF GetPedestalN() {return fPedestal->GetPedestalN(); }
	49	Float_t GetPedestalN(Int_t n) {return fPedestal->GetPedestalN(n); }
	50	void SetPedestal( AliITSPedestalSSD* pedestal) {fPedestal=pedestal;}
	51
fb4dfab9	52	void AddGainP(Int_t c, Float_t n) { fGain->AddGainP(c,n);}
	53	TArrayF GetGainP() {return fGain->GetGainP(); }
	54	Float_t GetGainP(Int_t n) {return fGain->GetGainP(n); }
	55	void AddGainN(Int_t c, Float_t n) { fGain->AddGainN(c,n);}
	56	TArrayF GetGainN() {return fGain->GetGainN(); }
	57	Float_t GetGainN(Int_t n) {return fGain->GetGainN(n); }
	58	void SetGain( AliITSGainSSD* gain) {fGain=gain;}
	59
df1d5598	60	void SetBad() {
	61	fIsBad = kTRUE;
	62	for(Int_t i=0;i<fgkChipsPerModule;i++) fIsChipBad[i]=kTRUE;
	63	}
	64	virtual Bool_t IsBad() const { return fIsBad; }
	65	void SetChipBad(Int_t nChip) {
	66	fIsChipBad[nChip] = kTRUE;
	67	}
	68	virtual Bool_t IsChipBad(Int_t nChip) const {
	69	return fIsChipBad[nChip];
	70	}
	71	Int_t ChipsPerModule() const{return fgkChipsPerModule;} // Number of chips/module
	72	Int_t ChannelsPerChip() const{ return fgkChannelsPerChip;}//Number of channels/chip
	73
ad2b76cb	74	TArrayI GetBadPChannelsList() { return fBadChannels->GetBadPChannelsList(); } const
ad2b76cb	75	TArrayI GetBadNChannelsList() { return fBadChannels->GetBadNChannelsList(); } const
fb4dfab9	76	void SetBadChannels( AliITSBadChannelsSSD* badchannels) {fBadChannels=badchannels;}
14a74335	77
	78	void SetNoisePThreshold(Int_t threshold) { fNoisePThreshold = threshold;}
	79	void AddNoisyPChannel(Int_t c, Int_t n) { fNoisyPChannelsList.AddAt(n,c);}
	80	TArrayI GetNoisyPChannelsList() const {return fNoisyPChannelsList; }
	81	void SetNoiseNThreshold(Int_t threshold) { fNoiseNThreshold = threshold;}
	82	void AddNoisyNChannel(Int_t c, Int_t n) { fNoisyNChannelsList.AddAt(n,c);}
	83	TArrayI GetNoisyNChannelsList() const {return fNoisyNChannelsList; }
	84
	85	void SetNDeadPChannelsList(Int_t n) { fDeadPChannelsList.Set(n); }
	86	void AddDeadPChannel(Int_t c, Int_t n) { fDeadPChannelsList.AddAt(n,c);}
	87	TArrayI GetDeadPChannelsList() const {return fDeadPChannelsList; }
	88	void SetNDeadNChannelsList(Int_t n) { fDeadNChannelsList.Set(n); }
	89	void AddDeadNChannel(Int_t c, Int_t n) { fDeadNChannelsList.AddAt(n,c);}
	90	TArrayI GetDeadNChannelsList() const {return fDeadNChannelsList; }
	91	//
	92
	93
fcf95fc7	94	virtual void SetNDetParam(Int_t npar) {
	95	// set number of param
	96	fNPar=npar;
	97	}
	98	virtual void SetDetParam(Double_t *par);
	99
	100	// Parameters options
	101	virtual Int_t NDetParam() const {
	102	// number of param
	103	return fNPar;
	104	}
	105	virtual void GetDetParam(Double_t *dpar) const;
	106
	107	virtual void SetSigmaSpread(Double_t p1, Double_t p2) {
14a74335	108	// Set sigmas of the charge spread function: Pside-Nside
	109	// square of (microns)
	110	fSigmaP=p1; fSigmaN=p2;
fcf95fc7	111	}
14a74335	112
fcf95fc7	113	virtual void SigmaSpread(Double_t &sP, Double_t &sN) const {
14a74335	114	// Get sigmas for the charge spread
14a74335	115	sP=fSigmaP; sN=fSigmaN;
fcf95fc7	116	}
	117
	118	virtual void SetThresholds(Double_t /* a /, Double_t / b */)
	119	{NotImplemented("SetThresholds");}
	120	virtual void Thresholds(Double_t & /* a /, Double_t & / b */) const
	121	{NotImplemented("Thresholds");}
	122
	123	virtual void SetParamOptions(const char opt1, const char opt2) {((AliITSresponseSSD*)fResponse)->SetParamOptions(opt1,opt2);}
	124	virtual void GetParamOptions(char opt1,char opt2) const {((AliITSresponseSSD*)fResponse)->ParamOptions(opt1,opt2);}
9acf2ecc	125	virtual void SetADCpereV(Double_t a=120./24888.9) {((AliITSresponseSSD*)fResponse)->SetADCpereV(a);}
fcf95fc7	126	virtual Double_t GetDEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->DEvToADC(eV);}
	127	virtual Int_t IEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->IEvToADC(eV);}
	128
9acf2ecc	129	virtual void SetKeVperADC(Double_t a=86.4/120.0) {((AliITSresponseSSD*)fResponse)->SetKeVperADC(a);}
	130	virtual Double_t ADCToKeV(Double_t adc) const {return ((AliITSresponseSSD*)fResponse)->ADCToKeV(adc);}
	131
14a74335	132	virtual Double_t GetCouplingPR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPR();}
	133	virtual Double_t GetCouplingPL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPL();}
	134	virtual Double_t GetCouplingNR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNR();}
	135	virtual Double_t GetCouplingNL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNL();}
	136	virtual void SetCouplings(Double_t pr, Double_t pl, Double_t nr, Double_t nl)
	137	{ ((AliITSresponseSSD*)fResponse)->SetCouplings(pr,pl,nr,nl);}
	138
	139	virtual Int_t GetZSThreshold() const {return ((AliITSresponseSSD*)fResponse)->GetZSThreshold();}
	140	virtual void SetZSThreshold(Int_t zsth)
	141	{ ((AliITSresponseSSD*)fResponse)->SetZSThreshold(zsth);}
	142
	143	protected:
df1d5598	144
	145	static const Int_t fgkChipsPerModule = 12; // Number of chips/module
	146	static const Int_t fgkChannelsPerChip = 128; // Number of channels/chip
	147
fcf95fc7	148	static const Double_t fgkNoiseNDefault; // default for fNoiseN
	149	static const Double_t fgkNoisePDefault; // default for fNoiseP
	150	static const Int_t fgkNParDefault; // default for fNPar
	151	static const Double_t fgkSigmaPDefault; //default for fSigmaP
	152	static const Double_t fgkSigmaNDefault; //default for fSigmaP
14a74335	153
fcf95fc7	154	Int_t fNPar; // Number of detector param
	155	Double_t *fDetPar; //[fNPar] Array of parameters
	156
	157	Double_t fNoiseP; // Noise on Pside
	158	Double_t fNoiseN; // Noise on Nside
fcf95fc7	159	Double_t fSigmaP; // Sigma charge spread on Pside
	160	Double_t fSigmaN; // Sigma charge spread on Nside
	161
fb4dfab9	162	AliITSNoiseSSD *fNoise;
b42cfa25	163	AliITSPedestalSSD *fPedestal;
fb4dfab9	164	AliITSGainSSD *fGain;
	165	AliITSBadChannelsSSD *fBadChannels;
	166
df1d5598	167	Bool_t fIsBad; // module is dead or alive ?
	168	Bool_t fIsChipBad[fgkChipsPerModule]; // chip is dead or alive ?
	169
569a17d8	170	TArrayF fGainP; // Gain for P side channels
569a17d8	171	TArrayF fGainN; // Gain for N side channels
14a74335	172
569a17d8	173	TArrayF fNoisP; // Noise for P side channels
569a17d8	174	TArrayF fNoisN; // Noise for N side channels
14a74335	175
b42cfa25	176	TArrayF fPedP;
	177	TArrayF fPedN;
	178
569a17d8	179	Float_t fNoisePThreshold; // need to decide if channel is noisy
14a74335	180	TArrayI fNoisyPChannelsList; // list of P side noisy channels
569a17d8	181	Float_t fNoiseNThreshold; // need to decide if channel is noisy
14a74335	182	TArrayI fNoisyNChannelsList; // list of N side noisy channels
	183
	184	TArrayI fDeadNChannelsList; // list of P side dead channels
	185	TArrayI fDeadPChannelsList; // list of N side dead channels
	186
fcf95fc7	187	private:
	188	AliITSCalibrationSSD(const AliITSCalibrationSSD &source); // copy constructor
	189	AliITSCalibrationSSD& operator=(const AliITSCalibrationSSD &source); // ass. op.
	190
df1d5598	191	ClassDef(AliITSCalibrationSSD,3) //Response class for SSD
fcf95fc7	192	};
fcf95fc7	193	#endif