Phase1PixelMaps Class Reference

#include <Phase1PixelMaps.h>

Public Member Functions
void	beautifyAllHistograms ()
void	book (const std::string &currentHistoName, const char *what, const char *zaxis)
void	bookBarrelHistograms (const std::string &currentHistoName, const char *what, const char *zaxis)
void	bookForwardHistograms (const std::string &currentHistoName, const char *what, const char *zaxis)
void	drawBarrelMaps (const std::string &currentHistoName, TCanvas &canvas, const char *drawOption=nullptr)
void	drawForwardMaps (const std::string &currentHistoName, TCanvas &canvas, const char *drawOption=nullptr)
void	drawSummaryMaps (const std::string &currentHistoName, TCanvas &canvas, const char *drawOption=nullptr)
void	fill (const std::string &currentHistoName, unsigned int id, double value)
void	fillBarrelBin (const std::string &currentHistoName, unsigned int id, double value)
void	fillForwardBin (const std::string &currentHistoName, unsigned int id, double value)
	Phase1PixelMaps (const char *option)
void	resetOption (const char *option)
void	setBarrelScale (const std::string &currentHistoName, std::pair< float, float > extrema)
void	setForwardScale (const std::string &currentHistoName, std::pair< float, float > extrema)
void	setNoRescale ()
	~Phase1PixelMaps ()=default

Private Member Functions
void	adjustCanvasMargins (TVirtualPad *pad, float top, float bottom, float left, float right)
void	bookBarrelBins (const std::string &currentHistoName)
void	bookForwardBins (const std::string &currentHistoName)
void	makeNicePlotStyle (TH1 *hist)
void	rescaleAllBarrel (const std::string &currentHistoName)
void	rescaleAllForward (const std::string &currentHistoName)
const indexedCorners	retrieveCorners (const std::vector< edm::FileInPath > &cornerFiles, const unsigned int reads)

Private Attributes
std::map< uint32_t, std::shared_ptr< TGraph > >	bins
std::map< uint32_t, std::shared_ptr< TGraph > >	binsSummary
bool	m_autorescale
std::vector< edm::FileInPath >	m_cornersBPIX
std::vector< edm::FileInPath >	m_cornersFPIX
std::pair< bool, bool >	m_isBooked
std::vector< std::string >	m_knownNames
Option_t *	m_option
TrackerTopology	m_trackerTopo
std::map< std::string, std::vector< std::shared_ptr< TH2Poly > > >	pxbTh2PolyBarrel
std::map< std::string, std::shared_ptr< TH2Poly > >	pxbTh2PolyBarrelSummary
std::map< std::string, std::vector< std::shared_ptr< TH2Poly > > >	pxfTh2PolyForward
std::map< std::string, std::shared_ptr< TH2Poly > >	pxfTh2PolyForwardSummary

Detailed Description

Definition at line 31 of file Phase1PixelMaps.h.

Constructor & Destructor Documentation

Phase1PixelMaps()

Phase1PixelMaps::Phase1PixelMaps ( const char *  option )

inline

Definition at line 33 of file Phase1PixelMaps.h.

References fileinputsource_cfi::option.

      : m_option{option},
        m_isBooked{std::make_pair(false, false)},
        m_knownNames{{}},
        m_trackerTopo{StandaloneTrackerTopology::fromTrackerParametersXMLFile(
            edm::FileInPath("Geometry/TrackerCommonData/data/PhaseI/trackerParameters.xml").fullPath())} {
    // set the rescale to true by default
    m_autorescale = true;

    // store the file in path for the corners (BPIX)
    for (unsigned int i = 1; i <= 4; i++) {
      m_cornersBPIX.push_back(edm::FileInPath(Form("DQM/SiStripMonitorClient/data/Geometry/vertices_barrel_%i", i)));
    }

    // store the file in path for the corners (BPIX)
    for (int j : {-3, -2, -1, 1, 2, 3}) {
      m_cornersFPIX.push_back(edm::FileInPath(Form("DQM/SiStripMonitorClient/data/Geometry/vertices_forward_%i", j)));
    }
  }

~Phase1PixelMaps()

Phase1PixelMaps::~Phase1PixelMaps ( )

default

Member Function Documentation

adjustCanvasMargins()

void Phase1PixelMaps::adjustCanvasMargins ( TVirtualPad *  pad, float  top, float  bottom, float  left, float  right  )

private

Definition at line 529 of file Phase1PixelMaps.cc.

Referenced by drawBarrelMaps(), drawForwardMaps(), and drawSummaryMaps().

                                                                                                            {
  if (top > 0)
    pad->SetTopMargin(top);
  if (bottom > 0)
    pad->SetBottomMargin(bottom);
  if (left > 0)
    pad->SetLeftMargin(left);
  if (right > 0)
    pad->SetRightMargin(right);
}

beautifyAllHistograms()

void Phase1PixelMaps::beautifyAllHistograms

(

)

Definition at line 289 of file Phase1PixelMaps.cc.

References m_isBooked, makeNicePlotStyle(), plotFactory::plot, pxbTh2PolyBarrel, and pxfTh2PolyForward.

