CaloGeometry Class Reference

#include <CaloGeometry.h>

Inheritance diagram for CaloGeometry:

Public Member Functions
	CaloGeometry ()
virtual	~CaloGeometry ()
virtual bool	PostProcessGeometry ()
virtual void Validate	ATLAS_NOT_THREAD_SAFE (int nrnd=100)
virtual const CaloDetDescrElement *	getDDE (Identifier identify) const
virtual const CaloDetDescrElement *	getDDE (int sampling, Identifier identify) const
virtual const CaloDetDescrElement *	getDDE (int sampling, float eta, float phi, float *distance=0, int *steps=0) const
virtual const CaloDetDescrElement *	getFCalDDE (int sampling, float x, float y, float z, float *distance=0, int *steps=0) const
bool	getClosestFCalCellIndex (int sampling, float x, float y, int &ieta, int &iphi, int *steps=0) const
double	deta (int sample, double eta) const
void	minmaxeta (int sample, double eta, double &mineta, double &maxeta) const
double	rzmid (int sample, double eta) const
double	rzent (int sample, double eta) const
double	rzext (int sample, double eta) const
double	rmid (int sample, double eta) const
double	rent (int sample, double eta) const
double	rext (int sample, double eta) const
double	zmid (int sample, double eta) const
double	zent (int sample, double eta) const
double	zext (int sample, double eta) const
double	rpos (int sample, double eta, int subpos=CaloSubPos::SUBPOS_MID) const
double	zpos (int sample, double eta, int subpos=CaloSubPos::SUBPOS_MID) const
double	rzpos (int sample, double eta, int subpos=CaloSubPos::SUBPOS_MID) const
bool	isCaloBarrel (int sample) const
TGraphErrors *	GetGraph (unsigned int sample) const
void	SetDoGraphs (bool dographs=true)
bool	DoGraphs () const
TGraph *	DrawGeoSampleForPhi0 (int sample, int calocol, bool print=false)
TCanvas *	DrawGeoForPhi0 ()
FCAL_ChannelMap *	GetFCAL_ChannelMap ()
void	SetFCal_ChannelMap (const FCAL_ChannelMap *fcal_ChannnelMap)
void	calculateFCalRminRmax ()
virtual bool	checkFCalGeometryConsistency ()
virtual void	PrintMapInfo (int i, int j)

Static Public Member Functions
static std::string	SamplingName (int sample)

Static Public Attributes
static const int	MAX_SAMPLING = CaloCell_ID_FCS::MaxSample
static const Identifier	m_debug_identify
static std::atomic< bool >	m_debug =false

Protected Member Functions
virtual void	addcell (const CaloDetDescrElement *cell)
virtual void	post_process (int layer)
virtual void	InitRZmaps ()

Protected Attributes
t_cellmap	m_cells
std::vector< t_cellmap >	m_cells_in_sampling
std::vector< t_eta_cellmap >	m_cells_in_sampling_for_phi0
std::vector< std::vector< CaloGeometryLookup * > >	m_cells_in_regions
std::vector< bool >	m_isCaloBarrel
std::vector< double >	m_min_eta_sample [2]
std::vector< double >	m_max_eta_sample [2]
std::vector< FSmap< double, double > >	m_rmid_map [2]
std::vector< FSmap< double, double > >	m_zmid_map [2]
std::vector< FSmap< double, double > >	m_rent_map [2]
std::vector< FSmap< double, double > >	m_zent_map [2]
std::vector< FSmap< double, double > >	m_rext_map [2]
std::vector< FSmap< double, double > >	m_zext_map [2]
bool	m_dographs
std::vector< TGraphErrors * >	m_graph_layers
FCAL_ChannelMap	m_FCal_ChannelMap
std::vector< double >	m_FCal_rmin
std::vector< double >	m_FCal_rmax

Detailed Description

Definition at line 27 of file CaloGeometry.h.

Constructor & Destructor Documentation

CaloGeometry()

CaloGeometry::CaloGeometry

(

)

Definition at line 41 of file CaloGeometry.cxx.

                           : m_cells_in_sampling(MAX_SAMPLING),m_cells_in_sampling_for_phi0(MAX_SAMPLING),m_cells_in_regions(MAX_SAMPLING),m_isCaloBarrel(MAX_SAMPLING),m_dographs(false),m_FCal_ChannelMap(0)
{
  //TMVA::Tools::Instance();
  for(int i=0;i<2;++i) {
    m_min_eta_sample[i].resize(MAX_SAMPLING); //[side][calosample]
    m_max_eta_sample[i].resize(MAX_SAMPLING); //[side][calosample]
    m_rmid_map[i].resize(MAX_SAMPLING); //[side][calosample]
    m_zmid_map[i].resize(MAX_SAMPLING); //[side][calosample]
    m_rent_map[i].resize(MAX_SAMPLING); //[side][calosample]
    m_zent_map[i].resize(MAX_SAMPLING); //[side][calosample]
    m_rext_map[i].resize(MAX_SAMPLING); //[side][calosample]
    m_zext_map[i].resize(MAX_SAMPLING); //[side][calosample]
  }
  m_graph_layers.resize(MAX_SAMPLING);
//  for(unsigned int i=CaloCell_ID_FCS::FirstSample;i<CaloCell_ID_FCS::MaxSample;++i) {
  for(unsigned int i=CaloSampling::PreSamplerB;i<=CaloSampling::FCAL2;++i) {
    m_graph_layers[i]=nullptr;
    CaloSampling::CaloSample s=static_cast<CaloSampling::CaloSample>(i);
    m_isCaloBarrel[i]=(CaloSampling::barrelPattern() & CaloSampling::getSamplingPattern(s))!=0;
  }
  m_isCaloBarrel[CaloCell_ID_FCS::TileGap3]=false; 
}

~CaloGeometry()

CaloGeometry::~CaloGeometry ( )

virtualdefault

Member Function Documentation

addcell()

void CaloGeometry::addcell ( const CaloDetDescrElement * cell )

protectedvirtual

Definition at line 67 of file CaloGeometry.cxx.

