#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
	not been done. More...

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

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

private

Definition at line 178 of file EcalHitMaker.h.

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

Definition at line 31 of file EcalHitMaker.h.

typedef math::XYZVector EcalHitMaker::XYZNormal

Definition at line 30 of file EcalHitMaker.h.

typedef math::XYZVector EcalHitMaker::XYZPoint

Definition at line 29 of file EcalHitMaker.h.

typedef math::XYZVector EcalHitMaker::XYZVector

Definition at line 28 of file EcalHitMaker.h.

Constructor & Destructor Documentation

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

override

Definition at line 126 of file EcalHitMaker.cc.

References myCrystalWindowMap_.

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

Member Function Documentation

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

overridevirtual

Implements CaloHitMaker.

Definition at line 195 of file EcalHitMaker.cc.

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

Referenced by HFShower::compute(), HDShower::compute(), EMShower::compute(), getX0back(), 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;
 }

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

Definition at line 135 of file EcalHitMaker.cc.

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

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

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

void EcalHitMaker::buildGeometry ( )

private

Definition at line 880 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(), and getCrystals().

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

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

private

Definition at line 664 of file EcalHitMaker.cc.

References CaloGeometryHelper::borderCrossing(), 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 getCrystals(), and 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.;
   
   unsigned ncrossedxtals = 0;
   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();
          ++ncrossedxtals;
          ++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;
     
 }

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

private

Definition at line 411 of file EcalHitMaker.cc.

References central_, gather_cfg::cout, 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(), FSimTrack::vertex(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by getCrystals(), and 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;
 //    }
 }

void EcalHitMaker::configureGeometry ( )

private

Definition at line 1116 of file EcalHitMaker.cc.

References bfactor_, CaloGeometryHelper::borderCrossing(), configuredGeometry_, DetIdMap_, 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 getCrystals(), and 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.;
 }

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

private

Definition at line 1208 of file EcalHitMaker.cc.

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

Referenced by fastInsideCell(), getCrystals(), 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;
 }

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

private

Definition at line 1305 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), TCMET_cfi::corner, dir, 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);  
 }

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

private

Definition at line 1493 of file EcalHitMaker.cc.

References crackPadProbability_, crackpadsatdepth_, EAST, CrystalPad::edge(), plotBeamSpotDB::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;
 }

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

private

Definition at line 1313 of file EcalHitMaker.cc.

References CaloDirectionOperations::neighbourDirection(), CaloDirectionOperations::oppositeSide(), padsatdepth_, 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;
 }

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

private

Definition at line 571 of file EcalHitMaker.cc.

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

Referenced by cellLine(), and getCrystals().

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

const XYZPoint& EcalHitMaker::ecalEntrance ( ) const

inline

used in FamosHcalHitMaker

Definition at line 128 of file EcalHitMaker.h.

References EcalEntrance_.

Referenced by HcalHitMaker::HcalHitMaker().

128 {return EcalEntrance_;};

EcalHitMaker::EcalEntrance_

XYZPoint EcalEntrance_

Definition: EcalHitMaker.h:225

double EcalHitMaker::ecalHcalGapTotalL0 ( ) const

inline

ECAL-HCAL transition.

Definition at line 95 of file EcalHitMaker.h.

References L0EHGAP_.

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

95 { return L0EHGAP_;}

EcalHitMaker::L0EHGAP_

double L0EHGAP_

Definition: EcalHitMaker.h:217

double EcalHitMaker::ecalHcalGapTotalX0 ( ) const

inline

ECAL-HCAL transition.

Definition at line 74 of file EcalHitMaker.h.

References X0EHGAP_.

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

74 { return X0EHGAP_;}

EcalHitMaker::X0EHGAP_

double X0EHGAP_

Definition: EcalHitMaker.h:210

double EcalHitMaker::ecalTotalL0 ( ) const

inline

in the ECAL

Definition at line 89 of file EcalHitMaker.h.

References L0ECAL_.

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

89 {return L0ECAL_;}

EcalHitMaker::L0ECAL_

double L0ECAL_

Definition: EcalHitMaker.h:215

double EcalHitMaker::ecalTotalX0 ( ) const

inline

in the ECAL

Definition at line 71 of file EcalHitMaker.h.

References X0ECAL_.

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

71 {return X0ECAL_;}

EcalHitMaker::X0ECAL_

double X0ECAL_

Definition: EcalHitMaker.h:209

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

