ana::SBNOsc::ROC Struct Reference

#include <ROC.h>

Classes
struct	NormalizedPrimitive
struct	Primitive

Public Member Functions
void	Initialize ()
void	Fill (const Cuts &cuts, const numu::RecoEvent &event, bool file_is_neutrino)
void	Write () const
void	BestCuts () const
void	Normalize (float neutrino_scale, float cosmic_scale)

Public Attributes
NormalizedPrimitive	crt_track_angle
NormalizedPrimitive	crt_hit_distance
std::vector < NormalizedPrimitive * >	fAllPrimitives

Detailed Description

Class to handle making of ROC plots for cut values

Definition at line 18 of file ROC.h.

Member Function Documentation

void ana::SBNOsc::ROC::BestCuts

(

)

const

Get the optimal cut values by S/sqrt(S+B)

Definition at line 24 of file ROC.cc.

                                  {
  for (const NormalizedPrimitive *prim: fAllPrimitives) {
    std::cout << "Cut (" << prim->name << ") value: " << prim->BestCut() << std::endl;
  }
}

void ana::SBNOsc::ROC::Fill	(	const Cuts &	cuts,
		const numu::RecoEvent &	event,
		bool	file_is_neutrino
	)

Fill the class with a reconstructed event

Parameters

event

The RecoEvent information

Definition at line 30 of file ROC.cc.

                                                                                                    {
  for (unsigned i = 0; i < event.reco.size(); i++) {
    const numu::RecoInteraction &reco = event.reco[i];
    bool is_signal = reco.match.mode == numu::mCC;
    bool is_bkg = reco.match.mode == numu::mCosmic || reco.match.mode == numu::mIntimeCosmic;
    if (!is_signal && !is_bkg) continue; // ignore NC stuff for now
    
    const numu::RecoTrack &track = event.tracks.at(reco.slice.primary_track_index);

    std::array<bool, ana::SBNOsc::Cuts::nCuts> results = cuts.ProcessRecoCuts(event, i, true);

    // only apply to events that already pass flash matching
    if (results[5] /* flashmatch */) {
      // CRT stuff
      if (track.crt_match.track.present) {
        if (file_is_neutrino) crt_track_angle.FillNeutrino(is_signal, track.crt_match.track.angle);
        else crt_track_angle.FillCosmic(is_signal, track.crt_match.track.angle);
      }
      else {
        if (file_is_neutrino) crt_track_angle.FillNeverNeutrino(is_signal);
        else crt_track_angle.FillNeverCosmic(is_signal);
      }
      if (track.crt_match.hit_match.present && !cuts.TimeInSpill(track.crt_match.hit_match.time)) {
        if (file_is_neutrino) crt_hit_distance.FillNeutrino(is_signal, track.crt_match.hit_match.distance);
        else crt_hit_distance.FillCosmic(is_signal, track.crt_match.hit_match.distance);
      }
      else {
        if (file_is_neutrino) crt_hit_distance.FillNeverNeutrino(is_signal);
        else crt_hit_distance.FillNeverCosmic(is_signal);
      }
    }
  }
    
}

void ana::SBNOsc::ROC::Initialize

(

)

Initialize the class

Definition at line 5 of file ROC.cc.

                              {
  crt_track_angle.Initialize("crt_track_angle", 0, 6, 60);
  fAllPrimitives.push_back(&crt_track_angle);
  crt_hit_distance.Initialize("crt_hit_distance", 0, 100, 100);
  fAllPrimitives.push_back(&crt_hit_distance);
}

void ana::SBNOsc::ROC::Normalize	(	float	neutrino_scale,
		float	cosmic_scale
	)

Definition at line 18 of file ROC.cc.

                                                                     {
  for (NormalizedPrimitive *prim: fAllPrimitives) {
    prim->Normalize(neutrino_scale, cosmic_scale);
  }
}

void ana::SBNOsc::ROC::Write

(

)

const

Write the ROC plots to disk

Definition at line 12 of file ROC.cc.

                               {
  for (const NormalizedPrimitive *prim: fAllPrimitives) {
    prim->Write();
  }
}

Member Data Documentation

NormalizedPrimitive ana::SBNOsc::ROC::crt_hit_distance

Definition at line 78 of file ROC.h.

NormalizedPrimitive ana::SBNOsc::ROC::crt_track_angle

Definition at line 77 of file ROC.h.

std::vector<NormalizedPrimitive *> ana::SBNOsc::ROC::fAllPrimitives

Definition at line 79 of file ROC.h.

---

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

• ROC.h
• ROC.cc