SupernovaAna Class Reference

Inheritance diagram for SupernovaAna:

Public Member Functions
	SupernovaAna (fhicl::ParameterSet const &p)
	SupernovaAna (SupernovaAna const &)=delete
	SupernovaAna (SupernovaAna &&)=delete
SupernovaAna &	operator= (SupernovaAna const &)=delete
SupernovaAna &	operator= (SupernovaAna &&)=delete
void	analyze (art::Event const &e) override
void	reconfigure (fhicl::ParameterSet const &p)
void	beginJob () 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
std::string	fHitLabel
TTree *	fNTuple
int	EventNum
int	EventNumCut
int	NHits
double	avgRMS
double	avgRMSloop [kMaxNumHits]
double	avgRMSloopsum
double	SummedADC
int	NMCTruths
int	NParticles [kMaxNumParticles]
double	MomentumYTruths [kMaxNumParticles]
double	MomentumZTruths [kMaxNumParticles]
double	NuEnergyTruths [kMaxNumParticles]
double	LeptonEnergyTruths [kMaxNumParticles]
double	ThetaTruths [kMaxNumParticles]
int	OriginTruths [kMaxNumParticles]
int	CCNCTruths [kMaxNumParticles]
int	ModeTruths [kMaxNumParticles]
int	NuPDGCodeTruths [kMaxNumParticles]
int	LeptonPDGCodeTruths [kMaxNumParticles]
double	xTruths [kMaxNumParticles]
double	yTruths [kMaxNumParticles]
double	zTruths [kMaxNumParticles]
double	tTruths [kMaxNumParticles]
double	EndxTruths [kMaxNumParticles]
double	EndyTruths [kMaxNumParticles]
double	EndzTruths [kMaxNumParticles]
double	EndtTruths [kMaxNumParticles]
double	DeltaXTruths [kMaxNumParticles]
double	DeltaYTruths [kMaxNumParticles]

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 42 of file SupernovaAna_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 108 of file SupernovaAna_module.cc.

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

SupernovaAna::SupernovaAna ( SupernovaAna const &  )

delete

SupernovaAna::SupernovaAna ( SupernovaAna &&  )

delete

Member Function Documentation

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

overridevirtual

Implements art::EDAnalyzer.

Definition at line 169 of file SupernovaAna_module.cc.

{

  
  //==============================HITS================================

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

  // Gets the number of hits per event and fills into a histogram
  NHits = hits_list->size();
 
  // Get summed ADC and average RMS of hit slope
  SummedADC = 0.0;
  for(unsigned int i = 0; i < hits_list->size();i++){
    recob::Hit const& hit = hits_list->at(i);  
    avgRMSloop[i] = hit.RMS();
    avgRMSloopsum += avgRMSloop[i];
    SummedADC += hit.SummedADC();
   }
  avgRMS = avgRMSloopsum/(hits_list->size());

  //============================MCTRUTHS==============================

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

  // Gets the number of MC truths per event and fills into a histogram
  NMCTruths = truths_list->size();
  
  // Loop through the truth arrays and sets the values of the arrays
  // to an arbitrary default (unfilled) value
  for(int i=0;i<kMaxNumParticles;i++){
    xTruths[i] = -1.0e9;
    yTruths[i] = -1.0e9;
    zTruths[i] = -1.0e9;
    tTruths[i] = -1.0e9;
    EndxTruths[i] = -1.0e9;
    EndyTruths[i] = -1.0e9;
    EndzTruths[i] = -1.0e9;
    EndtTruths[i] = -1.0e9;
    NParticles[i] = -999;
    MomentumYTruths[i] = -1.0e9;
    MomentumZTruths[i] = -1.0e9;
    NuEnergyTruths[i] = -999;
    LeptonEnergyTruths[i] = -999;
    ThetaTruths[i] = -999;
    OriginTruths[i] = -999;
    CCNCTruths[i] = -999;
    ModeTruths[i] = -999;
    NuPDGCodeTruths[i] = 0;
    LeptonPDGCodeTruths[i] = 0;
    DeltaXTruths[i] = -1.0e9;
    DeltaYTruths[i] = -1.0e9;

  }
  
  // Detect if the number of MC truths per event is greater than 
  // or equal to the maximum number of particles per event else 
  // we'll fall off the end of the array and cause a seg fault
  if(NMCTruths >= kMaxNumParticles) { 
    std::cerr << "ERROR: NMCTruths " << NMCTruths << 
    " >= kMaxNumParticles " << kMaxNumParticles << 
    " , will cause a segmentation fault" << std::endl;}
  // Detect if the number of hits per event is greater than 
  // or equal to the maximum number of particles per event else 
  // we'll fall off the end of the array and cause a seg fault
  if(NHits >= kMaxNumHits) { 
    std::cerr << "ERROR: NHits " << NHits << 
    " >= kMaxNumHits " << kMaxNumHits << 
    " , will cause a segmentation fault" << std::endl;}
  
  // Create nTuple data for any number of truths 
  for(unsigned int i = 0; i < truths_list->size();i++){
    simb::MCTruth const& truth = truths_list->at(i);

    NParticles[i] = truth.NParticles();

    // Get the ith particle's Origin truths
    OriginTruths[i] = truth.Origin();
           
    // Only get neutrino truths if they exist to avoid a segmentation fault
    if (truth.NeutrinoSet()){
      simb::MCNeutrino const& mc_neutrino = truth.GetNeutrino();
      // should change notation such that mc_neutrino = truth.GetNeutrino().Nu();     // and then remove .Nu() when called
      simb::MCParticle const& mc_lepton = mc_neutrino.Lepton();

      // Get data to be filled into nTuple       
      xTruths[i] = mc_lepton.Vx();
      yTruths[i] = mc_lepton.Vy();
      zTruths[i] = mc_lepton.Vz();
      tTruths[i] = mc_lepton.T();
      EndxTruths[i] = mc_lepton.EndX();
      EndyTruths[i] = mc_lepton.EndY();
      EndzTruths[i] = mc_lepton.EndZ();
      EndtTruths[i] = mc_lepton.EndT();
      MomentumYTruths[i] = mc_lepton.Py();
      MomentumZTruths[i] = mc_lepton.Pz();
      NuEnergyTruths[i] = mc_neutrino.Nu().E();
      LeptonEnergyTruths[i] = mc_lepton.E();
      ThetaTruths[i] = mc_neutrino.Theta();
      CCNCTruths[i] = mc_neutrino.CCNC();
      ModeTruths[i] = mc_neutrino.Mode();
      NuPDGCodeTruths[i] = mc_neutrino.Nu().PdgCode();
      LeptonPDGCodeTruths[i] = mc_lepton.PdgCode();
      DeltaXTruths[i] = EndxTruths[i] - xTruths[i];
      DeltaYTruths[i] = EndyTruths[i] - yTruths[i];
    }
  }
  
    art::EventNumber_t event = e.id().event();
    EventNum = event;
    if (NHits<2) {
      EventNumCut = event;
    }


  // Now we should have NMCTruths = truths_list->size() and 
  // PDGCodeTruths [0 - (NMCTruths - 1)] should contain PdgCodes, 
  // [ NMCTruths - 500 ] should be equal to -999, signifying that 
  // we did not fill this part
  
  // Fill nTuple
  fNTuple->Fill();
}

