MucGas Class Reference

#include <MucGas.h>

Inheritance diagram for MucGas:

Public Member Functions
	MucGas (int part, int segment, int layer, int upDown, int rpcId, int id)
	MucGas (const MucGas &other)
MucGas &	operator= (const MucGas &other)
	~MucGas ()
double	GetPhi ()
Public Member Functions inherited from MucEntity
	MucEntity (int part, int segment, int layer)
	MucEntity (int part, int segment, int layer, int id)
	MucEntity (int part, int segment, int layer, int upDown, int id)
	MucEntity (int part, int segment, int layer, int upDown, int rpcId, int id)
	MucEntity (const MucEntity &other)
MucEntity &	operator= (const MucEntity &other)
virtual	~MucEntity ()
int	GetPart ()
int	GetSegment ()
int	GetLayer ()
int	GetUpDown ()
int	GetRpcId ()
int	GetID ()
double	GetTheta ()
double	GetRin ()
double	GetRout ()
double	GetRc ()
double	GetThin ()
double	GetW ()
double	GetH ()
double	GetL ()
double	GetWu ()
double	GetWd ()
double	GetArea ()
double	GetLocOrgInBes (int i)
double	GetObjRotToMot (int i)
double	GetObjOrgInBes (int i)
double	GetObjOrgInLoc (int i)
double *	GetLocOrgInBes ()
double *	GetObjRotToMot ()
double *	GetObjOrgInLoc ()
double *	GetObjOrgInBes ()
void	TransBesToLoc (double LocCoord[], double BesCoord[], double LocOrgInBes[], double Rot_z)
void	TransLocToBes (double LocCoord[], double BesCoord[], double LocOrgInBes[], double Rot_z)

Protected Member Functions
virtual void	Init ()
virtual void	SetTheta ()
void	SetPhi ()
virtual void	SetRc ()
virtual void	SetRin ()
virtual void	SetRout ()
virtual void	SetThin ()
virtual void	SetW ()
virtual void	SetH ()
virtual void	SetL ()
virtual void	SetWu ()
virtual void	SetWd ()
int	SetEndcapRpcInBes ()
int	SetBarrelRpcInBes ()
virtual void	SetLocOrgInBes ()
virtual void	SetObjOrgInLoc ()
virtual void	SetObjOrgInBes ()
Protected Member Functions inherited from MucEntity
virtual void	SetArea ()
virtual void	SetObjRotToMot ()
virtual void	SetAlignment (double dx, double dy, double dz)

Additional Inherited Members
Protected Attributes inherited from MucEntity
int	m_Part
int	m_Segment
int	m_Layer
int	m_UpDown
int	m_RpcId
int	m_ID
double	m_Theta
double	m_Rin
double	m_Rout
double	m_Rc
double	m_Thin
double	m_W
double	m_H
double	m_L
double	m_Wu
double	m_Wd
double	m_Area
double	m_LocOrgInBes [3]
double	m_ObjRotToMot [3]
double	m_ObjOrgInBes [3]
double	m_ObjOrgInLoc [3]

Detailed Description

Definition at line 18 of file MucGas.h.

Constructor & Destructor Documentation

MucGas() [1/2]

MucGas::MucGas	(	int	part,
		int	segment,
		int	layer,
		int	upDown,
		int	rpcId,
		int	id )

Definition at line 17 of file MucGas.cxx.

  : MucEntity( part, segment, layer, upDown, rpcId, id )
{
  MucGas::Init();
}

MucGas() [2/2]

MucGas::MucGas

(

const MucGas &

other

)

Definition at line 24 of file MucGas.cxx.

                                    : MucEntity( other )
{
  m_Phi = other.m_Phi;
}

~MucGas()

MucGas::~MucGas

(

)

Definition at line 40 of file MucGas.cxx.

{ ; }

Member Function Documentation

GetPhi()

double MucGas::GetPhi

(

)

Definition at line 68 of file MucGas.cxx.

{ return m_Phi; }

Init()

void MucGas::Init ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 43 of file MucGas.cxx.

