ana::SBNOsc::Normalize Class Reference

#include <Normalize.h>

Public Member Functions
void	Initialize (const fhicl::ParameterSet &fcl)
void	AddCosmicEvent (const event::Event &event)
void	AddNeutrinoEvent (const event::Event &event)
void	AddNeutrinoSubRun (const SubRun &subrun)
void	AddCosmicFile (const FileMeta &meta)
double	ScaleNeutrino (double goal_pot) const
double	ScaleCosmic (double goal_pot) const

Private Attributes
int	fLastCosmicFileID
int	fLastCosmicEventID
int	fLastNeutrinoFileID
int	fLastNeutrinoEventID
double	fPotPerReadout
double	fNeutrinoPOT
unsigned	fCosmicEventsPerFile
unsigned	fNNeutrinoEvents
unsigned	fNCosmicEvents

Detailed Description

Class to handle normalization of overlay and intime-cosmic samples.

Definition at line 15 of file Normalize.h.

Member Function Documentation

void ana::SBNOsc::Normalize::AddCosmicEvent

(

const event::Event &

event

)

Add a cosmic event to the normalization

Parameters

event

The sbncode core event information

Definition at line 23 of file Normalize.cc.

{}

void ana::SBNOsc::Normalize::AddCosmicFile

(

const FileMeta &

meta

)

Definition at line 19 of file Normalize.cc.

                                                           {
  fNCosmicEvents += meta.n_gen_events;  
}

void ana::SBNOsc::Normalize::AddNeutrinoEvent

(

const event::Event &

event

)

Add a neutrino event to the normalization

Parameters

event

The sbncode core event information

Definition at line 25 of file Normalize.cc.

                                                                 {
  fNNeutrinoEvents ++;
}

void ana::SBNOsc::Normalize::AddNeutrinoSubRun

(

const SubRun &

subrun

)

Add a neutrino subrunto the normalization

Parameters

subrun

The sbncode core Subrun information

Definition at line 29 of file Normalize.cc.

                                                               {
  std::cout << "New subrun: " << subrun.totgoodpot << std::endl;
  fNeutrinoPOT += subrun.totgoodpot;
}

void ana::SBNOsc::Normalize::Initialize

(

const fhicl::ParameterSet &

fcl

)

Initialize the class

Parameters

fcl	The fhicl configuration

Definition at line 5 of file Normalize.cc.

                                                                {
  fPotPerReadout = fcl.get<double>("PotPerReadout", 5.e12);

  fLastCosmicFileID = -1;
  fLastCosmicEventID = -1;

  fLastNeutrinoFileID = -1;
  fLastNeutrinoEventID = -1;

  fNeutrinoPOT = 0.;
  fNNeutrinoEvents = 0;
  fNCosmicEvents = 0;
}

double ana::SBNOsc::Normalize::ScaleCosmic

(

double

goal_pot

)

const

Scale the cosmic events to a certain POT. To be called after all events are processed

Parameters

goal_pot

The POT to scale the sample to.

Returns

The scale factor for each cosmic event to the set POT

Definition at line 34 of file Normalize.cc.

                                                            {

  std::cout << "Total N cosmics: " << fNCosmicEvents << std::endl;
  std::cout << "Total N neutrinos: " << fNNeutrinoEvents << std::endl;
  std::cout << "Total neutrino POT: " << fNeutrinoPOT << std::endl;

  // fraction of all POT accounted for by the neutrino sample
  double neutrino_per_pot = 0.;
  if (fNeutrinoPOT > 1e-3 && fNNeutrinoEvents > 0) {
    neutrino_per_pot = fNNeutrinoEvents / (fNeutrinoPOT);
  }

  std::cout << "Neutrino per POT: " << neutrino_per_pot << std::endl;

  // POT per cosmic event
  double cosmic_per_pot = (1./fPotPerReadout) - neutrino_per_pot;

  std::cout << "Cosmic per POT: " << cosmic_per_pot << std::endl;

  // total POT equaivalent of cosmic
  double cosmic_pot = fNCosmicEvents / cosmic_per_pot;

  std::cout << "Cosmic POT: " << cosmic_pot << std::endl;
  
  // scale
  return goal_pot / cosmic_pot;
}

double ana::SBNOsc::Normalize::ScaleNeutrino

(

double

goal_pot

)

const

Scale the neutrino events to a certain POT. To be called after all events are processed

Parameters

goal_pot

The POT to scale the sample to.

Returns

The scale factor for each neutrino event to the set POT

Definition at line 62 of file Normalize.cc.

                                                              {
  if (fNeutrinoPOT < 1e-3 || fNNeutrinoEvents == 0) return 0;
  // simpler -- just scale neutrino POT to goal
  return goal_pot / (fNeutrinoPOT);
}

Member Data Documentation

unsigned ana::SBNOsc::Normalize::fCosmicEventsPerFile

private

Definition at line 68 of file Normalize.h.

int ana::SBNOsc::Normalize::fLastCosmicEventID

private

Definition at line 61 of file Normalize.h.

int ana::SBNOsc::Normalize::fLastCosmicFileID

private

Definition at line 60 of file Normalize.h.

int ana::SBNOsc::Normalize::fLastNeutrinoEventID

private

Definition at line 64 of file Normalize.h.

int ana::SBNOsc::Normalize::fLastNeutrinoFileID

private

Definition at line 63 of file Normalize.h.

unsigned ana::SBNOsc::Normalize::fNCosmicEvents

private

Definition at line 70 of file Normalize.h.

double ana::SBNOsc::Normalize::fNeutrinoPOT

private

Definition at line 67 of file Normalize.h.

unsigned ana::SBNOsc::Normalize::fNNeutrinoEvents

private

Definition at line 69 of file Normalize.h.

double ana::SBNOsc::Normalize::fPotPerReadout

private

Definition at line 66 of file Normalize.h.

---

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

• Normalize.h
• Normalize.cc