[u/mrichter/AliRoot.git] / TRD / AliTRDsim.h

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

/* $Id$ */

#include <TObject.h>

class TH1D;

class AliModule;

class AliTRDsim : public TObject {

 public:

  AliTRDsim();
  AliTRDsim(const AliTRDsim &s);
  AliTRDsim(AliModule *mod, Int_t foil, Int_t gap);
  virtual ~AliTRDsim();
  AliTRDsim &operator=(const AliTRDsim &s);

  virtual void          Copy(TObject &s);
  virtual void          Init();
  virtual Int_t         CreatePhotons(Int_t pdg, Float_t p
                                    , Int_t &nPhoton, Float_t *ePhoton);
  virtual Int_t         TrPhotons(Double_t gamma, Int_t &nPhoton, Float_t *ePhoton);
  virtual Double_t      Sigma(Double_t energykeV);
  virtual Double_t      Interpolate(Double_t energyMeV
                                  , Double_t *en, Double_t *mu, Int_t n);
  virtual Int_t         Locate(Double_t *xv, Int_t n, Double_t xval
                             , Int_t &kl, Double_t &dx);
  virtual Double_t      Omega(Float_t rho, Float_t z, Float_t a) 
                             { return (28.8 * TMath::Sqrt(rho * z / a)); };

          void          SetNFoils(Int_t n)      { fNFoils    = n;       };
          void          SetFoilThick(Float_t t) { fFoilThick = t;
                           SetSigma();                                  };
          void          SetGapThick(Float_t t)  { fGapThick  = t;
                           SetSigma();                                  };
          void          SetFoilDens(Float_t d)  { fFoilDens  = d; 
                           fFoilOmega = Omega(fFoilDens,fFoilZ,fFoilA);
                           SetSigma();                                  };
          void          SetFoilZ(Float_t z)     { fFoilZ     = z; 
                           fFoilOmega = Omega(fFoilDens,fFoilZ,fFoilA); };
          void          SetFoilA(Float_t a)     { fFoilA     = a; 
                           fFoilOmega = Omega(fFoilDens,fFoilZ,fFoilA); };
          void          SetGapDens(Float_t d)   { fGapDens   = d;
                           fGapOmega  = Omega(fGapDens ,fGapZ ,fGapA );
                           SetSigma();                                  };
          void          SetGapZ(Float_t z)      { fGapZ      = z;
                           fGapOmega  = Omega(fGapDens ,fGapZ ,fGapA ); };
          void          SetGapA(Float_t a)      { fGapA      = a;
                           fGapOmega  = Omega(fGapDens ,fGapZ ,fGapA ); };
          void          SetTemp(Float_t t)      { fTemp      = t; 
                           SetSigma();                                  };
          void          SetSigma();

  virtual Double_t      GetMuPo(Double_t energyMeV);
  virtual Double_t      GetMuCO(Double_t energyMeV);
  virtual Double_t      GetMuXe(Double_t energyMeV);
  virtual Double_t      GetMuBu(Double_t energyMeV);
  virtual Double_t      GetMuMy(Double_t energyMeV);
  virtual Double_t      GetMuN2(Double_t energyMeV);
  virtual Double_t      GetMuO2(Double_t energyMeV);
  virtual Double_t      GetMuHe(Double_t energyMeV);

          Int_t         GetNFoils() const       { return fNFoils;    };
          Float_t       GetFoilThick() const    { return fFoilThick; };
          Float_t       GetGapThick() const     { return fGapThick;  };
          Float_t       GetFoilDens() const     { return fFoilDens;  };
          Float_t       GetGapDens() const      { return fGapDens;   };
          Double_t      GetFoilOmega() const    { return fFoilOmega; };
          Double_t      GetGapOmega() const     { return fGapOmega;  };
          Float_t       GetTemp() const         { return fTemp;      };
          TH1D         *GetSpectrum() const     { return fSpectrum;  };

 protected:

  Int_t     fNFoils;               // Number of foils in the radiator stack
  Float_t   fFoilThick;            // Thickness of the foils (cm)
  Float_t   fGapThick;             // Thickness of the gaps between the foils (cm)

  Float_t   fFoilDens;             // Density of the radiator foils (g/cm^3) 
  Float_t   fGapDens;              // Density of the gas in the radiator gaps (g/cm^3)

  Double_t  fFoilOmega;            // Plasma frequency of the radiator foils
  Double_t  fGapOmega;             // Plasma frequency of the gas in the radiator gaps

  Float_t   fFoilZ;                // Z of the foil material
  Float_t   fGapZ;                 // Z of the gas in the gaps

  Float_t   fFoilA;                // A of the foil material
  Float_t   fGapA;                 // A of the gas in the gaps

  Float_t   fTemp;                 // Temperature of the radiator gas (Kelvin)

  Int_t     fSpNBins;              // Number of bins of the TR spectrum
  Float_t   fSpRange;              // Range of the TR spectrum
  Float_t   fSpBinWidth;           // Bin width of the TR spectrum
  Float_t   fSpLower;              // Lower border of the TR spectrum
  Float_t   fSpUpper;              // Upper border of the TR spectrum

  Double_t *fSigma;                //[fSpNBins] Array of sigma values

  TH1D     *fSpectrum;             //!TR photon energy spectrum