{
  SetPhi();
  SetTheta();
  SetRin();
  SetRout();
  SetRc();
  
  SetThin();
  SetW();
  SetWu();
  SetWd();
  SetH();
  SetL();
  
  SetEndcapRpcInBes();
  SetBarrelRpcInBes();
  
  SetLocOrgInBes();
  SetObjOrgInBes();
  SetObjOrgInLoc();
}

Referenced by MucGas().

operator=()

MucGas & MucGas::operator=

(

const MucGas &

other

)

Definition at line 30 of file MucGas.cxx.

{
  if( this == &other) return *this;
  MucEntity::operator =(other);
  m_Phi = other.m_Phi;
  
  return *this; 
}

SetBarrelRpcInBes()

int MucGas::SetBarrelRpcInBes ( )

protected

Definition at line 355 of file MucGas.cxx.

{
  if( m_Part != BRID ) return 0;
  
  double x, y, z;
  x = y = z = 0.;
  
  if( m_Segment != B_TOP )
  {
    // set x and y
    double r;
    r = (B_AS_RMIN[m_Layer] - AS_GAP/2.0) + (1-2*m_UpDown)*(STR_TH + RPC_TH)/2.0;
    x = r * cos( m_Theta );
    y = r * sin( m_Theta );
    
    // set z
    if( (m_UpDown == SL_UP && m_RpcId == 3) || (m_UpDown == SL_DOWN && m_RpcId == 0) )
    {
      if( m_Layer != 1 ) // not layer NO.2
        z = (2*m_UpDown - 1) * 3 * B_RPC_LTN/2.0;
      else
        z = (2*m_UpDown - 1) * 3 * B_RPC_LTS/2.0;
    }
    else if( m_UpDown == SL_UP )
      z = B_BOX_LT/2.0 - B_BOX_DZ[m_UpDown] - (1+2*m_RpcId)*B_RPC_LTN/2.0;
    else
      z = -B_BOX_LT/2.0 + B_BOX_DZ[m_UpDown] + (1+2*(3-m_RpcId))*B_RPC_LTN/2.0;
    
  } // end not top segment
  else // top segment
  {
    //---------- set x--------
    if( m_RpcId > 1 ) // NO.3, 4, 5
      x = 0.;
    else              // NO.1 and NO.2
      x = (1-2*m_RpcId) * ( (B_BOX_WT[m_Layer] + B_BOX_SLOT_WT)/4.0 + (B_BOX_DZ[m_UpDown] - B_BOX_DT[m_UpDown])/2.0 );
            
    //---------- set y--------
    y = (B_AS_RMIN[m_Layer] - AS_GAP / 2.0) + (1-2*m_UpDown)*(STR_TH + RPC_TH)/2.0;
    
    //---------- set z--------
    if( m_RpcId < 2 )             // rpc NO.1 and NO.2
      z = B_BOX_LT/2.0 - B_BOX_DZ[m_UpDown] - B_TOPRPC_LTS[ ((m_Layer==0)?1:(m_Layer%2)) ];
    else if( m_UpDown == SL_UP && m_RpcId == 4 )    // up NO.5
      z = -B_RPC_LTN - B_TOPRPC_LTS[ ((m_Layer==0)?1:(m_Layer%2)) ];
    else if( m_UpDown == SL_DOWN && m_RpcId == 2 )  // down NO.3
      z = B_RPC_LTN - B_TOPRPC_LTS[ ((m_Layer==0)?1:(m_Layer%2)) ];
    else if( m_UpDown == SL_UP )  // up NO.3 and NO.4
      z = B_BOX_LT/2.0 - B_BOX_DZ[m_UpDown] - B_TOPRPC_LTS[ ((m_Layer==0)?1:(m_Layer%2)) ] + (3-2*m_RpcId)*B_RPC_LTN/2.0;
    else                          // down NO.4 and NO.5
      z = -B_BOX_LT/2.0 - B_BOX_DZ[m_UpDown] + (1+2*(4-m_RpcId))*B_RPC_LTN/2.0;  
  } // end top segment
  
  m_LocOrgInBes[0] = x;
  m_LocOrgInBes[1] = y;
  m_LocOrgInBes[2] = z;
  
  // limit cut
  for( int i=0; i<2; i++ ) 
    if( fabs(m_LocOrgInBes[i]) < ERR_LIMIT ) m_LocOrgInBes[i] = 0.0;
        
  return 0;
}

