CSCDigiToPattern.cc

Go to the documentation of this file.

#include "EventFilter/CSCRawToDigi/src/CSCDigiToPattern.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"

CSCDigiToPattern::CSCDigiToPattern(edm::ParameterSet const& conf) {

  // If your module takes parameters, here is where you would define
  // their names and types, and access them to initialize internal
  // variables. Example as follows:
  //
}

void CSCDigiToPattern::analyze(edm::Event const& e, edm::EventSetup const& iSetup) {

  // These declarations create handles to the types of records that you want
  // to retrieve from event "e".
  //
  edm::Handle<CSCCorrelatedLCTDigiCollection> correlatedlcts;
  e.getByLabel("muonCSCDigis","MuonCSCCorrelatedLCTDigi",correlatedlcts);  
  for (CSCCorrelatedLCTDigiCollection::DigiRangeIterator j=correlatedlcts->begin(); j!=correlatedlcts->end(); j++) { 
    CSCDetId id=(*j).first;
    std::cout<<id<<std::endl;
    std::vector<CSCCorrelatedLCTDigi>::const_iterator digiIt = (*j).second.first;
    std::vector<CSCCorrelatedLCTDigi>::const_iterator last = (*j).second.second;
    for( ; digiIt != last; ++digiIt) {
      uint16_t wire       = digiIt->getKeyWG();    // 7 bits
      uint16_t pattern    = digiIt->getPattern();  // 4 bits
      uint16_t quality    = digiIt->getQuality();  // 4 bits
      uint16_t valid      = digiIt->isValid();     // 1 bit
      uint16_t strip      = digiIt->getStrip();    // 8 bits
      uint16_t bend       = digiIt->getBend();     // 1 bit
      uint16_t syncErr    = digiIt->getSyncErr();  // 1 bit
      uint16_t bx         = digiIt->getBX();       // 1 bit
      uint16_t bx0        = digiIt->getBX0();      // 1 bit
      uint16_t cscId      = digiIt->getCSCID();    // 4 bits
      //                                             __
      //                                             32 bits in total
      long unsigned int mpc =
    ((cscId&0xF)<<28) | ((bx0&0x1)<<27) | ((bx&0x1)<<26) |
    ((syncErr&0x1)<<25) | ((bend&0x1)<<24) | ((strip&0xFF)<<16) |
    ((valid&0x1)<<15) | ((quality&0xF)<<11) | ((pattern&0xF)<<7) |
    (wire&0x7F);
      std::cout <<"MPC"<<digiIt->getTrknmb()<< " " << std::hex << mpc <<std::dec <<std::endl;
    }
  }
}