HcalLayerDepthMap Class Reference

#include <HcalLayerDepthMap.h>

Public Member Functions
int	getDepth (const int subdet, const int ieta, const int iphi, const int zside, const int layer) const
int	getDepth16 (const int subdet, const int iphi, const int zside) const
int	getDepthMax (const int subdet, const int iphi, const int zside) const
int	getDepthMax (const int subdet, const int ieta, const int iphi, const int zside) const
int	getDepthMin (const int subdet, const int iphi, const int zside) const
std::pair< int, int >	getDepths (const int eta) const
double	getLayer0Wt (const int subdet, const int iphi, const int zside) const
int	getLayerBack (const int subdet, const int ieta, const int iphi, const int zside, const int depth) const
void	getLayerDepth (const int subdet, const int ieta, const int iphi, const int zside, std::map< int, int > &layers) const
void	getLayerDepth (const int ieta, std::map< int, int > &layers) const
int	getLayerFront (const int subdet, const int ieta, const int iphi, const int zside, const int depth) const
int	getMaxDepthLastHE (const int subdet, const int iphi, const int zside) const
const std::vector< int > &	getPhis () const
int	getSubdet () const
	HcalLayerDepthMap ()
void	initialize (const int subdet, const int ietaMax, const int dep16C, const int dep29C, const double wtl0C, std::vector< int > const &iphi, std::vector< int > const &ieta, std::vector< int > const &layer, std::vector< int > const &depth)
bool	isValid (const int det, const int phi, const int zside) const
int	validDet (std::vector< int > &phis) const
std::pair< int, int >	validEta () const
	~HcalLayerDepthMap ()

Private Attributes
int	dep16C_
int	dep29C_
std::map< std::pair< int, int >, int >	depth2LayerB_
std::map< std::pair< int, int >, int >	depth2LayerF_
int	depthMax_
int	depthMin_
std::map< int, std::pair< int, int > >	depthsEta_
int	ietaMax_
int	ietaMin_
std::vector< int >	iphi_
std::map< std::pair< int, int >, int >	layer2Depth_
int	subdet_
double	wtl0C_

Static Private Attributes
static const int	maxLayers_ = 18

Detailed Description

this class stores the map of layer to depth for special phi sections

Author

Sunanda Banerjee, SINP sunan.nosp@m.da.b.nosp@m.anerj.nosp@m.ee@c.nosp@m.ern.c.nosp@m.h

Definition at line 18 of file HcalLayerDepthMap.h.

Constructor & Destructor Documentation

HcalLayerDepthMap::HcalLayerDepthMap

(

)

Definition at line 9 of file HcalLayerDepthMap.cc.

References dep29C_, depthMax_, depthMin_, ietaMax_, ietaMin_, subdet_, and wtl0C_.

                                     {
  subdet_  = 0;
  ietaMin_ = ietaMax_ = 0;
  depthMin_= 99;
  depthMax_=-1;
  dep29C_  = 2;
  wtl0C_   = 1.;
}

HcalLayerDepthMap::~HcalLayerDepthMap

(

)

Definition at line 18 of file HcalLayerDepthMap.cc.

{}

Member Function Documentation

int HcalLayerDepthMap::getDepth	(	const int	subdet,
		const int	ieta,
		const int	iphi,
		const int	zside,
		const int	layer
	)		const

Definition at line 110 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, isValid(), and layer2Depth_.

