TcdsRawToDigi.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    ProdTutorial/TcdsRawToDigi
// Class:      TcdsRawToDigi
// 
//
// Original Author:  Chris Palmer
// Improved by    :  Salvatore Di Guida
//         Created:  Thu, 28 May 2015 19:54:56 GMT
//
//


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

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

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

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

#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDHeader.h"
#include "DataFormats/FEDRawData/interface/FEDTrailer.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "DataFormats/Scalers/interface/BSTRecord.h"

#include "EventFilter/FEDInterface/interface/GlobalEventNumber.h"
#include "EventFilter/FEDInterface/interface/FED1024.h"

using namespace std;

//
// class declaration
//


class TcdsRawToDigi : public edm::stream::EDProducer<> {
   public:
      explicit TcdsRawToDigi(const edm::ParameterSet&);
      ~TcdsRawToDigi();

      static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

   private:
      virtual void produce(edm::Event&, const edm::EventSetup&) override;
      
      edm::EDGetTokenT<FEDRawDataCollection> dataToken_;
      
};

TcdsRawToDigi::TcdsRawToDigi(const edm::ParameterSet& iConfig)
{
    edm::InputTag dataLabel = iConfig.getParameter<edm::InputTag>("InputLabel");
    dataToken_=consumes<FEDRawDataCollection>(dataLabel);
    produces<int>( "nibble" ).setBranchAlias( "nibble");
    produces<int64_t>( "triggerCount" ).setBranchAlias( "triggerCount");
    produces<BSTRecord>( "bstRecord" ).setBranchAlias( "bstRecord" );
}


TcdsRawToDigi::~TcdsRawToDigi()
{
}


//
// member functions
//

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

    edm::Handle<FEDRawDataCollection> rawdata;  
    iEvent.getByToken(dataToken_,rawdata);

    int nibble=-99;
    int triggerCount = -99;
    BSTRecord bstRecord;
    if( rawdata.isValid() ) {
        const FEDRawData& tcdsData = rawdata->FEDData(FEDNumbering::MINTCDSuTCAFEDID);
        if(tcdsData.size()>0){
            evf::evtn::TCDSRecord tcdsRecord(tcdsData.data());
            nibble = (int)tcdsRecord.getHeader().getData().header.nibble;
            triggerCount = (int64_t)tcdsRecord.getHeader().getData().header.triggerCount;
            uint64_t gpsh = tcdsRecord.getBST().getBST().bst.gpstimehigh;
            uint32_t gpsl = tcdsRecord.getBST().getBST().bst.gpstimelow;
            unsigned long long gpstime = static_cast<unsigned long long>((gpsh << 32) + gpsl);
            unsigned char bstMaster = static_cast<unsigned char>(tcdsRecord.getBST().getBST().bst.bstMaster_bireserved16 >> 8);
            uint32_t turnCountHigh = tcdsRecord.getBST().getBST().bst.turnCountHigh;
            uint16_t turnCountLow = tcdsRecord.getBST().getBST().bst.turnCountLow;
            unsigned int turnCount = static_cast<unsigned int>((turnCountHigh << 16) + turnCountLow);
            uint32_t lhcFillHigh = tcdsRecord.getBST().getBST().bst.lhcFillHigh;
            uint16_t lhcFillLow = tcdsRecord.getBST().getBST().bst.lhcFillLow;
            unsigned int lhcFill = static_cast<unsigned int>((lhcFillHigh << 16) + lhcFillLow);
            unsigned short beamMode = static_cast<unsigned short>(tcdsRecord.getBST().getBST().bst.beamMode);
            unsigned char particleType1 = static_cast<unsigned char>(tcdsRecord.getBST().getBST().bst.particleTypes & 0xFF);
            unsigned char particleType2 = static_cast<unsigned char>(tcdsRecord.getBST().getBST().bst.particleTypes >> 8);
            unsigned short beamMomentum = static_cast<unsigned short>(tcdsRecord.getBST().getBST().bst.beamMomentum);
            unsigned int intensityBeam1 = static_cast<unsigned int>(tcdsRecord.getBST().getBST().bst.intensityBeam1);
            unsigned int intensityBeam2 = static_cast<unsigned int>(tcdsRecord.getBST().getBST().bst.intensityBeam2);
            bstRecord.set(gpstime, bstMaster, turnCount, lhcFill, beamMode, particleType1, particleType2, beamMomentum, intensityBeam1, intensityBeam2);
        } else {
            nibble=-2;
            triggerCount=-2;
        }
    } else {
        nibble=-1;
        triggerCount=-1;
    }
    //std::cout<<"nibble is "<<nibble<<std::endl;

    iEvent.put(std::make_unique<int>(nibble), "nibble");
    iEvent.put(std::make_unique<int64_t>(triggerCount), "triggerCount");
    iEvent.put(std::make_unique<BSTRecord>(bstRecord), "bstRecord");
}

 
// ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
void
TcdsRawToDigi::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
    edm::ParameterSetDescription desc;
    desc.add<edm::InputTag>("InputLabel",edm::InputTag("rawDataCollector"));
    descriptions.add("tcdsRawToDigi", desc);
}

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