Gas mixture = mixture of gas components. More...

#include <TRT_PAI_gasMixture.h>

Inheritance diagram for TRT_PAI_gasMixture:

Collaboration diagram for TRT_PAI_gasMixture:

Public Member Functions
	TRT_PAI_gasMixture (const std::string &nm)
	Construct gas mixture. More...

void	addComponent (TRT_PAI_gasComponent *pgc, double frac)
	Add gas component to gas mixture. More...

const std::string &	getName ()
	Get name of gas mixture. More...

int	getNComponents ()
	Get number of different gas components in this gas mixture. More...

TRT_PAI_gasComponent *	getComponent (unsigned int n)
	Get gas component no. More...

double	getCompFraction (unsigned int n)
	Get fraction of gas component no. More...

int	getNElements ()
	Get number of different element in this gas mixture. More...

TRT_PAI_element *	getElement (unsigned int n)
	Get element no. More...

double	getElemWeight (unsigned int n)
	Get weight of element no. More...

void	showStructure ()
	Print out of structure of this gas mixture. More...

void	freezeGas ()
	Components can be added to gas mixture before freezeGas is called. More...

bool	msgLvl (const MSG::Level lvl) const
	Test the output level. More...

MsgStream &	msg () const
	The standard message stream. More...

MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...

void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
std::vector< TRT_PAI_gasComponent * >	m_pcomp

std::vector< double >	m_compFracs

std::vector< TRT_PAI_element * >	m_pelem

std::vector< double >	m_elemWeights

std::string	m_name

bool	m_gasFrozen

std::string	m_nm
	Message source name. More...

boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...

std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...

std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...

std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Detailed Description

Gas mixture = mixture of gas components.

Definition at line 18 of file TRT_PAI_gasMixture.h.

Constructor & Destructor Documentation

◆ TRT_PAI_gasMixture()

TRT_PAI_gasMixture::TRT_PAI_gasMixture ( const std::string & nm )

Construct gas mixture.

Parameters

nm	mixture name

Definition at line 13 of file TRT_PAI_gasMixture.cxx.

                                                           :
   AthMessaging(nm),
   m_name(nm),
   m_gasFrozen(0)
 {};

Member Function Documentation

◆ addComponent()

void TRT_PAI_gasMixture::addComponent	(	TRT_PAI_gasComponent *	pgc,
		double	frac
	)

Add gas component to gas mixture.

Parameters

pgc	pointer to gas component to be added
frac	fraction of this component in gas

Definition at line 21 of file TRT_PAI_gasMixture.cxx.

                                                                             {
   if ( m_gasFrozen ) {
     ATH_MSG_ERROR("gasMixture::addComponent: cannot add new gasComponent "
                   "- gas already frozen");
     return;
   }
  
   m_compFracs.push_back(frac);
   m_pcomp.push_back(pgc);
   return;
 }

◆ freezeGas()

void TRT_PAI_gasMixture::freezeGas ( )

Components can be added to gas mixture before freezeGas is called.

After call, gas is frozen, and no more components can be added.

Definition at line 35 of file TRT_PAI_gasMixture.cxx.

                                    {
   if ( m_gasFrozen ) return;
  
   int    nComp = m_pcomp.size();
   double wtot  = 0.;
   for (int j=0; j<nComp; j++) {
     if ( m_compFracs[j] <= 0. ) {
       ATH_MSG_ERROR("gasMixture::freezeGas: "
                     "A gasComponent has non-positive fraction");
       return;
     }
     wtot += m_compFracs[j];
   }
  
   if ( std::abs(wtot-1.) > 1e-5 ) {
     ATH_MSG_WARNING("gasMixture::freezeGas: "
                     "Gas fractions do not add to unity but " << wtot
                     << ". Re-normalizing!!");
   }
  
   for (int j=0; j<nComp; ++j ) {
     m_compFracs[j] /= wtot;
   }
  
   m_gasFrozen = 1;
  
   // Now, find the composition in terms of elements
  
   double w;
   TRT_PAI_element* pe;
   wtot = 0.;
   for (int j=0; j<nComp; j++) {
     for(int i=0; i<m_pcomp[j]->getNElementTypes(); i++) {
       pe = m_pcomp[j]->getElement(i);
       w  = m_compFracs[j] * m_pcomp[j]->getElementMultiplicity(i);
       wtot += w;
       for ( unsigned int k=0; k<m_pelem.size(); k++ ) {
     if ( pe == m_pelem[k] ) {
       m_elemWeights[k] += w;
       w = -999.;
       break;
     }
       }
       if ( w>0. ) {
     m_elemWeights.push_back(w);
     m_pelem.push_back(pe);
       }
     }
   }
   return;
 }

◆ getCompFraction()

double TRT_PAI_gasMixture::getCompFraction ( unsigned int n )

Get fraction of gas component no.

n of this gas mixture

Parameters

n	gas component number

Definition at line 124 of file TRT_PAI_gasMixture.cxx.

                                                          {
  
   if ( n >= m_compFracs.size() ) {
     ATH_MSG_ERROR("gasMixture::getCompFraction: out of bounds");
     return m_compFracs[0];
   }
  
   return m_compFracs[n];
 }

◆ getComponent()

TRT_PAI_gasComponent * TRT_PAI_gasMixture::getComponent ( unsigned int n )

Get gas component no.

n of this gas mixture

Parameters

n	gas component number

