A module to check the results from the Monte Carlo generator. More...

Inheritance diagram for rwgt::ReweightAna:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	ReweightAna (fhicl::ParameterSet const &pset)

virtual	~ReweightAna ()

void	analyze (art::Event const &evt)

void	beginSubRun (art::SubRun const &sr)

void	beginJob ()

void	endJob ()

void	endSubRun (art::SubRun const &sr)

void	reconfigure (const fhicl::ParameterSet &p)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Private Member Functions
void	LoadMCInfo (art::Event const &evt)

Private Attributes
TH1F *	fEnergyNeutrino
	Total number of events. More...

TH1F *	fNeventsSubrun
	Total number of events per subrun. More...

TH1F *	fWgtQE [3]

TH1F *	fWgtRES [3]

TH1F *	fWgtDIS [3]

rwgt::NuReweight *	fGrwgt [3]
	X-sec weight calculator. More...

std::string	fMCTruthModuleLabel
	label for module producing mc truth information More...

std::string	fPotLabel
	Module that produced the POTSum object. More...

Detailed Description

A module to check the results from the Monte Carlo generator.

Definition at line 50 of file ReweightAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

rwgt::ReweightAna::ReweightAna ( fhicl::ParameterSet const & pset )

explicit

Definition at line 85 of file ReweightAna_module.cc.

References reconfigure().

     : EDAnalyzer(p)
   {
     this->reconfigure(p);
   }

rwgt::ReweightAna::~ReweightAna ( )

virtual

Definition at line 93 of file ReweightAna_module.cc.

93 { }

Member Function Documentation

void rwgt::ReweightAna::analyze ( art::Event const & evt )

virtual

Implements art::EDAnalyzer.

Definition at line 138 of file ReweightAna_module.cc.

References rwgt::NuReweight::CalcWeight(), simb::MCNeutrino::CCNC(), simb::MCParticle::E(), fEnergyNeutrino, fGrwgt, fMCTruthModuleLabel, fWgtDIS, fWgtQE, fWgtRES, art::ProductRetriever::getByLabel(), simb::MCTruth::GetNeutrino(), MF_LOG_DEBUG, simb::MCNeutrino::Mode(), and simb::MCNeutrino::Nu().

   { 
 
     // // Pull the MC generator information out of the event
     mf::LogVerbatim("ReweightAna") << "Start analyze" ;
     art::Handle< std::vector<simb::MCTruth> > mclist;
     evt.getByLabel(fMCTruthModuleLabel, mclist);
     if (mclist->empty()) {
       mf::LogWarning("ReweightAna") << "Error retrieving MCTruth list" ;
       return;
     }
     
     art::Handle< std::vector<simb::GTruth> > gtlist;
     evt.getByLabel(fMCTruthModuleLabel, gtlist);
     if (gtlist->empty()) {
       mf::LogWarning("ReweightAna") << "Error retrieving GTruth list" ;
       return;
     }
     
    MF_LOG_DEBUG("ReweightAna")<<"MC List sizes:" << mclist->size() << " " << gtlist->size() << "\n";
     
     // // Loop over neutrino interactions
     for(size_t i_intx = 0; i_intx < mclist->size(); ++i_intx){
       MF_LOG_DEBUG("ReweightAna") << "start loop";
       
       //   // Link to the MCNeutrino class.
       //   // The class contains information not only about
       //   // the incoming neutrino, but about the products of the decay
       simb::MCTruth    const& truth       = mclist->at(i_intx);
       simb::GTruth     const& gtruth      = gtlist->at(i_intx);
       simb::MCNeutrino const& mc_neutrino = truth.GetNeutrino();
 
       fEnergyNeutrino->Fill(mc_neutrino.Nu().E());
       for(int i = 0; i < 3; i++) {
         double wgt = fGrwgt[i]->CalcWeight(truth, gtruth);
         //double wgt = 1.;
         if(mc_neutrino.Mode()==0 && mc_neutrino.CCNC()==0) {
           fWgtQE[i]->Fill(wgt);
         }
         else if(mc_neutrino.Mode()==1 && mc_neutrino.CCNC()==0) {
           fWgtRES[i]->Fill(wgt);
         }
         else if(mc_neutrino.Mode()==2 && mc_neutrino.CCNC()==0) {
           fWgtDIS[i]->Fill(wgt);
         }
       }
 
       MF_LOG_DEBUG("ReweightAna") << "end loop" ;
     }//end loop over interactions
 
  
     return;
   }

void rwgt::ReweightAna::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 96 of file ReweightAna_module.cc.

