#include <EcalHitMaker.h>

Inheritance diagram for EcalHitMaker:

Public Types
typedef ROOT::Math::Plane3D	Plane3D

typedef math::XYZVector	XYZNormal

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

bool	addHitDepth (double r, double phi, double depth=-1)

const XYZPoint &	ecalEntrance () const
	used in FamosHcalHitMaker More...

double	ecalHcalGapTotalL0 () const
	ECAL-HCAL transition. More...

double	ecalHcalGapTotalX0 () const
	ECAL-HCAL transition. More...

	EcalHitMaker (CaloGeometryHelper calo, const XYZPoint &ecalentrance, const DetId &cell, int onEcal, unsigned size, unsigned showertype, const RandomEngineAndDistribution engine)

double	ecalTotalL0 () const
	in the ECAL More...

double	ecalTotalX0 () const
	in the ECAL More...

unsigned	fastInsideCell (const CLHEP::Hep2Vector &point, double &sp, bool debug=false)

const std::vector< Crystal > &	getCrystals () const
	for debugging More...

const FSimTrack *	getFSimTrack () const
	To retrieve the track. More...

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

bool	getPads (double depth, bool inCm=false)

const std::vector< CaloSegment > &	getSegments () const

double	getX0back () const

double	hcalTotalL0 () const
	in the HCAL More...

double	hcalTotalX0 () const
	in the HCAL More...

double	ps1TotalL0 () const
	total number of L0 in the PS (Layer1). More...

double	ps1TotalX0 () const

double	ps2eeTotalL0 () const

double	ps2eeTotalX0 () const

double	ps2TotalL0 () const
	total number of L0 in the PS (Layer2). More...

double	ps2TotalX0 () const
	total number of X0 in the PS (Layer2). More...

void	setCrackPadSurvivalProbability (double val)

void	setPreshowerPresent (bool ps)

void	setPulledPadSurvivalProbability (double val)

void	setRadiusFactor (double r)

void	setSpotEnergy (double e) override

void	setTrackParameters (const XYZNormal &normal, double X0depthoffset, const FSimTrack &theTrack)

double	totalL0 () const
	Number of interaction length "seen" by the track. More...

double	totalX0 () const

double	x0DepthOffset () const
	get the offset (e.g the number of X0 after which the shower starts) More...

	~EcalHitMaker () 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 Types
typedef std::pair< CaloDirection, unsigned >	neighbour

Private Member Functions
void	buildGeometry ()

void	buildSegments (const std::vector< CaloPoint > &cp)

void	cellLine (std::vector< CaloPoint > &cp)

void	configureGeometry ()

void	convertIntegerCoordinates (double x, double y, unsigned &ix, unsigned &iy) const

CLHEP::Hep2Vector &	correspondingEdge (neighbour &myneighbour, CaloDirection dir2)

void	cracksPads (std::vector< neighbour > &cracks, unsigned iq)

bool	diagonalEdge (unsigned myPad, CaloDirection dir, CLHEP::Hep2Vector &point)

void	ecalCellLine (const XYZPoint &, const XYZPoint &, std::vector< CaloPoint > &cp)

void	gapsLifting (std::vector< neighbour > &gaps, unsigned iq)

void	hcalCellLine (std::vector< CaloPoint > &cp) const

bool	inside3D (const std::vector< XYZPoint > &, const XYZPoint &p) const

void	prepareCrystalNumberArray ()

void	preshowerCellLine (std::vector< CaloPoint > &cp) const

bool	pulled (const XYZPoint &origin, const XYZNormal &normal, XYZPoint &fPoint) const

void	reorganizePads ()

bool	unbalancedDirection (const std::vector< neighbour > &dirs, unsigned &unb, unsigned &dir1, unsigned &dir2)

Private Attributes
double	bfactor_

std::vector< DetId >	CellsWindow_

int	central_

bool	configuredGeometry_

std::vector< XYZPoint >	corners

double	crackPadProbability_

std::vector< CrystalPad >	crackpadsatdepth_

double	currentdepth_

bool	detailedShowerTail_

std::map< DetId, unsigned >	DetIdMap_

bool	doreorg_

XYZPoint	EcalEntrance_

int	ecalFirstSegment_

unsigned	etasize_

bool	hitmaphasbeencalculated_

std::vector< float >	hits_

std::vector< CaloPoint >	intersections_

double	L0ECAL_

double	L0EHGAP_

double	L0HCAL_

double	L0PS1_

double	L0PS2_

double	L0PS2EE_

double	maxX0_

std::vector< CLHEP::Hep2Vector >	mycorners

std::vector< std::vector< unsigned > >	myCrystalNumberArray_

CrystalWindowMap *	myCrystalWindowMap_

const FSimTrack *	myTrack_

unsigned	ncrackpadsatdepth_

unsigned	ncrystals_

XYZNormal	normal_

unsigned	npadsatdepth_

unsigned	nx_

unsigned	ny_

int	onEcal_

double	outsideWindowEnergy_

std::vector< CrystalPad >	padsatdepth_

unsigned	phisize_

Crystal	pivot_

Plane3D	plan_

double	pulledPadProbability_

double	radiusCorrectionFactor_

double	radiusFactor_

const RandomEngineAndDistribution *	random

double	rearleakage_

std::vector< Crystal >	regionOfInterest_

std::vector< CaloSegment >	segments_

bool	simulatePreshower_

double	sizex_

double	sizey_

double	totalL0_

double	totalX0_

bool	truncatedGrid_

std::vector< bool >	validPads_

double	X0depthoffset_

double	X0ECAL_

double	X0EHGAP_

double	X0HCAL_

double	X0PS1_

double	X0PS2_

double	X0PS2EE_

double	xmax_

double	xmin_

double	ymax_

double	ymin_

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 24 of file EcalHitMaker.h.

Member Typedef Documentation

◆ neighbour

typedef std::pair<CaloDirection, unsigned> EcalHitMaker::neighbour

private

Definition at line 168 of file EcalHitMaker.h.

◆ Plane3D

typedef ROOT::Math::Plane3D EcalHitMaker::Plane3D

Definition at line 29 of file EcalHitMaker.h.

◆ XYZNormal

typedef math::XYZVector EcalHitMaker::XYZNormal

Definition at line 28 of file EcalHitMaker.h.

◆ XYZPoint

typedef math::XYZVector EcalHitMaker::XYZPoint

Definition at line 27 of file EcalHitMaker.h.

◆ XYZVector

typedef math::XYZVector EcalHitMaker::XYZVector

Definition at line 26 of file EcalHitMaker.h.

Constructor & Destructor Documentation

◆ EcalHitMaker()

EcalHitMaker::EcalHitMaker	(	CaloGeometryHelper *	calo,
		const XYZPoint &	ecalentrance,
		const DetId &	cell,
		int	onEcal,
		unsigned	size,
		unsigned	showertype,
		const RandomEngineAndDistribution *	engine
	)

Definition at line 33 of file EcalHitMaker.cc.

References bfactor_, CaloGeometryHelper::buildCrystal(), buildGeometry(), CellsWindow_, central_, corners, gather_cfg::cout, doreorg_, EcalEntrance_, ecalFirstSegment_, etasize_, CaloGeometryHelper::getClosestCell(), Crystal::getDetId(), Calorimeter::getEcalEndcapGeometry(), CaloSubdetectorGeometry::getGeometry(), CaloGeometryHelper::getWindow(), hitmaphasbeencalculated_, Histos::instance(), L0ECAL_, L0EHGAP_, L0HCAL_, L0PS1_, L0PS2_, L0PS2EE_, maxX0_, CaloHitMaker::myCalorimeter, mycorners, myCrystalWindowMap_, ncrystals_, outsideWindowEnergy_, phisize_, pivot_, pulledPadProbability_, rearleakage_, simulatePreshower_, findQualityFiles::size, totalL0_, totalX0_, truncatedGrid_, X0depthoffset_, X0ECAL_, X0EHGAP_, X0HCAL_, X0PS1_, X0PS2_, and X0PS2EE_.

     : CaloHitMaker(theCalo, DetId::Ecal, ((onEcal == 1) ? EcalBarrel : EcalEndcap), onEcal, showertype),
       EcalEntrance_(ecalentrance),
       onEcal_(onEcal),
       myTrack_(nullptr),
       random(engine) {
 #ifdef FAMOSDEBUG
   myHistos = Histos::instance();
 #endif
   //  myHistos->debug("Constructeur EcalHitMaker");
   simulatePreshower_ = true;
   X0depthoffset_ = 0.;
   X0PS1_ = 0.;
   X0PS2_ = 0.;
   X0PS2EE_ = 0.;
   X0ECAL_ = 0.;
   X0EHGAP_ = 0.;
   X0HCAL_ = 0.;
   L0PS1_ = 0.;
   L0PS2_ = 0.;
   L0PS2EE_ = 0.;
   L0ECAL_ = 0.;
   L0EHGAP_ = 0.;
   L0HCAL_ = 0.;
   maxX0_ = 0.;
   totalX0_ = 0;
   totalL0_ = 0.;
   pulledPadProbability_ = 1.;
   outsideWindowEnergy_ = 0.;
   rearleakage_ = 0.;
   bfactor_ = 1.;
   ncrystals_ = 0;
 
   doreorg_ = !showertype;
 
   hitmaphasbeencalculated_ = false;
 
   if (onEcal)
     myCalorimeter->buildCrystal(cell, pivot_);
   else
     pivot_ = Crystal();
   central_ = onEcal == 1;
   ecalFirstSegment_ = -1;
 
   myCrystalWindowMap_ = nullptr;
   // In some cases, a "dummy" grid, not based on a cell, can be built. The previous variables
   // should however be initialized. In such a case onEcal=0
   if (!onEcal)
     return;
 
   // Same size in eta-phi
   etasize_ = size;
   phisize_ = size;
 
   // Build the grid
   // The result is put in CellsWindow and is ordered by distance to the pivot
   myCalorimeter->getWindow(pivot_.getDetId(), size, size, CellsWindow_);
 
   buildGeometry();
   //  std::cout << " Geometry built " << regionOfInterest_.size() << std::endl;
 
   truncatedGrid_ = CellsWindow_.size() != (etasize_ * phisize_);
 
   // A local vector of corners
   mycorners.resize(4);
   corners.resize(4);
 
 #ifdef DEBUGGW
   myHistos->fill("h10", EcalEntrance_.eta(), CellsWindow_.size());
   if (onEcal == 2) {
     myHistos->fill("h20", EcalEntrance_.perp(), CellsWindow_.size());
     if (EcalEntrance_.perp() > 70 && EcalEntrance_.perp() < 80 && CellsWindow_.size() < 35) {
       std::cout << " Truncated grid " << CellsWindow_.size() << " " << EcalEntrance_.perp() << std::endl;
       std::cout << " Pivot "
                 << myCalorimeter->getEcalEndcapGeometry()->getGeometry(pivot_.getDetId())->getPosition().perp();
       std::cout << EEDetId(pivot_.getDetId()) << std::endl;
 
       std::cout << " Test getClosestCell " << EcalEntrance_ << std::endl;
       DetId testcell = myCalorimeter->getClosestCell(EcalEntrance_, true, false);
       std::cout << " Result " << EEDetId(testcell) << std::endl;
       std::cout << " Position " << myCalorimeter->getEcalEndcapGeometry()->getGeometry(testcell)->getPosition()
                 << std::endl;
     }
   }
 
 #endif
 }

◆ ~EcalHitMaker()

EcalHitMaker::~EcalHitMaker ( )

override

Definition at line 127 of file EcalHitMaker.cc.

References myCrystalWindowMap_.

                             {
   if (myCrystalWindowMap_ != nullptr) {
     delete myCrystalWindowMap_;
   }
 }

Member Function Documentation

◆ addHit()

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

overridevirtual

Implements CaloHitMaker.

Definition at line 184 of file EcalHitMaker.cc.

References funct::cos(), fastInsideCell(), RandomEngineAndDistribution::flatShoot(), hits_, outsideWindowEnergy_, phi, point, alignCSCRings::r, radiusFactor_, random, funct::sin(), and CaloHitMaker::spotEnergy.

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

                                                               {
   //  std::cout <<" Addhit " << std::endl;
   //  std::cout << " Before insideCell " << std::endl;
   double sp(1.);
   //  std::cout << " Trying to add " << r << " " << phi << " " << radiusFactor_ << std::endl;
   r *= radiusFactor_;
   CLHEP::Hep2Vector point(r * std::cos(phi), r * std::sin(phi));
   //  std::cout << "point " << point << std::endl;
   //  CellID cellid=insideCell(point,sp);
   unsigned xtal = fastInsideCell(point, sp);
   //  if(cellid.isZero()) std::cout << " cell is Zero " << std::endl;
   if (xtal < 1000) {
     if (sp == 1.)
       hits_[xtal] += spotEnergy;
     else
       hits_[xtal] += (random->flatShoot() < sp) * spotEnergy;
     return true;
   }
 
   outsideWindowEnergy_ += spotEnergy;
   //  std::cout << " This hit ; r= " << point << " hasn't been added "<<std::endl;
   //  std::cout << " Xtal " << xtal << std::endl;
   //  for(unsigned ip=0;ip<npadsatdepth_;++ip)
   //    {
   //      std::cout << padsatdepth_[ip] << std::endl;
   //    }
 
   return false;
 }

