EmcStructure Class Reference

#include <EmcStructure.h>

Public Member Functions
	EmcStructure ()
	~EmcStructure ()
long	getIndex (unsigned int thetaIndex, unsigned int phiIndex) const
bool	isOutofAccep (unsigned int thetaIndex, unsigned int phiIndex) const
unsigned int	getNumberOfTheRings ()
unsigned int	getNumberOfXtals ()
unsigned int	startingTheta ()
unsigned int	crystalsInRing (unsigned int theta) const
void	setEmcStruc ()
long	getThisThetaMaxIndex (int Theta) const
unsigned int	getPartId (long Index) const
unsigned int	getTheta (long Index) const
unsigned int	getPhi (long Index) const
int	getGeomIndex (unsigned int PartId, unsigned int ThetaIndex, unsigned int PhiIndex) const

Detailed Description

Definition at line 5 of file EmcStructure.h.

Constructor & Destructor Documentation

EmcStructure()

EmcStructure::EmcStructure

(

)

Definition at line 11 of file EmcStructure.cxx.

{
 
}

~EmcStructure()

EmcStructure::~EmcStructure

(

)

Definition at line 19 of file EmcStructure.cxx.

{
}

Member Function Documentation

crystalsInRing()

unsigned int EmcStructure::crystalsInRing

(

unsigned int

theta

)

const

Definition at line 151 of file EmcStructure.cxx.

{
  //The theta index is defined by Endcap_east(0-5),Barrel(6-49),Endcap_west(50-55)
  //in Emc Bhabha Calibration
  unsigned int theXtalInRing;
 
  if (theta==0 || theta==1 || theta==55 || theta==54) 
    {
      theXtalInRing= aEmcId->getPHI_ENDCAP_MAX(0)+1;
    }
 
  if (theta==2 || theta==3 || theta==53 || theta==52) 
    {
      theXtalInRing=aEmcId->getPHI_ENDCAP_MAX(2)+1;
    }
 
  if (theta==4 || theta==5 || theta==51 || theta==50) 
    {
      theXtalInRing=aEmcId->getPHI_ENDCAP_MAX(4)+1;
    }
 
  if (theta>=6 && theta<=49) 
    {
      theXtalInRing=aEmcId->getPHI_BARREL_MAX()+1;
    }
 
  //  if (theta<0 || theta>55) theXtalInRing= 0;
  if (theta>55) 
    {
      theXtalInRing= 0;
    }
 
  return theXtalInRing;
}

Referenced by EmcSelBhaEvent::initGeom(), and setEmcStruc().

getGeomIndex()

int EmcStructure::getGeomIndex	(	unsigned int	PartId,
		unsigned int	ThetaIndex,
		unsigned int	PhiIndex ) const

Definition at line 53 of file EmcStructure.cxx.

{
  if (PartId>2){
    cout<<"PartId is out of EMC"<<endl;
    return -1;
  } else {
    unsigned int newThetaIndex;
    //The newThetaIndex is defined by Endcap_east(0-5),Barrel(6-49),Endcap_west(50-55)
    if (PartId==0) newThetaIndex = ThetaIndex;
    if (PartId==1) newThetaIndex = ThetaIndex + 6;
    if (PartId==2) newThetaIndex = 55 - ThetaIndex;
    
    if ( isOutofAccep(newThetaIndex, PhiIndex) ) 
      {
    return -1;
      } else {
      
    int index;   
    index = getIndex(newThetaIndex,PhiIndex);
      
    return index;
      }
  }
}

Referenced by BesEmcParameter::GetLightOutput().

getIndex()

long EmcStructure::getIndex	(	unsigned int	thetaIndex,
		unsigned int	phiIndex ) const

Definition at line 81 of file EmcStructure.cxx.

{
  long index = -1;
  if(thetaIndex==0){
    index=phiIndex;
  } else {
    index = getThisThetaMaxIndex( thetaIndex-1 ) +1 + phiIndex;
  }
  return index;
}

Referenced by getGeomIndex(), EmcSelBhaEvent::initGeom(), and setEmcStruc().

getNumberOfTheRings()

unsigned int EmcStructure::getNumberOfTheRings

(

)

Definition at line 124 of file EmcStructure.cxx.

{  
  unsigned int numberOfTheRings=0;
  numberOfTheRings  = (aEmcId->getTHETA_BARREL_MAX() + 1) 
    + (aEmcId->getTHETA_ENDCAP_MAX() + 1 ) * 2;
  
  return  numberOfTheRings;
  
}

Referenced by EmcSelBhaEvent::initGeom(), and setEmcStruc().

