PreshowerHitMaker Class Reference

#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

Transform3D

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

Definition at line 15 of file PreshowerHitMaker.h.

XYZPoint

typedef math::XYZVector PreshowerHitMaker::XYZPoint

Definition at line 14 of file PreshowerHitMaker.h.

XYZVector

typedef math::XYZVector PreshowerHitMaker::XYZVector

Definition at line 13 of file PreshowerHitMaker.h.

Constructor & Destructor Documentation

PreshowerHitMaker()

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

inlineoverride

Definition at line 25 of file PreshowerHitMaker.h.

{ ; }

Member Function Documentation

addHit()

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

overridevirtual

Implements CaloHitMaker.

Definition at line 88 of file PreshowerHitMaker.cc.

References anglecorrection1_, anglecorrection2_, funct::cos(), EcalPreshowerGeometry::getClosestCellInPlane(), Calorimeter::getEcalPreshowerGeometry(), CaloSubdetectorGeometry::getGeometry(), CaloHitMaker::hitMap_, LandauFluctuationGenerator::landau(), phase1PixelTopology::layer, layer1valid_, layer2valid_, locToGlobal1_, locToGlobal2_, mip1_, mip2_, CaloHitMaker::moliereRadius, CaloHitMaker::myCalorimeter, phi, point, CaloGeometryHelper::preshowerZPosition(), alignCSCRings::r, random, funct::sin(), digitizers_cfi::strip, theGenerator, totalLayer1_, totalLayer2_, and z.

Referenced by EMShower::compute().

                                                                   {
  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;
}

getHits()

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

inlineoverridevirtual

Implements CaloHitMaker.

Definition at line 30 of file PreshowerHitMaker.h.

References CaloHitMaker::hitMap_.

Referenced by CalorimetryManager::EMShowerSimulation().

{ return hitMap_; };

layer1Calibrated()

float PreshowerHitMaker::layer1Calibrated ( ) const

inline

Definition at line 39 of file PreshowerHitMaker.h.

References totalLayer1_.

{ return 0.024 / 81.1E-6 * totalLayer1_; }

layer2Calibrated()

float PreshowerHitMaker::layer2Calibrated ( ) const

inline

Definition at line 40 of file PreshowerHitMaker.h.

References totalLayer2_.

{ return 0.024 * 0.7 / 81.1E-6 * totalLayer2_; }

setMipEnergy()

void PreshowerHitMaker::setMipEnergy ( double  e1, double  e2  )

inline

Definition at line 32 of file PreshowerHitMaker.h.

References StorageManager_cfg::e1, mip1_, and mip2_.

Referenced by CalorimetryManager::EMShowerSimulation().

                                                 {
    mip1_ = e1;
    mip2_ = e2;
  }

setSpotEnergy()

void PreshowerHitMaker::setSpotEnergy ( double  e )

inlineoverridevirtual

Implements CaloHitMaker.

Definition at line 27 of file PreshowerHitMaker.h.

References MillePedeFileConverter_cfg::e, and CaloHitMaker::spotEnergy.

Referenced by EMShower::compute().

{ spotEnergy = e; }

totalCalibrated()

float PreshowerHitMaker::totalCalibrated ( ) const

inline

Definition at line 41 of file PreshowerHitMaker.h.

References totalLayer1_, and totalLayer2_.

{ return 0.024 / 81.1E-6 * (totalLayer1_ + 0.7 * totalLayer2_); }

totalLayer1()

float PreshowerHitMaker::totalLayer1 ( ) const

inline

Definition at line 37 of file PreshowerHitMaker.h.

References totalLayer1_.

{ return totalLayer1_; }

totalLayer2()

float PreshowerHitMaker::totalLayer2 ( ) const

inline

Definition at line 38 of file PreshowerHitMaker.h.

References totalLayer2_.

{ return totalLayer2_; }

Member Data Documentation

anglecorrection1_

float PreshowerHitMaker::anglecorrection1_

private

Definition at line 52 of file PreshowerHitMaker.h.

Referenced by addHit(), and PreshowerHitMaker().

anglecorrection2_

float PreshowerHitMaker::anglecorrection2_

private

Definition at line 53 of file PreshowerHitMaker.h.

Referenced by addHit(), and PreshowerHitMaker().

layer1valid_

bool PreshowerHitMaker::layer1valid_

private

Definition at line 48 of file PreshowerHitMaker.h.

Referenced by addHit(), and PreshowerHitMaker().

layer2valid_

bool PreshowerHitMaker::layer2valid_

private

Definition at line 49 of file PreshowerHitMaker.h.

Referenced by addHit(), and PreshowerHitMaker().

locToGlobal1_

Transform3D PreshowerHitMaker::locToGlobal1_

private

Definition at line 50 of file PreshowerHitMaker.h.

Referenced by addHit(), and PreshowerHitMaker().

locToGlobal2_

Transform3D PreshowerHitMaker::locToGlobal2_

private

Definition at line 51 of file PreshowerHitMaker.h.

Referenced by addHit(), and PreshowerHitMaker().

mip1_

double PreshowerHitMaker::mip1_

private

Definition at line 54 of file PreshowerHitMaker.h.

Referenced by addHit(), and setMipEnergy().

mip2_

double PreshowerHitMaker::mip2_

private

Definition at line 54 of file PreshowerHitMaker.h.

Referenced by addHit(), and setMipEnergy().

psLayer1Dir_

XYZVector PreshowerHitMaker::psLayer1Dir_

private

Definition at line 45 of file PreshowerHitMaker.h.

psLayer1Entrance_

XYZPoint PreshowerHitMaker::psLayer1Entrance_

private

Definition at line 44 of file PreshowerHitMaker.h.

Referenced by PreshowerHitMaker().

psLayer2Dir_

XYZVector PreshowerHitMaker::psLayer2Dir_

private

Definition at line 47 of file PreshowerHitMaker.h.

psLayer2Entrance_

XYZPoint PreshowerHitMaker::psLayer2Entrance_

private

Definition at line 46 of file PreshowerHitMaker.h.

Referenced by PreshowerHitMaker().

random

const RandomEngineAndDistribution* PreshowerHitMaker::random

private

Definition at line 59 of file PreshowerHitMaker.h.

Referenced by addHit().

theGenerator

const LandauFluctuationGenerator* PreshowerHitMaker::theGenerator

private

The Landau Fluctuation generator.

Definition at line 58 of file PreshowerHitMaker.h.

Referenced by addHit().

totalLayer1_

float PreshowerHitMaker::totalLayer1_

private

Definition at line 55 of file PreshowerHitMaker.h.

Referenced by addHit(), layer1Calibrated(), totalCalibrated(), and totalLayer1().

totalLayer2_

float PreshowerHitMaker::totalLayer2_

private

Definition at line 56 of file PreshowerHitMaker.h.

Referenced by addHit(), layer2Calibrated(), totalCalibrated(), and totalLayer2().