◆ addHitDepth()

bool EcalHitMaker::addHitDepth	(	double	r,
		double	phi,
		double	depth = `-1`
	)

Definition at line 133 of file EcalHitMaker.cc.

References funct::cos(), hcalRecHitTable_cff::depth, fastInsideCell(), hits_, outsideWindowEnergy_, phi, point, alignCSCRings::r, radiusFactor_, rearleakage_, regionOfInterest_, funct::sin(), CaloHitMaker::spotEnergy, and X0depthoffset_.

Referenced by EMShower::compute().

                                                                  {
   //  std::cout << " Add hit depth called; Current deph is "  << currentdepth_;
   //  std::cout << " Required depth is " << depth << std::endl;
   depth += X0depthoffset_;
   double sp(1.);
   r *= radiusFactor_;
   CLHEP::Hep2Vector point(r * std::cos(phi), r * std::sin(phi));
 
   unsigned xtal = fastInsideCell(point, sp);
   //  if(cellid.isZero()) std::cout << " cell is Zero " << std::endl;
   //  if(xtal<1000)
   //    {
   //      std::cout << "Result " << regionOfInterest_[xtal].getX0Back() << " " ;
   //      std::cout << depth << std::endl;
   //    }
   //  myHistos->fill("h5000",depth);
   if (xtal < 1000) {
     //      myHistos->fill("h5002",regionOfInterest_[xtal].getX0Back(),depth);
     //      myHistos->fill("h5003",ecalentrance_.eta(),maxX0_);
     if (regionOfInterest_[xtal].getX0Back() > depth) {
       hits_[xtal] += spotEnergy;
       //      myHistos->fill("h5005",r);
       return true;
     } else {
       rearleakage_ += spotEnergy;
     }
   } else {
     //      std::cout << " Return false " << std::endl;
     //      std::cout << " Add hit depth called; Current deph is "  << currentdepth_;
     //      std::cout << " Required depth is " << depth << std::endl;
     //      std::cout << " R = " << r << " " << radiusFactor_ << std::endl;
   }
 
   outsideWindowEnergy_ += spotEnergy;
   return false;
 }

◆ buildGeometry()

void EcalHitMaker::buildGeometry ( )

private

Definition at line 790 of file EcalHitMaker.cc.

References CaloGeometryHelper::buildCrystal(), CellsWindow_, configuredGeometry_, DetIdMap_, relativeConstraints::empty, hits_, CaloHitMaker::myCalorimeter, myCrystalNumberArray_, myCrystalWindowMap_, ncrystals_, nx_, ny_, padsatdepth_, phisize_, regionOfInterest_, and validPads_.

Referenced by EcalHitMaker().

                                  {
   configuredGeometry_ = false;
   ncrystals_ = CellsWindow_.size();
   // create the vector with of pads with the appropriate size
   padsatdepth_.resize(ncrystals_);
 
   // This is fully correct in the barrel.
   ny_ = phisize_;
   nx_ = ncrystals_ / ny_;
   std::vector<unsigned> empty;
   empty.resize(ny_, 0);
   myCrystalNumberArray_.reserve((unsigned)nx_);
   for (unsigned inx = 0; inx < (unsigned)nx_; ++inx) {
     myCrystalNumberArray_.push_back(empty);
   }
 
   hits_.resize(ncrystals_, 0.);
   regionOfInterest_.clear();
   regionOfInterest_.resize(ncrystals_);
   validPads_.resize(ncrystals_);
   for (unsigned ic = 0; ic < ncrystals_; ++ic) {
     myCalorimeter->buildCrystal(CellsWindow_[ic], regionOfInterest_[ic]);
     regionOfInterest_[ic].setNumber(ic);
     DetIdMap_.insert(std::pair<DetId, unsigned>(CellsWindow_[ic], ic));
   }
 
   // Computes the map of the neighbours
   myCrystalWindowMap_ = new CrystalWindowMap(myCalorimeter, regionOfInterest_);
 }

◆ buildSegments()

void EcalHitMaker::buildSegments ( const std::vector< CaloPoint > & cp )

private

Definition at line 583 of file EcalHitMaker.cc.

References CaloGeometryHelper::borderCrossing(), hgcalTopologyTester_cfi::cell2, gather_cfg::cout, CaloSegment::CRACK, DetId::Ecal, EcalBarrel, EcalEndcap, EcalEntrance_, ecalFirstSegment_, CaloSegment::ECALHCALGAP, EcalPreshower, CaloSegment::GAP, CaloSegment::HCAL, DetId::Hcal, L0ECAL_, L0EHGAP_, L0HCAL_, CaloSegment::L0length(), L0PS1_, L0PS2_, L0PS2EE_, CaloSegment::length(), CaloHitMaker::myCalorimeter, myTrack_, CaloSegment::PbWO4, FSimVertex::position(), CaloSegment::PS, CaloSegment::PSEEGAP, alignCSCRings::s, segments_, findQualityFiles::size, totalL0_, totalX0_, FSimTrack::vertex(), X0ECAL_, X0EHGAP_, X0HCAL_, CaloSegment::X0length(), X0PS1_, X0PS2_, and X0PS2EE_.

Referenced by setTrackParameters().

                                                                {
   //  myHistos->debug();
   //  TimeMe theT("FamosGrid::buildSegments");
   unsigned size = cp.size();
   if (size % 2 != 0) {
     //      std::cout << " There is a problem " << std::endl;
     return;
   }
   //  myHistos->debug();
   unsigned nsegments = size / 2;
   segments_.reserve(nsegments);
   if (size == 0)
     return;
   // curv abs
   double s = 0.;
   double sX0 = 0.;
   double sL0 = 0.;
 
 #ifdef DEBUGCELLLINE
   unsigned ncrossedxtals = 0;
 #endif
   unsigned is = 0;
   while (is < nsegments) {
     if (cp[2 * is].getDetId() != cp[2 * is + 1].getDetId() && cp[2 * is].whichDetector() != DetId::Hcal &&
         cp[2 * is + 1].whichDetector() != DetId::Hcal) {
       //      std::cout << " Problem with the segments " << std::endl;
       //      std::cout << cp[2*is].whichDetector() << " " << cp[2*is+1].whichDetector() << std::endl;
       //      std::cout << is << " " <<cp[2*is].getDetId().rawId() << std::endl;
       //      std::cout << (2*is+1) << " " <<cp[2*is+1].getDetId().rawId() << std::endl;
       ++is;
       continue;
     }
 
     // Check if it is a Preshower segment - Layer 1
     // One segment per layer, nothing between
     //      myHistos->debug("Just avant Preshower");
     if (cp[2 * is].whichDetector() == DetId::Ecal && cp[2 * is].whichSubDetector() == EcalPreshower &&
         cp[2 * is].whichLayer() == 1) {
       if (cp[2 * is + 1].whichDetector() == DetId::Ecal && cp[2 * is + 1].whichSubDetector() == EcalPreshower &&
           cp[2 * is + 1].whichLayer() == 1) {
         CaloSegment preshsegment(cp[2 * is], cp[2 * is + 1], s, sX0, sL0, CaloSegment::PS, myCalorimeter);
         segments_.push_back(preshsegment);
         //        std::cout << " Added (1-1)" << preshsegment << std::endl;
         s += preshsegment.length();
         sX0 += preshsegment.X0length();
         sL0 += preshsegment.L0length();
         X0PS1_ += preshsegment.X0length();
         L0PS1_ += preshsegment.L0length();
       } else {
         std::cout << " Strange segment between Preshower1 and " << cp[2 * is + 1].whichDetector();
         std::cout << std::endl;
       }
       ++is;
       continue;
     }
 
     // Check if it is a Preshower segment - Layer 2
     // One segment per layer, nothing between
     if (cp[2 * is].whichDetector() == DetId::Ecal && cp[2 * is].whichSubDetector() == EcalPreshower &&
         cp[2 * is].whichLayer() == 2) {
       if (cp[2 * is + 1].whichDetector() == DetId::Ecal && cp[2 * is + 1].whichSubDetector() == EcalPreshower &&
           cp[2 * is + 1].whichLayer() == 2) {
         CaloSegment preshsegment(cp[2 * is], cp[2 * is + 1], s, sX0, sL0, CaloSegment::PS, myCalorimeter);
         segments_.push_back(preshsegment);
         //        std::cout << " Added (1-2)" << preshsegment << std::endl;
         s += preshsegment.length();
         sX0 += preshsegment.X0length();
         sL0 += preshsegment.L0length();
         X0PS2_ += preshsegment.X0length();
         L0PS2_ += preshsegment.L0length();
 
         // material between preshower and EE
         if (is < (nsegments - 1) && cp[2 * is + 2].whichDetector() == DetId::Ecal &&
             cp[2 * is + 2].whichSubDetector() == EcalEndcap) {
           CaloSegment gapsef(cp[2 * is + 1], cp[2 * is + 2], s, sX0, sL0, CaloSegment::PSEEGAP, myCalorimeter);
           segments_.push_back(gapsef);
           s += gapsef.length();
           sX0 += gapsef.X0length();
           sL0 += gapsef.L0length();
           X0PS2EE_ += gapsef.X0length();
           L0PS2EE_ += gapsef.L0length();
           //          std::cout << " Created  a segment " << gapsef.length()<< " " << gapsef.X0length()<< std::endl;
         }
       } else {
         std::cout << " Strange segment between Preshower2 and " << cp[2 * is + 1].whichDetector();
         std::cout << std::endl;
       }
       ++is;
       continue;
     }
     // Now deal with the ECAL
     // One segment in each crystal. Segment corresponding to cracks/gaps are added
     //      myHistos->debug("Just avant ECAL");
     if (cp[2 * is].whichDetector() == DetId::Ecal &&
         (cp[2 * is].whichSubDetector() == EcalBarrel || cp[2 * is].whichSubDetector() == EcalEndcap)) {
       if (cp[2 * is + 1].whichDetector() == DetId::Ecal &&
           (cp[2 * is + 1].whichSubDetector() == EcalBarrel || cp[2 * is + 1].whichSubDetector() == EcalEndcap)) {
         DetId cell2 = cp[2 * is + 1].getDetId();
         // set the real entrance
         if (ecalFirstSegment_ < 0)
           ecalFirstSegment_ = segments_.size();
 
         // !! Approximatiom : the first segment is always in a crystal
         if (cp[2 * is].getDetId() == cell2) {
           CaloSegment segment(cp[2 * is], cp[2 * is + 1], s, sX0, sL0, CaloSegment::PbWO4, myCalorimeter);
           segments_.push_back(segment);
           // std::cout << " Added (2)" << segment << std::endl;
           s += segment.length();
           sX0 += segment.X0length();
           sL0 += segment.L0length();
           X0ECAL_ += segment.X0length();
           L0ECAL_ += segment.L0length();
 #ifdef DEBUGCELLLINE
           ++ncrossedxtals;
 #endif
           ++is;
         } else {
           std::cout << " One more bug in the segment " << std::endl;
           ++is;
         }
         // Now check if a gap or crack should be added
         if (is > 0 && is < nsegments) {
           DetId cell3 = cp[2 * is].getDetId();
           if (cp[2 * is].whichDetector() != DetId::Hcal) {
             // Crack inside the ECAL
             bool bordercrossing = myCalorimeter->borderCrossing(cell2, cell3);
             CaloSegment cracksegment(cp[2 * is - 1],
                                      cp[2 * is],
                                      s,
                                      sX0,
                                      sL0,
                                      (bordercrossing) ? CaloSegment::CRACK : CaloSegment::GAP,
                                      myCalorimeter);
             segments_.push_back(cracksegment);
             s += cracksegment.length();
             sX0 += cracksegment.X0length();
             sL0 += cracksegment.L0length();
             X0ECAL_ += cracksegment.X0length();
             L0ECAL_ += cracksegment.L0length();
             //   std::cout <<" Added(3) "<< cracksegment << std::endl;
           } else {
             // a segment corresponding to ECAL/HCAL transition should be
             // added here
             CaloSegment cracksegment(cp[2 * is - 1], cp[2 * is], s, sX0, sL0, CaloSegment::ECALHCALGAP, myCalorimeter);
             segments_.push_back(cracksegment);
             s += cracksegment.length();
             sX0 += cracksegment.X0length();
             sL0 += cracksegment.L0length();
             X0EHGAP_ += cracksegment.X0length();
             L0EHGAP_ += cracksegment.L0length();
           }
         }
         continue;
       } else {
         std::cout << " Strange segment between " << cp[2 * is].whichDetector();
         std::cout << " and " << cp[2 * is + 1].whichDetector() << std::endl;
         ++is;
         continue;
       }
     }
     //      myHistos->debug("Just avant HCAL");
     // HCAL
     if (cp[2 * is].whichDetector() == DetId::Hcal && cp[2 * is + 1].whichDetector() == DetId::Hcal) {
       CaloSegment segment(cp[2 * is], cp[2 * is + 1], s, sX0, sL0, CaloSegment::HCAL, myCalorimeter);
       segments_.push_back(segment);
       s += segment.length();
       sX0 += segment.X0length();
       sL0 += segment.L0length();
       X0HCAL_ += segment.X0length();
       L0HCAL_ += segment.L0length();
       //      std::cout <<" Added(4) "<< segment << std::endl;
       ++is;
     }
   }
   //  std::cout << " PS1 " << X0PS1_ << " " << L0PS1_ << std::endl;
   //  std::cout << " PS2 " << X0PS2_ << " " << L0PS2_ << std::endl;
   //  std::cout << " ECAL " << X0ECAL_ << " " << L0ECAL_ << std::endl;
   //  std::cout << " HCAL " << X0HCAL_ << " " << L0HCAL_ << std::endl;
 
   totalX0_ = X0PS1_ + X0PS2_ + X0PS2EE_ + X0ECAL_ + X0EHGAP_ + X0HCAL_;
   totalL0_ = L0PS1_ + L0PS2_ + L0PS2EE_ + L0ECAL_ + L0EHGAP_ + L0HCAL_;
   //  myHistos->debug("Just avant le fill");
 
 #ifdef DEBUGCELLLINE
   myHistos->fill("h200", fabs(EcalEntrance_.eta()), X0ECAL_);
   myHistos->fill("h210", EcalEntrance_.phi(), X0ECAL_);
   if (X0ECAL_ < 20)
     myHistos->fill("h212", EcalEntrance_.phi(), X0ECAL_);
   //  if(X0ECAL_<1.)
   //    {
   //      for(unsigned ii=0; ii<segments_.size() ; ++ii)
   //    {
   //      std::cout << segments_[ii] << std::endl;
   //    }
   //    }
   myHistos->fillByNumber("h30", ncrossedxtals, EcalEntrance_.eta(), X0ECAL_);
 
   double zvertex = myTrack_->vertex().position().z();
 
   myHistos->fill("h310", EcalEntrance_.eta(), X0ECAL_);
   if (X0ECAL_ < 22)
     myHistos->fill("h410", EcalEntrance_.phi());
   myHistos->fill("h400", zvertex, X0ECAL_);
 #endif
   //  std::cout << " Finished the segments " << std::endl;
 }