{
  int sampling=cell->getSampling();
  Identifier identify=cell->identify();
  
  m_cells[identify]=cell;
  m_cells_in_sampling[sampling][identify]=cell;
  
  //m_cells[identify]= new CaloDetDescrElement(*cell);
  //m_cells_in_sampling[sampling][identify]= new CaloDetDescrElement(*cell);
  
  CaloGeometryLookup* lookup=nullptr;
  for(unsigned int i=0;i<m_cells_in_regions[sampling].size();++i) {
    if(m_cells_in_regions[sampling][i]->IsCompatible(cell)) {
      lookup=m_cells_in_regions[sampling][i];
      //cout<<sampling<<":  found compatible map"<<endl;
      break;
    }
  }
  if(!lookup) {
    lookup=new CaloGeometryLookup(m_cells_in_regions[sampling].size());
    m_cells_in_regions[sampling].push_back(lookup);
  }
  lookup->add(cell);
}

ATLAS_NOT_THREAD_SAFE()

virtual void Validate CaloGeometry::ATLAS_NOT_THREAD_SAFE ( int nrnd = 100 )

virtual

Implements ICaloGeometry.

calculateFCalRminRmax()

void CaloGeometry::calculateFCalRminRmax

(

)

Definition at line 905 of file CaloGeometry.cxx.

                                         {
   
   m_FCal_rmin.resize(3,FLT_MAX);
   m_FCal_rmax.resize(3,0.);
 
   double x(0.),y(0.),r(0.);
   for(int imap=1;imap<=3;imap++)for(auto it=m_FCal_ChannelMap.begin(imap);it!=m_FCal_ChannelMap.end(imap);it++){
      x=it->second.x();
      y=it->second.y();
      r=sqrt(x*x+y*y);
      if(r<m_FCal_rmin[imap-1])m_FCal_rmin[imap-1]=r;
      if(r>m_FCal_rmax[imap-1])m_FCal_rmax[imap-1]=r;
   }
   
}

checkFCalGeometryConsistency()

bool CaloGeometry::checkFCalGeometryConsistency ( )

virtual

Reimplemented in CaloGeometryFromFile.

Definition at line 922 of file CaloGeometry.cxx.

                                               {
  
  unsigned long long phi_index,eta_index;
  float x,y,dx,dy;
  long id;
  
  long mask1[]{0x34,0x34,0x35};
  long mask2[]{0x36,0x36,0x37};
  
  m_FCal_rmin.resize(3,FLT_MAX);
  m_FCal_rmax.resize(3,0.);
  
  
  for(int imap=1;imap<=3;imap++){
    
    int sampling = imap+20;
    
    if((int)m_cells_in_sampling[sampling].size() != 2*std::distance(m_FCal_ChannelMap.begin(imap), m_FCal_ChannelMap.end(imap))){
      cout << "Error: Incompatibility between FCalChannel map and GEO file: Different number of cells in m_cells_in_sampling and FCal_ChannelMap" << endl;
      cout << "m_cells_in_sampling: " << m_cells_in_sampling[sampling].size() << endl;
      cout << "FCal_ChannelMap: " << 2*std::distance(m_FCal_ChannelMap.begin(imap), m_FCal_ChannelMap.end(imap)) << endl;
      return false;
    }
 
    for(auto it=m_FCal_ChannelMap.begin(imap);it!=m_FCal_ChannelMap.end(imap);it++){
 
 
      phi_index = it->first & 0xffff;
      eta_index = it->first >> 16;
      x=it->second.x();
      y=it->second.y();
      m_FCal_ChannelMap.tileSize(imap, eta_index, phi_index,dx,dy);
 
      id=(mask1[imap-1]<<12) + (eta_index << 5) +2*phi_index;
 
      if(imap==2) id+= (8<<8);
 
      Identifier id1((unsigned long long)(id<<44));
      const CaloDetDescrElement *DDE1 =getDDE(id1);
 
      id=(mask2[imap-1]<<12) + (eta_index << 5) +2*phi_index;
      if(imap==2) id+= (8<<8);
        Identifier id2((unsigned long long)(id<<44));
      const CaloDetDescrElement *DDE2=getDDE(id2);
 
      if(!TMath::AreEqualRel(x, DDE1->x(),1.E-8) || !TMath::AreEqualRel(y, DDE1->y(),1.E-8) || !TMath::AreEqualRel(x, DDE2->x(),1.E-8) || !TMath::AreEqualRel(y, DDE2->y(),1.E-8) ){
     cout << "Error: Incompatibility between FCalChannel map and GEO file \n" << x << " " << DDE1->x() << " " << DDE2->x() << y << " " << DDE1->y() << " " << DDE2->y() << endl;
     return false;   
      }
    }
 
 
 
  }
   
  return true;
}

deta()

double CaloGeometry::deta ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 752 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  double mineta=m_min_eta_sample[side][sample];
  double maxeta=m_max_eta_sample[side][sample];
 
  if(eta<mineta)
  {
    return fabs(eta-mineta);
  }
  else if(eta>maxeta)
  {
   return fabs(eta-maxeta);
    }
    else
    {
   double d1=fabs(eta-mineta);
   double d2=fabs(eta-maxeta);
   if(d1<d2) return -d1;
   else return -d2;
  }
}

DoGraphs()

bool CaloGeometry::DoGraphs ( ) const

inline

Definition at line 67 of file CaloGeometry.h.

{return m_dographs;};

DrawGeoForPhi0()

TCanvas * CaloGeometry::DrawGeoForPhi0

(

)

Definition at line 440 of file CaloGeometry.cxx.

{
  TCanvas* c=new TCanvas("CaloGeoForPhi0","Calo geometry for #phi~0");
  TH2D* hcalolayout=new TH2D("hcalolayoutPhi0","Reconstruction geometry: calorimeter layout;z [mm];r [mm]",50,-7000,7000,50,0,4000);
  hcalolayout->Draw();
  hcalolayout->SetStats(0);
  hcalolayout->GetYaxis()->SetTitleOffset(1.4);
  
  TLegend* leg=new TLegend(0.30,0.13,0.70,0.37);
  leg->SetFillStyle(0);
  leg->SetFillColor(10);
  leg->SetBorderSize(1);
  leg->SetNColumns(2);
 
  for(int sample=0;sample<MAX_SAMPLING;++sample) {
    TGraph* gr=DrawGeoSampleForPhi0(sample,CaloGeometry_calocol[sample],(sample==17));
    if(gr) {
      std::string sname=Form("Sampling %2d : ",sample);
      sname+=SamplingName(sample);
      leg->AddEntry(gr,sname.c_str(),"LF");
    }
  }
  leg->Draw();
  return c;
}