Referenced by templateHelper::SiPixelIDs< PayloadType, myType >::fill(), templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill(), and main().

                                            {
  if (!m_isBooked.first && !m_isBooked.second) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to beautify a histogram not booked";
    return;
  }

  // only if the barrel is booked
  if (m_isBooked.first) {
    for (const auto& vec : pxbTh2PolyBarrel) {
      for (const auto& plot : vec.second) {
        this->makeNicePlotStyle(plot.get());
        plot->GetXaxis()->SetTitleOffset(0.9);
        plot->GetYaxis()->SetTitleOffset(0.9);
        plot->GetZaxis()->SetTitleOffset(1.2);
        plot->GetZaxis()->SetTitleSize(0.05);
      }
    }
  }

  // only if the forwards are booked
  if (m_isBooked.second) {
    for (const auto& vec : pxfTh2PolyForward) {
      for (const auto& plot : vec.second) {
        this->makeNicePlotStyle(plot.get());
        plot->GetXaxis()->SetTitleOffset(0.9);
        plot->GetYaxis()->SetTitleOffset(0.9);
        plot->GetZaxis()->SetTitleOffset(1.2);
        plot->GetZaxis()->SetTitleSize(0.05);
      }
    }
  }
}

book()

void Phase1PixelMaps::book	(	const std::string &	currentHistoName,
		const char *	what,
		const char *	zaxis
	)

Definition at line 216 of file Phase1PixelMaps.cc.

References bookBarrelHistograms(), bookForwardHistograms(), m_isBooked, and LaserClient_cfi::zaxis.

Referenced by main().

                                                                                                 {
  bookBarrelHistograms(currentHistoName, what, zaxis);
  bookForwardHistograms(currentHistoName, what, zaxis);
  m_isBooked = std::make_pair(true, true);
}

bookBarrelBins()

void Phase1PixelMaps::bookBarrelBins ( const std::string &  currentHistoName )

private

Definition at line 124 of file Phase1PixelMaps.cc.

References bins, binsSummary, ALCARECOPPSCalTrackBasedSel_cff::detid, mps_splice::entry, Exception, f, ALCARECOEcalPhiSym_cff::float, EcalPhiSymFlatTableProducers_cfi::id, isotrackApplyRegressor::k, PVValHelper::ladder, nano_mu_digi_cff::layer, m_cornersBPIX, m_trackerTopo, PixelSubdetector::PixelBarrel, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), pxbTh2PolyBarrel, pxbTh2PolyBarrelSummary, and retrieveCorners().

Referenced by bookBarrelHistograms().

                                                                      {
  auto theIndexedCorners = this->retrieveCorners(m_cornersBPIX, 4);

  for (const auto& entry : theIndexedCorners) {
    auto id = entry.first;
    auto detid = DetId(id);
    if (detid.subdetId() != PixelSubdetector::PixelBarrel)
      continue;

    int layer = m_trackerTopo.pxbLayer(detid);
    int ladder = m_trackerTopo.pxbLadder(detid);

    auto theVectX = entry.second.first;
    auto theVectY = entry.second.second;

    float vertX[] = {theVectX[0], theVectX[1], theVectX[2], theVectX[3], theVectX[4]};
    float vertY[] = {(ladder - 1.0f), (ladder - 1.0f), (float)ladder, (float)ladder, (ladder - 1.0f)};

    bins[id] = std::make_shared<TGraph>(5, vertX, vertY);
    bins[id]->SetName(TString::Format("%u", id));

    // Summary plot
    for (unsigned k = 0; k < 5; ++k) {
      vertX[k] += ((layer == 2 || layer == 3) ? 0.0f : -60.0f);
      vertY[k] += ((layer > 2) ? 30.0f : 0.0f);
    }

    binsSummary[id] = std::make_shared<TGraph>(5, vertX, vertY);
    binsSummary[id]->SetName(TString::Format("%u", id));

    if (pxbTh2PolyBarrel.find(currentHistoName) != pxbTh2PolyBarrel.end()) {
      pxbTh2PolyBarrel[currentHistoName][layer - 1]->AddBin(bins[id]->Clone());
    } else {
      throw cms::Exception("LogicError") << currentHistoName << " is not found in the Barrel map! Aborting.";
    }

    if (pxbTh2PolyBarrelSummary.find(currentHistoName) != pxbTh2PolyBarrelSummary.end()) {
      pxbTh2PolyBarrelSummary[currentHistoName]->AddBin(binsSummary[id]->Clone());
    } else {
      throw cms::Exception("LocalError") << currentHistoName << " is not found in the Barrel Summary map! Aborting.";
    }
  }
}

bookBarrelHistograms()

void Phase1PixelMaps::bookBarrelHistograms	(	const std::string &	currentHistoName,
		const char *	what,
		const char *	zaxis
	)

Definition at line 31 of file Phase1PixelMaps.cc.

References bookBarrelBins(), spr::find(), mps_fire::i, m_isBooked, m_knownNames, m_option, pxbTh2PolyBarrel, pxbTh2PolyBarrelSummary, AlCaHLTBitMon_QueryRunRegistry::string, to_string(), and LaserClient_cfi::zaxis.