void SupernovaAna::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 127 of file SupernovaAna_module.cc.

{

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

  // Make nTuple
  fNTuple = tfs->make<TTree>("SupernovaAna","Supernonva analysis tree");

  // Add branches to nTuple
  fNTuple->Branch("EventNum",&EventNum);
  fNTuple->Branch("EventNumCut",&EventNumCut);
  fNTuple->Branch("NHits",&NHits);
  fNTuple->Branch("avgRMS",&avgRMS);
  fNTuple->Branch("SummedADC",&SummedADC);
  fNTuple->Branch("NMCTruths",&NMCTruths);
  fNTuple->Branch("NParticles",NParticles,"NParticles[NMCTruths]/I");
  fNTuple->Branch("MomentumYTruths",MomentumYTruths,"MomentumYTruths[NMCTruths]/D");
  fNTuple->Branch("MomentumZTruths",MomentumZTruths,"MomentumZTruths[NMCTruths]/D");
  fNTuple->Branch("NuEnergyTruths",NuEnergyTruths,"NuEnergyTruths[NMCTruths]/D");
  fNTuple->Branch("LeptonEnergyTruths",LeptonEnergyTruths,"LeptonEnergyTruths[NMCTruths]/D");
  fNTuple->Branch("ThetaTruths",ThetaTruths,"ThetaTruths[NMCTruths]/D");
  fNTuple->Branch("OriginTruths",OriginTruths,"OriginTruths[NMCTruths]/I");
  fNTuple->Branch("CCNCTruths",CCNCTruths,"CCNCTruths[NMCTruths]/I");
  fNTuple->Branch("ModeTruths",ModeTruths,"ModeTruths[NMCTruths]/I");
  fNTuple->Branch("NuPDGCodeTruths",NuPDGCodeTruths,"NuPDGCodeTruths[NMCTruths]/I");
  fNTuple->Branch("LeptonPDGCodeTruths",LeptonPDGCodeTruths,"LeptonPDGCodeTruths[NMCTruths]/I");
  fNTuple->Branch("xTruths",xTruths,"xTruths[NMCTruths]/D");
  fNTuple->Branch("yTruths",yTruths,"yTruths[NMCTruths]/D");
  fNTuple->Branch("zTruths",zTruths,"zTruths[NMCTruths]/D");
  fNTuple->Branch("tTruths",tTruths,"tTruths[NMCTruths]/D");
  fNTuple->Branch("EndxTruths",EndxTruths,"EndxTruths[NMCTruths]/D");
  fNTuple->Branch("EndyTruths",EndyTruths,"EndyTruths[NMCTruths]/D");
  fNTuple->Branch("EndzTruths",EndzTruths,"EndzTruths[NMCTruths]/D");
  fNTuple->Branch("EndtTruths",EndtTruths,"EndtTruths[NMCTruths]/D");
  fNTuple->Branch("DeltaXTruths",DeltaXTruths,"DeltaXTruths[NMCTruths]/D");
  fNTuple->Branch("DeltaYTruths",DeltaYTruths,"DeltaYTruths[NMCTruths]/D");

}