Definition at line 244 of file EcalHitMaker.cc.

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

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

 {
 
   //  debug = true;
   bool found=false;
   unsigned niter=0;
   // 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;
       ++niter;
     }
           }
   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;
       ++niter;
     }
   //  myHistos->fill("h1002",niter);
   if(!found&&debug) std::cout << " Not found in the cracks " << std::endl;
   return (found) ? crackpadsatdepth_[iquadinside].getNumber(): 9999;
 }

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

private

Definition at line 1357 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), constexpr, correspondingEdge(), gather_cfg::cout, debug, diagonalEdge(), GetRecoTauVFromDQM_MC_cff::dir2, EAST, CrystalPad::edge(), plotBeamSpotDB::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;
       }   
 }

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

inline

for debugging

Definition at line 140 of file EcalHitMaker.h.

References buildGeometry(), buildSegments(), cellLine(), configureGeometry(), convertIntegerCoordinates(), SimDataFormats::CaloAnalysis::cp, ecalCellLine(), hcalCellLine(), prepareCrystalNumberArray(), preshowerCellLine(), pulled(), regionOfInterest_, reorganizePads(), x, and y.

140 {return regionOfInterest_;}

EcalHitMaker::regionOfInterest_

std::vector< Crystal > regionOfInterest_

Definition: EcalHitMaker.h:237

const FSimTrack* EcalHitMaker::getFSimTrack ( ) const

inline

To retrieve the track.

Definition at line 125 of file EcalHitMaker.h.

References myTrack_.

Referenced by HcalHitMaker::HcalHitMaker(), and HcalHitMaker::setDepth().

125 {return myTrack_;}

EcalHitMaker::myTrack_

const FSimTrack * myTrack_

Definition: EcalHitMaker.h:284

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

overridevirtual

not been done.

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

Implements CaloHitMaker.

Definition at line 1217 of file EcalHitMaker.cc.

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

Referenced by CalorimetryManager::EMShowerSimulation(), and setSpotEnergy().

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

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 917 of file EcalHitMaker.cc.

References a, b, bfactor_, central_, configuredGeometry_, configureGeometry(), corners, gather_cfg::cout, crackpadsatdepth_, currentdepth_, egammaForCoreTracking_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(), HDRShower::computeShower(), and getSegments().

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

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

References egammaForCoreTracking_cff::depth, getPads(), and segments_.

Referenced by HcalHitMaker::setDepth().

99 {return segments_;};

EcalHitMaker::segments_

std::vector< CaloSegment > segments_

Definition: EcalHitMaker.h:289

double EcalHitMaker::getX0back ( ) const

inline

Definition at line 109 of file EcalHitMaker.h.

References addHit(), addHitDepth(), debug, fastInsideCell(), maxX0_, phi, point, and alignCSCRings::r.

Referenced by EMShower::compute().

109 {return maxX0_;}

EcalHitMaker::maxX0_

double maxX0_

Definition: EcalHitMaker.h:219

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

private

Definition at line 532 of file EcalHitMaker.cc.

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

Referenced by cellLine(), and getCrystals().

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

double EcalHitMaker::hcalTotalL0 ( ) const

inline

in the HCAL

Definition at line 92 of file EcalHitMaker.h.

References L0HCAL_.

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

92 {return L0HCAL_;}

EcalHitMaker::L0HCAL_

double L0HCAL_

Definition: EcalHitMaker.h:216

double EcalHitMaker::hcalTotalX0 ( ) const

inline

in the HCAL

Definition at line 77 of file EcalHitMaker.h.

References X0HCAL_.

Referenced by EMShower::prepareSteps().

77 {return X0HCAL_;}

EcalHitMaker::X0HCAL_

double X0HCAL_

Definition: EcalHitMaker.h:211

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

private

Definition at line 1549 of file EcalHitMaker.cc.

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

void EcalHitMaker::prepareCrystalNumberArray ( )

private

Definition at line 1197 of file EcalHitMaker.cc.

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

Referenced by getCrystals(), and 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;
     }
 }

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

private

Definition at line 474 of file EcalHitMaker.cc.

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

Referenced by cellLine(), and getCrystals().

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

double EcalHitMaker::ps1TotalL0 ( ) const

inline

total number of L0 in the PS (Layer1).

Definition at line 80 of file EcalHitMaker.h.

References L0PS1_.

80 {return L0PS1_;}

EcalHitMaker::L0PS1_

double L0PS1_

Definition: EcalHitMaker.h:212

double EcalHitMaker::ps1TotalX0 ( ) const

inline

Definition at line 62 of file EcalHitMaker.h.

