HcalTBQADCUnpacker.cc

Go to the documentation of this file.

#include "RecoTBCalo/HcalTBObjectUnpacker/interface/HcalTBQADCUnpacker.h"
#include "FWCore/Utilities/interface/Exception.h"

// QADC channels
static const int N_QADCS_ALLOWED = 6;

// Channel to logical unit
// TB04
static const int aScint1          = 52;  // 14x14 cm
static const int aScint2          = 53;  // 4x4 cm
static const int aScint3          = 54;  // 2x2 cm
static const int aScint4          = 66;  // 14x14 cm
static const int aMuonV            = 56;  // Behind the table
static const int aMuonV3           = 64;  // on HB1
static const int aMuonV6           = 62;  // on HB2
static const int aMuonVH1          = 67;  // Oct. muon veto wall
static const int aMuonVH2          = 68;  // Oct. muon veto wall
static const int aMuonVH3          = 69;  // Oct. muon veto wall
static const int aMuonVH4          = 70;  // Oct. muon veto wall
static const int aCerenkov2        = 49;  // el id
static const int aCerenkov3        = 59;  // pi/proton separation
static const int aSCI_VLE          = 65;  // VLE line 
static const int aSCI_521          = 57;
static const int aSCI_528          = 58;
// TB06
static const int bMuonV1           = 0;  // Muon veto wall
static const int bMuonV2           = 1;  // Muon veto wall
static const int bMuonV3           = 2;  // Muon veto wall
static const int bMuonV4           = 3;  // Muon veto wall
static const int bMuonV5           = 4;  // Muon veto wall
static const int bMuonV6           = 5;  // Muon veto wall
static const int bMuonV7           = 6;  // Muon veto wall
static const int bMuonV8           = 7;  // Muon veto wall
static const int bMuonVF           = 10;  // Behind the table
static const int bMuonVB           = 11;  // Behind a beam dump
static const int bScint1           = 12;  // 14x14 cm
static const int bScint2           = 13;  // 4x4 cm
static const int bScint3           = 14;  // 2x2 cm
static const int bScint4           = 15;  // 14x14 cm
static const int bCerenkov1        = 16;  // in HE line
static const int bCerenkov2        = 17;  // el id in VLE
static const int bCerenkov3        = 18;  // pi/proton id
static const int bTOF1S            = 20;  // TOF1S (S=Saleve side, TOF1 upstream)
static const int bTOF1J            = 21;  // TOF1J (J=Jura side)
static const int bTOF2S            = 22;  // TOF2S (TOF2 downstream)
static const int bTOF2J            = 23;  // TOF2J
static const int bSCI_521          = 24;  // In HE beam line for now
static const int bSCI_528          = 25;  // NC
static const int bVH2              = 27;  // beam halo left from particle view
static const int bVH4              = 28;  // beam halo down
static const int bVH3              = 29;  // beam halo up
static const int bVH1              = 30;  // beam halo right from particle view

namespace hcaltb {

HcalTBQADCUnpacker::HcalTBQADCUnpacker(){}

struct ClassicQADCDataFormat {
  unsigned int cdfHeader0,cdfHeader1,cdfHeader2,cdfHeader3;
  unsigned short data[N_QADCS_ALLOWED*32];
  unsigned int cdfTrailer0,cdfTrailer1;
};

