EcalDigiToRaw.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    EcalDigiToRaw
// Class:      EcalDigiToRaw
// 
//
// Original Author:  Emmanuelle Perez
//         Created:  Sat Nov 25 13:59:51 CET 2006
// $Id: EcalDigiToRaw.cc,v 1.18 2010/01/13 21:55:15 wmtan Exp $
//
//


// system include files


// user include files
#include "EventFilter/EcalDigiToRaw/interface/EcalDigiToRaw.h"

#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/EcalDetId/interface/EcalDetIdCollections.h"
#include "DataFormats/EcalDigi/interface/EcalSrFlag.h"

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


// #include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"




using namespace edm;
using namespace std;

EcalDigiToRaw::EcalDigiToRaw(const edm::ParameterSet& iConfig)
{

   doTCC_    = iConfig.getUntrackedParameter<bool>("WriteTCCBlock");
   doSR_     = iConfig.getUntrackedParameter<bool>("WriteSRFlags");
   doTower_  = iConfig.getUntrackedParameter<bool>("WriteTowerBlock");

   doBarrel_ = iConfig.getUntrackedParameter<bool>("DoBarrel");
   doEndCap_ = iConfig.getUntrackedParameter<bool>("DoEndCap");

   listDCCId_ = iConfig.getUntrackedParameter< std::vector<int32_t> >("listDCCId");
   label_= iConfig.getParameter<string>("Label");
   instanceNameEB_ = iConfig.getParameter<string>("InstanceEB");
   instanceNameEE_ = iConfig.getParameter<string>("InstanceEE");

   labelTT_ = iConfig.getParameter<edm::InputTag>("labelTT");

   labelEBSR_ = iConfig.getParameter<edm::InputTag>("labelEBSRFlags");
   labelEESR_ = iConfig.getParameter<edm::InputTag>("labelEESRFlags");

   counter_ = 0;
   debug_ = iConfig.getUntrackedParameter<bool>("debug");


   Towerblockformatter_ = new TowerBlockFormatter;
   TCCblockformatter_   = new TCCBlockFormatter();
   SRblockformatter_    = new SRBlockFormatter();
   Headerblockformatter_= new BlockFormatter;

   produces<FEDRawDataCollection>();


}


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

 delete Towerblockformatter_;
 delete TCCblockformatter_;
 delete SRblockformatter_;
 delete Headerblockformatter_;

}


//
// member functions
//