◆ cellLine()

void EcalHitMaker::cellLine ( std::vector< CaloPoint > & cp )

private

Definition at line 369 of file EcalHitMaker.cc.

References central_, gather_cfg::cout, DeadROC_duringRun::dir, ecalCellLine(), FSimTrack::ecalEntrance(), EcalEntrance_, hcalCellLine(), FSimTrack::hcalEntrance(), FSimTrack::layer1Entrance(), FSimTrack::layer2Entrance(), myTrack_, normal_, FSimTrack::onEcal(), onEcal_, FSimTrack::onHcal(), FSimTrack::onLayer1(), FSimTrack::onLayer2(), FSimTrack::onVFcal(), FSimVertex::position(), preshowerCellLine(), simulatePreshower_, jetUpdater_cfi::sort, RawParticle::Vect(), bphysicsOniaDQM_cfi::vertex, FSimTrack::vertex(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by setTrackParameters().

                                                     {
   cp.clear();
   //  if(myTrack->onVFcal()!=2)
   //    {
   if (!central_ && onEcal_ && simulatePreshower_)
     preshowerCellLine(cp);
   if (onEcal_)
     ecalCellLine(EcalEntrance_, EcalEntrance_ + normal_, cp);
   //    }
 
   XYZPoint vertex(myTrack_->vertex().position().Vect());
 
   //sort the points by distance (in the ECAL they are not necessarily ordered)
   XYZVector dir(0., 0., 0.);
   if (myTrack_->onLayer1()) {
     vertex = (myTrack_->layer1Entrance().vertex()).Vect();
     dir = myTrack_->layer1Entrance().Vect().Unit();
   } else if (myTrack_->onLayer2()) {
     vertex = (myTrack_->layer2Entrance().vertex()).Vect();
     dir = myTrack_->layer2Entrance().Vect().Unit();
   } else if (myTrack_->onEcal()) {
     vertex = (myTrack_->ecalEntrance().vertex()).Vect();
     dir = myTrack_->ecalEntrance().Vect().Unit();
   } else if (myTrack_->onHcal()) {
     vertex = (myTrack_->hcalEntrance().vertex()).Vect();
     dir = myTrack_->hcalEntrance().Vect().Unit();
   } else if (myTrack_->onVFcal() == 2) {
     vertex = (myTrack_->vfcalEntrance().vertex()).Vect();
     dir = myTrack_->vfcalEntrance().Vect().Unit();
   } else {
     std::cout << " Problem with the grid " << std::endl;
   }
 
   // Move the vertex for distance comparison (5cm)
   vertex -= 5. * dir;
   CaloPoint::DistanceToVertex myDistance(vertex);
   sort(cp.begin(), cp.end(), myDistance);
 
   // The intersections with the HCAL shouldn't need to be sorted
   // with the N.I it is actually a source of problems
   hcalCellLine(cp);
 
   //  std::cout << " Intersections ordered by distance to " << vertex << std::endl;
   //
   //  for (unsigned ic=0;ic<cp.size();++ic)
   //    {
   //      XYZVector t=cp[ic]-vertex;
   //      std::cout << cp[ic] << " " << t.mag() << std::endl;
   //    }
 }

◆ configureGeometry()

void EcalHitMaker::configureGeometry ( )

private

Definition at line 1010 of file EcalHitMaker.cc.

References bfactor_, CaloGeometryHelper::borderCrossing(), configuredGeometry_, DetIdMap_, DeadROC_duringRun::dir, EcalEntrance_, M_PI, M_PI_2, CaloGeometryHelper::magneticField(), CaloGeometryHelper::move(), CaloHitMaker::myCalorimeter, ncrystals_, CaloDirectionOperations::neighbourDirection(), CaloDirectionOperations::oppositeDirection(), regionOfInterest_, mps_update::status, and theta().

Referenced by getPads().

                                      {
   configuredGeometry_ = true;
   for (unsigned ic = 0; ic < ncrystals_; ++ic) {
     //      std::cout << " Building " << cellids_[ic] << std::endl;
     for (unsigned idir = 0; idir < 8; ++idir) {
       unsigned oppdir = CaloDirectionOperations::oppositeDirection(idir);
       // Is there something else to do ?
       // The relationship with the neighbour may have been set previously.
       if (regionOfInterest_[ic].crystalNeighbour(idir).status() >= 0) {
         //        std::cout << " Nothing to do " << std::endl;
         continue;
       }
 
       const DetId& oldcell(regionOfInterest_[ic].getDetId());
       CaloDirection dir = CaloDirectionOperations::neighbourDirection(idir);
       DetId newcell(oldcell);
       if (!myCalorimeter->move(newcell, dir)) {
         // no neighbour in this direction
         regionOfInterest_[ic].crystalNeighbour(idir).setStatus(-1);
         continue;
       }
       // Determine the number of this neighbour
       //      std::cout << " The neighbour is " << newcell << std::endl;
       std::map<DetId, unsigned>::const_iterator niter(DetIdMap_.find(newcell));
       if (niter == DetIdMap_.end()) {
         //        std::cout << " The neighbour is not in the map " << std::endl;
         regionOfInterest_[ic].crystalNeighbour(idir).setStatus(-1);
         continue;
       }
       // Now there is a neighbour
       //      std::cout << " The neighbour is " << niter->second << " " << cellids_[niter->second] << std::endl;
       regionOfInterest_[ic].crystalNeighbour(idir).setNumber(niter->second);
       //      std::cout << " Managed to set crystalNeighbour " << ic << " " << idir << std::endl;
       //      std::cout << " Trying  " << niter->second << " " << oppdir << std::endl;
       regionOfInterest_[niter->second].crystalNeighbour(oppdir).setNumber(ic);
       //      std::cout << " Crack/gap " << std::endl;
       if (myCalorimeter->borderCrossing(oldcell, newcell)) {
         regionOfInterest_[ic].crystalNeighbour(idir).setStatus(1);
         regionOfInterest_[niter->second].crystalNeighbour(oppdir).setStatus(1);
         //        std::cout << " Crack ! " << std::endl;
       } else {
         regionOfInterest_[ic].crystalNeighbour(idir).setStatus(0);
         regionOfInterest_[niter->second].crystalNeighbour(oppdir).setStatus(0);
         //        std::cout << " Gap" << std::endl;
       }
     }
   }
   // Magnetic field a la Charlot
   double theta = EcalEntrance_.theta();
   if (theta > M_PI_2)
     theta = M_PI - theta;
   bfactor_ = 1. / (1. + 0.133 * theta);
   if (myCalorimeter->magneticField() == 0.)
     bfactor_ = 1.;
 }

◆ convertIntegerCoordinates()

void EcalHitMaker::convertIntegerCoordinates	(	double	x,
		double	y,
		unsigned &	ix,
		unsigned &	iy
	)		const

private

Definition at line 1088 of file EcalHitMaker.cc.

References createfilelist::int, sizex_, sizey_, x, xmin_, y, and ymin_.

Referenced by fastInsideCell(), and prepareCrystalNumberArray().

                                                                                                  {
   int tix = (int)((x - xmin_) / sizex_);
   int tiy = (int)((y - ymin_) / sizey_);
   ix = iy = 9999;
   if (tix >= 0)
     ix = (unsigned)tix;
   if (tiy >= 0)
     iy = (unsigned)tiy;
 }

◆ correspondingEdge()

CLHEP::Hep2Vector & EcalHitMaker::correspondingEdge	(	neighbour &	myneighbour,
		CaloDirection	dir2
	)

private

Definition at line 1176 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), TCMET_cfi::corner, DeadROC_duringRun::dir, L1TDiffHarvesting_cfi::dir2, CaloDirectionOperations::oppositeSide(), and padsatdepth_.

Referenced by gapsLifting().

                                                                                            {
   CaloDirection dir = CaloDirectionOperations::oppositeSide(myneighbour.first);
   CaloDirection corner = CaloDirectionOperations::add2d(dir, dir2);
   //  std::cout << "Corresponding Edge " << dir<< " " << dir2 << " " << corner << std::endl;
   return padsatdepth_[myneighbour.second].edge(corner);
 }

◆ cracksPads()

void EcalHitMaker::cracksPads	(	std::vector< neighbour > &	cracks,
		unsigned	iq
	)

private

Definition at line 1336 of file EcalHitMaker.cc.

References crackPadProbability_, crackpadsatdepth_, EAST, CrystalPad::edge(), dqmdumpme::first, mycorners, NORTH, NORTHEAST, NORTHWEST, padsatdepth_, CrystalPad::setSurvivalProbability(), SOUTH, SOUTHEAST, SOUTHWEST, and WEST.

Referenced by reorganizePads().

                                                                        {
   //  std::cout << " myPad " << &myPad << std::endl;
   unsigned ncracks = cracks.size();
   CrystalPad& myPad = padsatdepth_[iq];
   for (unsigned ic = 0; ic < ncracks; ++ic) {
     //      std::vector<CLHEP::Hep2Vector> mycorners;
     //      mycorners.reserve(4);
     switch (cracks[ic].first) {
       case NORTH: {
         mycorners[0] = (padsatdepth_[cracks[ic].second].edge(SOUTHWEST));
         mycorners[1] = (padsatdepth_[cracks[ic].second].edge(SOUTHEAST));
         mycorners[2] = (myPad.edge(NORTHEAST));
         mycorners[3] = (myPad.edge(NORTHWEST));
       } break;
       case SOUTH: {
         mycorners[0] = (myPad.edge(SOUTHWEST));
         mycorners[1] = (myPad.edge(SOUTHEAST));
         mycorners[2] = (padsatdepth_[cracks[ic].second].edge(NORTHEAST));
         mycorners[3] = (padsatdepth_[cracks[ic].second].edge(NORTHWEST));
       } break;
       case EAST: {
         mycorners[0] = (myPad.edge(NORTHEAST));
         mycorners[1] = (padsatdepth_[cracks[ic].second].edge(NORTHWEST));
         mycorners[2] = (padsatdepth_[cracks[ic].second].edge(SOUTHWEST));
         mycorners[3] = (myPad.edge(SOUTHEAST));
       } break;
       case WEST: {
         mycorners[0] = (padsatdepth_[cracks[ic].second].edge(NORTHEAST));
         mycorners[1] = (myPad.edge(NORTHWEST));
         mycorners[2] = (myPad.edge(SOUTHWEST));
         mycorners[3] = (padsatdepth_[cracks[ic].second].edge(SOUTHEAST));
       } break;
       default: {
       }
     }
     CrystalPad crackpad(ic, mycorners);
     // to be tuned. A simpleconfigurable should be used
     crackpad.setSurvivalProbability(crackPadProbability_);
     crackpadsatdepth_.push_back(crackpad);
   }
   //  std::cout << " Finished cracksPads " << std::endl;
 }

◆ diagonalEdge()

bool EcalHitMaker::diagonalEdge	(	unsigned	myPad,
		CaloDirection	dir,
		CLHEP::Hep2Vector &	point
	)

private

Definition at line 1183 of file EcalHitMaker.cc.

References DeadROC_duringRun::dir, CaloDirectionOperations::neighbourDirection(), CaloDirectionOperations::oppositeSide(), padsatdepth_, point, regionOfInterest_, mps_update::status, and validPads_.

