CalibrationCluster.cc

Go to the documentation of this file.

#include "Calibration/Tools/interface/CalibrationCluster.h"
#include <iostream>
#include <string>

using namespace std;

CalibrationCluster::CalibrationCluster() {}

CalibrationCluster::~CalibrationCluster() {}


std::vector<EBDetId> CalibrationCluster::get5x5Id(EBDetId const& maxHitId) {
  Xtals5x5.clear();

  //   std::cout << "get5x5Id: max Containment crystal " << maxHitId.ic() << " eta " << maxHitId.ieta() << " phi " << maxHitId.iphi() << std::endl;

  for (unsigned int icry = 0; icry < 25; icry++) {
    unsigned int row = icry / 5;
    unsigned int column = icry % 5;
    //       std::cout << "CalibrationCluster::icry = " << icry << std::endl;

    int curr_eta = maxHitId.ieta() + column - (5 / 2);
    int curr_phi = maxHitId.iphi() + row - (5 / 2);

    if (curr_eta * maxHitId.ieta() <= 0) {
      if (maxHitId.ieta() > 0)
        curr_eta--;
      else
        curr_eta++;
    }  // JUMP over 0
    if (curr_phi < 1)
      curr_phi += 360;
    if (curr_phi > 360)
      curr_phi -= 360;

    try {
      //         Xtals5x5.push_back(EBDetId(maxHitId.ieta()+column-2,maxHitId.iphi()+row-2,EBDetId::ETAPHIMODE));
      Xtals5x5.push_back(EBDetId(curr_eta, curr_phi, EBDetId::ETAPHIMODE));
    } catch (...) {
      std::cout << "Cannot construct 5x5 matrix around EBDetId " << maxHitId << std::endl;
    }
  }

  return Xtals5x5;
}


std::vector<EBDetId> CalibrationCluster::get3x3Id(EBDetId const& maxHitId) {
  Xtals3x3.clear();

  for (unsigned int icry = 0; icry < 9; icry++) {
    unsigned int row = icry / 3;
    unsigned int column = icry % 3;

    try {
      Xtals3x3.push_back(EBDetId(maxHitId.ieta() + column - 1, maxHitId.iphi() + row - 1, EBDetId::ETAPHIMODE));
    } catch (...) {
      std::cout << "Cannot construct 3x3 matrix around EBDetId " << maxHitId << std::endl;
    }
  }

  return Xtals3x3;
}


CalibrationCluster::CalibMap CalibrationCluster::getMap(int minEta, int maxEta, int minPhi, int maxPhi) {
  calibRegion.clear();
  int rowSize = maxEta - minEta + 1;
  int columnSize = maxPhi - minPhi + 1;
  int reducedSize = rowSize * columnSize;

  for (int icry = 0; icry < reducedSize; icry++) {
    unsigned int eta = minEta + icry / columnSize;
    unsigned int phi = minPhi + icry % columnSize;

    try {
      calibRegion.insert(pippo(EBDetId(eta, phi, EBDetId::ETAPHIMODE), icry));
    } catch (...) {
      std::cout << "Cannot construct full matrix !!! " << std::endl;
    }
  }

  return calibRegion;
}


std::vector<float> CalibrationCluster::getEnergyVector(const EBRecHitCollection* hits,
                                                       CalibMap& ReducedMap,
                                                       std::vector<EBDetId>& XstalsNxN,
                                                       float& outBoundEnergy,
                                                       int& nXtalsOut) {
  energyVector.clear();
  std::vector<EBDetId>::iterator it;

  // std::cout << "Reduced Map Size =" << ReducedMap.size() << std::endl;
  // std::cout << "XstalsNxN Size =" << XstalsNxN.size() << std::endl;
  energyVector.resize(ReducedMap.size(), 0.);

  outBoundEnergy = 0.;
  nXtalsOut = 0;
  for (it = XstalsNxN.begin(); it != XstalsNxN.end(); ++it) {
    if (ReducedMap.find(*it) != ReducedMap.end()) {
      CalibMap::iterator it2 = ReducedMap.find(*it);

      int icry = it2->second;

      energyVector[icry] = (hits->find(*it))->energy();

    } else {
      //       std::cout << " Cell out of Reduced map: did you subtracted the cell energy from P ???" << std::endl;
      outBoundEnergy += (hits->find(*it))->energy();
      nXtalsOut++;
    }
  }

  return energyVector;
}