dump_event.C File Reference

#include "TBox.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH2.h"
#include "TSystem.h"
#include "stdio.h"
#include "LEM/func/Library.h"
#include <iostream>

Go to the source code of this file.

Functions
void	draw_event (TCanvas *theCanvas, TH2F *h[2], const MatchableEvent &sum)
void	printEventToJson (const MatchableEvent sum, FILE *, int idx=0)
void	dump_event (TString libPath, int idx=0, const char *outfile_name=NULL, bool draw=true, bool flipEven=false, bool flipOdd=false, bool idxs=false)
void	next_event ()
void	dump_lib (int start=1, int nEvts=1, const char *outfile_name=NULL, bool draw=false)

Variables
int	evt_idx =0
FILE *	outfile =NULL

Function Documentation

void draw_event	(	TCanvas *	theCanvas,
		TH2F *	h[2],
		const MatchableEvent &	sum
	)

Definition at line 76 of file dump_event.C.

References b, lem::kVertexCell, lem::kVertexPlane, lem::EventSummary::trueVtxCell, lem::EventSummary::trueVtxCellOther, lem::EventSummary::trueVtxPlane, and registry_explorer::v.

Referenced by adjusttick(), animation_step(), draw_whole_user_event(), dump_event(), and toggle_cum_ani().

                                                                          {
  for(int v = 0; v < 2; ++v){
    theCanvas->cd(v+1);
    h[v]->Draw("colz");

    TBox* b = new TBox(lem::kVertexPlane,   lem::kVertexCell,
                       lem::kVertexPlane+1, lem::kVertexCell+1);
    b->SetLineColor(kRed);
    b->SetFillStyle(0);
    b->SetLineWidth(2);
    b->SetLineStyle(7);
    b->Draw();

    const int trueCell = (sum.trueVtxPlane%2 == v) ? sum.trueVtxCellOther : sum.trueVtxCell;
    TBox* b2 = new TBox(sum.trueVtxPlane,   trueCell,
                        sum.trueVtxPlane+1, trueCell+1);
    b2->SetLineColor(kGreen+2);
    b2->SetFillStyle(0);
    b2->SetLineWidth(2);
    b2->SetLineStyle(7);
    b2->Draw();
  }
}

void dump_event	(	TString	libPath,
		int	idx = 0,
		const char *	outfile_name = NULL,
		bool	draw = true,
		bool	flipEven = false,
		bool	flipOdd = false,
		bool	idxs = false
	)

Definition at line 18 of file dump_event.C.

References draw_event(), lem::Library::Event(), evt_idx, lem::EventSummary::FillHists(), lem::Library::FromMMap(), make_syst_table_plots::h, MECModelEnuComparisons::i, compare_h5_caf::idx, allTimeWatchdog::index, outfile, printEventToJson(), make_template_knob_config::stdout, and sum.

Referenced by dump_lib(), and next_event().

{
  static Library* lib = NULL;
  static TCanvas* theCanvas = NULL;
  static int views = 2;
  static TH2F* h[2] = {NULL,NULL};
  long long int index=idx;
  evt_idx=idx;

  //FILE* outfile=NULL;

  if(outfile==NULL){
    if(outfile_name!=NULL){
      outfile = fopen(outfile_name,"w");
    }else{
      outfile = stdout;
    }
  }

  // Map the Library File if we haven't already mapped it
  if(lib==NULL){
    // Open the Library via Memory Map
    lib = Library::FromMMap(libPath.Data(), false);
  }

  // Create a Canvase for drawing the events
  if(theCanvas==NULL){
    theCanvas = new TCanvas("a", "b", 1200, 600);
    theCanvas->Divide(2, 1);
  }

  // Retrieve the template from the library
  const MatchableEvent& sum = lib->Event(idx);

  printEventToJson(sum,outfile,index); // Write the event out in JSon

  if(draw){

    // Book the histograms for drawing the event or clear them if they exist
    if(h[0]==NULL){
      for(int i=0; i<2; ++i){
        h[i] = new TH2F("", ";plane;cell", 128, 0, 128, 128, 64, 192);
      } // endfor
    }else{
      for(int i=0; i<2; ++i){
        // Clear/Reset the Histograms
        h[i]->Reset();
      } // endfor
    } //endif

    // Fill the histograms (even view first, then the odd view)
    sum.FillHists(h[0], h[1], flipEven, flipOdd, idxs);

    draw_event(theCanvas, h, sum);  // Draw the Event in question
  }

};

void dump_lib	(	int	start = 1,
		int	nEvts = 1,
		const char *	outfile_name = NULL,
		bool	draw = false
	)

Definition at line 138 of file dump_event.C.

References dump_event(), fclose(), MECModelEnuComparisons::i, nEvts, outfile, febshutoff_auto::start, and make_template_knob_config::stdout.

                                                                                       {

  // Close any open output files
  if(outfile!=NULL){
    fclose(outfile);
    outfile==NULL;
  }

  // Open a new file for output
  if(outfile_name!=NULL){
    outfile = fopen(outfile_name,"w");
  }else{
    outfile = stdout;
  }

  fprintf(outfile,"{ \"Library\": [");
  for(int i=start; i<(start+nEvts); ++i){
    dump_event("",i,outfile_name,false);
    if(i!=(start+nEvts-1)){
        fprintf(outfile,",\n");
      }
  }
  fprintf(outfile,"]\n }\n");
  fflush(outfile);
  fclose(outfile);
  outfile=NULL;
  return;
};

