PatternOptimizerBase Class Reference

#include <PatternOptimizerBase.h>

Inheritance diagram for PatternOptimizerBase:

Public Member Functions
void	endJob () override
void	observeEventEnd (const edm::Event &iEvent, std::unique_ptr< l1t::RegionalMuonCandBxCollection > &finalCandidates) override
	PatternOptimizerBase (const edm::ParameterSet &edmCfg, const OMTFConfiguration *omtfConfig, GoldenPatternVec< GoldenPatternWithStat > &gps)
	~PatternOptimizerBase () override
Public Member Functions inherited from EmulationObserverBase
	EmulationObserverBase (const edm::ParameterSet &edmCfg, const OMTFConfiguration *omtfConfig)
const SimTrack *	findSimMuon (const edm::Event &event, const SimTrack *previous=nullptr)
void	observeEventBegin (const edm::Event &iEvent) override
void	observeProcesorEmulation (unsigned int iProcessor, l1t::tftype mtfType, const std::shared_ptr< OMTFinput > &input, const AlgoMuons &algoCandidates, const AlgoMuons &gbCandidates, const std::vector< l1t::RegionalMuonCand > &candMuons) override
	~EmulationObserverBase () override
Public Member Functions inherited from IOMTFEmulationObserver
virtual void	beginRun (edm::EventSetup const &eventSetup)
	IOMTFEmulationObserver ()
virtual	~IOMTFEmulationObserver ()

Static Public Member Functions
static double	vxIntegMuRate (double pt_GeV, double dpt, double etaFrom, double etaTo)
static double	vxMuRate (double pt_GeV)

Protected Member Functions
virtual double	getEventRateWeight (double pt)
void	printPatterns ()
virtual void	saveHists (TFile &outfile)
void	savePatternsInRoot (std::string rootFileName)

Protected Attributes
GoldenPatternVec < GoldenPatternWithStat > &	goldenPatterns
TH1I *	simMuFoundByOmtfPt
TH1I *	simMuPt
TH1F *	simMuPtSpectrum
bool	writeLayerStat = false
Protected Attributes inherited from EmulationObserverBase
unsigned int	candProcIndx = 0
edm::ParameterSet	edmCfg
AlgoMuons::value_type	omtfCand
const OMTFConfiguration *	omtfConfig
l1t::RegionalMuonCand	regionalMuonCand
const SimTrack *	simMuon = nullptr

Detailed Description

Definition at line 25 of file PatternOptimizerBase.h.

Constructor & Destructor Documentation

PatternOptimizerBase::PatternOptimizerBase	(	const edm::ParameterSet &	edmCfg,
		const OMTFConfiguration *	omtfConfig,
		GoldenPatternVec< GoldenPatternWithStat > &	gps
	)

Definition at line 64 of file PatternOptimizerBase.cc.

References edm::ParameterSet::exists(), goldenPatterns, simMuFoundByOmtfPt, simMuPt, and simMuPtSpectrum.

    : EmulationObserverBase(edmCfg, omtfConfig), goldenPatterns(gps) {
  // TODO Auto-generated constructor stub

  simMuPt = new TH1I("simMuPt", "simMuPt", goldenPatterns.size(), -0.5, goldenPatterns.size() - 0.5);
  simMuFoundByOmtfPt =
      new TH1I("simMuFoundByOmtfPt", "simMuFoundByOmtfPt", goldenPatterns.size(), -0.5, goldenPatterns.size() - 0.5);

  simMuPtSpectrum = new TH1F("simMuPtSpectrum", "simMuPtSpectrum", 800, 0, 400);

  if (edmCfg.exists("simTracksTag") == false)
    edm::LogError("l1tOmtfEventPrint") << "simTracksTag not found !!!" << std::endl;
}

PatternOptimizerBase::~PatternOptimizerBase ( )

override

Definition at line 80 of file PatternOptimizerBase.cc.

{}

Member Function Documentation

void PatternOptimizerBase::endJob ( void  )

overridevirtual

Implements IOMTFEmulationObserver.

Definition at line 115 of file PatternOptimizerBase.cc.

References EmulationObserverBase::edmCfg, isotrackTrainRegressor::fName, edm::ParameterSet::getParameter(), goldenPatterns, EmulationObserverBase::omtfConfig, savePatternsInRoot(), AlCaHLTBitMon_QueryRunRegistry::string, and XMLConfigWriter::writeGPs().

Referenced by PatternGenerator::endJob().

                                  {
  std::string fName = edmCfg.getParameter<std::string>("optimisedPatsXmlFile");
  edm::LogImportant("PatternOptimizer") << " Writing optimized patterns to " << fName << std::endl;
  XMLConfigWriter xmlWriter(omtfConfig, false, false);
  xmlWriter.writeGPs(goldenPatterns, fName);

  fName.replace(fName.find('.'), fName.length(), ".root");
  savePatternsInRoot(fName);
}

virtual double PatternOptimizerBase::getEventRateWeight ( double  pt )

inlineprotectedvirtual

Definition at line 46 of file PatternOptimizerBase.h.

Referenced by observeEventEnd().

{ return 1; }

void PatternOptimizerBase::observeEventEnd ( const edm::Event &  iEvent, std::unique_ptr< l1t::RegionalMuonCandBxCollection > &  finalCandidates  )

overridevirtual

Reimplemented from IOMTFEmulationObserver.

Definition at line 100 of file PatternOptimizerBase.cc.

References funct::abs(), getEventRateWeight(), OMTFConfiguration::getPatternNum(), goldenPatterns, GoldenPattern::key(), CoreSimTrack::momentum(), EmulationObserverBase::omtfCand, EmulationObserverBase::omtfConfig, simMuFoundByOmtfPt, EmulationObserverBase::simMuon, simMuPtSpectrum, Key::theNumber, and CoreSimTrack::type().

Referenced by PatternGenerator::observeEventEnd().

                                                                                                            {
  if (simMuon == nullptr || omtfCand->getGoldenPatern() == nullptr)  //no sim muon or empty candidate
    return;

  double ptSim = simMuon->momentum().pt();
  int chargeSim = (abs(simMuon->type()) == 13) ? simMuon->type() / -13 : 0;

  unsigned int exptPatNum = omtfConfig->getPatternNum(ptSim, chargeSim);
  GoldenPatternWithStat* exptCandGp = goldenPatterns.at(exptPatNum).get();  // expected pattern
  simMuFoundByOmtfPt->Fill(exptCandGp->key().theNumber);                    //TODO add weight of the muons pt spectrum

  simMuPtSpectrum->Fill(ptSim, getEventRateWeight(ptSim));
}

void PatternOptimizerBase::printPatterns ( )

protected

Definition at line 82 of file PatternOptimizerBase.cc.

References OMTFConfiguration::getPatternPtRange(), goldenPatterns, RPCConst::iptFromPt(), submitPVResolutionJobs::key, EmulationObserverBase::omtfConfig, DiDispStaMuonMonitor_cfi::pt, and OMTFConfiguration::PatternPt::ptFrom.

                                         {
  edm::LogVerbatim("l1tOmtfEventPrint") << __FUNCTION__ << ": " << __LINE__ << " called!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! "
                                        << std::endl;
  for (int patNum = goldenPatterns.size() - 1; patNum >= 0; patNum--) {
    double pt = omtfConfig->getPatternPtRange(patNum).ptFrom;
    if (pt > 0) {
      edm::LogVerbatim("l1tOmtfEventPrint")
          << "cmsRun runThresholdCalc.py " << patNum << " " << (patNum + 1) << " _" << RPCConst::iptFromPt(pt) << "_";
      if (goldenPatterns[patNum]->key().theCharge == -1)
        edm::LogVerbatim("l1tOmtfEventPrint") << "m_";
      else
        edm::LogVerbatim("l1tOmtfEventPrint") << "p_";

      edm::LogVerbatim("l1tOmtfEventPrint") << " > out" << patNum << ".txt" << std::endl;
    }
  }
}

virtual void PatternOptimizerBase::saveHists ( TFile &  outfile )

inlineprotectedvirtual

Reimplemented in PatternGenerator.

Definition at line 48 of file PatternOptimizerBase.h.

Referenced by savePatternsInRoot().

{};

void PatternOptimizerBase::savePatternsInRoot ( std::string  rootFileName )

protected

Definition at line 125 of file PatternOptimizerBase.cc.

References svgfig::canvas(), OMTFConfiguration::getPatternPtRange(), OMTFConfiguration::getRefToLogicNumber(), goldenPatterns, runTauDisplay::gp, gpuVertexFinder::hist, OMTFConfiguration::nPdfBins(), EmulationObserverBase::omtfConfig, timingPdfMaker::outfile, OMTFConfiguration::pdfMaxValue(), OMTFConfiguration::PatternPt::ptFrom, OMTFConfiguration::PatternPt::ptTo, saveHists(), simMuFoundByOmtfPt, simMuPtSpectrum, and writeLayerStat.

