/data/refman/pasoursint/CMSSW_6_1_2_SLHC4_patch1/src/DQM/EcalBarrelMonitorTasks/src/TestPulseTask.cc

Go to the documentation of this file.

#include "../interface/TestPulseTask.h"

#include <algorithm>

#include "DataFormats/EcalRawData/interface/EcalDCCHeaderBlock.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDigi/interface/EcalDataFrame.h"

#include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"

namespace ecaldqm {

  TestPulseTask::TestPulseTask(const edm::ParameterSet &_params, const edm::ParameterSet& _paths) :
    DQWorkerTask(_params, _paths, "TestPulseTask"),
    MGPAGains_(),
    MGPAGainsPN_()
  {
    using namespace std;

    collectionMask_ = 
      (0x1 << kEBDigi) |
      (0x1 << kEEDigi) |
      (0x1 << kPnDiodeDigi) |
      (0x1 << kEBUncalibRecHit) |
      (0x1 << kEEUncalibRecHit);

    edm::ParameterSet const& commonParams(_params.getUntrackedParameterSet("Common"));
    MGPAGains_ = commonParams.getUntrackedParameter<std::vector<int> >("MGPAGains");
    MGPAGainsPN_ = commonParams.getUntrackedParameter<std::vector<int> >("MGPAGainsPN");

    for(int idcc(0); idcc < 54; idcc++){
      enable_[idcc] = false;
      gain_[idcc] = 12;
    }

    for(std::vector<int>::iterator gainItr(MGPAGains_.begin()); gainItr != MGPAGains_.end(); ++gainItr)
      if(*gainItr != 1 && *gainItr != 6 && *gainItr != 12) throw cms::Exception("InvalidConfiguration") << "MGPA gain" << std::endl;

    for(std::vector<int>::iterator gainItr(MGPAGainsPN_.begin()); gainItr != MGPAGainsPN_.end(); ++gainItr)
      if(*gainItr != 1 && *gainItr != 16) throw cms::Exception("InvalidConfiguration") << "PN diode gain" << std::endl; 

    map<string, string> replacements;
    stringstream ss;

    for(vector<int>::iterator gainItr(MGPAGains_.begin()); gainItr != MGPAGains_.end(); ++gainItr){
      ss.str("");
      ss << *gainItr;
      replacements["gain"] = ss.str();

      unsigned offset(0);
      switch(*gainItr){
      case 1: offset = 0; break;
      case 6: offset = 1; break;
      case 12: offset = 2; break;
      default: break;
      }

      MEs_[kOccupancy + offset]->name(replacements);
      MEs_[kShape + offset]->name(replacements);
      MEs_[kAmplitude + offset]->name(replacements);
    }

    for(vector<int>::iterator gainItr(MGPAGainsPN_.begin()); gainItr != MGPAGainsPN_.end(); ++gainItr){
      ss.str("");
      ss << *gainItr;
      replacements["pngain"] = ss.str();

      unsigned offset(0);
      switch(*gainItr){
      case 1: offset = 0; break;
      case 16: offset = 1; break;
      default: break;
      }

      MEs_[kPNOccupancy + offset]->name(replacements);
      MEs_[kPNAmplitude + offset]->name(replacements);
    }
  }

  TestPulseTask::~TestPulseTask()
  {
  }

  void
  TestPulseTask::bookMEs()
  {
    for(std::vector<int>::iterator gainItr(MGPAGains_.begin()); gainItr != MGPAGains_.end(); ++gainItr){
      unsigned offset(0);
      switch(*gainItr){
      case 1: offset = 0; break;
      case 6: offset = 1; break;
      case 12: offset = 2; break;
      default: break;
      }

      MEs_[kOccupancy + offset]->book();
      MEs_[kShape + offset]->book();
      MEs_[kAmplitude + offset]->book();
    }
    for(std::vector<int>::iterator gainItr(MGPAGainsPN_.begin()); gainItr != MGPAGainsPN_.end(); ++gainItr){
      unsigned offset(0);
      switch(*gainItr){
      case 1: offset = 0; break;
      case 16: offset = 1; break;
      default: break;
      }

      MEs_[kPNOccupancy + offset]->book();
      MEs_[kPNAmplitude + offset]->book();
    }
  }

  void
  TestPulseTask::beginRun(const edm::Run&, const edm::EventSetup&)
  {
    for(int idcc(0); idcc < 54; idcc++){
      enable_[idcc] = false;
      gain_[idcc] = 0;
    }
  }

  void
  TestPulseTask::endEvent(const edm::Event&, const edm::EventSetup&)
  {
    for(int idcc(0); idcc < 54; idcc++){
      enable_[idcc] = false;
      gain_[idcc] = 0;
    }
  }