getNumberOfXtals()

unsigned int EmcStructure::getNumberOfXtals

(

)

Definition at line 134 of file EmcStructure.cxx.

{
  unsigned int numberOfXtals=0;
  unsigned int numberOfBarrelXtals, numberOfEndcapXtals;
 
  numberOfEndcapXtals = numberOfOneEndcapRings * (aEmcId->getPHI_ENDCAP_MAX(0)+1)
    + numberOfTwoEndcapRings * (aEmcId->getPHI_ENDCAP_MAX(2)+1)
    + numberOfThreeEndcapRings * (aEmcId->getPHI_ENDCAP_MAX(4)+1);
 
  numberOfBarrelXtals = (aEmcId->getTHETA_BARREL_MAX()+1)*(aEmcId->getPHI_BARREL_MAX()+1);
  
  numberOfXtals = numberOfBarrelXtals + numberOfEndcapXtals * 2;
  
  return numberOfXtals;
}

Referenced by EmcSelBhaEvent::initGeom().

getPartId()

unsigned int EmcStructure::getPartId ( long Index ) const

inline

Definition at line 35 of file EmcStructure.h.

{return m_partID[Index];}

Referenced by main().

getPhi()

unsigned int EmcStructure::getPhi ( long Index ) const

inline

Definition at line 37 of file EmcStructure.h.

{return m_phiIndex[Index];}

Referenced by main().

getTheta()

unsigned int EmcStructure::getTheta ( long Index ) const

inline

Definition at line 36 of file EmcStructure.h.

{return m_thetaIndex[Index];}

Referenced by main().

getThisThetaMaxIndex()

long EmcStructure::getThisThetaMaxIndex ( int Theta ) const

inline

Definition at line 34 of file EmcStructure.h.

{return m_ThetaMaxIndex[Theta];}

Referenced by getIndex().

isOutofAccep()

bool EmcStructure::isOutofAccep	(	unsigned int	thetaIndex,
		unsigned int	phiIndex ) const

Definition at line 93 of file EmcStructure.cxx.

{
  if ( thetaIndex > 55 ) 
    {
      cout<<"warning " << " theta out of acceptance !"
      << endl;
      return true;  
      
    } 
  else 
    {
      if ( ((thetaIndex==0 || thetaIndex==1 || thetaIndex==55 || thetaIndex==54)
        && (phiIndex > 63))
       || ((thetaIndex==2 || thetaIndex==3 || thetaIndex==53 || thetaIndex==52)
           && (phiIndex > 79))
       || ((thetaIndex==4 || thetaIndex==5 || thetaIndex==51 || thetaIndex==50)
           && ( phiIndex > 95))
       || ((thetaIndex>=6 && thetaIndex<=49) && ( phiIndex > 119)) ) 
    {      
      cout<<"warning "<< " phi out of acceptance !" << endl;
      return true;    
    }  
      else 
    {
      return false;
    } 
    }
}

Referenced by EmcBhabhaEvent::getDepoMCShowerEnergy(), EmcBhabhaEvent::getDepoMCShowerEnergy_lab(), EmcBhabhaEvent::getErrorDepoMCShowerEnergy(), and getGeomIndex().

setEmcStruc()

void EmcStructure::setEmcStruc

(

)

Definition at line 23 of file EmcStructure.cxx.

{
  int index;
  int nrOfTheRings = getNumberOfTheRings();
  for ( int the = startingTheta(); the< nrOfTheRings; the++) {
       
    for ( int phi=0; phi < crystalsInRing((unsigned int) the ); phi++) {
      
      index = getIndex( (unsigned int)the , (unsigned int)phi);
 
      m_phiIndex[index]=phi;
      if (the>=0 && the<=5) {
    m_partID[index]=0;
    m_thetaIndex[index]=the;
      }
      if (the>=6 && the<=49) {
    m_partID[index]=1;
    m_thetaIndex[index]=the-6;
      }
      if (the>=50 && the<=55){ 
    m_partID[index]=2;
    m_thetaIndex[index]=55 - the;
      }
 
    }
    m_ThetaMaxIndex[the]=index;
  }
 
}

Referenced by BesEmcParameter::GetLightOutput(), EmcSelBhaEvent::initGeom(), and main().

startingTheta()

unsigned int EmcStructure::startingTheta ( )

inline

Definition at line 26 of file EmcStructure.h.

{return 0;}

Referenced by EmcSelBhaEvent::initGeom(), and setEmcStruc().

---

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

• EmcStructure.h
• EmcStructure.cxx