References rwgt::GENIEReweight::AddReweightValue(), rwgt::GENIEReweight::Configure(), fEnergyNeutrino, fGrwgt, fNeventsSubrun, rwgt::fReweightMaCCQE, rwgt::fReweightMaCCRES, rwgt::fReweightMaNCRES, rwgt::fReweightRvnCC1pi, rwgt::fReweightRvnCC2pi, rwgt::fReweightRvnNC1pi, rwgt::fReweightRvnNC2pi, rwgt::fReweightRvpCC1pi, rwgt::fReweightRvpCC2pi, rwgt::fReweightRvpNC1pi, rwgt::fReweightRvpNC2pi, fWgtDIS, fWgtQE, and fWgtRES.

   {
     art::ServiceHandle<art::TFileService> tfs;
    
     mf::LogVerbatim("ReweightAna") << "make histograms" ;
 
     fEnergyNeutrino   = tfs->make<TH1F>("fEnergyneutrino", "Total number of events", 50, 0., 25); 
     fNeventsSubrun    = tfs->make<TH1F>("fNeventsSubrun", "Total number of events", 1, 0., 1.); 
 
     char name[300];
     for(int i = 0; i < 3; i++) {
       sprintf(name, "fWgtQE_%dsigma", i+1);
       fWgtQE[i] = tfs->make<TH1F>(name, "Evt Wgts", 100, 0., 2.0); 
       sprintf(name, "fWgtRES_%dsigma", i+1);
       fWgtRES[i] = tfs->make<TH1F>(name, "Evt Wgts", 100, 0., 2.0); 
       sprintf(name, "fWgtDIS_%dsigma", i+1);
       fWgtDIS[i] = tfs->make<TH1F>(name, "Evt Wgts", 100, 0., 2.0); 
 
       double sigma = (double)(i+1);
       fGrwgt[i] = new rwgt::NuReweight();
       fGrwgt[i]->AddReweightValue(rwgt::fReweightMaCCQE, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightMaCCRES, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightMaNCRES, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvpCC1pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvnCC1pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvpCC2pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvnCC2pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvpNC1pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvnNC1pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvpNC2pi, sigma);
       fGrwgt[i]->AddReweightValue(rwgt::fReweightRvnNC2pi, sigma);
       
       fGrwgt[i]->Configure();
     }
   }

void rwgt::ReweightAna::beginSubRun ( art::SubRun const & sr )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 133 of file ReweightAna_module.cc.

                                                    {
     
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Analyzer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Analyzer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Analyzer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doEvent	(	EventPrincipal &	ep,
		ModuleContext const &	mc,
		std::atomic< std::size_t > &	counts_run,
		std::atomic< std::size_t > &	counts_passed,
		std::atomic< std::size_t > &	counts_failed
	)

inherited

Definition at line 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     auto const e = std::as_const(ep).makeEvent(mc);
     if (wantEvent(mc.scheduleID(), e)) {
       ++counts_run;
       ProcessingFrame const frame{mc.scheduleID()};
       analyzeWithFrame(e, frame);
       ++counts_passed;
     }
     return true;
   }

void art::detail::Analyzer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void rwgt::ReweightAna::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 216 of file ReweightAna_module.cc.

References DEFINE_ART_MODULE.

   {
 
   }

void rwgt::ReweightAna::endSubRun ( art::SubRun const & sr )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 208 of file ReweightAna_module.cc.

References rwgt::cntEvent, fNeventsSubrun, and art::SubRun::subRun().

   {
     fNeventsSubrun->Fill(sr.subRun(), cntEvent);
     cntEvent = 0;
 
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

Handle< TriggerResults > art::Observer::getTriggerResults ( Event const & e ) const

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (selectors_) {
       return selectors_->getOneTriggerResults(e);
     }
 
     // The following applies for cases where no SelectEvents entries
     // exist.
     Handle<TriggerResults> h;
     if (e.get(empty_process_name, h)) {
       return h;
     }
     return Handle<TriggerResults>{};
   }

void rwgt::ReweightAna::LoadMCInfo ( art::Event const & evt )

private

Definition at line 194 of file ReweightAna_module.cc.

   {
 
   }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

   {
     return process_name_;
   }

void rwgt::ReweightAna::reconfigure ( const fhicl::ParameterSet & p )

Definition at line 202 of file ReweightAna_module.cc.

References fMCTruthModuleLabel, and fhicl::ParameterSet::get().

Referenced by ReweightAna().

   {    
     fMCTruthModuleLabel = p.get< std::string>("MCTruthModuleLabel");
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

     {
       return wantAllEvents_;
     }

bool art::Observer::wantEvent	(	ScheduleID	id,
		Event const &	e
	)		const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

Referenced by art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (wantAllEvents_) {
       return true;
     }
     bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
     bool const reject_event =
       rejectors_ ? rejectors_->matchEvent(id, e) : false;
     return select_event and not reject_event;
   }

Member Data Documentation

TH1F* rwgt::ReweightAna::fEnergyNeutrino

private

Total number of events.

Definition at line 67 of file ReweightAna_module.cc.

Referenced by analyze(), and beginJob().

rwgt::NuReweight* rwgt::ReweightAna::fGrwgt[3]

private

X-sec weight calculator.

Definition at line 72 of file ReweightAna_module.cc.

Referenced by analyze(), and beginJob().

std::string rwgt::ReweightAna::fMCTruthModuleLabel

private

label for module producing mc truth information

Definition at line 74 of file ReweightAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1F* rwgt::ReweightAna::fNeventsSubrun

private

Total number of events per subrun.

Definition at line 68 of file ReweightAna_module.cc.

Referenced by beginJob(), and endSubRun().

std::string rwgt::ReweightAna::fPotLabel

private

Module that produced the POTSum object.

Definition at line 75 of file ReweightAna_module.cc.

TH1F* rwgt::ReweightAna::fWgtDIS[3]

private

Definition at line 71 of file ReweightAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* rwgt::ReweightAna::fWgtQE[3]

private

Definition at line 69 of file ReweightAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* rwgt::ReweightAna::fWgtRES[3]

private

Definition at line 70 of file ReweightAna_module.cc.

Referenced by analyze(), and beginJob().

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

nugen/NuReweight/art/ReweightAna_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation