HcalObjRepresent Namespace Reference

Classes
class	ADataRepr
class	HcalDataContainer

Functions
int	CalcEtaBin (int subdet, int ieta, int depth)
int	CalcIeta (int subdet, int eta, int depth)
int	CalcIeta (int eta, int depth)
void	drawTable (int nbRows, int nbColumns)
void	Fill (HcalDetId &id, double val, std::vector< TH2F > &depth)
void	fillOneGain (std::vector< TH2F > &graphData, HcalGains::tAllContWithNames &allContainers, std::string name, int id, std::string units="")
void	FillUnphysicalHEHFBins (std::vector< TH2F > &hh)
void	FillUnphysicalHEHFBins (TH2F &hh)
std::string	getBitsSummary (uint32_t bits, std::string statusBitArray[], short unsigned int bitMap[])
std::vector< std::string >	HcalEtaPhiHistNames ()
std::string	IntToBinary (unsigned int number)
const bool	isBitSet (unsigned int bitnumber, unsigned int status)
bool	isHB (int etabin, int depth)
bool	isHE (int etabin, int depth)
bool	isHF (int etabin, int depth)
bool	isHO (int etabin, int depth)
bool	isSiPM (int ieta, int iphi, int depth)
void	Reset (std::vector< TH2F > &depth)
std::string	SciNotatStr (float num)
void	setBinLabels (std::vector< TH2F > &depth)
void	SetEtaPhiLabels (TH2F &h)
void	setup (std::vector< TH2F > &depth, std::string name, std::string units="")
bool	validDetId (HcalSubdetector sd, int ies, int ip, int dp)

Variables
const int	binmapd2 []
const int	binmapd3 [] = {-28, -27, -9999, -16, -9999, 16, -9999, 27, 28}

Function Documentation

int HcalObjRepresent::CalcEtaBin ( int  subdet, int  ieta, int  depth  )

inline

Definition at line 621 of file HcalObjRepresent.h.

References funct::abs(), HcalEndcap, HcalForward, and HcalOuter.

                                                         {
    // This takes the eta value from a subdetector and return an eta counter value as used by eta-phi array
    // (ieta=-41 corresponds to bin 0, +41 to bin 85 -- there are two offsets to deal with the overlap at |ieta|=29).
    // For HO, ieta = -15 corresponds to bin 0, and ieta=15 is bin 30
    // For HE depth 3, things are more complicated, but feeding the ieta value will give back the corresponding counter eta value

    // The CalcEtaBin value is the value as used within our array counters, and thus starts at 0.
    // If you are using it with getBinContent or setBinContent, you will need to add +1 to the result of this function

    int etabin = -9999;  // default invalid value

    if (depth == 1) {
      // Depth 1 is fairly straightforward -- just shift HF-, HF+ by -/+1
      etabin = ieta + 42;
      if (subdet == HcalForward) {
        ieta < 0 ? etabin-- : etabin++;
      }
    }

    else if (depth == 2) {
      // Depth 2 is more complicated, given that there are no cells in the range |ieta|<15
      if (ieta < -14) {
        etabin = ieta + 42;
        if (subdet == HcalForward)
          etabin--;
      } else if (ieta > 14) {
        etabin = ieta + 14;
        if (subdet == HcalForward)
          etabin++;
      }

    }
    // HO is also straightforward; a simple offset to the ieta value is applied
    else if (subdet == HcalOuter && abs(ieta) < 16)
      etabin = ieta + 15;
    else if (subdet == HcalEndcap) {
      // HE depth 3 has spotty coverage; hard-code the bin response
      if (depth == 3) {
        if (ieta == -28)
          etabin = 0;
        else if (ieta == -27)
          etabin = 1;
        else if (ieta == -16)
          etabin = 3;
        else if (ieta == 16)
          etabin = 5;
        else if (ieta == 27)
          etabin = 7;
        else if (ieta == 28)
          etabin = 8;
      }
    }
    return etabin;
  }

int HcalObjRepresent::CalcIeta ( int  subdet, int  eta, int  depth  )

inline

Definition at line 676 of file HcalObjRepresent.h.

References funct::abs(), binmapd2, binmapd3, PVValHelper::eta, HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

