#include <HcalHitMaker.h>

Inheritance diagram for HcalHitMaker:

Public Types
typedef ROOT::Math::Transform3DPJ	Transform3D

typedef math::XYZVector	XYZPoint

typedef math::XYZVector	XYZVector

Public Types inherited from CaloHitMaker
typedef ROOT::Math::Plane3D	Plane3D

typedef math::XYZVector	XYZPoint

typedef math::XYZVector	XYZVector

Public Member Functions
bool	addHit (double r, double phi, unsigned layer=0) override
	add the hit in the HCAL in local coordinates More...

bool	addHit (const XYZPoint &point, unsigned layer=0)
	add the hit in the HCAL in global coordinates More...

const std::map< CaloHitID, float > &	getHits () override

	HcalHitMaker (EcalHitMaker &, unsigned)

bool	setDepth (double, bool inCm=false)
	set the depth in X0 or Lambda0 units depending on showerType More...

void	setSpotEnergy (double e) override
	Set the spot energy. More...

	~HcalHitMaker () override

Public Member Functions inherited from CaloHitMaker
	CaloHitMaker (const CaloGeometryHelper *calo, DetId::Detector det, int subdetn, int cal, unsigned sht=0)

const CaloGeometryHelper *	getCalorimeter () const

virtual	~CaloHitMaker ()

Static Public Member Functions
static int	getSubHcalDet (const FSimTrack *t)

Private Attributes
double	currentDepth_

XYZPoint	ecalEntrance_

Transform3D	locToGlobal_

bool	mapCalculated_

EcalHitMaker &	myGrid

const FSimTrack *	myTrack

int	onHcal

XYZVector	particleDirection

double	radiusFactor_

Additional Inherited Members
Static Protected Member Functions inherited from CaloHitMaker
static XYZPoint	intersect (const Plane3D &p, const XYZPoint &a, const XYZPoint &b, double &t, bool segment, bool debug=false)

Protected Attributes inherited from CaloHitMaker
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

Detailed Description

Definition at line 16 of file HcalHitMaker.h.

Member Typedef Documentation

typedef ROOT::Math::Transform3DPJ HcalHitMaker::Transform3D

Definition at line 20 of file HcalHitMaker.h.

typedef math::XYZVector HcalHitMaker::XYZPoint

Definition at line 19 of file HcalHitMaker.h.

typedef math::XYZVector HcalHitMaker::XYZVector

Definition at line 18 of file HcalHitMaker.h.

Constructor & Destructor Documentation

HcalHitMaker::HcalHitMaker	(	EcalHitMaker &	grid,
		unsigned	shower
	)

Definition at line 10 of file HcalHitMaker.cc.

References funct::abs(), gather_cfg::cout, EcalHitMaker::ecalEntrance(), FSimTrack::ecalEntrance(), ecalEntrance_, CaloHitMaker::EMSHOWER, EcalHitMaker::getFSimTrack(), CaloHitMaker::interactionLength, mapCalculated_, CaloHitMaker::moliereRadius, myGrid, myTrack, particleDirection, radiusFactor_, CoreSimTrack::type(), and RawParticle::Vect().

     : CaloHitMaker(grid.getCalorimeter(),
                    DetId::Hcal,
                    HcalHitMaker::getSubHcalDet(grid.getFSimTrack()),
                    grid.getFSimTrack()->onHcal() ? grid.getFSimTrack()->onHcal() : grid.getFSimTrack()->onVFcal() + 1,
                    shower),
       myGrid(grid),
       myTrack((grid.getFSimTrack())) {
   // normalize the direction
   ecalEntrance_ = myGrid.ecalEntrance();
   particleDirection = myTrack->ecalEntrance().Vect().Unit();
   radiusFactor_ = (EMSHOWER) ? moliereRadius : interactionLength;
   mapCalculated_ = false;
   //std::cout << " Famos HCAL " << grid.getTrack()->onHcal() << " " <<  grid.getTrack()->onVFcal() << " " << showerType << std::endl;
   if (EMSHOWER && (abs(grid.getFSimTrack()->type()) != 11 && grid.getFSimTrack()->type() != 22)) {
     std::cout << " FamosHcalHitMaker : Strange. The following shower has EM type" << std::endl
               << *grid.getFSimTrack() << std::endl;
   }
 }

HcalHitMaker::~HcalHitMaker ( )

inlineoverride

Definition at line 23 of file HcalHitMaker.h.

23 { ; }

Member Function Documentation

bool HcalHitMaker::addHit	(	double	r,
		double	phi,
		unsigned	layer = `0`
	)

overridevirtual

add the hit in the HCAL in local coordinates

Implements CaloHitMaker.

Definition at line 30 of file HcalHitMaker.cc.