  bool
  TestPulseTask::filterRunType(const std::vector<short>& _runType)
  {
    bool enable(false);

    for(int iFED(0); iFED < 54; iFED++){
      if(_runType[iFED] == EcalDCCHeaderBlock::TESTPULSE_MGPA ||
         _runType[iFED] == EcalDCCHeaderBlock::TESTPULSE_GAP){
        enable = true;
        enable_[iFED] = true;
      }
    }

    return enable;
  }

  void
  TestPulseTask::runOnDigis(const EcalDigiCollection &_digis)
  {
    for(EcalDigiCollection::const_iterator digiItr(_digis.begin()); digiItr != _digis.end(); ++digiItr){
      DetId id(digiItr->id());

      int iDCC(dccId(id) - 1);

      if(!enable_[iDCC]) continue;

      // EcalDataFrame is not a derived class of edm::DataFrame, but can take edm::DataFrame in the constructor
      EcalDataFrame dataFrame(*digiItr);

      unsigned offset(0);
      switch(dataFrame.sample(0).gainId()){
      case 1: offset = 2; gain_[iDCC] = 12; break;
      case 2: offset = 1; gain_[iDCC] = 6; break;
      case 3: offset = 0; gain_[iDCC] = 1; break;
      default: continue;
      }

      if(std::find(MGPAGains_.begin(), MGPAGains_.end(), gain_[iDCC]) == MGPAGains_.end()) continue;

      MEs_[kOccupancy + offset]->fill(id);

      for(int iSample(0); iSample < 10; iSample++)
        MEs_[kShape + offset]->fill(id, iSample + 0.5, float(dataFrame.sample(iSample).adc()));
    }
  }

  void
  TestPulseTask::runOnPnDigis(const EcalPnDiodeDigiCollection &_digis)
  {
    for(EcalPnDiodeDigiCollection::const_iterator digiItr(_digis.begin()); digiItr != _digis.end(); ++digiItr){
      EcalPnDiodeDetId id(digiItr->id());

      int iDCC(dccId(id) - 1);

      if(!enable_[iDCC]) continue;

      unsigned offset(0);
      int gain(0);
      switch(digiItr->sample(0).gainId()){
      case 0: offset = 0; gain = 1; break;
      case 1: offset = 1; gain = 16; break;
      default: continue;
      }

      if(std::find(MGPAGainsPN_.begin(), MGPAGainsPN_.end(), gain) == MGPAGainsPN_.end()) continue;

      MEs_[kPNOccupancy + offset]->fill(id);

      float pedestal(0.);
      for(int iSample(0); iSample < 4; iSample++)
        pedestal += digiItr->sample(iSample).adc();
      pedestal /= 4.;

      float max(0.);
      for(int iSample(0); iSample < 50; iSample++)
        if(digiItr->sample(iSample).adc() > max) max = digiItr->sample(iSample).adc();

      float amplitude(max - pedestal);

      MEs_[kPNAmplitude + offset]->fill(id, amplitude);
    }
  }

  void
  TestPulseTask::runOnUncalibRecHits(const EcalUncalibratedRecHitCollection &_uhits)
  {
    for(EcalUncalibratedRecHitCollection::const_iterator uhitItr(_uhits.begin()); uhitItr != _uhits.end(); ++uhitItr){
      DetId id(uhitItr->id());

      int iDCC(dccId(id) - 1);

      if(!enable_[iDCC]) continue;

      unsigned offset(0);
      switch(gain_[iDCC]){
      case 1: offset = 0; break;
      case 6: offset = 1; break;
      case 12: offset = 2; break;
      default: continue;
      }

      MEs_[kAmplitude + offset]->fill(id, uhitItr->amplitude());
    }
  }

  /*static*/
  void
  TestPulseTask::setMEData(std::vector<MEData>& _data)
  {
    BinService::AxisSpecs axis;
    axis.nbins = 10;
    axis.low = 0.;
    axis.high = 10.;

    for(unsigned iGain(0); iGain < nGain; iGain++){
      _data[kOccupancy + iGain] = MEData("Occupancy", BinService::kEcal2P, BinService::kSuperCrystal, MonitorElement::DQM_KIND_TH2F);
      _data[kShape + iGain] = MEData("Shape", BinService::kSM, BinService::kSuperCrystal, MonitorElement::DQM_KIND_TPROFILE2D, 0, &axis);
      _data[kAmplitude + iGain] = MEData("Amplitude", BinService::kSM, BinService::kCrystal, MonitorElement::DQM_KIND_TPROFILE2D);
    }
    for(unsigned iPNGain(0); iPNGain < nPNGain; iPNGain++){
      _data[kPNOccupancy + iPNGain] = MEData("PNOccupancy", BinService::kEcalMEM2P, BinService::kCrystal, MonitorElement::DQM_KIND_TH2F);
      _data[kPNAmplitude + iPNGain] = MEData("PNPedestal", BinService::kSMMEM, BinService::kCrystal, MonitorElement::DQM_KIND_TPROFILE);
    }
  }

  DEFINE_ECALDQM_WORKER(TestPulseTask);
}