Referenced by FillUnphysicalHEHFBins(), isHB(), isHE(), isHF(), and isHO().

                                                      {
    // This function returns the 'true' ieta value given subdet, eta, and depth
    // Here 'eta' is the index from our arrays (it starts at 0);
    // remember that histogram bins start with bin 1, so there's an offset of 1
    // to consider if using getBinContent(eta,phi)

    // eta runs from 0...X  (X depends on depth)
    int ieta = -9999;  // default value is nonsensical
    if (subdet == HcalBarrel) {
      if (depth == 1) {
        ieta = eta - 42;
        if (ieta == 0)
          return -9999;
        return ieta;
      } else if (depth == 2) {
        ieta = binmapd2[eta];
        if (ieta == 0)
          return -9999;
        if (ieta == 17 || ieta == -17)
          return -9999;  // no depth 2 cells at |ieta| = 17
        return ieta;
      } else
        return -9999;  // non-physical value
    } else if (subdet == HcalForward) {
      if (depth == 1) {
        ieta = eta - 42;
        if (eta < 13)
          ieta++;
        else if (eta > 71)
          ieta--;
        else
          return -9999;  // if outside forward range, return dummy
        return ieta;
      } else if (depth == 2) {
        ieta = binmapd2[eta];  // special map for depth 2
        if (ieta <= -30)
          ieta++;
        else if (ieta >= 30)
          ieta--;
        else
          return -9999;
        return ieta;
      } else
        return -9999;
    }

    else if (subdet == HcalEndcap) {
      if (depth == 1)
        ieta = eta - 42;
      else if (depth == 2) {
        ieta = binmapd2[eta];
        if (abs(ieta) > 29 || abs(ieta) < 18)
          return -9999;  // outside HE
        if (ieta == 0)
          return -9999;
        return ieta;
      } else if (depth == 3) {
        if (eta < 0 || eta > 8)
          return -9999;
        else
          ieta = binmapd3[eta];  // special map for depth 3
        if (ieta == 0)
          return -9999;
        return ieta;
      } else
        return -9999;
    }  // HcalEndcap
    else if (subdet == HcalOuter) {
      if (depth != 4)
        return -9999;
      else {
        ieta = eta - 15;  // bin 0 is ieta=-15, all bins increment normally from there
        if (abs(ieta) > 15)
          return -9999;
        if (ieta == 0)
          return -9999;
        return ieta;
      }
    }  // HcalOuter
    if (ieta == 0)
      return -9999;
    return ieta;
  }

int HcalObjRepresent::CalcIeta ( int  eta, int  depth  )

inline

Definition at line 760 of file HcalObjRepresent.h.

References funct::abs(), binmapd2, binmapd3, and PVValHelper::eta.

                                          {
    // This version of CalcIeta does the same as the function above,
    // but does not require that 'subdet' be specified.

    // returns ieta value give an eta counter.
    // eta runs from 0...X  (X depends on depth)
    int ieta = -9999;
    if (eta < 0)
      return ieta;
    if (depth == 1) {
      ieta =
          eta -
          42;  // default shift: bin 0 corresponds to a histogram ieta of -42 (which is offset by 1 from true HF value of -41)
      if (eta < 13)
        ieta++;
      else if (eta > 71)
        ieta--;
      if (ieta == 0)
        ieta = -9999;
      return ieta;
    } else if (depth == 2) {
      if (eta > 57)
        return -9999;
      else {
        ieta = binmapd2[eta];
        if (ieta == -9999)
          return ieta;
        if (ieta == 0)
          return -9999;
        if (ieta == 17 || ieta == -17)
          return -9999;  // no depth 2 cells at |ieta| = 17
        else if (ieta <= -30)
          ieta++;
        else if (ieta >= 30)
          ieta--;
        return ieta;
      }
    } else if (depth == 3) {
      if (eta > 8)
        return -9999;
      else
        ieta = binmapd3[eta];
      if (ieta == 0)
        return -9999;
      return ieta;
    } else if (depth == 4) {
      ieta = eta - 15;  // bin 0 is ieta=-15, all bins increment normally from there
      if (abs(ieta) > 15)
        return -9999;
      if (ieta == 0)
        return -9999;
      return ieta;
    }
    return ieta;  // avoids compilation warning
  }

void HcalObjRepresent::drawTable ( int  nbRows, int  nbColumns  )

inline

Definition at line 484 of file HcalObjRepresent.h.

References mps_fire::i, and cmsLHEtoEOSManager::l.

                                                   {
    TLine* l = new TLine;
    l->SetLineWidth(1);
    for (int i = 1; i < nbRows; i++) {
      double y = (double)i;
      l = new TLine(0., y, nbColumns, y);
      l->Draw();
    }

    for (int i = 1; i < nbColumns; i++) {
      double x = (double)i;
      double y = (double)nbRows;
      l = new TLine(x, 0., x, y);
      l->Draw();
    }
  }

void HcalObjRepresent::Fill ( HcalDetId &  id, double  val, std::vector< TH2F > &  depth  )

inline

Definition at line 1053 of file HcalObjRepresent.h.

References HcalForward.