References funct::cos(), locToGlobal_, point, radiusFactor_, and funct::sin().

Referenced by HFShower::compute(), HDShower::compute(), EMShower::compute(), CalorimetryManager::HDShowerSimulation(), CalorimetryManager::MuonMipSimulation(), and HDRShower::setHit().

                                                               {
   //  std::cout << " FamosHcalHitMaker::addHit - radiusFactor = " << radiusFactor
   //        << std::endl;
 
   XYZPoint point(r * radiusFactor_ * std::cos(phi), r * radiusFactor_ * std::sin(phi), 0.);
 
   // Watch out !!!! (Point) is a real point in the MathCore terminology (not a redefined a XYZPoint which
   // is actually a XYZVector in the MatchCore terminology). Therefore, the Transform3D is correctly applied
   point = locToGlobal_((Point)point);
   return addHit(point, layer);
 }

bool HcalHitMaker::addHit	(	const XYZPoint &	point,
		unsigned	layer = `0`
	)

add the hit in the HCAL in global coordinates

Definition at line 42 of file HcalHitMaker.cc.

References CaloHitMaker::EMSHOWER, CaloGeometryHelper::getClosestCell(), CaloHitMaker::HADSHOWER, HcalForward, CaloHitMaker::hitMap_, phase1PixelTopology::layer, CaloHitMaker::MIP, CaloHitMaker::myCalorimeter, CaloHitMaker::spotEnergy, and mathSSE::sqrt().

                                                                {
   // Temporary nasty hacks to avoid misbehaviour of not-intended-for-that
   //  getClosestCell in case of large (eta beyond HF ...)  and in EM showers
   if (fabs(point.Z()) > 2000 || fabs(point.X()) > 2000 || fabs(point.Y()) > 2000) {
     if (EMSHOWER)
       edm::LogWarning("HcalHitMaker") << "received a hit very far from the detector " << point
                                       << " coming from an electromagnetic shower. - Ignoring it" << std::endl;
     else if (HADSHOWER)
       edm::LogWarning("HcalHitMaker") << "received a hit very far from the detector " << point
                                       << " coming from a hadron shower. - Ignoring it" << std::endl;
     else if (MIP)
       edm::LogWarning("HcalHitMaker") << "received a hit very far from the detector " << point
                                       << " coming from a muon. - Ignoring it" << std::endl;
     return false;
   }
 
   double pointeta = fabs(point.eta());
   if (pointeta > 5.19)
     return false;
 
   //calculate time of flight
   double dist = std::sqrt(point.X() * point.X() + point.Y() * point.Y() + point.Z() * point.Z());
   double tof = dist / 29.98;  //speed of light
 
   DetId thecellID(myCalorimeter->getClosestCell(point, false, false));
 
   HcalDetId myDetId(thecellID);
 
   //   if ( myDetId.subdetId() == HcalForward ) {
   //     std::cout << "HcalHitMaker : " << point.Z() << " " << myDetId.depth()    << std::endl;
   //   }
 
   //   std::cout << "BEFORE" << std::endl;
   //   std::cout << "HcalHitMaker : subdetId : " << myDetId.subdetId() << std::endl;
   //   std::cout << "HcalHitMaker : depth    : " << myDetId.depth()    << std::endl;
   //   std::cout << "HcalHitMaker : ieta     : " << myDetId.ieta()     << std::endl;
   //   std::cout << "HcalHitMaker : iphi     : " << myDetId.iphi()     << std::endl;
   //   std::cout << "HcalHitMaker : spotE    : " << spotEnergy         << std::endl;
   //   std::cout << "HcalHitMaker : point.X  : " << point.X()          << std::endl;
   //   std::cout << "HcalHitMaker : point.Y  : " << point.Y()          << std::endl;
   //   std::cout << "HcalHitMaker : point.Z  : " << point.Z()          << std::endl;
 
   if (myDetId.subdetId() == HcalForward) {
     int mylayer = layer;
     if (myDetId.depth() == 2) {
       mylayer = (int)layer;
     } else {
       mylayer = 1;
     }
     HcalDetId myDetId2((HcalSubdetector)myDetId.subdetId(), myDetId.ieta(), myDetId.iphi(), mylayer);
     thecellID = myDetId2;
     myDetId = myDetId2;
   }
 
   if (!thecellID.null() && myDetId.depth() > 0) {
     CaloHitID current_id(thecellID.rawId(), tof, 0);  //no track yet
 
     //      std::cout << " FamosHcalHitMaker::addHit - the cell num " << cell
     //              << std::endl;
 
     std::map<CaloHitID, float>::iterator cellitr;
     cellitr = hitMap_.find(current_id);
     if (cellitr == hitMap_.end()) {
       hitMap_.insert(std::pair<CaloHitID, float>(current_id, spotEnergy));
     } else {
       cellitr->second += spotEnergy;
     }
     return true;
   }
   return false;
 }