Referenced by HcalDDDSimConstants::findDepth(), getDepthMax(), and getSubdet().

                                       {
  int depth(-1);
  if (isValid(subdet,iphi,zside)) {
    std::map<std::pair<int,int>,int>::const_iterator itr = layer2Depth_.find(std::pair<int,int>(ieta,layer));
    if (itr != layer2Depth_.end()) depth = itr->second;
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getDepth::Input " << subdet << ":" << ieta << ":" << iphi 
        << ":" << zside << ":" << layer << " Output " << depth <<std::endl;
#endif
  return depth;
}

int HcalLayerDepthMap::getDepth16	(	const int	subdet,
		const int	iphi,
		const int	zside
	)		const

Definition at line 125 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, dep16C_, particleFlowClusterECALTimeSelected_cfi::depth, and isValid().

Referenced by HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta16M(), and getSubdet().

                                         {
  int depth(-1);
  if (isValid(subdet,iphi,zside)) depth = dep16C_;
#ifdef EDM_ML_DEBUG
  std::cout << "getDepth16::Input " << subdet << ":" << iphi << ":" << zside
        << " Output "  << depth << std::endl;
#endif
  return depth;
}

int HcalLayerDepthMap::getDepthMax	(	const int	subdet,
		const int	iphi,
		const int	zside
	)		const

Definition at line 146 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, depthMax_, and isValid().

Referenced by getSubdet(), and HcalDDDSimConstants::initialize().

                                          {
  int depth = (isValid(subdet,iphi,zside)) ? depthMax_ : -1;
#ifdef EDM_ML_DEBUG
  std::cout << "getDepthMax::Input " << subdet << ":" << iphi << ":" << zside
        << " Output " << depth << std::endl;
#endif
  return depth;
}

int HcalLayerDepthMap::getDepthMax	(	const int	subdet,
		const int	ieta,
		const int	iphi,
		const int	zside
	)		const

Definition at line 156 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, getDepth(), isValid(), and maxLayers_.

                                                          {
  int depth = (isValid(subdet,iphi,zside)) ? getDepth(subdet,ieta,iphi,zside,maxLayers_) : -1;
#ifdef EDM_ML_DEBUG
  std::cout << "getDepthMax::Input " << subdet << ":" << iphi << ":" << zside
        << " Output " << depth << std::endl;
#endif
  return depth;
}

int HcalLayerDepthMap::getDepthMin	(	const int	subdet,
		const int	iphi,
		const int	zside
	)		const

Definition at line 136 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, depthMin_, and isValid().

Referenced by getSubdet().

                                          {
  int depth = (isValid(subdet,iphi,zside)) ? depthMin_ : -1;
#ifdef EDM_ML_DEBUG
  std::cout << "getDepthMin::Input " << subdet << ":" << iphi << ":" << zside
        << " Output " << depth << std::endl;
#endif
  return depth;
}

std::pair< int, int > HcalLayerDepthMap::getDepths

(

const int

eta

)

const

Definition at line 166 of file HcalLayerDepthMap.cc.

References depthsEta_.

Referenced by HcalDDDSimConstants::getMaxDepth(), HcalDDDSimConstants::getMinDepth(), getSubdet(), and HcalDDDRecConstants::initialize().

                                                                 {

  std::map<int,std::pair<int,int> >::const_iterator itr = depthsEta_.find(eta);
  if (itr == depthsEta_.end()) return std::pair<int,int>(-1,-1);
  else                         return itr->second;
}

double HcalLayerDepthMap::getLayer0Wt	(	const int	subdet,
		const int	iphi,
		const int	zside
	)		const

Definition at line 254 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, isValid(), and wtl0C_.

Referenced by HcalDDDSimConstants::getLayer0Wt(), and getSubdet().

                                             {
  double wt = isValid(subdet,iphi,zside) ? wtl0C_ : -1.0;
#ifdef EDM_ML_DEBUG
  std::cout << "getLayer0Wt::Input " << subdet << ":" << iphi << ":" << zside 
        << " Output " << wt << std::endl;
#endif
  return wt;
}

int HcalLayerDepthMap::getLayerBack	(	const int	subdet,
		const int	ieta,
		const int	iphi,
		const int	zside,
		const int	depth
	)		const

Definition at line 188 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, depth2LayerB_, and isValid().

Referenced by HcalDDDSimConstants::getLayerBack(), getSubdet(), and HcalDDDRecConstants::initialize().

                                           {
  int layer(-1);
  if (isValid(subdet,iphi,zside)) {
    std::map<std::pair<int,int>,int>::const_iterator itr = depth2LayerB_.find(std::pair<int,int>(ieta,depth));
    if (itr != depth2LayerB_.end()) layer = itr->second;
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getLayerBack::Input " << subdet << ":" << ieta << ":" << iphi 
        << ":" << zside << ":" << depth << " Output " << layer <<std::endl;
#endif
  return layer;
}

void HcalLayerDepthMap::getLayerDepth	(	const int	subdet,
		const int	ieta,
		const int	iphi,
		const int	zside,
		std::map< int, int > &	layers
	)		const

Definition at line 203 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, stringResolutionProvider_cfi::eta, isValid(), layer2Depth_, genParticles_cff::map, and edm::second().

Referenced by HcalDDDRecConstants::getDepth(), HcalDDDRecConstants::getEtaBins(), and getSubdet().

                                                     {
  layers.clear();
  if (isValid(subdet,phi,zside)) {
    for (std::map<std::pair<int,int>,int>::const_iterator itr=layer2Depth_.begin();
     itr != layer2Depth_.end(); ++itr) {
      if ((itr->first).first == eta) {
    layers[((itr->first).second)+1] = (itr->second);
      }
    }
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getLayerDepth::Input " << subdet << ":" << eta << ":" << phi
        << ":" << zside << " Output " << layers.size() << " entries\n";
  for (std::map<int,int>::iterator itr=layers.begin(); itr != layers.end();
       ++itr) std::cout << " [" << itr->first << "] " << itr->second;
  std::cout << std::endl;
#endif
}

void HcalLayerDepthMap::getLayerDepth	(	const int	ieta,
		std::map< int, int > &	layers
	)		const

Definition at line 224 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, stringResolutionProvider_cfi::eta, layer2Depth_, genParticles_cff::map, edm::second(), and subdet_.

                                                     {
  layers.clear();
  if (subdet_ > 0) {
    for (std::map<std::pair<int,int>,int>::const_iterator itr=layer2Depth_.begin();
     itr != layer2Depth_.end(); ++itr) {
      if ((itr->first).first == eta) {
    layers[((itr->first).second)+1] = (itr->second);
      }
    }
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getLayerDepth::Input " << eta << " Output " << layers.size()
        << " entries\n";
  for (std::map<int,int>::iterator itr=layers.begin(); itr != layers.end();
       ++itr) std::cout << " [" << itr->first << "] " << itr->second;
  std::cout << std::endl;
#endif
}

int HcalLayerDepthMap::getLayerFront	(	const int	subdet,
		const int	ieta,
		const int	iphi,
		const int	zside,
		const int	depth
	)		const

Definition at line 173 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, depth2LayerF_, and isValid().

Referenced by HcalDDDSimConstants::getLayerFront(), getSubdet(), and HcalDDDRecConstants::initialize().

                                            {
  int layer(-1);
  if (isValid(subdet,iphi,zside)) {
    std::map<std::pair<int,int>,int>::const_iterator itr = depth2LayerF_.find(std::pair<int,int>(ieta,depth));
    if (itr != depth2LayerF_.end()) layer = itr->second;
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getLayerFront::Input " << subdet << ":" << ieta << ":" << iphi 
        << ":" << zside << ":" << depth << " Output " << layer <<std::endl;
#endif
  return layer;
}

int HcalLayerDepthMap::getMaxDepthLastHE	(	const int	subdet,
		const int	iphi,
		const int	zside
	)		const

Definition at line 244 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, dep29C_, particleFlowClusterECALTimeSelected_cfi::depth, and isValid().

Referenced by HcalDDDSimConstants::getDepthEta29(), HcalDDDSimConstants::getDepthEta29M(), and getSubdet().

                                            {
  int depth = isValid(subdet,iphi,zside) ? dep29C_ : -1;
#ifdef EDM_ML_DEBUG
  std::cout << "getMaxDepthLastHE::Input " << subdet << ":" << iphi << ":"
        << zside << " Output " << depth << std::endl;
#endif
  return depth;
}

const std::vector<int>& HcalLayerDepthMap::getPhis ( ) const

inline

Definition at line 47 of file HcalLayerDepthMap.h.

References iphi_, isValid(), phi, and validDet().

Referenced by HcalDDDRecConstants::getPhiZOne().

{return iphi_;}

int HcalLayerDepthMap::getSubdet ( ) const

inline

Definition at line 28 of file HcalLayerDepthMap.h.

References stringResolutionProvider_cfi::eta, getDepth(), getDepth16(), getDepthMax(), getDepthMin(), getDepths(), getLayer0Wt(), getLayerBack(), getLayerDepth(), getLayerFront(), getMaxDepthLastHE(), LayerTriplets::layers(), subdet_, and HcalDigiParam_cfi::zside.

Referenced by HcalDDDRecConstants::getPhiZOne(), and HcalDDDRecConstants::withSpecialRBXHBHE().

{return subdet_;}

void HcalLayerDepthMap::initialize	(	const int	subdet,
		const int	ietaMax,
		const int	dep16C,
		const int	dep29C,
		const double	wtl0C,
		std::vector< int > const &	iphi,
		std::vector< int > const &	ieta,
		std::vector< int > const &	layer,
		std::vector< int > const &	depth
	)

Definition at line 20 of file HcalLayerDepthMap.cc.

References gather_cfg::cout, dep16C_, dep29C_, depth2LayerB_, depth2LayerF_, depthMax_, depthMin_, depthsEta_, stringResolutionProvider_cfi::eta, ietaMax_, ietaMin_, iphi_, gen::k, layer2Depth_, genParticles_cff::map, maxLayers_, edm::second(), subdet_, and wtl0C_.

Referenced by HcalDDDSimConstants::initialize().

                                                {
  
  subdet_  = subdet;
  ietaMin_ = ietaMax_ = ietaMax;
  dep16C_  = dep16C;
  dep29C_  = dep29C;
  wtl0C_   = wtl0C;
  iphi_.insert(iphi_.end(),iphi.begin(),iphi.end());
  layer2Depth_.clear(); depth2LayerF_.clear(); depth2LayerB_.clear();
  depthMin_= 99;
  depthMax_=-1;
  for (unsigned int k=0; k<ieta.size(); ++k) {
    if (ieta[k]  < ietaMin_)  ietaMin_  = ieta[k];
    if (depth[k] < depthMin_) depthMin_ = depth[k];
    if (depth[k] > depthMax_) depthMax_ = depth[k];
  }
  //Assume ieta, layer, depth are in increasing order of ieta and depth
  for (unsigned int k1=0; k1<ieta.size(); ++k1) {
    int ietaMin = ieta[k1];
    int ietaMax = ietaMax_;
    int layMin  = layer[k1];
    int layMax  = (k1+1 < ieta.size()) ? (layer[k1+1]-1) : maxLayers_;
    for (unsigned int k2=k1+1; k2<ieta.size(); ++k2) {
      if (ieta[k2] > ieta[k1]) {
    ietaMax = ieta[k2] - 1;
    if (k2 == k1+1) layMax = maxLayers_;
    break;
      }
    }
    for (int eta=ietaMin; eta<=ietaMax; ++eta) {
      depth2LayerF_[std::pair<int,int>(eta,depth[k1])] = layMin;
      depth2LayerB_[std::pair<int,int>(eta,depth[k1])] = layMax;
      for (int lay=layMin; lay<=layMax; ++lay)
    layer2Depth_[std::pair<int,int>(eta,lay)] = depth[k1];
    }
  }
  for (int eta=ietaMin_; eta<=ietaMax_; ++eta) {
    int dmin(99), dmax(-1);
    for (std::map<std::pair<int,int>,int>::iterator itr=layer2Depth_.begin();
     itr != layer2Depth_.end(); ++itr) {
      if ((itr->first).first == eta) {
    if ((itr->second) < dmin) dmin = (itr->second);
    if ((itr->second) > dmax) dmax = (itr->second);
      }
    }
    if (subdet == 2) {
      if      (eta == ietaMin_) dmin = dep16C_;
      else if (eta == ietaMax_) dmax = dep29C_;
    }
    depthsEta_[eta] = std::pair<int,int>(dmin,dmax);
  }
#ifdef EDM_ML_DEBUG
  std::cout << "HcalLayerDepthMap: Subdet " << subdet_ << " iEta " << ietaMin_
        << ":" << ietaMax_ << " depth " << depthMin_ << ":" << depthMax_
        << "\nMaximum Depth for last HE towers " << dep29C_  
        << " Layer 0 Weight " << wtl0C_ << " iPhi";
  for (unsigned int k=0; k<iphi_.size(); ++k) std::cout << ":" << iphi_[k];
  std::cout << std::endl << "Layer2Depth_ with " << layer2Depth_.size()
        << " elements" << std::endl;
  for (std::map<std::pair<int,int>,int>::iterator itr=layer2Depth_.begin();
       itr != layer2Depth_.end(); ++itr)
    std::cout << "iEta " << (itr->first).first << " Layer " 
          << (itr->first).second << " Depth " << itr->second << std::endl;
  std::cout << "Depth2LayerFront with " << depth2LayerF_.size() 
        << " elemsts" << std::endl;
  for (std::map<std::pair<int,int>,int>::iterator itr=depth2LayerF_.begin();
       itr != depth2LayerF_.end(); ++itr)
    std::cout << "iEta " << (itr->first).first << " Depth " 
          << (itr->first).second << " Layer " << itr->second << std::endl;
  std::cout << "Depth2LayerBack with " << depth2LayerB_.size()
        << " elemets" << std::endl;
  for (std::map<std::pair<int,int>,int>::iterator itr=depth2LayerB_.begin();
       itr != depth2LayerB_.end(); ++itr)
    std::cout << "iEta " << (itr->first).first << " Depth " 
          << (itr->first).second << " Layer " << itr->second << std::endl;
  std::cout << "DepthsEta_ with " << depthsEta_.size() << " elements\n";
  for (std::map<int,std::pair<int,int> >::iterator itr=depthsEta_.begin();
       itr != depthsEta_.end(); ++itr)
    std::cout << "iEta " << itr->first << " Depths " 
          << (itr->second).first << ":" << (itr->second).second 
          << std::endl;
#endif
}

bool HcalLayerDepthMap::isValid	(	const int	det,
		const int	phi,
		const int	zside
	)		const

Definition at line 264 of file HcalLayerDepthMap.cc.

References spr::find(), RemoveAddSevLevel::flag, iphi_, and subdet_.

Referenced by ntuple._Object::_checkIsValid(), HcalDDDSimConstants::findDepth(), getDepth(), getDepth16(), getDepthMax(), getDepthMin(), getLayer0Wt(), getLayerBack(), getLayerDepth(), getLayerFront(), HcalDDDSimConstants::getMaxDepth(), getMaxDepthLastHE(), HcalDDDSimConstants::getMinDepth(), getPhis(), and core.AutoHandle.AutoHandle::ReallyLoad().

                                       {
  bool flag(false);
  int  kphi = (zside > 0) ? iphi : -iphi;
  if (subdet == subdet_)
    flag = (std::find(iphi_.begin(),iphi_.end(),kphi) != iphi_.end());
  return flag;
}

int HcalLayerDepthMap::validDet

(

std::vector< int > &

phis

)

const

Definition at line 273 of file HcalLayerDepthMap.cc.

References iphi_, and subdet_.

Referenced by HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29M(), HcalDDDRecConstants::getEtaBins(), getPhis(), HcalDDDRecConstants::initialize(), and HcalDDDSimConstants::printTiles().

                                                         {

  phi.insert(phi.end(),iphi_.begin(),iphi_.end());
  return subdet_;
}

std::pair<int,int> HcalLayerDepthMap::validEta ( ) const

inline

Definition at line 50 of file HcalLayerDepthMap.h.

References ietaMax_, and ietaMin_.

Referenced by HcalDDDRecConstants::initialize().

{return std::pair<int,int>(ietaMin_,ietaMax_);}

Member Data Documentation

int HcalLayerDepthMap::dep16C_

private

Definition at line 59 of file HcalLayerDepthMap.h.

Referenced by getDepth16(), and initialize().

int HcalLayerDepthMap::dep29C_

private

Definition at line 60 of file HcalLayerDepthMap.h.

Referenced by getMaxDepthLastHE(), HcalLayerDepthMap(), and initialize().

std::map<std::pair<int,int>,int> HcalLayerDepthMap::depth2LayerB_

private

Definition at line 65 of file HcalLayerDepthMap.h.

Referenced by getLayerBack(), and initialize().

std::map<std::pair<int,int>,int> HcalLayerDepthMap::depth2LayerF_

private

Definition at line 64 of file HcalLayerDepthMap.h.

Referenced by getLayerFront(), and initialize().

int HcalLayerDepthMap::depthMax_

private

Definition at line 58 of file HcalLayerDepthMap.h.

Referenced by getDepthMax(), HcalLayerDepthMap(), and initialize().

int HcalLayerDepthMap::depthMin_

private

Definition at line 57 of file HcalLayerDepthMap.h.

Referenced by getDepthMin(), HcalLayerDepthMap(), and initialize().

std::map<int,std::pair<int,int> > HcalLayerDepthMap::depthsEta_

private

Definition at line 66 of file HcalLayerDepthMap.h.

Referenced by getDepths(), and initialize().

int HcalLayerDepthMap::ietaMax_

private

Definition at line 56 of file HcalLayerDepthMap.h.

Referenced by HcalLayerDepthMap(), initialize(), and validEta().

int HcalLayerDepthMap::ietaMin_

private

Definition at line 55 of file HcalLayerDepthMap.h.

Referenced by HcalLayerDepthMap(), initialize(), and validEta().

std::vector<int> HcalLayerDepthMap::iphi_

private

Definition at line 62 of file HcalLayerDepthMap.h.

Referenced by getPhis(), initialize(), isValid(), and validDet().

std::map<std::pair<int,int>,int> HcalLayerDepthMap::layer2Depth_

private

Definition at line 63 of file HcalLayerDepthMap.h.

Referenced by getDepth(), getLayerDepth(), and initialize().

const int HcalLayerDepthMap::maxLayers_ = 18

staticprivate

Definition at line 53 of file HcalLayerDepthMap.h.

Referenced by getDepthMax(), and initialize().

int HcalLayerDepthMap::subdet_

private

Definition at line 54 of file HcalLayerDepthMap.h.

Referenced by getLayerDepth(), getSubdet(), HcalLayerDepthMap(), initialize(), isValid(), and validDet().

double HcalLayerDepthMap::wtl0C_

private

Definition at line 61 of file HcalLayerDepthMap.h.

Referenced by getLayer0Wt(), HcalLayerDepthMap(), and initialize().