EmcStructure Class Reference

#include <EmcStructure.h>

List of all members.

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
Private Attributes
EmcID *	aEmcId
long	m_ThetaMaxIndex [56]
unsigned int	m_partID [6240]
unsigned int	m_thetaIndex [6240]
unsigned int	m_phiIndex [6240]
Static Private Attributes
static const unsigned int	numberOfOneEndcapRings = 2
static const unsigned int	numberOfTwoEndcapRings = 2
static const unsigned int	numberOfThreeEndcapRings = 2

---

Detailed Description

Definition at line 5 of file EmcStructure.h.

---

Constructor & Destructor Documentation

EmcStructure::EmcStructure

(

)

Definition at line 11 of file EmcStructure.cxx.

{

}

EmcStructure::~EmcStructure

(

)

Definition at line 19 of file EmcStructure.cxx.

{
}

---

Member Function Documentation

unsigned int EmcStructure::crystalsInRing

(

unsigned int

theta

)

const

Definition at line 151 of file EmcStructure.cxx.

References aEmcId, EmcID::getPHI_BARREL_MAX(), and EmcID::getPHI_ENDCAP_MAX().

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

{
  //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;
}

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

Definition at line 53 of file EmcStructure.cxx.

References getIndex(), and isOutofAccep().

Referenced by EmcCalibConstSvc::getIndex(), and BesEmcParameter::GetLightOutput().

{
  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;
      }
  }
}

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

Definition at line 81 of file EmcStructure.cxx.

References getThisThetaMaxIndex().

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

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

unsigned int EmcStructure::getNumberOfTheRings

(

)

Definition at line 124 of file EmcStructure.cxx.

References aEmcId, EmcID::getTHETA_BARREL_MAX(), and EmcID::getTHETA_ENDCAP_MAX().

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

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

unsigned int EmcStructure::getNumberOfXtals

(

)

Definition at line 134 of file EmcStructure.cxx.

References aEmcId, EmcID::getPHI_BARREL_MAX(), EmcID::getPHI_ENDCAP_MAX(), EmcID::getTHETA_BARREL_MAX(), numberOfOneEndcapRings, numberOfThreeEndcapRings, and numberOfTwoEndcapRings.

Referenced by EmcSelBhaEvent::initGeom().

{
  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;
}

unsigned int EmcStructure::getPartId

(

long

Index

)

const [inline]

Definition at line 35 of file EmcStructure.h.

References m_partID.

Referenced by EmcSelBhaEvent::expectedEnergy(), EmcCalibConstSvc::getPartID(), and main().

{return m_partID[Index];}

unsigned int EmcStructure::getPhi

(

long

Index

)

const [inline]

Definition at line 37 of file EmcStructure.h.

References m_phiIndex.

Referenced by EmcSelBhaEvent::expectedEnergy(), EmcCalibConstSvc::getPhiIndex(), and main().

{return m_phiIndex[Index];}

unsigned int EmcStructure::getTheta

(

long

Index

)

const [inline]

Definition at line 36 of file EmcStructure.h.

References m_thetaIndex.

Referenced by EmcSelBhaEvent::expectedEnergy(), EmcCalibConstSvc::getThetaIndex(), and main().

{return m_thetaIndex[Index];}

long EmcStructure::getThisThetaMaxIndex

(

int

Theta

)

const [inline]

Definition at line 34 of file EmcStructure.h.

References m_ThetaMaxIndex.

Referenced by getIndex().

{return m_ThetaMaxIndex[Theta];}

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

Definition at line 93 of file EmcStructure.cxx.

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

{
  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;
        } 
    }
}

void EmcStructure::setEmcStruc

(

)

Definition at line 23 of file EmcStructure.cxx.

References crystalsInRing(), getIndex(), getNumberOfTheRings(), m_partID, m_phiIndex, m_thetaIndex, m_ThetaMaxIndex, and startingTheta().

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

{
  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;
  }

}

unsigned int EmcStructure::startingTheta

(

)

[inline]

Definition at line 26 of file EmcStructure.h.

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

{return 0;}

---

Member Data Documentation

EmcID* EmcStructure::aEmcId [private]

Definition at line 44 of file EmcStructure.h.

Referenced by crystalsInRing(), getNumberOfTheRings(), and getNumberOfXtals().

unsigned int EmcStructure::m_partID[6240] [private]

Definition at line 49 of file EmcStructure.h.

Referenced by getPartId(), and setEmcStruc().

unsigned int EmcStructure::m_phiIndex[6240] [private]

Definition at line 51 of file EmcStructure.h.

Referenced by getPhi(), and setEmcStruc().

unsigned int EmcStructure::m_thetaIndex[6240] [private]

Definition at line 50 of file EmcStructure.h.

Referenced by getTheta(), and setEmcStruc().

long EmcStructure::m_ThetaMaxIndex[56] [private]

Definition at line 48 of file EmcStructure.h.

Referenced by getThisThetaMaxIndex(), and setEmcStruc().

const unsigned int EmcStructure::numberOfOneEndcapRings = 2 [static, private]

Definition at line 45 of file EmcStructure.h.

Referenced by getNumberOfXtals().

const unsigned int EmcStructure::numberOfThreeEndcapRings = 2 [static, private]

Definition at line 47 of file EmcStructure.h.

Referenced by getNumberOfXtals().

const unsigned int EmcStructure::numberOfTwoEndcapRings = 2 [static, private]

Definition at line 46 of file EmcStructure.h.

Referenced by getNumberOfXtals().

---