Logo ROOT   6.16/01
Reference Guide
TGeoMaterial.h
Go to the documentation of this file.
1// @(#)root/geom:$Id$
2// Author: Andrei Gheata 25/10/01
3
4/*************************************************************************
5 * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6 * All rights reserved. *
7 * *
8 * For the licensing terms see $ROOTSYS/LICENSE. *
9 * For the list of contributors see $ROOTSYS/README/CREDITS. *
10 *************************************************************************/
11
12#ifndef ROOT_TGeoMaterial
13#define ROOT_TGeoMaterial
14
15#include "TNamed.h"
16
17#include "TAttFill.h"
18
19#include "TGeoElement.h"
20
21
22// forward declarations
23class TGeoExtension;
24
25// Some units used in G4
26static const Double_t STP_temperature = 273.15; // [K]
27static const Double_t STP_pressure = 6.32420e+8; // [MeV/mm3]
28
29class TGeoMaterial : public TNamed,
30 public TAttFill
31{
32public:
36 };
42 };
43
44protected:
45 Int_t fIndex; // material index
46 Double_t fA; // A of material
47 Double_t fZ; // Z of material
48 Double_t fDensity; // density of material
49 Double_t fRadLen; // radiation length
50 Double_t fIntLen; // interaction length
51 Double_t fTemperature; // temperature
52 Double_t fPressure; // pressure
53 EGeoMaterialState fState; // material state
54 TObject *fShader; // shader with optical properties
55 TObject *fCerenkov; // pointer to class with Cerenkov properties
56 TGeoElement *fElement; // pointer to element composing the material
57 TGeoExtension *fUserExtension; //! Transient user-defined extension to materials
58 TGeoExtension *fFWExtension; //! Transient framework-defined extension to materials
59
60// methods
63
64
65public:
66 // constructors
68 TGeoMaterial(const char *name);
69 TGeoMaterial(const char *name, Double_t a, Double_t z,
70 Double_t rho, Double_t radlen=0, Double_t intlen=0);
71 TGeoMaterial(const char *name, Double_t a, Double_t z, Double_t rho,
73 TGeoMaterial(const char *name, TGeoElement *elem, Double_t rho);
74
75 // destructor
76 virtual ~TGeoMaterial();
77 // methods
78 static Double_t Coulomb(Double_t z);
79 // radioactive mixture evolution
80 virtual TGeoMaterial *DecayMaterial(Double_t time, Double_t precision=0.001);
81 virtual void FillMaterialEvolution(TObjArray *population, Double_t precision=0.001);
82 // getters & setters
83 virtual Int_t GetByteCount() const {return sizeof(*this);}
84 virtual Double_t GetA() const {return fA;}
85 virtual Double_t GetZ() const {return fZ;}
86 virtual Int_t GetDefaultColor() const;
87 virtual Double_t GetDensity() const {return fDensity;}
88 virtual Int_t GetNelements() const {return 1;}
89 virtual TGeoElement *GetElement(Int_t i=0) const;
90 virtual void GetElementProp(Double_t &a, Double_t &z, Double_t &w, Int_t i=0);
92 char *GetPointerName() const;
93 virtual Double_t GetRadLen() const {return fRadLen;}
94 virtual Double_t GetIntLen() const {return fIntLen;}
96 virtual TObject *GetCerenkovProperties() const {return fCerenkov;}
97 Char_t GetTransparency() const {return (fFillStyle<3000 || fFillStyle>3100)?0:Char_t(fFillStyle-3000);}
99 Double_t GetPressure() const {return fPressure;}
101 virtual Double_t GetSpecificActivity(Int_t) const {return 0.;}
106 virtual Bool_t IsEq(const TGeoMaterial *other) const;
108 virtual Bool_t IsMixture() const {return kFALSE;}
109 virtual void Print(const Option_t *option="") const;
110 virtual void SavePrimitive(std::ostream &out, Option_t *option = "");
111 virtual void SetA(Double_t a) {fA = a; SetRadLen(0);}
112 virtual void SetZ(Double_t z) {fZ = z; SetRadLen(0);}
113 virtual void SetDensity(Double_t density) {fDensity = density; SetRadLen(0);}
114 void SetIndex(Int_t index) {fIndex=index;}
115 virtual void SetCerenkovProperties(TObject* cerenkov) {fCerenkov = cerenkov;}
116 void SetRadLen(Double_t radlen, Double_t intlen=0.);
118 void SetTransparency(Char_t transparency=0) {fFillStyle = 3000+transparency;}
119 void SetTemperature(Double_t temperature) {fTemperature = temperature;}
120 void SetPressure(Double_t pressure) {fPressure = pressure;}
121 void SetState(EGeoMaterialState state) {fState = state;}
123 void SetFWExtension(TGeoExtension *ext);
125
126
127
128 ClassDef(TGeoMaterial, 5) // base material class
129
130//***** Need to add classes and globals to LinkDef.h *****
131};
132
133
135{
136protected :
137// data members
138 Int_t fNelements; // number of elements
139 Double_t *fZmixture; // [fNelements] array of Z of the elements
140 Double_t *fAmixture; // [fNelements] array of A of the elements
141 Double_t *fWeights; // [fNelements] array of relative proportions by mass
142 Int_t *fNatoms; // [fNelements] array of numbers of atoms
143 Double_t *fVecNbOfAtomsPerVolume; //[fNelements] array of numbers of atoms per unit volume
144 TObjArray *fElements; // array of elements composing the mixture
145// methods
146 TGeoMixture(const TGeoMixture&); // Not implemented
147 TGeoMixture& operator=(const TGeoMixture&); // Not implemented
148 void AverageProperties();
149
150public:
151 // constructors
152 TGeoMixture();
153 TGeoMixture(const char *name, Int_t nel, Double_t rho=-1);
154 // destructor
155 virtual ~TGeoMixture();
156 // methods for adding elements
157 void AddElement(Double_t a, Double_t z, Double_t weight);
158 void AddElement(TGeoMaterial *mat, Double_t weight);
159 void AddElement(TGeoElement *elem, Double_t weight);
160 void AddElement(TGeoElement *elem, Int_t natoms);
161 // backward compatibility for defining elements
162 void DefineElement(Int_t iel, Double_t a, Double_t z, Double_t weight);
163 void DefineElement(Int_t iel, TGeoElement *elem, Double_t weight);
164 void DefineElement(Int_t iel, Int_t z, Int_t natoms);
165 // radioactive mixture evolution
166 virtual TGeoMaterial *DecayMaterial(Double_t time, Double_t precision=0.001);
167 virtual void FillMaterialEvolution(TObjArray *population, Double_t precision=0.001);
168 // getters
169 virtual Int_t GetByteCount() const {return 48+12*fNelements;}
170 virtual TGeoElement *GetElement(Int_t i=0) const;
171 virtual void GetElementProp(Double_t &a, Double_t &z, Double_t &w, Int_t i=0) {a=fAmixture[i]; z=fZmixture[i]; w=fWeights[i];}
172 virtual Int_t GetNelements() const {return fNelements;}
173 Double_t *GetZmixt() const {return fZmixture;}
174 Double_t *GetAmixt() const {return fAmixture;}
175 Double_t *GetWmixt() const {return fWeights;}
176 Int_t *GetNmixt() const {return fNatoms;}
177 virtual Double_t GetSpecificActivity(Int_t i=-1) const;
178 // utilities
179 virtual Bool_t IsEq(const TGeoMaterial *other) const;
180 virtual Bool_t IsMixture() const {return kTRUE;}
181 virtual void Print(const Option_t *option="") const;
182 virtual void SavePrimitive(std::ostream &out, Option_t *option = "");
183 virtual void SetA(Double_t a) {fA = a;}
184 virtual void SetZ(Double_t z) {fZ = z;}
185 virtual void SetDensity(Double_t density) {fDensity = density; AverageProperties();}
189
190 ClassDef(TGeoMixture, 3) // material mixtures
191};
192
194 {return AddElement(a,z,weight);}
196 {return AddElement(elem,weight);}
197
198
199#endif
200
int Int_t
Definition: RtypesCore.h:41
char Char_t
Definition: RtypesCore.h:29
const Bool_t kFALSE
Definition: RtypesCore.h:88
bool Bool_t
Definition: RtypesCore.h:59
double Double_t
Definition: RtypesCore.h:55
const Bool_t kTRUE
Definition: RtypesCore.h:87
const char Option_t
Definition: RtypesCore.h:62
#define ClassDef(name, id)
Definition: Rtypes.h:324
#define BIT(n)
Definition: Rtypes.h:82
static const Double_t STP_temperature
Definition: TGeoMaterial.h:26
static const Double_t STP_pressure
Definition: TGeoMaterial.h:27
Fill Area Attributes class.
Definition: TAttFill.h:19
Style_t fFillStyle
Fill area style.
Definition: TAttFill.h:23
Base class for chemical elements.
Definition: TGeoElement.h:37
ABC for user objects attached to TGeoVolume or TGeoNode.
Definition: TGeoExtension.h:20
Base class describing materials.
Definition: TGeoMaterial.h:31
virtual ~TGeoMaterial()
Destructor.
virtual Double_t GetIntLen() const
Definition: TGeoMaterial.h:94
virtual void SetDensity(Double_t density)
Definition: TGeoMaterial.h:113
virtual void SetZ(Double_t z)
Definition: TGeoMaterial.h:112
void SetUserExtension(TGeoExtension *ext)
Connect user-defined extension to the material.
TGeoExtension * GetUserExtension() const
Definition: TGeoMaterial.h:102
TGeoExtension * GetFWExtension() const
Definition: TGeoMaterial.h:103
char * GetPointerName() const
Provide a pointer name containing uid.
virtual TObject * GetCerenkovProperties() const
Definition: TGeoMaterial.h:96
EGeoMaterialState fState
Definition: TGeoMaterial.h:53
virtual void SetCerenkovProperties(TObject *cerenkov)
Definition: TGeoMaterial.h:115
Bool_t IsUsed() const
Definition: TGeoMaterial.h:107
static Double_t ScreenFactor(Double_t z)
static function Compute screening factor for pair production and Bremsstrahlung REFERENCE : EGS MANUA...
virtual Bool_t IsMixture() const
Definition: TGeoMaterial.h:108
Double_t fPressure
Definition: TGeoMaterial.h:52
virtual Int_t GetByteCount() const
Definition: TGeoMaterial.h:83
void SetFWExtension(TGeoExtension *ext)
Connect framework defined extension to the material.
Double_t fTemperature
Definition: TGeoMaterial.h:51
virtual void GetElementProp(Double_t &a, Double_t &z, Double_t &w, Int_t i=0)
Single interface to get element properties.
virtual void Print(const Option_t *option="") const
print characteristics of this material
virtual TGeoMaterial * DecayMaterial(Double_t time, Double_t precision=0.001)
Create the material representing the decay product of this material at a given time.
void SetState(EGeoMaterialState state)
Definition: TGeoMaterial.h:121
EGeoMaterialState GetState() const
Definition: TGeoMaterial.h:100
virtual void SavePrimitive(std::ostream &out, Option_t *option="")
Save a primitive as a C++ statement(s) on output stream "out".
Double_t GetTemperature() const
Definition: TGeoMaterial.h:98
virtual Double_t GetSpecificActivity(Int_t) const
Definition: TGeoMaterial.h:101
virtual Int_t GetNelements() const
Definition: TGeoMaterial.h:88
Double_t fZ
Definition: TGeoMaterial.h:47
void SetRadLen(Double_t radlen, Double_t intlen=0.)
Set radiation/absorption lengths.
virtual Bool_t IsEq(const TGeoMaterial *other) const
return true if the other material has the same physical properties
TGeoElement * GetBaseElement() const
Definition: TGeoMaterial.h:91
Double_t fA
Definition: TGeoMaterial.h:46
TObject * fCerenkov
Definition: TGeoMaterial.h:55
void SetIndex(Int_t index)
Definition: TGeoMaterial.h:114
Double_t fDensity
Definition: TGeoMaterial.h:48
void SetUsed(Bool_t flag=kTRUE)
Definition: TGeoMaterial.h:117
virtual void FillMaterialEvolution(TObjArray *population, Double_t precision=0.001)
Fills a user array with all the elements deriving from the possible decay of the top element composin...
Int_t GetIndex()
Retrieve material index in the list of materials.
Double_t fIntLen
Definition: TGeoMaterial.h:50
Double_t GetPressure() const
Definition: TGeoMaterial.h:99
TGeoExtension * fUserExtension
Definition: TGeoMaterial.h:57
Char_t GetTransparency() const
Definition: TGeoMaterial.h:97
void SetPressure(Double_t pressure)
Definition: TGeoMaterial.h:120
TGeoElement * fElement
Definition: TGeoMaterial.h:56
void SetTransparency(Char_t transparency=0)
Definition: TGeoMaterial.h:118
TGeoMaterial & operator=(const TGeoMaterial &)
assignment operator
TObject * fShader
Definition: TGeoMaterial.h:54
TGeoExtension * GrabUserExtension() const
Get a copy of the user extension pointer.
Double_t fRadLen
Definition: TGeoMaterial.h:49
virtual TGeoElement * GetElement(Int_t i=0) const
Get a pointer to the element this material is made of.
TGeoExtension * GrabFWExtension() const
Get a copy of the framework extension pointer.
virtual Int_t GetDefaultColor() const
Get some default color related to this material.
TGeoMaterial()
Default constructor.
static Double_t Coulomb(Double_t z)
static function Compute Coulomb correction for pair production and Brem REFERENCE : EGS MANUAL SLAC 2...
virtual void SetA(Double_t a)
Definition: TGeoMaterial.h:111
virtual Double_t GetRadLen() const
Definition: TGeoMaterial.h:93
virtual Double_t GetA() const
Definition: TGeoMaterial.h:84
TGeoExtension * fFWExtension
Transient user-defined extension to materials.
Definition: TGeoMaterial.h:58
virtual Double_t GetDensity() const
Definition: TGeoMaterial.h:87
virtual Double_t GetZ() const
Definition: TGeoMaterial.h:85
void SetTemperature(Double_t temperature)
Definition: TGeoMaterial.h:119
Mixtures of elements.
Definition: TGeoMaterial.h:135
Int_t * GetNmixt() const
Definition: TGeoMaterial.h:176
TObjArray * fElements
Definition: TGeoMaterial.h:144
virtual TGeoElement * GetElement(Int_t i=0) const
Retrieve the pointer to the element corresponding to component I.
void ComputeNuclearInterLength()
Compute Nuclear Interaction Length based on Geant4 formula.
Double_t * GetZmixt() const
Definition: TGeoMaterial.h:173
void AddElement(Double_t a, Double_t z, Double_t weight)
add an element to the mixture using fraction by weight Check if the element is already defined
TGeoMixture & operator=(const TGeoMixture &)
assignment operator
Double_t * fZmixture
Definition: TGeoMaterial.h:139
virtual Double_t GetSpecificActivity(Int_t i=-1) const
Get specific activity (in Bq/gram) for the whole mixture (no argument) or for a given component.
TGeoMixture()
Default constructor.
virtual Int_t GetByteCount() const
Definition: TGeoMaterial.h:169
virtual ~TGeoMixture()
Destructor.
void ComputeDerivedQuantities()
Compute Derived Quantities as in Geant4.
virtual void Print(const Option_t *option="") const
print characteristics of this material
virtual void SetZ(Double_t z)
Definition: TGeoMaterial.h:184
void ComputeRadiationLength()
Compute Radiation Length based on Geant4 formula.
Double_t * fVecNbOfAtomsPerVolume
Definition: TGeoMaterial.h:143
virtual void SetDensity(Double_t density)
Definition: TGeoMaterial.h:185
virtual void FillMaterialEvolution(TObjArray *population, Double_t precision=0.001)
Fills a user array with all the elements deriving from the possible decay of the top elements composi...
Double_t * fAmixture
Definition: TGeoMaterial.h:140
void AverageProperties()
Compute effective A/Z and radiation length.
Int_t fNelements
Definition: TGeoMaterial.h:138
Double_t * GetWmixt() const
Definition: TGeoMaterial.h:175
virtual TGeoMaterial * DecayMaterial(Double_t time, Double_t precision=0.001)
Create the mixture representing the decay product of this material at a given time.
virtual void SavePrimitive(std::ostream &out, Option_t *option="")
Save a primitive as a C++ statement(s) on output stream "out".
virtual Int_t GetNelements() const
Definition: TGeoMaterial.h:172
virtual void GetElementProp(Double_t &a, Double_t &z, Double_t &w, Int_t i=0)
Single interface to get element properties.
Definition: TGeoMaterial.h:171
Int_t * fNatoms
Definition: TGeoMaterial.h:142
virtual void SetA(Double_t a)
Definition: TGeoMaterial.h:183
virtual Bool_t IsEq(const TGeoMaterial *other) const
Return true if the other material has the same physical properties.
virtual Bool_t IsMixture() const
Definition: TGeoMaterial.h:180
Double_t * fWeights
Definition: TGeoMaterial.h:141
void DefineElement(Int_t iel, Double_t a, Double_t z, Double_t weight)
Definition: TGeoMaterial.h:193
Double_t * GetAmixt() const
Definition: TGeoMaterial.h:174
The TNamed class is the base class for all named ROOT classes.
Definition: TNamed.h:29
An array of TObjects.
Definition: TObjArray.h:37
Mother of all ROOT objects.
Definition: TObject.h:37
R__ALWAYS_INLINE Bool_t TestBit(UInt_t f) const
Definition: TObject.h:172
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
Definition: TObject.cxx:694
auto * a
Definition: textangle.C:12