MakeSNeTProfileHistos Class Reference

Inheritance diagram for MakeSNeTProfileHistos:

Public Member Functions
	MakeSNeTProfileHistos (fhicl::ParameterSet const &p)
	MakeSNeTProfileHistos (MakeSNeTProfileHistos const &)=delete
	MakeSNeTProfileHistos (MakeSNeTProfileHistos &&)=delete
MakeSNeTProfileHistos &	operator= (MakeSNeTProfileHistos const &)=delete
MakeSNeTProfileHistos &	operator= (MakeSNeTProfileHistos &&)=delete
void	analyze (art::Event const &e) override
void	reconfigure (fhicl::ParameterSet const &p)
void	beginJob () override
void	endJob () override
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
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)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
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)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fTruthLabel
	label for process that made the MCTruth info More...
std::string	fHitLabel
	label for process that made the reco hits More...
TH1F *	fEnergySpec
	1D event energy spectrum (used only for backgrounds) More...
TH1F *	fTrigEff_numerator
	how many times a trigger was issued per true event energy More...
TH1F *	fTrigEff
	trigger efficiency per true event energy (fTrigEff_numerator/fEnergySpec) More...
TH2F *	fNHitsVsEnergy
	number of reconstructed hits per true event energy More...
TH2F *	fSumADCVsEnergy
	summed ADC for all reco hits per true event energy More...
TH1F *	fTPSplitProb
	probability of splitting a trigger primitive into two (or more) TPs per true event energy More...
TH1F *	fNMCTruths
	number of MCTruths per art::event More...
TH1F *	fCCNC
	CCNC for each MCTruth. More...
TH1F *	fOrigin
	mctruth->Origin() More...
TH1F *	fNuPDG
	mc_neutrino->Nu()->PdgCode() More...

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
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 ()

Detailed Description

Definition at line 43 of file MakeSNeTProfileHistos_module.cc.

Constructor & Destructor Documentation

MakeSNeTProfileHistos::MakeSNeTProfileHistos ( fhicl::ParameterSet const &  p )

explicit

Definition at line 95 of file MakeSNeTProfileHistos_module.cc.

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

MakeSNeTProfileHistos::MakeSNeTProfileHistos ( MakeSNeTProfileHistos const &  )

delete

MakeSNeTProfileHistos::MakeSNeTProfileHistos ( MakeSNeTProfileHistos &&  )

delete

Member Function Documentation

void MakeSNeTProfileHistos::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 212 of file MakeSNeTProfileHistos_module.cc.

{

  // TODO:
  // 
  // 1. implement dividing the TrigEff histos in the endJob function
  // 2. implement anything that has anything to do with the trigger data products
  //       (including the TrigEff histos, TP splitting histos, etc...???)
  // 3. figure out a good way to generalize this for multiple MCTruths...




  // ===== Notes and Assumptions ===== //
  //
  // This module assumes the following:
  //
  // 1. There is only one neutrino interaction per art::event. This is
  //    important because it assumes that all reco hits go with that one
  //    neutrino.





  //
  // Lift out the desired data products:
  //

  // Get the hits out of the event
  auto hits_list = e.getHandle< std::vector< recob::Hit > >(fHitLabel);

  // Get the MCTruths out of the event
  auto truths_list = e.getHandle< std::vector< simb::MCTruth > >(fTruthLabel);
  
  fNMCTruths->Fill(truths_list->size());

  // Assert assumption that there is one and only one MCTruth in the event
  // (note:   this needs to be fixed later so that it can handle multiple MCTruths.
  //          This only works for single generated SNe neutrinos for now...)
  if(truths_list->size() == 1) {

    // calculate Summed ADC
    float SummedADC = 0.0; // ...I know... this number is a float?! Who would have guessed that...
    for(unsigned int i = 0; i < hits_list->size();i++){
      recob::Hit const& hit = hits_list->at(i);  
      SummedADC += hit.SummedADC();
    }

    
    simb::MCTruth const& truth = truths_list->at(0);

    fOrigin->Fill(truth.Origin());

    if (truth.NeutrinoSet()){
      simb::MCNeutrino const& mc_nu = truth.GetNeutrino();
      double NuE = 1000.0*mc_nu.Nu().E(); // convert to MeV

      fEnergySpec->Fill(NuE);
      fNHitsVsEnergy->Fill(NuE,hits_list->size());
      fSumADCVsEnergy->Fill(NuE,SummedADC);

      fCCNC->Fill(mc_nu.CCNC());
      fNuPDG->Fill(mc_nu.Nu().PdgCode());

    }    

  }
  else {
    std::cerr << "\n\n\nWARNING: Number of MCTruths != 1\n\n\n\n";
  }



}

void MakeSNeTProfileHistos::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 114 of file MakeSNeTProfileHistos_module.cc.