Referenced by endJob().

                                                                    {
  gStyle->SetOptStat(111111);
  TFile outfile(rootFileName.c_str(), "RECREATE");
  edm::LogVerbatim("l1tOmtfEventPrint") << __FUNCTION__ << ": " << __LINE__ << " out fileName " << rootFileName
                                        << " outfile->GetName() " << outfile.GetName() << " writeLayerStat "
                                        << writeLayerStat << endl;

  outfile.cd();
  simMuFoundByOmtfPt->Write();

  simMuPtSpectrum->Write();

  outfile.mkdir("patternsPdfs")->cd();
  outfile.mkdir("patternsPdfs/canvases");
  outfile.mkdir("layerStats");
  ostringstream ostrName;
  ostringstream ostrTtle;
  vector<TH1F*> classProbHists;
  for (unsigned int iRefLayer = 0; iRefLayer < goldenPatterns[0]->getPdf()[0].size(); ++iRefLayer) {
    ostrName.str("");
    ostrName << "Neg_RefLayer_" << iRefLayer;
    ostrTtle.str("");
    ostrTtle << "Neg_RefLayer_" << iRefLayer;
    classProbHists.push_back(new TH1F(
        ostrName.str().c_str(), ostrTtle.str().c_str(), goldenPatterns.size(), -0.5, goldenPatterns.size() - 0.5));

    ostrName.str("");
    ostrName << "Pos_RefLayer_" << iRefLayer;
    ostrTtle.str("");
    ostrTtle << "Pos_RefLayer_" << iRefLayer;
    classProbHists.push_back(new TH1F(
        ostrName.str().c_str(), ostrTtle.str().c_str(), goldenPatterns.size(), -0.5, goldenPatterns.size() - 0.5));
  }

  for (auto& gp : goldenPatterns) {
    OMTFConfiguration::PatternPt patternPt = omtfConfig->getPatternPtRange(gp->key().theNumber);
    if (gp->key().thePt == 0)
      continue;
    //edm::LogVerbatim("l1tOmtfEventPrint") <<__FUNCTION__<<": "<<__LINE__<<" "<<gp->key()<<std::endl;
    ostrName.str("");
    ostrName << "PatNum_" << gp->key().theNumber;
    ostrTtle.str("");
    ostrTtle << "PatNum_" << gp->key().theNumber << "_ptCode_" << gp->key().thePt << "_Pt_" << patternPt.ptFrom << "_"
             << patternPt.ptTo << "_GeV";
    TCanvas* canvas = new TCanvas(ostrName.str().c_str(), ostrTtle.str().c_str(), 1200, 1000);
    canvas->Divide(gp->getPdf().size(), gp->getPdf()[0].size(), 0, 0);

    for (unsigned int iLayer = 0; iLayer < gp->getPdf().size(); ++iLayer) {
      for (unsigned int iRefLayer = 0; iRefLayer < gp->getPdf()[iLayer].size(); ++iRefLayer) {
        canvas->cd(1 + iLayer + iRefLayer * gp->getPdf().size());
        ostrName.str("");
        ostrName << "PatNum_" << gp->key().theNumber << "_refLayer_" << iRefLayer << "_Layer_" << iLayer;
        ostrTtle.str("");
        ostrTtle << "PatNum " << gp->key().theNumber << " ptCode " << gp->key().thePt << " refLayer " << iRefLayer
                 << " Layer " << iLayer << " meanDistPhi " << gp->meanDistPhi[iLayer][iRefLayer][0]
                 << " distPhiBitShift "
                 << gp->getDistPhiBitShift(iLayer, iRefLayer);  //"_Pt_"<<patternPt.ptFrom<<"_"<<patternPt.ptTo<<"_GeV
        //edm::LogVerbatim("l1tOmtfEventPrint") <<__FUNCTION__<<": "<<__LINE__<<" creating hist "<<ostrTtle.str()<<std::endl;
        TH1F* hist = new TH1F(
            ostrName.str().c_str(), ostrTtle.str().c_str(), omtfConfig->nPdfBins(), -0.5, omtfConfig->nPdfBins() - 0.5);
        for (unsigned int iPdf = 0; iPdf < gp->getPdf()[iLayer][iRefLayer].size(); iPdf++) {
          hist->Fill(iPdf, gp->pdfAllRef[iLayer][iRefLayer][iPdf]);
        }
        if ((int)iLayer == (omtfConfig->getRefToLogicNumber()[iRefLayer]))
          hist->SetLineColor(kGreen);

        hist->GetYaxis()->SetRangeUser(0, omtfConfig->pdfMaxValue() + 1);
        hist->Draw("hist");

        outfile.cd("patternsPdfs");
        hist->Write();

        if (writeLayerStat) {
          string histName = "histLayerStat_" + ostrName.str();
          unsigned int binCnt = gp->getStatistics()[iLayer][iRefLayer].size();
          TH1I* histLayerStat = new TH1I(histName.c_str(), histName.c_str(), binCnt, -0.5, binCnt - 0.5);
          for (unsigned int iBin = 0; iBin < binCnt; iBin++) {
            histLayerStat->Fill(iBin, gp->getStatistics()[iLayer][iRefLayer][iBin][0]);
          }

          outfile.cd("layerStats");
          histLayerStat->Write();
        }
      }
    }
    outfile.cd("patternsPdfs/canvases");
    canvas->Write();
    delete canvas;

    unsigned int iPdf = omtfConfig->nPdfBins() / 2;
    for (unsigned int iRefLayer = 0; iRefLayer < gp->getPdf()[0].size(); ++iRefLayer) {
      unsigned int refLayerLogicNumber = omtfConfig->getRefToLogicNumber()[iRefLayer];
      if (gp->key().theCharge == -1) {
        classProbHists[2 * iRefLayer]->Fill(gp->key().theNumber, gp->pdfAllRef[refLayerLogicNumber][iRefLayer][iPdf]);
      } else
        classProbHists[2 * iRefLayer + 1]->Fill(gp->key().theNumber,
                                                gp->pdfAllRef[refLayerLogicNumber][iRefLayer][iPdf]);
    }
  }

  outfile.cd();
  for (auto& classProbHist : classProbHists) {
    classProbHist->Write();
  }

  saveHists(outfile);

  outfile.Close();
}