DrawGeoSampleForPhi0()

TGraph * CaloGeometry::DrawGeoSampleForPhi0	(	int	sample,
		int	calocol,
		bool	print = false )

Definition at line 349 of file CaloGeometry.cxx.

{
  if (sample < CaloSampling::PreSamplerB || sample >= CaloSampling::Unknown) {
    return nullptr;
  }
  TGraph* firstgr=nullptr;
  cout<<"Start sample "<<sample<<" ("<<SamplingName(sample)<<")"<<endl;
  int ngr=0;
  for(t_eta_cellmap::iterator calo_iter=m_cells_in_sampling_for_phi0[sample].begin();calo_iter!=m_cells_in_sampling_for_phi0[sample].end();++calo_iter) {
    const CaloDetDescrElement* theDDE=(*calo_iter).second;
    if(theDDE) {
      TVector3 cv;
      TGraph* gr=new TGraph(5);
      gr->SetLineColor(TMath::Abs(calocol));
      gr->SetFillColor(TMath::Abs(calocol));
      if(calocol<0) {
        gr->SetFillStyle(1001);
      } else {
        gr->SetFillStyle(0);
      }
      gr->SetLineWidth(2);
      double r=theDDE->r();
      double dr=theDDE->dr();
      double x=theDDE->x();
      double dx=theDDE->dx();
      double y=theDDE->y();
      double dy=theDDE->dy();
      double z=theDDE->z();
      double dz=theDDE->dz()*2;
      double eta=theDDE->eta();
      double deta=theDDE->deta();
      
      if(CaloSampling::PreSamplerB<=sample && sample<=CaloSampling::EMB3) {
       dr*=2;
      }
      if(print) {
        cout<<"sample="<<sample<<" r="<<r<<" dr="<<dr<<" eta="<<eta<<" deta="<<deta<<" x="<<x<<" y="<<y<<" z="<<z<<" dz="<<dz<<endl;
      }
      if(isCaloBarrel(sample)) {
        cv.SetPtEtaPhi(r-dr/2,eta-deta/2,0);
        gr->SetPoint(0,cv.Z(),cv.Pt());
        gr->SetPoint(4,cv.Z(),cv.Pt());
        cv.SetPtEtaPhi(r-dr/2,eta+deta/2,0);
        gr->SetPoint(1,cv.Z(),cv.Pt());
        cv.SetPtEtaPhi(r+dr/2,eta+deta/2,0);
        gr->SetPoint(2,cv.Z(),cv.Pt());
        cv.SetPtEtaPhi(r+dr/2,eta-deta/2,0);
        gr->SetPoint(3,cv.Z(),cv.Pt());
      } else {
        if(sample<CaloSampling::FCAL0) {
          cv.SetPtEtaPhi(1,eta-deta/2,0);cv*=(z-dz/2)/cv.Z();
          gr->SetPoint(0,cv.Z(),cv.Pt());
          gr->SetPoint(4,cv.Z(),cv.Pt());
          cv.SetPtEtaPhi(1,eta+deta/2,0);cv*=(z-dz/2)/cv.Z();
          gr->SetPoint(1,cv.Z(),cv.Pt());
          cv.SetPtEtaPhi(1,eta+deta/2,0);cv*=(z+dz/2)/cv.Z();
          gr->SetPoint(2,cv.Z(),cv.Pt());
          cv.SetPtEtaPhi(1,eta-deta/2,0);cv*=(z+dz/2)/cv.Z();
          gr->SetPoint(3,cv.Z(),cv.Pt());
        } else {
          double minr=r;
          double maxr=r;
          for(double px=x-dx/2;px<=x+dx/2;px+=dx) {
            for(double py=y-dy/2;py<=y+dy/2;py+=dy) {
              double pr=TMath::Sqrt(px*px+py*py);
              minr=TMath::Min(minr,pr);
              maxr=TMath::Max(maxr,pr);
            }
          }
          cv.SetXYZ(minr,0,z-dz/2);
          gr->SetPoint(0,cv.Z(),cv.Pt());
          gr->SetPoint(4,cv.Z(),cv.Pt());
          cv.SetXYZ(maxr,0,z-dz/2);
          gr->SetPoint(1,cv.Z(),cv.Pt());
          cv.SetXYZ(maxr,0,z+dz/2);
          gr->SetPoint(2,cv.Z(),cv.Pt());
          cv.SetXYZ(minr,0,z+dz/2);
          gr->SetPoint(3,cv.Z(),cv.Pt());
        }  
      }
      //if(calocol[sample]>0) gr->Draw("Lsame");
      // else gr->Draw("LFsame");
      gr->Draw("LFsame");
      if(ngr==0) firstgr=gr;
      ++ngr;
    }
  }  
  cout<<"Done sample "<<sample<<" ("<<SamplingName(sample)<<")="<<ngr<<endl;
  return firstgr;
}

getClosestFCalCellIndex()

bool CaloGeometry::getClosestFCalCellIndex	(	int	sampling,
		float	x,
		float	y,
		int &	ieta,
		int &	iphi,
		int *	steps = 0 ) const