Referenced by SiStripGainCosmicCalculator::algoAnalyze(), L1TStage2EMTF::analyze(), RPCFEDIntegrity::analyze(), L1TPUM::analyze(), PrimaryVertexMonitor::analyze(), DTVDriftSegmentCalibration::analyze(), DTOccupancyEfficiency::analyze(), MuonTiming::analyze(), L1TDEMON::analyze(), DTPreCalibrationTask::analyze(), TrackingRecoMaterialAnalyser::analyze(), DTTTrigOffsetCalibration::analyze(), MuonRecoAnalyzer::analyze(), SiPixelRecHitsValid::analyze(), SiStripMonitorDigi::analyze(), CTPPSPixelDQMSource::analyze(), SiPixelHLTSource::analyze(), ECALpedestalPCLworker::analyze(), tadqm::TrackAnalyzer::analyze(), HGCalHitCalibration::analyze(), SiPixelTrackResidualSource::analyze(), SiStripMonitorCluster::analyze(), HcalLutAnalyzer::analyze(), DTChamberEfficiency::analyze(), PhiSymmetryCalibration::analyze(), APVCyclePhaseMonitor::analyze(), SiStripMonitorPedestals::analyze(), DTRunConditionVar::analyze(), HLTJetMETValidation::analyze(), EventTimeDistribution::analyze(), HcalRecHitsAnalyzer::analyze(), PhotonValidator::analyze(), SiStripQualityHistory::analyze(), SiPixelQualityHistory::analyze(), cms::Analyzer_minbias::analyze(), TrackToTrackComparisonHists::analyze(), DQMPFCandidateAnalyzer::analyze(), HLTMuonPlotter::analyze(), CommonModeAnalyzer::analyze(), CTPPSDiamondDQMSource::analyze(), analyzer::SiPixelLorentzAngle::analyze(), GeneralPurposeTrackAnalyzer::analyze(), DMRChecker::analyze(), DTNoiseComputation::beginRun(), plotTTrigFromHistos::compareDiff(), plotT0FromHistos::compareDiff(), CommissioningHistograms::copyCustomInformation(), SiStripCorrelateBadStripAndNoise::correlateWithNoise(), SiStripPlotGain::DoAnalysis(), SiStripCorrelateNoise::DoAnalysis(), PPSTimingCalibrationPCLWorker::dqmAnalyze(), PPSDiamondSampicTimingCalibrationPCLWorker::dqmAnalyze(), DTResolutionAnalysisTest::dqmEndJob(), DTRunConditionVarClient::dqmEndJob(), InvRingCalib::duringLoop(), EcalEleCalibLooper::duringLoop(), ESSummaryClient::endLumiAnalyze(), DigiInvestigatorHistogramMaker::fill(), BenchmarkTree::Fill(), DigiVertexCorrHistogramMaker::fill(), DigiLumiCorrHistogramMaker::fill(), DigiBXCorrHistogramMaker< T >::fill(), SiPixelTrackResidualModule::fill(), SiPixelRawDataErrorModule::fill(), Histos::fill(), SiPixelClusterModule::fill(), SiPixelDigiModule::fill(), XHistogram::fill(), DD4hep_XHistogram::fill(), Book::fill(), CSCValHists::fill2DHistByEvent(), TrackToTrackComparisonHists::fill_generic_tracks_histos(), TrackToTrackComparisonHists::fill_matching_tracks_histos(), SiPixelRecHitsValid::fillBarrel(), SiStripMonitorPedestals::fillCondDBMEs(), BPhysicsValidation::FillDaughters(), SiPixelRawDataErrorModule::fillFED(), tadqm::TrackAnalyzer::fillHistosForEfficiencyFromHitPatter(), SiStripMonitorCluster::fillModuleMEs(), METAnalyzer::fillMonitorElement(), GlobalDigisAnalyzer::fillMuon(), LRHelpFunctions::fillToBackgroundHists(), edm::RootOutputTree::fillTTree(), SiStripGainCosmicCalculator::getNewObject(), overlapValidationPlot::hist(), TMultiDimFet::MakeCoefficients(), TMultiDimFet::MakeNormalized(), CosmicGenFilterHelix::monitorEnd(), CosmicGenFilterHelix::monitorStart(), DTSegmentAnalysisTest::performClientDiagnostic(), bigModule::plot(), subModule::plot(), SiStripGainFromCalibTree::processEvent(), DTDataIntegrityTask::processuROS(), ecaldqm::MESetDet0D::reset(), ecaldqm::MESet::resetAll(), and PrimaryVertexMonitor::vertexPlots().

                                                                        {
    // If in HF, need to shift by 1 bin (-1 bin lower in -HF, +1 bin higher in +HF)
    if (id.subdet() == HcalForward)
      depth[id.depth() - 1].Fill(id.ieta() < 0 ? id.ieta() - 1 : id.ieta() + 1, id.iphi(), val);
    else
      depth[id.depth() - 1].Fill(id.ieta(), id.iphi(), val);
  }