{

  art::ServiceHandle<art::TFileService> tfs;

  // binning for the energy axes:
  int NbinsE = 200;
  float minE = 0.0;   // in MeV
  float maxE = 100.0; // in MeV

  // binning for the NHits axes:
  int NbinsNH = 5000;
  float minNH = 0.0;
  float maxNH = 5000.0; // a reasonable number determined from running over single nue events

  // binning for the ADC axes:
  int NbinsADC = 500;
  float minADC = 0.0;
  float maxADC = 500.0e3; // a reasonable number determined from running over single nue events



  //
  // Book required histograms for the Time-Profile simulator:
  //

  // 2D neutrino energy spectrum as a function of time:
  //   (Note: generated separately by Kate Schol.)

  // 1D event energy spectrum: (used only for backgrounds)
  fEnergySpec = tfs->make<TH1F>("EnergySpec",
                                "True energy spectra;true energy [MeV];",
                                NbinsE, minE, maxE);

  // Overall efficency for triggering as a function of true event energy:
  fTrigEff_numerator = tfs->make<TH1F>("TrigEff_numerator",
                                       "Total number of triggers issued;true energy [MeV];",
                                       NbinsE, minE, maxE);
  fTrigEff           = tfs->make<TH1F>("TrigEff",
                                       "Trigger efficiency;true energy [MeV];",
                                       NbinsE, minE, maxE);

  // Number of reco hits vs. true event energy:
  fNHitsVsEnergy = tfs->make<TH2F>("NHitsVsEnergy",
                                   "Number of reco hits per true energy;true energy [MeV];NHits",
                                   NbinsE, minE, maxE,
                                   NbinsNH, minNH, maxNH);

  // Summed ADC vs. true event energy:
  fSumADCVsEnergy = tfs->make<TH2F>("SumADCVsEnergy",
                                    "Summed ADC over all reco hits per true energy;true energy [MeV];Summed ADC",
                                    NbinsE, minE, maxE,
                                    NbinsADC, minADC, maxADC);

  // Probability of splitting a TP into 2 TPs as a function of true event energy:
  fTPSplitProb = tfs->make<TH1F>("TPSplitProb",
                                 "Probability of spliting a TP into 2 or more per true energy;true energy [MeV];",
                                 NbinsE, minE, maxE);



  //
  // Book a few extra check histos:
  //

  // Number of MCTruths per art::event:
  fNMCTruths = tfs->make<TH1F>("NMCTruths",
                               "Number of MCTruths per art::event;# of MCTruths;",
                               101,-0.5,100.5);

  // CCNC:
  fCCNC = tfs->make<TH1F>("CCNC",
                          "CCNC;CCNC;",
                          11,-5.5,5.5);

  // Origin:
  fOrigin = tfs->make<TH1F>("Origin",
                            "Neutrino origin;origin type;",
                            16,-5.5,10.5);

  // Neutrino pdg code:
  fNuPDG = tfs->make<TH1F>("NuPDG",
                           "Neutrino PDG code;pdg code;",
                           41,-20.5,20.5);

}

void MakeSNeTProfileHistos::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 204 of file MakeSNeTProfileHistos_module.cc.

{

}

MakeSNeTProfileHistos& MakeSNeTProfileHistos::operator= ( MakeSNeTProfileHistos const &  )

delete

MakeSNeTProfileHistos& MakeSNeTProfileHistos::operator= ( MakeSNeTProfileHistos &&  )

delete

void MakeSNeTProfileHistos::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 105 of file MakeSNeTProfileHistos_module.cc.

{
  fTruthLabel = p.get<std::string>("TruthLabel");
  fHitLabel   = p.get<std::string>("HitLabel");
}

Member Data Documentation

TH1F* MakeSNeTProfileHistos::fCCNC

private

CCNC for each MCTruth.

Definition at line 86 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fEnergySpec

private

1D event energy spectrum (used only for backgrounds)

Definition at line 75 of file MakeSNeTProfileHistos_module.cc.

std::string MakeSNeTProfileHistos::fHitLabel

private

label for process that made the reco hits

Definition at line 70 of file MakeSNeTProfileHistos_module.cc.

TH2F* MakeSNeTProfileHistos::fNHitsVsEnergy

private

number of reconstructed hits per true event energy

Definition at line 78 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fNMCTruths

private

number of MCTruths per art::event

Definition at line 85 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fNuPDG

private

mc_neutrino->Nu()->PdgCode()

Definition at line 88 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fOrigin

private

mctruth->Origin()

Definition at line 87 of file MakeSNeTProfileHistos_module.cc.

TH2F* MakeSNeTProfileHistos::fSumADCVsEnergy

private

summed ADC for all reco hits per true event energy

Definition at line 79 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fTPSplitProb

private

probability of splitting a trigger primitive into two (or more) TPs per true event energy

Definition at line 80 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fTrigEff

private

trigger efficiency per true event energy (fTrigEff_numerator/fEnergySpec)

Definition at line 77 of file MakeSNeTProfileHistos_module.cc.

TH1F* MakeSNeTProfileHistos::fTrigEff_numerator

private

how many times a trigger was issued per true event energy

Definition at line 76 of file MakeSNeTProfileHistos_module.cc.

std::string MakeSNeTProfileHistos::fTruthLabel

private

label for process that made the MCTruth info

Definition at line 69 of file MakeSNeTProfileHistos_module.cc.

---

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

• duneana/duneana/SupernovaAna/MakeSNeTProfileHistos_module.cc