const std::map<CaloHitID, float>& HcalHitMaker::getHits ( )

inlineoverridevirtual

Implements CaloHitMaker.

Definition at line 35 of file HcalHitMaker.h.

References CaloHitMaker::hitMap_.

Referenced by CalorimetryManager::EMShowerSimulation(), CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::MuonMipSimulation().

35 { return hitMap_; };

CaloHitMaker::hitMap_

std::map< CaloHitID, float > hitMap_

Definition: CaloHitMaker.h:61

static int HcalHitMaker::getSubHcalDet ( const FSimTrack * t )

inlinestatic

Definition at line 54 of file HcalHitMaker.h.

References HcalBarrel, HcalEndcap, HcalForward, FSimTrack::onHcal(), and FSimTrack::onVFcal().

                                                {
     //  std::cout << " getSubHcalDet " << std::endl;
     // According to  DataFormats/ HcalDetId/ interface/ HcalSubdetector.h
     //  std::cout << " onHcal " << t->onHcal() << " onVFcal " << t->onVFcal() << std::endl;
     if (t->onHcal() == 1)
       return HcalBarrel;
     if (t->onHcal() == 2)
       return HcalEndcap;
     if (t->onVFcal() == 2)
       return HcalForward;
     return -1;
   }

bool HcalHitMaker::setDepth	(	double	depth,
		bool	inCm = `false`
	)

set the depth in X0 or Lambda0 units depending on showerType

Definition at line 114 of file HcalHitMaker.cc.

References gather_cfg::cout, currentDepth_, HLT_FULL_cff::depth, CaloHitMaker::EMSHOWER, EcalHitMaker::getFSimTrack(), EcalHitMaker::getSegments(), CaloHitMaker::HADSHOWER, locToGlobal_, myGrid, and particleDirection.

Referenced by HFShower::compute(), HDShower::compute(), EMShower::compute(), CalorimetryManager::HDShowerSimulation(), and HDRShower::setHit().

                                                    {
   currentDepth_ = depth;
   std::vector<CaloSegment>::const_iterator segiterator;
   if (inCm) {
     segiterator =
         find_if(myGrid.getSegments().begin(), myGrid.getSegments().end(), CaloSegment::inSegment(currentDepth_));
   } else {
     if (EMSHOWER)
       segiterator =
           find_if(myGrid.getSegments().begin(), myGrid.getSegments().end(), CaloSegment::inX0Segment(currentDepth_));
 
     //Hadron shower
     if (HADSHOWER)
       segiterator =
           find_if(myGrid.getSegments().begin(), myGrid.getSegments().end(), CaloSegment::inL0Segment(currentDepth_));
   }
 
   if (segiterator == myGrid.getSegments().end()) {
     // Special trick  - As advised by Salavat, no leakage should be simulated
     if (depth > myGrid.getSegments().back().sL0Exit()) {
       segiterator = find_if(myGrid.getSegments().begin(),
                             myGrid.getSegments().end(),
                             CaloSegment::inL0Segment(myGrid.getSegments().back().sL0Exit() - 1.));
       depth = segiterator->sL0Exit() - 1.;
       currentDepth_ = depth;
       if (segiterator == myGrid.getSegments().end()) {
         std::cout << " Could not go at such depth " << EMSHOWER << "  " << currentDepth_ << std::endl;
         std::cout << " Track " << *(myGrid.getFSimTrack()) << std::endl;
         return false;
       }
     } else {
       std::cout << " Could not go at such depth " << EMSHOWER << "  " << currentDepth_ << " "
                 << myGrid.getSegments().back().sL0Exit() << std::endl;
       std::cout << " Track " << *(myGrid.getFSimTrack()) << std::endl;
       return false;
     }
   }
 
   XYZPoint origin;
   if (inCm) {
     origin = segiterator->positionAtDepthincm(currentDepth_);
   } else {
     if (EMSHOWER)
       origin = segiterator->positionAtDepthinX0(currentDepth_);
     if (HADSHOWER)
       origin = segiterator->positionAtDepthinL0(currentDepth_);
   }
   XYZVector zaxis(0, 0, 1);
   XYZVector planeVec1 = (zaxis.Cross(particleDirection)).Unit();
   locToGlobal_ = Transform3D(Point(0, 0, 0),
                              Point(0, 0, 1),
                              Point(1, 0, 0),
                              (Point)origin,
                              (Point)(origin + particleDirection),
                              (Point)(origin + planeVec1));
   return true;
 }