Analyze_MCTruth.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Analyze MCTruth Example\n",
    "\n",
    "OK, now that we have that out of the way, let's do a little work with art-ROOT files via gallery in a ROOT C++-kernel notebook.\n",
    "\n",
    "Remember, `gallery` is the set of utilities/functions that allow us to have a sort of \"easy access\" to art-ROOT files, which are the standard output format for our simulation and reconstruction. You can read more about gallery, and see some additional examples from the [art group's webpage](http://art.fnal.gov/gallery/).\n",
    "\n",
    "---\n",
    "## Finding data objects\n",
    "\n",
    "This is not a full art/gallery demonstration, so I'm going to assume that you are reasonably aware of how to use and navigate the data format. One of the most useful things though is to know where to see the data objects. In general, when you setup a ups product, the source code is available to peruse in the `$THATPRODUCT_DIR/source` area. So, important areas are:\n",
    "* `$LARDATAOBJ_DIR/source/lardataobj`\n",
    "* `$NUSIMDATA_DIR/source/nusimdata`\n",
    "\n",
    "The contents of those areas are here:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[?1034hlardataobj/AnalysisBase/BackTrackerMatchingData.h\n",
      "lardataobj/AnalysisBase/Calorimetry.h\n",
      "lardataobj/AnalysisBase/CosmicTag.h\n",
      "lardataobj/AnalysisBase/FlashMatch.h\n",
      "lardataobj/AnalysisBase/MVAOutput.h\n",
      "lardataobj/AnalysisBase/MVAPIDResult.h\n",
      "lardataobj/AnalysisBase/ParticleID.h\n",
      "lardataobj/AnalysisBase/T0.h\n",
      "lardataobj/AnalysisBase/classes.h\n",
      "lardataobj/MCBase/MCBaseException.h\n",
      "lardataobj/MCBase/MCDataHolder.h\n",
      "lardataobj/MCBase/MCHit.h\n",
      "lardataobj/MCBase/MCHitCollection.h\n",
      "lardataobj/MCBase/MCLimits.h\n",
      "lardataobj/MCBase/MCShower.h\n",
      "lardataobj/MCBase/MCStep.h\n",
      "lardataobj/MCBase/MCTrack.h\n",
      "lardataobj/MCBase/MCWire.h\n",
      "lardataobj/MCBase/MCWireCollection.h\n",
      "lardataobj/MCBase/classes.h\n",
      "lardataobj/OpticalDetectorData/ChannelData.h\n",
      "lardataobj/OpticalDetectorData/ChannelDataGroup.h\n",
      "lardataobj/OpticalDetectorData/FIFOChannel.h\n",
      "lardataobj/OpticalDetectorData/OpticalRawDigit.h\n",
      "lardataobj/OpticalDetectorData/OpticalTypes.h\n",
      "lardataobj/OpticalDetectorData/PMTTrigger.h\n",
      "lardataobj/OpticalDetectorData/classes.h\n",
      "lardataobj/RawData/AuxDetDigit.h\n",
      "lardataobj/RawData/BeamInfo.h\n",
      "lardataobj/RawData/DAQHeader.h\n",
      "lardataobj/RawData/ExternalTrigger.h\n",
      "lardataobj/RawData/OpDetPulse.h\n",
      "lardataobj/RawData/OpDetWaveform.h\n",
      "lardataobj/RawData/RawDigit.h\n",
      "lardataobj/RawData/TriggerData.h\n",
      "lardataobj/RawData/classes.h\n",
      "lardataobj/RawData/raw.h\n",
      "lardataobj/RecoBase/Cluster.h\n",
      "lardataobj/RecoBase/Edge.h\n",
      "lardataobj/RecoBase/EndPoint2D.h\n",
      "lardataobj/RecoBase/Event.h\n",
      "lardataobj/RecoBase/Hit.h\n",
      "lardataobj/RecoBase/MCSFitResult.h\n",
      "lardataobj/RecoBase/OpFlash.h\n",
      "lardataobj/RecoBase/OpHit.h\n",
      "lardataobj/RecoBase/PCAxis.h\n",
      "lardataobj/RecoBase/PFParticle.h\n",
      "lardataobj/RecoBase/PointCharge.h\n",
      "lardataobj/RecoBase/Seed.h\n",
      "lardataobj/RecoBase/Shower.h\n",
      "lardataobj/RecoBase/SpacePoint.h\n",
      "lardataobj/RecoBase/Track.h\n",
      "lardataobj/RecoBase/TrackFitHitInfo.h\n",
      "lardataobj/RecoBase/TrackHitMeta.h\n",
      "lardataobj/RecoBase/TrackTrajectory.h\n",
      "lardataobj/RecoBase/TrackingPlane.h\n",
      "lardataobj/RecoBase/TrackingTypes.h\n",
      "lardataobj/RecoBase/Trajectory.h\n",
      "lardataobj/RecoBase/TrajectoryPointFlags.h\n",
      "lardataobj/RecoBase/Vertex.h\n",
      "lardataobj/RecoBase/VertexAssnMeta.h\n",
      "lardataobj/RecoBase/Wire.h\n",
      "lardataobj/RecoBase/classes.h\n",
      "lardataobj/Simulation/AuxDetSimChannel.h\n",
      "lardataobj/Simulation/BeamGateInfo.h\n",
      "lardataobj/Simulation/BeamTypes.h\n",
      "lardataobj/Simulation/GeneratedParticleInfo.h\n",
      "lardataobj/Simulation/OpDetBacktrackerRecord.h\n",
      "lardataobj/Simulation/SimChannel.h\n",
      "lardataobj/Simulation/SimEnergyDeposit.h\n",
      "lardataobj/Simulation/SimPhotons.h\n",
      "lardataobj/Simulation/SupernovaTruth.h\n",
      "lardataobj/Simulation/classes.h\n",
      "lardataobj/Simulation/sim.h\n",
      "lardataobj/Utilities/BitMask.h\n",
      "lardataobj/Utilities/DataIOmanip.h\n",
      "lardataobj/Utilities/FlagSet.h\n",
      "lardataobj/Utilities/SumSecondFunction.h\n",
      "lardataobj/Utilities/VectorMap.h\n",
      "lardataobj/Utilities/sparse_vector.h\n"
     ]
    }
   ],
   "source": [
    ".! (cd $LARDATAOBJ_DIR/source/ && ls lardataobj/*/*.h)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "nusimdata/SimulationBase/GTruth.h\n",
      "nusimdata/SimulationBase/MCFlux.h\n",
      "nusimdata/SimulationBase/MCNeutrino.h\n",
      "nusimdata/SimulationBase/MCParticle.h\n",
      "nusimdata/SimulationBase/MCTrajectory.h\n",
      "nusimdata/SimulationBase/MCTruth.h\n",
      "nusimdata/SimulationBase/classes.h\n",
      "nusimdata/SimulationBase/simb.h\n"
     ]
    }
   ],
   "source": [
    ".! (cd $NUSIMDATA_DIR/source/ && ls nusimdata/*/*.h)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Necessary includes and linking libraries <a href=\"includes\"></a>\n",
    "\n",
    "There are a number of neccesary `include` statements to use gallery, and you'll also need to include headers for the objects you're going to want to analyze from the event. You can do a simple include statement directly in a cell in the macro:\n",
    "\n",
    "```\n",
    "//gallery includes\n",
    "#include \"canvas/Utilities/InputTag.h\"\n",
    "#include \"gallery/Event.h\"\n",
    "#include \"gallery/ValidHandle.h\"\n",
    "#include \"canvas/Persistency/Common/FindMany.h\"\n",
    "#include \"canvas/Persistency/Common/FindOne.h\"\n",
    "#include \"canvas/Persistency/Common/FindManyP.h\"\n",
    "#include \"canvas/Persistency/Common/FindOneP.h\"\n",
    "```\n",
    "\n",
    "That may often work, but I sometimes have trouble, especially when it comes to getting proper templates known by ROOT for associations. Instead, I find doing things in the following way rather useful.\n",
    "\n",
    "First, I setup a directory to contain a lot of \"common\" use files. These are collections of headers and maybe other functions/utilities I think I might use a lot. You can see that directory <a href=\"common/\">here</a>. For this example in particular, we're going to load in these two files:\n",
    "* <a href=\"common/gallery_includes.h\">common/gallery_includes.h</a>\n",
    "* <a href=\"common/larsoftobj_includes.h\">common/larsoftobj_includes.h</a>\n",
    "\n",
    "Now, when we load these in, we're going to tell ROOT to compile them (with an optimization flag!). ROOT is smart enough to know if there's been a change since the last compilation, so it will only do it as necessary. We use the `ProcessLine()` functionality that you can use just the same from the ROOT command-line:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "gROOT->ProcessLine(\".L common/gallery_includes.h+O\");\n",
    "gROOT->ProcessLine(\".L common/larsoftobj_includes.h+O\");"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "And, just like that, we should be good to go on to writing code!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Declarations\n",
    "\n",
    "OK, remeber what we talked about in the <a href=\"ROOT_CPP_Kernel_Example.ipynb\">C++-kernel tutorial</a>, where it makes sense to declare the variables we're going to be using in cells ahead of where we use them. Let's go ahead and do that here too.\n",
    "\n",
    "Note: this is a process of trial and error, of course. You're not going to perfectly set everything up in one go before you even write the bulk of your code doing the analysis on each event. Just remember a few things:\n",
    "- You can \"Restart & Run All\" for the kernel\n",
    "- You can execute single cells to declare something new, and then re-execute a cell where you've made a change using the thing in that cell.\n",
    "- You don't have to execute cells in order. And you can move cells around as you wish.\n",
    "\n",
    "That said, let's talk about what we'll do here, and then what we want to set up. Let's say we want to loop over the MCTruth objects in a file and store the neutrino and lepton (if charged current) kinematics in a tree, and then write that tree to a file, and maybe make a plot or two in this notebook too.\n",
    "\n",
    "So, let's list out what we'll need:\n",
    "- a vector of strings to hold an input file list;\n",
    "- an `art::InputTag` to specify the label for the MCTruth object in the event;\n",
    "- the output tree;\n",
    "- an output file; and,\n",
    "- a canvas.\n",
    "\n",
    "Let's even think ahead slightly and give ourselves a \"verbose\" flag that we can turn on and off for testing/debugging.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "// string to contain our input files\n",
    "std::vector<std::string> filenames;\n",
    "\n",
    "//input tag for the generator\n",
    "art::InputTag gen_tag;\n",
    "\n",
    "//placeholder for the output file\n",
    "TFile* f_output;\n",
    "\n",
    "//a canvas in case we want one\n",
    "TCanvas canvas(\"canvas\",\"MyCanvas\");\n",
    "\n",
    "//setup the ROOT output tree\n",
    "unsigned int run, event;\n",
    "int nu_pdg,ccnc,mode;\n",
    "float q2,nu_energy;\n",
    "float vtx_x,vtx_y,vtx_z;\n",
    "float lep_energy,lep_p,lep_px,lep_py,lep_pz;\n",
    "  \n",
    "TTree* fAnaTree = new TTree(\"ana\",\"Nu Interaction Tree\");\n",
    "fAnaTree->Branch(\"run\",&run,\"run/I\");\n",
    "fAnaTree->Branch(\"event\",&event,\"event/I\");\n",
    "fAnaTree->Branch(\"nu_pdg\",&nu_pdg,\"nu_pdg/I\");\n",
    "fAnaTree->Branch(\"ccnc\",&ccnc,\"ccnc/I\");\n",
    "fAnaTree->Branch(\"mode\",&mode,\"mode/I\");\n",
    "fAnaTree->Branch(\"q2\",&q2,\"q2/F\");\n",
    "fAnaTree->Branch(\"nu_energy\",&nu_energy,\"nu_energy/F\");\n",
    "fAnaTree->Branch(\"vtx_x\",&vtx_x,\"vtx_x/F\");\n",
    "fAnaTree->Branch(\"vtx_y\",&vtx_y,\"vtx_y/F\");\n",
    "fAnaTree->Branch(\"vtx_z\",&vtx_z,\"vtx_z/F\");\n",
    "fAnaTree->Branch(\"lep_energy\",&lep_energy,\"lep_energy/F\");\n",
    "fAnaTree->Branch(\"lep_p\",&lep_p,\"lep_p/F\");\n",
    "fAnaTree->Branch(\"lep_px\",&lep_px,\"lep_px/F\");\n",
    "fAnaTree->Branch(\"lep_py\",&lep_py,\"lep_py/F\");\n",
    "fAnaTree->Branch(\"lep_pz\",&lep_pz,\"lep_pz/F\");\n",
    "\n",
    "//verbosity mode flag\n",
    "bool verbose=false;"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Initializations\n",
    "\n",
    "Right, now that everything has been declared, we actually want to initialize a lot of it. Note: this initialization cell should be able to be run repeatedly now, given that we're not re-declaring anythng, but just re-assigning things.\n",
    "\n",
    "One quick trick: while you might imagine we don't need to really think about the output file until later, it's helpful to setup the output file early, as it gives ROOT a place to write things into temporarily. This is especially important if you run over a ton of events."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "//setup our output file here to avoid write buffer problems\n",
    "//and, you know, make sure you change this to something in your area :)\n",
    "\n",
    "f_output = new TFile(\"/uboone/data/users/wketchum/tmp.root\",\"RECREATE\");\n",
    "\n",
    "//setup our input file and input tag\n",
    "//the nice thing about doing this on the gpvms is having easy access to all your data areas!\n",
    "filenames = \n",
    "{\n",
    "  \"/uboone/data/users/wketchum/prodgenie_bnb_intrinsic_nue_uboone_0_20170413T202620_gen2_419cbfe3-cbf6-4a31-808b-ef9fbbae0254.root\"\n",
    "};\n",
    "\n",
    "gen_tag = { \"generator\" };"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## The gallery loop\n",
    "\n",
    "With that all set, we are ready to run the gallery loop! Again, not gonna give you the full gallery demo here, but you can use this as a working example of how to do it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Successfully opened file /uboone/data/users/wketchum/prodgenie_bnb_intrinsic_nue_uboone_0_20170413T202620_gen2_419cbfe3-cbf6-4a31-808b-ef9fbbae0254.root\n",
      "Processing Run 2, Event 469951\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.26009 0 1 149.549,20.2299,342.826\n",
      "Processing Run 2, Event 469952\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.72622 0 2 126.533,3.13313,903.273\n",
      "Processing Run 2, Event 469953\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.315355 0 0 103.413,63.6562,791.466\n",
      "Processing Run 2, Event 469954\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.726046 0 0 216.097,-88.0515,262.828\n",
      "Processing Run 2, Event 469955\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 2.51577 1 0 141.224,-13.5385,987.196\n",
      "Processing Run 2, Event 469956\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.35134 0 1 206.98,94.4059,961.351\n",
      "Processing Run 2, Event 469957\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.39069 0 0 97.6515,-21.6837,705.237\n",
      "Processing Run 2, Event 469958\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.33284 0 1 43.987,-80.697,937.657\n",
      "Processing Run 2, Event 469959\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.57557 1 0 16.956,99.1211,403.037\n",
      "Processing Run 2, Event 469960\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.776472 0 0 190.723,110.549,680.566\n",
      "Processing Run 2, Event 469961\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.471156 1 0 235.539,58.9542,929.069\n",
      "Processing Run 2, Event 469962\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.482173 0 0 150.609,70.1337,555.627\n",
      "Processing Run 2, Event 469963\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.763834 0 0 249.847,96.8296,987.532\n",
      "Processing Run 2, Event 469964\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 3.86244 1 1 179.203,75.5334,28.1037\n",
      "Processing Run 2, Event 469965\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 2.59224 1 2 115.695,97.0794,348.071\n",
      "Processing Run 2, Event 469966\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.867861 0 0 10.7453,-55.3415,340.216\n",
      "Processing Run 2, Event 469967\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.36114 0 10 89.0203,96.2009,213.446\n",
      "Processing Run 2, Event 469968\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.350947 1 0 116.139,83.9783,157.789\n",
      "Processing Run 2, Event 469969\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.809077 0 10 131.86,114.968,545.546\n",
      "Processing Run 2, Event 469970\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t-12 1.78078 0 0 193.039,40.7338,962.797\n",
      "Processing Run 2, Event 469971\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 2.66725 1 2 115.645,-52.4288,698.821\n",
      "Processing Run 2, Event 469972\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.278989 0 0 104.759,56.6213,320.405\n",
      "Processing Run 2, Event 469973\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 2.39887 0 0 59.1386,-110.702,937.961\n",
      "Processing Run 2, Event 469974\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.949425 0 10 170.904,-14.2987,368.879\n",
      "Processing Run 2, Event 469975\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 3.16613 0 2 218.316,-32.6422,402.549\n",
      "Processing Run 2, Event 469976\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.547508 0 10 134.663,82.3767,143.614\n",
      "Processing Run 2, Event 469977\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.494289 0 0 232.375,71.7754,52.9169\n",
      "Processing Run 2, Event 469978\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.37765 0 0 132.125,105.885,937.581\n",
      "Processing Run 2, Event 469979\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.58385 0 10 76.1297,33.9576,283.741\n",
      "Processing Run 2, Event 469980\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.779145 1 1 170.396,-56.1199,943.821\n",
      "Processing Run 2, Event 469981\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.833849 0 0 249.518,24.6997,875.202\n",
      "Processing Run 2, Event 469982\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.02608 0 1 236.472,111.256,236.356\n",
      "Processing Run 2, Event 469983\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.66148 1 1 143.648,-16.0451,426.845\n",
      "Processing Run 2, Event 469984\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.66905 0 10 237.967,-72.3359,991.578\n",
      "Processing Run 2, Event 469985\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.634481 0 1 25.3528,29.4411,216.24\n",
      "Processing Run 2, Event 469986\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.38608 0 0 187.214,99.6613,475.951\n",
      "Processing Run 2, Event 469987\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.26555 0 1 28.0521,-48.4065,21.1801\n",
      "Processing Run 2, Event 469988\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.721126 0 1 203.494,19.9135,12.922\n",
      "Processing Run 2, Event 469989\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.564142 0 0 98.1473,-29.0521,452.257\n",
      "Processing Run 2, Event 469990\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.608834 0 10 61.3799,45.5627,8.87405\n",
      "Processing Run 2, Event 469991\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.751741 1 0 29.7029,84.3073,777.526\n",
      "Processing Run 2, Event 469992\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.01475 1 1 177.77,-63.7864,616.986\n",
      "Processing Run 2, Event 469993\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.1531 0 1 12.4536,53.9129,735.855\n",
      "Processing Run 2, Event 469994\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.08346 0 1 165.232,-41.0941,383.754\n",
      "Processing Run 2, Event 469995\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.992315 0 0 149.194,-6.00802,76.9998\n",
      "Processing Run 2, Event 469996\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 0.950547 0 0 250.856,-76.4272,542.312\n",
      "Processing Run 2, Event 469997\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t-12 2.02016 0 10 15.4649,70.3546,387.135\n",
      "Processing Run 2, Event 469998\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 1.08528 0 1 202.687,-75.6907,750.093\n",
      "Processing Run 2, Event 469999\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 3.3566 0 2 144.42,-14.1759,181.669\n",
      "Processing Run 2, Event 470000\n",
      "\tThere are 1 MCTruths in this event.\n",
      "\t\t12 2.72869 0 0 129.93,-82.1686,226.123\n",
      "Closing file, read 596107 bytes in 366 transactions\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "%MSG-d ProductID: \n",
      "Product created with id: [3130632973] from canonical product name: \"simb::MCTruths_generator__GenieGen.\"\n",
      "%MSG\n",
      "%MSG-d ProductID: \n",
      "Product created with id: [4283826371] from canonical product name: \"simb::MCTruths_generator__G4.\"\n",
      "%MSG\n",
      "%MSG-d ProductID: \n",
      "Product created with id: [852458230] from canonical product name: \"simb::MCTruths_generator__Detsim.\"\n",
      "%MSG\n"
     ]
    }
   ],
   "source": [
    "verbose = true;\n",
    "fAnaTree->Reset();\n",
    "for(gallery::Event ev(filenames) ; !ev.atEnd(); ev.next()) {\n",
    "\n",
    "    if(verbose){\n",
    "        std::cout << \"Processing \" \n",
    "                  << \"Run \" << ev.eventAuxiliary().run() << \", \"\n",
    "                  << \"Event \" << ev.eventAuxiliary().event() << endl;\n",
    "    }\n",
    "    \n",
    "    auto const& mctruth_handle = ev.getValidHandle<vector<simb::MCTruth>>(gen_tag);\n",
    "    auto const& truth_vec(*mctruth_handle);\n",
    "    if(verbose) std::cout << \"\\tThere are \" << truth_vec.size() << \" MCTruths in this event.\" << std::endl;\n",
    "\n",
    "    //Loop over and print them out\n",
    "    for( simb::MCTruth const& mct : truth_vec){\n",
    "        \n",
    "        //we only want to look at neutrinos...\n",
    "        if(mct.Origin()!=simb::Origin_t::kBeamNeutrino) continue;\n",
    "      \n",
    "        //get the neutrino and start filling our output info\n",
    "        auto const& nu = mct.GetNeutrino();\n",
    "        ccnc = nu.CCNC();\n",
    "        mode = nu.Mode();\n",
    "        q2 = nu.QSqr();\n",
    "        nu_energy = nu.Nu().E();\n",
    "        nu_pdg = nu.Nu().PdgCode();\n",
    "        vtx_x = nu.Nu().Vx();\n",
    "        vtx_y = nu.Nu().Vy();\n",
    "        vtx_z = nu.Nu().Vz();\n",
    "      \n",
    "      if(verbose) \n",
    "          std::cout << \"\\t\\t\" << nu_pdg << \" \" << nu_energy << \" \" << ccnc << \" \" << mode << \" \"\n",
    "                    << nu.Nu().Vx() << \",\" << nu.Nu().Vy() << \",\" << nu.Nu().Vz() << std::endl;\n",
    "      \n",
    "      if(ccnc==0) { //CC\n",
    "          lep_energy = nu.Lepton().E();\n",
    "          lep_p = nu.Lepton().P();\n",
    "          lep_px = nu.Lepton().Px();\n",
    "          lep_py = nu.Lepton().Py();\n",
    "          lep_pz = nu.Lepton().Pz();\n",
    "      }\n",
    "      else{\n",
    "          lep_energy = -999;\n",
    "          lep_p = -999;\n",
    "          lep_px = -999;\n",
    "          lep_py = -999;\n",
    "          lep_pz = -999;\n",
    "      }\n",
    "      fAnaTree->Fill();\n",
    "\n",
    "    }//end loop over truth object\n",
    "\n",
    "}//end loop over events\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Analyzing and storing the output\n",
    "\n",
    "Since we've already set everything up, we're ready to just look at it all in this notebook if we'd like! Let's make a few quick plots, like ... CC vs. NC? And neutrino energy vs. lepton energy in CC interactions?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "<div id=\"root_plot_6\"\n",
       "     style=\"width: 696px; height: 472px\">\n",
       "</div>\n",
       "\n",
       "<script>\n",
       " requirejs.config({\n",
       "     paths: {\n",
       "       'JSRootCore' : 'https://root.cern.ch/js/notebook//scripts/JSRootCore',\n",
       "     }\n",
       "   });\n",
       " require(['JSRootCore'],\n",
       "     function(Core) {\n",
       "       var obj = Core.JSONR_unref({\"_typename\":\"TCanvas\",\"fUniqueID\":0,\"fBits\":36896776,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fLeftMargin\":0.1,\"fRightMargin\":0.1,\"fBottomMargin\":0.1,\"fTopMargin\":0.1,\"fXfile\":2,\"fYfile\":2,\"fAfile\":1,\"fXstat\":0.99,\"fYstat\":0.99,\"fAstat\":2,\"fFrameFillColor\":0,\"fFrameLineColor\":1,\"fFrameFillStyle\":1001,\"fFrameLineStyle\":1,\"fFrameLineWidth\":1,\"fFrameBorderSize\":1,\"fFrameBorderMode\":0,\"fX1\":-0.250000018626453,\"fY1\":-5.11875038137662,\"fX2\":2.25000001862645,\"fY2\":46.0687503813766,\"fXtoAbsPixelk\":69.6000541484835,\"fXtoPixelk\":69.6000541484835,\"fXtoPixel\":278.399995851517,\"fYtoAbsPixelk\":424.800047186661,\"fYtoPixelk\":424.800047186661,\"fYtoPixel\":-9.22100108359759,\"fUtoAbsPixelk\":5e-5,\"fUtoPixelk\":5e-5,\"fUtoPixel\":696,\"fVtoAbsPixelk\":472.00005,\"fVtoPixelk\":472,\"fVtoPixel\":-472,\"fAbsPixeltoXk\":-0.250000018626453,\"fPixeltoXk\":-0.250000018626453,\"fPixeltoX\":0.00359195407651279,\"fAbsPixeltoYk\":46.0687503813766,\"fPixeltoYk\":-5.11875038137662,\"fPixeltoY\":-0.108448094836342,\"fXlowNDC\":0,\"fYlowNDC\":0,\"fXUpNDC\":0,\"fYUpNDC\":0,\"fWNDC\":1,\"fHNDC\":1,\"fAbsXlowNDC\":0,\"fAbsYlowNDC\":0,\"fAbsWNDC\":1,\"fAbsHNDC\":1,\"fUxmin\":0,\"fUymin\":0,\"fUxmax\":2,\"fUymax\":40.95,\"fTheta\":30,\"fPhi\":30,\"fAspectRatio\":0,\"fNumber\":0,\"fTickx\":0,\"fTicky\":0,\"fLogx\":0,\"fLogy\":0,\"fLogz\":0,\"fPadPaint\":0,\"fCrosshair\":0,\"fCrosshairPos\":0,\"fBorderSize\":2,\"fBorderMode\":0,\"fModified\":false,\"fGridx\":false,\"fGridy\":false,\"fAbsCoord\":false,\"fEditable\":true,\"fFixedAspectRatio\":false,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TFrame\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":0,\"fY1\":0,\"fX2\":2,\"fY2\":40.95,\"fBorderSize\":1,\"fBorderMode\":0},{\"_typename\":\"TH1F\",\"fUniqueID\":0,\"fBits\":50331656,\"fName\":\"h\",\"fTitle\":\"ccnc\",\"fLineColor\":602,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fMarkerColor\":1,\"fMarkerStyle\":1,\"fMarkerSize\":1,\"fNcells\":4,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":54525952,\"fName\":\"xaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":2,\"fXmin\":0,\"fXmax\":2,\"fXbins\":[],\"fFirst\":1,\"fLast\":100,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":50,\"fTsumw\":50,\"fTsumw2\":50,\"fTsumwx\":11,\"fTsumwx2\":11,\"fMaximum\":-1111,\"fMinimum\":0,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TPaveStats\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":1.70000005066395,\"fY1\":34.5515630148584,\"fX2\":2.20000006556511,\"fY2\":42.7415629538382,\"fX1NDC\":0.780000016093254,\"fY1NDC\":0.775000005960464,\"fX2NDC\":0.980000019073486,\"fY2NDC\":0.935000002384186,\"fBorderSize\":1,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"brNDC\",\"fName\":\"stats\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":12,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":18,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"h\",\"fTextAngle\":0,\"fTextSize\":0.0368,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0368000008165836},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Entries = 50     \",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Mean  =   0.22\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04},{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"Std Dev   = 0.4142\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.04}],\"opt\":[\"\",\"\",\"\",\"\"]},\"fOptFit\":0,\"fOptStat\":1111,\"fFitFormat\":\"5.4g\",\"fStatFormat\":\"6.4g\",\"fParent\":{\"$ref\":3}}],\"opt\":[\"brNDC\"]},\"fBufferSize\":2001,\"fBuffer\":[-50,1,0,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,1,1,0,1,0,1,1,1,1,1,0,1,0,1,1,1,0,1,0,1,1,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,1,1,0,1,0,1,1,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,1,1,1,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],\"fBinStatErrOpt\":0,\"fArray\":[0,39,11,0]},{\"_typename\":\"TPaveText\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fX1\":0.878017239561625,\"fY1\":42.997500537741,\"fX2\":1.12198276043837,\"fY2\":45.8128131216439,\"fX1NDC\":0.451206896551724,\"fY1NDC\":0.940000003948808,\"fX2NDC\":0.548793103448276,\"fY2NDC\":0.995000004768372,\"fBorderSize\":0,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"blNDC\",\"fName\":\"title\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":22,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":4,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"ccnc\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0467500016093254}],\"opt\":[\"\"]}}],\"opt\":[\"\",\"\",\"blNDC\"]},\"fExecs\":null,\"fName\":\"canvas\",\"fTitle\":\"MyCanvas\",\"fNumPaletteColor\":0,\"fNextPaletteColor\":0,\"fDISPLAY\":\"$DISPLAY\",\"fDoubleBuffer\":0,\"fRetained\":true,\"fXsizeUser\":0,\"fYsizeUser\":0,\"fXsizeReal\":20,\"fYsizeReal\":14.28571,\"fWindowTopX\":0,\"fWindowTopY\":0,\"fWindowWidth\":0,\"fWindowHeight\":0,\"fCw\":696,\"fCh\":472,\"fCatt\":{\"_typename\":\"TAttCanvas\",\"fXBetween\":2,\"fYBetween\":2,\"fTitleFromTop\":1.2,\"fXdate\":0.2,\"fYdate\":0.3,\"fAdate\":1},\"kMoveOpaque\":true,\"kResizeOpaque\":true,\"fHighLightColor\":2,\"fBatch\":true,\"kShowEventStatus\":false,\"kAutoExec\":true,\"kMenuBar\":true});\n",
       "       Core.draw(\"root_plot_6\", obj, \"\");\n",
       "     }\n",
       " );\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%jsroot on\n",
    "canvas.cd();\n",
    "fAnaTree->Draw(\"ccnc>>h\");\n",
    "h->SetMinimum(0);\n",
    "canvas.Draw();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "<div id=\"root_plot_10\"\n",
       "     style=\"width: 696px; height: 472px\">\n",
       "</div>\n",
       "\n",
       "<script>\n",
       " requirejs.config({\n",
       "     paths: {\n",
       "       'JSRootCore' : 'https://root.cern.ch/js/notebook//scripts/JSRootCore',\n",
       "     }\n",
       "   });\n",
       " require(['JSRootCore'],\n",
       "     function(Core) {\n",
       "       var obj = Core.JSONR_unref({\"_typename\":\"TCanvas\",\"fUniqueID\":0,\"fBits\":36896776,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fLeftMargin\":0.1,\"fRightMargin\":0.1,\"fBottomMargin\":0.1,\"fTopMargin\":0.1,\"fXfile\":2,\"fYfile\":2,\"fAfile\":1,\"fXstat\":0.99,\"fYstat\":0.99,\"fAstat\":2,\"fFrameFillColor\":0,\"fFrameLineColor\":1,\"fFrameFillStyle\":1001,\"fFrameLineStyle\":1,\"fFrameLineWidth\":1,\"fFrameBorderSize\":1,\"fFrameBorderMode\":0,\"fX1\":-0.350000026077034,\"fY1\":-0.450000033527615,\"fX2\":3.15000002607703,\"fY2\":4.05000003352761,\"fXtoAbsPixelk\":69.6000541484835,\"fXtoPixelk\":69.6000541484835,\"fXtoPixel\":198.85713989394,\"fYtoAbsPixelk\":424.800047186661,\"fYtoPixelk\":424.800047186661,\"fYtoPixel\":-104.888887325923,\"fUtoAbsPixelk\":5e-5,\"fUtoPixelk\":5e-5,\"fUtoPixel\":696,\"fVtoAbsPixelk\":472.00005,\"fVtoPixelk\":472,\"fVtoPixel\":-472,\"fAbsPixeltoXk\":-0.350000026077034,\"fPixeltoXk\":-0.350000026077034,\"fPixeltoX\":0.00502873570711791,\"fAbsPixeltoYk\":4.05000003352761,\"fPixeltoYk\":-0.450000033527615,\"fPixeltoY\":-0.00953389844715091,\"fXlowNDC\":0,\"fYlowNDC\":0,\"fXUpNDC\":0,\"fYUpNDC\":0,\"fWNDC\":1,\"fHNDC\":1,\"fAbsXlowNDC\":0,\"fAbsYlowNDC\":0,\"fAbsWNDC\":1,\"fAbsHNDC\":1,\"fUxmin\":0,\"fUymin\":0,\"fUxmax\":2.8,\"fUymax\":3.6,\"fTheta\":30,\"fPhi\":30,\"fAspectRatio\":0,\"fNumber\":0,\"fTickx\":0,\"fTicky\":0,\"fLogx\":0,\"fLogy\":0,\"fLogz\":0,\"fPadPaint\":0,\"fCrosshair\":0,\"fCrosshairPos\":0,\"fBorderSize\":2,\"fBorderMode\":0,\"fModified\":false,\"fGridx\":false,\"fGridy\":false,\"fAbsCoord\":false,\"fEditable\":true,\"fFixedAspectRatio\":false,\"fPrimitives\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TFrame\",\"fUniqueID\":0,\"fBits\":50331656,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fX1\":0,\"fY1\":0,\"fX2\":2.8,\"fY2\":3.6,\"fBorderSize\":1,\"fBorderMode\":0},{\"_typename\":\"TH2F\",\"fUniqueID\":0,\"fBits\":50332169,\"fName\":\"htemp\",\"fTitle\":\"nu_energy:lep_energy {ccnc==0}\",\"fLineColor\":602,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fMarkerColor\":1,\"fMarkerStyle\":8,\"fMarkerSize\":1,\"fNcells\":1764,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"xaxis\",\"fTitle\":\"lep_energy\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":40,\"fXmin\":0,\"fXmax\":2.8,\"fXbins\":[],\"fFirst\":1,\"fLast\":40,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\"nu_energy\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":40,\"fXmin\":0,\"fXmax\":3.6,\"fXbins\":[],\"fFirst\":1,\"fLast\":40,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":0,\"fTsumw\":0,\"fTsumw2\":0,\"fTsumwx\":0,\"fTsumwx2\":0,\"fMaximum\":-1111,\"fMinimum\":-1111,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[],\"opt\":[]},\"fBufferSize\":3001,\"fBuffer\":[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],\"fBinStatErrOpt\":0,\"fScalefactor\":1,\"fTsumwy\":0,\"fTsumwy2\":0,\"fTsumwxy\":0,\"fArray\":[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]},{\"_typename\":\"TPaveText\",\"fUniqueID\":0,\"fBits\":50331657,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":0,\"fX1\":0.417083318686733,\"fY1\":3.75292377937206,\"fX2\":2.38291668131327,\"fY2\":4.02750005465001,\"fX1NDC\":0.219166666666667,\"fY1NDC\":0.933983055615829,\"fX2NDC\":0.780833333333333,\"fY2NDC\":0.995000004768372,\"fBorderSize\":0,\"fInit\":1,\"fShadowColor\":1,\"fCornerRadius\":0,\"fOption\":\"blNDC\",\"fName\":\"title\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":22,\"fTextColor\":1,\"fTextFont\":42,\"fLabel\":\"\",\"fLongest\":30,\"fMargin\":0.05,\"fLines\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[{\"_typename\":\"TLatex\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"\",\"fTitle\":\"nu_energy:lep_energy {ccnc==0}\",\"fTextAngle\":0,\"fTextSize\":0,\"fTextAlign\":0,\"fTextColor\":0,\"fTextFont\":0,\"fX\":0,\"fY\":0,\"fLineColor\":1,\"fLineStyle\":1,\"fLineWidth\":2,\"fLimitFactorSize\":3,\"fOriginSize\":0.0518644079566002}],\"opt\":[\"\"]}},{\"_typename\":\"TGraph\",\"fUniqueID\":0,\"fBits\":50594825,\"fName\":\"Graph\",\"fTitle\":\"Graph\",\"fLineColor\":602,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fMarkerColor\":1,\"fMarkerStyle\":8,\"fMarkerSize\":1,\"fNpoints\":39,\"fX\":[0.323625952005386,0.305804342031479,0.265897274017334,0.638446152210236,0.996000587940216,1.01262819766998,0.718929827213287,0.382989883422852,0.41228511929512,0.503746569156647,0.657756328582764,0.203627526760101,0.48620343208313,1.65346252918243,0.229688510298729,1.94391238689423,0.66587096452713,0.86928391456604,0.188297688961029,0.366372317075729,0.58515477180481,1.30306327342987,0.365816295146942,0.324208468198776,1.38854098320007,0.193710207939148,0.309802561998367,0.591222763061523,0.170312538743019,0.223700255155563,0.20206144452095,0.437045484781265,0.317095160484314,0.9421306848526,0.604887247085571,1.54373598098755,0.103969186544418,0.844232082366943,2.62651371955872],\"fY\":[1.26008629798889,1.72622060775757,0.315354973077774,0.726045966148376,1.35133504867554,1.39069426059723,1.332843542099,0.776472091674805,0.482172608375549,0.76383364200592,0.867861390113831,0.361139982938766,0.809077322483063,1.78077590465546,0.278988569974899,2.3988676071167,0.949425280094147,3.16613125801086,0.547507703304291,0.494289249181747,1.37764847278595,1.58384537696838,0.833849310874939,1.02608048915863,1.66905403137207,0.634480714797974,0.386080145835876,1.2655531167984,0.72112649679184,0.564142107963562,0.608834087848663,1.15310299396515,1.08345818519592,0.99231493473053,0.950547158718109,2.02015590667725,1.08528196811676,3.35659956932068,2.72869229316711],\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[],\"opt\":[]},\"fHistogram\":{\"_typename\":\"TH1F\",\"fUniqueID\":0,\"fBits\":50332168,\"fName\":\"Graph\",\"fTitle\":\"Graph\",\"fLineColor\":602,\"fLineStyle\":1,\"fLineWidth\":1,\"fFillColor\":0,\"fFillStyle\":1001,\"fMarkerColor\":1,\"fMarkerStyle\":1,\"fMarkerSize\":1,\"fNcells\":102,\"fXaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"xaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":100,\"fXmin\":0,\"fXmax\":2.87876817286015,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fYaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"yaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":0,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0.251089712977409,\"fXmax\":3.66436066925526,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fZaxis\":{\"_typename\":\"TAxis\",\"fUniqueID\":0,\"fBits\":50331648,\"fName\":\"zaxis\",\"fTitle\":\"\",\"fNdivisions\":510,\"fAxisColor\":1,\"fLabelColor\":1,\"fLabelFont\":42,\"fLabelOffset\":0.005,\"fLabelSize\":0.035,\"fTickLength\":0.03,\"fTitleOffset\":1,\"fTitleSize\":0.035,\"fTitleColor\":1,\"fTitleFont\":42,\"fNbins\":1,\"fXmin\":0,\"fXmax\":1,\"fXbins\":[],\"fFirst\":0,\"fLast\":0,\"fBits2\":0,\"fTimeDisplay\":false,\"fTimeFormat\":\"\",\"fLabels\":null,\"fModLabs\":null},\"fBarOffset\":0,\"fBarWidth\":1000,\"fEntries\":0,\"fTsumw\":0,\"fTsumw2\":0,\"fTsumwx\":0,\"fTsumwx2\":0,\"fMaximum\":3.66436066925526,\"fMinimum\":0.251089712977409,\"fNormFactor\":0,\"fContour\":[],\"fSumw2\":[],\"fOption\":\"\",\"fFunctions\":{\"_typename\":\"TList\",\"name\":\"TList\",\"arr\":[],\"opt\":[]},\"fBufferSize\":0,\"fBuffer\":[],\"fBinStatErrOpt\":0,\"fArray\":[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]},\"fMinimum\":-1111,\"fMaximum\":-1111}],\"opt\":[\"\",\"\",\"blNDC\",\"p\"]},\"fExecs\":null,\"fName\":\"canvas\",\"fTitle\":\"MyCanvas\",\"fNumPaletteColor\":0,\"fNextPaletteColor\":0,\"fDISPLAY\":\"$DISPLAY\",\"fDoubleBuffer\":0,\"fRetained\":true,\"fXsizeUser\":0,\"fYsizeUser\":0,\"fXsizeReal\":20,\"fYsizeReal\":14.28571,\"fWindowTopX\":0,\"fWindowTopY\":0,\"fWindowWidth\":0,\"fWindowHeight\":0,\"fCw\":696,\"fCh\":472,\"fCatt\":{\"_typename\":\"TAttCanvas\",\"fXBetween\":2,\"fYBetween\":2,\"fTitleFromTop\":1.2,\"fXdate\":0.2,\"fYdate\":0.3,\"fAdate\":1},\"kMoveOpaque\":true,\"kResizeOpaque\":true,\"fHighLightColor\":2,\"fBatch\":true,\"kShowEventStatus\":false,\"kAutoExec\":true,\"kMenuBar\":true});\n",
       "       Core.draw(\"root_plot_10\", obj, \"\");\n",
       "     }\n",
       " );\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fAnaTree->SetMarkerStyle(8);\n",
    "fAnaTree->Draw(\"nu_energy:lep_energy\",\"ccnc==0\");\n",
    "canvas.Draw();"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "OK, and then, so we can save our work and analyze it in a notebook later, let's save the tree to our output file."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "f_output->cd();\n",
    "fAnaTree->Write();\n",
    "f_output->Close();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ROOT C++",
   "language": "c++",
   "name": "root"
  },
  "language_info": {
   "codemirror_mode": "text/x-c++src",
   "file_extension": ".C",
   "mimetype": " text/x-c++src",
   "name": "c++"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}