#include <PreshowerHitMaker.h>

Inheritance diagram for PreshowerHitMaker:

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

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

float	layer1Calibrated () const

float	layer2Calibrated () const

	PreshowerHitMaker (CaloGeometryHelper calo, const XYZPoint &, const XYZVector &, const XYZPoint &, const XYZVector &, const LandauFluctuationGenerator aGenerator, const RandomEngineAndDistribution *engine)

void	setMipEnergy (double e1, double e2)

void	setSpotEnergy (double e) override

float	totalCalibrated () const

float	totalLayer1 () const

float	totalLayer2 () const

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

Private Attributes
float	anglecorrection1_

float	anglecorrection2_

bool	layer1valid_

bool	layer2valid_

Transform3D	locToGlobal1_

Transform3D	locToGlobal2_

double	mip1_

double	mip2_

XYZVector	psLayer1Dir_

XYZPoint	psLayer1Entrance_

XYZVector	psLayer2Dir_

XYZPoint	psLayer2Entrance_

const RandomEngineAndDistribution *	random

const LandauFluctuationGenerator *	theGenerator
	The Landau Fluctuation generator. More...

float	totalLayer1_

float	totalLayer2_

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 11 of file PreshowerHitMaker.h.

Member Typedef Documentation

typedef ROOT::Math::Transform3DPJ PreshowerHitMaker::Transform3D

Definition at line 17 of file PreshowerHitMaker.h.

typedef math::XYZVector PreshowerHitMaker::XYZPoint

Definition at line 16 of file PreshowerHitMaker.h.

typedef math::XYZVector PreshowerHitMaker::XYZVector

Definition at line 15 of file PreshowerHitMaker.h.

Constructor & Destructor Documentation

PreshowerHitMaker::PreshowerHitMaker	(	CaloGeometryHelper *	calo,
		const XYZPoint &	layer1entrance,
		const XYZVector &	layer1dir,
		const XYZPoint &	layer2entrance,
		const XYZVector &	layer2dir,
		const LandauFluctuationGenerator *	aGenerator,
		const RandomEngineAndDistribution *	engine
	)

Definition at line 18 of file PreshowerHitMaker.cc.

References anglecorrection1_, anglecorrection2_, CaloHitMaker::intersect(), layer1valid_, layer2valid_, locToGlobal1_, locToGlobal2_, CaloHitMaker::myCalorimeter, CaloGeometryHelper::preshowerZPosition(), psLayer1Entrance_, psLayer2Entrance_, z, and trackingPlots::zaxis.

                                               :
   CaloHitMaker(calo,DetId::Ecal,EcalPreshower,2),
   psLayer1Entrance_(layer1entrance),
   psLayer1Dir_(layer1dir),
   psLayer2Entrance_(layer2entrance),
   psLayer2Dir_(layer2dir),
   totalLayer1_(0.),totalLayer2_(0.),
   theGenerator(aGenerator),
   random(engine)
 {
   double dummyt;
   anglecorrection1_ = 0.;
   anglecorrection2_ = 0.;
    // Check if the entrance points are really on the wafers
   // Layer 1 
   layer1valid_ = (layer1entrance.Mag2()>0.);
   if(layer1valid_)
     {
       int z=(psLayer1Entrance_.z()>0)? 1:-1;
       Plane3D plan1(0.,0.,1.,-z*myCalorimeter->preshowerZPosition(1));
 
       psLayer1Entrance_ = intersect(plan1,layer1entrance,layer1entrance+layer1dir,dummyt,false);
 
       XYZVector zaxis(0,0,1);
       XYZVector planeVec1=(zaxis.Cross(layer1dir)).Unit();
       locToGlobal1_=Transform3D(Point(0,0,0),
                 Point(0,0,1),
                 Point(1,0,0),
                 (Point)psLayer1Entrance_,
                 (Point)(psLayer1Entrance_+layer1dir),
                 (Point)(psLayer1Entrance_+planeVec1));
 
       anglecorrection1_ = fabs(zaxis.Dot(layer1dir));
       if(anglecorrection1_!=0.) anglecorrection1_ = 1./anglecorrection1_;
       //      std::cout << " Layer 1 entrance " << psLayer1Entrance_ << std::endl;
       //      std::cout << " Layer 1 corr " << anglecorrection1_ << std::endl;
     }
 
   // Layer 2
   layer2valid_ = (layer2entrance.Mag2()>0.);
   if(layer2valid_)
     {
       int z=(psLayer2Entrance_.z()>0) ? 1:-1;
       Plane3D plan2(0.,0.,1.,-z*myCalorimeter->preshowerZPosition(2));
       
       psLayer2Entrance_ = intersect(plan2,layer2entrance,layer2entrance+layer2dir,dummyt,false);
 
       XYZVector zaxis(0,0,1);
       XYZVector planeVec2=(zaxis.Cross(layer2dir)).Unit();
       locToGlobal2_=Transform3D(Point(0,0,0),
                 Point(0,0,1),
                 Point(1,0,0),
                    (Point)psLayer2Entrance_,
                    (Point)(psLayer2Entrance_+layer2dir),
                    (Point)(psLayer2Entrance_+planeVec2));
       
       anglecorrection2_ = fabs(zaxis.Dot(layer2dir));
       if(anglecorrection2_!=0.) anglecorrection2_ = 1./anglecorrection2_;
       //      std::cout << " Layer 2 entrance " << psLayer2Entrance_ << std::endl;
       //      std::cout << " Layer 2 corr " << anglecorrection2_ << std::endl;
     }
   //  theGenerator=LandauFluctuationGenerator();
 }

