demo6.C File Reference

#include "CAFAna/Core/Binning.h"
#include "CAFAna/Cuts/Cuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "CAFAna/Core/EventList.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/Vars.h"
#include "StandardRecord/Proxy/SRProxy.h"
#include "TCanvas.h"
#include "TH2.h"
#include "TStyle.h"
#include "TLegend.h"
#include "TLatex.h"

Go to the source code of this file.

Functions
void	Preliminary ()
	Put NOvA Preliminary on plots. More...
void	legend (TH1 *data, TH1 *genie, TH1 *cosmic)
	legend to depict Data vs. GENIE MC vs. COSMIC MC More...
void	demo6 ()

Function Documentation

void demo6

(

)

Remove statistics from plots

Data, GENIE, Cosmic MC files

These are just place-holder histograms for ease of use with the TLegend

Spectrum to be filled from the loader List spectra we'd like to draw along with the cuts we'd like to apply Here we're defining cut to be on kCleanup

Definition at line 50 of file demo6.C.

References ana::bins, caf::Proxy< caf::SRSlice >::boxmax, caf::Proxy< caf::SRSlice >::boxmin, caf::Proxy< caf::SRTrackBranch >::cosmic, ana::DataMCComparisonAreaNormalized(), ana::SpectrumLoader::Go(), plot_xsec_1d::hgenie, Integral(), caf::Proxy< caf::SRTrackBranch >::kalman, kBlue, ana::kCleanup, ana::kCosmic, ana::kNContPlanes, ana::kNHit, ana::kNumuContainFD, kRed, demo0::kTrackLen, legend(), ana::MakeEventListFile(), caf::Proxy< caf::StandardRecord >::mc, caf::Proxy< caf::SRSlice >::ncontplanes, caf::Proxy< caf::SRSlice >::nhit, caf::Proxy< caf::SRTruthBranch >::nnu, caf::Proxy< caf::SRTrackBase >::ntracks, caf::Proxy< caf::SRKalman >::ntracks, ana::Spectrum::POT(), Preliminary(), ana::Binning::Simple(), caf::Proxy< caf::StandardRecord >::slc, sr, string, ana::Spectrum::ToTH1(), ana::Spectrum::ToTH2(), caf::Proxy< caf::SRTrackBase >::tracks, caf::Proxy< caf::SRKalman >::tracks, caf::Proxy< caf::StandardRecord >::trk, caf::Proxy< caf::SRVector3D >::X(), caf::Proxy< caf::SRVector3D >::Y(), and caf::Proxy< caf::SRVector3D >::Z().

