sbnd::crt::CRTSlimmer Class Reference

Inheritance diagram for sbnd::crt::CRTSlimmer:

Public Member Functions
	CRTSlimmer (fhicl::ParameterSet const &p)
	CRTSlimmer (CRTSlimmer const &)=delete
	CRTSlimmer (CRTSlimmer &&)=delete
CRTSlimmer &	operator= (CRTSlimmer const &)=delete
CRTSlimmer &	operator= (CRTSlimmer &&)=delete
void	produce (art::Event &e) override

Private Attributes
std::string	_feb_data_producer
uint16_t	_adc_threshold

Detailed Description

Definition at line 50 of file CRTSlimmer_module.cc.

Constructor & Destructor Documentation

sbnd::crt::CRTSlimmer::CRTSlimmer ( fhicl::ParameterSet const &  p )

explicit

Definition at line 73 of file CRTSlimmer_module.cc.

  : EDProducer{p}
  , _feb_data_producer(p.get<std::string>("FEBDataProducer"))
  , _adc_threshold(p.get<uint16_t>("ADCThreshold"))
{
  produces<std::vector<sbnd::crt::CRTData>>();
  produces<art::Assns<sbnd::crt::CRTData, sbnd::crt::FEBData>>();
  produces<art::Assns<sbnd::crt::CRTData, sim::AuxDetIDE>>();

  consumes<std::vector<sbnd::crt::FEBData>>(_feb_data_producer);
  consumes<art::Assns<sbnd::crt::FEBData, sim::AuxDetIDE, sbnd::crt::FEBTruthInfo>>(_feb_data_producer);
}

sbnd::crt::CRTSlimmer::CRTSlimmer ( CRTSlimmer const &  )

delete

sbnd::crt::CRTSlimmer::CRTSlimmer ( CRTSlimmer &&  )

delete

Member Function Documentation

CRTSlimmer& sbnd::crt::CRTSlimmer::operator= ( CRTSlimmer const &  )

delete

CRTSlimmer& sbnd::crt::CRTSlimmer::operator= ( CRTSlimmer &&  )

delete

void sbnd::crt::CRTSlimmer::produce ( art::Event &  e )

override

Definition at line 86 of file CRTSlimmer_module.cc.

{
  // Implementation of required member function here.
  std::unique_ptr<std::vector<sbnd::crt::CRTData>> crt_data_v(new std::vector<sbnd::crt::CRTData>);
  std::unique_ptr<art::Assns<sbnd::crt::CRTData, sbnd::crt::FEBData>> crtdata_to_febdata_assns
      (new art::Assns<sbnd::crt::CRTData, sbnd::crt::FEBData>);
  std::unique_ptr<art::Assns<sbnd::crt::CRTData, sim::AuxDetIDE>> crtdata_to_ide_assns
      (new art::Assns<sbnd::crt::CRTData, sim::AuxDetIDE>);

  art::Handle<std::vector<sbnd::crt::FEBData>> feb_data_h;
  e.getByLabel(_feb_data_producer, feb_data_h);

  // make sure hits look good
  if (!feb_data_h.isValid()) {
    throw art::Exception(art::errors::Configuration) << "could not locate FEBData." << std::endl;;
  }

  std::vector<art::Ptr<sbnd::crt::FEBData>> feb_data_v;
  art::fill_ptr_vector(feb_data_v, feb_data_h);

  art::FindManyP<sim::AuxDetIDE, sbnd::crt::FEBTruthInfo> febdata_to_ides (feb_data_h, e, _feb_data_producer);

  art::PtrMaker<sbnd::crt::CRTData> makeDataPtr(e);

  for (size_t feb_i = 0; feb_i < feb_data_v.size(); feb_i++) {

    auto const feb_data = feb_data_v[feb_i];
    mf::LogDebug("CRTSlimmer") << "FEB " << feb_i << " with mac " << feb_data->Mac5() << std::endl;

    auto ides = febdata_to_ides.at(feb_data.key());
    mf::LogDebug("CRTSlimmer") << "We have " << ides.size() << " IDEs." << std::endl;

    // Construct a map to go from SiPM ID to the index of the IDE
    // FEBTruthInfo stores, for each AuxDetIDE, a vector containing
    // the SiPM IDs contributing to that energy deposit
    std::map<int, std::vector<int>> sipm_to_ideids;
    for (size_t j = 0; j < 32; j++) sipm_to_ideids[j] = std::vector<int>();
    for (size_t ide_i = 0; ide_i < febdata_to_ides.data(feb_i).size(); ide_i++) {
      const sbnd::crt::FEBTruthInfo *fti = febdata_to_ides.data(feb_i)[ide_i];
      sipm_to_ideids[fti->GetChannel()].push_back(ide_i);
      mf::LogDebug("CRTSlimmer") << "ide_i " << ide_i << " ene " << ides[ide_i]->energyDeposited << " fti->GetChannel() " << fti->GetChannel() << std::endl;
    }

    auto adcs = feb_data->ADC();

    for (size_t i = 0; i < adcs.size(); i+=2) {

      // Both sipms need to be above threshold to save them
      if (not (adcs[i] >= _adc_threshold and adcs[i+1] >= _adc_threshold)) {
        continue;
      }

      for (size_t sipm = 0; sipm < 2; sipm++) {
        sbnd::crt::CRTData crt_data = sbnd::crt::CRTData(feb_data->Mac5() * 32 + i + sipm,
                                                         feb_data->Ts0(),
                                                         feb_data->Ts1(),
                                                         adcs[i+sipm]);

        mf::LogDebug("CRTSlimmer") << "Adding SiPM " << i + sipm << " with mac " << feb_data->Mac5()
                                   << " mapped to channel " << crt_data.Channel() << std::endl;

        crt_data_v->emplace_back(std::move(crt_data));

        art::Ptr<sbnd::crt::CRTData> crt_data_p = makeDataPtr(crt_data_v->size() - 1);

        // Create the association between CRTData and FEBData
        crtdata_to_febdata_assns->addSingle(crt_data_p, feb_data);

        // Create the association between CRTData and AuxDetIDEs
        auto & ide_ids = sipm_to_ideids[i];
        for (auto ide_id: ide_ids) {
          auto ide = ides[ide_id];
          mf::LogDebug("CRTSlimmer") << "Adding IDE with ID " << ide_id << std::endl;
          crtdata_to_ide_assns->addSingle(crt_data_p, ide);
        }
      }
    }
  }

  mf::LogInfo("CRTSlimmer") << "Creating " << crt_data_v->size()
                            << " CRTData products." << std::endl;

  e.put(std::move(crt_data_v));
  e.put(std::move(crtdata_to_febdata_assns));
  e.put(std::move(crtdata_to_ide_assns));
}

Member Data Documentation

uint16_t sbnd::crt::CRTSlimmer::_adc_threshold

private

Definition at line 68 of file CRTSlimmer_module.cc.

std::string sbnd::crt::CRTSlimmer::_feb_data_producer

private

Definition at line 67 of file CRTSlimmer_module.cc.

---

The documentation for this class was generated from the following file:

• CRTSlimmer_module.cc