EEHltTask.cc

Go to the documentation of this file.

/*
 * \file EEHltTask.cc
 *
 * \author G. Della Ricca
 *
*/

#include <iostream>
#include <fstream>

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"

#include "DQMServices/Core/interface/MonitorElement.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "DataFormats/FEDRawData/interface/FEDRawData.h"

#include "DataFormats/EcalDetId/interface/EEDetId.h"

#include "DQM/EcalEndcapMonitorTasks/interface/EEHltTask.h"

EEHltTask::EEHltTask(const edm::ParameterSet& ps){

  init_ = false;

  initGeometry_ = false;

  dqmStore_ = edm::Service<DQMStore>().operator->();

  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");

  folderName_ = ps.getUntrackedParameter<std::string>("folderName", "FEDIntegrity");

  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);

  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);

  EEDetIdCollection0_ =  consumes<EEDetIdCollection>(ps.getParameter<edm::InputTag>("EEDetIdCollection0"));
  EEDetIdCollection1_ =  consumes<EEDetIdCollection>(ps.getParameter<edm::InputTag>("EEDetIdCollection1"));
  EEDetIdCollection2_ =  consumes<EEDetIdCollection>(ps.getParameter<edm::InputTag>("EEDetIdCollection2"));
  EEDetIdCollection3_ =  consumes<EEDetIdCollection>(ps.getParameter<edm::InputTag>("EEDetIdCollection3"));
  EcalElectronicsIdCollection1_ = consumes<EcalElectronicsIdCollection>(ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection1"));
  EcalElectronicsIdCollection2_ = consumes<EcalElectronicsIdCollection>(ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection2"));
  EcalElectronicsIdCollection3_ = consumes<EcalElectronicsIdCollection>(ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection3"));
  EcalElectronicsIdCollection4_ = consumes<EcalElectronicsIdCollection>(ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection4"));
  EcalElectronicsIdCollection5_ = consumes<EcalElectronicsIdCollection>(ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection5"));
  EcalElectronicsIdCollection6_ = consumes<EcalElectronicsIdCollection>(ps.getParameter<edm::InputTag>("EcalElectronicsIdCollection6"));
  FEDRawDataCollection_ = consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("FEDRawDataCollection"));

  meEEFedsOccupancy_ = 0;
  meEEFedsSizeErrors_ = 0;
  meEEFedsIntegrityErrors_ = 0;

  map = 0;

}

EEHltTask::~EEHltTask(){

}

void EEHltTask::beginJob(void){

  ievt_ = 0;

  if ( dqmStore_ ) {
    dqmStore_->setCurrentFolder(prefixME_ + "/" + folderName_);
    dqmStore_->rmdir(prefixME_ + "/" + folderName_);
  }

}

void EEHltTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {

  initGeometry(c);

  if ( ! mergeRuns_ ) this->reset();

}

void EEHltTask::endRun(const edm::Run& r, const edm::EventSetup& c) {

}

void EEHltTask::reset(void) {

  if ( meEEFedsOccupancy_ ) meEEFedsOccupancy_->Reset();
  if ( meEEFedsSizeErrors_ ) meEEFedsSizeErrors_->Reset();
  if ( meEEFedsIntegrityErrors_ ) meEEFedsIntegrityErrors_->Reset();

}

void EEHltTask::setup(void){

  init_ = true;

  std::string name;

  if ( dqmStore_ ) {
    dqmStore_->setCurrentFolder(prefixME_ + "/" + folderName_);

    name = "FEDEntries";
    meEEFedsOccupancy_ = dqmStore_->book1D(name, name, 54, 601, 655);

    name = "FEDFatal";
    meEEFedsSizeErrors_ = dqmStore_->book1D(name, name, 54, 601, 655);

    name = "FEDNonFatal";
    meEEFedsIntegrityErrors_ = dqmStore_->book1D(name, name, 54, 601, 655);

  }

}

void EEHltTask::cleanup(void){

  if ( ! init_ ) return;

  if ( dqmStore_ ) {
    dqmStore_->setCurrentFolder(prefixME_ + "/" + folderName_);

    if ( meEEFedsOccupancy_ ) dqmStore_->removeElement( meEEFedsOccupancy_->getName() );
    meEEFedsOccupancy_ = 0;

    if ( meEEFedsSizeErrors_ ) dqmStore_->removeElement( meEEFedsSizeErrors_->getName() );
    meEEFedsSizeErrors_ = 0;

    if ( meEEFedsIntegrityErrors_ ) dqmStore_->removeElement( meEEFedsIntegrityErrors_->getName() );
    meEEFedsIntegrityErrors_ = 0;

  }

  init_ = false;

}

void EEHltTask::endJob(void){

  edm::LogInfo("EEHltTask") << "analyzed " << ievt_ << " events";

  if ( enableCleanup_ ) this->cleanup();

}

void EEHltTask::analyze(const edm::Event& e, const edm::EventSetup& c){

  if ( ! init_ ) this->setup();

  ievt_++;

  // ECAL endcap FEDs
  int EEFirstFED[2];
  EEFirstFED[0] = 601; // EE-
  EEFirstFED[1] = 646; // EE+

  int FedsSizeErrors[18];
  for ( int i=0; i<18; i++ ) FedsSizeErrors[i]=0;

  edm::Handle<EEDetIdCollection> ids0;

  if ( e.getByToken(EEDetIdCollection0_, ids0) ) {

    for ( EEDetIdCollection::const_iterator idItr = ids0->begin(); idItr != ids0->end(); ++idItr ) {

      int ism = iSM( *idItr );

      if ( ism > -1 ) FedsSizeErrors[ism-1]++;

    }

  } else {

//    edm::LogWarning("EEHltTask") << "EEDetIdCollection0 not available";

  }

  edm::Handle<FEDRawDataCollection> allFedRawData;

  if ( e.getByToken(FEDRawDataCollection_, allFedRawData) ) {

    for(int zside=0; zside<2; zside++) {

      int firstFedOnSide=EEFirstFED[zside];

      for ( int ism=1; ism<=9; ism++ ) {

    const FEDRawData& fedData = allFedRawData->FEDData( firstFedOnSide + ism - 1 );

    int length = fedData.size()/sizeof(uint64_t);

    if ( length > 0 ) {

      if ( meEEFedsOccupancy_ ) meEEFedsOccupancy_->Fill( firstFedOnSide + ism - 1 );

      uint64_t * pData = (uint64_t *)(fedData.data());
      uint64_t * fedTrailer = pData + (length - 1);
      bool crcError = (*fedTrailer >> 2 ) & 0x1;

      if (crcError) FedsSizeErrors[ism-1]++;

    }

      }

    }

  } else {
    edm::LogWarning("EEHltTask") << "FEDRawDataCollection not available";
  }


  for( int ism=1; ism<=18; ism++ ) {

    if ( FedsSizeErrors[ism-1] != 0 ) {

      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( meEEFedsSizeErrors_ ) meEEFedsSizeErrors_->Fill( fednumber );

    }

  }


  // Integrity errors
  edm::Handle<EEDetIdCollection> ids1;

  if ( e.getByToken(EEDetIdCollection1_, ids1) ) {

    for ( EEDetIdCollection::const_iterator idItr = ids1->begin(); idItr != ids1->end(); ++idItr ) {

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EEDetIdCollection1 not available";

  }

  edm::Handle<EEDetIdCollection> ids2;

  if ( e.getByToken(EEDetIdCollection2_, ids2) ) {

    for ( EEDetIdCollection::const_iterator idItr = ids2->begin(); idItr != ids2->end(); ++idItr ) {

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EEDetIdCollection2 not available";

  }

  edm::Handle<EEDetIdCollection> ids3;

  if ( e.getByToken(EEDetIdCollection3_, ids3) ) {

    for ( EEDetIdCollection::const_iterator idItr = ids3->begin(); idItr != ids3->end(); ++idItr ) {

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EEDetIdCollection3 not available";

  }

  edm::Handle<EcalElectronicsIdCollection> ids4;

  if ( e.getByToken(EcalElectronicsIdCollection1_, ids4) ) {

    for ( EcalElectronicsIdCollection::const_iterator idItr = ids4->begin(); idItr != ids4->end(); ++idItr ) {

      if ( subDet( *idItr ) != EcalEndcap ) continue;

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./34.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EcalElectronicsIdCollection1 not available";

  }

  edm::Handle<EcalElectronicsIdCollection> ids5;

  if ( e.getByToken(EcalElectronicsIdCollection2_, ids5) ) {

    for ( EcalElectronicsIdCollection::const_iterator idItr = ids5->begin(); idItr != ids5->end(); ++idItr ) {

      if ( subDet( *idItr ) != EcalEndcap ) continue;

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EcalElectronicsIdCollection2 not available";

  }

  edm::Handle<EcalElectronicsIdCollection> ids6;

  if ( e.getByToken(EcalElectronicsIdCollection3_, ids6) ) {

    for ( EcalElectronicsIdCollection::const_iterator idItr = ids6->begin(); idItr != ids6->end(); ++idItr ) {

      if ( subDet( *idItr ) != EcalEndcap ) continue;

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./34.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EcalElectronicsIdCollection3 not available";

  }

  edm::Handle<EcalElectronicsIdCollection> ids7;

  if ( e.getByToken(EcalElectronicsIdCollection4_, ids7) ) {

    for ( EcalElectronicsIdCollection::const_iterator idItr = ids7->begin(); idItr != ids7->end(); ++idItr ) {

      if ( subDet( *idItr ) != EcalEndcap ) continue;

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EcalElectronicsIdCollection4 not available";

  }

  edm::Handle<EcalElectronicsIdCollection> ids8;

  if ( e.getByToken(EcalElectronicsIdCollection5_, ids8) ) {

    for ( EcalElectronicsIdCollection::const_iterator idItr = ids8->begin(); idItr != ids8->end(); ++idItr ) {

      if ( subDet( *idItr ) != EcalEndcap ) continue;

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EcalElectronicsIdCollection5 not available";

  }

  edm::Handle<EcalElectronicsIdCollection> ids9;

  if ( e.getByToken(EcalElectronicsIdCollection6_, ids9) ) {

    for ( EcalElectronicsIdCollection::const_iterator idItr = ids9->begin(); idItr != ids9->end(); ++idItr ) {

      if ( subDet( *idItr ) != EcalEndcap ) continue;

      int ism = iSM( *idItr );
      int fednumber = ( ism < 10 ) ? 600 + ism : 636 + ism;

      if ( ism > -1 ) meEEFedsIntegrityErrors_->Fill( fednumber, 1./850.);

    }

  } else {

    edm::LogWarning("EEHltTask") << "EcalElectronicsIdCollection6 not available";

  }

}

//-------------------------------------------------------------------------

void EEHltTask::initGeometry( const edm::EventSetup& setup ) {

  if( initGeometry_ ) return;

  initGeometry_ = true;

  edm::ESHandle< EcalElectronicsMapping > handle;
  setup.get< EcalMappingRcd >().get(handle);
  map = handle.product();

  if( ! map ) edm::LogWarning("EEHltTask") << "EcalElectronicsMapping not available";

}

int EEHltTask::iSM( const EEDetId& id ) {

  if( ! map ) return -1;

  EcalElectronicsId eid = map->getElectronicsId(id);
  int idcc = eid.dccId();

  // EE-
  if( idcc >=  1 && idcc <=  9 ) return( idcc );

  // EE+
  if( idcc >= 46 && idcc <= 54 ) return( idcc - 45 + 9 );

  edm::LogWarning("EEHltTask") << "Wrong DCC id: dcc = " << idcc;
  return -1;

}

int EEHltTask::iSM( const EcalElectronicsId& id ) {

  int idcc = id.dccId();

  // EE-
  if( idcc >=  1 && idcc <=  9 ) return( idcc );

  // EE+
  if( idcc >= 46 && idcc <= 54 ) return( idcc - 45 + 9 );

  edm::LogWarning("EEHltTask") << "Wrong DCC id: dcc = " << idcc;
  return -1;

}