/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Emc/EmcGeneralClass/EmcGeneralClass-00-00-04/src/EmcStructure.cxx

Go to the documentation of this file.

#include "EmcGeneralClass/EmcStructure.h"
#include <vector>
#include <iostream>

using namespace std;

//----------------
// Constructors --
//----------------
EmcStructure::EmcStructure()
{

}

//--------------
// Destructor --
//--------------
EmcStructure::~EmcStructure()
{
}

void EmcStructure::setEmcStruc()
{
  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;
  }

}
int 
EmcStructure::getGeomIndex( unsigned int PartId, 
                            unsigned int ThetaIndex, 
                            unsigned int PhiIndex)  const
{
  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
{
  long index = -1;
  if(thetaIndex==0){
    index=phiIndex;
  } else {
    index = getThisThetaMaxIndex( thetaIndex-1 ) +1 + phiIndex;
  }
  return index;
}

bool 
EmcStructure::isOutofAccep(unsigned int thetaIndex, 
                           unsigned int phiIndex) const
{
  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;
        } 
    }
}

unsigned int 
EmcStructure::getNumberOfTheRings()
{  
  unsigned int numberOfTheRings=0;
  numberOfTheRings  = (aEmcId->getTHETA_BARREL_MAX() + 1) 
    + (aEmcId->getTHETA_ENDCAP_MAX() + 1 ) * 2;
  
  return  numberOfTheRings;
  
}
unsigned int 
EmcStructure::getNumberOfXtals()
{
  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::crystalsInRing( unsigned int theta ) const
{
  //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;
}

---