HIPMonitorConfig.cc

Go to the documentation of this file.

#include <Alignment/HIPAlignmentAlgorithm/interface/HIPMonitorConfig.h>

HIPMonitorConfig::HIPMonitorConfig(const edm::ParameterSet& cfg)
    : cfgMonitor(cfg.getParameter<edm::ParameterSet>("monitorConfig")),
      outfilecore(cfgMonitor.getParameter<std::string>("outfile")),
      maxEventsPerJob(cfgMonitor.getParameter<int>("maxEventsPerJob")),
      fillTrackMonitoring(cfgMonitor.getParameter<bool>("fillTrackMonitoring")),
      maxTracks(cfgMonitor.getParameter<int>("maxTracks")),
      trackmonitorvars(maxEventsPerJob, maxTracks),
      fillTrackHitMonitoring(cfgMonitor.getParameter<bool>("fillTrackHitMonitoring")),
      maxHits(cfgMonitor.getParameter<int>("maxHits")),
      hitmonitorvars(maxHits),
      eventCounter(0),
      hitCounter(0) {
  outfile = cfg.getParameter<std::string>("outpath") + outfilecore;
}

HIPMonitorConfig::HIPMonitorConfig(const HIPMonitorConfig& other)
    : cfgMonitor(other.cfgMonitor),
      outfilecore(other.outfilecore),
      maxEventsPerJob(other.maxEventsPerJob),
      fillTrackMonitoring(other.fillTrackMonitoring),
      maxTracks(other.maxTracks),
      fillTrackHitMonitoring(other.fillTrackHitMonitoring),
      maxHits(other.maxHits),
      outfile(other.outfile),
      eventCounter(other.eventCounter),
      hitCounter(other.hitCounter) {}

bool HIPMonitorConfig::checkNevents() {
  bool res = (maxEventsPerJob < 0 || maxEventsPerJob > eventCounter);
  eventCounter++;
  return res;
}
bool HIPMonitorConfig::checkNhits() {
  bool res = (maxHits < 0 || maxHits > hitCounter);
  hitCounter++;
  return res;
}

void HIPTrackMonitorVariables::bookBranches() {
  if (tree != nullptr) {
    tree->Branch("Ntracks", &m_Ntracks);
    tree->Branch("Nhits", &m_Nhits);
    //tree->Branch("DataType", &m_datatype); This is done in HIPAlignmentAlgorithm
    tree->Branch("nhPXB", &m_nhPXB);
    tree->Branch("nhPXF", &m_nhPXF);
    tree->Branch("nhTIB", &m_nhTIB);
    tree->Branch("nhTOB", &m_nhTOB);
    tree->Branch("nhTID", &m_nhTID);
    tree->Branch("nhTEC", &m_nhTEC);
    tree->Branch("Pt", &m_Pt);
    tree->Branch("P", &m_P);
    tree->Branch("Eta", &m_Eta);
    tree->Branch("Phi", &m_Phi);
    tree->Branch("Chi2n", &m_Chi2n);
    tree->Branch("d0", &m_d0);
    tree->Branch("dz", &m_dz);
    tree->Branch("wt", &m_wt);
  }
}
void HIPHitMonitorVariables::bookBranches() {
  if (tree != nullptr) {
    tree->Branch("Id", &m_detId, "Id/i");
    tree->Branch("sinTheta", &m_sinTheta);
    tree->Branch("impactAngle", &m_angle);
    tree->Branch("wt", &m_hitwt);
    tree->Branch("probPresent", &m_hasHitProb);
    tree->Branch("probXY", &m_probXY);
    tree->Branch("probQ", &m_probQ);
    tree->Branch("qualityWord", &m_rawQualityWord);
  }
}
void HIPTrackMonitorVariables::fill() {
  if (tree == nullptr)
    return;
  if (maxNEvents >= 0 && nEvents >= maxNEvents)
    return;

  bool doFill = false;
  m_Ntracks = m_Pt.size();
  if (m_Ntracks == 0)
    return;

  if (maxTracksRcd < 0)
    doFill = true;
  else if (nTracks < maxTracksRcd) {
    if ((nTracks + m_Ntracks) < maxTracksRcd)
      doFill = true;
    else {
      int NewSize = maxTracksRcd - nTracks;
      if ((int)m_Pt.size() < NewSize)
        NewSize = m_Pt.size();

      // Do not touch m_Ntracks, just resize these vectors
      resizeVectors(NewSize);

      doFill = true;
    }
  } else {
    // Record only m_Ntracks
    int tmpNtracks = m_Ntracks;
    resetPerEvent();
    m_Ntracks = tmpNtracks;

    doFill = true;
  }

  if (doFill) {
    tree->Fill();
    nEvents++;
    nTracks += m_Ntracks;
  }

  resetPerEvent();
}
void HIPHitMonitorVariables::fill() {
  if (tree == nullptr)
    return;

  bool doFill = (maxHitsRcd < 0 || nHits < maxHitsRcd);

  if (doFill) {
    tree->Fill();
    nHits++;
  }

  resetPerHit();
}