SimFastTiming/FastTimingCommon/test/MTDDigiDump.cc - master - githubmirror/cmssw

SimFastTiming/FastTimingCommon/test/MTDDigiDump.cc
131 строка4 KB
Fabio Cossutti
Update debug printouts for MTD in several packages, fix simulation LogDebug
21 фев 2023, 20:05
21 фев 2023, 20:0544d647a
#include <memory>
#include <iostream>
#include <vector>

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/FTLDigi/interface/FTLDigiCollections.h"

class MTDDigiDump : public edm::one::EDAnalyzer<edm::one::SharedResources> {
public:
  explicit MTDDigiDump(const edm::ParameterSet&);
  ~MTDDigiDump() override;

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

private:
  void beginJob() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void endJob() override;

  // ----------member data ---------------------------

  edm::EDGetTokenT<BTLDigiCollection> tok_BTL_digi;
  edm::EDGetTokenT<ETLDigiCollection> tok_ETL_digi;
};

MTDDigiDump::MTDDigiDump(const edm::ParameterSet& iConfig)

{
  tok_BTL_digi = consumes<BTLDigiCollection>(edm::InputTag("mix", "FTLBarrel"));
  tok_ETL_digi = consumes<ETLDigiCollection>(edm::InputTag("mix", "FTLEndcap"));
}

MTDDigiDump::~MTDDigiDump() {}

//
// member functions
//

// ------------ method called for each event ------------
void MTDDigiDump::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace std;

  edm::Handle<BTLDigiCollection> h_BTL_digi;
  iEvent.getByToken(tok_BTL_digi, h_BTL_digi);

  edm::Handle<ETLDigiCollection> h_ETL_digi;
  iEvent.getByToken(tok_ETL_digi, h_ETL_digi);

  // --- BTL DIGIs:

  if (!h_BTL_digi->empty()) {
    std::cout << " ----------------------------------------" << std::endl;
    std::cout << " BTL DIGI collection: \n" << std::endl;

    for (const auto& dataFrame : *h_BTL_digi) {
      // --- detector element ID:
      std::cout << "\n BTL DIGI:  row = " << dataFrame.row() << " col = " << dataFrame.column()
                << " BTLDetId = " << dataFrame.id();

      // --- loop over the dataFrame samples
      for (int isample = 0; isample < dataFrame.size(); ++isample) {
        const auto& sample = dataFrame.sample(isample);

        std::cout << "       sample " << isample << ":";
        if (sample.data() == 0 && sample.toa() == 0) {
          std::cout << std::endl;
          continue;
        }
        std::cout << "  amplitude = " << sample.data() << "  time1 = " << sample.toa() << "  time2 = " << sample.toa2()
                  << " r/c = " << (uint32_t)sample.row() << " / " << (uint32_t)sample.column()
                  << " th = " << sample.threshold() << " mode = " << sample.mode() << std::endl;

      }  // isaple loop

    }  // digi loop

  }  // if ( h_BTL_digi->size() > 0 )

  // --- ETL DIGIs:

  if (!h_ETL_digi->empty()) {
    std::cout << "\n ----------------------------------------" << std::endl;
    std::cout << " ETL DIGI collection: \n" << std::endl;

    for (const auto& dataFrame : *h_ETL_digi) {
      // --- detector element ID:
      std::cout << "\n ETL DIGI:  row = " << dataFrame.row() << " col = " << dataFrame.column()
                << " ETLDetId = " << dataFrame.id();

      // --- loop over the dataFrame samples
      for (int isample = 0; isample < dataFrame.size(); ++isample) {
        const auto& sample = dataFrame.sample(isample);

        std::cout << "       sample " << isample << ":";
        if (sample.data() == 0 && sample.toa() == 0) {
          std::cout << std::endl;
          continue;
        }
        std::cout << "  amplitude = " << sample.data() << "  time = " << sample.toa() << " r/c = " << sample.row()
                  << " / " << sample.column() << " th = " << sample.threshold() << " mode = " << sample.mode()
                  << std::endl;

      }  // isample loop

    }  // digi loop

  }  // if ( h_ETL_digi->size() > 0 )
}

// ------------ method called once each job just before starting event loop  ------------
void MTDDigiDump::beginJob() {}

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

// ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
void MTDDigiDump::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

//define this as a plug-in
DEFINE_FWK_MODULE(MTDDigiDump);