Referenced by Init(), and SetLocOrgInBes().

SetEndcapRpcInBes()

int MucGas::SetEndcapRpcInBes ( )

protected

Definition at line 419 of file MucGas.cxx.

{
  if( m_Part == BRID ) return 0;
  
  double x, y, z;
  x = y = z = 0.;
  
  //------------ set x and y ---------------
  if( m_UpDown == SL_UP )
  {
    if( m_RpcId == 0 ) {
      x = ( VALUE*m_Rc + E_GP_DX + E_STRPLN_DA + RPC_SL_OFFSET ) / 2.0;
      y = m_Rc;
    }   
    else if( m_RpcId == 1 ) {
      x = m_Rc * cos( MUC_PI/4.0 );
      y = m_Rc * sin( MUC_PI/4.0 );
    }
    else {
      x = m_Rc;
      y = ( VALUE*m_Rc + E_GP_DY + E_STRPLN_DA + RPC_SL_OFFSET ) / 2.0;
    }       
  }
  else
  {
    switch( m_RpcId )
    {
      case 0:
              x = ( VALUE*m_Rc + E_GP_DX + E_STRPLN_DA )/2.0 + RPC_DEATH_EDGE/2.0;
              y = m_Rc;
              break;
      case 1:
              x = VALUE*m_Rc + sin(MUC_PI/4.0)*m_Wd/4.0;
              x += RPC_DEATH_EDGE * ( 3 - 2*sqrt(2.0) )/4.0;
              y = m_Rc - cos(MUC_PI/4.0)*m_Wd/4.0;
              y += RPC_DEATH_EDGE * ( 1 + 4.0/sqrt(2.0) )/4.0;
              break;
      case 2:
              x = m_Rc * cos( MUC_PI/4.0 ); 
              y = m_Rc * sin( MUC_PI/4.0 );
              break;
      case 3:
              x = m_Rc - cos(MUC_PI/4.0)*m_Wd/4.0;
              x += RPC_DEATH_EDGE * ( 1 + 4.0/sqrt(2.0) )/4.0;
              y = VALUE*m_Rc + sin(MUC_PI/4.0)*m_Wd/4.0;
              y += RPC_DEATH_EDGE * ( 3 - 2*sqrt(2.0) )/4.0;
              break;
      case 4:
              x = m_Rc;
              y = ( VALUE*m_Rc + E_GP_DY + E_STRPLN_DA )/2.0 + RPC_DEATH_EDGE/2.0;
    }
  }
 
  // sign different by segment
  if      ( m_Segment == 0 )  { ;               }
  else if ( m_Segment == 1 )  { x = -x;         }
  else if ( m_Segment == 2 )  { x = -x; y = -y; }                           
  else                        { y = -y;         }
  
  //------------ set z ---------------------
  for( int i=0; i<m_Layer+1; i++ ) z += E_AS_TH[i]; 
  z += m_Layer * AS_GAP;
  z += ( (E_AS_ZMAX - E_AS_TOTAL_TH) + AS_GAP/2.0 );  
  z += (1-2*m_UpDown) * ( STR_TH + RPC_TH ) / 2.0;  
  z *= cos( m_Part*MUC_PI/2.0 );
  
  m_LocOrgInBes[0] = x;
  m_LocOrgInBes[1] = y;
  m_LocOrgInBes[2] = z;
  
  // limit cut
  for( int i=0; i<2; i++ ) 
    if( fabs(m_LocOrgInBes[i]) < ERR_LIMIT ) m_LocOrgInBes[i] = 0.;
    
  return 0;
}

Referenced by Init(), and SetLocOrgInBes().

SetH()

void MucGas::SetH ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 194 of file MucGas.cxx.