Referenced by gapsLifting().

                                                                                          {
   unsigned idir = CaloDirectionOperations::neighbourDirection(dir);
   if (regionOfInterest_[myPad].crystalNeighbour(idir).status() < 0)
     return false;
   unsigned nneighbour = regionOfInterest_[myPad].crystalNeighbour(idir).number();
   if (!validPads_[nneighbour]) {
     //      std::cout << " Wasn't able to move " << std::endl;
     return false;
   }
   point = padsatdepth_[nneighbour].edge(CaloDirectionOperations::oppositeSide(dir));
   return true;
 }

◆ ecalCellLine()

void EcalHitMaker::ecalCellLine	(	const XYZPoint &	a,
		const XYZPoint &	b,
		std::vector< CaloPoint > &	cp
	)

private

Definition at line 497 of file EcalHitMaker.cc.

References a, angle(), b, corners, DOWN, inside3D(), CaloHitMaker::intersect(), mag2(), ncrystals_, normal_, regionOfInterest_, CaloDirectionOperations::Side(), mathSSE::sqrt(), submitPVValidationJobs::t, funct::tan(), and UP.

Referenced by cellLine().

                                                                                               {
   //  std::vector<XYZPoint> corners;
   //  corners.resize(4);
   unsigned ic = 0;
   double t;
   XYZPoint xp;
   DetId c_entrance, c_exit;
   bool entrancefound(false), exitfound(false);
   //  std::cout << " Look for intersections " << ncrystals_ << std::endl;
   //  std::cout << " regionOfInterest_ " << truncatedGrid_ << " " << regionOfInterest_.size() << std::endl;
   // try to determine the number of crystals to test
   // First determine the incident angle
   double angle = std::acos(normal_.Dot(regionOfInterest_[0].getAxis().Unit()));
 
   //  std::cout << " Normal " << normal_<< " Axis " << regionOfInterest_[0].getAxis().Unit() << std::endl;
   double backdistance = std::sqrt(regionOfInterest_[0].getAxis().mag2()) * std::tan(angle);
   // 1/2.2cm = 0.45
   //   std::cout << " Angle " << angle << std::endl;
   //   std::cout << " Back distance " << backdistance << std::endl;
   unsigned ncrystals = (unsigned)(backdistance * 0.45);
   unsigned highlim = (ncrystals + 4);
   highlim *= highlim;
   if (highlim > ncrystals_)
     highlim = ncrystals_;
   //  unsigned lowlim=(ncrystals>2)? (ncrystals-2):0;
   //  std::cout << " Ncrys " << ncrystals << std::endl;
 
   while (ic < ncrystals_ && (ic < highlim || !exitfound)) {
     // Check front side
     //      if(!entrancefound)
     {
       const Plane3D& plan = regionOfInterest_[ic].getFrontPlane();
       //      XYZVector axis1=(plan.Normal());
       //      XYZVector axis2=regionOfInterest_[ic].getFirstEdge();
       xp = intersect(plan, a, b, t, false);
       regionOfInterest_[ic].getFrontSide(corners);
       //      CrystalPad pad(9999,onEcal_,corners,regionOfInterest_[ic].getCorner(0),axis1,axis2);
       //      if(pad.globalinside(xp))
       if (inside3D(corners, xp)) {
         cp.push_back(CaloPoint(regionOfInterest_[ic].getDetId(), UP, xp));
         entrancefound = true;
         c_entrance = regionOfInterest_[ic].getDetId();
         //      myHistos->fill("j12",highlim,ic);
       }
     }
 
     // check rear side
     //  if(!exitfound)
     {
       const Plane3D& plan = regionOfInterest_[ic].getBackPlane();
       //      XYZVector axis1=(plan.Normal());
       //      XYZVector axis2=regionOfInterest_[ic].getFifthEdge();
       xp = intersect(plan, a, b, t, false);
       regionOfInterest_[ic].getBackSide(corners);
       //      CrystalPad pad(9999,onEcal_,corners,regionOfInterest_[ic].getCorner(4),axis1,axis2);
       //      if(pad.globalinside(xp))
       if (inside3D(corners, xp)) {
         cp.push_back(CaloPoint(regionOfInterest_[ic].getDetId(), DOWN, xp));
         exitfound = true;
         c_exit = regionOfInterest_[ic].getDetId();
         //        std::cout << " Crystal : " << ic << std::endl;
         //        myHistos->fill("j10",highlim,ic);
       }
     }
 
     if (entrancefound && exitfound && c_entrance == c_exit)
       return;
     // check lateral sides
     for (unsigned iside = 0; iside < 4; ++iside) {
       const Plane3D& plan = regionOfInterest_[ic].getLateralPlane(iside);
       xp = intersect(plan, a, b, t, false);
       //      XYZVector axis1=(plan.Normal());
       //      XYZVector axis2=regionOfInterest_[ic].getLateralEdge(iside);
       regionOfInterest_[ic].getLateralSide(iside, corners);
       //      CrystalPad pad(9999,onEcal_,corners,regionOfInterest_[ic].getCorner(iside),axis1,axis2);
       //      if(pad.globalinside(xp))
       if (inside3D(corners, xp)) {
         cp.push_back(CaloPoint(regionOfInterest_[ic].getDetId(), CaloDirectionOperations::Side(iside), xp));
         //        std::cout << cp[cp.size()-1] << std::endl;
       }
     }
     // Go to next crystal
     ++ic;
   }
 }

◆ ecalEntrance()

const XYZPoint& EcalHitMaker::ecalEntrance ( ) const

inline

used in FamosHcalHitMaker

Definition at line 123 of file EcalHitMaker.h.

References EcalEntrance_.

Referenced by HcalHitMaker::HcalHitMaker().

123 { return EcalEntrance_; };

EcalHitMaker::EcalEntrance_

XYZPoint EcalEntrance_

Definition: EcalHitMaker.h:213

◆ ecalHcalGapTotalL0()

double EcalHitMaker::ecalHcalGapTotalL0 ( ) const

inline

ECAL-HCAL transition.

Definition at line 91 of file EcalHitMaker.h.

References L0EHGAP_.

Referenced by HDRShower::computeShower(), HDShower::HDShower(), and HFShower::HFShower().

91 { return L0EHGAP_; }

EcalHitMaker::L0EHGAP_

double L0EHGAP_

Definition: EcalHitMaker.h:205

◆ ecalHcalGapTotalX0()

double EcalHitMaker::ecalHcalGapTotalX0 ( ) const

inline

ECAL-HCAL transition.

Definition at line 70 of file EcalHitMaker.h.

References X0EHGAP_.

Referenced by HDShower::HDShower(), and HFShower::HFShower().

70 { return X0EHGAP_; }

EcalHitMaker::X0EHGAP_

double X0EHGAP_

Definition: EcalHitMaker.h:198

◆ ecalTotalL0()

double EcalHitMaker::ecalTotalL0 ( ) const

inline

in the ECAL

Definition at line 85 of file EcalHitMaker.h.

References L0ECAL_.

Referenced by HDRShower::computeShower(), HDShower::HDShower(), and HFShower::HFShower().

85 { return L0ECAL_; }

EcalHitMaker::L0ECAL_

double L0ECAL_

Definition: EcalHitMaker.h:203

◆ ecalTotalX0()

double EcalHitMaker::ecalTotalX0 ( ) const

inline

in the ECAL

Definition at line 67 of file EcalHitMaker.h.

References X0ECAL_.

Referenced by EMShower::prepareSteps(), and setTrackParameters().

67 { return X0ECAL_; }

EcalHitMaker::X0ECAL_

double X0ECAL_

Definition: EcalHitMaker.h:197

◆ fastInsideCell()

unsigned EcalHitMaker::fastInsideCell	(	const CLHEP::Hep2Vector &	point,
		double &	sp,
		bool	debug = `false`
	)

Definition at line 229 of file EcalHitMaker.cc.

References convertIntegerCoordinates(), gather_cfg::cout, crackpadsatdepth_, d1, debug, newFWLiteAna::found, CrystalWindowMap::getCrystalWindow(), myCrystalNumberArray_, myCrystalWindowMap_, ncrackpadsatdepth_, nx_, ny_, padsatdepth_, point, findQualityFiles::size, mps_update::status, and validPads_.

Referenced by addHit(), and addHitDepth().

                                                                                           {
   //  debug = true;
   bool found = false;
   // something clever has to be implemented here
   unsigned d1, d2;
   convertIntegerCoordinates(point.x(), point.y(), d1, d2);
   //  std::cout << "Fastinside cell " << point.x() << " " <<  point.y() << " " << d1 << " "<< d2 << " " << nx_ << " " << ny_ << std::endl;
   if (d1 >= nx_ || d2 >= ny_) {
     //      std::cout << " Not in the map " <<std::endl;
     return 9999;
   }
   unsigned cell = myCrystalNumberArray_[d1][d2];
   // We are likely to be lucky
   //  std::cout << " Got the cell " << cell << std::endl;
   if (validPads_[cell] && padsatdepth_[cell].inside(point)) {
     //      std::cout << " We are lucky " << cell << std::endl;
     sp = padsatdepth_[cell].survivalProbability();
     return cell;
   }
 
   //  std::cout << "Starting the loop " << std::endl;
   bool status(true);
   const std::vector<unsigned>& localCellVector(myCrystalWindowMap_->getCrystalWindow(cell, status));
   if (status) {
     unsigned size = localCellVector.size();
     //      std::cout << " Starting from " << EBDetId(regionOfInterest_[cell].getDetId()) << std::endl;
     //      const std::vector<DetId>& neighbours=myCalorimeter->getNeighbours(regionOfInterest_[cell].getDetId());
     //      std::cout << " The neighbours are " << std::endl;
     //      for(unsigned ic=0;ic<neighbours.size(); ++ic)
     //  {
     //    std::cout << EBDetId(neighbours[ic]) << std::endl;
     //  }
     //      std::cout << " Done " << std::endl;
     for (unsigned ic = 0; ic < 8 && ic < size; ++ic) {
       unsigned iq = localCellVector[ic];
       //      std::cout << " Testing " << EBDetId(regionOfInterest_[iq].getDetId()) << std::endl; ;
       //      std::cout << " " << iq << std::endl;
       //      std::cout << padsatdepth_[iq] ;
       if (validPads_[iq] && padsatdepth_[iq].inside(point)) {
         //        std::cout << " Yes " << std::endl;
         //        myHistos->fill("h1000",niter);
         sp = padsatdepth_[iq].survivalProbability();
         //        std::cout << "Finished the loop " << niter << std::endl;
         //        std::cout << "Inside " << std::endl;
         return iq;
       }
       //      std::cout << " Not inside " << std::endl;
       //      std::cout << "No " << std::endl;
     }
   }
   if (debug)
     std::cout << " not found in a quad, let's check the " << ncrackpadsatdepth_ << " cracks " << std::endl;
   //  std::cout << "Finished the loop " << niter << std::endl;
   // Let's check the cracks
   //  std::cout << " Let's check the cracks " << ncrackpadsatdepth_ << " " << crackpadsatdepth_.size() << std::endl;
   unsigned iquad = 0;
   unsigned iquadinside = 999;
   while (iquad < ncrackpadsatdepth_ && !found) {
     //      std::cout << " Inside the while " << std::endl;
     if (crackpadsatdepth_[iquad].inside(point)) {
       iquadinside = iquad;
       found = true;
       sp = crackpadsatdepth_[iquad].survivalProbability();
     }
     ++iquad;
   }
   //  myHistos->fill("h1002",niter);
   if (!found && debug)
     std::cout << " Not found in the cracks " << std::endl;
   return (found) ? crackpadsatdepth_[iquadinside].getNumber() : 9999;
 }

◆ gapsLifting()

void EcalHitMaker::gapsLifting	(	std::vector< neighbour > &	gaps,
		unsigned	iq
	)

private

Definition at line 1225 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), correspondingEdge(), gather_cfg::cout, debug, diagonalEdge(), L1TDiffHarvesting_cfi::dir1, L1TDiffHarvesting_cfi::dir2, EAST, CrystalPad::edge(), dqmdumpme::first, NONE, NORTH, NORTHEAST, NORTHWEST, CaloDirectionOperations::oppositeSide(), padsatdepth_, point, SOUTH, SOUTHEAST, SOUTHWEST, unbalancedDirection(), and WEST.

