HOWTO simulate and analyze b1pi events
From GlueXWiki
Contents
Introduction
This page gives a very short tutorial on how one might simulate and analyze a specific physics channel using the Hall-D source code. The channel represented here is the so-called "b1π" channel in which a 2GeV/c2 meson is produced which subsequently decays into a b1 and π. (Note: The plots on this page were made using rev. 5911 of the Hall-D source code with JANA 0.6.0.)
Event Generation
We generate the events with the genr8 program. This uses an isobar model with an exponential t distribution. It is likely not exactly how the final distributions will look, but it does give something with a final state that is complex and should roughly overlap the true distributions. There already exists a configuration file for the b1π reaction for genr8 in the following location within the HAll-D source code tree:
trunk/sim-recon/src/programs/Simulation/genr8/InputFiles/b1_pi.input
For this exercise, I edited the resonance widths to be artificially narrow so as to get cleaner signals in the end which indicate directly the reconstructed resolutions. In particular, I used:
| Meson | Γ |
|---|---|
| X(2000) | 1 MeV |
| b1(1235) | 1 MeV |
| ω(783) | 1 MeV |
| ρ(776) | 35 MeV |
| πo | 0 MeV |
(n.b. the ρ width is larger because genr8 seems to take much longer to run if it is made narrower.)
Run the genr8 program to create an ASCII file of events and then convert the ASCII file to HDDM with the genr8_2_hddm program. Finally, you'll want to rename the default name of output.hddm of the output file to match the other base filenames.
genr8 -Ab1_pi.ascii -M1000 < b1_pi.input genr8_2_hddm b1_pi.ascii mv output.hddm b1_pi.hddm
Simulation
You'll need to edit a copy of a control.in file to use to configure hdgeant to run. A minimal file would be the following:
INFILE 'b1_pi.hddm' TRIG 1000 OUTFILE 'hdgeant.hddm' RNDM 121 c The following line controls the cutoffs for tracking of particles. c CUTS cutgam cutele cutneu cuthad cutmuo bcute bcutm dcute dcutm ppcutm tofmax CUTS 1e-4 1e-4 1e-3 1e-3 1e-4 c The following card enables the GelHad package (from BaBar) c on/off ecut scale mode thresh GELH 1 0.2 1.0 4 0.160 DCAY 0 no decay in flight END
Not all detector resolution is included in the hdgeant program so the output must be further processed with the mcsmear program. In fact, calorimetry reconstruction will not work at all without this. Run the simulation and smear (digitize) the output:
hdgeant mcsmear -n hdgeant.hddm
We use the "-n" option with mcsmear to prevent it from adding random noise hits. This should result in a file named "hdgeant_smeared.hddm" which has all of the detector hits with expected resolution. Note that if you used the control.in file above then no extra beam background was included in the simulation.
Reconstruction (DEPRECATED)
For reconstruction we use the hd_root program. This is a program that has all reconstruction modules already compiled in. The program is only useful if plugins are attached that will activate the reconstruction and produce ROOT objects (trees or histograms) from the reconstructed values. For this exercise we use the somewhat generic phys_tree plugin. This, like many plugins, are not compiled as part of the default Hall-D build so you may have to go in and explicitly build it by invoking "make" in the src/programs/Analysis/plugins/phys_tree directory.
Run hd_root with the phys_tree plugin on your hdgeant_smeared.hddm file with the following:
hd_root -PPLUGINS=phys_tree hdgeant_smeared.hddm
This will produce a file called hd_root.root that has a TDirectory in it that contains a couple of trees. There are actually 2 trees (thrown and recon) plus a third that simply adds the other 2 as friends. The thrown and recon trees have exactly the same format, but the former holds generated values and the latter reconstructed ones.
Analysis (DEPRECATED)
One can do simple plotting within ROOT of the reconstructed values. For example, to plot the reconstructed π+ momentum distribution, one can do as illustrated in the the following session:
>root hd_root.root
*******************************************
* *
* W E L C O M E to R O O T *
* *
* Version 5.22/00 17 December 2008 *
* *
* You are welcome to visit our Web site *
* http://root.cern.ch *
* *
*******************************************
ROOT 5.22/00 (trunk@26997, De 18 2008, 10:17:00 on macosx)
CINT/ROOT C/C++ Interpreter version 5.16.29, Jan 08, 2008
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.
Executing ~/.rootlogon.C
root [0]
Attaching file hd_root.root as _file0...
Warning in <TClass::TClass>: no dictionary for class Event is available
Warning in <TClass::TClass>: no dictionary for class Particle is available
root [1] PHYSICS->cd()
(Bool_t)1
root [2] recon->Draw("R.pip.p.P()", "R.pip.p.P()<12")
<TCanvas::MakeDefCanvas>: created default TCanvas with name c1
(Long64_t)2244
root [3]
|
The plot to the right shows the result. The meaning of the "R.pip.p.P()" is as follows:
| R : Branch name of Reconstructed tree |
| pip : List of π+ Particles. By not specifying an index (e.g. pip[2]) we are plotting all π+s |
| p : The 4-momentum of the particle as a TLorentzVector |
| P() : The "P()" method of TLorentzVector which just returns the magnitude of the 3-momentum |
Similarly, one can look at the multiplicity of say, reconstructed photons in the event with the Nphoton member:
root [3] recon->Draw("Nphoton")
|
Notice here that I didn't need to specify the branch name "R" explicitly since Nphoton is not a complex object.
TSelector
To do an actual analysis with the tree produced by phys_tree, one would need to use a TSelector object. An appropriate template can be created by invoking the MakeSelector("") method of the recon tree. Describing this in detail is beyond the scope of this HOWTO. However, I've copied the header and source files I used to make the plots below into the Appendix at the end of this document. To invoke the selector and make invariant mass plots, you need to copied the header and source files as "b1_pi.h" and "b1_pi.C" respectively into the working directory and then do the following:
>root hd_root.root
root [1] PHYSICS->cd()
root [2] recon->Process("b1_pi.C")
X(2000) meson
ρ(776) mass from π+π-
πo mass from γγ
ω(783) mass from ρπo
b1(1235) mass from ωπ+
The above plots contain no fiducial cuts or kinematic fitting. Nor do they include any information from reconstructed protons. They do, however indicate about 52% of the generated events are reconstructed fully. A more more detailed analysis would include beam backgrounds, in-flight decays, etc. This should at least get you started.
Please direct any questions or comments to David Lawrence
Appendix
- Note that the phys_tree plugin and the below TSelector are DEPRECATED.
Script to generate, simulate, obfuscate, and "recreate" b1 π events
#!/bin/tcsh -f
set NEVENTS=1000
set VERTEX="0 0 50 80" # x y zmin zmax
set NTHREADS=Ncores
set b1pi_input=${HALLD_HOME}/src/programs/Simulation/genr8/InputFiles/b1_pi.input
if ( ! -e ${b1pi_input} ) then
echo " "
echo "No file: ${b1pi_input}"
echo "Please make sure your HALLD_HOME env. var. is set."
echo " "
exit -1
endif
echo "Copying ${b1pi_input} ..."
cp $b1pi_input b1_pi.input
echo "Running genr8 ..."
genr8 -M${NEVENTS} -Ab1_pi.ascii < b1_pi.input >& /dev/null
echo "Converting generated events to HDDM ..."
genr8_2_hddm -V"${VERTEX}" b1_pi.ascii >& /dev/null
mv output.hddm b1_pi.hddm
echo "Creating control.in file ..."
cat - << EOF > control.in
INFILE 'b1_pi.hddm'
TRIG ${NEVENTS}
OUTFILE 'hdgeant.hddm'
RNDM 123
HADR 0
EOF
echo "Running hdgeant ..."
hdgeant
echo "Running mcsmear ..."
mcsmear hdgeant.hddm
echo "Running hd_root ..."
hd_root -PNTHREADS=${NTHREADS} -PPLUGINS=phys_tree hdgeant_smeared.hddm
b1_pi.h
//////////////////////////////////////////////////////////
// This class has been automatically generated on
// Mon May 9 09:29:11 2011 by ROOT version 5.26/00
// from TTree recon/Reconstructed Event parameters
// found on file: hd_root.root
//////////////////////////////////////////////////////////
#ifndef b1_pi_h
#define b1_pi_h
#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include <TSelector.h>
const Int_t kMaxpip = 4;
const Int_t kMaxpim = 3;
const Int_t kMaxKp = 1;
const Int_t kMaxKm = 1;
const Int_t kMaxproton = 3;
const Int_t kMaxphoton = 6;
const Int_t kMaxpip_match = 4;
const Int_t kMaxpim_match = 3;
const Int_t kMaxKp_match = 1;
const Int_t kMaxKm_match = 1;
const Int_t kMaxproton_match = 3;
const Int_t kMaxphoton_match = 6;
class b1_pi : public TSelector {
public :
TTree *fChain; //!pointer to the analyzed TTree or TChain
// Declaration of leaf types
//Event *R;
UInt_t fUniqueID;
UInt_t fBits;
Int_t event;
UInt_t Npip;
UInt_t Npim;
UInt_t NKp;
UInt_t NKm;
UInt_t Nproton;
UInt_t Nphoton;
Int_t pip_;
UInt_t pip_fUniqueID[kMaxpip]; //[pip_]
UInt_t pip_fBits[kMaxpip]; //[pip_]
UInt_t pip_p_fUniqueID[kMaxpip]; //[pip_]
UInt_t pip_p_fBits[kMaxpip]; //[pip_]
UInt_t pip_p_fP_fUniqueID[kMaxpip]; //[pip_]
UInt_t pip_p_fP_fBits[kMaxpip]; //[pip_]
Double_t pip_p_fP_fX[kMaxpip]; //[pip_]
Double_t pip_p_fP_fY[kMaxpip]; //[pip_]
Double_t pip_p_fP_fZ[kMaxpip]; //[pip_]
Double_t pip_p_fE[kMaxpip]; //[pip_]
UInt_t pip_x_fUniqueID[kMaxpip]; //[pip_]
UInt_t pip_x_fBits[kMaxpip]; //[pip_]
Double_t pip_x_fX[kMaxpip]; //[pip_]
Double_t pip_x_fY[kMaxpip]; //[pip_]
Double_t pip_x_fZ[kMaxpip]; //[pip_]
Bool_t pip_is_fiducial[kMaxpip]; //[pip_]
Double_t pip_chisq[kMaxpip]; //[pip_]
Int_t pip_Ndof[kMaxpip]; //[pip_]
Double_t pip_FOM_pid[kMaxpip]; //[pip_]
Int_t pim_;
UInt_t pim_fUniqueID[kMaxpim]; //[pim_]
UInt_t pim_fBits[kMaxpim]; //[pim_]
UInt_t pim_p_fUniqueID[kMaxpim]; //[pim_]
UInt_t pim_p_fBits[kMaxpim]; //[pim_]
UInt_t pim_p_fP_fUniqueID[kMaxpim]; //[pim_]
UInt_t pim_p_fP_fBits[kMaxpim]; //[pim_]
Double_t pim_p_fP_fX[kMaxpim]; //[pim_]
Double_t pim_p_fP_fY[kMaxpim]; //[pim_]
Double_t pim_p_fP_fZ[kMaxpim]; //[pim_]
Double_t pim_p_fE[kMaxpim]; //[pim_]
UInt_t pim_x_fUniqueID[kMaxpim]; //[pim_]
UInt_t pim_x_fBits[kMaxpim]; //[pim_]
Double_t pim_x_fX[kMaxpim]; //[pim_]
Double_t pim_x_fY[kMaxpim]; //[pim_]
Double_t pim_x_fZ[kMaxpim]; //[pim_]
Bool_t pim_is_fiducial[kMaxpim]; //[pim_]
Double_t pim_chisq[kMaxpim]; //[pim_]
Int_t pim_Ndof[kMaxpim]; //[pim_]
Double_t pim_FOM_pid[kMaxpim]; //[pim_]
Int_t Kp_;
UInt_t Kp_fUniqueID[kMaxKp]; //[Kp_]
UInt_t Kp_fBits[kMaxKp]; //[Kp_]
UInt_t Kp_p_fUniqueID[kMaxKp]; //[Kp_]
UInt_t Kp_p_fBits[kMaxKp]; //[Kp_]
UInt_t Kp_p_fP_fUniqueID[kMaxKp]; //[Kp_]
UInt_t Kp_p_fP_fBits[kMaxKp]; //[Kp_]
Double_t Kp_p_fP_fX[kMaxKp]; //[Kp_]
Double_t Kp_p_fP_fY[kMaxKp]; //[Kp_]
Double_t Kp_p_fP_fZ[kMaxKp]; //[Kp_]
Double_t Kp_p_fE[kMaxKp]; //[Kp_]
UInt_t Kp_x_fUniqueID[kMaxKp]; //[Kp_]
UInt_t Kp_x_fBits[kMaxKp]; //[Kp_]
Double_t Kp_x_fX[kMaxKp]; //[Kp_]
Double_t Kp_x_fY[kMaxKp]; //[Kp_]
Double_t Kp_x_fZ[kMaxKp]; //[Kp_]
Bool_t Kp_is_fiducial[kMaxKp]; //[Kp_]
Double_t Kp_chisq[kMaxKp]; //[Kp_]
Int_t Kp_Ndof[kMaxKp]; //[Kp_]
Double_t Kp_FOM_pid[kMaxKp]; //[Kp_]
Int_t Km_;
UInt_t Km_fUniqueID[kMaxKm]; //[Km_]
UInt_t Km_fBits[kMaxKm]; //[Km_]
UInt_t Km_p_fUniqueID[kMaxKm]; //[Km_]
UInt_t Km_p_fBits[kMaxKm]; //[Km_]
UInt_t Km_p_fP_fUniqueID[kMaxKm]; //[Km_]
UInt_t Km_p_fP_fBits[kMaxKm]; //[Km_]
Double_t Km_p_fP_fX[kMaxKm]; //[Km_]
Double_t Km_p_fP_fY[kMaxKm]; //[Km_]
Double_t Km_p_fP_fZ[kMaxKm]; //[Km_]
Double_t Km_p_fE[kMaxKm]; //[Km_]
UInt_t Km_x_fUniqueID[kMaxKm]; //[Km_]
UInt_t Km_x_fBits[kMaxKm]; //[Km_]
Double_t Km_x_fX[kMaxKm]; //[Km_]
Double_t Km_x_fY[kMaxKm]; //[Km_]
Double_t Km_x_fZ[kMaxKm]; //[Km_]
Bool_t Km_is_fiducial[kMaxKm]; //[Km_]
Double_t Km_chisq[kMaxKm]; //[Km_]
Int_t Km_Ndof[kMaxKm]; //[Km_]
Double_t Km_FOM_pid[kMaxKm]; //[Km_]
Int_t proton_;
UInt_t proton_fUniqueID[kMaxproton]; //[proton_]
UInt_t proton_fBits[kMaxproton]; //[proton_]
UInt_t proton_p_fUniqueID[kMaxproton]; //[proton_]
UInt_t proton_p_fBits[kMaxproton]; //[proton_]
UInt_t proton_p_fP_fUniqueID[kMaxproton]; //[proton_]
UInt_t proton_p_fP_fBits[kMaxproton]; //[proton_]
Double_t proton_p_fP_fX[kMaxproton]; //[proton_]
Double_t proton_p_fP_fY[kMaxproton]; //[proton_]
Double_t proton_p_fP_fZ[kMaxproton]; //[proton_]
Double_t proton_p_fE[kMaxproton]; //[proton_]
UInt_t proton_x_fUniqueID[kMaxproton]; //[proton_]
UInt_t proton_x_fBits[kMaxproton]; //[proton_]
Double_t proton_x_fX[kMaxproton]; //[proton_]
Double_t proton_x_fY[kMaxproton]; //[proton_]
Double_t proton_x_fZ[kMaxproton]; //[proton_]
Bool_t proton_is_fiducial[kMaxproton]; //[proton_]
Double_t proton_chisq[kMaxproton]; //[proton_]
Int_t proton_Ndof[kMaxproton]; //[proton_]
Double_t proton_FOM_pid[kMaxproton]; //[proton_]
Int_t photon_;
UInt_t photon_fUniqueID[kMaxphoton]; //[photon_]
UInt_t photon_fBits[kMaxphoton]; //[photon_]
UInt_t photon_p_fUniqueID[kMaxphoton]; //[photon_]
UInt_t photon_p_fBits[kMaxphoton]; //[photon_]
UInt_t photon_p_fP_fUniqueID[kMaxphoton]; //[photon_]
UInt_t photon_p_fP_fBits[kMaxphoton]; //[photon_]
Double_t photon_p_fP_fX[kMaxphoton]; //[photon_]
Double_t photon_p_fP_fY[kMaxphoton]; //[photon_]
Double_t photon_p_fP_fZ[kMaxphoton]; //[photon_]
Double_t photon_p_fE[kMaxphoton]; //[photon_]
UInt_t photon_x_fUniqueID[kMaxphoton]; //[photon_]
UInt_t photon_x_fBits[kMaxphoton]; //[photon_]
Double_t photon_x_fX[kMaxphoton]; //[photon_]
Double_t photon_x_fY[kMaxphoton]; //[photon_]
Double_t photon_x_fZ[kMaxphoton]; //[photon_]
Bool_t photon_is_fiducial[kMaxphoton]; //[photon_]
Double_t photon_chisq[kMaxphoton]; //[photon_]
Int_t photon_Ndof[kMaxphoton]; //[photon_]
Double_t photon_FOM_pid[kMaxphoton]; //[photon_]
Int_t pip_match_;
UInt_t pip_match_fUniqueID[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_fBits[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_p_fUniqueID[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_p_fBits[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_p_fP_fUniqueID[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_p_fP_fBits[kMaxpip_match]; //[pip_match_]
Double_t pip_match_p_fP_fX[kMaxpip_match]; //[pip_match_]
Double_t pip_match_p_fP_fY[kMaxpip_match]; //[pip_match_]
Double_t pip_match_p_fP_fZ[kMaxpip_match]; //[pip_match_]
Double_t pip_match_p_fE[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_x_fUniqueID[kMaxpip_match]; //[pip_match_]
UInt_t pip_match_x_fBits[kMaxpip_match]; //[pip_match_]
Double_t pip_match_x_fX[kMaxpip_match]; //[pip_match_]
Double_t pip_match_x_fY[kMaxpip_match]; //[pip_match_]
Double_t pip_match_x_fZ[kMaxpip_match]; //[pip_match_]
Bool_t pip_match_is_fiducial[kMaxpip_match]; //[pip_match_]
Double_t pip_match_chisq[kMaxpip_match]; //[pip_match_]
Int_t pip_match_Ndof[kMaxpip_match]; //[pip_match_]
Double_t pip_match_FOM_pid[kMaxpip_match]; //[pip_match_]
Int_t pim_match_;
UInt_t pim_match_fUniqueID[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_fBits[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_p_fUniqueID[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_p_fBits[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_p_fP_fUniqueID[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_p_fP_fBits[kMaxpim_match]; //[pim_match_]
Double_t pim_match_p_fP_fX[kMaxpim_match]; //[pim_match_]
Double_t pim_match_p_fP_fY[kMaxpim_match]; //[pim_match_]
Double_t pim_match_p_fP_fZ[kMaxpim_match]; //[pim_match_]
Double_t pim_match_p_fE[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_x_fUniqueID[kMaxpim_match]; //[pim_match_]
UInt_t pim_match_x_fBits[kMaxpim_match]; //[pim_match_]
Double_t pim_match_x_fX[kMaxpim_match]; //[pim_match_]
Double_t pim_match_x_fY[kMaxpim_match]; //[pim_match_]
Double_t pim_match_x_fZ[kMaxpim_match]; //[pim_match_]
Bool_t pim_match_is_fiducial[kMaxpim_match]; //[pim_match_]
Double_t pim_match_chisq[kMaxpim_match]; //[pim_match_]
Int_t pim_match_Ndof[kMaxpim_match]; //[pim_match_]
Double_t pim_match_FOM_pid[kMaxpim_match]; //[pim_match_]
Int_t Kp_match_;
UInt_t Kp_match_fUniqueID[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_fBits[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_p_fUniqueID[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_p_fBits[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_p_fP_fUniqueID[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_p_fP_fBits[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_p_fP_fX[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_p_fP_fY[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_p_fP_fZ[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_p_fE[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_x_fUniqueID[kMaxKp_match]; //[Kp_match_]
UInt_t Kp_match_x_fBits[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_x_fX[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_x_fY[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_x_fZ[kMaxKp_match]; //[Kp_match_]
Bool_t Kp_match_is_fiducial[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_chisq[kMaxKp_match]; //[Kp_match_]
Int_t Kp_match_Ndof[kMaxKp_match]; //[Kp_match_]
Double_t Kp_match_FOM_pid[kMaxKp_match]; //[Kp_match_]
Int_t Km_match_;
UInt_t Km_match_fUniqueID[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_fBits[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_p_fUniqueID[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_p_fBits[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_p_fP_fUniqueID[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_p_fP_fBits[kMaxKm_match]; //[Km_match_]
Double_t Km_match_p_fP_fX[kMaxKm_match]; //[Km_match_]
Double_t Km_match_p_fP_fY[kMaxKm_match]; //[Km_match_]
Double_t Km_match_p_fP_fZ[kMaxKm_match]; //[Km_match_]
Double_t Km_match_p_fE[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_x_fUniqueID[kMaxKm_match]; //[Km_match_]
UInt_t Km_match_x_fBits[kMaxKm_match]; //[Km_match_]
Double_t Km_match_x_fX[kMaxKm_match]; //[Km_match_]
Double_t Km_match_x_fY[kMaxKm_match]; //[Km_match_]
Double_t Km_match_x_fZ[kMaxKm_match]; //[Km_match_]
Bool_t Km_match_is_fiducial[kMaxKm_match]; //[Km_match_]
Double_t Km_match_chisq[kMaxKm_match]; //[Km_match_]
Int_t Km_match_Ndof[kMaxKm_match]; //[Km_match_]
Double_t Km_match_FOM_pid[kMaxKm_match]; //[Km_match_]
Int_t proton_match_;
UInt_t proton_match_fUniqueID[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_fBits[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_p_fUniqueID[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_p_fBits[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_p_fP_fUniqueID[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_p_fP_fBits[kMaxproton_match]; //[proton_match_]
Double_t proton_match_p_fP_fX[kMaxproton_match]; //[proton_match_]
Double_t proton_match_p_fP_fY[kMaxproton_match]; //[proton_match_]
Double_t proton_match_p_fP_fZ[kMaxproton_match]; //[proton_match_]
Double_t proton_match_p_fE[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_x_fUniqueID[kMaxproton_match]; //[proton_match_]
UInt_t proton_match_x_fBits[kMaxproton_match]; //[proton_match_]
Double_t proton_match_x_fX[kMaxproton_match]; //[proton_match_]
Double_t proton_match_x_fY[kMaxproton_match]; //[proton_match_]
Double_t proton_match_x_fZ[kMaxproton_match]; //[proton_match_]
Bool_t proton_match_is_fiducial[kMaxproton_match]; //[proton_match_]
Double_t proton_match_chisq[kMaxproton_match]; //[proton_match_]
Int_t proton_match_Ndof[kMaxproton_match]; //[proton_match_]
Double_t proton_match_FOM_pid[kMaxproton_match]; //[proton_match_]
Int_t photon_match_;
UInt_t photon_match_fUniqueID[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_fBits[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_p_fUniqueID[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_p_fBits[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_p_fP_fUniqueID[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_p_fP_fBits[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_p_fP_fX[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_p_fP_fY[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_p_fP_fZ[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_p_fE[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_x_fUniqueID[kMaxphoton_match]; //[photon_match_]
UInt_t photon_match_x_fBits[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_x_fX[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_x_fY[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_x_fZ[kMaxphoton_match]; //[photon_match_]
Bool_t photon_match_is_fiducial[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_chisq[kMaxphoton_match]; //[photon_match_]
Int_t photon_match_Ndof[kMaxphoton_match]; //[photon_match_]
Double_t photon_match_FOM_pid[kMaxphoton_match]; //[photon_match_]
UInt_t target_fUniqueID;
UInt_t target_fBits;
UInt_t target_fP_fUniqueID;
UInt_t target_fP_fBits;
Double_t target_fP_fX;
Double_t target_fP_fY;
Double_t target_fP_fZ;
Double_t target_fE;
UInt_t beam_fUniqueID;
UInt_t beam_fBits;
UInt_t beam_fP_fUniqueID;
UInt_t beam_fP_fBits;
Double_t beam_fP_fX;
Double_t beam_fP_fY;
Double_t beam_fP_fZ;
Double_t beam_fE;
UInt_t vertex_fUniqueID;
UInt_t vertex_fBits;
Double_t vertex_fX;
Double_t vertex_fY;
Double_t vertex_fZ;
UInt_t W_fUniqueID;
UInt_t W_fBits;
UInt_t W_fP_fUniqueID;
UInt_t W_fP_fBits;
Double_t W_fP_fX;
Double_t W_fP_fY;
Double_t W_fP_fZ;
Double_t W_fE;
Bool_t all_fiducial;
Bool_t all_mesons_fiducial;
Bool_t all_photons_fiducial;
Bool_t all_protons_fiducial;
// List of branches
TBranch *b_R_fUniqueID; //!
TBranch *b_R_fBits; //!
TBranch *b_R_event; //!
TBranch *b_R_Npip; //!
TBranch *b_R_Npim; //!
TBranch *b_R_NKp; //!
TBranch *b_R_NKm; //!
TBranch *b_R_Nproton; //!
TBranch *b_R_Nphoton; //!
TBranch *b_R_pip_; //!
TBranch *b_pip_fUniqueID; //!
TBranch *b_pip_fBits; //!
TBranch *b_pip_p_fUniqueID; //!
TBranch *b_pip_p_fBits; //!
TBranch *b_pip_p_fP_fUniqueID; //!
TBranch *b_pip_p_fP_fBits; //!
TBranch *b_pip_p_fP_fX; //!
TBranch *b_pip_p_fP_fY; //!
TBranch *b_pip_p_fP_fZ; //!
TBranch *b_pip_p_fE; //!
TBranch *b_pip_x_fUniqueID; //!
TBranch *b_pip_x_fBits; //!
TBranch *b_pip_x_fX; //!
TBranch *b_pip_x_fY; //!
TBranch *b_pip_x_fZ; //!
TBranch *b_pip_is_fiducial; //!
TBranch *b_pip_chisq; //!
TBranch *b_pip_Ndof; //!
TBranch *b_pip_FOM_pid; //!
TBranch *b_R_pim_; //!
TBranch *b_pim_fUniqueID; //!
TBranch *b_pim_fBits; //!
TBranch *b_pim_p_fUniqueID; //!
TBranch *b_pim_p_fBits; //!
TBranch *b_pim_p_fP_fUniqueID; //!
TBranch *b_pim_p_fP_fBits; //!
TBranch *b_pim_p_fP_fX; //!
TBranch *b_pim_p_fP_fY; //!
TBranch *b_pim_p_fP_fZ; //!
TBranch *b_pim_p_fE; //!
TBranch *b_pim_x_fUniqueID; //!
TBranch *b_pim_x_fBits; //!
TBranch *b_pim_x_fX; //!
TBranch *b_pim_x_fY; //!
TBranch *b_pim_x_fZ; //!
TBranch *b_pim_is_fiducial; //!
TBranch *b_pim_chisq; //!
TBranch *b_pim_Ndof; //!
TBranch *b_pim_FOM_pid; //!
TBranch *b_R_Kp_; //!
TBranch *b_Kp_fUniqueID; //!
TBranch *b_Kp_fBits; //!
TBranch *b_Kp_p_fUniqueID; //!
TBranch *b_Kp_p_fBits; //!
TBranch *b_Kp_p_fP_fUniqueID; //!
TBranch *b_Kp_p_fP_fBits; //!
TBranch *b_Kp_p_fP_fX; //!
TBranch *b_Kp_p_fP_fY; //!
TBranch *b_Kp_p_fP_fZ; //!
TBranch *b_Kp_p_fE; //!
TBranch *b_Kp_x_fUniqueID; //!
TBranch *b_Kp_x_fBits; //!
TBranch *b_Kp_x_fX; //!
TBranch *b_Kp_x_fY; //!
TBranch *b_Kp_x_fZ; //!
TBranch *b_Kp_is_fiducial; //!
TBranch *b_Kp_chisq; //!
TBranch *b_Kp_Ndof; //!
TBranch *b_Kp_FOM_pid; //!
TBranch *b_R_Km_; //!
TBranch *b_Km_fUniqueID; //!
TBranch *b_Km_fBits; //!
TBranch *b_Km_p_fUniqueID; //!
TBranch *b_Km_p_fBits; //!
TBranch *b_Km_p_fP_fUniqueID; //!
TBranch *b_Km_p_fP_fBits; //!
TBranch *b_Km_p_fP_fX; //!
TBranch *b_Km_p_fP_fY; //!
TBranch *b_Km_p_fP_fZ; //!
TBranch *b_Km_p_fE; //!
TBranch *b_Km_x_fUniqueID; //!
TBranch *b_Km_x_fBits; //!
TBranch *b_Km_x_fX; //!
TBranch *b_Km_x_fY; //!
TBranch *b_Km_x_fZ; //!
TBranch *b_Km_is_fiducial; //!
TBranch *b_Km_chisq; //!
TBranch *b_Km_Ndof; //!
TBranch *b_Km_FOM_pid; //!
TBranch *b_R_proton_; //!
TBranch *b_proton_fUniqueID; //!
TBranch *b_proton_fBits; //!
TBranch *b_proton_p_fUniqueID; //!
TBranch *b_proton_p_fBits; //!
TBranch *b_proton_p_fP_fUniqueID; //!
TBranch *b_proton_p_fP_fBits; //!
TBranch *b_proton_p_fP_fX; //!
TBranch *b_proton_p_fP_fY; //!
TBranch *b_proton_p_fP_fZ; //!
TBranch *b_proton_p_fE; //!
TBranch *b_proton_x_fUniqueID; //!
TBranch *b_proton_x_fBits; //!
TBranch *b_proton_x_fX; //!
TBranch *b_proton_x_fY; //!
TBranch *b_proton_x_fZ; //!
TBranch *b_proton_is_fiducial; //!
TBranch *b_proton_chisq; //!
TBranch *b_proton_Ndof; //!
TBranch *b_proton_FOM_pid; //!
TBranch *b_R_photon_; //!
TBranch *b_photon_fUniqueID; //!
TBranch *b_photon_fBits; //!
TBranch *b_photon_p_fUniqueID; //!
TBranch *b_photon_p_fBits; //!
TBranch *b_photon_p_fP_fUniqueID; //!
TBranch *b_photon_p_fP_fBits; //!
TBranch *b_photon_p_fP_fX; //!
TBranch *b_photon_p_fP_fY; //!
TBranch *b_photon_p_fP_fZ; //!
TBranch *b_photon_p_fE; //!
TBranch *b_photon_x_fUniqueID; //!
TBranch *b_photon_x_fBits; //!
TBranch *b_photon_x_fX; //!
TBranch *b_photon_x_fY; //!
TBranch *b_photon_x_fZ; //!
TBranch *b_photon_is_fiducial; //!
TBranch *b_photon_chisq; //!
TBranch *b_photon_Ndof; //!
TBranch *b_photon_FOM_pid; //!
TBranch *b_R_pip_match_; //!
TBranch *b_pip_match_fUniqueID; //!
TBranch *b_pip_match_fBits; //!
TBranch *b_pip_match_p_fUniqueID; //!
TBranch *b_pip_match_p_fBits; //!
TBranch *b_pip_match_p_fP_fUniqueID; //!
TBranch *b_pip_match_p_fP_fBits; //!
TBranch *b_pip_match_p_fP_fX; //!
TBranch *b_pip_match_p_fP_fY; //!
TBranch *b_pip_match_p_fP_fZ; //!
TBranch *b_pip_match_p_fE; //!
TBranch *b_pip_match_x_fUniqueID; //!
TBranch *b_pip_match_x_fBits; //!
TBranch *b_pip_match_x_fX; //!
TBranch *b_pip_match_x_fY; //!
TBranch *b_pip_match_x_fZ; //!
TBranch *b_pip_match_is_fiducial; //!
TBranch *b_pip_match_chisq; //!
TBranch *b_pip_match_Ndof; //!
TBranch *b_pip_match_FOM_pid; //!
TBranch *b_R_pim_match_; //!
TBranch *b_pim_match_fUniqueID; //!
TBranch *b_pim_match_fBits; //!
TBranch *b_pim_match_p_fUniqueID; //!
TBranch *b_pim_match_p_fBits; //!
TBranch *b_pim_match_p_fP_fUniqueID; //!
TBranch *b_pim_match_p_fP_fBits; //!
TBranch *b_pim_match_p_fP_fX; //!
TBranch *b_pim_match_p_fP_fY; //!
TBranch *b_pim_match_p_fP_fZ; //!
TBranch *b_pim_match_p_fE; //!
TBranch *b_pim_match_x_fUniqueID; //!
TBranch *b_pim_match_x_fBits; //!
TBranch *b_pim_match_x_fX; //!
TBranch *b_pim_match_x_fY; //!
TBranch *b_pim_match_x_fZ; //!
TBranch *b_pim_match_is_fiducial; //!
TBranch *b_pim_match_chisq; //!
TBranch *b_pim_match_Ndof; //!
TBranch *b_pim_match_FOM_pid; //!
TBranch *b_R_Kp_match_; //!
TBranch *b_Kp_match_fUniqueID; //!
TBranch *b_Kp_match_fBits; //!
TBranch *b_Kp_match_p_fUniqueID; //!
TBranch *b_Kp_match_p_fBits; //!
TBranch *b_Kp_match_p_fP_fUniqueID; //!
TBranch *b_Kp_match_p_fP_fBits; //!
TBranch *b_Kp_match_p_fP_fX; //!
TBranch *b_Kp_match_p_fP_fY; //!
TBranch *b_Kp_match_p_fP_fZ; //!
TBranch *b_Kp_match_p_fE; //!
TBranch *b_Kp_match_x_fUniqueID; //!
TBranch *b_Kp_match_x_fBits; //!
TBranch *b_Kp_match_x_fX; //!
TBranch *b_Kp_match_x_fY; //!
TBranch *b_Kp_match_x_fZ; //!
TBranch *b_Kp_match_is_fiducial; //!
TBranch *b_Kp_match_chisq; //!
TBranch *b_Kp_match_Ndof; //!
TBranch *b_Kp_match_FOM_pid; //!
TBranch *b_R_Km_match_; //!
TBranch *b_Km_match_fUniqueID; //!
TBranch *b_Km_match_fBits; //!
TBranch *b_Km_match_p_fUniqueID; //!
TBranch *b_Km_match_p_fBits; //!
TBranch *b_Km_match_p_fP_fUniqueID; //!
TBranch *b_Km_match_p_fP_fBits; //!
TBranch *b_Km_match_p_fP_fX; //!
TBranch *b_Km_match_p_fP_fY; //!
TBranch *b_Km_match_p_fP_fZ; //!
TBranch *b_Km_match_p_fE; //!
TBranch *b_Km_match_x_fUniqueID; //!
TBranch *b_Km_match_x_fBits; //!
TBranch *b_Km_match_x_fX; //!
TBranch *b_Km_match_x_fY; //!
TBranch *b_Km_match_x_fZ; //!
TBranch *b_Km_match_is_fiducial; //!
TBranch *b_Km_match_chisq; //!
TBranch *b_Km_match_Ndof; //!
TBranch *b_Km_match_FOM_pid; //!
TBranch *b_R_proton_match_; //!
TBranch *b_proton_match_fUniqueID; //!
TBranch *b_proton_match_fBits; //!
TBranch *b_proton_match_p_fUniqueID; //!
TBranch *b_proton_match_p_fBits; //!
TBranch *b_proton_match_p_fP_fUniqueID; //!
TBranch *b_proton_match_p_fP_fBits; //!
TBranch *b_proton_match_p_fP_fX; //!
TBranch *b_proton_match_p_fP_fY; //!
TBranch *b_proton_match_p_fP_fZ; //!
TBranch *b_proton_match_p_fE; //!
TBranch *b_proton_match_x_fUniqueID; //!
TBranch *b_proton_match_x_fBits; //!
TBranch *b_proton_match_x_fX; //!
TBranch *b_proton_match_x_fY; //!
TBranch *b_proton_match_x_fZ; //!
TBranch *b_proton_match_is_fiducial; //!
TBranch *b_proton_match_chisq; //!
TBranch *b_proton_match_Ndof; //!
TBranch *b_proton_match_FOM_pid; //!
TBranch *b_R_photon_match_; //!
TBranch *b_photon_match_fUniqueID; //!
TBranch *b_photon_match_fBits; //!
TBranch *b_photon_match_p_fUniqueID; //!
TBranch *b_photon_match_p_fBits; //!
TBranch *b_photon_match_p_fP_fUniqueID; //!
TBranch *b_photon_match_p_fP_fBits; //!
TBranch *b_photon_match_p_fP_fX; //!
TBranch *b_photon_match_p_fP_fY; //!
TBranch *b_photon_match_p_fP_fZ; //!
TBranch *b_photon_match_p_fE; //!
TBranch *b_photon_match_x_fUniqueID; //!
TBranch *b_photon_match_x_fBits; //!
TBranch *b_photon_match_x_fX; //!
TBranch *b_photon_match_x_fY; //!
TBranch *b_photon_match_x_fZ; //!
TBranch *b_photon_match_is_fiducial; //!
TBranch *b_photon_match_chisq; //!
TBranch *b_photon_match_Ndof; //!
TBranch *b_photon_match_FOM_pid; //!
TBranch *b_R_target_fUniqueID; //!
TBranch *b_R_target_fBits; //!
TBranch *b_R_target_fP_fUniqueID; //!
TBranch *b_R_target_fP_fBits; //!
TBranch *b_R_target_fP_fX; //!
TBranch *b_R_target_fP_fY; //!
TBranch *b_R_target_fP_fZ; //!
TBranch *b_R_target_fE; //!
TBranch *b_R_beam_fUniqueID; //!
TBranch *b_R_beam_fBits; //!
TBranch *b_R_beam_fP_fUniqueID; //!
TBranch *b_R_beam_fP_fBits; //!
TBranch *b_R_beam_fP_fX; //!
TBranch *b_R_beam_fP_fY; //!
TBranch *b_R_beam_fP_fZ; //!
TBranch *b_R_beam_fE; //!
TBranch *b_R_vertex_fUniqueID; //!
TBranch *b_R_vertex_fBits; //!
TBranch *b_R_vertex_fX; //!
TBranch *b_R_vertex_fY; //!
TBranch *b_R_vertex_fZ; //!
TBranch *b_R_W_fUniqueID; //!
TBranch *b_R_W_fBits; //!
TBranch *b_R_W_fP_fUniqueID; //!
TBranch *b_R_W_fP_fBits; //!
TBranch *b_R_W_fP_fX; //!
TBranch *b_R_W_fP_fY; //!
TBranch *b_R_W_fP_fZ; //!
TBranch *b_R_W_fE; //!
TBranch *b_R_all_fiducial; //!
TBranch *b_R_all_mesons_fiducial; //!
TBranch *b_R_all_photons_fiducial; //!
TBranch *b_R_all_protons_fiducial; //!
b1_pi(TTree * /*tree*/ =0) { }
virtual ~b1_pi() { }
virtual Int_t Version() const { return 2; }
virtual void Begin(TTree *tree);
virtual void SlaveBegin(TTree *tree);
virtual void Init(TTree *tree);
virtual Bool_t Notify();
virtual Bool_t Process(Long64_t entry);
virtual Int_t GetEntry(Long64_t entry, Int_t getall = 0) { return fChain ? fChain->GetTree()->GetEntry(entry, getall) : 0; }
virtual void SetOption(const char *option) { fOption = option; }
virtual void SetObject(TObject *obj) { fObject = obj; }
virtual void SetInputList(TList *input) { fInput = input; }
virtual TList *GetOutputList() const { return fOutput; }
virtual void SlaveTerminate();
virtual void Terminate();
//-------------------------------------
// Add pointers to my histograms as members of this class
TH1D *mass_X;
TH1D *mass_b1;
TH1D *mass_omega;
TH1D *mass_rho;
TH1D *mass_pi0;
//-------------------------------------
ClassDef(b1_pi,0);
};
#endif
#ifdef b1_pi_cxx
void b1_pi::Init(TTree *tree)
{
// The Init() function is called when the selector needs to initialize
// a new tree or chain. Typically here the branch addresses and branch
// pointers of the tree will be set.
// It is normally not necessary to make changes to the generated
// code, but the routine can be extended by the user if needed.
// Init() will be called many times when running on PROOF
// (once per file to be processed).
// Set branch addresses and branch pointers
if (!tree) return;
fChain = tree;
fChain->SetMakeClass(1);
fChain->SetBranchAddress("fUniqueID", &fUniqueID, &b_R_fUniqueID);
fChain->SetBranchAddress("fBits", &fBits, &b_R_fBits);
fChain->SetBranchAddress("event", &event, &b_R_event);
fChain->SetBranchAddress("Npip", &Npip, &b_R_Npip);
fChain->SetBranchAddress("Npim", &Npim, &b_R_Npim);
fChain->SetBranchAddress("NKp", &NKp, &b_R_NKp);
fChain->SetBranchAddress("NKm", &NKm, &b_R_NKm);
fChain->SetBranchAddress("Nproton", &Nproton, &b_R_Nproton);
fChain->SetBranchAddress("Nphoton", &Nphoton, &b_R_Nphoton);
fChain->SetBranchAddress("pip", &pip_, &b_R_pip_);
fChain->SetBranchAddress("pip.fUniqueID", pip_fUniqueID, &b_pip_fUniqueID);
fChain->SetBranchAddress("pip.fBits", pip_fBits, &b_pip_fBits);
fChain->SetBranchAddress("pip.p.fUniqueID", pip_p_fUniqueID, &b_pip_p_fUniqueID);
fChain->SetBranchAddress("pip.p.fBits", pip_p_fBits, &b_pip_p_fBits);
fChain->SetBranchAddress("pip.p.fP.fUniqueID", pip_p_fP_fUniqueID, &b_pip_p_fP_fUniqueID);
fChain->SetBranchAddress("pip.p.fP.fBits", pip_p_fP_fBits, &b_pip_p_fP_fBits);
fChain->SetBranchAddress("pip.p.fP.fX", pip_p_fP_fX, &b_pip_p_fP_fX);
fChain->SetBranchAddress("pip.p.fP.fY", pip_p_fP_fY, &b_pip_p_fP_fY);
fChain->SetBranchAddress("pip.p.fP.fZ", pip_p_fP_fZ, &b_pip_p_fP_fZ);
fChain->SetBranchAddress("pip.p.fE", pip_p_fE, &b_pip_p_fE);
fChain->SetBranchAddress("pip.x.fUniqueID", pip_x_fUniqueID, &b_pip_x_fUniqueID);
fChain->SetBranchAddress("pip.x.fBits", pip_x_fBits, &b_pip_x_fBits);
fChain->SetBranchAddress("pip.x.fX", pip_x_fX, &b_pip_x_fX);
fChain->SetBranchAddress("pip.x.fY", pip_x_fY, &b_pip_x_fY);
fChain->SetBranchAddress("pip.x.fZ", pip_x_fZ, &b_pip_x_fZ);
fChain->SetBranchAddress("pip.is_fiducial", pip_is_fiducial, &b_pip_is_fiducial);
fChain->SetBranchAddress("pip.chisq", pip_chisq, &b_pip_chisq);
fChain->SetBranchAddress("pip.Ndof", pip_Ndof, &b_pip_Ndof);
fChain->SetBranchAddress("pip.FOM_pid", pip_FOM_pid, &b_pip_FOM_pid);
fChain->SetBranchAddress("pim", &pim_, &b_R_pim_);
fChain->SetBranchAddress("pim.fUniqueID", pim_fUniqueID, &b_pim_fUniqueID);
fChain->SetBranchAddress("pim.fBits", pim_fBits, &b_pim_fBits);
fChain->SetBranchAddress("pim.p.fUniqueID", pim_p_fUniqueID, &b_pim_p_fUniqueID);
fChain->SetBranchAddress("pim.p.fBits", pim_p_fBits, &b_pim_p_fBits);
fChain->SetBranchAddress("pim.p.fP.fUniqueID", pim_p_fP_fUniqueID, &b_pim_p_fP_fUniqueID);
fChain->SetBranchAddress("pim.p.fP.fBits", pim_p_fP_fBits, &b_pim_p_fP_fBits);
fChain->SetBranchAddress("pim.p.fP.fX", pim_p_fP_fX, &b_pim_p_fP_fX);
fChain->SetBranchAddress("pim.p.fP.fY", pim_p_fP_fY, &b_pim_p_fP_fY);
fChain->SetBranchAddress("pim.p.fP.fZ", pim_p_fP_fZ, &b_pim_p_fP_fZ);
fChain->SetBranchAddress("pim.p.fE", pim_p_fE, &b_pim_p_fE);
fChain->SetBranchAddress("pim.x.fUniqueID", pim_x_fUniqueID, &b_pim_x_fUniqueID);
fChain->SetBranchAddress("pim.x.fBits", pim_x_fBits, &b_pim_x_fBits);
fChain->SetBranchAddress("pim.x.fX", pim_x_fX, &b_pim_x_fX);
fChain->SetBranchAddress("pim.x.fY", pim_x_fY, &b_pim_x_fY);
fChain->SetBranchAddress("pim.x.fZ", pim_x_fZ, &b_pim_x_fZ);
fChain->SetBranchAddress("pim.is_fiducial", pim_is_fiducial, &b_pim_is_fiducial);
fChain->SetBranchAddress("pim.chisq", pim_chisq, &b_pim_chisq);
fChain->SetBranchAddress("pim.Ndof", pim_Ndof, &b_pim_Ndof);
fChain->SetBranchAddress("pim.FOM_pid", pim_FOM_pid, &b_pim_FOM_pid);
fChain->SetBranchAddress("Kp", &Kp_, &b_R_Kp_);
fChain->SetBranchAddress("Kp.fUniqueID", &Kp_fUniqueID, &b_Kp_fUniqueID);
fChain->SetBranchAddress("Kp.fBits", &Kp_fBits, &b_Kp_fBits);
fChain->SetBranchAddress("Kp.p.fUniqueID", &Kp_p_fUniqueID, &b_Kp_p_fUniqueID);
fChain->SetBranchAddress("Kp.p.fBits", &Kp_p_fBits, &b_Kp_p_fBits);
fChain->SetBranchAddress("Kp.p.fP.fUniqueID", &Kp_p_fP_fUniqueID, &b_Kp_p_fP_fUniqueID);
fChain->SetBranchAddress("Kp.p.fP.fBits", &Kp_p_fP_fBits, &b_Kp_p_fP_fBits);
fChain->SetBranchAddress("Kp.p.fP.fX", &Kp_p_fP_fX, &b_Kp_p_fP_fX);
fChain->SetBranchAddress("Kp.p.fP.fY", &Kp_p_fP_fY, &b_Kp_p_fP_fY);
fChain->SetBranchAddress("Kp.p.fP.fZ", &Kp_p_fP_fZ, &b_Kp_p_fP_fZ);
fChain->SetBranchAddress("Kp.p.fE", &Kp_p_fE, &b_Kp_p_fE);
fChain->SetBranchAddress("Kp.x.fUniqueID", &Kp_x_fUniqueID, &b_Kp_x_fUniqueID);
fChain->SetBranchAddress("Kp.x.fBits", &Kp_x_fBits, &b_Kp_x_fBits);
fChain->SetBranchAddress("Kp.x.fX", &Kp_x_fX, &b_Kp_x_fX);
fChain->SetBranchAddress("Kp.x.fY", &Kp_x_fY, &b_Kp_x_fY);
fChain->SetBranchAddress("Kp.x.fZ", &Kp_x_fZ, &b_Kp_x_fZ);
fChain->SetBranchAddress("Kp.is_fiducial", &Kp_is_fiducial, &b_Kp_is_fiducial);
fChain->SetBranchAddress("Kp.chisq", &Kp_chisq, &b_Kp_chisq);
fChain->SetBranchAddress("Kp.Ndof", &Kp_Ndof, &b_Kp_Ndof);
fChain->SetBranchAddress("Kp.FOM_pid", &Kp_FOM_pid, &b_Kp_FOM_pid);
fChain->SetBranchAddress("Km", &Km_, &b_R_Km_);
fChain->SetBranchAddress("Km.fUniqueID", &Km_fUniqueID, &b_Km_fUniqueID);
fChain->SetBranchAddress("Km.fBits", &Km_fBits, &b_Km_fBits);
fChain->SetBranchAddress("Km.p.fUniqueID", &Km_p_fUniqueID, &b_Km_p_fUniqueID);
fChain->SetBranchAddress("Km.p.fBits", &Km_p_fBits, &b_Km_p_fBits);
fChain->SetBranchAddress("Km.p.fP.fUniqueID", &Km_p_fP_fUniqueID, &b_Km_p_fP_fUniqueID);
fChain->SetBranchAddress("Km.p.fP.fBits", &Km_p_fP_fBits, &b_Km_p_fP_fBits);
fChain->SetBranchAddress("Km.p.fP.fX", &Km_p_fP_fX, &b_Km_p_fP_fX);
fChain->SetBranchAddress("Km.p.fP.fY", &Km_p_fP_fY, &b_Km_p_fP_fY);
fChain->SetBranchAddress("Km.p.fP.fZ", &Km_p_fP_fZ, &b_Km_p_fP_fZ);
fChain->SetBranchAddress("Km.p.fE", &Km_p_fE, &b_Km_p_fE);
fChain->SetBranchAddress("Km.x.fUniqueID", &Km_x_fUniqueID, &b_Km_x_fUniqueID);
fChain->SetBranchAddress("Km.x.fBits", &Km_x_fBits, &b_Km_x_fBits);
fChain->SetBranchAddress("Km.x.fX", &Km_x_fX, &b_Km_x_fX);
fChain->SetBranchAddress("Km.x.fY", &Km_x_fY, &b_Km_x_fY);
fChain->SetBranchAddress("Km.x.fZ", &Km_x_fZ, &b_Km_x_fZ);
fChain->SetBranchAddress("Km.is_fiducial", &Km_is_fiducial, &b_Km_is_fiducial);
fChain->SetBranchAddress("Km.chisq", &Km_chisq, &b_Km_chisq);
fChain->SetBranchAddress("Km.Ndof", &Km_Ndof, &b_Km_Ndof);
fChain->SetBranchAddress("Km.FOM_pid", &Km_FOM_pid, &b_Km_FOM_pid);
fChain->SetBranchAddress("proton", &proton_, &b_R_proton_);
fChain->SetBranchAddress("proton.fUniqueID", proton_fUniqueID, &b_proton_fUniqueID);
fChain->SetBranchAddress("proton.fBits", proton_fBits, &b_proton_fBits);
fChain->SetBranchAddress("proton.p.fUniqueID", proton_p_fUniqueID, &b_proton_p_fUniqueID);
fChain->SetBranchAddress("proton.p.fBits", proton_p_fBits, &b_proton_p_fBits);
fChain->SetBranchAddress("proton.p.fP.fUniqueID", proton_p_fP_fUniqueID, &b_proton_p_fP_fUniqueID);
fChain->SetBranchAddress("proton.p.fP.fBits", proton_p_fP_fBits, &b_proton_p_fP_fBits);
fChain->SetBranchAddress("proton.p.fP.fX", proton_p_fP_fX, &b_proton_p_fP_fX);
fChain->SetBranchAddress("proton.p.fP.fY", proton_p_fP_fY, &b_proton_p_fP_fY);
fChain->SetBranchAddress("proton.p.fP.fZ", proton_p_fP_fZ, &b_proton_p_fP_fZ);
fChain->SetBranchAddress("proton.p.fE", proton_p_fE, &b_proton_p_fE);
fChain->SetBranchAddress("proton.x.fUniqueID", proton_x_fUniqueID, &b_proton_x_fUniqueID);
fChain->SetBranchAddress("proton.x.fBits", proton_x_fBits, &b_proton_x_fBits);
fChain->SetBranchAddress("proton.x.fX", proton_x_fX, &b_proton_x_fX);
fChain->SetBranchAddress("proton.x.fY", proton_x_fY, &b_proton_x_fY);
fChain->SetBranchAddress("proton.x.fZ", proton_x_fZ, &b_proton_x_fZ);
fChain->SetBranchAddress("proton.is_fiducial", proton_is_fiducial, &b_proton_is_fiducial);
fChain->SetBranchAddress("proton.chisq", proton_chisq, &b_proton_chisq);
fChain->SetBranchAddress("proton.Ndof", proton_Ndof, &b_proton_Ndof);
fChain->SetBranchAddress("proton.FOM_pid", proton_FOM_pid, &b_proton_FOM_pid);
fChain->SetBranchAddress("photon", &photon_, &b_R_photon_);
fChain->SetBranchAddress("photon.fUniqueID", photon_fUniqueID, &b_photon_fUniqueID);
fChain->SetBranchAddress("photon.fBits", photon_fBits, &b_photon_fBits);
fChain->SetBranchAddress("photon.p.fUniqueID", photon_p_fUniqueID, &b_photon_p_fUniqueID);
fChain->SetBranchAddress("photon.p.fBits", photon_p_fBits, &b_photon_p_fBits);
fChain->SetBranchAddress("photon.p.fP.fUniqueID", photon_p_fP_fUniqueID, &b_photon_p_fP_fUniqueID);
fChain->SetBranchAddress("photon.p.fP.fBits", photon_p_fP_fBits, &b_photon_p_fP_fBits);
fChain->SetBranchAddress("photon.p.fP.fX", photon_p_fP_fX, &b_photon_p_fP_fX);
fChain->SetBranchAddress("photon.p.fP.fY", photon_p_fP_fY, &b_photon_p_fP_fY);
fChain->SetBranchAddress("photon.p.fP.fZ", photon_p_fP_fZ, &b_photon_p_fP_fZ);
fChain->SetBranchAddress("photon.p.fE", photon_p_fE, &b_photon_p_fE);
fChain->SetBranchAddress("photon.x.fUniqueID", photon_x_fUniqueID, &b_photon_x_fUniqueID);
fChain->SetBranchAddress("photon.x.fBits", photon_x_fBits, &b_photon_x_fBits);
fChain->SetBranchAddress("photon.x.fX", photon_x_fX, &b_photon_x_fX);
fChain->SetBranchAddress("photon.x.fY", photon_x_fY, &b_photon_x_fY);
fChain->SetBranchAddress("photon.x.fZ", photon_x_fZ, &b_photon_x_fZ);
fChain->SetBranchAddress("photon.is_fiducial", photon_is_fiducial, &b_photon_is_fiducial);
fChain->SetBranchAddress("photon.chisq", photon_chisq, &b_photon_chisq);
fChain->SetBranchAddress("photon.Ndof", photon_Ndof, &b_photon_Ndof);
fChain->SetBranchAddress("photon.FOM_pid", photon_FOM_pid, &b_photon_FOM_pid);
fChain->SetBranchAddress("pip_match", &pip_match_, &b_R_pip_match_);
fChain->SetBranchAddress("pip_match.fUniqueID", pip_match_fUniqueID, &b_pip_match_fUniqueID);
fChain->SetBranchAddress("pip_match.fBits", pip_match_fBits, &b_pip_match_fBits);
fChain->SetBranchAddress("pip_match.p.fUniqueID", pip_match_p_fUniqueID, &b_pip_match_p_fUniqueID);
fChain->SetBranchAddress("pip_match.p.fBits", pip_match_p_fBits, &b_pip_match_p_fBits);
fChain->SetBranchAddress("pip_match.p.fP.fUniqueID", pip_match_p_fP_fUniqueID, &b_pip_match_p_fP_fUniqueID);
fChain->SetBranchAddress("pip_match.p.fP.fBits", pip_match_p_fP_fBits, &b_pip_match_p_fP_fBits);
fChain->SetBranchAddress("pip_match.p.fP.fX", pip_match_p_fP_fX, &b_pip_match_p_fP_fX);
fChain->SetBranchAddress("pip_match.p.fP.fY", pip_match_p_fP_fY, &b_pip_match_p_fP_fY);
fChain->SetBranchAddress("pip_match.p.fP.fZ", pip_match_p_fP_fZ, &b_pip_match_p_fP_fZ);
fChain->SetBranchAddress("pip_match.p.fE", pip_match_p_fE, &b_pip_match_p_fE);
fChain->SetBranchAddress("pip_match.x.fUniqueID", pip_match_x_fUniqueID, &b_pip_match_x_fUniqueID);
fChain->SetBranchAddress("pip_match.x.fBits", pip_match_x_fBits, &b_pip_match_x_fBits);
fChain->SetBranchAddress("pip_match.x.fX", pip_match_x_fX, &b_pip_match_x_fX);
fChain->SetBranchAddress("pip_match.x.fY", pip_match_x_fY, &b_pip_match_x_fY);
fChain->SetBranchAddress("pip_match.x.fZ", pip_match_x_fZ, &b_pip_match_x_fZ);
fChain->SetBranchAddress("pip_match.is_fiducial", pip_match_is_fiducial, &b_pip_match_is_fiducial);
fChain->SetBranchAddress("pip_match.chisq", pip_match_chisq, &b_pip_match_chisq);
fChain->SetBranchAddress("pip_match.Ndof", pip_match_Ndof, &b_pip_match_Ndof);
fChain->SetBranchAddress("pip_match.FOM_pid", pip_match_FOM_pid, &b_pip_match_FOM_pid);
fChain->SetBranchAddress("pim_match", &pim_match_, &b_R_pim_match_);
fChain->SetBranchAddress("pim_match.fUniqueID", pim_match_fUniqueID, &b_pim_match_fUniqueID);
fChain->SetBranchAddress("pim_match.fBits", pim_match_fBits, &b_pim_match_fBits);
fChain->SetBranchAddress("pim_match.p.fUniqueID", pim_match_p_fUniqueID, &b_pim_match_p_fUniqueID);
fChain->SetBranchAddress("pim_match.p.fBits", pim_match_p_fBits, &b_pim_match_p_fBits);
fChain->SetBranchAddress("pim_match.p.fP.fUniqueID", pim_match_p_fP_fUniqueID, &b_pim_match_p_fP_fUniqueID);
fChain->SetBranchAddress("pim_match.p.fP.fBits", pim_match_p_fP_fBits, &b_pim_match_p_fP_fBits);
fChain->SetBranchAddress("pim_match.p.fP.fX", pim_match_p_fP_fX, &b_pim_match_p_fP_fX);
fChain->SetBranchAddress("pim_match.p.fP.fY", pim_match_p_fP_fY, &b_pim_match_p_fP_fY);
fChain->SetBranchAddress("pim_match.p.fP.fZ", pim_match_p_fP_fZ, &b_pim_match_p_fP_fZ);
fChain->SetBranchAddress("pim_match.p.fE", pim_match_p_fE, &b_pim_match_p_fE);
fChain->SetBranchAddress("pim_match.x.fUniqueID", pim_match_x_fUniqueID, &b_pim_match_x_fUniqueID);
fChain->SetBranchAddress("pim_match.x.fBits", pim_match_x_fBits, &b_pim_match_x_fBits);
fChain->SetBranchAddress("pim_match.x.fX", pim_match_x_fX, &b_pim_match_x_fX);
fChain->SetBranchAddress("pim_match.x.fY", pim_match_x_fY, &b_pim_match_x_fY);
fChain->SetBranchAddress("pim_match.x.fZ", pim_match_x_fZ, &b_pim_match_x_fZ);
fChain->SetBranchAddress("pim_match.is_fiducial", pim_match_is_fiducial, &b_pim_match_is_fiducial);
fChain->SetBranchAddress("pim_match.chisq", pim_match_chisq, &b_pim_match_chisq);
fChain->SetBranchAddress("pim_match.Ndof", pim_match_Ndof, &b_pim_match_Ndof);
fChain->SetBranchAddress("pim_match.FOM_pid", pim_match_FOM_pid, &b_pim_match_FOM_pid);
fChain->SetBranchAddress("Kp_match", &Kp_match_, &b_R_Kp_match_);
fChain->SetBranchAddress("Kp_match.fUniqueID", &Kp_match_fUniqueID, &b_Kp_match_fUniqueID);
fChain->SetBranchAddress("Kp_match.fBits", &Kp_match_fBits, &b_Kp_match_fBits);
fChain->SetBranchAddress("Kp_match.p.fUniqueID", &Kp_match_p_fUniqueID, &b_Kp_match_p_fUniqueID);
fChain->SetBranchAddress("Kp_match.p.fBits", &Kp_match_p_fBits, &b_Kp_match_p_fBits);
fChain->SetBranchAddress("Kp_match.p.fP.fUniqueID", &Kp_match_p_fP_fUniqueID, &b_Kp_match_p_fP_fUniqueID);
fChain->SetBranchAddress("Kp_match.p.fP.fBits", &Kp_match_p_fP_fBits, &b_Kp_match_p_fP_fBits);
fChain->SetBranchAddress("Kp_match.p.fP.fX", &Kp_match_p_fP_fX, &b_Kp_match_p_fP_fX);
fChain->SetBranchAddress("Kp_match.p.fP.fY", &Kp_match_p_fP_fY, &b_Kp_match_p_fP_fY);
fChain->SetBranchAddress("Kp_match.p.fP.fZ", &Kp_match_p_fP_fZ, &b_Kp_match_p_fP_fZ);
fChain->SetBranchAddress("Kp_match.p.fE", &Kp_match_p_fE, &b_Kp_match_p_fE);
fChain->SetBranchAddress("Kp_match.x.fUniqueID", &Kp_match_x_fUniqueID, &b_Kp_match_x_fUniqueID);
fChain->SetBranchAddress("Kp_match.x.fBits", &Kp_match_x_fBits, &b_Kp_match_x_fBits);
fChain->SetBranchAddress("Kp_match.x.fX", &Kp_match_x_fX, &b_Kp_match_x_fX);
fChain->SetBranchAddress("Kp_match.x.fY", &Kp_match_x_fY, &b_Kp_match_x_fY);
fChain->SetBranchAddress("Kp_match.x.fZ", &Kp_match_x_fZ, &b_Kp_match_x_fZ);
fChain->SetBranchAddress("Kp_match.is_fiducial", &Kp_match_is_fiducial, &b_Kp_match_is_fiducial);
fChain->SetBranchAddress("Kp_match.chisq", &Kp_match_chisq, &b_Kp_match_chisq);
fChain->SetBranchAddress("Kp_match.Ndof", &Kp_match_Ndof, &b_Kp_match_Ndof);
fChain->SetBranchAddress("Kp_match.FOM_pid", &Kp_match_FOM_pid, &b_Kp_match_FOM_pid);
fChain->SetBranchAddress("Km_match", &Km_match_, &b_R_Km_match_);
fChain->SetBranchAddress("Km_match.fUniqueID", &Km_match_fUniqueID, &b_Km_match_fUniqueID);
fChain->SetBranchAddress("Km_match.fBits", &Km_match_fBits, &b_Km_match_fBits);
fChain->SetBranchAddress("Km_match.p.fUniqueID", &Km_match_p_fUniqueID, &b_Km_match_p_fUniqueID);
fChain->SetBranchAddress("Km_match.p.fBits", &Km_match_p_fBits, &b_Km_match_p_fBits);
fChain->SetBranchAddress("Km_match.p.fP.fUniqueID", &Km_match_p_fP_fUniqueID, &b_Km_match_p_fP_fUniqueID);
fChain->SetBranchAddress("Km_match.p.fP.fBits", &Km_match_p_fP_fBits, &b_Km_match_p_fP_fBits);
fChain->SetBranchAddress("Km_match.p.fP.fX", &Km_match_p_fP_fX, &b_Km_match_p_fP_fX);
fChain->SetBranchAddress("Km_match.p.fP.fY", &Km_match_p_fP_fY, &b_Km_match_p_fP_fY);
fChain->SetBranchAddress("Km_match.p.fP.fZ", &Km_match_p_fP_fZ, &b_Km_match_p_fP_fZ);
fChain->SetBranchAddress("Km_match.p.fE", &Km_match_p_fE, &b_Km_match_p_fE);
fChain->SetBranchAddress("Km_match.x.fUniqueID", &Km_match_x_fUniqueID, &b_Km_match_x_fUniqueID);
fChain->SetBranchAddress("Km_match.x.fBits", &Km_match_x_fBits, &b_Km_match_x_fBits);
fChain->SetBranchAddress("Km_match.x.fX", &Km_match_x_fX, &b_Km_match_x_fX);
fChain->SetBranchAddress("Km_match.x.fY", &Km_match_x_fY, &b_Km_match_x_fY);
fChain->SetBranchAddress("Km_match.x.fZ", &Km_match_x_fZ, &b_Km_match_x_fZ);
fChain->SetBranchAddress("Km_match.is_fiducial", &Km_match_is_fiducial, &b_Km_match_is_fiducial);
fChain->SetBranchAddress("Km_match.chisq", &Km_match_chisq, &b_Km_match_chisq);
fChain->SetBranchAddress("Km_match.Ndof", &Km_match_Ndof, &b_Km_match_Ndof);
fChain->SetBranchAddress("Km_match.FOM_pid", &Km_match_FOM_pid, &b_Km_match_FOM_pid);
fChain->SetBranchAddress("proton_match", &proton_match_, &b_R_proton_match_);
fChain->SetBranchAddress("proton_match.fUniqueID", proton_match_fUniqueID, &b_proton_match_fUniqueID);
fChain->SetBranchAddress("proton_match.fBits", proton_match_fBits, &b_proton_match_fBits);
fChain->SetBranchAddress("proton_match.p.fUniqueID", proton_match_p_fUniqueID, &b_proton_match_p_fUniqueID);
fChain->SetBranchAddress("proton_match.p.fBits", proton_match_p_fBits, &b_proton_match_p_fBits);
fChain->SetBranchAddress("proton_match.p.fP.fUniqueID", proton_match_p_fP_fUniqueID, &b_proton_match_p_fP_fUniqueID);
fChain->SetBranchAddress("proton_match.p.fP.fBits", proton_match_p_fP_fBits, &b_proton_match_p_fP_fBits);
fChain->SetBranchAddress("proton_match.p.fP.fX", proton_match_p_fP_fX, &b_proton_match_p_fP_fX);
fChain->SetBranchAddress("proton_match.p.fP.fY", proton_match_p_fP_fY, &b_proton_match_p_fP_fY);
fChain->SetBranchAddress("proton_match.p.fP.fZ", proton_match_p_fP_fZ, &b_proton_match_p_fP_fZ);
fChain->SetBranchAddress("proton_match.p.fE", proton_match_p_fE, &b_proton_match_p_fE);
fChain->SetBranchAddress("proton_match.x.fUniqueID", proton_match_x_fUniqueID, &b_proton_match_x_fUniqueID);
fChain->SetBranchAddress("proton_match.x.fBits", proton_match_x_fBits, &b_proton_match_x_fBits);
fChain->SetBranchAddress("proton_match.x.fX", proton_match_x_fX, &b_proton_match_x_fX);
fChain->SetBranchAddress("proton_match.x.fY", proton_match_x_fY, &b_proton_match_x_fY);
fChain->SetBranchAddress("proton_match.x.fZ", proton_match_x_fZ, &b_proton_match_x_fZ);
fChain->SetBranchAddress("proton_match.is_fiducial", proton_match_is_fiducial, &b_proton_match_is_fiducial);
fChain->SetBranchAddress("proton_match.chisq", proton_match_chisq, &b_proton_match_chisq);
fChain->SetBranchAddress("proton_match.Ndof", proton_match_Ndof, &b_proton_match_Ndof);
fChain->SetBranchAddress("proton_match.FOM_pid", proton_match_FOM_pid, &b_proton_match_FOM_pid);
fChain->SetBranchAddress("photon_match", &photon_match_, &b_R_photon_match_);
fChain->SetBranchAddress("photon_match.fUniqueID", photon_match_fUniqueID, &b_photon_match_fUniqueID);
fChain->SetBranchAddress("photon_match.fBits", photon_match_fBits, &b_photon_match_fBits);
fChain->SetBranchAddress("photon_match.p.fUniqueID", photon_match_p_fUniqueID, &b_photon_match_p_fUniqueID);
fChain->SetBranchAddress("photon_match.p.fBits", photon_match_p_fBits, &b_photon_match_p_fBits);
fChain->SetBranchAddress("photon_match.p.fP.fUniqueID", photon_match_p_fP_fUniqueID, &b_photon_match_p_fP_fUniqueID);
fChain->SetBranchAddress("photon_match.p.fP.fBits", photon_match_p_fP_fBits, &b_photon_match_p_fP_fBits);
fChain->SetBranchAddress("photon_match.p.fP.fX", photon_match_p_fP_fX, &b_photon_match_p_fP_fX);
fChain->SetBranchAddress("photon_match.p.fP.fY", photon_match_p_fP_fY, &b_photon_match_p_fP_fY);
fChain->SetBranchAddress("photon_match.p.fP.fZ", photon_match_p_fP_fZ, &b_photon_match_p_fP_fZ);
fChain->SetBranchAddress("photon_match.p.fE", photon_match_p_fE, &b_photon_match_p_fE);
fChain->SetBranchAddress("photon_match.x.fUniqueID", photon_match_x_fUniqueID, &b_photon_match_x_fUniqueID);
fChain->SetBranchAddress("photon_match.x.fBits", photon_match_x_fBits, &b_photon_match_x_fBits);
fChain->SetBranchAddress("photon_match.x.fX", photon_match_x_fX, &b_photon_match_x_fX);
fChain->SetBranchAddress("photon_match.x.fY", photon_match_x_fY, &b_photon_match_x_fY);
fChain->SetBranchAddress("photon_match.x.fZ", photon_match_x_fZ, &b_photon_match_x_fZ);
fChain->SetBranchAddress("photon_match.is_fiducial", photon_match_is_fiducial, &b_photon_match_is_fiducial);
fChain->SetBranchAddress("photon_match.chisq", photon_match_chisq, &b_photon_match_chisq);
fChain->SetBranchAddress("photon_match.Ndof", photon_match_Ndof, &b_photon_match_Ndof);
fChain->SetBranchAddress("photon_match.FOM_pid", photon_match_FOM_pid, &b_photon_match_FOM_pid);
fChain->SetBranchAddress("target.fUniqueID", &target_fUniqueID, &b_R_target_fUniqueID);
fChain->SetBranchAddress("target.fBits", &target_fBits, &b_R_target_fBits);
fChain->SetBranchAddress("target.fP.fUniqueID", &target_fP_fUniqueID, &b_R_target_fP_fUniqueID);
fChain->SetBranchAddress("target.fP.fBits", &target_fP_fBits, &b_R_target_fP_fBits);
fChain->SetBranchAddress("target.fP.fX", &target_fP_fX, &b_R_target_fP_fX);
fChain->SetBranchAddress("target.fP.fY", &target_fP_fY, &b_R_target_fP_fY);
fChain->SetBranchAddress("target.fP.fZ", &target_fP_fZ, &b_R_target_fP_fZ);
fChain->SetBranchAddress("target.fE", &target_fE, &b_R_target_fE);
fChain->SetBranchAddress("beam.fUniqueID", &beam_fUniqueID, &b_R_beam_fUniqueID);
fChain->SetBranchAddress("beam.fBits", &beam_fBits, &b_R_beam_fBits);
fChain->SetBranchAddress("beam.fP.fUniqueID", &beam_fP_fUniqueID, &b_R_beam_fP_fUniqueID);
fChain->SetBranchAddress("beam.fP.fBits", &beam_fP_fBits, &b_R_beam_fP_fBits);
fChain->SetBranchAddress("beam.fP.fX", &beam_fP_fX, &b_R_beam_fP_fX);
fChain->SetBranchAddress("beam.fP.fY", &beam_fP_fY, &b_R_beam_fP_fY);
fChain->SetBranchAddress("beam.fP.fZ", &beam_fP_fZ, &b_R_beam_fP_fZ);
fChain->SetBranchAddress("beam.fE", &beam_fE, &b_R_beam_fE);
fChain->SetBranchAddress("vertex.fUniqueID", &vertex_fUniqueID, &b_R_vertex_fUniqueID);
fChain->SetBranchAddress("vertex.fBits", &vertex_fBits, &b_R_vertex_fBits);
fChain->SetBranchAddress("vertex.fX", &vertex_fX, &b_R_vertex_fX);
fChain->SetBranchAddress("vertex.fY", &vertex_fY, &b_R_vertex_fY);
fChain->SetBranchAddress("vertex.fZ", &vertex_fZ, &b_R_vertex_fZ);
fChain->SetBranchAddress("W.fUniqueID", &W_fUniqueID, &b_R_W_fUniqueID);
fChain->SetBranchAddress("W.fBits", &W_fBits, &b_R_W_fBits);
fChain->SetBranchAddress("W.fP.fUniqueID", &W_fP_fUniqueID, &b_R_W_fP_fUniqueID);
fChain->SetBranchAddress("W.fP.fBits", &W_fP_fBits, &b_R_W_fP_fBits);
fChain->SetBranchAddress("W.fP.fX", &W_fP_fX, &b_R_W_fP_fX);
fChain->SetBranchAddress("W.fP.fY", &W_fP_fY, &b_R_W_fP_fY);
fChain->SetBranchAddress("W.fP.fZ", &W_fP_fZ, &b_R_W_fP_fZ);
fChain->SetBranchAddress("W.fE", &W_fE, &b_R_W_fE);
fChain->SetBranchAddress("all_fiducial", &all_fiducial, &b_R_all_fiducial);
fChain->SetBranchAddress("all_mesons_fiducial", &all_mesons_fiducial, &b_R_all_mesons_fiducial);
fChain->SetBranchAddress("all_photons_fiducial", &all_photons_fiducial, &b_R_all_photons_fiducial);
fChain->SetBranchAddress("all_protons_fiducial", &all_protons_fiducial, &b_R_all_protons_fiducial);
}
Bool_t b1_pi::Notify()
{
// The Notify() function is called when a new file is opened. This
// can be either for a new TTree in a TChain or when when a new TTree
// is started when using PROOF. It is normally not necessary to make changes
// to the generated code, but the routine can be extended by the
// user if needed. The return value is currently not used.
return kTRUE;
}
#endif // #ifdef b1_pi_cxx
b1_pi.C
#define b1_pi_cxx
// The class definition in b1_pi.h has been generated automatically
// by the ROOT utility TTree::MakeSelector(). This class is derived
// from the ROOT class TSelector. For more information on the TSelector
// framework see $ROOTSYS/README/README.SELECTOR or the ROOT User Manual.
// The following methods are defined in this file:
// Begin(): called every time a loop on the tree starts,
// a convenient place to create your histograms.
// SlaveBegin(): called after Begin(), when on PROOF called only on the
// slave servers.
// Process(): called for each event, in this function you decide what
// to read and fill your histograms.
// SlaveTerminate: called at the end of the loop on the tree, when on PROOF
// called only on the slave servers.
// Terminate(): called at the end of the loop on the tree,
// a convenient place to draw/fit your histograms.
//
// To use this file, try the following session on your Tree T:
//
// Root > T->Process("b1_pi.C")
// Root > T->Process("b1_pi.C","some options")
// Root > T->Process("b1_pi.C+")
//
#include "b1_pi.h"
#include <TH2.h>
#include <TStyle.h>
void b1_pi::Begin(TTree * /*tree*/)
{
// The Begin() function is called at the start of the query.
// When running with PROOF Begin() is only called on the client.
// The tree argument is deprecated (on PROOF 0 is passed).
TString option = GetOption();
//-------------------------------------
// Instantiate my histograms
mass_X = new TH1D("mass_X", "Mass of Exotic Hybrid", 500, 0.0, 4.0);
mass_b1 = new TH1D("mass_b1", "Mass of b1 meson", 200, 0.0, 2.0);
mass_omega = new TH1D("mass_omega", "Mass of #omega meson", 500, 0.0, 2.0);
mass_rho = new TH1D("mass_rho", "Mass of #rho meson", 200, 0.0, 2.0);
mass_pi0 = new TH1D("mass_pi0", "Mass of #pi^{o} meson", 500, 0.0, 0.5);
//-------------------------------------
}
void b1_pi::SlaveBegin(TTree * /*tree*/)
{
// The SlaveBegin() function is called after the Begin() function.
// When running with PROOF SlaveBegin() is called on each slave server.
// The tree argument is deprecated (on PROOF 0 is passed).
TString option = GetOption();
}
Bool_t b1_pi::Process(Long64_t entry)
{
// The Process() function is called for each entry in the tree (or possibly
// keyed object in the case of PROOF) to be processed. The entry argument
// specifies which entry in the currently loaded tree is to be processed.
// It can be passed to either b1_pi::GetEntry() or TBranch::GetEntry()
// to read either all or the required parts of the data. When processing
// keyed objects with PROOF, the object is already loaded and is available
// via the fObject pointer.
//
// This function should contain the "body" of the analysis. It can contain
// simple or elaborate selection criteria, run algorithms on the data
// of the event and typically fill histograms.
//
// The processing can be stopped by calling Abort().
//
// Use fStatus to set the return value of TTree::Process().
//
// The return value is currently not used.
//-------------------------------------
// Read the event into memory
GetEntry(entry);
//-------------------------------------
// For some reason, ROOT no longer makes the objects available as objects (???)
// so we have to explicitly re-create out TLorentzVector objects
// (Not sure why I'm forced to do this!!)
if(Nphoton>kMaxphoton)Nphoton = kMaxphoton;
if(Nproton>kMaxproton)Nproton = kMaxproton;
if(Npip>kMaxpip)Npip = kMaxpip;
if(Npim>kMaxpim)Npim = kMaxpim;
TLorentzVector *photon = new TLorentzVector[Nphoton];
TLorentzVector *pip = new TLorentzVector[Npip];
TLorentzVector *pim = new TLorentzVector[Npim];
TLorentzVector *proton = new TLorentzVector[Nproton];
for(Int_t i=0; i<Nphoton; i++)photon[i].SetPxPyPzE(photon_p_fP_fX[i], photon_p_fP_fY[i], photon_p_fP_fZ[i], photon_p_fE[i]);
for(Int_t i=0; i<Npip; i++)pip[i].SetPxPyPzE(pip_p_fP_fX[i], pip_p_fP_fY[i], pip_p_fP_fZ[i], pip_p_fE[i]);
for(Int_t i=0; i<Npim; i++)pim[i].SetPxPyPzE(pim_p_fP_fX[i], pim_p_fP_fY[i], pim_p_fP_fZ[i], pim_p_fE[i]);
for(Int_t i=0; i<Nproton; i++)proton[i].SetPxPyPzE(proton_p_fP_fX[i], proton_p_fP_fY[i], proton_p_fP_fZ[i], proton_p_fE[i]);
TLorentzVector target(target_fP_fX, target_fP_fY, target_fP_fZ, target_fE);
TLorentzVector beam(beam_fP_fX, beam_fP_fY, beam_fP_fZ, beam_fE);
TVector3 vertex(vertex_fX, vertex_fY, vertex_fZ);
TLorentzVector W(W_fP_fX, W_fP_fY, W_fP_fZ, W_fE);
//-------------------------------------
// Fill histograms
// Note: Below we loop over all combinations but keep only the
// one closest the relevant meson mass.
// Find the 2 gammas that best reconstruct a pi0
TLorentzVector pi0(0,0,0,0);
if(Nphoton>0){
for(Int_t i=0; i<Nphoton-1; i++){
for(Int_t j=i+1; j<Nphoton; j++){
TLorentzVector my_pi0 = photon[i]+photon[j];
if(fabs(my_pi0.M()-0.135) < fabs(pi0.M()-0.135))pi0 = my_pi0;
}
}
}
// Uncomment the following to force the pi0 to correct mass
//pi0.SetE(sqrt(pow(pi0.P(),2) + pow(0.135,2)));
// Find pi+pi+pi-pi-pi0 particles that give us the best X(2000) mass
TLorentzVector X(0,0,0,0);
Int_t ipip1=-1, ipip2=-1;
Int_t ipim1=-1, ipim2=-1;
if(pi0.P()!=0.0 && Npip>0 && Npim>0){
for(Int_t k=0; k<Npip-1; k++){
for(Int_t m=k+1; m<Npip; m++){
for(Int_t n=0; n<Npim-1; n++){
for(Int_t p=n+1; p<Npim; p++){
TLorentzVector my_X = pi0+pip[k]+pip[m]+pim[n]+pim[p];
if(fabs(my_X.M()-2.000) < fabs(X.M()-2.000)){
X = my_X;
ipip1 = k;
ipip2 = m;
ipim1 = n;
ipim2 = p;
}
}
}
}
}
}
// Fill all histograms if an "X" was found
if(X.P()!=0.0){
mass_pi0->Fill(pi0.M());
mass_X->Fill(X.M());
// b1
// (find the pi- that did NOT come from the X decay directly)
TLorentzVector b1;
TLorentzVector b1_1 = X - pim[ipim1];
TLorentzVector b1_2 = X - pim[ipim2];
Int_t index_rho_pim; // which ever is not from the X must be from the rho
if(fabs(b1_1.M()-1.235) < fabs(b1_2.M()-1.235)){
b1 = b1_1;
index_rho_pim = ipim1;
}else{
b1 = b1_2;
index_rho_pim = ipim2;
}
mass_b1->Fill(b1.M());
// omega
// (find the pi+ that did NOT come from the b1 decay directly)
TLorentzVector omega;
TLorentzVector omega_1 = b1 - pip[ipip1];
TLorentzVector omega_2 = b1 - pip[ipip2];
Int_t index_rho_pip; // which ever is not from the b1 must be from the rho
if(fabs(omega_1.M()-0.783) < fabs(omega_2.M()-0.783)){
omega = omega_1;
index_rho_pip = ipip1;
}else{
omega = omega_2;
index_rho_pip = ipip2;
}
mass_omega->Fill(omega.M());
// rho
//TLorentzVector rho = pip[index_rho_pip] + pim[index_rho_pim];
TLorentzVector rho = omega - pi0;
mass_rho->Fill(rho.M());
} // if(X.P()!=0.0)
//-------------------------------------
// Delete our TLorentzVector objects
if(Nphoton>0)delete[] photon;
if(Npip>0)delete[] pip;
if(Npim>0)delete[] pim;
if(Nproton>0)delete[] proton;
return kTRUE;
}
void b1_pi::SlaveTerminate()
{
// The SlaveTerminate() function is called after all entries or objects
// have been processed. When running with PROOF SlaveTerminate() is called
// on each slave server.
}
void b1_pi::Terminate()
{
// The Terminate() function is the last function to be called during
// a query. It always runs on the client, it can be used to present
// the results graphically or save the results to file.
}