Definition at line 592 of file CaloGeometry.cxx.

                                                                                                             {
  
  double rmin = m_FCal_rmin.at(sampling-21);
  double rmax = m_FCal_rmax.at(sampling-21);
  int isam=sampling-20;
  double a=1.;
  const double b=0.01;
  const int nmax=100;
  int i=0;
  
  const double r = sqrt(x*x +y*y);
  if(r==0.) return false;
  const double r_inverse=1./r;
  
  if((r/rmax)>(rmin*r_inverse)){
    x=x*rmax*r_inverse;
    y=y*rmax*r_inverse;
    while((!m_FCal_ChannelMap.getTileID(isam, a*x, a*y, ieta, iphi)) && i<nmax){
      a-=b;
      i++;
    }
  }
  else {
    x=x*rmin*r_inverse;
    y=y*rmin*r_inverse;
    while((!m_FCal_ChannelMap.getTileID(isam, a*x, a*y, ieta, iphi)) && i<nmax){
      a+=b;
      i++;
    }
    
  }
  if(steps)*steps=i+1;
  return i<nmax;
}

getDDE() [1/3]

const CaloDetDescrElement * CaloGeometry::getDDE ( Identifier identify ) const

virtual

Implements ICaloGeometry.

Definition at line 466 of file CaloGeometry.cxx.

{
  return m_cells.at(identify);
}

getDDE() [2/3]

const CaloDetDescrElement * CaloGeometry::getDDE ( int sampling, float eta, float phi, float * distance = 0, int * steps = 0 ) const

virtual

Implements ICaloGeometry.

Definition at line 475 of file CaloGeometry.cxx.

{
  if(sampling<0) return nullptr;
  if(sampling>=MAX_SAMPLING) return nullptr;
  if(m_cells_in_regions.at(sampling).empty()) return nullptr;
  
  float dist = 0;
  const CaloDetDescrElement* bestDDE=nullptr;
  if(!distance) distance=&dist;
  *distance=+10000000;
  int intsteps;
  int beststeps;
  if(steps) beststeps=(*steps);
   else beststeps=0;
  
  if(sampling<21) {
    for(int skip_range_check=0;skip_range_check<=1;++skip_range_check) {
      for(unsigned int j=0;j<m_cells_in_regions.at(sampling).size();++j) {
        if(!skip_range_check) {
          if(eta<m_cells_in_regions.at(sampling).at(j)->mineta()) continue;
          if(eta>m_cells_in_regions.at(sampling).at(j)->maxeta()) continue;
        }  
        if(steps) intsteps=(*steps);
         else intsteps=0;
        if(m_debug) {
          cout<<"CaloGeometry::getDDE : check map"<<j<<" skip_range_check="<<skip_range_check<<endl;
        }
        float newdist;
        const CaloDetDescrElement* newDDE=m_cells_in_regions.at(sampling).at(j)->getDDE(eta,phi,&newdist,&intsteps);
        if(m_debug) {
          cout<<"CaloGeometry::getDDE : map"<<j<<" dist="<<newdist<<" best dist="<<*distance<<" steps="<<intsteps<<endl;
        }
        if(newdist<*distance) {
          bestDDE=newDDE;
          *distance=newdist;
          if(steps) beststeps=intsteps;
          if(newdist<-0.1) break; //stop, we are well within the hit cell
        }
      }
      if(bestDDE) break;
    }
  } else {
        return nullptr;
    //for(int skip_range_check=0;skip_range_check<=1;++skip_range_check) {
      //for(unsigned int j=0;j<m_cells_in_regions[sampling].size();++j) {
        //if(!skip_range_check) {
          //if(eta<m_cells_in_regions[sampling][j]->minx()) continue;
          //if(eta>m_cells_in_regions[sampling][j]->maxx()) continue;
          //if(phi<m_cells_in_regions[sampling][j]->miny()) continue;
          //if(phi>m_cells_in_regions[sampling][j]->maxy()) continue;
        //}  
        //if(steps) intsteps=*steps;
         //else intsteps=0;
        //if(m_debug) {
          //cout<<"CaloGeometry::getDDE : check map"<<j<<" skip_range_check="<<skip_range_check<<endl;
        //}
        //float newdist;
        //const CaloGeoDetDescrElement* newDDE=m_cells_in_regions[sampling][j]->getDDE(eta,phi,&newdist,&intsteps);
        //if(m_debug) {
          //cout<<"CaloGeometry::getDDE : map"<<j<<" dist="<<newdist<<" best dist="<<*distance<<" steps="<<intsteps<<endl;
        //}
        //if(newdist<*distance) {
          //bestDDE=newDDE;
          //*distance=newdist;
          //if(steps) beststeps=intsteps;
          //if(newdist<-0.1) break; //stop, we are well within the hit cell
        //}
      //}
      //if(bestDDE) break;
    //}  
  }
  if(steps) *steps=beststeps;
  return bestDDE;
}

getDDE() [3/3]

const CaloDetDescrElement * CaloGeometry::getDDE ( int sampling, Identifier identify ) const

virtual

Definition at line 470 of file CaloGeometry.cxx.

{
  return m_cells_in_sampling.at(sampling).at(identify);
}

GetFCAL_ChannelMap()

FCAL_ChannelMap * CaloGeometry::GetFCAL_ChannelMap ( )

inline

Definition at line 72 of file CaloGeometry.h.

{return &m_FCal_ChannelMap;}

getFCalDDE()

const CaloDetDescrElement * CaloGeometry::getFCalDDE ( int sampling, float x, float y, float z, float * distance = 0, int * steps = 0 ) const

virtual

Implements ICaloGeometry.

Definition at line 550 of file CaloGeometry.cxx.

                                                                                                                        {
  int isam = sampling - 20;
  int iphi(-100000),ieta(-100000);
  Long64_t mask1[]{0x34,0x34,0x35};
  Long64_t mask2[]{0x36,0x36,0x37};
  bool found = m_FCal_ChannelMap.getTileID(isam, x, y, ieta, iphi);
  if(steps && found) *steps=0;
  if(!found) {
    //cout << "Warning: Hit is not matched with any FCal cell! Looking for the closest cell" << endl;
    found = getClosestFCalCellIndex(sampling, x, y, ieta, iphi,steps);
  }
  if(!found) {
    cout << "Error: Unable to find the closest FCal cell!" << endl;
    return nullptr;
  }
  
  
  //cout << "CaloGeometry::getFCalDDE: x:" << x << " y: " << y << " ieta: " << ieta << " iphi: " << iphi << " cmap->x(): " << m_FCal_ChannelMap.x(isam,ieta,iphi) << " cmap->y(): " << m_FCal_ChannelMap.y(isam,ieta,iphi) << endl;
  Long64_t id = (ieta << 5) + 2*iphi;
  if(isam==2)id+= (8<<8);
  
  
  
  if(z>0) id+=(mask2[isam-1] << 12);
  else id+=(mask1[isam-1] << 12);
  
  id = id << 44; 
  Identifier identify((unsigned long long)id);
  const CaloDetDescrElement* foundcell=m_cells.at(identify);
 
  // Report the shortest distance in x or y to the face of the cell if outside the cell.
  // If inside the cell this will be negative and be the shortest distance to a cell face.
  float dist = 0;
  if(!distance) distance=&dist;
  float distanceX = abs(foundcell->x()-x) - foundcell->dx()/2;
  float distanceY = abs(foundcell->y()-y) - foundcell->dy()/2;
  *distance = max(distanceX,distanceY);
  
  return foundcell;
}

