SumRunsCalib_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// \file    SumRunsCalib_module.cc
// \brief   Sum up AttenProfilesMap objects
// \author  Christopher Backhouse - bckhouse@caltech.edu
////////////////////////////////////////////////////////////////////////

#include <sstream>

// Framework includes
#include "art/Framework/Core/EDFilter.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Handle.h"
#include "art/Framework/Principal/Run.h"
#include "art/Framework/Principal/SubRun.h"

// NOvAsoft includes
#include "CalibrationDataProducts/AttenProfiles.h"

namespace calib
{
  /// Sum up AttenProfilesMap objects
  class SumRunsCalib: public art::EDFilter
  {
  public:
    explicit SumRunsCalib(const fhicl::ParameterSet& pset);
    ~SumRunsCalib();

    void reconfigure(const fhicl::ParameterSet& pset);
    virtual bool filter(art::Event& evt) {return false;}
    virtual bool endRun(art::Run& run);

  protected:
    std::list<std::pair<art::InputTag, int> > fInputTags;
    int                                       fMinBlock;
    int                                       fMaxBlock;
    unsigned int                              fLastRun;

    std::map<int, std::unique_ptr<caldp::AttenProfilesMap>> fChannelMapProf;
  };

  //......................................................................
  SumRunsCalib::SumRunsCalib(const fhicl::ParameterSet& pset)
  {
    reconfigure(pset);

    for(auto tag : fInputTags){
      produces<caldp::AttenProfilesMap, art::InRun>(tag.first.instance());
    }
  }

  //......................................................................
  SumRunsCalib::~SumRunsCalib()
  {
  }

  //......................................................................
  void SumRunsCalib::reconfigure(const fhicl::ParameterSet& pset)
  {
    fMinBlock          = pset.get<int         >("MinBlock");
    fMaxBlock          = pset.get<int         >("MaxBlock");
    fLastRun           = pset.get<unsigned int>("LastRun" );

    assert(fMinBlock >= 0);
    assert(fMaxBlock < 0 || fMaxBlock >= fMinBlock);
    
    // I don't think we can determine the detector in time
    if(fMaxBlock < 0) fMaxBlock = 27;
    
    // make the input tags
    fInputTags.clear();
    auto histsMapLabel = pset.get<std::string>("AttenHistsMapLabel");
    
    for(int block = fMinBlock; block < fMaxBlock + 1; ++block){
      std::ostringstream blockLabel;
      blockLabel << "block" << block;
      fInputTags.push_back(std::make_pair(art::InputTag(histsMapLabel, blockLabel.str()),
                                          block));
    }

  }

  //......................................................................
  bool SumRunsCalib::endRun(art::Run& run)
  {
    for(auto tag : fInputTags){
      art::Handle<caldp::AttenProfilesMap> profsmap;
      
      LOG_VERBATIM("SumRunsCalib")
      << run.run() << " " << tag.first.label() << " " << tag.first.instance();
      
      run.getByLabel(tag.first, profsmap);

      if(profsmap.failedToGet()) continue;

      if(!fChannelMapProf[tag.second]) fChannelMapProf[tag.second] = std::unique_ptr<caldp::AttenProfilesMap>(new caldp::AttenProfilesMap(profsmap->MinW(), profsmap->MaxW()));

      *fChannelMapProf[tag.second] += *profsmap;

      // No one else will ever want this. Save the memory
      run.removeCachedProduct(profsmap);
    }

    if(run.run() == fLastRun){
      for(auto& it: fChannelMapProf){
        std::ostringstream blockLabel;
        blockLabel << "block" << it.first;

        run.put(std::move(it.second), blockLabel.str());
      }

      fChannelMapProf.clear();

      return true;
    }

    return false;
  }

  DEFINE_ART_MODULE(SumRunsCalib)

} // namespace