← Back to branch summary

~lbno-nd/lbno-nd/trunk

« back to all changes in this revision

Viewing changes to workdir/analysis/scripts/reconHists.cxx

  • Committer: tstainer
  • Date: 2014-06-16 14:50:45 UTC
  • Revision ID: tstainer@liv.ac.uk-20140616145045-6hb0uouxkh8pcp2r
added readme and updates

Show diffs side-by-side

added added

removed removed

Lines of Context:
workdir/analysis/scripts/reconHists.cxx
17
17
#include "TPad.h"
18
18
#include "TInterpreter.h"
19
19
#include "TColor.h"
 
20
#include "TGraph.h"
20
21
 
21
22
int binning = 4;
22
23
double maxenergy = 30.;
34
35
 
35
36
void reconHists(){
36
37
  TString loadfile = filename + fidCut + "cmFidCut" + extension;
 
38
  loadfile = "../../test0_recon.root";
37
39
 
38
40
  TFile file(loadfile);
39
41
 
69
71
  TH2F * tracking_eDepVsMom = (TH2F*)file.Get("tracking_EDepVsMom_h2d");
70
72
  TH2F * tracking_fullyContained_truthVsRecon_mom = (TH2F*)file.Get("tracking_fullyContained_Mom_truthVsRecon_h2d");
71
73
 
 
74
  TH1F * tracking_piZero_reconEnergy_all = (TH1F*)file.Get("tracking_piZero_reconEnergy_all_h1d");
 
75
  TH1F * tracking_piZero_reconEnergy_1pi0 = (TH1F*)file.Get("tracking_piZero_reconEnergy_1pi0_h1d");
 
76
  TH1F * tracking_piZero_reconEnergy_2pi0 = (TH1F*)file.Get("tracking_piZero_reconEnergy_2pi0_h1d");
 
77
  TH1F * tracking_piZero_reconEnergy_3pi0 = (TH1F*)file.Get("tracking_piZero_reconEnergy_3pi0_h1d");
 
78
 
 
79
  TH1F * tracking_piZero_truthCount_all = (TH1F*)file.Get("tracking_piZero_truthCount_all_h1d");
 
80
  TH1F * tracking_piZero_reconCount_all = (TH1F*)file.Get("tracking_piZero_reconCount_all_h1d");
 
81
 
 
82
  if(tracking_piZero_reconEnergy_all)tracking_piZero_reconEnergy_all->Rebin(binning);
 
83
  tracking_piZero_reconEnergy_1pi0->Rebin(binning);
 
84
  tracking_piZero_reconEnergy_2pi0->Rebin(binning);
 
85
  tracking_piZero_reconEnergy_3pi0->Rebin(binning);
 
86
 
 
87
  tracking_piZero_truthCount_all->Rebin(binning);
 
88
  tracking_piZero_reconCount_all->Rebin(binning);
 
89
 
72
90
  setupHistTH1(neutrino_truth_energy,kBlack);
73
91
  setupHistTH1(neutrino_recon_mom,kRed);
74
92
  setupHistTH1(neutrino_recon_mom_ccqe,kGreen);
344
362
  tracking_rangeVsEDep->DrawCopy("colz");
345
363
 
346
364
  c8->cd();
347
 
  c8->Divide(1,2);
 
365
  c8->Divide(2,2);
348
366
  if(logZ)gPad->SetLogz();
349
367
  gPad->SetGridx();
350
368
  gPad->SetGridy();
366
384
  tracking_fullyContained_truthVsRecon_mom->GetXaxis()->SetTitle("Truth Momentum [GeV/c]");
367
385
  tracking_fullyContained_truthVsRecon_mom->GetYaxis()->SetTitle("Reconstructed Momentum [GeV/c]");
368
386
  tracking_fullyContained_truthVsRecon_mom->DrawCopy("colz");
369
 
 
370
 
 
371
 
  c9->cd();
 
387
  c8->cd(3);
372
388
  gStyle->SetOptFit(1111);
373
 
  c9->Divide(1,2);
374
 
  c9->cd(1);
375
389
  if(logZ)gPad->SetLogz();
376
390
  gPad->SetGridx();
377
391
  gPad->SetGridy();
378
392
  gPad->SetTickx();
379
393
  gPad->SetTicky();
380
394
  gPad->Update();
381
 
  c9->cd(1);
382
395
  neutrino_reconDiff_mom->GetXaxis()->SetTitle("(Truth - Recon)/Truth Momentum [GeV/c]");
383
396
  neutrino_reconDiff_mom->GetYaxis()->SetTitle("Normalised units");
384
397
  neutrino_reconDiff_mom->SetTitle("(Truth - recon)/truth neutrino momentum");
385
398
  //neutrino_reconDiff_mom->Fit("gaus","","",-0.5,0.2);
386
399
  neutrino_reconDiff_mom->DrawCopy();
387
 
  c9->cd(2);
 
400
  c8->cd(4);
388
401
  if(logZ)gPad->SetLogz(); 
389
402
  gPad->SetGridx();
390
403
  gPad->SetGridy();
417
430
  neutrino_reconDiff_mom_ccqe->DrawCopy(); 
418
431
  neutrino_reconDiff_mom_ccqe_calc->DrawCopy("SAME");
419
432
 
 
433
 
 
434
  c9->cd();
 
435
  c9->Divide(2,2);
 
436
  gPad->SetGridx();
 
437
  gPad->SetGridy();
 
438
  gPad->SetTickx();
 
439
  gPad->SetTicky();
 
440
  gPad->Update();
 
441
  c9->cd(1);
 
442
  if(tracking_piZero_reconEnergy_all)tracking_piZero_reconEnergy_all->DrawCopy();
 
443
  c9->cd(2);
 
444
  gPad->SetGridx();
 
445
  gPad->SetGridy();
 
446
  gPad->SetTickx();
 
447
  gPad->SetTicky();
 
448
  gPad->Update();
 
449
 
 
450
  tracking_piZero_reconEnergy_1pi0->GetXaxis()->SetTitle("recon energy [GeV]");
 
451
  tracking_piZero_reconEnergy_1pi0->GetYaxis()->SetTitle("count");
 
452
 
 
453
  tracking_piZero_reconEnergy_1pi0->SetLineWidth(2);
 
454
  tracking_piZero_reconEnergy_2pi0->SetLineWidth(2);
 
455
  tracking_piZero_reconEnergy_3pi0->SetLineWidth(2);
 
456
  tracking_piZero_reconEnergy_1pi0->SetLineColor(kBlack);
 
457
  tracking_piZero_reconEnergy_2pi0->SetLineColor(kGreen);
 
458
  tracking_piZero_reconEnergy_3pi0->SetLineColor(kRed);
 
459
 
 
460
  tracking_piZero_reconEnergy_1pi0->DrawCopy();
 
461
  tracking_piZero_reconEnergy_2pi0->DrawCopy("SAME");
 
462
  tracking_piZero_reconEnergy_3pi0->DrawCopy("SAME");
 
463
 
 
464
  c9->cd(3);
 
465
  gPad->SetGridx();
 
466
  gPad->SetGridy();
 
467
  gPad->SetTickx();
 
468
  gPad->SetTicky();
 
469
  gPad->Update();
 
470
  tracking_piZero_truthCount_all->GetXaxis()->SetTitle("neutrino energy [GeV]");
 
471
  tracking_piZero_truthCount_all->GetYaxis()->SetTitle("count");
 
472
  tracking_piZero_truthCount_all->SetLineWidth(2);
 
473
  tracking_piZero_reconCount_all->SetLineWidth(2);
 
474
  tracking_piZero_truthCount_all->SetLineColor(kBlack);
 
475
  tracking_piZero_reconCount_all->SetLineColor(kRed);
 
476
 
 
477
  tracking_piZero_truthCount_all->DrawCopy();
 
478
  tracking_piZero_reconCount_all->DrawCopy("SAME");
 
479
 
 
480
  c9->cd(4);
 
481
  gPad->SetGridx();
 
482
  gPad->SetGridy();
 
483
  gPad->SetTickx();
 
484
  gPad->SetTicky();
 
485
  gPad->Update();  
 
486
  TGraphAsymmErrors * piZero_eff_grp = new TGraphAsymmErrors(tracking_piZero_reconCount_all,tracking_piZero_truthCount_all);
 
487
  piZero_eff_grp->SetLineColor(kBlack);
 
488
  piZero_eff_grp->SetLineWidth(1);
 
489
  piZero_eff_grp->Draw("AP");
420
490
}
421
491
 
422
492
void setupHistTH1(TH1 * hist,Color_t colour){