EventWithHistoryProducerFromL1ABC.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    EventWithHistoryProducerFromL1ABC
// Class:      EventWithHistoryProducerFromL1ABC
// 
//
// Original Author:  Andrea Venturi
//         Created:  Tue Jun 30 15:26:20 CET 2009
//
//


// system include files
#include <memory>
#include <map>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/Utilities/interface/InputTag.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/Scalers/interface/L1AcceptBunchCrossing.h"
#include "DPGAnalysis/SiStripTools/interface/EventWithHistory.h"
//
// class decleration
//

class EventWithHistoryProducerFromL1ABC : public edm::EDProducer {
   public:
      explicit EventWithHistoryProducerFromL1ABC(const edm::ParameterSet&);
      ~EventWithHistoryProducerFromL1ABC();

   private:
      virtual void beginJob() ;
      virtual void beginRun(edm::Run&, const edm::EventSetup&) ;
      virtual void produce(edm::Event&, const edm::EventSetup&);
      virtual void endRun(edm::Run&, const edm::EventSetup&) ;
      virtual void endJob() ;
      
      // ----------member data ---------------------------

  edm::InputTag _l1abccollection;
  const bool _forceNoOffset;
  std::map<unsigned int, long long> _offsets;
  long long _curroffset;
  unsigned int _curroffevent;
};

//
// constants, enums and typedefs
//


//
// static data member definitions
//

//
// constructors and destructor
//
EventWithHistoryProducerFromL1ABC::EventWithHistoryProducerFromL1ABC(const edm::ParameterSet& iConfig):
  _l1abccollection(iConfig.getParameter<edm::InputTag>("l1ABCCollection")),
  _forceNoOffset(iConfig.getUntrackedParameter<bool>("forceNoOffset",false)),
  _offsets(), _curroffset(0), _curroffevent(0)
{

  if(_forceNoOffset) edm::LogWarning("NoOffsetComputation") << "Orbit and BX offset will NOT be computed: Be careful!";

  produces<EventWithHistory>();
   
   //now do what ever other initialization is needed
  
}


EventWithHistoryProducerFromL1ABC::~EventWithHistoryProducerFromL1ABC()
{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
EventWithHistoryProducerFromL1ABC::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
   using namespace edm;

   if(iEvent.run() < 110878 ) {

     std::auto_ptr<EventWithHistory> pOut(new EventWithHistory(iEvent));
     iEvent.put(pOut);

   }
   else {

     Handle<L1AcceptBunchCrossingCollection > pIn;
     iEvent.getByLabel(_l1abccollection,pIn);
     
     // offset computation
     
     long long orbitoffset = 0;
     int bxoffset = 0;
     if(!_forceNoOffset) {
       for(L1AcceptBunchCrossingCollection::const_iterator l1abc=pIn->begin();l1abc!=pIn->end();++l1abc) {
     if(l1abc->l1AcceptOffset()==0) {
       orbitoffset = (long long)l1abc->orbitNumber() - (long long)iEvent.orbitNumber();
       bxoffset = l1abc->bunchCrossing() - iEvent.bunchCrossing();
     }
       }
     }
     
     
     std::auto_ptr<EventWithHistory> pOut(new EventWithHistory(iEvent,*pIn,orbitoffset,bxoffset));
     iEvent.put(pOut);
     
     // monitor offset
     
     long long absbxoffset = orbitoffset*3564 + bxoffset;
     
     if(_offsets.size()==0) {
       _curroffset = absbxoffset;
       _curroffevent = iEvent.id().event();
       _offsets[iEvent.id().event()] = absbxoffset;
     }
     else {
       if(_curroffset != absbxoffset || iEvent.id().event() < _curroffevent ) {

     if( _curroffset != absbxoffset) {
       edm::LogInfo("AbsoluteBXOffsetChanged") << "Absolute BX offset changed from " 
                           << _curroffset << " to "
                           << absbxoffset << " at orbit "
                           << iEvent.orbitNumber() << " and BX "
                           << iEvent.bunchCrossing();
       for(L1AcceptBunchCrossingCollection::const_iterator l1abc=pIn->begin();l1abc!=pIn->end();++l1abc) {
         edm::LogVerbatim("AbsoluteBXOffsetChanged") << *l1abc;
       }
     }

     _curroffset = absbxoffset;
     _curroffevent = iEvent.id().event();
     _offsets[iEvent.id().event()] = absbxoffset;
       }
     }
   }
}

// ------------ method called once each job just before starting event loop  ------------
void 
EventWithHistoryProducerFromL1ABC::beginJob()
{
}

void 
EventWithHistoryProducerFromL1ABC::beginRun(edm::Run&, const edm::EventSetup&)
{
  // reset offset vector

  _offsets.clear();
  edm::LogInfo("AbsoluteBXOffsetReset") << "Absolute BX offset map reset";

}

void 
EventWithHistoryProducerFromL1ABC::endRun(edm::Run&, const edm::EventSetup&)
{
  // summary of absolute bx offset vector

  edm::LogInfo("AbsoluteBXOffsetSummary") << "Absolute BX offset summary:";
  for(std::map<unsigned int, long long>::const_iterator offset=_offsets.begin();offset!=_offsets.end();++offset) {
    edm::LogVerbatim("AbsoluteBXOffsetSummary") << offset->first << " " << offset->second;
  }

}

// ------------ method called once each job just after ending the event loop  ------------
void 
EventWithHistoryProducerFromL1ABC::endJob() {
}

//define this as a plug-in
DEFINE_FWK_MODULE(EventWithHistoryProducerFromL1ABC);