{
  if( m_Part == BRID )
    m_H = GAS_TH;
  else
  {
    m_H = m_Rout - m_Rin;
    
    // cut rpc chamber death edge
    if( m_UpDown == SL_DOWN && (m_ID == 1 || m_ID == 3) )
      m_H -= RPC_DEATH_EDGE * (1 + 1.0/VALUE);
    else
      m_H -= RPC_DEATH_EDGE * 2.0;
  }
}

Referenced by Init().

SetL()

void MucGas::SetL ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 210 of file MucGas.cxx.

{
  if( m_Part == BRID ) // barrel
  {
    if( m_Segment != B_TOP ) // not top segment
    {
      // the first rpc of SL_DOWN is constructed as the same as the last rpc of SL_UP
      if( (m_UpDown == SL_UP && m_RpcId == 3) || (m_UpDown == SL_DOWN && m_RpcId == 0) )
      {
        if( m_Layer != 1 )  // not layer NO.2
          m_L = B_BOX_LT - 2*B_BOX_DZ[m_UpDown] - 3*B_RPC_LTN;
        else                // layer NO.2 is special
          m_L = B_BOX_LT - 2*B_BOX_DZ[m_UpDown] - 3*B_RPC_LTS;
      }
      else
      {
        if( m_Layer != 1 )  // not layer NO.2
          m_L = B_RPC_LTN;
        else                // layer NO.2 is special
          m_L = B_RPC_LTS;
      }
    }
    else // top segment
    {
      if( m_RpcId < 2 ) // RPC NO.1 and NO.2
        m_L = B_TOPRPC_LTS[ ((m_Layer == 0)?1:(m_Layer%2)) ];
      else if( (m_UpDown == SL_UP && m_RpcId == 4) || (m_UpDown == SL_DOWN && m_RpcId == 2) )
        m_L = B_BOX_LT - B_TOPRPC_LTS[ ((m_Layer == 0)?1:(m_Layer%2)) ] - 2*(B_BOX_DZ[m_UpDown]+B_RPC_LTN);
      else
        m_L = B_RPC_LTN;
    } // else, top segment  
    
    m_L -= RPC_DEATH_EDGE * 2.0;    // cut rpc chamber death edge
  
  } // End barrel 
  else // endcap
    m_L = GAS_TH;
}

Referenced by Init().

SetLocOrgInBes()

void MucGas::SetLocOrgInBes ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 497 of file MucGas.cxx.

{
  if( m_Part == BRID )  SetBarrelRpcInBes();
  else                  SetEndcapRpcInBes();
}

Referenced by Init().

SetObjOrgInBes()

void MucGas::SetObjOrgInBes ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 509 of file MucGas.cxx.

{
  MucEntity::TransLocToBes( m_ObjOrgInBes, m_ObjOrgInLoc, m_LocOrgInBes, m_Phi );
  // limit cut
  for( int i=0; i<3; i++ ) 
    if( fabs(m_ObjOrgInBes[i]) < ERR_LIMIT ) m_ObjOrgInBes[i] = 0.0;
}

Referenced by Init().

SetObjOrgInLoc()

void MucGas::SetObjOrgInLoc ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 504 of file MucGas.cxx.

{
  for( int i=0; i<3; i++ )  m_ObjOrgInLoc[i] = 0.;
}

Referenced by Init().

SetPhi()

void MucGas::SetPhi ( )

protected

Definition at line 73 of file MucGas.cxx.

{
  if( m_Part == BRID )  m_Phi = B_PHI[m_Segment] * (MUC_PI/4.0);
  else                  m_Phi = 0.;
}

Referenced by Init().

SetRc()

void MucGas::SetRc ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 113 of file MucGas.cxx.

{
  m_Rc = (m_Rin + m_Rout)/2.0;
}

Referenced by Init().

SetRin()

void MucGas::SetRin ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 97 of file MucGas.cxx.

{
  if( m_Part == BRID )  // Rin of gas( same as which of RPC )
    m_Rin = B_AS_RMIN[m_Layer] - AS_GAP/2.0 + (1-2*m_UpDown) * (STR_TH + RPC_TH)/2.0 - GAS_TH/2.0;
  else                  // Rin of RPC
    m_Rin = E_GP_RMIN[m_Layer] + E_STRPLN_DR + (1-m_UpDown)*RPC_SL_OFFSET;
}

