HtrXmlPatternTool Class Reference

#include <HtrXmlPatternTool.h>

Public Member Functions
void	createHists ()
void	Fill (const HcalElectronicsId HEID, HBHEDigiCollection::const_iterator data)
void	Fill (const HcalElectronicsId HEID, HFDigiCollection::const_iterator data)
void	Fill (const HcalElectronicsId HEID, HODigiCollection::const_iterator data)
HtrXmlPatternSet *	GetPatternSet ()
	HtrXmlPatternTool (HtrXmlPatternToolParameters *m_params)
void	prepareDirs ()
void	writeXML ()
	~HtrXmlPatternTool ()

Private Attributes
HtrXmlPatternToolParameters *	m_params
HtrXmlPatternSet *	m_patternSet
HtrXmlPatternWriter	m_xmlWriter

Detailed Description

Definition at line 9 of file HtrXmlPatternTool.h.

Constructor & Destructor Documentation

HtrXmlPatternTool::HtrXmlPatternTool

(

HtrXmlPatternToolParameters *

m_params

)

Definition at line 11 of file HtrXmlPatternTool.cc.

References mps_fire::i, HtrXmlPatternToolParameters::m_file_tag, m_params, m_patternSet, m_xmlWriter, ChannelPattern::NUM_CRATES, ChannelPattern::NUM_SLOTS, and HtrXmlPatternWriter::setTagName().

                                                                        {
  m_params = params;

  //slotsOn array index corresponds to physical number (e.g., slotsOn[1] for slot 1)
  int slotsOn[ChannelPattern::NUM_SLOTS];
  int cratesOn[ChannelPattern::NUM_CRATES];

  //turn off all slots
  for (int i = 0; i < ChannelPattern::NUM_SLOTS; i++)
    slotsOn[i] = 0;
  //turn on slots 2,3,4,5,6,7,8
  for (int i = 2; i < 9; i++)
    slotsOn[i] = 1;
  //turn on slots 13,14,15,16,17,18
  for (int i = 13; i < 19; i++)
    slotsOn[i] = 1;

  //turn on all crates
  for (int i = 0; i < ChannelPattern::NUM_CRATES; i++)
    cratesOn[i] = 1;
  //turn off two unused crates
  cratesOn[8] = 0;
  cratesOn[16] = 0;

  m_patternSet = new HtrXmlPatternSet(cratesOn, slotsOn);
  m_xmlWriter.setTagName(m_params->m_file_tag);
}

HtrXmlPatternTool::~HtrXmlPatternTool

(

)

Definition at line 39 of file HtrXmlPatternTool.cc.

References m_patternSet.

{ delete m_patternSet; }

Member Function Documentation

void HtrXmlPatternTool::createHists

(

)

Definition at line 194 of file HtrXmlPatternTool.cc.

References hippyaddtobaddatafiles::cd(), MessageLogger_cfi::cerr, officialStyle::chan, gather_cfg::cout, SimDataFormats::CaloAnalysis::cp, ChannelPattern::getCoded(), HtrXmlPatternSet::getCrate(), CrateData::getHalfHtrData(), HalfHtrData::getPattern(), ChannelPattern::getQuantized(), AnalysisDataFormats_SUSYBSMObjects::hp, mps_fire::i, HtrXmlPatternToolParameters::m_file_tag, HtrXmlPatternToolParameters::m_output_directory, m_params, m_patternSet, dataset::name, ChannelPattern::NUM_CRATES, ChannelPattern::NUM_SLOTS, ChannelPattern::SAMPLES, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by HtrXmlPattern::endJob().

                                    {
  std::cout << "Writing root file..." << std::endl;
  std::string name = m_params->m_output_directory + "/" + (m_params->m_file_tag) + ".root";
  TFile of(name.c_str(), "RECREATE");
  if (of.IsZombie()) {
    std::cerr << "root output file " << name << " is bad." << std::endl;
    return;
  }

  of.mkdir("adc");

  for (int crate = 0; crate < ChannelPattern::NUM_CRATES; crate++) {
    CrateData* cd = m_patternSet->getCrate(crate);
    if (cd == nullptr)
      continue;
    for (int slot = 0; slot < ChannelPattern::NUM_SLOTS; slot++) {
      for (int tb = 0; tb <= 1; tb++) {
        HalfHtrData* hd = cd->getHalfHtrData(slot, tb);
        if (hd == nullptr)
          continue;
        for (int chan = 1; chan <= 24; chan++) {
          ChannelPattern* cp = hd->getPattern(chan);
          char hname[128];
          sprintf(hname, "Exact fC Cr%d,%d%s-%d", crate, slot, ((tb == 1) ? ("t") : ("b")), chan);
          TH1* hp = new TH1F(hname, hname, ChannelPattern::SAMPLES, -0.5, ChannelPattern::SAMPLES - 0.5);
          hp->SetDirectory(nullptr);
          sprintf(hname, "Quantized fC Cr%d,%d%s-%d", crate, slot, ((tb == 1) ? ("t") : ("b")), chan);
          TH1* hq = new TH1F(hname, hname, ChannelPattern::SAMPLES, -0.5, ChannelPattern::SAMPLES - 0.5);
          hp->SetDirectory(nullptr);
          sprintf(hname, "Encoded fC Cr%d,%d%s-%d", crate, slot, ((tb == 1) ? ("t") : ("b")), chan);
          TH1* ha = new TH1F(hname, hname, ChannelPattern::SAMPLES, -0.5, ChannelPattern::SAMPLES - 0.5);
          ha->SetDirectory(nullptr);
          for (int i = 0; i < ChannelPattern::SAMPLES; i++) {
            hp->Fill(i * 1.0, (*cp)[i]);
            hq->Fill(i * 1.0, cp->getQuantized(i));
            ha->Fill(i * 1.0, cp->getCoded(i));
          }
          //of.cd("perfect");   hp->Write();
          //of.cd("quantized"); hq->Write();
          of.cd("adc");
          ha->Write();
          delete hp;
          delete hq;
          delete ha;
        }
      }
    }
  }
  of.Close();
}

void HtrXmlPatternTool::Fill	(	const HcalElectronicsId	HEID,
		HBHEDigiCollection::const_iterator	data
	)

Definition at line 41 of file HtrXmlPatternTool.cc.

References hippyaddtobaddatafiles::cd(), MessageLogger_cfi::cerr, SimDataFormats::CaloAnalysis::cp, data, HcalElectronicsId::dccid(), ChannelPattern::Fill(), HtrXmlPatternSet::getCrate(), CrateData::getHalfHtrData(), HalfHtrData::getPattern(), HcalElectronicsId::htrChanId(), HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), m_params, m_patternSet, HtrXmlPatternToolParameters::m_show_errors, HcalElectronicsId::readoutVMECrateId(), HalfHtrData::setDCC(), HalfHtrData::setSpigot(), and HcalElectronicsId::spigot().

Referenced by HtrXmlPattern::analyze().

                                                                                                {
  CrateData* cd = m_patternSet->getCrate(HEID.readoutVMECrateId());
  HalfHtrData* hd = nullptr;
  ChannelPattern* cp = nullptr;

  if (cd) {
    hd = cd->getHalfHtrData(HEID.htrSlot(), HEID.htrTopBottom());
    if (hd) {
      hd->setSpigot(HEID.spigot());
      hd->setDCC(HEID.dccid());

      cp = hd->getPattern(HEID.htrChanId());
      if (cp)
        cp->Fill(m_params, data);
      else if (m_params->m_show_errors)
        std::cerr << "Bad (crate,slot,channel): (" << HEID.readoutVMECrateId() << "," << HEID.htrSlot() << ","
                  << HEID.htrChanId() << ")" << std::endl;
    } else if (m_params->m_show_errors)
      std::cerr << "Bad (crate,slot): (" << HEID.readoutVMECrateId() << "," << HEID.htrSlot() << ")" << std::endl;
  } else if (m_params->m_show_errors)
    std::cerr << "Bad (crate): (" << HEID.readoutVMECrateId() << ")" << std::endl;
}

void HtrXmlPatternTool::Fill	(	const HcalElectronicsId	HEID,
		HFDigiCollection::const_iterator	data
	)

void HtrXmlPatternTool::Fill	(	const HcalElectronicsId	HEID,
		HODigiCollection::const_iterator	data
	)

HtrXmlPatternSet* HtrXmlPatternTool::GetPatternSet ( )

inline

Definition at line 19 of file HtrXmlPatternTool.h.

References m_patternSet.

Referenced by HtrXmlPattern::do_hand_fill().

{ return m_patternSet; }

void HtrXmlPatternTool::prepareDirs

(

)

Definition at line 110 of file HtrXmlPatternTool.cc.

References HtrXmlPatternToolParameters::m_output_directory, and m_params.

Referenced by HtrXmlPattern::endJob().

{ boost::filesystem::create_directory(m_params->m_output_directory); }

void HtrXmlPatternTool::writeXML

(

)

Definition at line 112 of file HtrXmlPatternTool.cc.

