#include <CaloHitMaker.h>

Inheritance diagram for CaloHitMaker:

Public Types
typedef ROOT::Math::Plane3D	Plane3D
typedef math::XYZVector	XYZPoint
typedef math::XYZVector	XYZVector
Public Member Functions
virtual bool	addHit (double r, double phi, unsigned layer=0)=0
	CaloHitMaker (const CaloGeometryHelper *calo, DetId::Detector det, int subdetn, int cal, unsigned sht=0)
const CaloGeometryHelper *	getCalorimeter () const
virtual const std::map < CaloHitID, float > &	getHits ()=0
virtual void	setSpotEnergy (double e)=0
virtual	~CaloHitMaker ()
Static Protected Member Functions
static XYZPoint	intersect (const Plane3D &p, const XYZPoint &a, const XYZPoint &b, double &t, bool segment, bool debug=false)
Protected Attributes
bool	EMSHOWER
bool	HADSHOWER
std::map< CaloHitID, float >	hitMap_
double	interactionLength
bool	MIP
double	moliereRadius
const CaloGeometryHelper *	myCalorimeter
unsigned	showerType_
double	spotEnergy
const CalorimeterProperties *	theCaloProperties
Private Attributes
DetId::Detector	base_
int	onCal_
int	subdetn_

Detailed Description

Definition at line 21 of file CaloHitMaker.h.

Member Typedef Documentation

typedef ROOT::Math::Plane3D CaloHitMaker::Plane3D

Reimplemented in EcalHitMaker.

Definition at line 27 of file CaloHitMaker.h.

typedef math::XYZVector CaloHitMaker::XYZPoint

Reimplemented in EcalHitMaker, HcalHitMaker, and PreshowerHitMaker.

Definition at line 26 of file CaloHitMaker.h.

typedef math::XYZVector CaloHitMaker::XYZVector

Reimplemented in EcalHitMaker, HcalHitMaker, and PreshowerHitMaker.

Definition at line 25 of file CaloHitMaker.h.

Constructor & Destructor Documentation

CaloHitMaker::CaloHitMaker	(	const CaloGeometryHelper *	calo,
		DetId::Detector	det,
		int	subdetn,
		int	cal,
		unsigned	sht = `0`
	)

Definition at line 10 of file CaloHitMaker.cc.

References base_, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, Calorimeter::ecalProperties(), EMSHOWER, HADSHOWER, DetId::Hcal, Calorimeter::hcalProperties(), CalorimeterProperties::interactionLength(), interactionLength, Calorimeter::layer1Properties(), MIP, CalorimeterProperties::moliereRadius(), moliereRadius, myCalorimeter, onCal_, subdetn_, and theCaloProperties.

  :myCalorimeter(theCalo),theCaloProperties(NULL),base_(basedet),subdetn_(subdetn),onCal_(cal),showerType_(sht)
{
  //  std::cout << " FamosCalorimeter " << basedet << " " << cal << std::endl;
  EMSHOWER=(sht==0);
  HADSHOWER=(sht==1);
  MIP=(sht==2);
  if(base_==DetId::Ecal&&(subdetn_==EcalBarrel||subdetn==EcalEndcap)&&onCal_)
    theCaloProperties = (CalorimeterProperties*)myCalorimeter->ecalProperties(onCal_);
  // is it really necessary to cast here ? 
  if(base_==DetId::Ecal&&subdetn_==EcalPreshower&&onCal_)
    theCaloProperties = (PreshowerProperties*)myCalorimeter->layer1Properties(onCal_);
  if(base_==DetId::Hcal&&cal) theCaloProperties = myCalorimeter->hcalProperties(onCal_);

  if(theCaloProperties)
    {
      moliereRadius=theCaloProperties->moliereRadius();
      interactionLength=theCaloProperties->interactionLength();
    }
  else
    {
      moliereRadius=999;
      interactionLength=999;
    }
}

virtual CaloHitMaker::~CaloHitMaker ( ) [inline, virtual]

Definition at line 30 of file CaloHitMaker.h.

{;}

Member Function Documentation

virtual bool CaloHitMaker::addHit	(	double	r,
		double	phi,
		unsigned	layer = `0`
	)		`[pure virtual]`

Implemented in EcalHitMaker, HcalHitMaker, and PreshowerHitMaker.

const CaloGeometryHelper* CaloHitMaker::getCalorimeter ( ) const [inline]

Definition at line 36 of file CaloHitMaker.h.

References myCalorimeter.

    {
      //      std::cout << "CaloHitMaker is returning myCalorimeter " << myCalorimeter << std::endl;
      return myCalorimeter;
    }

virtual const std::map<CaloHitID,float>& CaloHitMaker::getHits ( ) [pure virtual]

Implemented in EcalHitMaker, HcalHitMaker, and PreshowerHitMaker.

CaloHitMaker::XYZPoint CaloHitMaker::intersect	(	const Plane3D &	p,
		const XYZPoint &	a,
		const XYZPoint &	b,
		double &	t,
		bool	segment,
		bool	debug = `false`
	)		`[static, protected]`

computes the intersection between a straight line defined by a & b and a plan

Definition at line 38 of file CaloHitMaker.cc.

References gather_cfg::cout.

Referenced by EcalHitMaker::ecalCellLine(), EcalHitMaker::getPads(), and PreshowerHitMaker::PreshowerHitMaker().

{
  t=-9999.;
  // En Attendant //
  XYZVector normal = p.Normal();
  double AAA = normal.X();
  double BBB = normal.Y();
  double CCC = normal.Z();
  //  double DDD = p.Distance(Point(0.,0.,0.));
  double DDD = p.HesseDistance();
  //  double denom = p.A()*(b.X()-a.X()) + p.B()*(b.Y()-a.Y()) + p.C()*(b.Z()-a.Z());
  double denom = AAA*(b.X()-a.X()) + BBB*(b.Y()-a.Y()) + CCC*(b.Z()-a.Z());
  if(denom!=0.)
    {
      // t=-(p.A()*a.X()+p.B()*a.Y()+p.C()*a.Z()+p.D());
      t=-(AAA*a.X()+BBB*a.Y()+CCC*a.Z()+DDD);
      t/=denom;
      if(debug) std::cout << " T = " << t <<std::endl; 
      if(segment)
        {
          if(t>=0&&t<=1)
            return XYZPoint(a.X()+(b.X()-a.X())*t,
                            a.Y()+(b.Y()-a.Y())*t,
                            a.Z()+(b.Z()-a.Z())*t);      
        }
      else
        {
          return XYZPoint(a.X()+(b.X()-a.X())*t,
                          a.Y()+(b.Y()-a.Y())*t,
                          a.Z()+(b.Z()-a.Z())*t);      
        }
          
     
    }

  return XYZPoint(0.,0.,0.);
}