Referenced by book(), templateHelper::SiPixelIDs< PayloadType, myType >::fill(), and templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                                                                                 {
  std::string histName;
  std::shared_ptr<TH2Poly> th2p;

  // check if the passed histogram name already exists, if not store it
  if (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) == m_knownNames.end()) {
    m_knownNames.emplace_back(currentHistoName);
  }

  for (unsigned i = 0; i < 4; ++i) {
    histName = "barrel_layer_";

    th2p = std::make_shared<TH2Poly>(
        (histName + std::to_string(i + 1)).c_str(), Form("PXBMap of %s - Layer %i", what, i + 1), -15.0, 15.0, 0.0, 5.0);

    th2p->SetFloat();

    th2p->GetXaxis()->SetTitle("z [cm]");
    th2p->GetYaxis()->SetTitle("ladder");
    th2p->GetZaxis()->SetTitle(zaxis);
    th2p->GetZaxis()->CenterTitle();
    th2p->SetStats(false);
    th2p->SetOption(m_option);
    pxbTh2PolyBarrel[currentHistoName].push_back(th2p);
  }

  th2p = std::make_shared<TH2Poly>("barrel_summary", Form("Barrel Pixel Map of %s", what), -5.0, 5.0, 0.0, 15.0);
  th2p->SetFloat();

  th2p->GetXaxis()->SetTitle("");
  th2p->GetYaxis()->SetTitle("");
  th2p->GetZaxis()->SetTitle(zaxis);
  th2p->GetZaxis()->CenterTitle();
  th2p->SetStats(false);
  th2p->SetOption(m_option);
  pxbTh2PolyBarrelSummary[currentHistoName] = th2p;

  // book the bins
  bookBarrelBins(currentHistoName);

  // set the isBooked bit to true;
  m_isBooked.first = true;
}

bookForwardBins()

void Phase1PixelMaps::bookForwardBins ( const std::string &  currentHistoName )

private

Definition at line 169 of file Phase1PixelMaps.cc.

References bins, binsSummary, ALCARECOPPSCalTrackBasedSel_cff::detid, mps_splice::entry, Exception, f, ALCARECOEcalPhiSym_cff::float, EcalPhiSymFlatTableProducers_cfi::id, isotrackApplyRegressor::k, m_cornersFPIX, m_trackerTopo, PixelSubdetector::PixelEndcap, TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), pxfTh2PolyForward, pxfTh2PolyForwardSummary, retrieveCorners(), and ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side().

Referenced by bookForwardHistograms().

                                                                       {
  auto theIndexedCorners = this->retrieveCorners(m_cornersFPIX, 3);

  for (const auto& entry : theIndexedCorners) {
    auto id = entry.first;
    auto detid = DetId(id);
    if (detid.subdetId() != PixelSubdetector::PixelEndcap)
      continue;

    int disk = m_trackerTopo.pxfDisk(detid);
    int side = m_trackerTopo.pxfSide(detid);

    unsigned mapIdx = disk + (side - 1) * 3 - 1;

    auto theVectX = entry.second.first;
    auto theVectY = entry.second.second;

    float vertX[] = {theVectX[0], theVectX[1], theVectX[2], theVectX[3]};
    float vertY[] = {theVectY[0], theVectY[1], theVectY[2], theVectY[3]};

    bins[id] = std::make_shared<TGraph>(4, vertX, vertY);
    bins[id]->SetName(TString::Format("%u", id));

    // Summary plot
    for (unsigned k = 0; k < 4; ++k) {
      vertX[k] += (float(side) - 1.5f) * 40.0f;
      vertY[k] += (disk - 1) * 35.0f;
    }

    binsSummary[id] = std::make_shared<TGraph>(4, vertX, vertY);
    binsSummary[id]->SetName(TString::Format("%u", id));

    if (pxfTh2PolyForward.find(currentHistoName) != pxfTh2PolyForward.end()) {
      pxfTh2PolyForward[currentHistoName][mapIdx]->AddBin(bins[id]->Clone());
    } else {
      throw cms::Exception("LogicError") << currentHistoName << " is not found in the Forward map! Aborting.";
    }

    if (pxfTh2PolyForwardSummary.find(currentHistoName) != pxfTh2PolyForwardSummary.end()) {
      pxfTh2PolyForwardSummary[currentHistoName]->AddBin(binsSummary[id]->Clone());
    } else {
      throw cms::Exception("LogicError") << currentHistoName << " is not found in the Forward Summary map! Aborting.";
    }
  }
}

bookForwardHistograms()

void Phase1PixelMaps::bookForwardHistograms	(	const std::string &	currentHistoName,
		const char *	what,
		const char *	zaxis
	)

Definition at line 76 of file Phase1PixelMaps.cc.

References bookForwardBins(), spr::find(), m_isBooked, m_knownNames, m_option, pxfTh2PolyForward, pxfTh2PolyForwardSummary, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), AlCaHLTBitMon_QueryRunRegistry::string, to_string(), and LaserClient_cfi::zaxis.