GetGraph()

TGraphErrors * CaloGeometry::GetGraph ( unsigned int sample ) const

inline

Definition at line 65 of file CaloGeometry.h.

{return m_graph_layers[sample];};

InitRZmaps()

void CaloGeometry::InitRZmaps ( )

protectedvirtual

Definition at line 197 of file CaloGeometry.cxx.

{
  FSmap< double , double > rz_map_eta [2][MAX_SAMPLING];
  FSmap< double , double > rz_map_rmid[2][MAX_SAMPLING];
  FSmap< double , double > rz_map_zmid[2][MAX_SAMPLING];
  FSmap< double , double > rz_map_rent[2][MAX_SAMPLING];
  FSmap< double , double > rz_map_zent[2][MAX_SAMPLING];
  FSmap< double , double > rz_map_rext[2][MAX_SAMPLING];
  FSmap< double , double > rz_map_zext[2][MAX_SAMPLING];
  FSmap< double , int    > rz_map_n   [2][MAX_SAMPLING];
 
 
  for(unsigned int side=0;side<=1;++side) for(unsigned int sample=0;sample<MAX_SAMPLING;++sample)
  {
    m_min_eta_sample[side][sample]=+1000;
    m_max_eta_sample[side][sample]=-1000;
  }  
  
  
  for(t_cellmap::iterator calo_iter=m_cells.begin();calo_iter!=m_cells.end();++calo_iter)
  {
    const CaloDetDescrElement* theDDE=(*calo_iter).second;
    if(theDDE)
    {
      unsigned int sample=theDDE->getSampling();
 
      int side=0;
      int sign_side=-1;
      double eta_raw=theDDE->eta_raw();
      if(eta_raw>0) {
        side=1;
        sign_side=+1;
      }
      
      if(!m_cells_in_sampling_for_phi0[sample][eta_raw]) {
        m_cells_in_sampling_for_phi0[sample][eta_raw]=theDDE;
      } else {
        if(TMath::Abs(theDDE->phi()) < TMath::Abs(m_cells_in_sampling_for_phi0[sample][eta_raw]->phi())) {
          m_cells_in_sampling_for_phi0[sample][eta_raw]=theDDE;
        }
      }  
      
      double min_eta=theDDE->eta()-theDDE->deta()/2;
      double max_eta=theDDE->eta()+theDDE->deta()/2;
      if(min_eta<m_min_eta_sample[side][sample]) m_min_eta_sample[side][sample]=min_eta;
      if(max_eta>m_max_eta_sample[side][sample]) m_max_eta_sample[side][sample]=max_eta;
      
      if(rz_map_eta[side][sample].find(eta_raw)==rz_map_eta[side][sample].end()) {
        rz_map_eta [side][sample][eta_raw]=0;
        rz_map_rmid[side][sample][eta_raw]=0;
        rz_map_zmid[side][sample][eta_raw]=0;
        rz_map_rent[side][sample][eta_raw]=0;
        rz_map_zent[side][sample][eta_raw]=0;
        rz_map_rext[side][sample][eta_raw]=0;
        rz_map_zext[side][sample][eta_raw]=0;
        rz_map_n   [side][sample][eta_raw]=0;
      }
      rz_map_eta [side][sample][eta_raw]+=theDDE->eta();
      rz_map_rmid[side][sample][eta_raw]+=theDDE->r();
      rz_map_zmid[side][sample][eta_raw]+=theDDE->z();
      double drh=theDDE->dr()/2;
      double dzh=theDDE->dz();
      if(sample<=CaloSampling::EMB3) { // ensure we have a valid sampling
        drh=theDDE->dr();
      } 
      // An `else` would be good here since we can't continue with a *bad* sampling...
      // Should most likely handle the situation nicely rather than with a crash.
 
      rz_map_rent[side][sample][eta_raw]+=theDDE->r()-drh;
      rz_map_zent[side][sample][eta_raw]+=theDDE->z()-dzh*sign_side;
      rz_map_rext[side][sample][eta_raw]+=theDDE->r()+drh;
      rz_map_zext[side][sample][eta_raw]+=theDDE->z()+dzh*sign_side;
      rz_map_n   [side][sample][eta_raw]++;
      
    }
  }
  
  for(int side=0;side<=1;++side)
  {
   for(int sample=0;sample<MAX_SAMPLING;++sample)
   {
 
    if(!rz_map_n[side][sample].empty())
    {
      for(FSmap< double , int >::iterator iter=rz_map_n[side][sample].begin();iter!=rz_map_n[side][sample].end();++iter)
      {
        double eta_raw=iter->first;
        if(iter->second<1)
        {
          //ATH_MSG_WARNING("rz-map for side="<<side<<" sample="<<sample<<" eta_raw="<<eta_raw<<" : #cells="<<iter->second<<" !!!");
        }
        else
        {
          double eta =rz_map_eta[side][sample][eta_raw]/iter->second;
          double rmid=rz_map_rmid[side][sample][eta_raw]/iter->second;
          double zmid=rz_map_zmid[side][sample][eta_raw]/iter->second;
          double rent=rz_map_rent[side][sample][eta_raw]/iter->second;
          double zent=rz_map_zent[side][sample][eta_raw]/iter->second;
          double rext=rz_map_rext[side][sample][eta_raw]/iter->second;
          double zext=rz_map_zext[side][sample][eta_raw]/iter->second;
          
          m_rmid_map[side][sample][eta]=rmid;
          m_zmid_map[side][sample][eta]=zmid;
          m_rent_map[side][sample][eta]=rent;
          m_zent_map[side][sample][eta]=zent;
          m_rext_map[side][sample][eta]=rext;
          m_zext_map[side][sample][eta]=zext;
        }
      }
      //ATH_MSG_DEBUG("rz-map for side="<<side<<" sample="<<sample<<" #etas="<<m_rmid_map[side][sample].size());
     }
    else
    {
     std::cout<<"rz-map for side="<<side<<" sample="<<sample<<" is empty!!!"<<std::endl;
    }
   } //sample
  } //side
  
  if(DoGraphs()) {
    for(int sample=0;sample<MAX_SAMPLING;++sample) {
      m_graph_layers[sample]=new TGraphErrors(rz_map_n[0][sample].size()+rz_map_n[1][sample].size());
      m_graph_layers[sample]->SetMarkerColor(TMath::Abs(CaloGeometry_calocol[sample]));
      m_graph_layers[sample]->SetLineColor(TMath::Abs(CaloGeometry_calocol[sample]));
      int np=0;
      for(int side=0;side<=1;++side) {
        for(FSmap< double , int >::iterator iter=rz_map_n[side][sample].begin();iter!=rz_map_n[side][sample].end();++iter) {
          double eta_raw=iter->first;
          int sign_side=-1;
          if(eta_raw>0) sign_side=+1;
          //double eta =rz_map_eta[side][sample][eta_raw]/iter->second;
          double rmid=rz_map_rmid[side][sample][eta_raw]/iter->second;
          double zmid=rz_map_zmid[side][sample][eta_raw]/iter->second;
          //double rent=rz_map_rent[side][sample][eta_raw]/iter->second;
          //double zent=rz_map_zent[side][sample][eta_raw]/iter->second;
          double rext=rz_map_rext[side][sample][eta_raw]/iter->second;
          double zext=rz_map_zext[side][sample][eta_raw]/iter->second;
          m_graph_layers[sample]->SetPoint(np,zmid,rmid);
          /*
          if(isCaloBarrel(sample)) {
            m_graph_layers[sample]->SetPointError(np,0,rext-rmid);
          } else {
            m_graph_layers[sample]->SetPointError(np,(zext-zent)*sign_side,0);
          }
          */
          m_graph_layers[sample]->SetPointError(np,(zext-zmid)*sign_side,rext-rmid);
          ++np;
        }
      }  
    }
  }
}