Referenced by reorganizePads().

                                                                       {
   //  std::cout << " Entering gapsLifting "  << std::endl;
   CrystalPad& myPad = padsatdepth_[iq];
   unsigned ngaps = gaps.size();
   constexpr bool debug = false;
   if (ngaps == 1) {
     if (debug) {
       std::cout << " Avant " << ngaps << " " << gaps[0].first << std::endl;
       std::cout << myPad << std::endl;
     }
     if (gaps[0].first == NORTH || gaps[0].first == SOUTH) {
       CaloDirection dir1 = CaloDirectionOperations::add2d(gaps[0].first, EAST);
       CaloDirection dir2 = CaloDirectionOperations::add2d(gaps[0].first, WEST);
       myPad.edge(dir1) = correspondingEdge(gaps[0], EAST);
       myPad.edge(dir2) = correspondingEdge(gaps[0], WEST);
     } else {
       CaloDirection dir1 = CaloDirectionOperations::add2d(gaps[0].first, NORTH);
       CaloDirection dir2 = CaloDirectionOperations::add2d(gaps[0].first, SOUTH);
       myPad.edge(dir1) = correspondingEdge(gaps[0], NORTH);
       myPad.edge(dir2) = correspondingEdge(gaps[0], SOUTH);
     }
     if (debug) {
       std::cout << " Apres " << std::endl;
       std::cout << myPad << std::endl;
     }
   } else if (ngaps == 2) {
     if (debug) {
       std::cout << " Avant " << ngaps << " " << gaps[0].first << " " << gaps[1].first << std::endl;
       std::cout << myPad << std::endl;
       std::cout << " Voisin 1 " << (gaps[0].second) << std::endl;
       std::cout << " Voisin 2 " << (gaps[1].second) << std::endl;
     }
     CaloDirection corner0 = CaloDirectionOperations::add2d(gaps[0].first, gaps[1].first);
 
     CLHEP::Hep2Vector point(0., 0.);
     if (corner0 != NONE && diagonalEdge(iq, corner0, point)) {
       CaloDirection corner1 =
           CaloDirectionOperations::add2d(CaloDirectionOperations::oppositeSide(gaps[0].first), gaps[1].first);
       CaloDirection corner2 =
           CaloDirectionOperations::add2d(gaps[0].first, CaloDirectionOperations::oppositeSide(gaps[1].first));
       myPad.edge(corner0) = point;
       myPad.edge(corner1) = correspondingEdge(gaps[1], CaloDirectionOperations::oppositeSide(gaps[0].first));
       myPad.edge(corner2) = correspondingEdge(gaps[0], CaloDirectionOperations::oppositeSide(gaps[1].first));
     } else if (corner0 == NONE) {
       if (gaps[0].first == EAST || gaps[0].first == WEST) {
         CaloDirection corner1 = CaloDirectionOperations::add2d(gaps[0].first, NORTH);
         CaloDirection corner2 = CaloDirectionOperations::add2d(gaps[0].first, SOUTH);
         myPad.edge(corner1) = correspondingEdge(gaps[0], NORTH);
         myPad.edge(corner2) = correspondingEdge(gaps[0], SOUTH);
 
         corner1 = CaloDirectionOperations::add2d(gaps[1].first, NORTH);
         corner2 = CaloDirectionOperations::add2d(gaps[1].first, SOUTH);
         myPad.edge(corner1) = correspondingEdge(gaps[1], NORTH);
         myPad.edge(corner2) = correspondingEdge(gaps[1], SOUTH);
       } else {
         CaloDirection corner1 = CaloDirectionOperations::add2d(gaps[0].first, EAST);
         CaloDirection corner2 = CaloDirectionOperations::add2d(gaps[0].first, WEST);
         myPad.edge(corner1) = correspondingEdge(gaps[0], EAST);
         myPad.edge(corner2) = correspondingEdge(gaps[0], WEST);
 
         corner1 = CaloDirectionOperations::add2d(gaps[1].first, EAST);
         corner2 = CaloDirectionOperations::add2d(gaps[1].first, WEST);
         myPad.edge(corner1) = correspondingEdge(gaps[1], EAST);
         myPad.edge(corner2) = correspondingEdge(gaps[1], WEST);
       }
     }
     if (debug) {
       std::cout << " Apres " << std::endl;
       std::cout << myPad << std::endl;
     }
   } else if (ngaps == 3) {
     // in this case the four corners have to be changed
     unsigned iubd, idir1, idir2;
     CaloDirection diag;
     CLHEP::Hep2Vector point(0., 0.);
     //    std::cout << " Yes : 3 gaps" << std::endl;
     if (unbalancedDirection(gaps, iubd, idir1, idir2)) {
       CaloDirection ubd(gaps[iubd].first), dir1(gaps[idir1].first);
       CaloDirection dir2(gaps[idir2].first);
 
       //          std::cout << " Avant " << std::endl << myPad << std::endl;
       //          std::cout << ubd << " " << dir1 << " " << dir2 << std::endl;
       diag = CaloDirectionOperations::add2d(ubd, dir1);
       if (diagonalEdge(iq, diag, point))
         myPad.edge(diag) = point;
       diag = CaloDirectionOperations::add2d(ubd, dir2);
       if (diagonalEdge(iq, diag, point))
         myPad.edge(diag) = point;
       CaloDirection oppside = CaloDirectionOperations::oppositeSide(ubd);
       myPad.edge(CaloDirectionOperations::add2d(oppside, dir1)) = correspondingEdge(gaps[idir1], oppside);
       myPad.edge(CaloDirectionOperations::add2d(oppside, dir2)) = correspondingEdge(gaps[idir2], oppside);
       //          std::cout << " Apres " << std::endl << myPad << std::endl;
     }
   } else if (ngaps == 4) {
     //      std::cout << " Waouh :4 gaps" << std::endl;
     //      std::cout << " Avant " << std::endl;
     //      std::cout << myPad<< std::endl;
     CLHEP::Hep2Vector point(0., 0.);
     if (diagonalEdge(iq, NORTHEAST, point))
       myPad.edge(NORTHEAST) = point;
     if (diagonalEdge(iq, NORTHWEST, point))
       myPad.edge(NORTHWEST) = point;
     if (diagonalEdge(iq, SOUTHWEST, point))
       myPad.edge(SOUTHWEST) = point;
     if (diagonalEdge(iq, SOUTHEAST, point))
       myPad.edge(SOUTHEAST) = point;
     //      std::cout << " Apres " << std::endl;
     //      std::cout << myPad<< std::endl;
   }
 }

◆ getCrystals()

const std::vector<Crystal>& EcalHitMaker::getCrystals ( ) const

inline

for debugging

Definition at line 135 of file EcalHitMaker.h.

References regionOfInterest_.

135 { return regionOfInterest_; }

EcalHitMaker::regionOfInterest_

std::vector< Crystal > regionOfInterest_

Definition: EcalHitMaker.h:225

◆ getFSimTrack()

const FSimTrack* EcalHitMaker::getFSimTrack ( ) const

inline

To retrieve the track.

Definition at line 120 of file EcalHitMaker.h.

References myTrack_.

Referenced by HcalHitMaker::setDepth().

120 { return myTrack_; }

EcalHitMaker::myTrack_

const FSimTrack * myTrack_

Definition: EcalHitMaker.h:271

◆ getHits()

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

overridevirtual

get the map of the stored hits. Triggers the calculation of the grid if it has not been done.

Implements CaloHitMaker.

Definition at line 1098 of file EcalHitMaker.cc.

References Calorimeter::getEcalGeometry(), CaloSubdetectorGeometry::getGeometry(), ecalpyutils::hashedIndex(), CaloHitMaker::hitMap_, hitmaphasbeencalculated_, hits_, CaloHitMaker::myCalorimeter, ncrystals_, onEcal_, nano_mu_digi_cff::rawId, and regionOfInterest_.

Referenced by CalorimetryManager::EMShowerSimulation().

                                                       {
   if (hitmaphasbeencalculated_)
     return hitMap_;
   for (unsigned ic = 0; ic < ncrystals_; ++ic) {
     //calculate time of flight
     float tof = 0.0;
     if (onEcal_ == 1 || onEcal_ == 2)
       tof =
           (myCalorimeter->getEcalGeometry(onEcal_)->getGeometry(regionOfInterest_[ic].getDetId())->getPosition().mag()) /
           29.98;  //speed of light
 
     if (onEcal_ == 1) {
       CaloHitID current_id(EBDetId(regionOfInterest_[ic].getDetId().rawId()).hashedIndex(), tof, 0);  //no track yet
       hitMap_.insert(std::pair<CaloHitID, float>(current_id, hits_[ic]));
     } else if (onEcal_ == 2) {
       CaloHitID current_id(EEDetId(regionOfInterest_[ic].getDetId().rawId()).hashedIndex(), tof, 0);  //no track yet
       hitMap_.insert(std::pair<CaloHitID, float>(current_id, hits_[ic]));
     }
   }
   hitmaphasbeencalculated_ = true;
   return hitMap_;
 }

◆ getPads()

bool EcalHitMaker::getPads	(	double	depth,
		bool	inCm = `false`
	)

computes the crystals-plan intersection at depth (in X0 or L0 depending on the shower type) if it is not possible to go at such a depth, the result is false

Definition at line 821 of file EcalHitMaker.cc.

References a, b, bfactor_, central_, configuredGeometry_, configureGeometry(), corners, gather_cfg::cout, crackpadsatdepth_, currentdepth_, hcalRecHitTable_cff::depth, detailedShowerTail_, doreorg_, DetId::Ecal, CaloHitMaker::EMSHOWER, Crystal::getFirstEdge(), CaloHitMaker::HADSHOWER, cuy::ii, CaloHitMaker::interactionLength, CaloHitMaker::intersect(), CaloHitMaker::moliereRadius, myTrack_, ncrackpadsatdepth_, ncrystals_, normal_, npadsatdepth_, nx_, ny_, padsatdepth_, pivot_, plan_, prepareCrystalNumberArray(), pulled(), pulledPadProbability_, radiusCorrectionFactor_, radiusFactor_, regionOfInterest_, reorganizePads(), segments_, Validation_hcalonly_cfi::sign, sizex_, sizey_, validPads_, X0depthoffset_, xmax_, xmin_, geometryCSVtoXML::xx, ymax_, and ymin_.

Referenced by HFShower::compute(), HDShower::compute(), EMShower::compute(), and HDRShower::computeShower().

                                                   {
   //std::cout << " New depth " << depth << std::endl;
   // The first time, the relationship between crystals must be calculated
   // but only in the case of EM showers
 
   if (EMSHOWER && !configuredGeometry_)
     configureGeometry();
 
   radiusFactor_ = (EMSHOWER) ? moliereRadius * radiusCorrectionFactor_ : interactionLength;
   detailedShowerTail_ = false;
   if (EMSHOWER)
     currentdepth_ = depth + X0depthoffset_;
   else
     currentdepth_ = depth;
 
   //  if(currentdepth_>maxX0_+ps1TotalX0()+ps2TotalX0())
   //    {
   //      currentdepth_=maxX0_+ps1TotalX0()+ps2TotalX0()-1.; // the -1 is for safety
   //      detailedShowerTail_=true;
   //    }
 
   //  std::cout << " FamosGrid::getQuads " << currentdepth_ << " " << maxX0_ << std::endl;
 
   ncrackpadsatdepth_ = 0;
 
   xmin_ = ymin_ = 999;
   xmax_ = ymax_ = -999;
   double locxmin, locxmax, locymin, locymax;
 
   // Get the depth of the pivot
   std::vector<CaloSegment>::const_iterator segiterator;
   // First identify the correct segment
 
   if (inCm)  // centimeter
   {
     segiterator = find_if(segments_.begin(), segments_.end(), CaloSegment::inSegment(currentdepth_));
   } else {
     // EM shower
     if (EMSHOWER)
       segiterator = find_if(segments_.begin(), segments_.end(), CaloSegment::inX0Segment(currentdepth_));
 
     //Hadron shower
     if (HADSHOWER)
       segiterator = find_if(segments_.begin(), segments_.end(), CaloSegment::inL0Segment(currentdepth_));
   }
   if (segiterator == segments_.end()) {
     std::cout << " FamosGrid: Could not go at such depth " << depth << std::endl;
     std::cout << " EMSHOWER " << EMSHOWER << std::endl;
     std::cout << " Track " << *myTrack_ << std::endl;
     std::cout << " Segments " << segments_.size() << std::endl;
     for (unsigned ii = 0; ii < segments_.size(); ++ii) {
       std::cout << segments_[ii] << std::endl;
     }
 
     return false;
   }
   //  std::cout << *segiterator << std::endl;
 
   if (segiterator->whichDetector() != DetId::Ecal) {
     std::cout << " In  " << segiterator->whichDetector() << std::endl;
     //      buildPreshower();
     return false;
   }
 
   //  std::cout << *segiterator << std::endl;
   // get the position of the origin
 
   XYZPoint origin;
   if (inCm) {
     origin = segiterator->positionAtDepthincm(currentdepth_);
   } else {
     if (EMSHOWER)
       origin = segiterator->positionAtDepthinX0(currentdepth_);
     if (HADSHOWER)
       origin = segiterator->positionAtDepthinL0(currentdepth_);
   }
   //  std::cout << " currentdepth_ " << currentdepth_ << " " << origin << std::endl;
   XYZVector newaxis = pivot_.getFirstEdge().Cross(normal_);
 
   //  std::cout  << "Normal " << normal_ << std::endl;
   //  std::cout << " New axis " << newaxis << std::endl;
 
   //  std::cout << " ncrystals  " << ncrystals << std::endl;
   plan_ = Plane3D((Vector)normal_, (Point)origin);
 
   unsigned nquads = 0;
   double sign = (central_) ? -1. : 1.;
   Transform3DR trans((Point)origin,
                      (Point)(origin + normal_),
                      (Point)(origin + newaxis),
                      Point(0, 0, 0),
                      Point(0., 0., sign),
                      Point(0., 1., 0.));
   for (unsigned ic = 0; ic < ncrystals_; ++ic) {
     //      std::cout << " Building geometry for " << regionOfInterest_[ic].getCellID() << std::endl;
     XYZPoint a, b;
 
     //      std::cout << " Origin " << origin << std::endl;
 
     //      std::vector<XYZPoint> corners;
     //      corners.reserve(4);
     double dummyt;
     bool hasbeenpulled = false;
     bool behindback = false;
     for (unsigned il = 0; il < 4; ++il) {
       // a is the il-th front corner of the crystal. b is the corresponding rear corner
       regionOfInterest_[ic].getLateralEdges(il, a, b);
 
       // pull the surface if necessary (only in the front of the crystals)
       XYZPoint aprime = a;
       if (pulled(origin, normal_, a)) {
         b = aprime;
         hasbeenpulled = true;
       }
 
       // compute the intersection.
       // Check that the intersection is in the [a,b] segment  if HADSHOWER
       // if EMSHOWER the intersection is calculated as if the crystals were infinite
       XYZPoint xx = (EMSHOWER) ? intersect(plan_, a, b, dummyt, false) : intersect(plan_, a, b, dummyt, true);
 
       if (dummyt > 1)
         behindback = true;
       //      std::cout << " Intersect " << il << " " << a << " " << b << " " << plan_ << " " << xx << std::endl;
       // check that the intersection actually exists
       if (xx.mag2() != 0) {
         corners[il] = xx;
       }
     }
     //      std::cout << " ncorners " << corners.size() << std::endl;
     if (behindback && EMSHOWER)
       detailedShowerTail_ = true;
     // If the quad is completly defined. Store it !
     if (corners.size() == 4) {
       padsatdepth_[ic] = CrystalPad(ic, corners, trans, bfactor_, !central_);
       // Parameter to be tuned
       if (hasbeenpulled)
         padsatdepth_[ic].setSurvivalProbability(pulledPadProbability_);
       validPads_[ic] = true;
       ++nquads;
       // In principle, this should be done after the quads reorganization. But it would cost one more loop
       //  quadsatdepth_[ic].extrems(locxmin,locxmax,locymin,locymax);
       //  if(locxmin<xmin_) xmin_=locxmin;
       //  if(locymin<ymin_) ymin_=locymin;
       //  if(locxmax>xmax_) xmax_=locxmax;
       //  if(locymax>ymax_) ymax_=locymax;
     } else {
       padsatdepth_[ic] = CrystalPad();
       validPads_[ic] = false;
     }
   }
   //  std::cout << " Number of quads " << quadsatdepth_.size() << std::endl;
   if (doreorg_)
     reorganizePads();
   //  std::cout << "Finished to reorganize " << std::endl;
   npadsatdepth_ = nquads;
   //  std::cout << " prepareCellIDMap " << std::endl;
 
   // Resize the Quads to allow for some numerical inaccuracy
   // in the "inside" function
   for (unsigned ic = 0; ic < ncrystals_; ++ic) {
     if (!validPads_[ic])
       continue;
 
     if (EMSHOWER)
       padsatdepth_[ic].resetCorners();
 
     padsatdepth_[ic].extrems(locxmin, locxmax, locymin, locymax);
     if (locxmin < xmin_)
       xmin_ = locxmin;
     if (locymin < ymin_)
       ymin_ = locymin;
     if (locxmax > xmax_)
       xmax_ = locxmax;
     if (locymax > ymax_)
       ymax_ = locymax;
   }
 
   sizex_ = (xmax_ - xmin_) / nx_;
   sizey_ = (ymax_ - ymin_) / ny_;
 
   // Make sure that sizex_ and sizey_ are set before running prepareCellIDMap
   prepareCrystalNumberArray();
 
   //  std::cout << " Finished  prepareCellIDMap " << std::endl;
   ncrackpadsatdepth_ = crackpadsatdepth_.size();
 
   return true;
 }

