EcalHitMaker Class Reference

#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)
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 ()
	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)
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 ()
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(), CaloCellGeometry::getPosition(), 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_(NULL),
  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_ = 0; 
  // 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

(

)

Definition at line 126 of file EcalHitMaker.cc.

References myCrystalWindowMap_.

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

Member Function Documentation

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

virtual

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

{
//  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().

{
  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 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_, python.multivaluedict::sort(), 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;
//    }
}

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_, ntuplemaker::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.;
}

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

private

Definition at line 1208 of file EcalHitMaker.cc.

References sizex_, sizey_, xmin_, 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;
}

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

private

Definition at line 1305 of file EcalHitMaker.cc.

References CaloDirectionOperations::add2d(), 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_, ntuplemaker::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(), tree::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;
    }    
}

const XYZPoint& EcalHitMaker::ecalEntrance ( ) const

inline

used in FamosHcalHitMaker

Definition at line 128 of file EcalHitMaker.h.

References EcalEntrance_.

Referenced by HcalHitMaker::HcalHitMaker().

{return EcalEntrance_;};

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

{  return L0EHGAP_;}

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

{  return X0EHGAP_;}

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

{return L0ECAL_;}

double EcalHitMaker::ecalTotalX0 ( ) const

inline

in the ECAL

Definition at line 71 of file EcalHitMaker.h.

References X0ECAL_.

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

{return X0ECAL_;}

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_, newFWLiteAna::found, CrystalWindowMap::getCrystalWindow(), myCrystalNumberArray_, myCrystalWindowMap_, ncrackpadsatdepth_, nx_, ny_, padsatdepth_, findQualityFiles::size, ntuplemaker::status, and validPads_.

Referenced by addHit(), and addHitDepth().

{

  //  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 regionOfInterest_.

{return regionOfInterest_;}

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

{return myTrack_;}

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

virtual

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(), CaloCellGeometry::getPosition(), ecalpyutils::hashedIndex(), CaloHitMaker::hitMap_, hitmaphasbeencalculated_, hits_, PV3DBase< T, PVType, FrameType >::mag(), CaloHitMaker::myCalorimeter, ncrystals_, onEcal_, 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_;
}

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_, HLT_25ns14e33_v1_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_, jetcorrextractor::sign(), sizex_, sizey_, validPads_, X0depthoffset_, xmax_, xmin_, 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;
}

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 segments_.

Referenced by HcalHitMaker::setDepth().

{return segments_;};

double EcalHitMaker::getX0back ( ) const

inline

Definition at line 109 of file EcalHitMaker.h.

References maxX0_.

Referenced by EMShower::compute().

{return maxX0_;}

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

private

Definition at line 532 of file EcalHitMaker.cc.

References dir, eta, DetId::Hcal, FSimTrack::hcalEntrance(), Calorimeter::hcalProperties(), CaloHitMaker::myCalorimeter, myTrack_, FSimTrack::onHcal(), FSimTrack::onVFcal(), CalorimeterProperties::thickness(), 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));
    }

    }
}

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

{return L0HCAL_;}

double EcalHitMaker::hcalTotalX0 ( ) const

inline

in the HCAL

Definition at line 77 of file EcalHitMaker.h.

References X0HCAL_.

Referenced by EMShower::prepareSteps().

{return X0HCAL_;}

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 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(), 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;
}

double EcalHitMaker::ps1TotalL0 ( ) const

inline

total number of L0 in the PS (Layer1).

Definition at line 80 of file EcalHitMaker.h.

References L0PS1_.

{return L0PS1_;}

double EcalHitMaker::ps1TotalX0 ( ) const

inline

Definition at line 62 of file EcalHitMaker.h.

References X0PS1_.

Referenced by EMShower::prepareSteps().

{return X0PS1_;}

double EcalHitMaker::ps2eeTotalL0 ( ) const

inline

Definition at line 86 of file EcalHitMaker.h.

References L0PS2EE_.

{return L0PS2EE_;}

double EcalHitMaker::ps2eeTotalX0 ( ) const

inline

Definition at line 68 of file EcalHitMaker.h.

References X0PS2EE_.

Referenced by EMShower::prepareSteps().

{return X0PS2EE_;}

double EcalHitMaker::ps2TotalL0 ( ) const

inline

total number of L0 in the PS (Layer2).

Definition at line 83 of file EcalHitMaker.h.

References L0PS2_.

{return L0PS2_;}

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

{return X0PS2_;}

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

private

Definition at line 1182 of file EcalHitMaker.cc.

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

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 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_.

Referenced by CalorimetryManager::EMShowerSimulation().

{crackPadProbability_ = val ;};

void EcalHitMaker::setPreshowerPresent ( bool  ps )

inline

Definition at line 137 of file EcalHitMaker.h.

References simulatePreshower_.

Referenced by CalorimetryManager::EMShowerSimulation().

{simulatePreshower_=ps;};

void EcalHitMaker::setPulledPadSurvivalProbability ( double  val )

inline

Definition at line 132 of file EcalHitMaker.h.

References pulledPadProbability_.

Referenced by CalorimetryManager::EMShowerSimulation().

{pulledPadProbability_ = val;};

void EcalHitMaker::setRadiusFactor ( double  r )

inline

Definition at line 130 of file EcalHitMaker.h.

References alignCSCRings::r, and radiusCorrectionFactor_.

Referenced by CalorimetryManager::EMShowerSimulation().

{radiusCorrectionFactor_ = r;}

void EcalHitMaker::setSpotEnergy ( double  e )

inlinevirtual

Implements CaloHitMaker.

Definition at line 117 of file EcalHitMaker.h.

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

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

{ spotEnergy=e;}

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_.

{return totalL0_;}; 

double EcalHitMaker::totalX0 ( ) const

inline

Definition at line 53 of file EcalHitMaker.h.

References totalX0_.

Referenced by EMShower::prepareSteps().

{return totalX0_;}; 

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

{return X0depthoffset_;}

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