References X0PS1_.

Referenced by EMShower::prepareSteps().

62 {return X0PS1_;}

EcalHitMaker::X0PS1_

double X0PS1_

Definition: EcalHitMaker.h:206

double EcalHitMaker::ps2eeTotalL0 ( ) const

inline

Definition at line 86 of file EcalHitMaker.h.

References L0PS2EE_.

86 {return L0PS2EE_;}

EcalHitMaker::L0PS2EE_

double L0PS2EE_

Definition: EcalHitMaker.h:214

double EcalHitMaker::ps2eeTotalX0 ( ) const

inline

Definition at line 68 of file EcalHitMaker.h.

References X0PS2EE_.

Referenced by EMShower::prepareSteps().

68 {return X0PS2EE_;}

EcalHitMaker::X0PS2EE_

double X0PS2EE_

Definition: EcalHitMaker.h:208

double EcalHitMaker::ps2TotalL0 ( ) const

inline

total number of L0 in the PS (Layer2).

Definition at line 83 of file EcalHitMaker.h.

References L0PS2_.

83 {return L0PS2_;}

EcalHitMaker::L0PS2_

double L0PS2_

Definition: EcalHitMaker.h:213

double EcalHitMaker::ps2TotalX0 ( ) const

inline

total number of X0 in the PS (Layer2).

Definition at line 65 of file EcalHitMaker.h.

References X0PS2_.

Referenced by EMShower::prepareSteps().

65 {return X0PS2_;}

EcalHitMaker::X0PS2_

double X0PS2_

Definition: EcalHitMaker.h:207

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

private

Definition at line 1182 of file EcalHitMaker.cc.

Referenced by getCrystals(), and 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;
 }

void EcalHitMaker::reorganizePads ( )

private

Definition at line 1243 of file EcalHitMaker.cc.

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

Referenced by getCrystals(), and 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);
     }
 
 }

void EcalHitMaker::setCrackPadSurvivalProbability ( double val )

inline

Definition at line 134 of file EcalHitMaker.h.

References crackPadProbability_, and heppy_batch::val.

Referenced by CalorimetryManager::EMShowerSimulation().

134 {crackPadProbability_ = val ;};

EcalHitMaker::crackPadProbability_

double crackPadProbability_

Definition: EcalHitMaker.h:279

heppy_batch.val

val

Definition: heppy_batch.py:351

void EcalHitMaker::setPreshowerPresent ( bool ps )

inline

Definition at line 137 of file EcalHitMaker.h.

References simulatePreshower_.

Referenced by CalorimetryManager::EMShowerSimulation().

137 {simulatePreshower_=ps;};

EcalHitMaker::simulatePreshower_

bool simulatePreshower_

Definition: EcalHitMaker.h:269

void EcalHitMaker::setPulledPadSurvivalProbability ( double val )

inline

Definition at line 132 of file EcalHitMaker.h.

References pulledPadProbability_, and heppy_batch::val.

Referenced by CalorimetryManager::EMShowerSimulation().

132 {pulledPadProbability_ = val;};

EcalHitMaker::pulledPadProbability_

double pulledPadProbability_

Definition: EcalHitMaker.h:277

heppy_batch.val

val

Definition: heppy_batch.py:351

void EcalHitMaker::setRadiusFactor ( double r )

inline

Definition at line 130 of file EcalHitMaker.h.

References alignCSCRings::r, and radiusCorrectionFactor_.

Referenced by CalorimetryManager::EMShowerSimulation().

130 {radiusCorrectionFactor_ = r;}

EcalHitMaker::radiusCorrectionFactor_

double radiusCorrectionFactor_

Definition: EcalHitMaker.h:259

alignCSCRings.r

r

Definition: alignCSCRings.py:93

void EcalHitMaker::setSpotEnergy ( double e )

inlineoverridevirtual

Implements CaloHitMaker.

Definition at line 117 of file EcalHitMaker.h.

References MillePedeFileConverter_cfg::e, getHits(), and CaloHitMaker::spotEnergy.

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

117 { spotEnergy=e;}

MillePedeFileConverter_cfg.e

e

Definition: MillePedeFileConverter_cfg.py:37

CaloHitMaker::spotEnergy

double spotEnergy

Definition: CaloHitMaker.h:51

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

Definition at line 328 of file EcalHitMaker.cc.

References buildSegments(), cellLine(), dir, DetId::Ecal, ecalFirstSegment_, ecalTotalX0(), CaloHitMaker::EMSHOWER, intersections_, maxX0_, myTrack_, ncrystals_, 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
 }

double EcalHitMaker::totalL0 ( ) const

inline

Number of interaction length "seen" by the track.

Definition at line 56 of file EcalHitMaker.h.

References totalL0_.

56 {return totalL0_;};

EcalHitMaker::totalL0_

double totalL0_

Definition: EcalHitMaker.h:204

double EcalHitMaker::totalX0 ( ) const

inline

Definition at line 53 of file EcalHitMaker.h.

References totalX0_.

Referenced by EMShower::prepareSteps().

53 {return totalX0_;};

EcalHitMaker::totalX0_

double totalX0_

Definition: EcalHitMaker.h:203

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

private

Definition at line 1328 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), plotBeamSpotDB::first, NONE, and 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;
 }

double EcalHitMaker::x0DepthOffset ( ) const

inline

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

Definition at line 59 of file EcalHitMaker.h.

References X0depthoffset_.

Referenced by EMShower::prepareSteps().

59 {return X0depthoffset_;}

EcalHitMaker::X0depthoffset_

double X0depthoffset_

Definition: EcalHitMaker.h:205

Member Data Documentation

double EcalHitMaker::bfactor_

private

Definition at line 267 of file EcalHitMaker.h.

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

std::vector<DetId> EcalHitMaker::CellsWindow_

private

Definition at line 236 of file EcalHitMaker.h.

Referenced by buildGeometry(), and EcalHitMaker().

int EcalHitMaker::central_

private

Definition at line 227 of file EcalHitMaker.h.

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

bool EcalHitMaker::configuredGeometry_

private

Definition at line 230 of file EcalHitMaker.h.

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

std::vector<XYZPoint> EcalHitMaker::corners

private

Definition at line 301 of file EcalHitMaker.h.

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

double EcalHitMaker::crackPadProbability_

private

Definition at line 279 of file EcalHitMaker.h.

Referenced by cracksPads(), and setCrackPadSurvivalProbability().

std::vector<CrystalPad> EcalHitMaker::crackpadsatdepth_

private

Definition at line 295 of file EcalHitMaker.h.

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

double EcalHitMaker::currentdepth_

private

Definition at line 265 of file EcalHitMaker.h.

Referenced by getPads().

bool EcalHitMaker::detailedShowerTail_

private

Definition at line 263 of file EcalHitMaker.h.

Referenced by getPads().

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

private

Definition at line 243 of file EcalHitMaker.h.

Referenced by buildGeometry(), and configureGeometry().

bool EcalHitMaker::doreorg_

private

Definition at line 291 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and getPads().

XYZPoint EcalHitMaker::EcalEntrance_

private

Definition at line 225 of file EcalHitMaker.h.

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

int EcalHitMaker::ecalFirstSegment_

private

Definition at line 248 of file EcalHitMaker.h.

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

unsigned EcalHitMaker::etasize_

private

Definition at line 251 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and reorganizePads().

bool EcalHitMaker::hitmaphasbeencalculated_

private

Definition at line 297 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and getHits().

std::vector<float> EcalHitMaker::hits_

private

Definition at line 238 of file EcalHitMaker.h.

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

std::vector<CaloPoint> EcalHitMaker::intersections_

private

Definition at line 287 of file EcalHitMaker.h.

Referenced by setTrackParameters().

double EcalHitMaker::L0ECAL_

private

Definition at line 215 of file EcalHitMaker.h.

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

double EcalHitMaker::L0EHGAP_

private

Definition at line 217 of file EcalHitMaker.h.

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

double EcalHitMaker::L0HCAL_

private

Definition at line 216 of file EcalHitMaker.h.

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

double EcalHitMaker::L0PS1_

private

Definition at line 212 of file EcalHitMaker.h.

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

double EcalHitMaker::L0PS2_

private

Definition at line 213 of file EcalHitMaker.h.

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

double EcalHitMaker::L0PS2EE_

private

Definition at line 214 of file EcalHitMaker.h.

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

double EcalHitMaker::maxX0_

private

Definition at line 219 of file EcalHitMaker.h.

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

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

private

Definition at line 300 of file EcalHitMaker.h.

Referenced by cracksPads(), and EcalHitMaker().

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

private

Definition at line 200 of file EcalHitMaker.h.

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

CrystalWindowMap* EcalHitMaker::myCrystalWindowMap_

private

Definition at line 245 of file EcalHitMaker.h.

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

const FSimTrack* EcalHitMaker::myTrack_

private

Definition at line 284 of file EcalHitMaker.h.

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

