SystematicUncertaintyScaleAna_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// \file    SystematicUncertaintyScaleAna_module.cc
// \brief
// \author  Brian Rebel brebel@fnal.gov
////////////////////////////////////////////////////////////////////////

// Framework includes
#include "art/Framework/Core/EDAnalyzer.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art_root_io/TFileService.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Handle.h"

#include "Calibrator/Calibrator.h"

#include "TH1.h"

namespace calib
{
  
  class SystematicUncertaintyScaleAna: public art::EDAnalyzer
  {
  public:
    explicit SystematicUncertaintyScaleAna(fhicl::ParameterSet const& pset);
    ~SystematicUncertaintyScaleAna();
    
    void analyze(art::Event const& evt);
    
    void reconfigure(fhicl::ParameterSet const& pset);
    
    void endJob();

  private:

    TH1F *fMuCXViewW; ///< W for the X view Muon Catcher
    TH1F *fXViewW;    ///< W for the X view
    TH1F *fYViewW;    ///< W for the Y view
    
  };
  
  //......................................................................
  SystematicUncertaintyScaleAna::SystematicUncertaintyScaleAna(const fhicl::ParameterSet& pset)
  : EDAnalyzer(pset)
  {
    this->reconfigure(pset);
    
    return;
  }
  
  //......................................................................
  SystematicUncertaintyScaleAna::~SystematicUncertaintyScaleAna()
  {
  }
  
  //......................................................................
  void SystematicUncertaintyScaleAna::reconfigure(fhicl::ParameterSet const& /*pset*/)
  {
  }
  
  //......................................................................
  void SystematicUncertaintyScaleAna::endJob()
  {
    art::ServiceHandle<art::TFileService> tfs;
    
    art::ServiceHandle<geo::Geometry> geo;
    
    float halfWidth = geo->DetHalfWidth();
    int   numBins   = 200;
    
    fMuCXViewW = tfs->make<TH1F>("MuCxViewW", ";w (cm);Systematic Scale Factor", numBins, -1.*halfWidth, halfWidth);
    fXViewW    = tfs->make<TH1F>("xViewW", ";w (cm);Systematic Scale Factor", numBins, -1.*halfWidth, halfWidth);
    fYViewW    = tfs->make<TH1F>("yViewW", ";w (cm);Systematic Scale Factor", numBins, -1.*halfWidth, halfWidth);
    
    // this module tests the calibrator as configured for the current job
    // one could create new parameter sets to do multiple different tests
    // of the calibrator by calling the Calibrator::reconfigure() method
    // for each new parameter set
    art::ServiceHandle<calib::Calibrator> cal;

    for(int b = 1; b < fXViewW->GetNbinsX(); ++b){
      const double w = fXViewW->GetBinCenter(b);
      if( geo->DetId() == novadaq::cnv::kNEARDET )
        fMuCXViewW->Fill(w, cal->SystematicUncertaintyScale(geo::kX, geo->FirstPlaneInMuonCatcher()+1, w));
      fXViewW   ->Fill(w, cal->SystematicUncertaintyScale(geo::kX, 0,   w)); // plane moot
      fYViewW   ->Fill(w, cal->SystematicUncertaintyScale(geo::kY, 0,   w)); // plane moot
    }

    
    return;
  }
  
  //......................................................................
  void SystematicUncertaintyScaleAna::analyze(art::Event const& /*evt*/)
  {
    return;
  }
  
  DEFINE_ART_MODULE(SystematicUncertaintyScaleAna)
  
} // end namespace