EcalTBCrystalMap.cc

Go to the documentation of this file.

#include "Geometry/EcalTestBeam/interface/EcalTBCrystalMap.h"

EcalTBCrystalMap::EcalTBCrystalMap(std::string const & MapFileName) {
                                                                                                                                          
  // check if File exists
  std::ifstream * input= new std::ifstream(MapFileName.c_str(), std::ios::in);
  if(! (*input))
    {
      throw cms::Exception("FileNotFound", "Ecal TB Crystal map file not found") 
        <<  "\n" << MapFileName << " could not be opened.\n";
   }
  if (input){

    int nCrysCount = 0;

    while (nCrysCount <= NCRYSTAL ) {
      
      (*input) >> crysIndex >> crysEta >>  crysPhi;
      map_[std::pair<double,double>(crysEta,crysPhi)] = crysIndex;

      nCrysCount++;
      
    }

    input->close();
  }

}

EcalTBCrystalMap::~EcalTBCrystalMap() {
}

int EcalTBCrystalMap::CrystalIndex(double thisEta, double thisPhi) {

  int thisCrysIndex = 0;

  CrystalTBIndexMap::const_iterator mapItr = map_.find(std::make_pair(thisEta,thisPhi));
  if ( mapItr != map_.end() ) {
    thisCrysIndex = mapItr->second;
  }

  return thisCrysIndex;

}

void EcalTBCrystalMap::findCrystalAngles(const int thisCrysIndex, double & thisEta, double & thisPhi) {
  
  thisEta = thisPhi = 0.;

  if ( thisCrysIndex < 1 || thisCrysIndex > NCRYSTAL ) {
    edm::LogError("OutOfBounds") << "Required crystal number " << thisCrysIndex << " outside range";
    return;
  }

  for ( CrystalTBIndexMap::const_iterator mapItr = map_.begin() ; mapItr != map_.end() ; ++mapItr) {
    int theCrysIndex = mapItr->second;
    if ( theCrysIndex == thisCrysIndex ) {
      thisEta = (mapItr->first).first;
      thisPhi = (mapItr->first).second;
      return;
    }
  }

}