{
  /// Remove statistics from plots
  gStyle->SetOptStat(0);

  /// Data, GENIE, Cosmic MC files
  const std::string fData = "$NOVA_PROD/data/FarDet/S13-10-11/numi/hadd/fardet_r11761*v3a*caf.root";
  const std::string fGENIE = "$NOVA_PROD/mc/B13-10-23/genie/fd/ImprovedTrans_newMultiplexReadoutSim/fd_r01_s??_*fhc*caf.root";
  const std::string fCOSMICS = "$NOVA_PROD/mc/B13-10-23/cosmics/fd/ImprovedTrans_newMultiplexReadoutSim/fd_r01_s??_*cosmics*caf.root";

  //-------------------------------------------------------------------------
  /// These are just place-holder histograms for ease of use with the TLegend
  TH1 *hdata = new TH1("hdata", "hdata", 10, 0, 1);
  TH1 *hgenie = new TH1("hgenie", "hgenie", 10, 0, 1);
  TH1 *hcosmic = new TH1("hcosmic", "hcosmic", 10, 0, 1);
  
  hdata->SetLineWidth(2);
  hgenie->SetLineWidth(2);
  hcosmic->SetLineWidth(2);

  hdata->SetLineColor(kBlack);
  hgenie->SetLineColor(kRed);
  hcosmic->SetLineColor(kBlue);
  
  hdata->SetMarkerStyle(kFullCircle);
  //------------------------------------------------------------------------
  

  //------------------------------------------------------------------------
  // Specify variables needed and arbitrary code to extract value from
  // SRProxy
  const Var kBoxMinX( [](const caf::SRProxy* sr)
                      {
                        return sr->slc.boxmin.X();
                      });
  // Don't use this one, but is another example of a Var definition
  // with the lambda function
  const Var kNContPlanes([](const caf::SRProxy* sr) 
                         {
                           if(sr->slc.ncontplanes == 0) return (unsigned int)0.0;
                           return sr->slc.ncontplanes;
                         });

  const Var kHitRatio( [](const caf::SRProxy* sr)
                       {
                         if(sr->trk.cosmic.ntracks == 0) return 0.0;
                         return ((sr->trk.cosmic.tracks[0].nhit*1.0)/(1.0*sr->slc.nhit));
                       });

  // Cut variable on hit ratio
  const Cut kHitRatioCut([](const caf::SRProxy* sr)
                         {
                           return ((sr->trk.cosmic.tracks[0].nhit*1.0)/(1.0*sr->slc.nhit)) > 0.88;
                         });

  const Cut kCleanup([](const caf::SRProxy* sr)
                     {
                       return (sr->slc.nhit > 20 &&
                               sr->trk.cosmic.tracks[0].view == 2);
                     });
  
  const Cut kNNu([](const caf::SRProxy* sr)
                 {
                   return (sr->mc.nnu > 0);
                 });
  
  const Cut k2DBContain(  // If we pretend there's no Cosmics, this works well
                      [](const caf::SRProxy* sr)
                      {
                        return (sr->slc.boxmin.X() > -700 &&
                                sr->slc.boxmin.Y() > -700 &&
                                sr->slc.boxmax.X() < 700 &&
                                sr->slc.boxmax.Y() < 660 &&
                                sr->slc.boxmin.Z() > 100 &&
                                sr->slc.boxmax.Z() < 800);
                      });

  const Var kTrackLen([](const caf::SRProxy* sr)
                      {
                        if(sr->trk.kalman.ntracks == 0) return 0.0f;
                        return sr->trk.kalman.tracks[0].len;
                      });

  const Var kNHit([](const caf::SRProxy* sr){return sr->slc.nhit;});

  
  //SpectrumLoader definitions
  SpectrumLoader loadData(fData);
  SpectrumLoader loadGENIE(fGENIE);
  SpectrumLoader loadCOSMICS(fCOSMICS,kCosmic);

  //Any particular binning types for spectra
  const Binning bins   = Binning::Simple(100, 0, 1000);

  const Binning binsX  = Binning::Simple(425, -850, 850);

  //--------------------------------------------------------------------------
  /// Spectrum to be filled from the loader
  /// List spectra we'd like to draw along with the cuts we'd like to apply
  /// Here we're defining cut to be on kCleanup
  Spectrum boxminxD("Slice Box Min X (cm)", binsX, loadData, kBoxMinX,kCleanup); // DATA
  Spectrum boxminxG("Slice Box Min X (cm)", binsX, loadGENIE, kBoxMinX, kNNu && kCleanup); // GENIE MC
  Spectrum boxminxC("Slice Box Min X", binsX, loadCOSMICS, kBoxMinX, kCleanup); // COSMIC MC

  // Nhits vs Track Length 2D Spectrum
  Spectrum nhittracklen("Nhits vs Track Length", loadCOSMICS, 
                        bins, kTrackLen, bins, kNHit,kCleanup && kNumuContainFD); // COSMIC MC

  // Do it!
  loadGENIE.Go();
  loadData.Go();
  loadCOSMICS.Go();

  //Pick data set we'd like to print run/subrun/event information to text file
  //Cut on genie events with a neutrino that pass kCleanup and k2DBEnter
  MakeEventListFile(fGENIE,kNNu && kCleanup,"event-list-genie.txt"); 
  //Same cuts to a file for Data - obviously can't rely on truth for a neutrino
  MakeEventListFile(fData,kCleanup,"event-list-data.txt"); 
//Same cuts to a file for Data - obviously can't rely on truth for a neutrino
  MakeEventListFile(fCOSMICS,kCleanup,"event-list-cosmics.txt"); 

  //--------------------------------------------------------------------------
  //We have histograms

  // Nhits vs TrackLength

  
  TH2* hnht = nhittracklen.ToTH2(nhittracklen.POT());
  hnht->Draw("colz");
  Preliminary();
  
  // Slice Box Min X
  TCanvas *minx = new TCanvas();
  //Make Data vs. MC area normalized
  DataMCComparisonAreaNormalized(boxminxD,boxminxG);
  //Add Cosmic MC for comparison
  TH1* boxminxc = boxminxC.ToTH1(boxminxC.POT());
  boxminxc->SetLineColor(kBlue);
  //Scale Cosmic MC appropriately
  boxminxc->Scale(boxminxD.ToTH1(boxminxD.POT())->Integral()/boxminxc->Integral());
  boxminxc->Draw("hist same");
  minx->SetLogy();
  //Add Preliminary tag and legend
  Preliminary();
  legend(hdata,hgenie,hcosmic);

}

void legend	(	TH1 *	data,
		TH1 *	genie,
		TH1 *	cosmic
	)

legend to depict Data vs. GENIE MC vs. COSMIC MC

Definition at line 37 of file demo6.C.

References MECModelEnuComparisons::leg.

Referenced by demo6().

                                               {
  
  TLegend *leg = new TLegend(.50,.57,.84,.84);
  leg->AddEntry(data,"Far Detector Data","lpe"); 
  leg->AddEntry(genie,"GENIE Simulation","l");   
  leg->AddEntry(cosmic,"Cosmic Simulation","l");
  leg->SetBorderSize(0); //no border for legend
  leg->SetFillColor(0);  //fill colour is white
  leg->Draw();

}

void Preliminary

(

)

Put NOvA Preliminary on plots.

Definition at line 26 of file demo6.C.

References kBlue, and prelim.

Referenced by demo6().

{
  TLatex* prelim = new TLatex(.9, .95, "NO#nuA Preliminary");
  prelim->SetTextColor(kBlue);
  prelim->SetNDC();
  prelim->SetTextSize(2/30.);
  prelim->SetTextAlign(32);
  prelim->Draw();
}