 struct CombinedTDCQDCDataFormat {
    unsigned int cdfHeader0,cdfHeader1,cdfHeader2,cdfHeader3;
    unsigned int n_qdc_hits; // Count of QDC channels
    unsigned int n_tdc_hits; // upper/lower TDC counts    
    unsigned short qdc_values[4];
  };

// Sets the pedestal and gain
void HcalTBQADCUnpacker::setCalib(const std::vector<std::vector<std::string> >& calibLines_) {
// The default pedestal and gain
    for(int i=0;i<N_QADCS_ALLOWED*32;i++)
     {
      qdc_ped[i]=0.;qdc_gain[i]=1.;
     }
// Pedestal and gains from configuration file.
    for(unsigned int ii=0;ii<calibLines_.size();ii++)
     {
      if(calibLines_[ii][0]=="QDC")
        {
        if(calibLines_[ii].size()==4)
          {
          int channel=atoi(calibLines_[ii][1].c_str());
          qdc_ped[channel]=atof(calibLines_[ii][2].c_str());
          qdc_gain[channel]=atof(calibLines_[ii][3].c_str());
    //    printf("Got QDC %i ped %f , gain %f\n",channel, qdc_ped[channel],qdc_gain[channel]);
          }
         else
          {
                  throw cms::Exception("Incomplete configuration") << 
          "Wrong QADC configuration format: expected 3 parameters, got "<<calibLines_[ii].size()-1;
          }
        }
     } // End of calibLines.
    }

void HcalTBQADCUnpacker::unpack(const FEDRawData& raw,
                   HcalTBBeamCounters& beamadc, bool is04) const {

  if (raw.size()<3*8) {
    throw cms::Exception("Missing Data") << "No data in the QDC block";
  }


  if(is04){ 
    const ClassicQADCDataFormat* qadc=(const ClassicQADCDataFormat*)raw.data();
    double qdc_calib_hits[N_QADCS_ALLOWED*32];
    // Applying mask, pedestal subtraction and gain.
    for (unsigned int i=0;i<N_QADCS_ALLOWED*32;i++)
      qdc_calib_hits[i]=((qadc->data[i]&0xFFF)-qdc_ped[i])/qdc_gain[i];

    // Ecal energy sum should go here.
    double Ecal7x7=0.;
    for(int i=0;i<49;i++)Ecal7x7+=qdc_calib_hits[i];


    beamadc.setADCs04(qdc_calib_hits[aMuonV],qdc_calib_hits[aMuonV3],qdc_calib_hits[aMuonV6],
            qdc_calib_hits[aMuonVH1],qdc_calib_hits[aMuonVH2],qdc_calib_hits[aMuonVH3],
            qdc_calib_hits[aMuonVH4],qdc_calib_hits[aCerenkov2],qdc_calib_hits[aCerenkov3],
            qdc_calib_hits[aSCI_VLE],qdc_calib_hits[aSCI_521],qdc_calib_hits[aSCI_528],
            qdc_calib_hits[aScint1],qdc_calib_hits[aScint2],qdc_calib_hits[aScint3],
            qdc_calib_hits[aScint4],Ecal7x7);
       }
  else{ 
    const CombinedTDCQDCDataFormat* qdctdc=(const CombinedTDCQDCDataFormat*)raw.data();
    double qdc_calib_hits[32];
    for (unsigned int i=0;i<qdctdc->n_qdc_hits;i++)
      qdc_calib_hits[i]=((qdctdc->qdc_values[i]&0xFFF)-qdc_ped[i])/qdc_gain[i];

    beamadc.setADCs06( qdc_calib_hits[bMuonVF], qdc_calib_hits[bMuonVB],
                    qdc_calib_hits[bMuonV1],qdc_calib_hits[bMuonV2],qdc_calib_hits[bMuonV3], 
                    qdc_calib_hits[bMuonV4],qdc_calib_hits[bMuonV5],qdc_calib_hits[bMuonV6], 
                    qdc_calib_hits[bMuonV7], qdc_calib_hits[bMuonV8],
                    qdc_calib_hits[bCerenkov1], qdc_calib_hits[bCerenkov2], qdc_calib_hits[bCerenkov3],
                    qdc_calib_hits[bScint1], qdc_calib_hits[bScint2], 
            qdc_calib_hits[bScint3], qdc_calib_hits[bScint4],
                    qdc_calib_hits[bTOF1S], qdc_calib_hits[bTOF1J],
                    qdc_calib_hits[bTOF2S], qdc_calib_hits[bTOF2J],
            qdc_calib_hits[bSCI_521],qdc_calib_hits[bSCI_528],
            qdc_calib_hits[bVH1],qdc_calib_hits[bVH2],
            qdc_calib_hits[bVH3],qdc_calib_hits[bVH4]);

       }
 }
  
}