isCaloBarrel()

bool CaloGeometry::isCaloBarrel ( int sample ) const

inlinevirtual

Implements ICaloGeometry.

Definition at line 62 of file CaloGeometry.h.

{return m_isCaloBarrel[sample];};

minmaxeta()

void CaloGeometry::minmaxeta ( int sample, double eta, double & mineta, double & maxeta ) const

virtual

Implements ICaloGeometry.

Definition at line 778 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
  
  mineta=m_min_eta_sample[side][sample];
  maxeta=m_max_eta_sample[side][sample];
}

post_process()

void CaloGeometry::post_process ( int layer )

protectedvirtual

Definition at line 121 of file CaloGeometry.cxx.

{
  //cout<<"post processing sampling "<<sampling<<endl;
  bool found_overlap=false;
  for(unsigned int j=0;j<m_cells_in_regions[sampling].size();++j) {
    /*
    cout<<"Sample "<<sampling<<": checking map "<<j<<"/"<<m_cells_in_regions[sampling].size();
    for(unsigned int k=0;k<m_cells_in_regions[sampling].size();++k) {
      cout<<", "<<k<<":"<<m_cells_in_regions[sampling][k]->size();
    }
    cout<<endl;
    */
    for(unsigned int i=j+1;i<m_cells_in_regions[sampling].size();++i) {
      if(m_cells_in_regions[sampling][i]->has_overlap(m_cells_in_regions[sampling][j])) {
        if(!found_overlap) {
          cout<<"Sample "<<sampling<<", starting maps : "<<m_cells_in_regions[sampling].size();
          for(unsigned int k=0;k<m_cells_in_regions[sampling].size();++k) {
            cout<<", "<<k<<":"<<m_cells_in_regions[sampling][k]->size();
          }
          cout<<endl;
        }
        found_overlap=true;
        /*
        cout<<"Sample "<<sampling<<": Found overlap between map "<<j<<" and "<<i<<", "
            <<m_cells_in_regions[sampling].size()<<" total maps"<<endl;
        cout<<"  current configuration map "<<j<<"/"<<m_cells_in_regions[sampling].size();
        for(unsigned int k=0;k<m_cells_in_regions[sampling].size();++k) {
          cout<<", "<<k<<":"<<m_cells_in_regions[sampling][k]->size();
        }
        cout<<endl;
 
        PrintMapInfo(sampling,j);
        PrintMapInfo(sampling,i);
        */
        
        CaloGeometryLookup* toremove=m_cells_in_regions[sampling][i];
        toremove->merge_into_ref(m_cells_in_regions[sampling][j]);
        
        /*
        cout<<"  New ";
        PrintMapInfo(sampling,j);
        */
        
        for(unsigned int k=i;k<m_cells_in_regions[sampling].size()-1;++k) {
          m_cells_in_regions[sampling][k]=m_cells_in_regions[sampling][k+1];
          m_cells_in_regions[sampling][k]->set_index(k);
        }
        m_cells_in_regions[sampling].resize(m_cells_in_regions[sampling].size()-1);
 
        /*
        cout<<"  new configuration map "<<j<<"/"<<m_cells_in_regions[sampling].size();
        for(unsigned int k=0;k<m_cells_in_regions[sampling].size();++k) {
          cout<<", "<<k<<":"<<m_cells_in_regions[sampling][k]->size();
        }
        cout<<endl;
        */
        
        --i;
      }
    }
  }
  
  if(found_overlap) {
    cout<<"Sample "<<sampling<<", final maps : "<<m_cells_in_regions[sampling].size();
    for(unsigned int k=0;k<m_cells_in_regions[sampling].size();++k) {
      cout<<", "<<k<<":"<<m_cells_in_regions[sampling][k]->size();
    }
    cout<<endl;
  }
 
  if(found_overlap) {
    cout<<"Run another round of ";
    post_process(sampling);
  }  
}

