#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

Definition at line 27 of file CaloHitMaker.h.

typedef math::XYZVector CaloHitMaker::XYZPoint

Definition at line 26 of file CaloHitMaker.h.

typedef math::XYZVector CaloHitMaker::XYZVector

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 ( )

inlinevirtual

Definition at line 30 of file CaloHitMaker.h.

30 {;}

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`
	)

staticprotected

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.);
 }