void HcalObjRepresent::fillOneGain ( std::vector< TH2F > &  graphData, HcalGains::tAllContWithNames &  allContainers, std::string  name, int  id, std::string  units = ""  )

inline

Definition at line 1135 of file HcalObjRepresent.h.

References ztail::d, HcalDetId::depth(), HLT_FULL_cff::depth, HcalForward, mps_fire::i, gpuClustering::id, HcalDetId::ieta(), HcalDetId::iphi(), setup(), and HcalDetId::subdet().

                                                {
    setup(graphData, name);

    std::stringstream x;
    // Change the titles of each individual histogram
    for (unsigned int d = 0; d < graphData.size(); ++d) {
      graphData[d].Reset();
      x << "Gain " << id << " for HCAL depth " << d + 1;

      //BUG CAN BE HERE:
      //if (ChannelStatus->depth[d])
      graphData[d].SetTitle(
          x.str()
              .c_str());  // replace "setTitle" with "SetTitle", since you are using TH2F objects instead of MonitorElements
      x.str("");
    }

    HcalDetId hcal_id;
    int ieta, depth, iphi;

    //main loop
    // get all containers with names
    //HcalGains::tAllContWithNames allContainers = object().getAllContainers();

    //ITERATORS AND VALUES:
    HcalGains::tAllContWithNames::const_iterator iter;
    std::vector<HcalGain>::const_iterator contIter;
    float gain = 0.0;

    //Run trough given id gain:
    int i = id;

    //run trough all pair containers
    for (iter = allContainers.begin(); iter != allContainers.end(); ++iter) {
      //Run trough all values:
      for (contIter = (*iter).second.begin(); contIter != (*iter).second.end(); ++contIter) {
        hcal_id = HcalDetId((uint32_t)(*contIter).rawId());

        depth = hcal_id.depth();
        if (depth < 1 || depth > 4)
          continue;

        ieta = hcal_id.ieta();
        iphi = hcal_id.iphi();

        if (hcal_id.subdet() == HcalForward)
          ieta > 0 ? ++ieta : --ieta;

        //GET VALUE:
        gain = (*contIter).getValue(i);
        //logstatus = log2(1.*channelBits)+1;

        //FILLING GOES HERE:
        graphData[depth - 1].Fill(ieta, iphi, gain);

        //FOR DEBUGGING:
        //std::cout << "ieta: " << ieta << "; iphi: " << iphi << "; logstatus: " << logstatus << "; channelBits: " << channelBits<< std::endl;
      }
    }
  }

void HcalObjRepresent::FillUnphysicalHEHFBins ( std::vector< TH2F > &  hh )

inline

Definition at line 979 of file HcalObjRepresent.h.

References funct::abs(), CalcIeta(), ztail::d, and PVValHelper::eta.

Referenced by HcalObjRepresent::ADataRepr::fillOneGain().

                                                          {
    int ieta = 0;
    int iphi = 0;
    // First 2 depths have 5-10-20 degree corrections
    for (unsigned int d = 0; d < 3; ++d) {
      //BUG CAN BE HERE:
      //if (hh[d] != 0) continue;

      for (int eta = 0; eta < hh[d].GetNbinsX(); ++eta) {
        ieta = CalcIeta(eta, d + 1);
        if (ieta == -9999 || abs(ieta) < 21)
          continue;
        for (int phi = 0; phi < hh[d].GetNbinsY(); ++phi) {
          iphi = phi + 1;
          if (iphi % 2 == 1 && abs(ieta) < 40 && iphi < 73) {
            hh[d].SetBinContent(eta + 1, iphi + 1, hh[d].GetBinContent(eta + 1, iphi));
          }
          // last two eta strips span 20 degrees in phi
          // Fill the phi cell above iphi, and the 2 below it
          else if (abs(ieta) > 39 && iphi % 4 == 3 && iphi < 73) {
            //ieta=40, iphi=3 covers iphi 3,4,5,6
            hh[d].SetBinContent(eta + 1, (iphi) % 72 + 1, hh[d].GetBinContent(eta + 1, iphi));
            hh[d].SetBinContent(eta + 1, (iphi + 1) % 72 + 1, hh[d].GetBinContent(eta + 1, iphi));
            hh[d].SetBinContent(eta + 1, (iphi + 2) % 72 + 1, hh[d].GetBinContent(eta + 1, iphi));
          }
        }  // for (int phi...)
      }    // for (int eta...)
    }      // for (int d=0;...)
    // no corrections needed for HO (depth 4)
    return;
  }  // FillUnphysicalHEHFBins(MonitorElement* hh)