double PatternOptimizerBase::vxIntegMuRate ( double  pt_GeV, double  dpt, double  etaFrom, double  etaTo  )

static

Definition at line 55 of file PatternOptimizerBase.cc.

References RPCpg::rate(), and vxMuRate().

Referenced by PatternGenerator::modifyClassProb().

                                                                                                  {
  //integration using trapeze method - not exact but good enough
  double rate = 0.5 * (vxMuRate(pt_GeV) + vxMuRate(pt_GeV + dpt)) * dpt;

  rate = rate * (etaTo - etaFrom);
  //edm::LogError("l1tOmtfEventPrint")<<ptCode<<" "<<rate;//<<<<<<<<<<<<<<<<<<<<<<<<
  return rate;
}

double PatternOptimizerBase::vxMuRate ( double  pt_GeV )

static

Definition at line 35 of file PatternOptimizerBase.cc.

References a, funct::exp(), RPCpg::mu, and RPCpg::rate().

Referenced by vxIntegMuRate().

                                                   {
  if (pt_GeV == 0)
    return 0.0;
  const double lum = 2.0e34;  //defoult is 1.0e34;
  const double dabseta = 1.0;
  const double dpt = 1.0;
  const double afactor = 1.0e-34 * lum * dabseta * dpt;
  const double a = 2 * 1.3084E6;
  const double mu = -0.725;
  const double sigma = 0.4333;
  const double s2 = 2 * sigma * sigma;

  double ptlog10;
  ptlog10 = log10(pt_GeV);
  double ex = (ptlog10 - mu) * (ptlog10 - mu) / s2;
  double rate = (a * exp(-ex) * afactor);
  //edm::LogError("l1tOmtfEventPrint")<<ptCode<<" "<<rate;//<<<<<<<<<<<<<<<<<<<<<<<<
  return rate;
}

Member Data Documentation

GoldenPatternVec<GoldenPatternWithStat>& PatternOptimizerBase::goldenPatterns

protected

Definition at line 48 of file PatternOptimizerBase.h.

Referenced by PatternGenerator::endJob(), endJob(), PatternGenerator::groupPatterns(), PatternGenerator::initPatternGen(), PatternGenerator::modifyClassProb(), observeEventEnd(), PatternOptimizerBase(), printPatterns(), PatternGenerator::reCalibratePt(), savePatternsInRoot(), PatternGenerator::upadatePdfs(), and PatternGenerator::updateStat().

TH1I* PatternOptimizerBase::simMuFoundByOmtfPt

protected

Definition at line 53 of file PatternOptimizerBase.h.

Referenced by observeEventEnd(), PatternOptimizerBase(), and savePatternsInRoot().

TH1I* PatternOptimizerBase::simMuPt

protected

Definition at line 52 of file PatternOptimizerBase.h.

Referenced by PatternOptimizerBase().

TH1F* PatternOptimizerBase::simMuPtSpectrum

protected

Definition at line 55 of file PatternOptimizerBase.h.

Referenced by observeEventEnd(), PatternOptimizerBase(), and savePatternsInRoot().

bool PatternOptimizerBase::writeLayerStat = false

protected

Definition at line 57 of file PatternOptimizerBase.h.

Referenced by PatternGenerator::endJob(), and savePatternsInRoot().