◆ getSegments()

const std::vector<CaloSegment>& EcalHitMaker::getSegments ( ) const

inline

retrieve the segments (the path in the crystal crossed by the extrapolation of the track. Debugging only

Definition at line 95 of file EcalHitMaker.h.

References segments_.

Referenced by HcalHitMaker::setDepth().

95 { return segments_; };

EcalHitMaker::segments_

std::vector< CaloSegment > segments_

Definition: EcalHitMaker.h:276

◆ getX0back()

double EcalHitMaker::getX0back ( ) const

inline

Definition at line 104 of file EcalHitMaker.h.

References maxX0_.

Referenced by EMShower::compute().

104 { return maxX0_; }

EcalHitMaker::maxX0_

double maxX0_

Definition: EcalHitMaker.h:207

◆ hcalCellLine()

void EcalHitMaker::hcalCellLine ( std::vector< CaloPoint > & cp ) const

private

Definition at line 465 of file EcalHitMaker.cc.

References DeadROC_duringRun::dir, PVValHelper::eta, DetId::Hcal, FSimTrack::hcalEntrance(), Calorimeter::hcalProperties(), CaloHitMaker::myCalorimeter, myTrack_, FSimTrack::onHcal(), FSimTrack::onVFcal(), CalorimeterProperties::thickness(), Calorimetry_cff::thickness, RawParticle::Vect(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by cellLine().

                                                               {
   //  FSimEvent& mySimEvent = myEventMgr->simSignal();
   //  FSimTrack myTrack = mySimEvent.track(fsimtrack_);
   int onHcal = myTrack_->onHcal();
 
   if (onHcal <= 2 && onHcal > 0) {
     XYZPoint point1 = (myTrack_->hcalEntrance().vertex()).Vect();
 
     double eta = point1.eta();
     // HCAL thickness in cm (assuming that the particle is coming from 000)
     double thickness = myCalorimeter->hcalProperties(onHcal)->thickness(eta);
     cp.push_back(CaloPoint(DetId::Hcal, point1));
     XYZVector dir = myTrack_->hcalEntrance().Vect().Unit();
     XYZPoint point2 = point1 + dir * thickness;
 
     cp.push_back(CaloPoint(DetId::Hcal, point2));
   }
   int onVFcal = myTrack_->onVFcal();
   if (onVFcal == 2) {
     XYZPoint point1 = (myTrack_->vfcalEntrance().vertex()).Vect();
     double eta = point1.eta();
     // HCAL thickness in cm (assuming that the particle is coming from 000)
     double thickness = myCalorimeter->hcalProperties(3)->thickness(eta);
     cp.push_back(CaloPoint(DetId::Hcal, point1));
     XYZVector dir = myTrack_->vfcalEntrance().Vect().Unit();
     if (thickness > 0) {
       XYZPoint point2 = point1 + dir * thickness;
       cp.push_back(CaloPoint(DetId::Hcal, point2));
     }
   }
 }

◆ hcalTotalL0()

double EcalHitMaker::hcalTotalL0 ( ) const

inline

in the HCAL

Definition at line 88 of file EcalHitMaker.h.

References L0HCAL_.

Referenced by HDRShower::computeShower(), HDShower::HDShower(), and HFShower::HFShower().

88 { return L0HCAL_; }

EcalHitMaker::L0HCAL_

double L0HCAL_

Definition: EcalHitMaker.h:204

◆ hcalTotalX0()

double EcalHitMaker::hcalTotalX0 ( ) const

inline

in the HCAL

Definition at line 73 of file EcalHitMaker.h.

References X0HCAL_.

Referenced by EMShower::prepareSteps().

73 { return X0HCAL_; }

EcalHitMaker::X0HCAL_

double X0HCAL_

Definition: EcalHitMaker.h:199

◆ inside3D()

bool EcalHitMaker::inside3D	(	const std::vector< XYZPoint > &	corners,
		const XYZPoint &	p
	)		const

private

Definition at line 1379 of file EcalHitMaker.cc.

References corners, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by ecalCellLine().

                                                                                        {
   // corners and p are in the same plane
   // p is inside "corners" if the four crossproducts (corners[i]xcorners[i+1])
   // are in the same direction
 
   XYZVector crossproduct(0., 0., 0.), previouscrossproduct(0., 0., 0.);
 
   for (unsigned ip = 0; ip < 4; ++ip) {
     crossproduct = (corners[ip] - p).Cross(corners[(ip + 1) % 4] - p);
     if (ip == 0)
       previouscrossproduct = crossproduct;
     else if (crossproduct.Dot(previouscrossproduct) < 0.)
       return false;
   }
 
   return true;
 }

◆ prepareCrystalNumberArray()

void EcalHitMaker::prepareCrystalNumberArray ( )

private

Definition at line 1078 of file EcalHitMaker.cc.

References convertIntegerCoordinates(), d1, myCrystalNumberArray_, npadsatdepth_, padsatdepth_, validPads_, x, and y.

Referenced by getPads().

                                              {
   for (unsigned iq = 0; iq < npadsatdepth_; ++iq) {
     if (!validPads_[iq])
       continue;
     unsigned d1, d2;
     convertIntegerCoordinates(padsatdepth_[iq].center().x(), padsatdepth_[iq].center().y(), d1, d2);
     myCrystalNumberArray_[d1][d2] = iq;
   }
 }

◆ preshowerCellLine()

void EcalHitMaker::preshowerCellLine ( std::vector< CaloPoint > & cp ) const

private

Definition at line 420 of file EcalHitMaker.cc.

References DeadROC_duringRun::dir, DetId::Ecal, EcalPreshower, RawParticle::eta(), FSimTrack::layer1Entrance(), Calorimeter::layer1Properties(), FSimTrack::layer2Entrance(), Calorimeter::layer2Properties(), CaloHitMaker::myCalorimeter, myTrack_, FSimTrack::onLayer1(), FSimTrack::onLayer2(), PreshowerLayer1Properties::thickness(), PreshowerLayer2Properties::thickness(), RawParticle::Vect(), and RawParticle::vertex().

Referenced by cellLine().

                                                                    {
   //  FSimEvent& mySimEvent = myEventMgr->simSignal();
   //  FSimTrack myTrack = mySimEvent.track(fsimtrack_);
   //  std::cout << "FsimTrack " << fsimtrack_<< std::endl;
   //  std::cout << " On layer 1 " << myTrack.onLayer1() << std::endl;
   // std::cout << " preshowerCellLine " << std::endl;
   if (myTrack_->onLayer1()) {
     XYZPoint point1 = (myTrack_->layer1Entrance().vertex()).Vect();
     double phys_eta = myTrack_->layer1Entrance().eta();
     double cmthickness = myCalorimeter->layer1Properties(1)->thickness(phys_eta);
 
     if (cmthickness > 0) {
       XYZVector dir = myTrack_->layer1Entrance().Vect().Unit();
       XYZPoint point2 = point1 + dir * cmthickness;
 
       CaloPoint cp1(DetId::Ecal, EcalPreshower, 1, point1);
       CaloPoint cp2(DetId::Ecal, EcalPreshower, 1, point2);
       cp.push_back(cp1);
       cp.push_back(cp2);
     } else {
       //      std::cout << "Track on ECAL " << myTrack.EcalEntrance_().vertex()*0.1<< std::endl;
     }
   }
 
   //  std::cout << " On layer 2 " << myTrack.onLayer2() << std::endl;
   if (myTrack_->onLayer2()) {
     XYZPoint point1 = (myTrack_->layer2Entrance().vertex()).Vect();
     double phys_eta = myTrack_->layer2Entrance().eta();
     double cmthickness = myCalorimeter->layer2Properties(1)->thickness(phys_eta);
     if (cmthickness > 0) {
       XYZVector dir = myTrack_->layer2Entrance().Vect().Unit();
       XYZPoint point2 = point1 + dir * cmthickness;
 
       CaloPoint cp1(DetId::Ecal, EcalPreshower, 2, point1);
       CaloPoint cp2(DetId::Ecal, EcalPreshower, 2, point2);
 
       cp.push_back(cp1);
       cp.push_back(cp2);
     } else {
       //      std::cout << "Track on ECAL " << myTrack.EcalEntrance_().vertex()*0.1 << std::endl;
     }
   }
   //  std::cout << " Exit preshower CellLine " << std::endl;
 }

◆ ps1TotalL0()

double EcalHitMaker::ps1TotalL0 ( ) const

inline

total number of L0 in the PS (Layer1).

Definition at line 76 of file EcalHitMaker.h.

References L0PS1_.

76 { return L0PS1_; }

EcalHitMaker::L0PS1_

double L0PS1_

Definition: EcalHitMaker.h:200

◆ ps1TotalX0()

double EcalHitMaker::ps1TotalX0 ( ) const

inline

Definition at line 58 of file EcalHitMaker.h.

References X0PS1_.

Referenced by EMShower::prepareSteps().

58 { return X0PS1_; }

EcalHitMaker::X0PS1_

double X0PS1_

Definition: EcalHitMaker.h:194

◆ ps2eeTotalL0()

double EcalHitMaker::ps2eeTotalL0 ( ) const

inline

Definition at line 82 of file EcalHitMaker.h.

References L0PS2EE_.

82 { return L0PS2EE_; }

EcalHitMaker::L0PS2EE_

double L0PS2EE_

Definition: EcalHitMaker.h:202

◆ ps2eeTotalX0()

double EcalHitMaker::ps2eeTotalX0 ( ) const

inline

Definition at line 64 of file EcalHitMaker.h.

References X0PS2EE_.

Referenced by EMShower::prepareSteps().

64 { return X0PS2EE_; }

EcalHitMaker::X0PS2EE_

double X0PS2EE_

Definition: EcalHitMaker.h:196

◆ ps2TotalL0()

double EcalHitMaker::ps2TotalL0 ( ) const

inline

total number of L0 in the PS (Layer2).

Definition at line 79 of file EcalHitMaker.h.

References L0PS2_.

79 { return L0PS2_; }

EcalHitMaker::L0PS2_

double L0PS2_

Definition: EcalHitMaker.h:201

◆ ps2TotalX0()

double EcalHitMaker::ps2TotalX0 ( ) const

inline

total number of X0 in the PS (Layer2).

Definition at line 61 of file EcalHitMaker.h.

References X0PS2_.

Referenced by EMShower::prepareSteps().

61 { return X0PS2_; }

EcalHitMaker::X0PS2_

double X0PS2_

Definition: EcalHitMaker.h:195

◆ pulled()

bool EcalHitMaker::pulled	(	const XYZPoint &	origin,
		const XYZNormal &	normal,
		XYZPoint &	fPoint
	)		const

private

Definition at line 1067 of file EcalHitMaker.cc.

References normal.

Referenced by getPads().

                                                                                                  {
   // check if fPoint is behind the origin
   double dotproduct = normal.Dot(fPoint - origin);
   if (dotproduct <= 0.)
     return false;
 
   //norm of normal is 1
   fPoint -= (1 + dotproduct) * normal;
   return true;
 }

◆ reorganizePads()

void EcalHitMaker::reorganizePads ( )

private

Definition at line 1123 of file EcalHitMaker.cc.

References crackpadsatdepth_, cracksPads(), etasize_, gapsLifting(), ncrackpadsatdepth_, ncrystals_, phisize_, regionOfInterest_, CaloDirectionOperations::Side(), and validPads_.

Referenced by getPads().

                                   {
   // Some cleaning first
   //  std::cout << " Starting reorganize " << std::endl;
   crackpadsatdepth_.clear();
   crackpadsatdepth_.reserve(etasize_ * phisize_);
   ncrackpadsatdepth_ = 0;
   std::vector<neighbour> gaps;
   std::vector<std::vector<neighbour> > cracks;
   //  cracks.clear();
   cracks.resize(ncrystals_);
 
   for (unsigned iq = 0; iq < ncrystals_; ++iq) {
     if (!validPads_[iq])
       continue;
 
     gaps.clear();
     //   std::cout << padsatdepth_[iq] << std::endl;
     //check all the directions
     for (unsigned iside = 0; iside < 4; ++iside) {
       //      std::cout << " To be glued " << iside << " " << regionOfInterest_[iq].crystalNeighbour(iside).toBeGlued() << std::endl;
       CaloDirection thisside = CaloDirectionOperations::Side(iside);
       if (regionOfInterest_[iq].crystalNeighbour(iside).toBeGlued()) {
         // look for the neighbour and check that it exists
         int neighbourstatus = regionOfInterest_[iq].crystalNeighbour(iside).status();
         if (neighbourstatus < 0)
           continue;
 
         unsigned neighbourNumber = regionOfInterest_[iq].crystalNeighbour(iside).number();
         if (!validPads_[neighbourNumber])
           continue;
         // there is a crack between
         if (neighbourstatus == 1) {
           //          std::cout << " 1 Crack : " << thisside << " " << cellids_[iq]<< " " << cellids_[neighbourNumber] << std::endl;
           cracks[iq].push_back(neighbour(thisside, neighbourNumber));
         }  // else it is a gap
         else {
           gaps.push_back(neighbour(thisside, neighbourNumber));
         }
       }
     }
     // Now lift the gaps
     gapsLifting(gaps, iq);
   }
 
   unsigned ncracks = cracks.size();
   //  std::cout << " Cracks treatment : " << cracks.size() << std::endl;
   for (unsigned icrack = 0; icrack < ncracks; ++icrack) {
     //      std::cout << " Crack number " << crackiter->first << std::endl;
     cracksPads(cracks[icrack], icrack);
   }
 }

◆ setCrackPadSurvivalProbability()

void EcalHitMaker::setCrackPadSurvivalProbability ( double val )

inline

Definition at line 129 of file EcalHitMaker.h.

References crackPadProbability_, and heppy_batch::val.

Referenced by CalorimetryManager::EMShowerSimulation().

129 { crackPadProbability_ = val; };

EcalHitMaker::crackPadProbability_

double crackPadProbability_

Definition: EcalHitMaker.h:266

heppy_batch.val

val

Definition: heppy_batch.py:351

◆ setPreshowerPresent()

void EcalHitMaker::setPreshowerPresent ( bool ps )

inline

Definition at line 132 of file EcalHitMaker.h.

References simulatePreshower_.

Referenced by CalorimetryManager::EMShowerSimulation().

132 { simulatePreshower_ = ps; };

EcalHitMaker::simulatePreshower_

bool simulatePreshower_

Definition: EcalHitMaker.h:256

◆ setPulledPadSurvivalProbability()

void EcalHitMaker::setPulledPadSurvivalProbability ( double val )

inline

Definition at line 127 of file EcalHitMaker.h.

References pulledPadProbability_, and heppy_batch::val.

Referenced by CalorimetryManager::EMShowerSimulation().

127 { pulledPadProbability_ = val; };

EcalHitMaker::pulledPadProbability_

double pulledPadProbability_

Definition: EcalHitMaker.h:264

heppy_batch.val

val

Definition: heppy_batch.py:351

◆ setRadiusFactor()

void EcalHitMaker::setRadiusFactor ( double r )

inline

Definition at line 125 of file EcalHitMaker.h.

References alignCSCRings::r, and radiusCorrectionFactor_.

Referenced by CalorimetryManager::EMShowerSimulation().

125 { radiusCorrectionFactor_ = r; }

EcalHitMaker::radiusCorrectionFactor_

double radiusCorrectionFactor_

Definition: EcalHitMaker.h:246

alignCSCRings.r

r

Definition: alignCSCRings.py:93

◆ setSpotEnergy()

void EcalHitMaker::setSpotEnergy ( double e )

inlineoverridevirtual

Implements CaloHitMaker.

Definition at line 112 of file EcalHitMaker.h.

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

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

112 { spotEnergy = e; }

MillePedeFileConverter_cfg.e

e

Definition: MillePedeFileConverter_cfg.py:37

CaloHitMaker::spotEnergy

double spotEnergy

Definition: CaloHitMaker.h:48

◆ setTrackParameters()

void EcalHitMaker::setTrackParameters	(	const XYZNormal &	normal,
		double	X0depthoffset,
		const FSimTrack &	theTrack
	)

Definition at line 301 of file EcalHitMaker.cc.

References buildSegments(), cellLine(), DeadROC_duringRun::dir, DetId::Ecal, ecalFirstSegment_, ecalTotalX0(), CaloHitMaker::EMSHOWER, intersections_, maxX0_, myTrack_, ncrystals_, normal, normal_, onEcal_, regionOfInterest_, segments_, CaloDirectionOperations::Side(), and X0depthoffset_.

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

                                                                                                               {
   //  myHistos->debug("setTrackParameters");
   //  std::cout << " Track " << theTrack << std::endl;
   intersections_.clear();
   // This is certainly enough
   intersections_.reserve(50);
   myTrack_ = &theTrack;
   normal_ = normal.Unit();
   X0depthoffset_ = X0depthoffset;
   cellLine(intersections_);
   buildSegments(intersections_);
   //  std::cout << " Segments " << segments_.size() << std::endl;
   //  for(unsigned ii=0; ii<segments_.size() ; ++ii)
   //    {
   //      std::cout << segments_[ii] << std::endl;
   //    }
 
   // This is only needed in case of electromagnetic showers
   if (EMSHOWER && onEcal_ && ecalTotalX0() > 0.) {
     //      std::cout << "Total X0  " << ecalTotalX0() << std::endl;
     for (unsigned ic = 0; ic < ncrystals_; ++ic) {
       for (unsigned idir = 0; idir < 4; ++idir) {
         XYZVector norm = regionOfInterest_[ic].exitingNormal(CaloDirectionOperations::Side(idir));
         regionOfInterest_[ic].crystalNeighbour(idir).setToBeGlued((norm.Dot(normal_) < 0.));
       }
       // Now calculate the distance in X0 of the back sides of the crystals
       // (only for EM showers)
       if (EMSHOWER) {
         XYZVector dir = regionOfInterest_[ic].getBackCenter() - segments_[ecalFirstSegment_].entrance();
         double dist = dir.Dot(normal_);
         double absciss = dist + segments_[ecalFirstSegment_].sEntrance();
         std::vector<CaloSegment>::const_iterator segiterator;
         // First identify the correct segment
         //      std::cout << " Crystal " << ic  << regionOfInterest_[ic].getBackCenter() ;
         //      std::cout << " Entrance : " << segments_[ecalFirstSegment_].entrance()<< std::endl;
         //      std::cout << " Looking for the segment " << dist << std::endl;
         segiterator = find_if(segments_.begin(), segments_.end(), CaloSegment::inSegment(absciss));
         //      std::cout << " Done " << std::endl;
         if (segiterator == segments_.end()) {
           // in this case, we won't have any problem. No need to
           // calculate the real depth.
           regionOfInterest_[ic].setX0Back(9999);
         } else {
           DetId::Detector det(segiterator->whichDetector());
           if (det != DetId::Ecal) {
             regionOfInterest_[ic].setX0Back(9999);
           } else {
             double x0 = segiterator->x0FromCm(dist);
             if (x0 < maxX0_)
               maxX0_ = x0;
             regionOfInterest_[ic].setX0Back(x0);
           }
         }
         //  myHistos->fill("h4000",ecalentrance_.eta(), regionOfInterest_[ic].getX0Back());
         //}
         //      else
         //{
         //   // in this case, we won't have any problem. No need to
         //  // calculate the real depth.
         //  regionOfInterest_[ic].setX0Back(9999);
         //}
       }  //EMSHOWER
     }    // ndir
     //      myHistos->fill("h6000",segments_[ecalFirstSegment_].entrance().eta(),maxX0_);
   }
   //  std::cout << "Leaving setTrackParameters" << std::endl
 }

◆ totalL0()

double EcalHitMaker::totalL0 ( ) const

inline

Number of interaction length "seen" by the track.

Definition at line 52 of file EcalHitMaker.h.

References totalL0_.

52 { return totalL0_; };

EcalHitMaker::totalL0_

double totalL0_

Definition: EcalHitMaker.h:192

◆ totalX0()

double EcalHitMaker::totalX0 ( ) const

inline

Definition at line 49 of file EcalHitMaker.h.

References totalX0_.

Referenced by EMShower::prepareSteps().

49 { return totalX0_; };

EcalHitMaker::totalX0_

double totalX0_

Definition: EcalHitMaker.h:191

◆ unbalancedDirection()

bool EcalHitMaker::unbalancedDirection	(	const std::vector< neighbour > &	dirs,
		unsigned &	unb,
		unsigned &	dir1,
		unsigned &	dir2
	)

private

Definition at line 1196 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), L1TDiffHarvesting_cfi::dir1, L1TDiffHarvesting_cfi::dir2, heppy_check::dirs, dqmdumpme::first, NONE, and createJobs::tmp.

Referenced by gapsLifting().

                                                        {
   if (dirs.size() == 1)
     return false;
   if (dirs.size() % 2 == 0)
     return false;
   CaloDirection tmp;
   tmp = CaloDirectionOperations::add2d(dirs[0].first, dirs[1].first);
   if (tmp == NONE) {
     unb = 2;
     dir1 = 0;
     dir2 = 1;
     return true;
   }
   tmp = CaloDirectionOperations::add2d(dirs[0].first, dirs[2].first);
   if (tmp == NONE) {
     unb = 1;
     dir1 = 0;
     dir2 = 2;
     return true;
   }
   unb = 0;
   dir1 = 1;
   dir2 = 2;
   return true;
 }

◆ x0DepthOffset()

double EcalHitMaker::x0DepthOffset ( ) const

inline

get the offset (e.g the number of X0 after which the shower starts)

Definition at line 55 of file EcalHitMaker.h.

References X0depthoffset_.

Referenced by EMShower::prepareSteps().

55 { return X0depthoffset_; }

EcalHitMaker::X0depthoffset_

double X0depthoffset_

Definition: EcalHitMaker.h:193

Member Data Documentation

◆ bfactor_

double EcalHitMaker::bfactor_

private

Definition at line 254 of file EcalHitMaker.h.

Referenced by configureGeometry(), EcalHitMaker(), and getPads().

◆ CellsWindow_

std::vector<DetId> EcalHitMaker::CellsWindow_

private

Definition at line 224 of file EcalHitMaker.h.

Referenced by buildGeometry(), and EcalHitMaker().

◆ central_

int EcalHitMaker::central_

private

Definition at line 215 of file EcalHitMaker.h.

Referenced by cellLine(), EcalHitMaker(), and getPads().

◆ configuredGeometry_

bool EcalHitMaker::configuredGeometry_

private

Definition at line 218 of file EcalHitMaker.h.

Referenced by buildGeometry(), configureGeometry(), and getPads().

◆ corners

std::vector<XYZPoint> EcalHitMaker::corners

private

Definition at line 288 of file EcalHitMaker.h.

Referenced by ecalCellLine(), EcalHitMaker(), getPads(), and inside3D().

◆ crackPadProbability_

double EcalHitMaker::crackPadProbability_

private

Definition at line 266 of file EcalHitMaker.h.

Referenced by cracksPads(), and setCrackPadSurvivalProbability().

◆ crackpadsatdepth_

std::vector<CrystalPad> EcalHitMaker::crackpadsatdepth_

private

Definition at line 282 of file EcalHitMaker.h.

Referenced by cracksPads(), fastInsideCell(), getPads(), and reorganizePads().

◆ currentdepth_

double EcalHitMaker::currentdepth_

private

Definition at line 252 of file EcalHitMaker.h.

Referenced by getPads().

◆ detailedShowerTail_

bool EcalHitMaker::detailedShowerTail_

private

Definition at line 250 of file EcalHitMaker.h.

Referenced by getPads().

◆ DetIdMap_

std::map<DetId, unsigned> EcalHitMaker::DetIdMap_

private

Definition at line 231 of file EcalHitMaker.h.

Referenced by buildGeometry(), and configureGeometry().

◆ doreorg_

bool EcalHitMaker::doreorg_

private

Definition at line 278 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and getPads().

◆ EcalEntrance_

XYZPoint EcalHitMaker::EcalEntrance_

private

Definition at line 213 of file EcalHitMaker.h.

Referenced by buildSegments(), cellLine(), configureGeometry(), ecalEntrance(), and EcalHitMaker().

◆ ecalFirstSegment_

int EcalHitMaker::ecalFirstSegment_

private

Definition at line 236 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and setTrackParameters().

◆ etasize_

unsigned EcalHitMaker::etasize_

private

Definition at line 239 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and reorganizePads().

◆ hitmaphasbeencalculated_

bool EcalHitMaker::hitmaphasbeencalculated_

private

Definition at line 284 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and getHits().

◆ hits_

std::vector<float> EcalHitMaker::hits_

private

Definition at line 226 of file EcalHitMaker.h.

Referenced by addHit(), addHitDepth(), buildGeometry(), and getHits().

◆ intersections_

std::vector<CaloPoint> EcalHitMaker::intersections_

private

Definition at line 274 of file EcalHitMaker.h.

Referenced by setTrackParameters().

◆ L0ECAL_

double EcalHitMaker::L0ECAL_

private

Definition at line 203 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ecalTotalL0().

◆ L0EHGAP_

double EcalHitMaker::L0EHGAP_

private

Definition at line 205 of file EcalHitMaker.h.

Referenced by buildSegments(), ecalHcalGapTotalL0(), and EcalHitMaker().

◆ L0HCAL_

double EcalHitMaker::L0HCAL_

private

Definition at line 204 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and hcalTotalL0().

◆ L0PS1_

double EcalHitMaker::L0PS1_

private

Definition at line 200 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ps1TotalL0().

◆ L0PS2_

double EcalHitMaker::L0PS2_

private

Definition at line 201 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ps2TotalL0().

◆ L0PS2EE_

double EcalHitMaker::L0PS2EE_

private

Definition at line 202 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ps2eeTotalL0().

◆ maxX0_

double EcalHitMaker::maxX0_

private

Definition at line 207 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), getX0back(), and setTrackParameters().