PostProcessGeometry()

bool CaloGeometry::PostProcessGeometry ( )

virtual

Implements ICaloGeometry.

Definition at line 627 of file CaloGeometry.cxx.

{
  for(int i=0;i<MAX_SAMPLING;++i) {
    post_process(i);
    for(unsigned int j=0;j<m_cells_in_regions[i].size();++j) {
      m_cells_in_regions[i][j]->post_process();
    }
    //if(i>=21) break;
  }
  
  InitRZmaps(); 
  
  /*
  cout<<"all : "<<m_cells.size()<<endl;
  for(int sampling=0;sampling<MAX_SAMPLING;++sampling) {
    cout<<"cells sampling "<<sampling<<" : "<<m_cells_in_sampling[sampling].size()<<" cells";
    cout<<", "<<m_cells_in_regions[sampling].size()<<" lookup map(s)"<<endl;
    for(unsigned int j=0;j<m_cells_in_regions[sampling].size();++j) {
      PrintMapInfo(sampling,j);
      //break;
    }
    //if(i>0) break;
  } 
  */ 
 
  return true;
}

PrintMapInfo()

void CaloGeometry::PrintMapInfo ( int i, int j )

virtual

Definition at line 93 of file CaloGeometry.cxx.

{
  cout<<"  map "<<j<<": "<<m_cells_in_regions[i][j]->size()<<" cells";
  if(i<21) {
    cout<<", "<<m_cells_in_regions[i][j]->cell_grid_eta()<<"*"<<m_cells_in_regions[i][j]->cell_grid_phi();
    cout<<", deta="<<m_cells_in_regions[i][j]->deta().mean()<<"+-"<<m_cells_in_regions[i][j]->deta().rms();
    cout<<", dphi="<<m_cells_in_regions[i][j]->dphi().mean()<<"+-"<<m_cells_in_regions[i][j]->dphi().rms();
    cout<<", mineta="<<m_cells_in_regions[i][j]->mineta()<<", maxeta="<<m_cells_in_regions[i][j]->maxeta();
    cout<<", minphi="<<m_cells_in_regions[i][j]->minphi()<<", maxphi="<<m_cells_in_regions[i][j]->maxphi();
    cout<<endl<<"         ";  
    cout<<", mineta_raw="<<m_cells_in_regions[i][j]->mineta_raw()<<", maxeta_raw="<<m_cells_in_regions[i][j]->maxeta_raw();
    cout<<", minphi_raw="<<m_cells_in_regions[i][j]->minphi_raw()<<", maxphi_raw="<<m_cells_in_regions[i][j]->maxphi_raw();
    cout<<endl;
  } else {
    cout<<", "<<m_cells_in_regions[i][j]->cell_grid_eta()<<"*"<<m_cells_in_regions[i][j]->cell_grid_phi();
    cout<<", dx="<<m_cells_in_regions[i][j]->dx().mean()<<"+-"<<m_cells_in_regions[i][j]->dx().rms();
    cout<<", dy="<<m_cells_in_regions[i][j]->dy().mean()<<"+-"<<m_cells_in_regions[i][j]->dy().rms();
    cout<<", mindx="<<m_cells_in_regions[i][j]->mindx();
    cout<<", mindy="<<m_cells_in_regions[i][j]->mindy();
    cout<<", minx="<<m_cells_in_regions[i][j]->minx()<<", maxx="<<m_cells_in_regions[i][j]->maxx();
    cout<<", miny="<<m_cells_in_regions[i][j]->miny()<<", maxy="<<m_cells_in_regions[i][j]->maxy();
    cout<<endl<<"         ";  
    cout<<", minx_raw="<<m_cells_in_regions[i][j]->minx_raw()<<", maxx_raw="<<m_cells_in_regions[i][j]->maxx_raw();
    cout<<", miny_raw="<<m_cells_in_regions[i][j]->miny_raw()<<", maxy_raw="<<m_cells_in_regions[i][j]->maxy_raw();
    cout<<endl;
  }
}

rent()

double CaloGeometry::rent ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 812 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  return m_rent_map[side][sample].find_closest(eta)->second;
}

rext()

double CaloGeometry::rext ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 837 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  return m_rext_map[side][sample].find_closest(eta)->second;
}

rmid()

double CaloGeometry::rmid ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 787 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  return m_rmid_map[side][sample].find_closest(eta)->second;
}

rpos()

double CaloGeometry::rpos ( int sample, double eta, int subpos = CaloSubPos::SUBPOS_MID ) const

virtual

Implements ICaloGeometry.

Definition at line 862 of file CaloGeometry.cxx.

{
    
  int side=0;
  if(eta>0) side=1;
   
  if(subpos==CaloSubPos::SUBPOS_ENT) return m_rent_map[side][sample].find_closest(eta)->second;
  if(subpos==CaloSubPos::SUBPOS_EXT) return m_rext_map[side][sample].find_closest(eta)->second;
    
  return m_rmid_map[side][sample].find_closest(eta)->second;
}

rzent()

double CaloGeometry::rzent ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 828 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  if(isCaloBarrel(sample)) return m_rent_map[side][sample].find_closest(eta)->second;
   else                    return m_zent_map[side][sample].find_closest(eta)->second;
}

rzext()

double CaloGeometry::rzext ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 853 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  if(isCaloBarrel(sample)) return m_rext_map[side][sample].find_closest(eta)->second;
   else                    return m_zext_map[side][sample].find_closest(eta)->second;
}

rzmid()

double CaloGeometry::rzmid ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 803 of file CaloGeometry.cxx.

{
 int side=0;
 if(eta>0) side=1;
    
 if(isCaloBarrel(sample)) return m_rmid_map[side][sample].find_closest(eta)->second;
 else                     return m_zmid_map[side][sample].find_closest(eta)->second;
}