void HcalObjRepresent::FillUnphysicalHEHFBins ( TH2F &  hh )

inline

Definition at line 1012 of file HcalObjRepresent.h.

References funct::abs(), CalcIeta(), PVValHelper::eta, and defaults_cfi::etabins.

                                               {
    // Fills unphysical HE/HF bins for Summary Histogram
    // Summary Histogram is binned with the same binning as the Depth 1 EtaPhiHists

    //CAN BE BUG HERE:
    //if (hh==0) return;

    int ieta = 0;
    int iphi = 0;
    int etabins = hh.GetNbinsX();
    int phibins = hh.GetNbinsY();
    float binval = 0;
    for (int eta = 0; eta < etabins; ++eta)  // loop over eta bins
    {
      ieta = CalcIeta(eta, 1);
      if (ieta == -9999 || abs(ieta) < 21)
        continue;  // ignore etas that don't exist, or that have 5 degree phi binning

      for (int phi = 0; phi < phibins; ++phi) {
        iphi = phi + 1;
        if (iphi % 2 == 1 && abs(ieta) < 40 && iphi < 73)  // 10 degree phi binning condition
        {
          binval = hh.GetBinContent(eta + 1, iphi);
          hh.SetBinContent(eta + 1, iphi + 1, binval);
        }                                                       // if (iphi%2==1...)
        else if (abs(ieta) > 39 && iphi % 4 == 3 && iphi < 73)  // 20 degree phi binning condition
        {
          // Set last two eta strips where each cell spans 20 degrees in phi
          // Set next phi cell above iphi, and 2 cells below the actual cell
          hh.SetBinContent(eta + 1, (iphi) % 72 + 1, hh.GetBinContent(eta + 1, iphi));
          hh.SetBinContent(eta + 1, (iphi + 1) % 72 + 1, hh.GetBinContent(eta + 1, iphi));
          hh.SetBinContent(eta + 1, (iphi + 2) % 72 + 1, hh.GetBinContent(eta + 1, iphi));
        }  // else if (abs(ieta)>39 ...)
      }    // for (int phi=0;phi<72;++phi)

    }  // for (int eta=0; eta< (etaBins_-2);++eta)

    return;
  }  // FillUnphysicalHEHFBins(std::vector<MonitorElement*> &hh)

std::string HcalObjRepresent::getBitsSummary ( uint32_t  bits, std::string  statusBitArray[], short unsigned int  bitMap[]  )

inline

Definition at line 536 of file HcalObjRepresent.h.

References mps_fire::i, isBitSet(), and contentValuesCheck::ss.

                                                                                                          {
    std::stringstream ss;
    for (unsigned int i = 0; i < 9; ++i) {
      if (isBitSet(bitMap[i], bits)) {
        ss << "[" << bitMap[i] << "]" << statusBitArray[bitMap[i]] << "; ";
      }
    }
    ss << std::endl;
    return ss.str();
  }

std::vector<std::string> HcalObjRepresent::HcalEtaPhiHistNames ( )

inline

Definition at line 818 of file HcalObjRepresent.h.

References mergeVDriftHistosByStation::name.

                                                    {
    std::vector<std::string> name;
    name.push_back("HB HE HF Depth 1 ");
    name.push_back("HB HE HF Depth 2 ");
    name.push_back("HE Depth 3 ");
    name.push_back("HO Depth 4 ");
    return name;
  }

std::string HcalObjRepresent::IntToBinary ( unsigned int  number )

inline

Definition at line 516 of file HcalObjRepresent.h.

References mps_fire::i, and contentValuesCheck::ss.

                                                    {
    std::stringstream ss;
    unsigned int mask = 1 << 31;
    for (unsigned short int i = 0; i < 32; ++i) {
      //if (!(i % 4))
      //    ss << "_";
      if (mask & number)
        ss << "1";
      else
        ss << "0";
      mask = mask >> 1;
    }
    return ss.str();
  }

const bool HcalObjRepresent::isBitSet ( unsigned int  bitnumber, unsigned int  status  )

inline

Definition at line 531 of file HcalObjRepresent.h.

Referenced by getBitsSummary().

                                                                          {
    unsigned int statadd = 0x1 << (bitnumber);
    return (status & statadd) ? (true) : (false);
  }

bool HcalObjRepresent::isHB ( int  etabin, int  depth  )

inline

Definition at line 827 of file HcalObjRepresent.h.

References funct::abs(), and CalcIeta().