◆ mycorners

std::vector<CLHEP::Hep2Vector> EcalHitMaker::mycorners

private

Definition at line 287 of file EcalHitMaker.h.

Referenced by cracksPads(), and EcalHitMaker().

◆ myCrystalNumberArray_

std::vector<std::vector<unsigned> > EcalHitMaker::myCrystalNumberArray_

private

Definition at line 188 of file EcalHitMaker.h.

Referenced by buildGeometry(), fastInsideCell(), and prepareCrystalNumberArray().

◆ myCrystalWindowMap_

CrystalWindowMap* EcalHitMaker::myCrystalWindowMap_

private

Definition at line 233 of file EcalHitMaker.h.

Referenced by buildGeometry(), EcalHitMaker(), fastInsideCell(), and ~EcalHitMaker().

◆ myTrack_

const FSimTrack* EcalHitMaker::myTrack_

private

Definition at line 271 of file EcalHitMaker.h.

Referenced by buildSegments(), cellLine(), getFSimTrack(), getPads(), hcalCellLine(), preshowerCellLine(), and setTrackParameters().

◆ ncrackpadsatdepth_

unsigned EcalHitMaker::ncrackpadsatdepth_

private

Definition at line 259 of file EcalHitMaker.h.

Referenced by fastInsideCell(), getPads(), and reorganizePads().

◆ ncrystals_

unsigned EcalHitMaker::ncrystals_

private

Definition at line 219 of file EcalHitMaker.h.

Referenced by buildGeometry(), configureGeometry(), ecalCellLine(), EcalHitMaker(), getHits(), getPads(), reorganizePads(), and setTrackParameters().

◆ normal_

XYZNormal EcalHitMaker::normal_

private

Definition at line 214 of file EcalHitMaker.h.

Referenced by cellLine(), ecalCellLine(), getPads(), and setTrackParameters().

◆ npadsatdepth_

unsigned EcalHitMaker::npadsatdepth_

private

Definition at line 260 of file EcalHitMaker.h.

Referenced by getPads(), and prepareCrystalNumberArray().

◆ nx_

unsigned EcalHitMaker::nx_

private

Definition at line 221 of file EcalHitMaker.h.

Referenced by buildGeometry(), fastInsideCell(), and getPads().

◆ ny_

unsigned EcalHitMaker::ny_

private

Definition at line 221 of file EcalHitMaker.h.

Referenced by buildGeometry(), fastInsideCell(), and getPads().

◆ onEcal_

int EcalHitMaker::onEcal_

private

Definition at line 216 of file EcalHitMaker.h.

Referenced by cellLine(), getHits(), and setTrackParameters().

◆ outsideWindowEnergy_

double EcalHitMaker::outsideWindowEnergy_

private

Definition at line 209 of file EcalHitMaker.h.

Referenced by addHit(), addHitDepth(), and EcalHitMaker().

◆ padsatdepth_

std::vector<CrystalPad> EcalHitMaker::padsatdepth_

private

Definition at line 281 of file EcalHitMaker.h.

Referenced by buildGeometry(), correspondingEdge(), cracksPads(), diagonalEdge(), fastInsideCell(), gapsLifting(), getPads(), and prepareCrystalNumberArray().

◆ phisize_

unsigned EcalHitMaker::phisize_

private

Definition at line 240 of file EcalHitMaker.h.

Referenced by buildGeometry(), EcalHitMaker(), and reorganizePads().

◆ pivot_

Crystal EcalHitMaker::pivot_

private

Definition at line 212 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and getPads().

◆ plan_

Plane3D EcalHitMaker::plan_

private

Definition at line 261 of file EcalHitMaker.h.

Referenced by getPads().

◆ pulledPadProbability_

double EcalHitMaker::pulledPadProbability_

private

Definition at line 264 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), getPads(), and setPulledPadSurvivalProbability().

◆ radiusCorrectionFactor_

double EcalHitMaker::radiusCorrectionFactor_

private

Definition at line 246 of file EcalHitMaker.h.

Referenced by getPads(), and setRadiusFactor().

◆ radiusFactor_

double EcalHitMaker::radiusFactor_

private

Definition at line 248 of file EcalHitMaker.h.

Referenced by addHit(), addHitDepth(), and getPads().

◆ random

const RandomEngineAndDistribution* EcalHitMaker::random

private

Definition at line 290 of file EcalHitMaker.h.

Referenced by addHit().

◆ rearleakage_

double EcalHitMaker::rearleakage_

private

Definition at line 208 of file EcalHitMaker.h.

Referenced by addHitDepth(), and EcalHitMaker().

◆ regionOfInterest_

std::vector<Crystal> EcalHitMaker::regionOfInterest_

private

Definition at line 225 of file EcalHitMaker.h.

Referenced by addHitDepth(), buildGeometry(), configureGeometry(), diagonalEdge(), ecalCellLine(), getCrystals(), getHits(), getPads(), reorganizePads(), and setTrackParameters().

◆ segments_

std::vector<CaloSegment> EcalHitMaker::segments_

private

Definition at line 276 of file EcalHitMaker.h.

Referenced by buildSegments(), getPads(), getSegments(), and setTrackParameters().

◆ simulatePreshower_

bool EcalHitMaker::simulatePreshower_

private

Definition at line 256 of file EcalHitMaker.h.

Referenced by cellLine(), EcalHitMaker(), and setPreshowerPresent().

◆ sizex_

double EcalHitMaker::sizex_

private

Definition at line 268 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

◆ sizey_

double EcalHitMaker::sizey_

private

Definition at line 268 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

◆ totalL0_

double EcalHitMaker::totalL0_

private

Definition at line 192 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and totalL0().

◆ totalX0_

double EcalHitMaker::totalX0_

private

Definition at line 191 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and totalX0().

◆ truncatedGrid_

bool EcalHitMaker::truncatedGrid_

private

Definition at line 242 of file EcalHitMaker.h.

Referenced by EcalHitMaker().

◆ validPads_

std::vector<bool> EcalHitMaker::validPads_

private

Definition at line 228 of file EcalHitMaker.h.

Referenced by buildGeometry(), diagonalEdge(), fastInsideCell(), getPads(), prepareCrystalNumberArray(), and reorganizePads().

◆ X0depthoffset_

double EcalHitMaker::X0depthoffset_

private

Definition at line 193 of file EcalHitMaker.h.

Referenced by addHitDepth(), EcalHitMaker(), getPads(), setTrackParameters(), and x0DepthOffset().

◆ X0ECAL_

double EcalHitMaker::X0ECAL_

private

Definition at line 197 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ecalTotalX0().

◆ X0EHGAP_

double EcalHitMaker::X0EHGAP_

private

Definition at line 198 of file EcalHitMaker.h.

Referenced by buildSegments(), ecalHcalGapTotalX0(), and EcalHitMaker().

◆ X0HCAL_

double EcalHitMaker::X0HCAL_

private

Definition at line 199 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and hcalTotalX0().

◆ X0PS1_

double EcalHitMaker::X0PS1_

private

Definition at line 194 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ps1TotalX0().

◆ X0PS2_

double EcalHitMaker::X0PS2_

private

Definition at line 195 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ps2TotalX0().

◆ X0PS2EE_

double EcalHitMaker::X0PS2EE_

private

Definition at line 196 of file EcalHitMaker.h.

Referenced by buildSegments(), EcalHitMaker(), and ps2eeTotalX0().

◆ xmax_

double EcalHitMaker::xmax_

private

Definition at line 222 of file EcalHitMaker.h.

Referenced by getPads().

◆ xmin_

double EcalHitMaker::xmin_

private

Definition at line 222 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

◆ ymax_

double EcalHitMaker::ymax_

private

Definition at line 222 of file EcalHitMaker.h.

Referenced by getPads().

◆ ymin_

double EcalHitMaker::ymin_

private

Definition at line 222 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ neighbour

◆ Plane3D

◆ XYZNormal

◆ XYZPoint

◆ XYZVector

Constructor & Destructor Documentation

◆ EcalHitMaker()

◆ ~EcalHitMaker()

Member Function Documentation

◆ addHit()

◆ addHitDepth()

◆ buildGeometry()

◆ buildSegments()

◆ cellLine()

◆ configureGeometry()

◆ convertIntegerCoordinates()

◆ correspondingEdge()

◆ cracksPads()

◆ diagonalEdge()

◆ ecalCellLine()

◆ ecalEntrance()

◆ ecalHcalGapTotalL0()

◆ ecalHcalGapTotalX0()

◆ ecalTotalL0()

◆ ecalTotalX0()

◆ fastInsideCell()

◆ gapsLifting()

◆ getCrystals()

◆ getFSimTrack()

◆ getHits()

◆ getPads()

◆ getSegments()

◆ getX0back()

◆ hcalCellLine()

◆ hcalTotalL0()

◆ hcalTotalX0()

◆ inside3D()

◆ prepareCrystalNumberArray()

◆ preshowerCellLine()

◆ ps1TotalL0()

◆ ps1TotalX0()

◆ ps2eeTotalL0()

◆ ps2eeTotalX0()

◆ ps2TotalL0()

◆ ps2TotalX0()

◆ pulled()

◆ reorganizePads()

◆ setCrackPadSurvivalProbability()

◆ setPreshowerPresent()

◆ setPulledPadSurvivalProbability()

◆ setRadiusFactor()

◆ setSpotEnergy()

◆ setTrackParameters()

◆ totalL0()

◆ totalX0()

◆ unbalancedDirection()

◆ x0DepthOffset()

Member Data Documentation

◆ bfactor_

◆ CellsWindow_

◆ central_

◆ configuredGeometry_

◆ corners

◆ crackPadProbability_

◆ crackpadsatdepth_

◆ currentdepth_

◆ detailedShowerTail_

◆ DetIdMap_

◆ doreorg_

◆ EcalEntrance_

◆ ecalFirstSegment_

◆ etasize_