Referenced by book(), templateHelper::SiPixelIDs< PayloadType, myType >::fill(), and templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                                                                                  {
  std::string histName;
  std::shared_ptr<TH2Poly> th2p;

  // check if the passed histogram name already exists, if not store it
  if (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) == m_knownNames.end()) {
    m_knownNames.emplace_back(currentHistoName);
  }

  for (unsigned side = 1; side <= 2; ++side) {
    for (unsigned disk = 1; disk <= 3; ++disk) {
      histName = "forward_disk_";

      th2p = std::make_shared<TH2Poly>((histName + std::to_string((side == 1 ? -(int(disk)) : (int)disk))).c_str(),
                                       Form("PXFMap of %s - Side %i Disk %i", what, side, disk),
                                       -15.0,
                                       15.0,
                                       -15.0,
                                       15.0);
      th2p->SetFloat();
      th2p->GetXaxis()->SetTitle("x [cm]");
      th2p->GetYaxis()->SetTitle("y [cm]");
      th2p->GetZaxis()->SetTitle(zaxis);
      th2p->GetZaxis()->CenterTitle();
      th2p->SetStats(false);
      th2p->SetOption(m_option);
      pxfTh2PolyForward[currentHistoName].push_back(th2p);
    }
  }

  th2p = std::make_shared<TH2Poly>("forward_summary", Form("Forward Pixel Map of %s", what), -40.0, 40.0, -20.0, 90.0);
  th2p->SetFloat();

  th2p->GetXaxis()->SetTitle("");
  th2p->GetYaxis()->SetTitle("");
  th2p->GetZaxis()->SetTitle(zaxis);
  th2p->GetZaxis()->CenterTitle();
  th2p->SetStats(false);
  th2p->SetOption(m_option);
  pxfTh2PolyForwardSummary[currentHistoName] = th2p;

  // book the bins
  bookForwardBins(currentHistoName);

  m_isBooked.second = true;
}

drawBarrelMaps()

void Phase1PixelMaps::drawBarrelMaps	(	const std::string &	currentHistoName,
		TCanvas &	canvas,
		const char *	drawOption = nullptr
	)

Definition at line 337 of file Phase1PixelMaps.cc.

References adjustCanvasMargins(), histoStyle::drawOption, spr::find(), newFWLiteAna::found, mps_fire::i, m_autorescale, m_isBooked, m_knownNames, m_option, compareTotals::pad1, compareTotals::pad2, pxbTh2PolyBarrel, and rescaleAllBarrel().

Referenced by templateHelper::SiPixelIDs< PayloadType, myType >::fill(), and templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                                                                               {
  auto found = (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) != m_knownNames.end());

  if (!m_isBooked.first || !found) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to draw a histogram not booked";
    return;
  }

  TPad* pad1 = new TPad("pad1", "pad1", 0.0, 0.025, 1.0, 1.0);
  TPad* pad2 = new TPad("pad2", "pad2", 0.0, 0.00, 1.0, 0.025);
  pad1->Divide(2, 2);
  pad1->Draw();
  pad2->Draw();
  for (int i = 1; i <= 4; i++) {
    pad1->cd(i);
    if (strcmp(m_option, "text") == 0) {
      pad1->cd(i)->SetRightMargin(0.02);
      pxbTh2PolyBarrel[currentHistoName].at(i - 1)->SetMarkerColor(kRed);
    } else {
      if (m_autorescale)
        rescaleAllBarrel(currentHistoName);
      adjustCanvasMargins(pad1->cd(i), 0.07, 0.12, 0.10, 0.18);
    }
    if (drawOption) {
      pxbTh2PolyBarrel[currentHistoName].at(i - 1)->Draw("L");
      pxbTh2PolyBarrel[currentHistoName].at(i - 1)->Draw(fmt::sprintf("%s%ssame", m_option, drawOption).c_str());
    } else {
      pxbTh2PolyBarrel[currentHistoName].at(i - 1)->Draw("L");
      pxbTh2PolyBarrel[currentHistoName].at(i - 1)->Draw(fmt::sprintf("%ssame", m_option).c_str());
    }
  }
}

drawForwardMaps()

void Phase1PixelMaps::drawForwardMaps	(	const std::string &	currentHistoName,
		TCanvas &	canvas,
		const char *	drawOption = nullptr
	)

Definition at line 371 of file Phase1PixelMaps.cc.

References adjustCanvasMargins(), histoStyle::drawOption, spr::find(), newFWLiteAna::found, mps_fire::i, m_autorescale, m_isBooked, m_knownNames, m_option, compareTotals::pad1, compareTotals::pad2, pxfTh2PolyForward, and rescaleAllForward().

Referenced by templateHelper::SiPixelIDs< PayloadType, myType >::fill(), and templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                                                                                {
  auto found = (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) != m_knownNames.end());

  if (!m_isBooked.second || !found) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to draw a histogram not booked";
    return;
  }

  TPad* pad1 = new TPad("pad1", "pad1", 0.0, 0.025, 1.0, 1.0);
  TPad* pad2 = new TPad("pad2", "pad2", 0.0, 0.00, 1.0, 0.025);
  pad1->Divide(3, 2);
  pad1->Draw();
  pad2->Draw();

  for (int i = 1; i <= 6; i++) {
    pad1->cd(i);
    if (strcmp(m_option, "text") == 0) {
      pad1->cd(i)->SetRightMargin(0.02);
      pxfTh2PolyForward[currentHistoName].at(i - 1)->SetMarkerColor(kRed);
    } else {
      if (m_autorescale)
        rescaleAllForward(currentHistoName);
      adjustCanvasMargins(pad1->cd(i), 0.07, 0.12, 0.10, 0.18);
    }
    if (drawOption) {
      pxfTh2PolyForward[currentHistoName].at(i - 1)->Draw("L");
      pxfTh2PolyForward[currentHistoName].at(i - 1)->Draw(fmt::sprintf("%s%ssame", m_option, drawOption).c_str());
    } else {
      pxfTh2PolyForward[currentHistoName].at(i - 1)->Draw("L");
      pxfTh2PolyForward[currentHistoName].at(i - 1)->Draw(fmt::sprintf("%ssame", m_option).c_str());
    }
  }
}