Referenced by Init().

SetRout()

void MucGas::SetRout ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 105 of file MucGas.cxx.

{
  if( m_Part == BRID )
    m_Rout = B_AS_RMIN[m_Layer] - AS_GAP/2.0 + (1-2*m_UpDown) * (STR_TH + RPC_TH)/2.0 + GAS_TH/2.0;
  else
    m_Rout = E_AS_RMAX - E_BOX_DR - E_STRPLN_DB[(m_Layer==0)?0:1] + (m_UpDown-1)*RPC_SL_OFFSET;
}

Referenced by Init().

SetTheta()

void MucGas::SetTheta ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 79 of file MucGas.cxx.

{
  if( m_Part == BRID )      m_Theta = m_Segment * (MUC_PI/4.0);
  else
  {
    if( m_UpDown == SL_UP ) m_Theta = ( MUC_PI/4.0 ) + (1 - m_RpcId) * MUC_PI / 8.0 ;
    else
    {
      if( m_RpcId < 2 ) // RPC NO.1 and NO.2 
        m_Theta = ( MUC_PI/4.0 ) + MUC_PI / 8.0 ;
      else if ( m_RpcId == 2 ) // center RPC, NO.3
        m_Theta = ( MUC_PI/4.0 ) ;
      else
        m_Theta = ( MUC_PI/4.0 ) - MUC_PI / 8.0 ;
    }
  }
}

Referenced by Init().

SetThin()

void MucGas::SetThin ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 119 of file MucGas.cxx.

{
  m_Thin = GAS_TH;
}

Referenced by Init().

SetW()

void MucGas::SetW ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 124 of file MucGas.cxx.

{
  if( m_Part == BRID)
  {
    if( m_Segment == B_TOP && m_RpcId < 2 ) // the first 2 gases in top segment
    {
      if( (m_Layer+1)%2 == 1 )  // odd layer
        m_W = (B_BOX_WT[m_Layer] - 2*B_BOX_DT[m_UpDown] - 2*B_BOX_DZ[m_UpDown] - B_BOX_SLOT_WT)/2.0;
      else                      // even layer
        m_W = (B_BOX_WT[m_Layer] - 4*B_BOX_DZ[m_UpDown] - B_BOX_SLOT_WT)/2.0;
    }   
    else                                    // other gases
    {
      if( (m_Layer+1)%2 ==1 ) // odd layer
        m_W = B_BOX_WT[m_Layer] - 2 * B_BOX_DT[m_UpDown];
      else                    // even layer
        m_W = B_BOX_WT[m_Layer] - 2 * B_BOX_DZ[m_UpDown];
    }   
  } // End barrel
  else
  {
    if( m_UpDown == SL_UP ) // up layer
    {   
      if( m_RpcId == 0 )
      m_W = VALUE * m_Rc - E_GP_DX - E_STRPLN_DA - RPC_SL_OFFSET;
      else if ( m_RpcId == 1 )
      m_W = 2 * VALUE * m_Rc;
      else 
      m_W = VALUE * m_Rc - E_GP_DY - E_STRPLN_DA - RPC_SL_OFFSET;
    }
    else // down layer
    {   
      switch( m_RpcId )
      {
        case 0:
                m_W = VALUE * m_Rc - E_GP_DX - E_STRPLN_DA;
                break;
        case 1:
        case 3:
                m_W = VALUE * ( m_Rc - m_Rin );
                break;
        case 2:
                m_W = 2 * VALUE * m_Rin;
                break;
        case 4:
                m_W = VALUE * m_Rc - E_GP_DY - E_STRPLN_DA; 
                break;
        default: ;
      }
    }
  } // else, endcap
  
  // cut rpc chamber death edge
  if( m_Part == BRID )  m_W -= RPC_DEATH_EDGE * 2.0;
  else
  {
    if( m_UpDown == SL_UP ) {
      if( m_RpcId == 1 )  m_W -= RPC_DEATH_EDGE * 2.0; // center rpc
      else                m_W -= ( 1 + 1.0/cos(MUC_PI/8.0) ) * RPC_DEATH_EDGE;
    }
    else
    {
      if( m_RpcId == 2 )  // center rpc
        m_W -= RPC_DEATH_EDGE * 2.0;
      else                // DOWN 0,1,3,4 RPC
        m_W -= RPC_DEATH_EDGE;
    }
  }
}

