EcalFloatCondObjectContainerUtils.h

Go to the documentation of this file.

#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "CondFormats/EcalObjects/interface/EcalCondObjectContainer.h"

#include "TH2F.h"
#include <string>

enum { kEBChannels = 61200, kEEChannels = 14648 };
enum {
  MIN_IETA = 1,
  MIN_IPHI = 1,
  MAX_IETA = 85,
  MAX_IPHI = 360,
  EBhistEtaMax = 171
};  // barrel lower and upper bounds on eta and phi
enum {
  IX_MIN = 1,
  IY_MIN = 1,
  IX_MAX = 100,
  IY_MAX = 100,
  EEhistXMax = 220
};  // endcaps lower and upper bounds on x and y

template <class floatCondObj>
void fillEBMap_SingleIOV(std::shared_ptr<floatCondObj> payload, TH2F*& barrel) {
  for (int cellid = EBDetId::MIN_HASH; cellid < EBDetId::kSizeForDenseIndexing; ++cellid) {
    uint32_t rawid = EBDetId::unhashIndex(cellid);
    EcalCondObjectContainer<float>::const_iterator value_ptr = payload->find(rawid);
    if (value_ptr == payload->end())
      continue;  // cell absent from payload

    float weight = (float)(*value_ptr);
    Double_t phi = (Double_t)(EBDetId(rawid)).iphi() - 0.5;
    Double_t eta = (Double_t)(EBDetId(rawid)).ieta();

    if (eta > 0.)
      eta = eta - 0.5;  //   0.5 to 84.5
    else
      eta = eta + 0.5;  //  -84.5 to -0.5

    barrel->Fill(phi, eta, weight);
  }  // loop over cellid
}

template <class floatCondObj>
void fillEEMap_SingleIOV(std::shared_ptr<floatCondObj> payload, TH2F*& endc_m, TH2F*& endc_p) {
  // looping over the EE channels
  for (int iz = -1; iz < 2; iz = iz + 2)  // -1 or +1
    for (int iy = IY_MIN; iy < IY_MAX + IY_MIN; iy++)
      for (int ix = IX_MIN; ix < IX_MAX + IX_MIN; ix++)
        if (EEDetId::validDetId(ix, iy, iz)) {
          EEDetId myEEId = EEDetId(ix, iy, iz, EEDetId::XYMODE);
          uint32_t rawid = myEEId.rawId();
          EcalCondObjectContainer<float>::const_iterator value_ptr = payload->find(rawid);
          if (value_ptr == payload->end())
            continue;  // cell absent from payload

          float weight = (float)(*value_ptr);
          if (iz == 1)
            endc_p->Fill(ix, iy, weight);
          else
            endc_m->Fill(ix, iy, weight);

        }  // validDetId
}

template <class floatCondObj>
void fillEBMap_TwoIOVs(std::shared_ptr<floatCondObj> payload,
                       TH2F*& barrel,
                       int irun,
                       float pEB[],
                       float& pEBmin,
                       float& pEBmax,
                       int method) {
  for (int cellid = EBDetId::MIN_HASH; cellid < EBDetId::kSizeForDenseIndexing; ++cellid) {
    uint32_t rawid = EBDetId::unhashIndex(cellid);
    EcalCondObjectContainer<float>::const_iterator value_ptr = payload->find(rawid);
    if (value_ptr == payload->end())
      continue;  // cell absent from payload

    float weight = (float)(*value_ptr);
    if (irun == 0)
      pEB[cellid] = weight;
    else {
      Double_t phi = (Double_t)(EBDetId(rawid)).iphi() - 0.5;
      Double_t eta = (Double_t)(EBDetId(rawid)).ieta();

      if (eta > 0.)
        eta = eta - 0.5;  //   0.5 to 84.5
      else
        eta = eta + 0.5;  //  -84.5 to -0.5
      double dr;
      if (method == 0) {
        dr = weight - pEB[cellid];
      }  // diff
      else {
        if (pEB[cellid] != 0.)
          dr = weight / pEB[cellid];
        else {
          if (weight == 0.)
            dr = 1.;
          else
            dr = 9999.;  // use a large value
        }
      }  // ratio
      if (dr < pEBmin)
        pEBmin = dr;
      if (dr > pEBmax)
        pEBmax = dr;
      barrel->Fill(phi, eta, dr);
    }
  }  // loop over cellid
}

template <class floatCondObj>
void fillEEMap_TwoIOVs(std::shared_ptr<floatCondObj> payload,
                       TH2F*& endc_m,
                       TH2F*& endc_p,
                       int irun,
                       float pEE[],
                       float& pEEmin,
                       float& pEEmax,
                       int method) {
  // looping over the EE channels
  for (int iz = -1; iz < 2; iz = iz + 2)  // -1 or +1
    for (int iy = IY_MIN; iy < IY_MAX + IY_MIN; iy++)
      for (int ix = IX_MIN; ix < IX_MAX + IX_MIN; ix++)
        if (EEDetId::validDetId(ix, iy, iz)) {
          EEDetId myEEId = EEDetId(ix, iy, iz, EEDetId::XYMODE);
          uint32_t cellid = myEEId.hashedIndex();
          uint32_t rawid = myEEId.rawId();
          EcalCondObjectContainer<float>::const_iterator value_ptr = payload->find(rawid);

          if (value_ptr == payload->end())
            continue;  // cell absent from payload
          float weight = (float)(*value_ptr);

          if (irun == 0)
            pEE[cellid] = weight;
          else {
            double dr;
            if (method == 0) {
              dr = weight - pEE[cellid];
            }  // diff
            else {
              if (pEE[cellid] != 0.)
                dr = weight / pEE[cellid];
              else {
                if (weight == 0.)
                  dr = 1.;
                else
                  dr = 9999.;  // use a large value
              }
            }  // ratio
            if (dr < pEEmin)
              pEEmin = dr;
            if (dr > pEEmax)
              pEEmax = dr;
            if (iz == 1)
              endc_p->Fill(ix, iy, dr);
            else
              endc_m->Fill(ix, iy, dr);
          }

        }  // validDetId
}

inline void fillTableWithSummary(TH2F*& align,
                                 std::string title,
                                 const float& mean_x_EB,
                                 const float& rms_EB,
                                 const int& num_x_EB,
                                 const float& mean_x_EE,
                                 const float& rms_EE,
                                 const int& num_x_EE) {
  int NbRows = 2;
  align = new TH2F(title.c_str(), "EB/EE      mean_x      rms        num_x", 4, 0, 4, NbRows, 0, NbRows);

  double row = NbRows - 0.5;

  align->Fill(0.5, row, 1);
  align->Fill(1.5, row, mean_x_EB);
  align->Fill(2.5, row, rms_EB);
  align->Fill(3.5, row, num_x_EB);

  row--;

  align->Fill(0.5, row, 2);
  align->Fill(1.5, row, mean_x_EE);
  align->Fill(2.5, row, rms_EE);
  align->Fill(3.5, row, num_x_EE);

  align->GetXaxis()->SetTickLength(0.);
  align->GetXaxis()->SetLabelSize(0.);
  align->GetYaxis()->SetTickLength(0.);
  align->GetYaxis()->SetLabelSize(0.);
}