drawSummaryMaps()

void Phase1PixelMaps::drawSummaryMaps	(	const std::string &	currentHistoName,
		TCanvas &	canvas,
		const char *	drawOption = nullptr
	)

Definition at line 406 of file Phase1PixelMaps.cc.

References adjustCanvasMargins(), histoStyle::drawOption, spr::find(), newFWLiteAna::found, m_isBooked, m_knownNames, m_option, compareTotals::pad1, compareTotals::pad2, pxbTh2PolyBarrelSummary, pxfTh2PolyForwardSummary, AlCaHLTBitMon_QueryRunRegistry::string, and groupFilesInBlocks::temp.

Referenced by templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill(), and main().

                                                                                                                {
  auto found = (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) != m_knownNames.end());

  if (!m_isBooked.second || !m_isBooked.first || !found) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to draw a histogram not booked";
    return;
  }

  TPad* pad1 = new TPad("pad1", "pad1", 0.0, 0.025, 1.0, 1.0);
  TPad* pad2 = new TPad("pad2", "pad2", 0.0, 0.00, 1.0, 0.025);
  pad1->Divide(2, 1);
  pad1->Draw();
  pad2->Draw();

  pad1->cd(1);
  std::string temp(m_option);  // create a std string
  auto isText = (temp.find("text") != std::string::npos);
  adjustCanvasMargins(pad1->cd(1), 0.07, 0.02, 0.01, isText ? 0.05 : 0.15);
  if (isText) {
    pxbTh2PolyBarrelSummary[currentHistoName]->SetMarkerColor(kRed);
    pxbTh2PolyBarrelSummary[currentHistoName]->SetMarkerSize(0.5);
  }
  pxbTh2PolyBarrelSummary[currentHistoName]->GetZaxis()->SetTitleOffset(1.4);
  pxbTh2PolyBarrelSummary[currentHistoName]->Draw("AL");
  pxbTh2PolyBarrelSummary[currentHistoName]->Draw(fmt::sprintf("%s%ssame", m_option, drawOption).c_str());

  pad1->cd(2);
  adjustCanvasMargins(pad1->cd(2), 0.07, 0.02, 0.01, isText ? 0.05 : 0.15);
  if (isText) {
    pxfTh2PolyForwardSummary[currentHistoName]->SetMarkerColor(kRed);
    pxfTh2PolyForwardSummary[currentHistoName]->SetMarkerSize(0.5);
  }
  pxfTh2PolyForwardSummary[currentHistoName]->GetZaxis()->SetTitleOffset(1.4);
  pxfTh2PolyForwardSummary[currentHistoName]->Draw("AL");
  pxfTh2PolyForwardSummary[currentHistoName]->Draw(fmt::sprintf("%s%ssame", m_option, drawOption).c_str());
}

fill()

void Phase1PixelMaps::fill	(	const std::string &	currentHistoName,
		unsigned int	id,
		double	value
	)

Definition at line 223 of file Phase1PixelMaps.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, nano_mu_digi_cff::layer, LogDebug, m_isBooked, m_trackerTopo, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TrackerTopology::pxbLayer(), pxbTh2PolyBarrel, pxbTh2PolyBarrelSummary, TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), pxfTh2PolyForward, pxfTh2PolyForwardSummary, and ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side().

Referenced by main().

                                                                                           {
  auto detid = DetId(id);
  if (detid.subdetId() == PixelSubdetector::PixelBarrel) {
    int layer = m_trackerTopo.pxbLayer(id);
    if (!m_isBooked.first) {
      edm::LogError("Phase1PixelMaps") << __func__ << ": trying to fill a histogram not booked";
      return;
    }

    LogDebug("Phase1PixelMaps") << __func__ << " filling barrel with value: " << value << std::endl;

    pxbTh2PolyBarrel[currentHistoName][layer - 1]->Fill(TString::Format("%u", id), value);
    pxbTh2PolyBarrelSummary[currentHistoName]->Fill(TString::Format("%u", id), value);
  } else if (detid.subdetId() == PixelSubdetector::PixelEndcap) {
    int disk = m_trackerTopo.pxfDisk(id);
    int side = m_trackerTopo.pxfSide(id);
    unsigned mapIdx = disk + (side - 1) * 3 - 1;
    if (!m_isBooked.second) {
      edm::LogError("Phase1PixelMaps") << __func__ << ": trying to fill a histogram not booked";
      return;
    }

    LogDebug("Phase1PixelMaps") << __func__ << " filling endcaps with value: " << value << std::endl;

    pxfTh2PolyForward[currentHistoName][mapIdx]->Fill(TString::Format("%u", id), value);
    pxfTh2PolyForwardSummary[currentHistoName]->Fill(TString::Format("%u", id), value);
  }
}

fillBarrelBin()