Definition at line 112 of file TRT_PAI_gasMixture.cxx.

                                                                      {
  
   if ( n >= m_pcomp.size() ) {
     ATH_MSG_ERROR("gasMixture::getComponent: out of bounds");
     return m_pcomp[0];
   };
  
   return m_pcomp[n];
 }

◆ getElement()

TRT_PAI_element * TRT_PAI_gasMixture::getElement ( unsigned int n )

Get element no.

n in this gas mixture

Parameters

n	element number

Definition at line 136 of file TRT_PAI_gasMixture.cxx.

                                                               {
  
   if ( n >= m_pelem.size() ) {
     ATH_MSG_ERROR("TRT_PAI_gasMixture::getElement: out of bounds");
     return m_pelem[0];
   };
  
   return m_pelem[n];
 }

◆ getElemWeight()

double TRT_PAI_gasMixture::getElemWeight ( unsigned int n )

Get weight of element no.

n in this gas mixture

Parameters

n	element number

Definition at line 148 of file TRT_PAI_gasMixture.cxx.

                                                        {
  
   if ( n >= m_elemWeights.size() ) {
     ATH_MSG_ERROR("gasMixture::getElemFraction:Error:out of bounds");
     return 0;
   };
  
   return m_elemWeights[n];
 }

◆ getName()

const std::string& TRT_PAI_gasMixture::getName ( )

inline

Get name of gas mixture.

Definition at line 38 of file TRT_PAI_gasMixture.h.

38 { return m_name; };

◆ getNComponents()

int TRT_PAI_gasMixture::getNComponents ( )

inline

Get number of different gas components in this gas mixture.

Definition at line 43 of file TRT_PAI_gasMixture.h.

43 {return m_compFracs.size();}

◆ getNElements()

int TRT_PAI_gasMixture::getNElements ( )

inline

Get number of different element in this gas mixture.

Definition at line 60 of file TRT_PAI_gasMixture.h.

60 {return m_elemWeights.size();}

◆ initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

 {
   m_imsg = Athena::getMessageSvc();
   m_lvl = m_imsg ?
     static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
     MSG::INFO;
 }

◆ msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 164 of file AthMessaging.h.

 {
   MsgStream* ms = m_msg_tls.get();
   if (!ms) {
     if (!m_initialized.test_and_set()) initMessaging();
     ms = new MsgStream(m_imsg,m_nm);
     m_msg_tls.reset( ms );
   }
  
   ms->setLevel (m_lvl);
   return *ms;
 }

◆ msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 179 of file AthMessaging.h.

180 { return msg() << lvl; }

◆ msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns: boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

 {
   if (!m_initialized.test_and_set()) initMessaging();
   if (m_lvl <= lvl) {
     msg() << lvl;
     return true;
   } else {
     return false;
   }
 }

◆ setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

 {
   m_lvl = lvl;
 }

◆ showStructure()

void TRT_PAI_gasMixture::showStructure ( )

Print out of structure of this gas mixture.

Definition at line 89 of file TRT_PAI_gasMixture.cxx.

                                        {
  
   if ( !m_gasFrozen ) {
     ATH_MSG_WARNING("gasMixture::showStructure: Showing structure of non-frozen gas");
   }
  
   ATH_MSG_INFO("The gas named '" << m_name << "' has the following components:");
  
   for (unsigned int i=0; i<m_compFracs.size(); i++) {
     msg(MSG::INFO) << " - " << m_compFracs[i]*100. << " percent "
                    << m_pcomp[i]->getName() << " consisting of: ";
     for (int j=0; j<m_pcomp[i]->getNElementTypes(); j++) {
       if ( j>0 ) msg(MSG::INFO) << ",";
       msg(MSG::INFO) << " " << m_pcomp[i]->getElementMultiplicity(j)
                      << " atoms " << m_pcomp[i]->getElement(j)->getName();
     }
     msg(MSG::INFO) << endmsg;
   }
   return;
 }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

◆ m_compFracs

std::vector<double> TRT_PAI_gasMixture::m_compFracs

private

Definition at line 87 of file TRT_PAI_gasMixture.h.

◆ m_elemWeights

std::vector<double> TRT_PAI_gasMixture::m_elemWeights

private

Definition at line 89 of file TRT_PAI_gasMixture.h.

◆ m_gasFrozen

bool TRT_PAI_gasMixture::m_gasFrozen

private

Definition at line 91 of file TRT_PAI_gasMixture.h.

◆ m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

◆ m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

◆ m_name

std::string TRT_PAI_gasMixture::m_name

private

Definition at line 90 of file TRT_PAI_gasMixture.h.

◆ m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_pcomp

std::vector<TRT_PAI_gasComponent*> TRT_PAI_gasMixture::m_pcomp

private

Definition at line 86 of file TRT_PAI_gasMixture.h.

◆ m_pelem

std::vector<TRT_PAI_element*> TRT_PAI_gasMixture::m_pelem

private

Definition at line 88 of file TRT_PAI_gasMixture.h.

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TRT_PAI_gasMixture()

Member Function Documentation

◆ addComponent()

◆ freezeGas()

◆ getCompFraction()

◆ getComponent()

◆ getElement()

◆ getElemWeight()

◆ getName()

◆ getNComponents()

◆ getNElements()

◆ initMessaging()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ setLevel()

◆ showStructure()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_compFracs

◆ m_elemWeights

◆ m_gasFrozen

◆ m_imsg

◆ m_lvl

◆ m_msg_tls

◆ m_name

◆ m_nm

◆ m_pcomp

◆ m_pelem