// ------------ method called to for each event  ------------
void
EcalDigiToRaw::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{

   if (debug_) cout << "Enter in EcalDigiToRaw::produce ... " << endl;

   ESHandle< EcalElectronicsMapping > ecalmapping;
   iSetup.get< EcalMappingRcd >().get(ecalmapping);
   const EcalElectronicsMapping* TheMapping = ecalmapping.product();

   Towerblockformatter_ -> StartEvent();
   SRblockformatter_ -> StartEvent();

  runnumber_ = iEvent.id().run();

  // bx_ = (counter_ % BXMAX);
  // orbit_number_ = counter_ / BXMAX;
  // counter_ ++;

  counter_ = iEvent.id().event();
  bx_ = iEvent.bunchCrossing();
  orbit_number_ = iEvent.orbitNumber();

  lv1_ = counter_ % (0x1<<24);

  auto_ptr<FEDRawDataCollection> productRawData(new FEDRawDataCollection);


  Headerblockformatter_ -> DigiToRaw(productRawData.get());


// ---------   Now the Trigger Block part

  Handle<EcalTrigPrimDigiCollection> ecalTrigPrim;
  
  Handle<EBSrFlagCollection> ebSrFlags;
  Handle<EESrFlagCollection> eeSrFlags;


  if (doTCC_) {

     if (debug_) cout << "Creation of the TCC block  " << endl;
     // iEvent.getByType(ecalTrigPrim);
    iEvent.getByLabel(labelTT_, ecalTrigPrim);

     // loop on TP's and add one by one to the block
     for (EcalTrigPrimDigiCollection::const_iterator it = ecalTrigPrim -> begin();
               it != ecalTrigPrim -> end(); it++) {

       const EcalTriggerPrimitiveDigi& trigprim = *it;
       const EcalTrigTowerDetId& detid = it -> id();

       if ( (detid.subDet() == EcalBarrel) && (! doBarrel_) ) continue;
           if ( (detid.subDet() == EcalEndcap) && (! doEndCap_) ) continue;

       int iDCC = TheMapping -> DCCid(detid);
           int FEDid = FEDNumbering::MINECALFEDID + iDCC;

           FEDRawData& rawdata = productRawData.get() -> FEDData(FEDid);
       
       // adding the primitive to the block
       TCCblockformatter_ -> DigiToRaw(trigprim, rawdata, TheMapping);

     }   // end loop on ecalTrigPrim

   }  // endif doTCC


   if (doSR_) { 
    if (debug_) cout << " Process the SR flags " << endl;

    if (doBarrel_) {

        // iEvent.getByType(ebSrFlags);
       iEvent.getByLabel(labelEBSR_, ebSrFlags);
                                                                                                                                                
           for (EBSrFlagCollection::const_iterator it = ebSrFlags -> begin();
                        it != ebSrFlags -> end(); it++) {
                const EcalSrFlag& srflag = *it;
        int flag = srflag.value();

        EcalTrigTowerDetId id = srflag.id();
        int Dccid = TheMapping -> DCCid(id);
        int DCC_Channel = TheMapping -> iTT(id);
        int FEDid = FEDNumbering::MINECALFEDID + Dccid;
        // if (Dccid == 10) cout << "Dcc " << Dccid << " DCC_Channel " << DCC_Channel << " flag " << flag << endl;
        if (debug_) cout << "will process SRblockformatter_ for FEDid " << dec << FEDid << endl;
        FEDRawData& rawdata = productRawData.get() -> FEDData(FEDid);
        if (debug_) Headerblockformatter_ -> print(rawdata);
        SRblockformatter_ -> DigiToRaw(Dccid,DCC_Channel,flag, rawdata);

           }
    }  // end DoBarrel


    if (doEndCap_) {
    // iEvent.getByType(eeSrFlags);
    iEvent.getByLabel(labelEESR_, eeSrFlags);

           for (EESrFlagCollection::const_iterator it = eeSrFlags -> begin();
                        it != eeSrFlags -> end(); it++) {
                const EcalSrFlag& srflag = *it;
                int flag = srflag.value();
        EcalScDetId id = srflag.id();
        pair<int, int> ind = TheMapping -> getDCCandSC(id);
        int Dccid = ind.first;
        int DCC_Channel = ind.second;

                int FEDid = FEDNumbering::MINECALFEDID + Dccid;
                FEDRawData& rawdata = productRawData.get() -> FEDData(FEDid);
                SRblockformatter_ -> DigiToRaw(Dccid,DCC_Channel,flag, rawdata);
           }
    }  // end doEndCap

   }   // endif doSR


// ---------  Now the Tower Block part

  Handle<EBDigiCollection> ebDigis;
  Handle<EEDigiCollection> eeDigis;

  if (doTower_) {

    if (doBarrel_) {
    if (debug_) cout << "Creation of the TowerBlock ... Barrel case " << endl;
        iEvent.getByLabel(label_,instanceNameEB_,ebDigis);
        for (EBDigiCollection::const_iterator it=ebDigis -> begin();
                                it != ebDigis->end(); it++) {
                const EBDataFrame& dataframe = *it;
                const EBDetId& ebdetid = it -> id();
        int DCCid = TheMapping -> DCCid(ebdetid);
            int FEDid = FEDNumbering::MINECALFEDID + DCCid ;
                FEDRawData& rawdata = productRawData.get() -> FEDData(FEDid);
                Towerblockformatter_ -> DigiToRaw(dataframe, rawdata, TheMapping);
        }

    }

    if (doEndCap_) {
    if (debug_) cout << "Creation of the TowerBlock ... EndCap case " << endl;
        iEvent.getByLabel(label_,instanceNameEE_,eeDigis);
        for (EEDigiCollection::const_iterator it=eeDigis -> begin();
                                it != eeDigis->end(); it++) {
                const EEDataFrame& dataframe = *it;
                const EEDetId& eedetid = it -> id();
        EcalElectronicsId elid = TheMapping -> getElectronicsId(eedetid);
                int DCCid = elid.dccId() ;   
                int FEDid = FEDNumbering::MINECALFEDID + DCCid;
                FEDRawData& rawdata = productRawData.get() -> FEDData(FEDid);
                Towerblockformatter_ -> DigiToRaw(dataframe, rawdata, TheMapping);
        }
    }

  }  // endif doTower_



// -------- Clean up things ...

  map<int, map<int,int> >* FEDorder = Towerblockformatter_ -> GetFEDorder();

  Headerblockformatter_ -> CleanUp(productRawData.get(), FEDorder);


/*
   cout << "For FED 633 " << endl;
         FEDRawData& rawdata = productRawData -> FEDData(633);
         Headerblockformatter_ -> print(rawdata);
*/

 // Headerblockformatter_ -> PrintSizes(productRawData.get());



 Towerblockformatter_ -> EndEvent(productRawData.get());

 iEvent.put(productRawData);


 return;

}


// ------------ method called once each job just before starting event loop  ------------
void 
EcalDigiToRaw::beginJob()
{
    Headerblockformatter_ -> SetParam(this);
    Towerblockformatter_  -> SetParam(this);
    TCCblockformatter_  -> SetParam(this);
    SRblockformatter_   -> SetParam(this);
}

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