void Phase1PixelMaps::fillBarrelBin	(	const std::string &	currentHistoName,
		unsigned int	id,
		double	value
	)

Definition at line 253 of file Phase1PixelMaps.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, nano_mu_digi_cff::layer, m_isBooked, m_trackerTopo, PixelSubdetector::PixelBarrel, TrackerTopology::pxbLayer(), pxbTh2PolyBarrel, and pxbTh2PolyBarrelSummary.

Referenced by templateHelper::SiPixelIDs< PayloadType, myType >::fill(), and templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                                                                    {
  auto detid = DetId(id);
  if (detid.subdetId() != PixelSubdetector::PixelBarrel) {
    edm::LogError("Phase1PixelMaps") << "fillBarrelBin() The following detid " << id << " is not Pixel Barrel!"
                                     << std::endl;
    return;
  }
  if (!m_isBooked.first) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to fill a histogram not booked";
    return;
  }
  int layer = m_trackerTopo.pxbLayer(id);
  pxbTh2PolyBarrel[currentHistoName][layer - 1]->Fill(TString::Format("%u", id), value);
  pxbTh2PolyBarrelSummary[currentHistoName]->Fill(TString::Format("%u", id), value);
}

fillForwardBin()

void Phase1PixelMaps::fillForwardBin	(	const std::string &	currentHistoName,
		unsigned int	id,
		double	value
	)

Definition at line 270 of file Phase1PixelMaps.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, m_isBooked, m_trackerTopo, PixelSubdetector::PixelEndcap, TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), pxfTh2PolyForward, pxfTh2PolyForwardSummary, and ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side().

Referenced by templateHelper::SiPixelIDs< PayloadType, myType >::fill(), and templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                                                                     {
  auto detid = DetId(id);
  if (detid.subdetId() != PixelSubdetector::PixelEndcap) {
    edm::LogError("Phase1PixelMaps") << "fillForwardBin() The following detid " << id << " is not Pixel Forward!"
                                     << std::endl;
    return;
  }
  if (!m_isBooked.second) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to fill a histogram not booked";
    return;
  }
  int disk = m_trackerTopo.pxfDisk(id);
  int side = m_trackerTopo.pxfSide(id);
  unsigned mapIdx = disk + (side - 1) * 3 - 1;
  pxfTh2PolyForward[currentHistoName][mapIdx]->Fill(TString::Format("%u", id), value);
  pxfTh2PolyForwardSummary[currentHistoName]->Fill(TString::Format("%u", id), value);
}

makeNicePlotStyle()

void Phase1PixelMaps::makeNicePlotStyle ( TH1 *  hist )

private

Definition at line 504 of file Phase1PixelMaps.cc.

References HLT_2024v14_cff::Class, and compareTotals::hist.

Referenced by beautifyAllHistograms().

                                                 {
  hist->SetStats(kFALSE);
  hist->SetLineWidth(2);
  hist->GetXaxis()->CenterTitle(true);
  hist->GetYaxis()->CenterTitle(true);
  hist->GetXaxis()->SetTitleFont(42);
  hist->GetYaxis()->SetTitleFont(42);
  hist->GetXaxis()->SetTitleSize(0.05);
  hist->GetYaxis()->SetTitleSize(0.05);
  hist->GetXaxis()->SetTitleOffset(1.1);
  hist->GetYaxis()->SetTitleOffset(1.3);
  hist->GetXaxis()->SetLabelFont(42);
  hist->GetYaxis()->SetLabelFont(42);
  hist->GetYaxis()->SetLabelSize(.05);
  hist->GetXaxis()->SetLabelSize(.05);

  if (hist->InheritsFrom(TH2::Class())) {
    hist->GetZaxis()->SetLabelFont(42);
    hist->GetZaxis()->SetLabelFont(42);
    hist->GetZaxis()->SetLabelSize(.05);
    hist->GetZaxis()->SetLabelSize(.05);
  }
}

rescaleAllBarrel()

void Phase1PixelMaps::rescaleAllBarrel ( const std::string &  currentHistoName )

private

Definition at line 541 of file Phase1PixelMaps.cc.

References SplitLinear::begin, mps_fire::end, spr::find(), timingPdfMaker::histo, m_knownNames, pxbTh2PolyBarrel, and HcalDetIdTransform::transform().

Referenced by drawBarrelMaps().

                                                                        {
  if (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) == m_knownNames.end()) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to manipulate a histogram not booked";
    return;
  }

  std::vector<float> maxima;
  std::transform(pxbTh2PolyBarrel[currentHistoName].begin(),
                 pxbTh2PolyBarrel[currentHistoName].end(),
                 std::back_inserter(maxima),
                 [](std::shared_ptr<TH2Poly> thp) -> float { return thp->GetMaximum(); });
  std::vector<float> minima;
  std::transform(pxbTh2PolyBarrel[currentHistoName].begin(),
                 pxbTh2PolyBarrel[currentHistoName].end(),
                 std::back_inserter(minima),
                 [](std::shared_ptr<TH2Poly> thp) -> float { return thp->GetMinimum(); });

  auto globalMax = *std::max_element(maxima.begin(), maxima.end());
  auto globalMin = *std::min_element(minima.begin(), minima.end());

  // in case the two coincide do not rescale
  if (globalMax == globalMin)
    return;

  for (auto& histo : pxbTh2PolyBarrel[currentHistoName]) {
    histo->GetZaxis()->SetRangeUser(globalMin, globalMax);
  }
}

