CaloVertexedCell.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id$
 
 
#include "CaloUtils/CaloVertexedCell.h"
#include "CaloEvent/CaloCell.h"
 

CaloVertexedCell::CaloVertexedCell (const CaloCell& cell)
  : m_cell (&cell),
    m_eta (cell.eta()),
    m_phi (cell.phi()),
    m_sinth (cell.sinTh()),
    m_costh (cell.cosTh()),
    m_sinph (cell.sinPhi()),
    m_cosph (cell.cosPhi())
{
}
 

CaloVertexedCell::CaloVertexedCell (const CaloCell& cell,
                                    const Amg::Vector3D&  vx)
  : m_cell (&cell),
    m_eta (cell.eta()),
    m_phi (cell.phi()),
    m_sinth (cell.sinTh()),
    m_costh (cell.cosTh()),
    m_sinph (cell.sinPhi()),
    m_cosph (cell.cosPhi())
{
  computeVertexCorr (vx);
}
 

void
CaloVertexedCell::computeVertexCorr (const Amg::Vector3D& vx)
{
  const CaloDetDescrElement& dd = *m_cell->caloDDE();
  m_eta += ((vx.x()*dd.cosPhi()+vx.y()*dd.sinPhi())*dd.cosTh()-vx.z()*dd.sinTh())/dd.r();
  m_phi += (vx.x()*dd.sinPhi()-vx.y()*dd.cosPhi())/dd.r() ;
 
  m_sinth += (vx.z()*dd.sinTh() - (vx.x()*dd.cosPhi() + vx.y()*dd.sinPhi())*dd.cosTh()) * dd.cosTh()*dd.sinTh()/dd.r();
  m_costh += - (vx.z()*dd.sinTh() - (vx.x()*dd.cosPhi() + vx.y()*dd.sinPhi())*dd.cosTh()) * dd.sinTh()*dd.sinTh()/dd.r();
 
  m_sinph += - ( vx.y()*dd.cosPhi() - vx.x()*dd.sinPhi())*dd.cosPhi()/dd.r();
  m_cosph += + ( vx.y()*dd.cosPhi() - vx.x()*dd.sinPhi())*dd.sinPhi()/dd.r();
}