SupernovaAna& SupernovaAna::operator= ( SupernovaAna const &  )

delete

SupernovaAna& SupernovaAna::operator= ( SupernovaAna &&  )

delete

void SupernovaAna::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 118 of file SupernovaAna_module.cc.

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

Member Data Documentation

double SupernovaAna::avgRMS

private

Definition at line 76 of file SupernovaAna_module.cc.

double SupernovaAna::avgRMSloop[kMaxNumHits]

private

Definition at line 77 of file SupernovaAna_module.cc.

double SupernovaAna::avgRMSloopsum

private

Definition at line 78 of file SupernovaAna_module.cc.

int SupernovaAna::CCNCTruths[kMaxNumParticles]

private

Definition at line 88 of file SupernovaAna_module.cc.

double SupernovaAna::DeltaXTruths[kMaxNumParticles]

private

Definition at line 100 of file SupernovaAna_module.cc.

double SupernovaAna::DeltaYTruths[kMaxNumParticles]

private

Definition at line 101 of file SupernovaAna_module.cc.

double SupernovaAna::EndtTruths[kMaxNumParticles]

private

Definition at line 99 of file SupernovaAna_module.cc.

double SupernovaAna::EndxTruths[kMaxNumParticles]

private

Definition at line 96 of file SupernovaAna_module.cc.

double SupernovaAna::EndyTruths[kMaxNumParticles]

private

Definition at line 97 of file SupernovaAna_module.cc.

double SupernovaAna::EndzTruths[kMaxNumParticles]

private

Definition at line 98 of file SupernovaAna_module.cc.

int SupernovaAna::EventNum

private

Definition at line 73 of file SupernovaAna_module.cc.

int SupernovaAna::EventNumCut

private

Definition at line 74 of file SupernovaAna_module.cc.

std::string SupernovaAna::fHitLabel

private

Definition at line 67 of file SupernovaAna_module.cc.

TTree* SupernovaAna::fNTuple

private

Definition at line 70 of file SupernovaAna_module.cc.

std::string SupernovaAna::fTruthLabel

private

Definition at line 66 of file SupernovaAna_module.cc.

double SupernovaAna::LeptonEnergyTruths[kMaxNumParticles]

private

Definition at line 85 of file SupernovaAna_module.cc.

int SupernovaAna::LeptonPDGCodeTruths[kMaxNumParticles]

private

Definition at line 91 of file SupernovaAna_module.cc.

int SupernovaAna::ModeTruths[kMaxNumParticles]

private

Definition at line 89 of file SupernovaAna_module.cc.

double SupernovaAna::MomentumYTruths[kMaxNumParticles]

private

Definition at line 82 of file SupernovaAna_module.cc.

double SupernovaAna::MomentumZTruths[kMaxNumParticles]

private

Definition at line 83 of file SupernovaAna_module.cc.

int SupernovaAna::NHits

private

Definition at line 75 of file SupernovaAna_module.cc.

int SupernovaAna::NMCTruths

private

Definition at line 80 of file SupernovaAna_module.cc.

int SupernovaAna::NParticles[kMaxNumParticles]

private

Definition at line 81 of file SupernovaAna_module.cc.

double SupernovaAna::NuEnergyTruths[kMaxNumParticles]

private

Definition at line 84 of file SupernovaAna_module.cc.

int SupernovaAna::NuPDGCodeTruths[kMaxNumParticles]

private

Definition at line 90 of file SupernovaAna_module.cc.

int SupernovaAna::OriginTruths[kMaxNumParticles]

private

Definition at line 87 of file SupernovaAna_module.cc.

double SupernovaAna::SummedADC

private

Definition at line 79 of file SupernovaAna_module.cc.

double SupernovaAna::ThetaTruths[kMaxNumParticles]

private

Definition at line 86 of file SupernovaAna_module.cc.

double SupernovaAna::tTruths[kMaxNumParticles]

private

Definition at line 95 of file SupernovaAna_module.cc.

double SupernovaAna::xTruths[kMaxNumParticles]

private

Definition at line 92 of file SupernovaAna_module.cc.

double SupernovaAna::yTruths[kMaxNumParticles]

private

Definition at line 93 of file SupernovaAna_module.cc.

double SupernovaAna::zTruths[kMaxNumParticles]

private

Definition at line 94 of file SupernovaAna_module.cc.

---

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

• duneana/duneana/SupernovaAna/SupernovaAna_module.cc