rescaleAllForward()

void Phase1PixelMaps::rescaleAllForward ( const std::string &  currentHistoName )

private

Definition at line 571 of file Phase1PixelMaps.cc.

References SplitLinear::begin, mps_fire::end, spr::find(), timingPdfMaker::histo, m_knownNames, pxfTh2PolyForward, and HcalDetIdTransform::transform().

Referenced by drawForwardMaps().

                                                                         {
  if (std::find(m_knownNames.begin(), m_knownNames.end(), currentHistoName) == m_knownNames.end()) {
    edm::LogError("Phase1PixelMaps") << __func__ << ": trying to manipulate a histogram not booked";
    return;
  }

  std::vector<float> maxima;
  std::transform(pxfTh2PolyForward[currentHistoName].begin(),
                 pxfTh2PolyForward[currentHistoName].end(),
                 std::back_inserter(maxima),
                 [](std::shared_ptr<TH2Poly> thp) -> float { return thp->GetMaximum(); });
  std::vector<float> minima;
  std::transform(pxfTh2PolyForward[currentHistoName].begin(),
                 pxfTh2PolyForward[currentHistoName].end(),
                 std::back_inserter(minima),
                 [](std::shared_ptr<TH2Poly> thp) -> float { return thp->GetMinimum(); });

  auto globalMax = *std::max_element(maxima.begin(), maxima.end());
  auto globalMin = *std::min_element(minima.begin(), minima.end());

  // in case the two coincide do not rescale
  if (globalMax == globalMin)
    return;

  for (auto& histo : pxfTh2PolyForward[currentHistoName]) {
    histo->GetZaxis()->SetRangeUser(globalMin, globalMax);
  }
}

resetOption()

void Phase1PixelMaps::resetOption

(

const char *

option

)

Definition at line 21 of file Phase1PixelMaps.cc.

References m_option, and fileinputsource_cfi::option.

Referenced by templateHelper::SiPixelTemplateHeaderInfo< PayloadType, StoreType, TransientType, myType, myParam >::fill().

                                                    {
  if (m_option != nullptr && !m_option[0]) {
    m_option = option;
  } else {
    edm::LogError("Phase1PixelMaps") << "Option has already been set to " << m_option
                                     << ". It's not possible to reset it.";
  }
}

retrieveCorners()

const indexedCorners Phase1PixelMaps::retrieveCorners ( const std::vector< edm::FileInPath > &  cornerFiles, const unsigned int  reads  )

private

Definition at line 444 of file Phase1PixelMaps.cc.

References SplitLinear::begin, mps_fire::end, Exception, geometryDiff::file, mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, mps_splice::line, LOGDEBUG, mergeVDriftHistosByStation::name, MatrixUtil::remove(), AlCaHLTBitMon_QueryRunRegistry::string, and submitPVValidationJobs::t.

Referenced by bookBarrelBins(), and bookForwardBins().

                                                                                {
  indexedCorners theOutMap;

  for (const auto& file : cornerFiles) {
    auto cornerFileName = file.fullPath();
    std::ifstream cornerFile(cornerFileName.c_str());
    if (!cornerFile.good()) {
      throw cms::Exception("FileError") << "Problem opening corner file: " << cornerFileName;
    }
    std::string line;
    while (std::getline(cornerFile, line)) {
      if (!line.empty()) {
        std::istringstream iss(line);
        unsigned int id;
        std::string name;
        std::vector<std::string> corners(reads, "");
        std::vector<float> xP, yP;

        iss >> id >> name;
        for (unsigned int i = 0; i < reads; ++i) {
          iss >> corners.at(i);
        }

        LOGDEBUG("Phase1PixelMaps") << id << " : ";
        for (unsigned int i = 0; i < reads; i++) {
          // remove the leading and trailing " signs in the corners list
          (corners[i]).erase(std::remove(corners[i].begin(), corners[i].end(), '"'), corners[i].end());
          LOGDEBUG("Phase1PixelMaps") << corners.at(i) << " ";
          typedef boost::tokenizer<boost::char_separator<char>> tokenizer;
          boost::char_separator<char> sep{","};
          tokenizer tok{corners.at(i), sep};
          for (const auto& t : tok | boost::adaptors::indexed(0)) {
            if (t.index() == 0) {
              xP.push_back(atof((t.value()).c_str()));
            } else if (t.index() == 1) {
              yP.push_back(atof((t.value()).c_str()));
            } else {
              edm::LogError("LogicError") << "There should not be any token with index " << t.index() << std::endl;
            }
          }
        }
        LOGDEBUG("Phase1PixelMaps") << std::endl;

        xP.push_back(xP.front());
        yP.push_back(yP.front());

        for (unsigned int i = 0; i < xP.size(); i++) {
          LOGDEBUG("Phase1PixelMaps") << "x[" << i << "]=" << xP[i] << " y[" << i << "]" << yP[i] << std::endl;
        }

        theOutMap[id] = std::make_pair(xP, yP);

      }  // if line is empty
    }  // loop on lines
  }  // loop on files
  return theOutMap;
}

setBarrelScale()

