HGCalTBRawToDigi.cc

Go to the documentation of this file.

#include <iostream>
#include "HGCal/RawToDigi/plugins/HGCalTBRawToDigi.h"
#include "HGCal/Geometry/interface/HGCalTBGeometryParameters.h"
using namespace std;

unsigned int gray_to_binary (unsigned int gray);
int counter=0;
HGCalTBRawToDigi::HGCalTBRawToDigi(edm::ParameterSet const& conf):
    dataTag_(conf.getParameter<edm::InputTag>("InputLabel")),
    fedId_(conf.getUntrackedParameter<int>("fedId")),
    mapfile_(conf.getUntrackedParameter<std::string>("electronicsMap"))
{
    produces<SKIROC2DigiCollection>();
    consumes<FEDRawDataCollection>(dataTag_);
}

void HGCalTBRawToDigi::beginJob()
{
    ///\todo this should be done by an ESProducer
    HGCalCondObjectTextIO io(0); // don't need a numbering scheme for this

    edm::FileInPath fip(mapfile_);
    if (!io.load(fip.fullPath(), essource_.emap_)) {
        throw cms::Exception("Unable to load electronics map");
    }
}

void HGCalTBRawToDigi::produce(edm::Event& e, const edm::EventSetup& c)
{
    // Step A: Get Inputs
    edm::Handle<FEDRawDataCollection> rawraw;
    e.getByLabel(dataTag_, rawraw);

    std::auto_ptr<SKIROC2DigiCollection> digis(0);
    //
    const FEDRawData& fed = rawraw->FEDData(fedId_);
    if (fed.size() != 0) { /// \todo Exception if 0????
         int ski_up = 2;
         int ski_down = 1;
        // we can figure out the number of samples from the size of the raw data
        int nsamples = fed.size() / (sizeof(uint16_t) * SKIROC::NCHANNELS * 2); // 2 is for ADC and TDC
        digis = std::auto_ptr<SKIROC2DigiCollection>(new SKIROC2DigiCollection(nsamples));
        const uint16_t* pdata = (const uint16_t*)(fed.data());

        // we start from the back...
        int ptr = fed.size() / sizeof(uint16_t) - 1;
/*
        printf("Starting on SKIROC %x\n", pdata[ptr]);
        ptr--; // now we are pointing at a relatively-useless header word
        ptr--; // now we are pointing at the first TDC word
*/
         counter++;
         if(counter > 150){
            ski_up = 4;
            ski_down = 3;
           }  
 
         if(counter > 150*2 ){
            ski_up = 6;
            ski_down = 5;  
           }
         if(counter > 150*3){
            ski_up = 8;
            ski_down = 7;
           }
         if(counter > 150*4){
            ski_up = 10;
            ski_down = 9;
           }
         if(counter > 150*5){
            ski_up = 12;
            ski_down = 11;
           }
         if(counter > 150*6){
            ski_up = 14;
            ski_down = 13;
           }
         if(counter > 150*7){
            ski_up = 16;
            ski_down = 15;
           }
         if(counter > 150*8 ){
            ski_up = 18;
            ski_down = 17;
           }
         if(counter > 150*9){
            ski_up = 20;
            ski_down = 19;
           }
         if(counter > 150*10){
            ski_up = 22;
            ski_down = 21;
           }
         if(counter > 150*11 ){
            ski_up = 24;
            ski_down = 23;
           }
         if(counter > 150*12){
            ski_up = 26;
            ski_down = 25;
           }
     if(counter > 150*13){
            ski_up = 28;
        ski_down = 27;      
       }
         if(counter > 150*14){
            ski_up = 30;
            ski_down = 29;
           }
         if(counter > 150*15){
            ski_up = 32;
            ski_down = 31;
           }

         if(counter == 150*MAXLAYERS) counter = 0;   
//                for (int ski = 2; ski >= 1; ski--) {
            for (int ichan = 0; ichan < SKIROC::NCHANNELS; ichan++) {
                            for (int ski = ski_down; ski <= ski_up; ski++) {
                HGCalTBElectronicsId eid(ski, ichan);
                if (!essource_.emap_.existsEId(eid.rawId())) {
//                  std::cout << "We do not have a mapping for " << eid;
                } else {
                    HGCalTBDetId did = essource_.emap_.eid2detId(eid);
                    digis->addDataFrame(did);
                                        if(ski <=32){
                                    int ptradc1 = ptr - ichan*2 - (ski - ski_down);
                    int ptradc2 = ptr - ichan*2 - (ski - ski_down) -  128;
                                        digis->backDataFrame().setSample(0, gray_to_binary(pdata[ptradc1] & 0xFFF),gray_to_binary( pdata[ptradc2] & 0xFFF),0); 
//                                        cout<<endl<<dec<<"SKI= "<<ski<<" chan= "<<ichan<<" "<<ptradc1<<" "<<ptradc2<<" "<<" High= "<<hex<<(pdata[ptradc1])<<dec<<"  "<<gray_to_binary(pdata[ptradc1] & 0xFFF)<<" Low= "<<hex<<( pdata[ptradc2] & 0xFFF)<<"  "<<dec<<gray_to_binary( pdata[ptradc2] & 0xFFF)<<endl; 
                                          }

/*
                                        else{
                                              int ptradc1 = ptr - ichan*2 + (ski_up - ski) ;
                                              int ptradc2 = ptr - ichan*2 + (ski_up - ski) -  128;
//                      cout<<endl<<dec<<" Layer= "<<did.layer()<<" SKI= "<<eid.iskiroc()<<" chan= "<<ichan<<" "<<ptradc1<<" "<<ptradc2<<" "<<" High= "<<hex<<(pdata[ptradc1])<<dec<<"  "<<gray_to_binary(pdata[ptradc1] & 0xFFF)<<" Low= "<<hex<<( pdata[ptradc2] & 0xFFF)<<"  "<<dec<<gray_to_binary( pdata[ptradc2] & 0xFFF)<<endl;
                                              digis->backDataFrame().setSample(0, gray_to_binary(pdata[ptradc1] & 0xFFF),gray_to_binary( pdata[ptradc2] & 0xFFF),0); 

                                             }  
*/ 

                }
                             }//loop over skirocs
            }
//                 }

    }// fed size > 0

    // put it into the event
    e.put(digis);
}

unsigned int gray_to_binary (unsigned int gray){

  unsigned int result = gray & (1 << 11);
  result |= (gray ^ (result >> 1)) & (1 << 10);
  result |= (gray ^ (result >> 1)) & (1 << 9);
  result |= (gray ^ (result >> 1)) & (1 << 8);
  result |= (gray ^ (result >> 1)) & (1 << 7);
  result |= (gray ^ (result >> 1)) & (1 << 6);
  result |= (gray ^ (result >> 1)) & (1 << 5);
  result |= (gray ^ (result >> 1)) & (1 << 4);
  result |= (gray ^ (result >> 1)) & (1 << 3);
  result |= (gray ^ (result >> 1)) & (1 << 2);
  result |= (gray ^ (result >> 1)) & (1 << 1);
  result |= (gray ^ (result >> 1)) & (1 << 0);
  return result;
}


#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/Framework/interface/MakerMacros.h"

DEFINE_FWK_MODULE(HGCalTBRawToDigi);