Referenced by HcalGeometry::alignmentTransformIndexLocal().

                                          {
    if (depth > 2)
      return false;
    else if (depth < 1)
      return false;
    else {
      int ieta = CalcIeta(etabin, depth);
      if (ieta == -9999)
        return false;
      if (depth == 1) {
        if (abs(ieta) <= 16)
          return true;
        else
          return false;
      } else if (depth == 2) {
        if (abs(ieta) == 15 || abs(ieta) == 16)
          return true;
        else
          return false;
      }
    }
    return false;
  }

bool HcalObjRepresent::isHE ( int  etabin, int  depth  )

inline

Definition at line 851 of file HcalObjRepresent.h.

References funct::abs(), and CalcIeta().

Referenced by HcalGeometry::alignmentTransformIndexLocal(), and l1t::HGCalClusterT< l1t::HGCalCluster >::hOverE().

                                          {
    if (depth > 3)
      return false;
    else if (depth < 1)
      return false;
    else {
      int ieta = CalcIeta(etabin, depth);
      if (ieta == -9999)
        return false;
      if (depth == 1) {
        if (abs(ieta) >= 17 && abs(ieta) <= 28)
          return true;
        if (ieta == -29 && etabin == 13)
          return true;  // HE -29
        if (ieta == 29 && etabin == 71)
          return true;  // HE +29
      } else if (depth == 2) {
        if (abs(ieta) >= 17 && abs(ieta) <= 28)
          return true;
        if (ieta == -29 && etabin == 13)
          return true;  // HE -29
        if (ieta == 29 && etabin == 43)
          return true;  // HE +29
      } else if (depth == 3)
        return true;
    }
    return false;
  }

bool HcalObjRepresent::isHF ( int  etabin, int  depth  )

inline

Definition at line 880 of file HcalObjRepresent.h.

References funct::abs(), and CalcIeta().

Referenced by HcalGeometry::alignmentTransformIndexLocal(), hcaldqm::hashfunctions::hash_FED(), and reco::CentralityProducer::produce().

                                          {
    if (depth > 2)
      return false;
    else if (depth < 1)
      return false;
    else {
      int ieta = CalcIeta(etabin, depth);
      if (ieta == -9999)
        return false;
      if (depth == 1) {
        if (ieta == -29 && etabin == 13)
          return false;  // HE -29
        else if (ieta == 29 && etabin == 71)
          return false;  // HE +29
        else if (abs(ieta) >= 29)
          return true;
      } else if (depth == 2) {
        if (ieta == -29 && etabin == 13)
          return false;  // HE -29
        else if (ieta == 29 && etabin == 43)
          return false;  // HE +29
        else if (abs(ieta) >= 29)
          return true;
      }
    }
    return false;
  }

bool HcalObjRepresent::isHO ( int  etabin, int  depth  )

inline

Definition at line 908 of file HcalObjRepresent.h.

References CalcIeta().

Referenced by hcaldqm::hashfunctions::hash_FED().

                                          {
    if (depth != 4)
      return false;
    int ieta = CalcIeta(etabin, depth);
    if (ieta != -9999)
      return true;
    return false;
  }

bool HcalObjRepresent::isSiPM ( int  ieta, int  iphi, int  depth  )

inline

Definition at line 919 of file HcalObjRepresent.h.

                                                    {
    if (depth != 4)
      return false;
    // HOP1
    if (ieta >= 5 && ieta <= 10 && iphi >= 47 && iphi <= 58)
      return true;
    // HOP2
    if (ieta >= 11 && ieta <= 15 && iphi >= 59 && iphi <= 70)
      return true;
    return false;
  }  // bool isSiPM

void HcalObjRepresent::Reset ( std::vector< TH2F > &  depth )

inline

Definition at line 1061 of file HcalObjRepresent.h.

References ztail::d.

Referenced by plotTTrigFromHistos::compareDiff(), plotT0FromHistos::compareDiff(), TH2Analyzer::Eval(), and RootSig2XTReqHandler::XTReqTimer::FireAway().

                                            {
    for (unsigned int d = 0; d < depth.size(); d++)
      //BUG CAN BE HERE:
      //if(depth[d])
      depth[d].Reset();
  }  // void Reset(void)

std::string HcalObjRepresent::SciNotatStr ( float  num )

inline

Definition at line 501 of file HcalObjRepresent.h.

References pileupDistInMC::num, and AlCaHLTBitMon_QueryRunRegistry::string.

                                          {
    // Create an output string stream
    std::ostringstream streamObj2;

    if (num == -1)
      return "NOT FOUND";

    //Add double to stream
    streamObj2 << num;
    // Get string from output string stream
    std::string strObj2 = streamObj2.str();

    return strObj2;
  }

void HcalObjRepresent::setBinLabels ( std::vector< TH2F > &  depth )