unsigned EcalHitMaker::ncrackpadsatdepth_

private

Definition at line 272 of file EcalHitMaker.h.

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

unsigned EcalHitMaker::ncrystals_

private

Definition at line 231 of file EcalHitMaker.h.

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

XYZNormal EcalHitMaker::normal_

private

Definition at line 226 of file EcalHitMaker.h.

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

unsigned EcalHitMaker::npadsatdepth_

private

Definition at line 273 of file EcalHitMaker.h.

Referenced by getPads(), and prepareCrystalNumberArray().

unsigned EcalHitMaker::nx_

private

Definition at line 233 of file EcalHitMaker.h.

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

unsigned EcalHitMaker::ny_

private

Definition at line 233 of file EcalHitMaker.h.

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

int EcalHitMaker::onEcal_

private

Definition at line 228 of file EcalHitMaker.h.

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

double EcalHitMaker::outsideWindowEnergy_

private

Definition at line 221 of file EcalHitMaker.h.

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

std::vector<CrystalPad> EcalHitMaker::padsatdepth_

private

Definition at line 294 of file EcalHitMaker.h.

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

unsigned EcalHitMaker::phisize_

private

Definition at line 252 of file EcalHitMaker.h.

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

Crystal EcalHitMaker::pivot_

private

Definition at line 224 of file EcalHitMaker.h.

Referenced by EcalHitMaker(), and getPads().

Plane3D EcalHitMaker::plan_

private

Definition at line 274 of file EcalHitMaker.h.

Referenced by getPads().

double EcalHitMaker::pulledPadProbability_

private

Definition at line 277 of file EcalHitMaker.h.

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

double EcalHitMaker::radiusCorrectionFactor_

private

Definition at line 259 of file EcalHitMaker.h.

Referenced by getPads(), and setRadiusFactor().

double EcalHitMaker::radiusFactor_

private

Definition at line 261 of file EcalHitMaker.h.

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

const RandomEngineAndDistribution* EcalHitMaker::random

private

Definition at line 303 of file EcalHitMaker.h.

Referenced by addHit().

double EcalHitMaker::rearleakage_

private

Definition at line 220 of file EcalHitMaker.h.

Referenced by addHitDepth(), and EcalHitMaker().

std::vector<Crystal> EcalHitMaker::regionOfInterest_

private

Definition at line 237 of file EcalHitMaker.h.

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

std::vector<CaloSegment> EcalHitMaker::segments_

private

Definition at line 289 of file EcalHitMaker.h.

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

bool EcalHitMaker::simulatePreshower_

private

Definition at line 269 of file EcalHitMaker.h.

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

double EcalHitMaker::sizex_

private

Definition at line 281 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

double EcalHitMaker::sizey_

private

Definition at line 281 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

double EcalHitMaker::totalL0_

private

Definition at line 204 of file EcalHitMaker.h.

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

double EcalHitMaker::totalX0_

private

Definition at line 203 of file EcalHitMaker.h.

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

bool EcalHitMaker::truncatedGrid_

private

Definition at line 254 of file EcalHitMaker.h.

Referenced by EcalHitMaker().

std::vector<bool> EcalHitMaker::validPads_

private

Definition at line 240 of file EcalHitMaker.h.

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

double EcalHitMaker::X0depthoffset_

private

Definition at line 205 of file EcalHitMaker.h.

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

double EcalHitMaker::X0ECAL_

private

Definition at line 209 of file EcalHitMaker.h.

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

double EcalHitMaker::X0EHGAP_

private

Definition at line 210 of file EcalHitMaker.h.

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

double EcalHitMaker::X0HCAL_

private

Definition at line 211 of file EcalHitMaker.h.

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

double EcalHitMaker::X0PS1_

private

Definition at line 206 of file EcalHitMaker.h.

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

double EcalHitMaker::X0PS2_

private

Definition at line 207 of file EcalHitMaker.h.

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

double EcalHitMaker::X0PS2EE_

private

Definition at line 208 of file EcalHitMaker.h.

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

double EcalHitMaker::xmax_

private

Definition at line 234 of file EcalHitMaker.h.

Referenced by getPads().

double EcalHitMaker::xmin_

private

Definition at line 234 of file EcalHitMaker.h.

Referenced by convertIntegerCoordinates(), and getPads().

double EcalHitMaker::ymax_

private

Definition at line 234 of file EcalHitMaker.h.

Referenced by getPads().

double EcalHitMaker::ymin_

private

Definition at line 234 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

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation