Merge pull request #21 from DUNE/feature_FDVariablePlotting

Re-introduce executable functionality

.gitignore

@@ -2,6 +2,9 @@ build/
 *.root
 *.tbz2
 *.tgz
+*.pdf
+*.png
+*~
 logs/
 jupyter_env
 .history

configs/EventRates_Atmospherics.yaml

@@ -6,8 +6,11 @@ General:
   #Nu-FIT
   #OscillationParameters: [0.310, 0.582, 0.224, 7.39E-5, 2.5254E-3, -2.498, 25]
 
+  #Nu-FIT 5.2 w. SK atm
+  OscillationParameters: [0.303, 0.452, 0.0223, 7.41E-5, 2.51E-3, -2.233, 25]
+
   #T2K-like best-fit
-  OscillationParameters: [0.307, 0.528, 0.0218, 7.53e-5, 2.509e-3, -1.601, 25]
+  #OscillationParameters: [0.307, 0.528, 0.0218, 7.53e-5, 2.509e-3, -1.601, 25]
 
   #Unoscillated
   #OscillationParameters: [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 2.5]
@@ -34,3 +37,54 @@ General:
   ProcessMCMC: No
   Seed: 0
   Debug: No
+
+"Projections": [
+  {
+    "Name": "RecoNuEnergy", 
+    "VarString": "RecoNeutrinoEnergy",
+    "VarBins": [0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.2,1.4],
+    "KinematicCuts": [
+      {
+        "Range": [0.1,0.8],
+        "Name": "TrueNuEnergy", 
+        "VarString": "TrueNeutrinoEnergy"
+      }, 
+      {
+        "Range": [0.4,1.0],
+        "Name": "TrueCosZ", 
+        "VarString": "TrueCosineZ"
+      }
+    ], 
+    "CategoryCuts": [
+      {
+        "Name": "OscillationChannel_Single", 
+        "VarString": "OscChannel",
+        "Breakdown": [[0.0],  [1.0],  [2.0],  [3.0],  [4.0],  [5.0],  [6.0],  [7.0],  [8.0],  [9.0], [10.0], [11.0]],
+	"Names":     ["nue_x_nue","nue_x_numu","nue_x_nutau","numu_x_nue","numu_x_numu","numu_x_nutau","nuebar_x_nuebar","nuebar_x_numubar","nuebar_x_nutaubar","numubar_x_nuebar","numubar_x_numubar","numubar_x_nutaubar"],
+      },
+      {
+        "Name": "OscillationChannel_Group", 
+        "VarString": "OscChannel",
+        "Breakdown": [[0.0, 1.0, 2.0, 3.0, 4.0, 5.0],  [6.0, 7.0, 8.0, 9.0, 10.0, 11.0]],
+	"Names":     ["Nu","Nubar"],
+      },
+      {
+        "Name": "Mode_Single", 
+        "VarString": "Mode",
+        "Breakdown": [[0.0] , [1.0], [2.0], [3.0], [4.0], [5.0], [6.0], [7.0], [8.0], [9.0], [10.0], [11.0], [12.0], [13.0], [14.0], [15.0], [16.0], [17.0], [18.0], [19.0], [20.0], [21.0], [22.0], [23.0], [24.0], [25.0], [26.0]],
+      },
+      {
+        "Name": "Mode_Group", 
+        "VarString": "Mode",
+        "Breakdown": [ [0.0],   [2.0],   [3.0],   [9.0],  [15.0],  [16.0], [1.0, 4.0, 5.0, 6.0, 7.0, 8.0, 10.0, 11.0, 12.0, 13.0], [14.0, 17.0, 18.0, 19.0, 20.0, 21.0, 22.0, 23.0, 24.0, 25.0, 26.0]],
+	"Names":     ["CCQE", "CCDIS", "CCRES", "CCMEC", "NCDIS", "NCRES",                                                "CCOth",                                                            "NCOth"],
+      }
+    ],
+  },
+
+  {
+    "Name": "TrueNuEnergy",
+    "VarString": "TrueNeutrinoEnergy",
+    "VarBins": [0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.2,1.4,1.6,1.8,2.0,2.2,2.4,2.6,2.8,3.0,3.5,4.0,4.5,5.0,6.0,7.0,8.0,9.0,10.0],
+  }
+]

configs/Samples/AtmSample_AllMC.yaml

@@ -0,0 +1,108 @@
+SampleName: "numuselec"
+SampleType: "Atm"
+# SelectionCuts:
+#   - KinematicStr: "TrueXPos"
+#     Bounds: [ -310.0, 310.0 ]
+#   - KinematicStr: "TrueYPos"
+#     Bounds: [ -550.0, 550.0 ]
+#   - KinematicStr: "TrueZPos"
+#     Bounds: [ 50.0, 1244.0 ]
+Binning:
+  XVarStr: "TrueNeutrinoEnergy"
+  XVarBins: [0., 0.25, 0.5, 0.75, 1.,  1.25, 1.5, 1.75, 2., 2.25, 2.5, 2.75, 3., 3.25, 3.5, 3.75, 4., 4.25, 4.5, 4.75, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 100000.]
+SampleBools:
+  IsELike: no
+InputFiles:
+  mtupleprefix: "inputs/Atmospherics/CAFs/atm_hd"
+  mtuplesuffix: ".root"
+  splineprefix: ""
+  splinesuffix: ""
+DetID: 24
+NuOsc:
+  NuOscConfigFile: "configs/NuOsc/CUDAProb3.yaml"
+NSubSamples: 12
+SubSamples:
+  - name: "nue_x_nue"
+    mtuplefile: "nue_x_nue"
+    splinefile: "nue_x_nue"
+    samplevecno: 0
+    nutype: 12
+    oscnutype: 12
+    signal: false
+  - name: "nue_x_numu"
+    mtuplefile: "nue_x_numu"
+    splinefile: "nue_x_numu"
+    samplevecno: 1
+    nutype: 12
+    oscnutype: 14
+    signal: true
+  - name: "nue_x_nutau"
+    mtuplefile: "nue_x_nutau"
+    splinefile: "nue_x_nutau"
+    samplevecno: 2
+    nutype: 12
+    oscnutype: 16
+    signal: true
+  - name: "numu_x_nue"
+    mtuplefile: "numu_x_nue"
+    splinefile: "numu_x_nue"
+    samplevecno: 3
+    nutype: 14
+    oscnutype: 12
+    signal: true
+  - name: "numu_x_numu"
+    mtuplefile: "numu_x_numu"
+    splinefile: "numu_x_numu"
+    samplevecno: 4
+    nutype: 14
+    oscnutype: 14
+    signal: false
+  - name: "numu_x_nutau"
+    mtuplefile: "numu_x_nutau"
+    splinefile: "numu_x_nutau"
+    samplevecno: 5
+    nutype: 14
+    oscnutype: 16
+    signal: true
+  - name: "nuebar_x_nuebar"
+    mtuplefile: "nuebar_x_nuebar"
+    splinefile: "nuebar_x_nuebar"
+    samplevecno: 6
+    nutype: -12
+    oscnutype: -12
+    signal: false
+  - name: "nuebar_x_numubar"
+    mtuplefile: "nuebar_x_numubar"
+    splinefile: "nuebar_x_numubar"
+    samplevecno: 7
+    nutype: -12
+    oscnutype: -14
+    signal: true
+  - name: "nuebar_x_nutaubar"
+    mtuplefile: "nuebar_x_nutaubar"
+    splinefile: "nuebar_x_nutaubar"
+    samplevecno: 8
+    nutype: -12
+    oscnutype: -16
+    signal: true
+  - name: "numubar_x_nuebar"
+    mtuplefile: "numubar_x_nuebar"
+    splinefile: "numubar_x_nuebar"
+    samplevecno: 9
+    nutype: -14
+    oscnutype: -12
+    signal: true
+  - name: "numubar_x_numubar"
+    mtuplefile: "numubar_x_numubar"
+    splinefile: "numubar_x_numubar"
+    samplevecno: 10
+    nutype: -14
+    oscnutype: -14
+    signal: false
+  - name: "numubar_x_nutaubar"
+    mtuplefile: "numubar_x_nutaubar"
+    splinefile: "numubar_x_nutaubar"
+    samplevecno: 11
+    nutype: -14
+    oscnutype: -16
+    signal: true

configs/Samples/AtmSample_ncselec.yaml

@@ -0,0 +1,108 @@
+SampleName: "numuselec"
+SampleType: "Atm"
+# SelectionCuts:
+#   - KinematicStr: "TrueXPos"
+#     Bounds: [ -310.0, 310.0 ]
+#   - KinematicStr: "TrueYPos"
+#     Bounds: [ -550.0, 550.0 ]
+#   - KinematicStr: "TrueZPos"
+#     Bounds: [ 50.0, 1244.0 ]
+Binning:
+  XVarStr: "TrueNeutrinoEnergy"
+  XVarBins: [0., 0.25, 0.5, 0.75, 1.,  1.25, 1.5, 1.75, 2., 2.25, 2.5, 2.75, 3., 3.25, 3.5, 3.75, 4., 4.25, 4.5, 4.75, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 100000.]
+SampleBools:
+  IsELike: no
+InputFiles:
+  mtupleprefix: "inputs/Atmospherics/CAFs/atm_hd"
+  mtuplesuffix: "_ncselec.root"
+  splineprefix: ""
+  splinesuffix: ""
+DetID: 24
+NuOsc:
+  NuOscConfigFile: "configs/NuOsc/CUDAProb3.yaml"
+NSubSamples: 12
+SubSamples:
+  - name: "nue_x_nue"
+    mtuplefile: "nue_x_nue"
+    splinefile: "nue_x_nue"
+    samplevecno: 0
+    nutype: 12
+    oscnutype: 12
+    signal: false
+  - name: "nue_x_numu"
+    mtuplefile: "nue_x_numu"
+    splinefile: "nue_x_numu"
+    samplevecno: 1
+    nutype: 12
+    oscnutype: 14
+    signal: true
+  - name: "nue_x_nutau"
+    mtuplefile: "nue_x_nutau"
+    splinefile: "nue_x_nutau"
+    samplevecno: 2
+    nutype: 12
+    oscnutype: 16
+    signal: true
+  - name: "numu_x_nue"
+    mtuplefile: "numu_x_nue"
+    splinefile: "numu_x_nue"
+    samplevecno: 3
+    nutype: 14
+    oscnutype: 12
+    signal: true
+  - name: "numu_x_numu"
+    mtuplefile: "numu_x_numu"
+    splinefile: "numu_x_numu"
+    samplevecno: 4
+    nutype: 14
+    oscnutype: 14
+    signal: false
+  - name: "numu_x_nutau"
+    mtuplefile: "numu_x_nutau"
+    splinefile: "numu_x_nutau"
+    samplevecno: 5
+    nutype: 14
+    oscnutype: 16
+    signal: true
+  - name: "nuebar_x_nuebar"
+    mtuplefile: "nuebar_x_nuebar"
+    splinefile: "nuebar_x_nuebar"
+    samplevecno: 6
+    nutype: -12
+    oscnutype: -12
+    signal: false
+  - name: "nuebar_x_numubar"
+    mtuplefile: "nuebar_x_numubar"
+    splinefile: "nuebar_x_numubar"
+    samplevecno: 7
+    nutype: -12
+    oscnutype: -14
+    signal: true
+  - name: "nuebar_x_nutaubar"
+    mtuplefile: "nuebar_x_nutaubar"
+    splinefile: "nuebar_x_nutaubar"
+    samplevecno: 8
+    nutype: -12
+    oscnutype: -16
+    signal: true
+  - name: "numubar_x_nuebar"
+    mtuplefile: "numubar_x_nuebar"
+    splinefile: "numubar_x_nuebar"
+    samplevecno: 9
+    nutype: -14
+    oscnutype: -12
+    signal: true
+  - name: "numubar_x_numubar"
+    mtuplefile: "numubar_x_numubar"
+    splinefile: "numubar_x_numubar"
+    samplevecno: 10
+    nutype: -14
+    oscnutype: -14
+    signal: false
+  - name: "numubar_x_nutaubar"
+    mtuplefile: "numubar_x_nutaubar"
+    splinefile: "numubar_x_nutaubar"
+    samplevecno: 11
+    nutype: -14
+    oscnutype: -16
+    signal: true