inline

Definition at line 548 of file HcalObjRepresent.h.

References mps_fire::i, and label.

Referenced by MTVHistoProducerAlgoForTracker::bookRecoHistos(), and setup().

                                                   {
    // Set labels for all depth histograms
    for (unsigned int i = 0; i < depth.size(); ++i) {
      depth[i].SetXTitle("i#eta");
      depth[i].SetYTitle("i#phi");
    }

    std::stringstream label;

    // set label on every other bin
    for (int i = -41; i <= -29; i = i + 2) {
      label << i;
      depth[0].GetXaxis()->SetBinLabel(i + 42, label.str().c_str());
      depth[1].GetXaxis()->SetBinLabel(i + 42, label.str().c_str());
      label.str("");
    }
    depth[0].GetXaxis()->SetBinLabel(14, "-29HE");
    depth[1].GetXaxis()->SetBinLabel(14, "-29HE");

    // offset by one for HE
    for (int i = -27; i <= 27; i = i + 2) {
      label << i;
      depth[0].GetXaxis()->SetBinLabel(i + 43, label.str().c_str());
      label.str("");
    }
    depth[0].GetXaxis()->SetBinLabel(72, "29HE");
    for (int i = 29; i <= 41; i = i + 2) {
      label << i;
      depth[0].GetXaxis()->SetBinLabel(i + 44, label.str().c_str());
      label.str("");
    }
    for (int i = 16; i <= 28; i = i + 2) {
      label << i - 43;
      depth[1].GetXaxis()->SetBinLabel(i, label.str().c_str());
      label.str("");
    }
    depth[1].GetXaxis()->SetBinLabel(29, "NULL");
    for (int i = 15; i <= 27; i = i + 2) {
      label << i;
      depth[1].GetXaxis()->SetBinLabel(i + 15, label.str().c_str());
      label.str("");
    }

    depth[1].GetXaxis()->SetBinLabel(44, "29HE");
    for (int i = 29; i <= 41; i = i + 2) {
      label << i;
      depth[1].GetXaxis()->SetBinLabel(i + 16, label.str().c_str());
      label.str("");
    }

    // HE depth 3 labels;
    depth[2].GetXaxis()->SetBinLabel(1, "-28");
    depth[2].GetXaxis()->SetBinLabel(2, "-27");
    depth[2].GetXaxis()->SetBinLabel(3, "Null");
    depth[2].GetXaxis()->SetBinLabel(4, "-16");
    depth[2].GetXaxis()->SetBinLabel(5, "Null");
    depth[2].GetXaxis()->SetBinLabel(6, "16");
    depth[2].GetXaxis()->SetBinLabel(7, "Null");
    depth[2].GetXaxis()->SetBinLabel(8, "27");
    depth[2].GetXaxis()->SetBinLabel(9, "28");
  }

void HcalObjRepresent::SetEtaPhiLabels ( TH2F &  h )

inline

Definition at line 954 of file HcalObjRepresent.h.

References mps_fire::i, and label.

                                       {
    std::stringstream label;
    for (int i = -41; i <= -29; i = i + 2) {
      label << i;
      h.GetXaxis()->SetBinLabel(i + 42, label.str().c_str());
      label.str("");
    }
    h.GetXaxis()->SetBinLabel(14, "-29HE");

    // offset by one for HE
    for (int i = -27; i <= 27; i = i + 2) {
      label << i;
      h.GetXaxis()->SetBinLabel(i + 43, label.str().c_str());
      label.str("");
    }
    h.GetXaxis()->SetBinLabel(72, "29HE");
    for (int i = 29; i <= 41; i = i + 2) {
      label << i;
      h.GetXaxis()->SetBinLabel(i + 44, label.str().c_str());
      label.str("");
    }
    return;
  }

void HcalObjRepresent::setup ( std::vector< TH2F > &  depth, std::string  name, std::string  units = ""  )

inline

Definition at line 1068 of file HcalObjRepresent.h.

References mps_fire::i, setBinLabels(), AlCaHLTBitMon_QueryRunRegistry::string, units(), and fw3dlego::xbins.

