SiStripFEDEmulator.h

Go to the documentation of this file.

#ifndef DQM_SiStripMonitorHardware_SiStripFEDEmulator_H
#define DQM_SiStripMonitorHardware_SiStripFEDEmulator_H

// Created 2010-01-20 by A.-M. Magnan
// Class intended to mimic the data path in the FED firmware in software
// steps: pedestal subtraction, CM subtraction, clustering and zero-suppression.
#include <sstream>
#include <fstream>
#include <iostream>

#include "FWCore/Framework/interface/ESHandle.h"

#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/SiStripDigi/interface/SiStripRawDigi.h"
#include "DataFormats/SiStripDigi/interface/SiStripProcessedRawDigi.h"
#include "DataFormats/SiStripDigi/interface/SiStripDigi.h"

#include "DataFormats/Common/interface/Wrapper.h"

#include "CondFormats/SiStripObjects/interface/SiStripFedCabling.h"
#include "CondFormats/SiStripObjects/interface/SiStripPedestals.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"

#include "EventFilter/SiStripRawToDigi/interface/SiStripFEDBuffer.h"

//for the zero suppression algorithm(s)
#include "RecoLocalTracker/SiStripZeroSuppression/interface/SiStripRawProcessingFactory.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/SiStripRawProcessingAlgorithms.h"

namespace sistrip {

  class FEDEmulator {

  public:
    FEDEmulator();
    ~FEDEmulator();

    void initialise(const bool byModule);

    void initialiseModule(const uint32_t aDetId,
              const uint32_t aNPairs,
              const uint32_t aPair);

    void retrievePedestals(const edm::ESHandle<SiStripPedestals> & aHandle);
    void retrieveNoises(const edm::ESHandle<SiStripNoises> & aHandle);

    void subtractPedestals(const edm::DetSetVector<SiStripRawDigi>::const_iterator & inputChannel,
               std::vector<SiStripRawDigi> & pedsDetSetData,
               std::vector<SiStripProcessedRawDigi> & noiseDetSetData,
               std::vector<SiStripRawDigi> & pedSubtrDetSetData,
               std::vector<uint32_t> & medsDetSetData,
               const bool fillApvsForCM);

    void subtractCM(const std::vector<SiStripRawDigi> & pedSubtrDetSetData,
            std::vector<SiStripRawDigi> & cmSubtrDetSetData);

    void zeroSuppress(const std::vector<SiStripRawDigi> & cmSubtrDetSetData,
              edm::DetSet<SiStripDigi>    & zsDetSetData,
              const std::auto_ptr<SiStripRawProcessingAlgorithms> & algorithms);

    uint32_t fedIndex(const uint16_t aFedChannel);

    void fillPeds(const edm::DetSetVector<SiStripRawDigi>::const_iterator & peds);
    void fillNoises(const edm::DetSetVector<SiStripProcessedRawDigi>::const_iterator & noise);

    void fillMedians(const std::map<uint32_t,std::vector<uint32_t> >::const_iterator & meds);

    void print(std::ostream & aOs);
    void printPeds(std::ostream & aOs);
    void printNoises(std::ostream & aOs);
    void printMeds(std::ostream & aOs);

  private:

    static const char* messageLabel_;

    bool byModule_;

    uint32_t detId_;
    uint32_t nPairs_;
    uint32_t pair_;

    uint32_t minStrip_;
    uint32_t maxStrip_;

    std::vector<int> pedestals_;
    std::vector<float> noises_;
    std::vector<uint32_t> medians_;



    
  };//class FEDEmulator

}//namespace sistrip
#endif//DQM_SiStripMonitorHardware_SiStripFEDEmulator_H