void Phase1PixelMaps::setBarrelScale	(	const std::string &	currentHistoName,
		std::pair< float, float >	extrema
	)

Definition at line 323 of file Phase1PixelMaps.cc.

References timingPdfMaker::histo, and pxbTh2PolyBarrel.

                                                                                                     {
  for (auto& histo : pxbTh2PolyBarrel[currentHistoName]) {
    histo->GetZaxis()->SetRangeUser(extrema.first, extrema.second);
  }
}

setForwardScale()

void Phase1PixelMaps::setForwardScale	(	const std::string &	currentHistoName,
		std::pair< float, float >	extrema
	)

Definition at line 330 of file Phase1PixelMaps.cc.

References timingPdfMaker::histo, and pxfTh2PolyForward.

                                                                                                      {
  for (auto& histo : pxfTh2PolyForward[currentHistoName]) {
    histo->GetZaxis()->SetRangeUser(extrema.first, extrema.second);
  }
}

setNoRescale()

void Phase1PixelMaps::setNoRescale ( )

inline

Definition at line 56 of file Phase1PixelMaps.h.

References m_autorescale.

{ m_autorescale = false; }

Member Data Documentation

bins

std::map<uint32_t, std::shared_ptr<TGraph> > Phase1PixelMaps::bins

private

Definition at line 93 of file Phase1PixelMaps.h.

Referenced by bookBarrelBins(), and bookForwardBins().

binsSummary

std::map<uint32_t, std::shared_ptr<TGraph> > Phase1PixelMaps::binsSummary

private

Definition at line 93 of file Phase1PixelMaps.h.

Referenced by bookBarrelBins(), and bookForwardBins().

m_autorescale

bool Phase1PixelMaps::m_autorescale

private

Definition at line 87 of file Phase1PixelMaps.h.

Referenced by drawBarrelMaps(), drawForwardMaps(), and setNoRescale().

m_cornersBPIX

std::vector<edm::FileInPath> Phase1PixelMaps::m_cornersBPIX

private

Definition at line 99 of file Phase1PixelMaps.h.

Referenced by bookBarrelBins().

m_cornersFPIX

std::vector<edm::FileInPath> Phase1PixelMaps::m_cornersFPIX

private

Definition at line 100 of file Phase1PixelMaps.h.

Referenced by bookForwardBins().

m_isBooked

std::pair<bool, bool> Phase1PixelMaps::m_isBooked

private

Definition at line 88 of file Phase1PixelMaps.h.

Referenced by beautifyAllHistograms(), book(), bookBarrelHistograms(), bookForwardHistograms(), drawBarrelMaps(), drawForwardMaps(), drawSummaryMaps(), fill(), fillBarrelBin(), and fillForwardBin().

m_knownNames

std::vector<std::string> Phase1PixelMaps::m_knownNames

private

Definition at line 89 of file Phase1PixelMaps.h.

Referenced by bookBarrelHistograms(), bookForwardHistograms(), drawBarrelMaps(), drawForwardMaps(), drawSummaryMaps(), rescaleAllBarrel(), and rescaleAllForward().

m_option

Option_t* Phase1PixelMaps::m_option

private

Definition at line 86 of file Phase1PixelMaps.h.

Referenced by bookBarrelHistograms(), bookForwardHistograms(), drawBarrelMaps(), drawForwardMaps(), drawSummaryMaps(), and resetOption().

m_trackerTopo

TrackerTopology Phase1PixelMaps::m_trackerTopo

private

Definition at line 91 of file Phase1PixelMaps.h.

Referenced by bookBarrelBins(), bookForwardBins(), fill(), fillBarrelBin(), and fillForwardBin().

pxbTh2PolyBarrel

std::map<std::string, std::vector<std::shared_ptr<TH2Poly> > > Phase1PixelMaps::pxbTh2PolyBarrel

private

Definition at line 94 of file Phase1PixelMaps.h.

Referenced by beautifyAllHistograms(), bookBarrelBins(), bookBarrelHistograms(), drawBarrelMaps(), fill(), fillBarrelBin(), rescaleAllBarrel(), and setBarrelScale().

pxbTh2PolyBarrelSummary

std::map<std::string, std::shared_ptr<TH2Poly> > Phase1PixelMaps::pxbTh2PolyBarrelSummary

private

Definition at line 95 of file Phase1PixelMaps.h.

Referenced by bookBarrelBins(), bookBarrelHistograms(), drawSummaryMaps(), fill(), and fillBarrelBin().

pxfTh2PolyForward

std::map<std::string, std::vector<std::shared_ptr<TH2Poly> > > Phase1PixelMaps::pxfTh2PolyForward

private

Definition at line 96 of file Phase1PixelMaps.h.

Referenced by beautifyAllHistograms(), bookForwardBins(), bookForwardHistograms(), drawForwardMaps(), fill(), fillForwardBin(), rescaleAllForward(), and setForwardScale().

pxfTh2PolyForwardSummary

std::map<std::string, std::shared_ptr<TH2Poly> > Phase1PixelMaps::pxfTh2PolyForwardSummary

private

Definition at line 97 of file Phase1PixelMaps.h.

Referenced by bookForwardBins(), bookForwardHistograms(), drawSummaryMaps(), fill(), and fillForwardBin().