Referenced by fillOneGain().

                                                                                  {
    std::string unittitle, unitname;
    if (units.empty()) {
      unitname = units;
      unittitle = "No Units";
    } else {
      unitname = " " + units;
      unittitle = units;
    }

    // Push back depth plots
    //1. create first plot
    depth.push_back(TH2F(("HB HE HF Depth 1 " + name + unitname).c_str(),
                         (name + " Depth 1 -- HB HE HF (" + unittitle + ")").c_str(),
                         85,
                         -42.5,
                         42.5,
                         72,
                         0.5,
                         72.5));

    //2.1 prepare second plot
    float ybins[73];
    for (int i = 0; i <= 72; i++)
      ybins[i] = (float)(i + 0.5);
    float xbinsd2[] = {-42.5, -41.5, -40.5, -39.5, -38.5, -37.5, -36.5, -35.5, -34.5, -33.5, -32.5, -31.5,
                       -30.5, -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5, -20.5, -19.5,
                       -18.5, -17.5, -16.5, -15.5, -14.5, 14.5,  15.5,  16.5,  17.5,  18.5,  19.5,  20.5,
                       21.5,  22.5,  23.5,  24.5,  25.5,  26.5,  27.5,  28.5,  29.5,  30.5,  31.5,  32.5,
                       33.5,  34.5,  35.5,  36.5,  37.5,  38.5,  39.5,  40.5,  41.5,  42.5};

    //2.2 create second plot
    depth.push_back(TH2F(("HB HE HF Depth 2 " + name + unitname).c_str(),
                         (name + " Depth 2 -- HB HE HF (" + unittitle + ")").c_str(),
                         57,
                         xbinsd2,
                         72,
                         ybins));

    //3.1 Set up variable-sized bins for HE depth 3 (MonitorElement also requires phi bins to be entered in array format)
    float xbins[] = {-28.5, -27.5, -26.5, -16.5, -15.5, 15.5, 16.5, 26.5, 27.5, 28.5};
    //3.2
    depth.push_back(TH2F(("HE Depth 3 " + name + unitname).c_str(),
                         (name + " Depth 3 -- HE (" + unittitle + ")").c_str(),
                         // Use variable-sized eta bins
                         9,
                         xbins,
                         72,
                         ybins));

    //4.1 HO bins are fixed width, but cover a smaller eta range (-15 -> 15)
    depth.push_back(TH2F(("HO Depth 4 " + name + unitname).c_str(),
                         (name + " Depth 4 -- HO (" + unittitle + ")").c_str(),
                         31,
                         -15.5,
                         15.5,
                         72,
                         0.5,
                         72.5));

    for (unsigned int i = 0; i < depth.size(); ++i)
      depth[i].Draw("colz");

    setBinLabels(depth);  // set axis titles, special bins
  }

bool HcalObjRepresent::validDetId ( HcalSubdetector  sd, int  ies, int  ip, int  dp  )

inline

Definition at line 933 of file HcalObjRepresent.h.

References funct::abs(), HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

                                                                      {
    // inputs are (subdetector, ieta, iphi, depth)
    // stolen from latest version of DataFormats/HcalDetId/src/HcalDetId.cc (not yet available in CMSSW_2_1_9)

    const int ie(abs(ies));

    return (
        (ip >= 1) && (ip <= 72) && (dp >= 1) && (ie >= 1) &&
        (((sd == HcalBarrel) && (((ie <= 14) && (dp == 1)) || (((ie == 15) || (ie == 16)) && (dp <= 2)))) ||
         ((sd == HcalEndcap) &&
          (((ie == 16) && (dp == 3)) || ((ie == 17) && (dp == 1)) || ((ie >= 18) && (ie <= 20) && (dp <= 2)) ||
           ((ie >= 21) && (ie <= 26) && (dp <= 2) && (ip % 2 == 1)) ||
           ((ie >= 27) && (ie <= 28) && (dp <= 3) && (ip % 2 == 1)) || ((ie == 29) && (dp <= 2) && (ip % 2 == 1)))) ||
         ((sd == HcalOuter) && (ie <= 15) && (dp == 4)) ||
         ((sd == HcalForward) && (dp <= 2) &&
          (((ie >= 29) && (ie <= 39) && (ip % 2 == 1)) || ((ie >= 40) && (ie <= 41) && (ip % 4 == 3))))));

  }  // bool validDetId(HcalSubdetector sd, int ies, int ip, int dp)

Variable Documentation

const int HcalObjRepresent::binmapd2[]

Initial value:

= {-42, -41, -40, -39, -38, -37, -36, -35, -34, -33, -32, -31, -30, -29,   -28,
                          -27, -26, -25, -24, -23, -22, -21, -20, -19, -18, -17, -16, -15, -9999, 15,
                          16,  17,  18,  19,  20,  21,  22,  23,  24,  25,  26,  27,  28,  29,    30,
                          31,  32,  33,  34,  35,  36,  37,  38,  39,  40,  41,  42}

Definition at line 612 of file HcalObjRepresent.h.

Referenced by CalcIeta().

const int HcalObjRepresent::binmapd3[] = {-28, -27, -9999, -16, -9999, 16, -9999, 27, 28}

Definition at line 619 of file HcalObjRepresent.h.

Referenced by CalcIeta().