BTransitionAnalyzer.h

Go to the documentation of this file.

#ifndef BTRANSITIONANALYZER_H
#define BTRANSITIONANALYZER_H

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CondCore/CondDB/interface/Utils.h"
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
#include "CondFormats/RunInfo/interface/RunInfo.h"
#include "CondFormats/DataRecord/interface/RunSummaryRcd.h"

namespace cond {
  template<class T, class R>
  class BTransitionAnalyzer: public edm::one::EDAnalyzer<edm::one::WatchRuns> {
  public:
    BTransitionAnalyzer( const edm::ParameterSet& pset ):
      m_currentThreshold( pset.getUntrackedParameter<double>( "currentThreshold", 18000. ) ) {
    }
#ifdef __INTEL_COMPILER
    virtual ~BTransitionAnalyzer() = default;
#endif
    // implicit copy constructor
    // implicit assignment operator
    // implicit destructor
    void beginJob() final {};
    void beginRun(edm::Run const&, edm::EventSetup const&) final {};
    void analyze(edm::Event const&, edm::EventSetup const&) final {};
    void endRun(edm::Run const& run , edm::EventSetup const& eventSetup) final {
      edm::ESHandle<RunInfo> runInfoHandle;
      edm::ESHandle<T> payloadHandle, payloadRefHandle;
      eventSetup.get<RunInfoRcd>().get(runInfoHandle);
      double avg_current = (double)runInfoHandle->m_avg_current;
      double current_default = -1;
      std::string bOnLabel = std::string( "38T" );
      std::string bOffLabel = std::string( "0T" );
      std::string bFieldLabel = bOnLabel;
      LogDebug( "BTransitionAnalyzer" ) << "Comparing value of magnet current: " << avg_current << " A for run: " << run.run()
                                        << " with the corresponding threshold: "<< m_currentThreshold << " A." << std::endl;
      if( avg_current != current_default && avg_current <= m_currentThreshold ) bFieldLabel = bOffLabel;
      edm::LogInfo( "BTransitionAnalyzer" ) << "The magnet was " 
                                            << ( bFieldLabel == bOnLabel ? "ON" : "OFF" ) << " during run " << run.run()
                                            << ".\nLoading the product for the corrisponding label " << bFieldLabel << std::endl;
      eventSetup.get<R>().get( bFieldLabel, payloadHandle );
      eventSetup.get<R>().get( payloadRefHandle );
      edm::Service<cond::service::PoolDBOutputService> mydbservice;
      if( mydbservice.isAvailable() ) {
        if( !equalPayloads( payloadHandle, payloadRefHandle ) ) {
          edm::LogInfo( "BTransitionAnalyzer" ) << "Exporting payload corresponding to the calibrations for magnetic field "
                                                << ( bFieldLabel == bOnLabel ? "ON" : "OFF" )
                                                << " starting from run number: " << run.run() << std::endl;
          mydbservice->writeOne( payloadHandle.product(), run.run(), demangledName( typeid(R) ) );
        } else {
          edm::LogInfo( "BTransitionAnalyzer" ) << "The payload corresponding to the calibrations for magnetic field "
                                                << ( bFieldLabel == bOnLabel ? "ON" : "OFF" )
                                                << " is still valid for run " << run.run()
                                                << ".\nNo transfer needed." << std::endl;
        }
      } else {
        edm::LogError( "BTransitionAnalyzer" ) << "PoolDBOutputService unavailable";
      }
    }
    void endJob() final {};
    virtual bool equalPayloads( edm::ESHandle<T> const & payloadHandle, edm::ESHandle<T> const & payloadRefHandle ) = 0;
  private:
    double m_currentThreshold;
  };
} //namespace cond
#endif //BTRANSITIONANALYZER_H