[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$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  TRD simulation - multimodule (regular rad.)                              //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

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