PreshowerHitMaker::~PreshowerHitMaker ( )

inlineoverride

Definition at line 27 of file PreshowerHitMaker.h.

27 {;}

Member Function Documentation

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

overridevirtual

Implements CaloHitMaker.

Definition at line 91 of file PreshowerHitMaker.cc.

References anglecorrection1_, anglecorrection2_, funct::cos(), EcalPreshowerGeometry::getClosestCellInPlane(), Calorimeter::getEcalPreshowerGeometry(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloHitMaker::hitMap_, LandauFluctuationGenerator::landau(), layer1valid_, layer2valid_, locToGlobal1_, locToGlobal2_, PV3DBase< T, PVType, FrameType >::mag(), mip1_, mip2_, CaloHitMaker::moliereRadius, CaloHitMaker::myCalorimeter, point, CaloGeometryHelper::preshowerZPosition(), random, funct::sin(), digi_MixPreMix_cfi::strip, theGenerator, totalLayer1_, totalLayer2_, and z.

Referenced by EMShower::compute(), and setSpotEnergy().

 {
   if((layer==1&&!layer1valid_)||((layer==2&&!layer2valid_))) return false;
 
   r*=moliereRadius;
   XYZPoint point (r*std::cos(phi),r*std::sin(phi),0.);
   point =  (layer==1) ? locToGlobal1_((Point)point) : locToGlobal2_((Point)point);
   //  std::cout << "  Point " << point  << std::endl;
   int z=(point.z()>0) ? 1: -1;
   point = XYZPoint(point.x(),point.y(),z*myCalorimeter->preshowerZPosition(layer));
   //  std::cout << "r " << r << "  Point after " << point  << std::endl;
   //  std::cout << " Layer " << layer << " " << point << std::endl;
   DetId strip = myCalorimeter->getEcalPreshowerGeometry()->getClosestCellInPlane(GlobalPoint(point.x(),point.y(),point.z()),layer);
 
   float meanspot=(layer==1) ? mip1_ : mip2_; 
   float spote = meanspot + 0.000021*theGenerator->landau(random);
   spote *= ( (layer==1) ? anglecorrection1_ : anglecorrection2_ );
 
   if(!strip.null())
     {
       //calculate time of flight
       double tof = (myCalorimeter->getEcalPreshowerGeometry()->getGeometry(strip)->getPosition().mag())/29.98;//speed of light
       CaloHitID current_id(strip.rawId(),tof,0); //no track yet
       std::map<CaloHitID,float>::iterator cellitr;
       cellitr = hitMap_.find(current_id);
       if( cellitr==hitMap_.end())
     {
       hitMap_.insert(std::pair<CaloHitID,float>(current_id,spote));
     }
       else
     {
       cellitr->second+=spote;
     }  
       //      std::cout << " found " << stripNumber << " " << spote <<std::endl;
       if(layer==1){
     totalLayer1_+=spote;
       }
       else if (layer==2) {
     totalLayer2_+=spote;
       }
       return true;
     }
   //  std::cout << "  Could not find a cell " << point << std::endl;
   return false;
 }