RctDigiToSourceCardText.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    RctDigiToSourceCardText
// Class:      RctDigiToSourceCardText
//
//
// Original Author:  Alex Tapper
//         Created:  Fri Feb 16 14:52:19 CET 2007
//
//

#include "FWCore/MessageLogger/interface/MessageLogger.h" // Logger
#include "FWCore/ServiceRegistry/interface/Service.h"     // Framework services
#include "RctDigiToSourceCardText.h"

using namespace edm;
using namespace std;

// Set constant
const static int NUM_RCT_CRATES = 18;

RctDigiToSourceCardText::RctDigiToSourceCardText(
    const edm::ParameterSet &iConfig)
    : m_rctInputLabel(iConfig.getParameter<edm::InputTag>("RctInputLabel")),
      m_textFileName(iConfig.getParameter<std::string>("TextFileName")),
      m_nevt(0) {
  // Open the output file
  m_file.open(m_textFileName.c_str(), std::ios::out);

  if (!m_file.good()) {
    throw cms::Exception("RctDigiToSourceCardTextFileOpenError")
        << "RctDigiToSourceCardText::RctDigiToSourceCardText : "
        << " couldn't open the file " << m_textFileName << " for writing"
        << std::endl;
  }

  // Make a SC routing object
  SourceCardRouting m_scRouting;
}

RctDigiToSourceCardText::~RctDigiToSourceCardText() {
  // Close the output file
  m_file.close();
}