  ClassDef(AliTRDsim,1)            // Simulates TR photons

};
#endif
Commit	Line	Data
46d29e70	1	#ifndef ALITRDSIM_H
	2	#define ALITRDSIM_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7
	8	#include <TObject.h>
	9
	10	class TH1D;
	11
	12	class AliModule;
	13
	14	class AliTRDsim : public TObject {
	15
	16	public:
	17
	18	AliTRDsim();
	19	AliTRDsim(const AliTRDsim &s);
	20	AliTRDsim(AliModule *mod, Int_t foil, Int_t gap);
	21	virtual ~AliTRDsim();
	22	AliTRDsim &operator=(const AliTRDsim &s);
	23
	24	virtual void Copy(TObject &s);
	25	virtual void Init();
	26	virtual Int_t CreatePhotons(Int_t pdg, Float_t p
	27	, Int_t &nPhoton, Float_t *ePhoton);
	28	virtual Int_t TrPhotons(Double_t gamma, Int_t &nPhoton, Float_t *ePhoton);
	29	virtual Double_t Sigma(Double_t energykeV);
	30	virtual Double_t Interpolate(Double_t energyMeV
	31	, Double_t en, Double_t mu, Int_t n);
	32	virtual Int_t Locate(Double_t *xv, Int_t n, Double_t xval
	33	, Int_t &kl, Double_t &dx);
	34	virtual Double_t Omega(Float_t rho, Float_t z, Float_t a)
	35	{ return (28.8 * TMath::Sqrt(rho * z / a)); };
	36
	37	void SetNFoils(Int_t n) { fNFoils = n; };
	38	void SetFoilThick(Float_t t) { fFoilThick = t;
	39	SetSigma(); };
	40	void SetGapThick(Float_t t) { fGapThick = t;
	41	SetSigma(); };
	42	void SetFoilDens(Float_t d) { fFoilDens = d;
	43	fFoilOmega = Omega(fFoilDens,fFoilZ,fFoilA);
	44	SetSigma(); };
	45	void SetFoilZ(Float_t z) { fFoilZ = z;
	46	fFoilOmega = Omega(fFoilDens,fFoilZ,fFoilA); };
	47	void SetFoilA(Float_t a) { fFoilA = a;
	48	fFoilOmega = Omega(fFoilDens,fFoilZ,fFoilA); };
	49	void SetGapDens(Float_t d) { fGapDens = d;
	50	fGapOmega = Omega(fGapDens ,fGapZ ,fGapA );
	51	SetSigma(); };
	52	void SetGapZ(Float_t z) { fGapZ = z;
	53	fGapOmega = Omega(fGapDens ,fGapZ ,fGapA ); };
	54	void SetGapA(Float_t a) { fGapA = a;
	55	fGapOmega = Omega(fGapDens ,fGapZ ,fGapA ); };
	56	void SetTemp(Float_t t) { fTemp = t;
	57	SetSigma(); };
	58	void SetSigma();
	59
	60	virtual Double_t GetMuPo(Double_t energyMeV);
	61	virtual Double_t GetMuCO(Double_t energyMeV);
	62	virtual Double_t GetMuXe(Double_t energyMeV);
	63	virtual Double_t GetMuBu(Double_t energyMeV);
	64	virtual Double_t GetMuMy(Double_t energyMeV);
65	virtual Double_t GetMuN2(Double_t energyMeV);
66	virtual Double_t GetMuO2(Double_t energyMeV);
67	virtual Double_t GetMuHe(Double_t energyMeV);
68
69	Int_t GetNFoils() const { return fNFoils; };
70	Float_t GetFoilThick() const { return fFoilThick; };
71	Float_t GetGapThick() const { return fGapThick; };
72	Float_t GetFoilDens() const { return fFoilDens; };
73	Float_t GetGapDens() const { return fGapDens; };
74	Double_t GetFoilOmega() const { return fFoilOmega; };
75	Double_t GetGapOmega() const { return fGapOmega; };
76	Float_t GetTemp() const { return fTemp; };
77	TH1D *GetSpectrum() const { return fSpectrum; };
78
79	protected:
80
81	Int_t fNFoils; // Number of foils in the radiator stack
82	Float_t fFoilThick; // Thickness of the foils (cm)
83	Float_t fGapThick; // Thickness of the gaps between the foils (cm)
84
85	Float_t fFoilDens; // Density of the radiator foils (g/cm^3)
86	Float_t fGapDens; // Density of the gas in the radiator gaps (g/cm^3)
87
88	Double_t fFoilOmega; // Plasma frequency of the radiator foils
89	Double_t fGapOmega; // Plasma frequency of the gas in the radiator gaps
90
91	Float_t fFoilZ; // Z of the foil material
92	Float_t fGapZ; // Z of the gas in the gaps
93
94	Float_t fFoilA; // A of the foil material
95	Float_t fGapA; // A of the gas in the gaps
96
97	Float_t fTemp; // Temperature of the radiator gas (Kelvin)
98
99	Int_t fSpNBins; // Number of bins of the TR spectrum
100	Float_t fSpRange; // Range of the TR spectrum
101	Float_t fSpBinWidth; // Bin width of the TR spectrum
102	Float_t fSpLower; // Lower border of the TR spectrum
103	Float_t fSpUpper; // Upper border of the TR spectrum
104
d1b06c24	105	Double_t *fSigma; //[fSpNBins] Array of sigma values
46d29e70	106
a2b90f83	107	TH1D *fSpectrum; //!TR photon energy spectrum
46d29e70	108
	109	ClassDef(AliTRDsim,1) // Simulates TR photons
	110
	111	};
	112	#endif