anadst.cxx

Go to the documentation of this file.

#include <TFile.h>
#include <TH2F.h>
#include "DSTAnalysis/NuMI/NuMIRunQAna.h"
#include "DSTAnalysis/NuMI/NuMIBeamAna.h"
#include "DSTAnalysis/NuMI/NuMITgtAna.h"
#include "MIPPEventSummary/MIPPEventSummary.h"
#include "MIPPEventSummary/MIPPTrackSummary.h"

static int firstRun = 14497;
static int lastRun  = 15250;

static NuMIRunQAna* fNuMIRunQAna = 0;
static NuMIBeamAna* fNuMIBeamAna = 0;
static NuMITgtAna*  fNuMITgtAna  = 0;

// ROOT histograms
static TFile* gsHistoFile = 0;

extern "C" 
{
  void InitDST(const char* fname) 
  {
    gsHistoFile = new TFile(fname,"RECREATE");
    gsHistoFile->cd();
    
    fNuMIRunQAna = new NuMIRunQAna();
    fNuMIBeamAna = new NuMIBeamAna();
    fNuMITgtAna  = new NuMITgtAna();
 }

  //....................................................................

  void EndDST() { gsHistoFile->Write(); }

  //....................................................................

  bool AnaDST(const MIPPEventSummary* dst) 
  {
    // Select only NuMI target data
    //    if (dst->tgta!=1201) return false;

    // Fill histograms related to run quality
    bool goodrun = fNuMIRunQAna->Ana(dst);
    if (goodrun==false) return false;

    // Look at the beam track quality
    bool goodbeam = fNuMIBeamAna->Ana(dst);
    if (goodbeam==false) return false;
    
    // Now fill analysis histograms
    bool goodevent = fNuMITgtAna->Ana(dst);
    if (goodevent==false) return false;

    return true;
  }

}
////////////////////////////////////////////////////////////////////////

---