// ------------ method called to for each event  ------------
void RctDigiToSourceCardText::analyze(const edm::Event &iEvent,
                                      const edm::EventSetup &iSetup) {

  // get the RCT data
  Handle<L1CaloEmCollection> em;
  Handle<L1CaloRegionCollection> rgn;
  iEvent.getByLabel(m_rctInputLabel, em);
  iEvent.getByLabel(m_rctInputLabel, rgn);

  // General variables
  int RoutingMode;
  unsigned short logicalCardID;
  std::string dataString;

  // Arrays to hold electron variables
  unsigned short eIsoRank[18][4] = {{0}};
  unsigned short eIsoCardId[18][4] = {{0}};
  unsigned short eIsoRegionId[18][4] = {{0}};
  unsigned short eNonIsoRank[18][4] = {{0}};
  unsigned short eNonIsoCardId[18][4] = {{0}};
  unsigned short eNonIsoRegionId[18][4] = {{0}};

  // Fill electrons
  unsigned numIsoEM[18] = {0};
  unsigned numNonIsoEM[18] = {0};

  for (L1CaloEmCollection::const_iterator iem = em->begin(); iem != em->end();
       iem++) {
    if (iem->isolated()) {
      eIsoRank[iem->rctCrate()][numIsoEM[iem->rctCrate()]] = iem->rank();
      eIsoCardId[iem->rctCrate()][numIsoEM[iem->rctCrate()]] = iem->rctCard();
      eIsoRegionId[iem->rctCrate()][numIsoEM[iem->rctCrate()]] =
          iem->rctRegion();
      numIsoEM[iem->rctCrate()]++;
    } else {
      eNonIsoRank[iem->rctCrate()][numNonIsoEM[iem->rctCrate()]] = iem->rank();
      eNonIsoCardId[iem->rctCrate()][numNonIsoEM[iem->rctCrate()]] =
          iem->rctCard();
      eNonIsoRegionId[iem->rctCrate()][numNonIsoEM[iem->rctCrate()]] =
          iem->rctRegion();
      numNonIsoEM[iem->rctCrate()]++;
    }
    // Debug info
    LogDebug("Electrons") << "Rank=" << iem->rank()
                          << " Card=" << iem->rctCard()
                          << " Region=" << iem->rctRegion()
                          << " Crate=" << iem->rctCrate()
                          << " Isolated=" << iem->isolated();
  }

  // Arrays to hold region variables
  unsigned short RC[18][7][2] = {{{0}}};
  unsigned short RCof[18][7][2] = {{{0}}};
  unsigned short RCtau[18][7][2] = {{{0}}};
  unsigned short HF[18][4][2] = {{{0}}};
  unsigned short HFQ[18][4][2] = {{{0}}};
  unsigned short MIPbits[18][7][2] = {{{0}}};
  unsigned short Qbits[18][7][2] = {{{0}}};

  // Fill regions
  for (L1CaloRegionCollection::const_iterator irgn = rgn->begin();
       irgn != rgn->end(); irgn++) {
    if (irgn->id().isHf()) {
      HF[irgn->rctCrate()][irgn->id().rctEta() - 7][irgn->id().rctPhi()] =
          irgn->et();
      HFQ[irgn->rctCrate()][irgn->id().rctEta() - 7][irgn->id().rctPhi()] =
          irgn->fineGrain();
      // Debug info
      LogDebug("HFRegions")
          << "Et=" << irgn->et() << " FineGrain=" << irgn->fineGrain()
          << " Eta=" << irgn->id().rctEta() << " Phi=" << irgn->id().rctPhi()
          << " Crate=" << irgn->rctCrate();
    } else {
      RC[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()] =
          irgn->et();
      RCof[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()] =
          irgn->overFlow();
      RCtau[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()] =
          irgn->tauVeto();
      MIPbits[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()] =
          irgn->mip();
      Qbits[irgn->rctCrate()][irgn->rctCard()][irgn->rctRegionIndex()] =
          irgn->quiet();
      // Debug info
      LogDebug("Regions") << "Et=" << irgn->et()
                          << " OverFlow=" << irgn->overFlow()
                          << " tauVeto=" << irgn->tauVeto()
                          << " mip=" << irgn->mip()
                          << " quiet=" << irgn->quiet()
                          << " Card=" << irgn->rctCard()
                          << " Region=" << irgn->rctRegionIndex()
                          << " Crate=" << irgn->rctCrate();
    }
  }

  for (int crate = 0; crate < NUM_RCT_CRATES; crate++) {

    // Logical Card ID = Source Card number
    RoutingMode = 0;
    m_scRouting.RoutingModetoLogicalCardID(logicalCardID, RoutingMode, crate);

    // Convert electrons to SC format
    m_scRouting.EMUtoSTRING(
        logicalCardID, m_nevt, eIsoRank[crate], eIsoCardId[crate],
        eIsoRegionId[crate], eNonIsoRank[crate], eNonIsoCardId[crate],
        eNonIsoRegionId[crate], MIPbits[crate], Qbits[crate], dataString);

    // Write electrons
    m_file << dataString;

    // Logical Card ID = Source Card number
    RoutingMode = 1;
    m_scRouting.RoutingModetoLogicalCardID(logicalCardID, RoutingMode, crate);

    // Convert regions to SC format
    m_scRouting.RC56HFtoSTRING(logicalCardID, m_nevt, RC[crate], RCof[crate],
                               RCtau[crate], HF[crate], HFQ[crate], dataString);

    // Write regions
    m_file << dataString;

    // Logical Card ID = Source Card number
    RoutingMode = 2;
    m_scRouting.RoutingModetoLogicalCardID(logicalCardID, RoutingMode, crate);

    // Convert regions to SC format
    m_scRouting.RC012toSTRING(logicalCardID, m_nevt, RC[crate], RCof[crate],
                              RCtau[crate], dataString);

    // Write regions
    m_file << dataString;

    // This is to 9 only as this is the shared source card
    if (crate < 9) {
      // Logical Card ID = Source Card number
      RoutingMode = 3;
      m_scRouting.RoutingModetoLogicalCardID(logicalCardID, RoutingMode, crate);

      // Convert regions to SC format
      m_scRouting.RC234toSTRING(logicalCardID, m_nevt, RC[crate], RCof[crate],
                                RCtau[crate], RC[crate + 9], RCof[crate + 9],
                                RCtau[crate + 9], dataString);

      // Write regions
      m_file << dataString;
    }
  }

  // Force write to file
  m_file << flush;

  m_nevt++;
}