void next_event

(

)

Definition at line 132 of file dump_event.C.

References dump_event(), and evt_idx.

                 {
  evt_idx++;
  dump_event("",evt_idx);
};

void printEventToJson	(	const MatchableEvent	sum,
		FILE *	theOutput,
		int	idx = 0
	)

Definition at line 100 of file dump_event.C.

References lem::EventSummary::ccnc, getBrightness::cell, lem::LiteHit::Cell(), energy, lem::EventSummary::hits, MECModelEnuComparisons::i, lem::EventSummary::mode, lem::EventSummary::nhits, lem::EventSummary::nParts, lem::EventSummary::partPdgs, lem::EventSummary::pdg, lem::LiteHit::pecorr, lem::EventSummary::photE0, lem::EventSummary::photE1, lem::EventSummary::photL0, lem::EventSummary::photL1, lem::LiteHit::Plane(), NDAPDHVSetting::plane, lem::EventSummary::totalPE, lem::EventSummary::trueEVis, and POTSpillRate::view.

Referenced by dump_event().

                                                                           {
  
  // Print out the event's information
  fprintf(theOutput,"{\"Event\": {   \n\"EventID\": %lld,\n",index);
  fprintf(theOutput,"  \"Truth\": {\n");
  fprintf(theOutput,"       \"PE\": %2.2f ,\n"   ,sum.totalPE);
  fprintf(theOutput,"       \"VisE\": %2.2f ,\n" ,sum.trueEVis);
  fprintf(theOutput,"       \"CC\": %d ,\n"   ,sum.ccnc);
  fprintf(theOutput,"       \"PDG\": %d ,\n"   ,sum.pdg);
  fprintf(theOutput,"       \"Mode\": %d ,\n"   ,sum.mode);
  fprintf(theOutput,"       \"PhotonL\": [ %2.2f, %2.2f ],\n",sum.photL0,sum.photL1);
  fprintf(theOutput,"       \"PhotonE\": [ %2.2f, %2.2f ],\n",sum.photE0,sum.photE1);

  fprintf(theOutput,"       \"ParticlesPDG\": [ ");
  for(int i=0; i< sum.nParts; ++i){
    fprintf(theOutput,"%d,",sum.partPdgs[i]);
  }
  fprintf(theOutput," ]\n      },\n");
  fprintf(theOutput,"   \"NHits\": %d,\n",sum.nhits);
  fprintf(theOutput,"   \"HitList\": [\n");
  for(unsigned int i =0; i<sum.nhits; ++i){
    if(i!=0){fprintf(theOutput,",\n");}
    const int plane  = sum.hits[i].Plane();
    const int cell   = sum.hits[i].Cell();
    const int view   = plane%2;
    const float energy = sum.hits[i].pecorr;
    fprintf(theOutput,"      [ %3d, %3d, %3d, %7.4f ]", plane, cell, view, energy);
  }
  fprintf(theOutput,"\n      ]\n");  
  fprintf(theOutput,"   }\n}\n");
};

Variable Documentation

int evt_idx =0

Definition at line 12 of file dump_event.C.

Referenced by dump_event(), and next_event().

FILE* outfile =NULL

Definition at line 13 of file dump_event.C.

Referenced by add_and_normalize_LLhistos(), novaddt::Evaluation::analyze(), air::AirKalmanAna::analyze(), BestFit(), BuildMetricsTree(), BuildMetricsTree_OnMon(), BuildOldTree(), CAF_makeCAFSensitivities_for_FNEX(), cafe_FD_predictions(), CalcFiberLoopCorr(), compare_cos_numi(), ConsolidateSystMaker(), CreateKNNTrainingTrees(), d3sigma_calc(), d4sigma_calc(), d4sigma_hist(), d4sigma_plot(), demoBarChart(), demoFitContours(), demoFitSlices(), doUnfolding(), dump_event(), dump_lib(), calib::AttenFit::endJob(), energyResolution(), generate_fake_data(), get_cosmic_sample(), getHists_FNEX(), joint_fit_2017_contours(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_bestfit_univ(), joint_fit_future_contour_univ(), jointsensitivity(), main(), make_BPF_Energy_Estimator_KNN_regression(), make_BPF_muon_PID_KNN(), make_DataMCComp_numu(), make_files_for_decomp_plots(), MakeCAFSensitivities_for_FNEX(), MakeISysts(), MakeNus17CosBkgd(), MakeNus18CosBkgd(), MakeNus2019CovMxSysts(), MakeSurface(), MakeSurfaceNoNDOsc(), MakeSysts(), makeXSecPlots_TemplateFit(), NDFDCVNBDTCutSelector(), neut_code_from_rootracker(), nue_fnex_vs_caf(), nue_fnex_vs_caf_noextrapolation(), evdb::ScanWindow::OpenOutputFile(), INukeNucleonCorr::OutputFiles(), plot_BEN(), readTree(), evdb::ScanFrame::Record(), reduce_the_predfile_nue_numu(), run_joint_fit_2020_bestfit(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), sensitivity2018(), sensitivity2020(), specprod_systematics(), Tutorial2019FitContours(), Tutorial2019FitSlices(), validation(), calib::AttenuationFit::writeResults(), and fnex::Correlations::~Correlations().