Referenced by Init().

SetWd()

void MucGas::SetWd ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 303 of file MucGas.cxx.

{
  if( m_Part == BRID )  m_Wd = m_W;
  else
  {
    if( m_UpDown == SL_UP )
    {
      if( m_RpcId == 0 )
        m_Wd = VALUE * m_Rout - E_GP_DX - E_STRPLN_DA - RPC_SL_OFFSET;
      else if( m_RpcId == 1 )
        m_Wd = 2.0* VALUE * m_Rout;
      else
        m_Wd = VALUE * m_Rout - E_GP_DY - E_STRPLN_DA - RPC_SL_OFFSET;
    }
    else
    {
      switch( m_RpcId )
      {
        case 0:
                m_Wd = VALUE * m_Rout - E_GP_DX - E_STRPLN_DA;
                break;
        case 1:
        case 3:
                m_Wd = VALUE * (m_Rout - m_Rin);
                break;
        case 2:
                m_Wd = 2.0 * VALUE * m_Rin;
                break;
        case 4:
                m_Wd = VALUE * m_Rout - E_GP_DY - E_STRPLN_DA;
                break;
        default: ;
      }
    }
    
    // cut rpc chamber death edge
    if( m_UpDown == SL_UP )
    {
      if( m_RpcId == 1 )  m_Wd -= RPC_DEATH_EDGE * 2.0; // center rpc
      else                m_Wd -= ( 1 + 1.0/cos(MUC_PI/8.0) ) * RPC_DEATH_EDGE;
    }
    else
    {
      if( m_RpcId == 2 )  // center rpc
        m_Wd -= RPC_DEATH_EDGE * 2.0;
      else                // DOWN 0,4 RPC
        m_Wd -= RPC_DEATH_EDGE;
    }  
  } // else, endcap
}

Referenced by Init().

SetWu()

void MucGas::SetWu ( )

protectedvirtual

Reimplemented from MucEntity.

Definition at line 249 of file MucGas.cxx.

{
  if( m_Part == BRID )  m_Wu = m_W;
  else
  { 
    if( m_UpDown == SL_UP )
    {
      if( m_RpcId == 0 )
        m_Wu = VALUE * m_Rin - E_GP_DX - E_STRPLN_DA - RPC_SL_OFFSET;
      else if( m_RpcId == 1 )
        m_Wu = 2.0* VALUE * m_Rin;
      else 
        m_Wu = VALUE * m_Rin - E_GP_DY - E_STRPLN_DA - RPC_SL_OFFSET;
    }
    else
    {
      switch( m_RpcId )
      {
        case 0:
                m_Wu = VALUE * m_Rin - E_GP_DX - E_STRPLN_DA;
                break;
        case 1:
        case 3:
                m_Wu = 0.;      
                break;
        case 2:
                m_Wu = 2.0* VALUE * m_Rin;
                break;
        case 4:
                m_Wu = VALUE * m_Rin - E_GP_DY - E_STRPLN_DA;
                break;
        default: ;
      }
    }
  
    // cut rpc chamber death edge
    if( m_Wu > 0.0 )
    {
      if( m_UpDown == SL_UP )
      {
        if( m_RpcId == 1 )  m_Wu -= RPC_DEATH_EDGE * 2.0; // center rpc
        else                m_Wu -= ( 1 + 1.0/cos(MUC_PI/8.0) ) * RPC_DEATH_EDGE;
      }
      else
      {
        if( m_RpcId == 2 ) // center rpc
          m_Wu -= RPC_DEATH_EDGE * 2.0;
        else // DOWN 0,4 RPC
          m_W -= RPC_DEATH_EDGE;
      }
    }
  } // else, endcap
}

Referenced by Init().

---

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

• MucGas.h
• MucGas.cxx