References hippyaddtobaddatafiles::cd(), MessageLogger_cfi::cerr, gather_cfg::cout, HtrXmlPatternSet::getCrate(), CrateData::getHalfHtrData(), HtrXmlPatternToolParameters::m_file_tag, HtrXmlPatternToolParameters::m_output_directory, m_params, m_patternSet, HtrXmlPatternToolParameters::m_XML_file_mode, m_xmlWriter, dataset::name, ChannelPattern::NUM_CRATES, ChannelPattern::NUM_SLOTS, MillePedeFileConverter_cfg::out, AlCaHLTBitMon_QueryRunRegistry::string, estimatePileup_makeJSON::trunc, and HtrXmlPatternWriter::writePattern().

Referenced by HtrXmlPattern::endJob().

                                 {
  std::cout << "Writing XML..." << std::endl;
  std::ofstream* of = nullptr;

  if (m_params->m_XML_file_mode == 1) {
    std::string name = m_params->m_output_directory + (m_params->m_file_tag) + "_all.xml";
    of = new std::ofstream(name.c_str(), std::ios_base::out | std::ios_base::trunc);
    if (!of->good()) {
      std::cerr << "XML output file " << name << " is bad." << std::endl;
      return;
    }
    (*of) << "<?xml version='1.0' encoding='UTF-8'?>" << std::endl;
    (*of) << "<CFGBrickSet name='" << m_params->m_file_tag << "'>" << std::endl;
  }

  for (int crate = 0; crate < ChannelPattern::NUM_CRATES; crate++) {
    CrateData* cd = m_patternSet->getCrate(crate);
    if (cd == nullptr)
      continue;

    if (m_params->m_XML_file_mode == 2) {
      std::string name = m_params->m_output_directory + (m_params->m_file_tag);
      char cr_name[256];
      snprintf(cr_name, 256, "_crate_%d.xml", crate);
      name += cr_name;
      of = new std::ofstream(name.c_str(), std::ios_base::out | std::ios_base::trunc);
      if (!of->good()) {
        std::cerr << "XML output file " << name << " is bad." << std::endl;
        delete of;
        return;
      }
      (*of) << "<?xml version='1.0' encoding='UTF-8'?>" << std::endl;
      (*of) << "<CFGBrickSet name='" << m_params->m_file_tag << "'>" << std::endl;
    }

    for (int slot = 0; slot < ChannelPattern::NUM_SLOTS; slot++) {
      for (int tb = 0; tb <= 1; tb++) {
        HalfHtrData* hd = cd->getHalfHtrData(slot, tb);
        if (hd == nullptr)
          continue;
        for (int fiber = 1; fiber <= 8; fiber++) {
          if (m_params->m_XML_file_mode == 3) {
            std::string name = m_params->m_output_directory + (m_params->m_file_tag);
            char cr_name[256];
            snprintf(cr_name, 256, "_crate_%d_slot_%d_tb_%d_fiber_%d.xml", crate, slot, tb, fiber);
            name += cr_name;
            of = new std::ofstream(name.c_str(), std::ios_base::out | std::ios_base::trunc);
            if (!of->good()) {
              std::cerr << "XML output file " << name << " is bad." << std::endl;
              delete of;
              return;
            }
            (*of) << "<?xml version='1.0' encoding='UTF-8'?>" << std::endl;
          }
          m_xmlWriter.writePattern(hd, fiber, *of, 1);
          if (m_params->m_XML_file_mode == 3) {
            of->close();
            delete of;
            of = nullptr;
          }

        }  //end fiber loop
      }    // end tb loop
    }      //end slot loop

    if (m_params->m_XML_file_mode == 2) {
      (*of) << "</CFGBrickSet>" << std::endl;
      of->close();
      delete of;
      of = nullptr;
    }

  }  //end crate loop

  if (m_params->m_XML_file_mode == 1) {
    (*of) << "</CFGBrickSet>" << std::endl;
    of->close();
    delete of;
    of = nullptr;
  }
}

Member Data Documentation

HtrXmlPatternToolParameters* HtrXmlPatternTool::m_params

private

Definition at line 23 of file HtrXmlPatternTool.h.

Referenced by createHists(), Fill(), HtrXmlPatternTool(), prepareDirs(), and writeXML().

HtrXmlPatternSet* HtrXmlPatternTool::m_patternSet

private

Definition at line 22 of file HtrXmlPatternTool.h.

Referenced by createHists(), Fill(), GetPatternSet(), HtrXmlPatternTool(), writeXML(), and ~HtrXmlPatternTool().

HtrXmlPatternWriter HtrXmlPatternTool::m_xmlWriter

private

Definition at line 24 of file HtrXmlPatternTool.h.

Referenced by HtrXmlPatternTool(), and writeXML().