CrystalWindowMap.cc

Go to the documentation of this file.

//FAMOS headers
#include "FastSimulation/CaloGeometryTools/interface/CrystalWindowMap.h"
#include "FastSimulation/CaloGeometryTools/interface/CaloGeometryHelper.h"

#include <algorithm>

CrystalWindowMap::CrystalWindowMap(const CaloGeometryHelper *calo,const std::vector<Crystal> & cw):
  myCalorimeter_(calo),originalVector_(cw)
{
  size_=cw.size();
  if(size_==0) return;
  
  // Loop over the crystals in the grid and stores the index number of the 8 neighbours
  myNeighbours_.resize(size_);

  for(unsigned ic=0; ic<size_; ++ic)
    {
      const std::vector<DetId>& neighbours=myCalorimeter_->getNeighbours(cw[ic].getDetId());
      myNeighbours_[ic].reserve(8);
      for(unsigned in=0; in<8;++in)
    {
      // Check that the crystal is in the grid

      Crystal::crystalEqual myCrystal(neighbours[in]);
      std::vector<Crystal>::const_iterator itcheck;
      itcheck=find_if(cw.begin(),cw.end(),myCrystal);
      // The neighbour might not be in the grid
      if(itcheck==cw.end())
        {
//        std::cout << " Ouh la " << std::endl;
//        for(unsigned ic=0;ic<size_;++ic)
//      {
//        std::cout << cw[ic].getDetId().rawId()<<  " " ;
//      }
//        std::cout << std::endl ; 
//        std::cout << " We are looking for " << neighbours[in].rawId() << std::endl;
//        edm::LogWarning("CrystalWindowMap") << " Inconsistency in the CellWindow " << std::endl;      
        }
      else
        {
          myNeighbours_[ic].push_back(itcheck-cw.begin());
          //          std::cout << " index " << itcheck-cw.begin() << std::endl;
        }
    }
    }
}

bool 
CrystalWindowMap::getCrystalWindow(unsigned iq,std::vector<unsigned>&  cw ) const
{
  if(iq<size_) // iq >= 0, since iq is unsigned
    {
      cw=myNeighbours_[iq];
      return true;
    }
  else
    return false;
}

bool 
CrystalWindowMap::getCrystalWindow(unsigned iq,const std::vector<unsigned>*  cw) const
{
  
  if(iq<size_) // iq >= 0, since iq is unsigned
    {
      cw=&myNeighbours_[iq];
      return true;
    }
  else
    return false;
    //  std::map<CrystalID,CrystalWindow>::const_iterator itcheck=myMap.find(cell);
}

const std::vector<unsigned>& 
CrystalWindowMap::getCrystalWindow(unsigned iq, bool& status) const
{
  return myNeighbours_[iq];
}