DTTPDeadWriter.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author S. Bolognesi
 */

#include "CalibMuon/DTCalibration/plugins/DTTPDeadWriter.h"
#include "CalibMuon/DTCalibration/interface/DTCalibDBUtils.h"

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"

#include "CondFormats/DTObjects/interface/DTT0.h"
#include "CondFormats/DataRecord/interface/DTT0Rcd.h"
#include "CondFormats/DataRecord/interface/DTDeadFlagRcd.h"
#include "CondFormats/DTObjects/interface/DTDeadFlag.h"

/* C++ Headers */
#include <vector>
#include <set>
#include <iostream>
#include <fstream>
#include <string>
#include <sstream>
#include "TFile.h"
#include "TH1.h"

using namespace std;
using namespace edm;

// Constructor
DTTPDeadWriter::DTTPDeadWriter(const ParameterSet& pset) : dtGeomToken_(esConsumes<edm::Transition::BeginRun>()) {
  // get selected debug option
  debug = pset.getUntrackedParameter<bool>("debug", false);
  t0Token_ = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", pset.getParameter<string>("debug")));

  // Create the object to be written to DB
  tpDeadList = new DTDeadFlag();

  if (debug)
    cout << "[DTTPDeadWriter]Constructor called!" << endl;
}

// Destructor
DTTPDeadWriter::~DTTPDeadWriter() {
  if (debug)
    cout << "[DTTPDeadWriter]Destructor called!" << endl;
}

void DTTPDeadWriter::beginRun(const edm::Run&, const EventSetup& setup) {
  // Get the t0 map
  tZeroMap = &setup.getData(t0Token_);

  // Get the muon Geometry
  muonGeom = setup.getHandle(dtGeomToken_);
}

// Do the job
void DTTPDeadWriter::analyze(const Event& event, const EventSetup& eventSetup) {
  set<DTLayerId> analyzedLayers;

  //Loop on tzero map
  for (DTT0::const_iterator tzero = tZeroMap->begin(); tzero != tZeroMap->end(); ++tzero) {
    //Consider what layers have been already considered
    // @@@ NEW DTT0 FORMAT
    //    DTLayerId layerId = (DTWireId((*tzero).first.wheelId,
    //                (*tzero).first.stationId,
    //                (*tzero).first.sectorId,
    //                (*tzero).first.slId,
    //                (*tzero).first.layerId,
    //                (*tzero).first.cellId)).layerId();
    int channelId = tzero->channelId;
    if (channelId == 0)
      continue;
    DTLayerId layerId = (DTWireId(channelId)).layerId();
    // @@@ NEW DTT0 END
    if (analyzedLayers.find(layerId) == analyzedLayers.end()) {
      analyzedLayers.insert(layerId);

      //Take the layer topology
      const DTTopology& dtTopo = muonGeom->layer(layerId)->specificTopology();
      const int firstWire = dtTopo.firstChannel();
      //const int lastWire = dtTopo.lastChannel();
      const int nWires = muonGeom->layer(layerId)->specificTopology().channels();

      //Loop on wires
      for (int wire = firstWire; wire <= nWires; wire++) {
        DTWireId wireId(layerId, wire);
        float t0 = 0;
        float t0rms = 0;
        tZeroMap->get(wireId, t0, t0rms, DTTimeUnits::ns);

        //If no t0 stored then is a tp dead channel
        if (!t0) {
          tpDeadList->setCellDead_TP(wireId, true);
          cout << "Wire id " << wireId << " is TP dead" << endl;
        }
      }
    }
  }
}

// Write objects to DB
void DTTPDeadWriter::endJob() {
  if (debug)
    cout << "[DTTPDeadWriter]Writing ttrig object to DB!" << endl;

  // FIXME: to be read from cfg?
  string deadRecord = "DTDeadFlagRcd";

  // Write the object to DB
  DTCalibDBUtils::writeToDB(deadRecord, tpDeadList);
}