EcalFloatCondObjectContainerUtils.h File Reference

#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>

Go to the source code of this file.

Enumerations
enum	{ kEBChannels = 61200, kEEChannels = 14648 }
enum	{   MIN_IETA = 1, MIN_IPHI = 1, MAX_IETA = 85, MAX_IPHI = 360,   EBhistEtaMax = 171 }
enum	{   IX_MIN = 1, IY_MIN = 1, IX_MAX = 100, IY_MAX = 100,   EEhistXMax = 220 }

Functions
template<class floatCondObj >
void	fillEBMap_SingleIOV (std::shared_ptr< floatCondObj > payload, TH2F *&barrel)
template<class floatCondObj >
void	fillEBMap_TwoIOVs (std::shared_ptr< floatCondObj > payload, TH2F *&barrel, int irun, float pEB[], float &pEBmin, float &pEBmax, int method)
template<class floatCondObj >
void	fillEEMap_SingleIOV (std::shared_ptr< floatCondObj > payload, TH2F *&endc_m, TH2F *&endc_p)
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)
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)

Enumeration Type Documentation

anonymous enum

anonymous enum

Enumerator
kEBChannels
kEEChannels

Definition at line 9 of file EcalFloatCondObjectContainerUtils.h.

{ kEBChannels = 61200, kEEChannels = 14648 };

anonymous enum

anonymous enum

Enumerator
MIN_IETA
MIN_IPHI
MAX_IETA
MAX_IPHI
EBhistEtaMax

Definition at line 10 of file EcalFloatCondObjectContainerUtils.h.

     {
  MIN_IETA = 1,
  MIN_IPHI = 1,
  MAX_IETA = 85,
  MAX_IPHI = 360,
  EBhistEtaMax = 171
};  // barrel lower and upper bounds on eta and phi

anonymous enum

anonymous enum

Enumerator
IX_MIN
IY_MIN
IX_MAX
IY_MAX
EEhistXMax

Definition at line 17 of file EcalFloatCondObjectContainerUtils.h.

     {
  IX_MIN = 1,
  IY_MIN = 1,
  IX_MAX = 100,
  IY_MAX = 100,
  EEhistXMax = 220
};  // endcaps lower and upper bounds on x and y

Function Documentation

fillEBMap_SingleIOV()

template<class floatCondObj >

void fillEBMap_SingleIOV	(	std::shared_ptr< floatCondObj >	payload,
		TH2F *&	barrel
	)

Definition at line 26 of file EcalFloatCondObjectContainerUtils.h.

References Reference_intrackfit_cff::barrel, PVValHelper::eta, nano_mu_digi_cff::float, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, EBDetId::kSizeForDenseIndexing, EBDetId::MIN_HASH, jetsAK4_Puppi_cff::payload, PVValHelper::phi, and EBDetId::unhashIndex().

                                                                             {
  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
}

fillEBMap_TwoIOVs()

template<class floatCondObj >

void fillEBMap_TwoIOVs	(	std::shared_ptr< floatCondObj >	payload,
		TH2F *&	barrel,
		int	irun,
		float	pEB[],
		float &	pEBmin,
		float &	pEBmax,
		int	method
	)

Definition at line 69 of file EcalFloatCondObjectContainerUtils.h.

References Reference_intrackfit_cff::barrel, l1ctLayer1_cff::dr, PVValHelper::eta, nano_mu_digi_cff::float, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, EBDetId::kSizeForDenseIndexing, AlcaSiPixelAliHarvester0T_cff::method, EBDetId::MIN_HASH, jetsAK4_Puppi_cff::payload, PVValHelper::phi, EBDetId::unhashIndex(), and mps_merge::weight.

                                   {
  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
}

fillEEMap_SingleIOV()

template<class floatCondObj >

void fillEEMap_SingleIOV	(	std::shared_ptr< floatCondObj >	payload,
		TH2F *&	endc_m,
		TH2F *&	endc_p
	)

Definition at line 47 of file EcalFloatCondObjectContainerUtils.h.

References nano_mu_digi_cff::float, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), IX_MAX, IX_MIN, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), IY_MAX, IY_MIN, jetsAK4_Puppi_cff::payload, DetId::rawId(), EEDetId::validDetId(), and EEDetId::XYMODE.

                                                                                            {
  // 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
}

fillEEMap_TwoIOVs()

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
	)

Definition at line 117 of file EcalFloatCondObjectContainerUtils.h.

References l1ctLayer1_cff::dr, nano_mu_digi_cff::float, EEDetId::hashedIndex(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), IX_MAX, IX_MIN, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), IY_MAX, IY_MIN, AlcaSiPixelAliHarvester0T_cff::method, jetsAK4_Puppi_cff::payload, DetId::rawId(), EEDetId::validDetId(), mps_merge::weight, and EEDetId::XYMODE.

                                   {
  // 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
}

fillTableWithSummary()

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  )

inline

Definition at line 169 of file EcalFloatCondObjectContainerUtils.h.

References runGCPTkAlMap::title.

                                                      {
  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.);
}