rzpos()

double CaloGeometry::rzpos ( int sample, double eta, int subpos = CaloSubPos::SUBPOS_MID ) const

virtual

Implements ICaloGeometry.

Definition at line 884 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  if(isCaloBarrel(sample)) {
    if(subpos==CaloSubPos::SUBPOS_ENT) return m_rent_map[side][sample].find_closest(eta)->second;
    if(subpos==CaloSubPos::SUBPOS_EXT) return m_rext_map[side][sample].find_closest(eta)->second;
    return m_rmid_map[side][sample].find_closest(eta)->second;
  } else {
    if(subpos==CaloSubPos::SUBPOS_ENT) return m_zent_map[side][sample].find_closest(eta)->second;
    if(subpos==CaloSubPos::SUBPOS_EXT) return m_zext_map[side][sample].find_closest(eta)->second;
    return m_zmid_map[side][sample].find_closest(eta)->second;
  }  
}

SamplingName()

std::string CaloGeometry::SamplingName ( int sample )

static

Definition at line 900 of file CaloGeometry.cxx.

{
  return CaloSampling::getSamplingName(sample);
}

SetDoGraphs()

void CaloGeometry::SetDoGraphs ( bool dographs = true )

inline

Definition at line 66 of file CaloGeometry.h.

{m_dographs=dographs;};

SetFCal_ChannelMap()

void CaloGeometry::SetFCal_ChannelMap ( const FCAL_ChannelMap * fcal_ChannnelMap )

inline

Definition at line 73 of file CaloGeometry.h.

{m_FCal_ChannelMap=*fcal_ChannnelMap;}

zent()

double CaloGeometry::zent ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 820 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  return m_zent_map[side][sample].find_closest(eta)->second;
}

zext()

double CaloGeometry::zext ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 845 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  return m_zext_map[side][sample].find_closest(eta)->second;
}

zmid()

double CaloGeometry::zmid ( int sample, double eta ) const

virtual

Implements ICaloGeometry.

Definition at line 795 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  return m_zmid_map[side][sample].find_closest(eta)->second;
}

zpos()

double CaloGeometry::zpos ( int sample, double eta, int subpos = CaloSubPos::SUBPOS_MID ) const

virtual

Implements ICaloGeometry.

Definition at line 874 of file CaloGeometry.cxx.

{
  int side=0;
  if(eta>0) side=1;
 
  if(subpos==CaloSubPos::SUBPOS_ENT) return m_zent_map[side][sample].find_closest(eta)->second;
  if(subpos==CaloSubPos::SUBPOS_EXT) return m_zext_map[side][sample].find_closest(eta)->second;
  return m_zmid_map[side][sample].find_closest(eta)->second;
}

Member Data Documentation

m_cells

t_cellmap CaloGeometry::m_cells

protected

Definition at line 100 of file CaloGeometry.h.

m_cells_in_regions

std::vector< std::vector< CaloGeometryLookup* > > CaloGeometry::m_cells_in_regions

protected

Definition at line 103 of file CaloGeometry.h.

m_cells_in_sampling

std::vector< t_cellmap > CaloGeometry::m_cells_in_sampling

protected

Definition at line 101 of file CaloGeometry.h.

m_cells_in_sampling_for_phi0

std::vector< t_eta_cellmap > CaloGeometry::m_cells_in_sampling_for_phi0

protected

Definition at line 102 of file CaloGeometry.h.

m_debug

std::atomic< bool > CaloGeometry::m_debug =false

static

Definition at line 32 of file CaloGeometry.h.

m_debug_identify

const Identifier CaloGeometry::m_debug_identify

static

Definition at line 31 of file CaloGeometry.h.

m_dographs

bool CaloGeometry::m_dographs

protected

Definition at line 115 of file CaloGeometry.h.

m_FCal_ChannelMap

FCAL_ChannelMap CaloGeometry::m_FCal_ChannelMap

protected

Definition at line 117 of file CaloGeometry.h.

m_FCal_rmax

std::vector<double> CaloGeometry::m_FCal_rmax

protected

Definition at line 118 of file CaloGeometry.h.

m_FCal_rmin

std::vector<double> CaloGeometry::m_FCal_rmin

protected

Definition at line 118 of file CaloGeometry.h.

m_graph_layers

std::vector< TGraphErrors* > CaloGeometry::m_graph_layers

protected

Definition at line 116 of file CaloGeometry.h.

m_isCaloBarrel

std::vector< bool > CaloGeometry::m_isCaloBarrel

protected

Definition at line 105 of file CaloGeometry.h.

m_max_eta_sample

std::vector< double > CaloGeometry::m_max_eta_sample[2]

protected

Definition at line 107 of file CaloGeometry.h.

m_min_eta_sample

std::vector< double > CaloGeometry::m_min_eta_sample[2]

protected

Definition at line 106 of file CaloGeometry.h.

m_rent_map

std::vector< FSmap< double , double > > CaloGeometry::m_rent_map[2]

protected

Definition at line 110 of file CaloGeometry.h.

m_rext_map

std::vector< FSmap< double , double > > CaloGeometry::m_rext_map[2]

protected

Definition at line 112 of file CaloGeometry.h.

m_rmid_map

std::vector< FSmap< double , double > > CaloGeometry::m_rmid_map[2]

protected

Definition at line 108 of file CaloGeometry.h.

m_zent_map

std::vector< FSmap< double , double > > CaloGeometry::m_zent_map[2]

protected

Definition at line 111 of file CaloGeometry.h.

m_zext_map

std::vector< FSmap< double , double > > CaloGeometry::m_zext_map[2]

protected

Definition at line 113 of file CaloGeometry.h.

m_zmid_map

std::vector< FSmap< double , double > > CaloGeometry::m_zmid_map[2]

protected

Definition at line 109 of file CaloGeometry.h.

MAX_SAMPLING

const int CaloGeometry::MAX_SAMPLING = CaloCell_ID_FCS::MaxSample

static

Definition at line 29 of file CaloGeometry.h.

---

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

• CaloGeometry.h
• CaloGeometry.cxx