#include <MillRecoTree.h>

Inheritance diagram for garana::MillRecoTree:

Public Member Functions
	MillRecoTree ()

	MillRecoTree (TTree treeIn, TTree treeOut)

void	MillTrees () override
	main method More...

Public Member Functions inherited from garana::Mill
virtual	~Mill ()

bool	IsVerified () const

std::string	CharStarToString (const char *cstr)

Public Member Functions inherited from garana::FlatRecoTree
	FlatRecoTree ()

	FlatRecoTree (TTree *tree)

	FlatRecoTree (TTree *tree, char opt)

const size_t	NTrack () const override
	number of tracks in this event More...

const size_t	NVertex () const override
	number of vertices in this event More...

const size_t	NVee () const override
	number of vees in this event More...

const size_t	NCalCluster () const override
	number of ECal clusters in this event More...

const TLorentzVector *	TrackVertex (const size_t &itrack) const override
	4-position of track's assumed start point More...

const TLorentzVector *	TrackEnd (const size_t &itrack) const override
	4-position of track's assumed end point More...

const size_t	NTrackHit (const size_t &itrack) const override
	number of reconstructed hits used in track fit More...

const TVector3 *	TrackMomBeg (const size_t &itrack) const override
	momentum 3-vector as track's assumed start point More...

const TVector3 *	TrackMomEnd (const size_t &itrack) const override
	momentum 3-vector as track's assumed end point More...

const float	TrackVtxDirectionX (const size_t &itrack) const override
	x-direction cosine at track vertex More...

const float	TrackVtxDirectionY (const size_t &itrack) const override
	y-direction cosine at track vertex More...

const float	TrackVtxDirectionZ (const size_t &itrack) const override
	z-direction cosine at track vertex More...

const float	TrackEndDirectionX (const size_t &itrack) const override
	x-direction cosine at track end More...

const float	TrackEndDirectionY (const size_t &itrack) const override
	y-direction cosine at track end More...

const float	TrackEndDirectionZ (const size_t &itrack) const override
	z-direction cosine at track end More...

const float	TrackLenFwd (const size_t &itrack) const override
	track length from forward fit More...

const float	TrackLenBkd (const size_t &itrack) const override
	track length from backward fit More...

const float	TrackIonizFwd (const size_t &itrack) const override
	track average ionization rate from forward fit More...

const float	TrackIonizBkd (const size_t &itrack) const override
	track average ionization rate from backward fit More...

const int	TrackChiSqrFwd (const size_t &itrack) const override
	chi-squared of track fit in forward direction More...

const int	TrackChiSqrBkd (const size_t &itrack) const override
	chi-squared of track fit in backward direction More...

const int	TrackChgFwd (const size_t &itrack) const override
	charge sign of track if current hypothesis correct More...

const int	TrackChgBkd (const size_t &itrack) const override
	charge sign of track if current hypothesis wrong More...

void	TrackParBeg (const size_t &itrack, float pars[5]) const override
	track fit parameters at the track's assumed start More...

void	TrackParEnd (const size_t &itrack, float pars[5]) const override
	track fit parameters at the track's assumed end More...

void	TrackCovarBeg (const size_t &itrack, float pars[15]) const override
	track fit covariance matrix (assume symmetry) at track's assumed start More...

void	TrackCovarEnd (const size_t &itrack, float pars[15]) const override
	track fit covariance matrix (assume symmetry) at track's assumed end More...

const TLorentzVector *	TrackTruePosBeg (const size_t &itrack) const override
	true 4-position at track vertex [cm,ns] More...

const TLorentzVector *	TrackTruePosEnd (const size_t &itrack) const override
	true 4-position at track end [cm,ns] More...

const TLorentzVector *	TrackTrueMomBeg (const size_t &itrack) const override
	true 4-momentum at track vertex [GeV/c,GeV] More...

const TLorentzVector *	TrackTrueMomEnd (const size_t &itrack) const override
	true 4-momentum at track end [GeV/c,GeV] More...

const float	TrackTrueEnergy (const size_t &icluster) const override
	total associated true energy deposited with ith cluster More...

const size_t	TrackNTrueTrack (const size_t &icluster) const override
	number of MCParticles associated with ith cluster More...

const int	TrackTrkIdMaxDeposit (const size_t &icluster) const override
	trackID of the MCParticle depositing the most energy in ith cluster More...

const float	TrackMaxDeposit (const size_t &icluster) const override
	maximum true deposited energy from a single MCParticle More...

const pair< int, float > *	TrackTrueDeposit (const size_t &itrack, size_t &itrue) const override

const vector< pair< int, float > > *	TrackTrueDeposits (const size_t &itrack) const override

const TLorentzVector *	GetVertex (const size_t &ivertex) const override
	================ vertex ========================== More...

void	VertexCovariance (const size_t &ivertex, float covar[][3]) const override
	given a vertex index, fill given position covariance matrix More...

const TLorentzVector *	VeeVertex (const size_t &ivee) const override
	=============== Vee ======================= More...

void	VeeCovariance (const size_t &ivee, float covar[][3]) const override
	given a vee index, fill given position covariance matrix More...

const vector< TLorentzVector > *	VeeMomentumPerHypothesis (const size_t &ivee) const override
	4-momentum for vee if it matches {Kshort, Lambda1, Lambda2} More...

const float	VeeChiSquared (const size_t &ivee) const override

const TLorentzVector *	CalClustPosition (const size_t &icluster) const override
	================ ECal cluster ====================== More...

const float	CalClustEnergy (const size_t &icluster) const override
	reconstructed ECal cluster energy More...

const float	CalClustEnergyError (const size_t &icluster) const override

const float	CalClustTrueEnergy (const size_t &icluster) const override
	total associated true energy deposited with ith cluster More...

const size_t	CalClustNTrueTrack (const size_t &icluster) const override
	number of MCParticles associated with ith cluster More...

const int	CalClustTrkIdMaxDeposit (const size_t &icluster) const override
	trackID of the MCParticle depositing the most energy in ith cluster More...

const float	CalClustMaxDeposit (const size_t &icluster) const override
	maximum true deposited energy from a single MCParticle More...

const std::pair< int, float > *	CalClustTrueDeposit (const size_t &icluster, const size_t &itrack) const override

const float	CalClustTimeDifference (const size_t &icluster) const override

const float *	CalClustShape (const size_t &icluster) const override

const float	CalClustTheta (const size_t &icluster) const override

const float	CalClustPhi (const size_t &icluster) const override

const vector< TVector3 > *	CalClustEigenVecs (const size_t &icluster) const override

Public Member Functions inherited from garana::RecoTree
virtual	~RecoTree ()

const float	TrackMaxDepositFrac (const size_t &itrack) const
	largest fraction of total energy contributed by single MCParticle More...

const float	CalClustMaxDepositFrac (const size_t &icluster) const
	largest fraction of total energy contributed by single MCParticle More...

void	GetTrackG4PIndices (const size_t &itrk, vector< UInt_t > &ig4ps) const
	given a track index, fill a given vector with matched G4 particle indices More...

void	GetVertexTrackIndices (const size_t &ivtx, vector< UInt_t > &itracks) const
	given a vertex index, fill a given vector with matched track indices More...

void	GetVeeTrackIndices (const size_t &ivee, vector< UInt_t > &itracks) const
	given a vee index, fill a given vector with matched G4 particle indices More...

void	GetCalClusterTrackIndices (const size_t &iclust, vector< UInt_t > &itracks) const
	given a calocluster index, fill a given vector with matched track indices More...

void	GetCalClusterG4Indices (const size_t &iclust, vector< UInt_t > &ig4ps) const
	given a calocluster index, fill a given vector with matched G4Particle indices More...

Public Member Functions inherited from garana::TreeReader
virtual	~TreeReader ()

void	SetupRead (TTree *tree)

TTree *	GetInputTree ()

size_t	NEntries () const

virtual void	GetEntry (const UInt_t &ientry)

Int_t	Event () const

const TObjArray *	GetBranchList () const

void	Fill ()

void	Write ()

void	CheckOpt (char opt)

bool	BlockWrite () const

Private Types
enum	RecoBranch { kEvent, kTracks, kVees, kVertices, kCalClusters, kVertTrackEnds, kVeeTrackEnds }

Private Member Functions
bool	VerifyBranches () override

void	FillTrack ()

void	FillVertex ()

void	FillVee ()

void	FillCaloCluster ()

Private Attributes
StructuredRecoTree *	fRecoIn = nullptr
	points to input structured recoTree we are flattening More...

std::map< std::string, RecoBranch >	nameToRecoBranch

std::map< RecoBranch, bool >	branchToDrawOpt

Additional Inherited Members
Protected Member Functions inherited from garana::FlatRecoTree
bool	SetBranchAddresses () override

void	SetVecs ()

void	ClearVecs ()

Protected Attributes inherited from garana::Mill
bool	fIsVerified = false

TTree *	fTreeOut = nullptr

TFile *	fOutfile = nullptr

Protected Attributes inherited from garana::FlatRecoTree
vector< float > *	fTrkLenFwd = nullptr
	length of the track in cm from forwards fit More...

vector< float > *	fTrkLenBac = nullptr
	length of the track in cm from backwards fit More...

vector< float > *	fTrkMomBeg = nullptr
	momentum of the track at the vertex in GeV/c More...

vector< float > *	fTrkMomEnd = nullptr
	momentum of the track at the end in GeV/c More...

vector< float > *	fTrkVtxX = nullptr
	track vertex position 4-vector cm,ns More...

vector< float > *	fTrkVtxY = nullptr

vector< float > *	fTrkVtxZ = nullptr

vector< float > *	fTrkVtxT = nullptr

vector< float > *	fTrkEndX = nullptr
	track end position 4-vector cm,ns More...

vector< float > *	fTrkEndY = nullptr

vector< float > *	fTrkEndZ = nullptr

vector< float > *	fTrkEndT = nullptr

vector< float > *	fTrkVtxDirX = nullptr
	track vertex direction More...

vector< float > *	fTrkVtxDirY = nullptr

vector< float > *	fTrkVtxDirZ = nullptr

vector< float > *	fTrkEndDirX = nullptr
	track end direction More...

vector< float > *	fTrkEndDirY = nullptr

vector< float > *	fTrkEndDirZ = nullptr

vector< float > *	fTrkChiFwd = nullptr
	chisquared forward fit More...

vector< float > *	fTrkChiBac = nullptr
	chisquared backward fit More...

vector< size_t > *	fTrkNHits = nullptr
	number of hits More...

vector< float > *	fTrackXBeg = nullptr
	Track x-coord at beginning of track [cm]. More...

vector< float > *	fTrackYBeg = nullptr
	Track y-coord (par[0]) at beginning of track [cm]. More...

vector< float > *	fTrackZBeg = nullptr
	Track z-coord (par[1]) at beginning of track [cm]. More...

vector< float > *	fTrackCurvBeg = nullptr
	Track curvature (par[2]) at beginning of track [cm^-1]. More...

vector< float > *	fTrackPhiBeg = nullptr
	Track phi (par[3]) at beginning of track [rad]. More...

vector< float > *	fTrackLambBeg = nullptr
	Track lambda (par[4]) at beginning of track [rad]. More...

vector< float > *	fTrackXEnd = nullptr
	Track x-coord at end of track [cm]. More...

vector< float > *	fTrackYEnd = nullptr
	Track y-coord (par[0]) at end of track [cm]. More...

vector< float > *	fTrackZEnd = nullptr
	Track z-coord (par[1]) at end of track [cm]. More...

vector< float > *	fTrackCurvEnd = nullptr
	Track curvature (par[2]) at end of track [cm^-1]. More...

vector< float > *	fTrackPhiEnd = nullptr
	Track phi (par[3]) at end of track [rad]. More...

vector< float > *	fTrackLambEnd = nullptr
	Track lambda (par[4]) at end of track [rad]. More...

vector< float > *	fTrkCovMat1Beg = nullptr
	covariance matrix at beginning of track – assuming symmetry More...

vector< float > *	fTrkCovMat2Beg = nullptr

vector< float > *	fTrkCovMat3Beg = nullptr

vector< float > *	fTrkCovMat4Beg = nullptr

vector< float > *	fTrkCovMat5Beg = nullptr

vector< float > *	fTrkCovMat6Beg = nullptr

vector< float > *	fTrkCovMat7Beg = nullptr

vector< float > *	fTrkCovMat8Beg = nullptr

vector< float > *	fTrkCovMat9Beg = nullptr

vector< float > *	fTrkCovMat10Beg = nullptr

vector< float > *	fTrkCovMat11Beg = nullptr

vector< float > *	fTrkCovMat12Beg = nullptr

vector< float > *	fTrkCovMat13Beg = nullptr

vector< float > *	fTrkCovMat14Beg = nullptr

vector< float > *	fTrkCovMat15Beg = nullptr

vector< float > *	fTrkCovMat1End = nullptr
	covariance matrix at end of track – assuming symmetry More...

vector< float > *	fTrkCovMat2End = nullptr

vector< float > *	fTrkCovMat3End = nullptr

vector< float > *	fTrkCovMat4End = nullptr

vector< float > *	fTrkCovMat5End = nullptr

vector< float > *	fTrkCovMat6End = nullptr

vector< float > *	fTrkCovMat7End = nullptr

vector< float > *	fTrkCovMat8End = nullptr

vector< float > *	fTrkCovMat9End = nullptr

vector< float > *	fTrkCovMat10End = nullptr

vector< float > *	fTrkCovMat11End = nullptr

vector< float > *	fTrkCovMat12End = nullptr

vector< float > *	fTrkCovMat13End = nullptr

vector< float > *	fTrkCovMat14End = nullptr

vector< float > *	fTrkCovMat15End = nullptr

vector< float > *	fTrkIonizFwd = nullptr

vector< float > *	fTrkIonizBac = nullptr

vector< int > *	fTrkChgFwd = nullptr

vector< int > *	fTrkChgBac = nullptr

vector< float > *	fVeeVtxX = nullptr
	x-position of vertex [cm] More...

vector< float > *	fVeeVtxY = nullptr
	y-position of vertex [cm] More...

vector< float > *	fVeeVtxZ = nullptr
	z-position of vertex [cm] More...

vector< float > *	fVeeVtxT = nullptr
	time of vertex [ns] More...

vector< float > *	fVeeVtxCovar11 = nullptr
	uncertainties in vertex position [cm] More...

vector< float > *	fVeeVtxCovar12 = nullptr

vector< float > *	fVeeVtxCovar13 = nullptr

vector< float > *	fVeeVtxCovar21 = nullptr

vector< float > *	fVeeVtxCovar22 = nullptr

vector< float > *	fVeeVtxCovar23 = nullptr

vector< float > *	fVeeVtxCovar31 = nullptr

vector< float > *	fVeeVtxCovar32 = nullptr

vector< float > *	fVeeVtxCovar33 = nullptr

vector< float > *	fVeePxKshort = nullptr
	four-momentum of vee, for Kshort hypothesis [GeV] More...

vector< float > *	fVeePyKshort = nullptr

vector< float > *	fVeePzKshort = nullptr

vector< float > *	fVeeEKshort = nullptr

vector< float > *	fVeePxLamb1 = nullptr
	four-momentum of vee, for Kshort hypothesis [GeV] More...

vector< float > *	fVeePyLamb1 = nullptr

vector< float > *	fVeePzLamb1 = nullptr

vector< float > *	fVeeELamb1 = nullptr

vector< float > *	fVeePxLamb2 = nullptr
	four-momentum of vee, for Kshort hypothesis [GeV] More...

vector< float > *	fVeePyLamb2 = nullptr

vector< float > *	fVeePzLamb2 = nullptr

vector< float > *	fVeeELamb2 = nullptr

vector< float > *	fVeeChisqr = nullptr
	chi-squared More...

vector< float > *	fVtxX = nullptr
	x-position of vertex [cm] More...

vector< float > *	fVtxY = nullptr
	y-position of vertex [cm] More...

vector< float > *	fVtxZ = nullptr
	z-position of vertex [cm] More...

vector< float > *	fVtxT = nullptr
	time of vertex [ns] More...

vector< float > *	fVtxCovar11 = nullptr
	uncertainties in vertex position [cm] More...

vector< float > *	fVtxCovar12 = nullptr

vector< float > *	fVtxCovar13 = nullptr

vector< float > *	fVtxCovar21 = nullptr

vector< float > *	fVtxCovar22 = nullptr

vector< float > *	fVtxCovar23 = nullptr

vector< float > *	fVtxCovar31 = nullptr

vector< float > *	fVtxCovar32 = nullptr

vector< float > *	fVtxCovar33 = nullptr

vector< float > *	fCalClustX = nullptr
	x-position of the cluster [cm] More...

vector< float > *	fCalClustY = nullptr
	y-position of the cluster [cm] More...

vector< float > *	fCalClustZ = nullptr
	z-position of the cluster [cm] More...

vector< float > *	fCalClustT = nullptr
	time of the cluster [ns] More...

vector< float > *	fCalClustEnergy = nullptr
	energy of the ECal cluster [GeV] More...

vector< float > *	fCalClustEnergyError = nullptr
	uncertainty of the cluster energy [GeV] More...

vector< float > *	fCalClustTdiff = nullptr
	time difference between the first and last layer of the cluster [ns] More...

vector< float > *	fCalClustEllip = nullptr
	cluster shape parameter ellipsoid More...

vector< float > *	fCalClustR1 = nullptr
	cluster shape parameter r1 More...

vector< float > *	fCalClustR2 = nullptr
	cluster shape parameter r2 More...

vector< float > *	fCalClustR3 = nullptr
	cluster shape parameter r3 More...

vector< float > *	fCalClustVol = nullptr
	cluster shape parameter volume More...

vector< float > *	fCalClustWidth = nullptr
	cluster shape parameter width More...

vector< float > *	fCalClustTheta = nullptr
	intrasic direction of the cluster theta More...

vector< float > *	fCalClustPhi = nullptr
	intrasic direction of the cluster phi More...

vector< float > *	fCalClustEigen1X = nullptr
	Eigenvectors of the cluster in ascending order (principal axis with smallest inertial mass) normalised to 1. More...

vector< float > *	fCalClustEigen1Y = nullptr

vector< float > *	fCalClustEigen1Z = nullptr

vector< float > *	fCalClustEigen2X = nullptr
	Eigenvectors of the cluster in ascending order (principal axis with smallest inertial mass) normalised to 1. More...

vector< float > *	fCalClustEigen2Y = nullptr

vector< float > *	fCalClustEigen2Z = nullptr

vector< float > *	fCalClustEigen3X = nullptr
	Eigenvectors of the cluster in ascending order (principal axis with smallest inertial mass) normalised to 1. More...

vector< float > *	fCalClustEigen3Y = nullptr

vector< float > *	fCalClustEigen3Z = nullptr

vector< vector< UInt_t > > *	fTrackG4PIndices = nullptr

vector< vector< UInt_t > > *	fVertTrackIndices = nullptr

vector< vector< UInt_t > > *	fVeeTrackIndices = nullptr

vector< vector< UInt_t > > *	fCalTrackIndices = nullptr

vector< vector< Int_t > > *	fVertTrackEnds = nullptr

vector< vector< Int_t > > *	fVeeTrackEnds = nullptr

TBranch *	b_TrkLenFwd = nullptr

TBranch *	b_TrkLenBac = nullptr

TBranch *	b_TrkMomBeg = nullptr

TBranch *	b_TrkMomEnd = nullptr

TBranch *	b_TrkVtxX = nullptr

TBranch *	b_TrkVtxY = nullptr

TBranch *	b_TrkVtxZ = nullptr

TBranch *	b_TrkVtxT = nullptr

TBranch *	b_TrkEndX = nullptr

TBranch *	b_TrkEndY = nullptr

TBranch *	b_TrkEndZ = nullptr

TBranch *	b_TrkEndT = nullptr

TBranch *	b_TrkVtxDirX = nullptr

TBranch *	b_TrkVtxDirY = nullptr

TBranch *	b_TrkVtxDirZ = nullptr

TBranch *	b_TrkEndDirX = nullptr

TBranch *	b_TrkEndDirY = nullptr

TBranch *	b_TrkEndDirZ = nullptr

TBranch *	b_TrkChiFwd = nullptr

TBranch *	b_TrkChiBac = nullptr

TBranch *	b_TrkNHits = nullptr

TBranch *	b_TrackXBeg = nullptr

TBranch *	b_TrackYBeg = nullptr

TBranch *	b_TrackZBeg = nullptr

TBranch *	b_TrackCurvBeg = nullptr

TBranch *	b_TrackPhiBeg = nullptr

TBranch *	b_TrackLambBeg = nullptr

TBranch *	b_TrackParEnd = nullptr

TBranch *	b_TrackXEnd = nullptr

TBranch *	b_TrackYEnd = nullptr

TBranch *	b_TrackZEnd = nullptr

TBranch *	b_TrackCurvEnd = nullptr

TBranch *	b_TrackPhiEnd = nullptr

TBranch *	b_TrackLambEnd = nullptr

TBranch *	b_TrkCovMat1Beg = nullptr

TBranch *	b_TrkCovMat2Beg = nullptr

TBranch *	b_TrkCovMat3Beg = nullptr

TBranch *	b_TrkCovMat4Beg = nullptr

TBranch *	b_TrkCovMat5Beg = nullptr

TBranch *	b_TrkCovMat6Beg = nullptr

TBranch *	b_TrkCovMat7Beg = nullptr

TBranch *	b_TrkCovMat8Beg = nullptr

TBranch *	b_TrkCovMat9Beg = nullptr

TBranch *	b_TrkCovMat10Beg = nullptr

TBranch *	b_TrkCovMat11Beg = nullptr

TBranch *	b_TrkCovMat12Beg = nullptr

TBranch *	b_TrkCovMat13Beg = nullptr

TBranch *	b_TrkCovMat14Beg = nullptr

TBranch *	b_TrkCovMat15Beg = nullptr

TBranch *	b_TrkCovMat1End = nullptr

TBranch *	b_TrkCovMat2End = nullptr

TBranch *	b_TrkCovMat3End = nullptr

TBranch *	b_TrkCovMat4End = nullptr

TBranch *	b_TrkCovMat5End = nullptr

TBranch *	b_TrkCovMat6End = nullptr

TBranch *	b_TrkCovMat7End = nullptr

TBranch *	b_TrkCovMat8End = nullptr

TBranch *	b_TrkCovMat9End = nullptr

TBranch *	b_TrkCovMat10End = nullptr

TBranch *	b_TrkCovMat11End = nullptr

TBranch *	b_TrkCovMat12End = nullptr

TBranch *	b_TrkCovMat13End = nullptr

TBranch *	b_TrkCovMat14End = nullptr

TBranch *	b_TrkCovMat15End = nullptr

TBranch *	b_VeeVtxX = nullptr

TBranch *	b_VeeVtxY = nullptr

TBranch *	b_VeeVtxZ = nullptr

TBranch *	b_VeeVtxT = nullptr

TBranch *	b_VeeVtxCovar11 = nullptr

TBranch *	b_VeeVtxCovar12 = nullptr

TBranch *	b_VeeVtxCovar13 = nullptr

TBranch *	b_VeeVtxCovar21 = nullptr

TBranch *	b_VeeVtxCovar22 = nullptr

TBranch *	b_VeeVtxCovar23 = nullptr

TBranch *	b_VeeVtxCovar31 = nullptr

TBranch *	b_VeeVtxCovar32 = nullptr

TBranch *	b_VeeVtxCovar33 = nullptr

TBranch *	b_VeePxKshort = nullptr

TBranch *	b_VeePyKshort = nullptr

TBranch *	b_VeePzKshort = nullptr

TBranch *	b_VeeEKshort = nullptr

TBranch *	b_VeePxLamb1 = nullptr

TBranch *	b_VeePyLamb1 = nullptr

TBranch *	b_VeePzLamb1 = nullptr

TBranch *	b_VeeELamb1 = nullptr

TBranch *	b_VeePxLamb2 = nullptr

TBranch *	b_VeePyLamb2 = nullptr

TBranch *	b_VeePzLamb2 = nullptr

TBranch *	b_VeeELamb2 = nullptr

TBranch *	b_VeeChisqr = nullptr

TBranch *	b_VtxX = nullptr

TBranch *	b_VtxY = nullptr

TBranch *	b_VtxZ = nullptr

TBranch *	b_VtxT = nullptr

TBranch *	b_VtxCovar11 = nullptr

TBranch *	b_VtxCovar12 = nullptr

TBranch *	b_VtxCovar13 = nullptr

TBranch *	b_VtxCovar21 = nullptr

TBranch *	b_VtxCovar22 = nullptr

TBranch *	b_VtxCovar23 = nullptr

TBranch *	b_VtxCovar31 = nullptr

TBranch *	b_VtxCovar32 = nullptr

TBranch *	b_VtxCovar33 = nullptr

TBranch *	b_CalClustX = nullptr

TBranch *	b_CalClustY = nullptr

TBranch *	b_CalClustZ = nullptr

TBranch *	b_CalClustT = nullptr

TBranch *	b_CalClustEnergy = nullptr

TBranch *	b_CalClustEnergyError = nullptr

TBranch *	b_CalClustTdiff = nullptr

TBranch *	b_CalClustEllip = nullptr

TBranch *	b_CalClustR1 = nullptr

TBranch *	b_CalClustR2 = nullptr

TBranch *	b_CalClustR3 = nullptr

TBranch *	b_CalClustVol = nullptr

TBranch *	b_CalClustWidth = nullptr

TBranch *	b_CalClustTheta = nullptr

TBranch *	b_CalClustPhi = nullptr

TBranch *	b_CalClustEigen1X = nullptr

TBranch *	b_CalClustEigen1Y = nullptr

TBranch *	b_CalClustEigen1Z = nullptr

TBranch *	b_CalClustEigen2X = nullptr

TBranch *	b_CalClustEigen2Y = nullptr

TBranch *	b_CalClustEigen2Z = nullptr

TBranch *	b_CalClustEigen3X = nullptr

TBranch *	b_CalClustEigen3Y = nullptr

TBranch *	b_CalClustEigen3Z = nullptr

TBranch *	b_TrackG4PIndices = nullptr

TBranch *	b_VertTrackIndices = nullptr

TBranch *	b_VeeTrackIndices = nullptr

TBranch *	b_CalTrackIndices = nullptr

TBranch *	b_VertTrackEnds = nullptr

TBranch *	b_VeeTrackEnds = nullptr

TBranch *	b_TrkIonizFwd = nullptr

TBranch *	b_TrkIonizBac = nullptr

TBranch *	b_TrkChgFwd = nullptr

TBranch *	b_TrkChgBac = nullptr

Protected Attributes inherited from garana::RecoTree
vector< vector< UInt_t > > *	fTrackG4PIndices = nullptr

vector< vector< UInt_t > > *	fVertTrackIndices = nullptr

vector< vector< UInt_t > > *	fVeeTrackIndices = nullptr

vector< vector< UInt_t > > *	fCalClusterTrackIndices = nullptr

vector< vector< UInt_t > > *	fCalClusterG4Indices = nullptr

TBranch *	b_TrackG4PIndices = nullptr

TBranch *	b_VertTrackIndices = nullptr

TBranch *	b_VeeTrackIndices = nullptr

TBranch *	b_CalClusterTrackIndices = nullptr

TBranch *	b_CalClusterG4Indices = nullptr

Protected Attributes inherited from garana::TreeReader
const std::string	treename

char	fOpt = 'r'

TTree *	fTreeIn = nullptr
	pointer to the analyzed TTree or TChain More...

UInt_t	fCurrentEntry = UINT_MAX

TBranch *	b_Event = nullptr

Int_t	fEvent = -1
	event number for tree entry More...

Detailed Description

Definition at line 16 of file MillRecoTree.h.

Member Enumeration Documentation

enum garana::MillRecoTree::RecoBranch

private

Enumerator
kEvent
kTracks
kVees
kVertices
kCalClusters
kVertTrackEnds
kVeeTrackEnds

Definition at line 37 of file MillRecoTree.h.

                                 {
                         kEvent,
                         kTracks,
                         kVees,
                         kVertices,
                         kCalClusters,
                         kVertTrackEnds,
                         kVeeTrackEnds
                 };//

Constructor & Destructor Documentation

garana::MillRecoTree::MillRecoTree ( )

inline

Definition at line 20 of file MillRecoTree.h.

20 {};

MillRecoTree::MillRecoTree	(	TTree *	treeIn,
		TTree *	treeOut
	)

Definition at line 20 of file MillRecoTree.cxx.

 {
     fRecoIn = new StructuredRecoTree(treeIn);
     fTreeIn = treeOut;
     fOpt='w';
         this->VerifyBranches();
         this->SetVecs();
         this->SetBranchAddresses();
         this->MillTrees();
 
 }//constructor()

Member Function Documentation

void MillRecoTree::FillCaloCluster ( )

private

Definition at line 292 of file MillRecoTree.cxx.

                                    {
 
         for(size_t iclust=0; iclust<fRecoIn->NCalCluster(); iclust++){
                 fCalClustX->push_back(fRecoIn->CalClustPosition(iclust)->X());
                 fCalClustY->push_back(fRecoIn->CalClustPosition(iclust)->Y());
                 fCalClustZ->push_back(fRecoIn->CalClustPosition(iclust)->Z());
                 fCalClustT->push_back(fRecoIn->CalClustPosition(iclust)->T());
                 fCalClustEnergy->push_back(fRecoIn->CalClustEnergy(iclust));
                 fCalClustEnergyError->push_back(fRecoIn->CalClustEnergyError(iclust));
                 fCalClustTdiff->push_back(fRecoIn->CalClustTimeDifference(iclust));
                 fCalClustEllip->push_back(fRecoIn->CalClustShape(iclust)[0]);
                 fCalClustR1->push_back(fRecoIn->CalClustShape(iclust)[1]);
                 fCalClustR2->push_back(fRecoIn->CalClustShape(iclust)[2]);
                 fCalClustR3->push_back(fRecoIn->CalClustShape(iclust)[3]);
                 fCalClustVol->push_back(fRecoIn->CalClustShape(iclust)[4]);
                 fCalClustWidth->push_back(fRecoIn->CalClustShape(iclust)[5]);
                 fCalClustTheta->push_back(fRecoIn->CalClustTheta(iclust));
                 fCalClustPhi->push_back(fRecoIn->CalClustPhi(iclust));
                 fCalClustEigen1X->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(0).X());
                 fCalClustEigen1Y->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(0).Y());
                 fCalClustEigen1Z->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(0).Z());
                 fCalClustEigen2X->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(1).X());
                 fCalClustEigen2Y->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(1).Y());
                 fCalClustEigen2Z->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(1).Z());
                 fCalClustEigen3X->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(2).X());
                 fCalClustEigen3Y->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(2).Y());
                 fCalClustEigen3Z->push_back(fRecoIn->CalClustEigenVecs(iclust)->at(2).Z());
         }
 }

void MillRecoTree::FillTrack ( )

private

Definition at line 139 of file MillRecoTree.cxx.

                              {
         for(size_t itrk=0; itrk<fRecoIn->NTrack(); itrk++){
 
                 fTrkLenFwd      ->push_back(fRecoIn->TrackLenFwd(itrk));
                 fTrkLenBac      ->push_back(fRecoIn->TrackLenBkd(itrk));
                 fTrkMomBeg      ->push_back(fRecoIn->TrackMomBeg(itrk)->Mag());
                 fTrkMomEnd      ->push_back(fRecoIn->TrackMomEnd(itrk)->Mag());
 
                 fTrkVtxX        ->push_back(fRecoIn->TrackVertex(itrk)->X());
                 fTrkVtxY        ->push_back(fRecoIn->TrackVertex(itrk)->Y());
                 fTrkVtxZ        ->push_back(fRecoIn->TrackVertex(itrk)->Z());
                 fTrkVtxT        ->push_back(fRecoIn->TrackVertex(itrk)->T());
                 fTrkEndX        ->push_back(fRecoIn->TrackEnd(itrk)->X());
                 fTrkEndY        ->push_back(fRecoIn->TrackEnd(itrk)->Y());
                 fTrkEndZ        ->push_back(fRecoIn->TrackEnd(itrk)->Z());
                 fTrkEndT        ->push_back(fRecoIn->TrackEnd(itrk)->T());
 
                 fTrkVtxDirX     ->push_back(fRecoIn->TrackVtxDirectionX(itrk));
                 fTrkVtxDirY     ->push_back(fRecoIn->TrackVtxDirectionY(itrk));
                 fTrkVtxDirZ     ->push_back(fRecoIn->TrackVtxDirectionZ(itrk));
                 fTrkEndDirX     ->push_back(fRecoIn->TrackEndDirectionX(itrk));
                 fTrkEndDirY     ->push_back(fRecoIn->TrackEndDirectionY(itrk));
                 fTrkEndDirZ     ->push_back(fRecoIn->TrackEndDirectionZ(itrk));
 
                 fTrkChiFwd      ->push_back(fRecoIn->TrackChiSqrFwd(itrk) );
                 fTrkChiBac      ->push_back(fRecoIn->TrackChiSqrBkd(itrk));
                 fTrkNHits       ->push_back(fRecoIn->NTrackHit(itrk));
 
                 float trackpars[5], trackcovar[15];
                 fRecoIn->TrackParBeg(itrk,trackpars);
 
                 fTrackXBeg   ->push_back(fRecoIn->TrackVertex(itrk)->X());
                 fTrackYBeg   ->push_back(trackpars[0]);
                 fTrackZBeg   ->push_back(trackpars[1]);
                 fTrackCurvBeg->push_back(trackpars[2]);
                 fTrackPhiBeg ->push_back(trackpars[3]);
                 fTrackLambBeg->push_back(trackpars[4]);
 
                 fRecoIn->TrackParEnd(itrk,trackpars);
                 fTrackXEnd      ->push_back(fRecoIn->TrackEnd(itrk)->X());
                 fTrackYEnd      ->push_back(trackpars[0]);
                 fTrackZEnd      ->push_back(trackpars[1]);
                 fTrackCurvEnd   ->push_back(trackpars[2]);
                 fTrackPhiEnd    ->push_back(trackpars[3]);
                 fTrackLambEnd   ->push_back(trackpars[4]);
 
                 fRecoIn->TrackCovarBeg(itrk,trackcovar);
                 fTrkCovMat1Beg ->push_back(trackcovar[0]);
                 fTrkCovMat2Beg ->push_back(trackcovar[1]);
                 fTrkCovMat3Beg ->push_back(trackcovar[2]);
                 fTrkCovMat4Beg ->push_back(trackcovar[3]);
                 fTrkCovMat5Beg ->push_back(trackcovar[4]);
                 fTrkCovMat6Beg ->push_back(trackcovar[5]);
                 fTrkCovMat7Beg ->push_back(trackcovar[6]);
                 fTrkCovMat8Beg ->push_back(trackcovar[7]);
                 fTrkCovMat9Beg ->push_back(trackcovar[8]);
                 fTrkCovMat10Beg->push_back(trackcovar[9]);
                 fTrkCovMat11Beg->push_back(trackcovar[10]);
                 fTrkCovMat12Beg->push_back(trackcovar[11]);
                 fTrkCovMat13Beg->push_back(trackcovar[12]);
                 fTrkCovMat14Beg->push_back(trackcovar[13]);
                 fTrkCovMat15Beg->push_back(trackcovar[14]);
 
                 fRecoIn->TrackCovarEnd(itrk,trackcovar);
                 fTrkCovMat1End ->push_back(trackcovar[0]);
                 fTrkCovMat2End ->push_back(trackcovar[1]);
                 fTrkCovMat3End ->push_back(trackcovar[2]);
                 fTrkCovMat4End ->push_back(trackcovar[3]);
                 fTrkCovMat5End ->push_back(trackcovar[4]);
                 fTrkCovMat6End ->push_back(trackcovar[5]);
                 fTrkCovMat7End ->push_back(trackcovar[6]);
                 fTrkCovMat8End ->push_back(trackcovar[7]);
                 fTrkCovMat9End ->push_back(trackcovar[8]);
                 fTrkCovMat10End->push_back(trackcovar[9]);
                 fTrkCovMat11End->push_back(trackcovar[10]);
                 fTrkCovMat12End->push_back(trackcovar[11]);
                 fTrkCovMat13End->push_back(trackcovar[12]);
                 fTrkCovMat14End->push_back(trackcovar[13]);
                 fTrkCovMat15End->push_back(trackcovar[14]);
 
                 fTrkIonizFwd->push_back(fRecoIn->TrackIonizFwd(itrk));
                 fTrkIonizBac->push_back(fRecoIn->TrackIonizBkd(itrk));
                 fTrkChgFwd->push_back(fRecoIn->TrackChgFwd(itrk));
                 fTrkChgBac->push_back(fRecoIn->TrackChgBkd(itrk));
 
                 //FixMe fill associations
                 //fRecoIn->GetTrackG4PIndices(itrk,fTrackG4PIndices);
 
         }//for track entries
 }

void MillRecoTree::FillVee ( )

private

Definition at line 252 of file MillRecoTree.cxx.

                            {
 
         for(size_t ivee=0; ivee<fRecoIn->NVee(); ivee++){
 
                 fVeeVtxX->push_back(fRecoIn->VeeVertex(ivee)->X());
                 fVeeVtxY->push_back(fRecoIn->VeeVertex(ivee)->Y());
                 fVeeVtxZ->push_back(fRecoIn->VeeVertex(ivee)->Z());
                 fVeeVtxT->push_back(fRecoIn->VeeVertex(ivee)->T());
                 fVeeChisqr->push_back(fRecoIn->VeeChiSquared(ivee));
 
                 float covar[3][3];
                 fRecoIn->VeeCovariance(ivee,covar);
                 fVeeVtxCovar11->push_back(covar[0][0]);
                 fVeeVtxCovar12->push_back(covar[0][1]);
                 fVeeVtxCovar13->push_back(covar[0][2]);
                 fVeeVtxCovar21->push_back(covar[1][0]);
                 fVeeVtxCovar22->push_back(covar[1][1]);
                 fVeeVtxCovar23->push_back(covar[1][2]);
                 fVeeVtxCovar31->push_back(covar[2][0]);
                 fVeeVtxCovar32->push_back(covar[2][1]);
                 fVeeVtxCovar33->push_back(covar[2][2]);
 
                 fVeePxKshort->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(0).Px());
                 fVeePyKshort->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(0).Py());
                 fVeePzKshort->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(0).Pz());
                 fVeeEKshort ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(0).E());
 
                 fVeePxLamb1 ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(1).Px());
                 fVeePyLamb1 ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(1).Py());
                 fVeePzLamb1 ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(1).Pz());
                 fVeeELamb1  ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(1).E());
 
                 fVeePxLamb2 ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(2).Px());
                 fVeePyLamb2 ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(2).Py());
                 fVeePzLamb2 ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(2).Pz());
                 fVeeELamb2  ->push_back(fRecoIn->VeeMomentumPerHypothesis(ivee)->at(2).E());
 
         }
 }

void MillRecoTree::FillVertex ( )

private

Definition at line 230 of file MillRecoTree.cxx.

                               {
 
         for(size_t ivtx=0; ivtx<fRecoIn->NVertex(); ivtx++) {
                 fVtxX->push_back(fRecoIn->GetVertex(ivtx)->X());
                 fVtxY->push_back(fRecoIn->GetVertex(ivtx)->Y());
                 fVtxZ->push_back(fRecoIn->GetVertex(ivtx)->Z());
                 fVtxT->push_back(fRecoIn->GetVertex(ivtx)->T());
 
                 float covar[3][3];
                 fRecoIn->VertexCovariance(ivtx,covar);
                 fVtxCovar11->push_back(covar[0][0]);
                 fVtxCovar12->push_back(covar[0][1]);
                 fVtxCovar13->push_back(covar[0][2]);
                 fVtxCovar21->push_back(covar[1][0]);
                 fVtxCovar22->push_back(covar[1][1]);
                 fVtxCovar23->push_back(covar[1][2]);
                 fVtxCovar31->push_back(covar[2][0]);
                 fVtxCovar32->push_back(covar[2][1]);
                 fVtxCovar33->push_back(covar[2][2]);
         }
 }

void MillRecoTree::MillTrees ( )

overridevirtual

main method

Implements garana::Mill.

Definition at line 33 of file MillRecoTree.cxx.

                              {
 
     if(!IsVerified()){
         std::cerr << "WARNING(MillRecoTree::MillTrees): "
                           << "trying to mill trees that have not been verified"
                           << std::endl;
         return;
     }
 
     std::cout << "MillRecoTree: loop over " << fRecoIn->NEntries() << " gen entries" << std::endl;
 
     //loop over input genTree events (entries)
     for(size_t ientry=0; ientry<fRecoIn->NEntries(); ientry++){
 
         ClearVecs(); //clear previous entry's data (if any)
         fRecoIn->GetEntry(ientry);
 
         if(branchToDrawOpt[kTracks]) {
                 fEvent = fRecoIn->Event();
         }
 
                 if(branchToDrawOpt[kTracks]) {
                         FillTrack();
                 }
 
                 if(branchToDrawOpt[kVertices]) {
                         FillVertex();
                 }
 
                 if(branchToDrawOpt[kVees]) {
                         FillVee();
                 }
 
                 if(branchToDrawOpt[kCalClusters]) {
                         FillCaloCluster();
                 }
 
         fTreeIn->Fill(); //actually our fTreeOut, but it's called fTreeIn in FlatRecoTree.h, fill once per event
     }//for genTree entries
 
     fTreeIn->Write();
 
 }//Mill()

bool MillRecoTree::VerifyBranches ( )

overrideprivatevirtual

Implements garana::Mill.

Definition at line 77 of file MillRecoTree.cxx.

                                   {
 
     // get list of branches and check it matches what we expect
     const TObjArray* branches = fRecoIn->GetBranchList();
     std::cout << "got ObjArray of branches from fTreeIn" << std::endl;
     try {
         if(!branches || branches->GetEntries()==0 )
                 throw branches;
     }
     catch(TObjArray* branches){
         std::cerr << "ERROR(MillRecoTree::VerifyBranches): no branches found in passed input tree"
                           << std::endl;
         return false;
     }
 
     if(branches->GetEntriesFast() != (Int_t)nameToRecoBranch.size())
         std::cout << "WARNING(MillRecoTree::VerifyBranches): Mismatch in number of branches (expected "
                   << nameToRecoBranch.size() << " but found " << branches->GetEntriesFast()
                   << ")" << std::endl;
     else
         std::cout << "found genTree with " << branches->GetEntriesFast() << " branches" << std::endl;
 
     // loop over branches
     TIter next(branches);
     TBranch* branch = nullptr;
     while( (branch=(TBranch*)next())) {
 
         // check if branch is expected
         if(nameToRecoBranch.find(CharStarToString(branch->GetFullName()))!=nameToRecoBranch.end()) {
 
                 std::cout << "   chopping branch '" << branch->GetFullName() << "'" << std::endl;
 
                 if(nameToRecoBranch[CharStarToString(branch->GetFullName())] == kEvent)
                         branchToDrawOpt[kEvent] = true;
 
                 if(nameToRecoBranch[CharStarToString(branch->GetFullName())] == kTracks)
                         branchToDrawOpt[kTracks] = true;
 
                 if(nameToRecoBranch[CharStarToString(branch->GetFullName())] == kVertices)
                         branchToDrawOpt[kVertices] = true;
 
                 if(nameToRecoBranch[CharStarToString(branch->GetFullName())] == kVees)
                         branchToDrawOpt[kVees] = true;
 
                 if(nameToRecoBranch[CharStarToString(branch->GetFullName())] == kCalClusters)
                         branchToDrawOpt[kCalClusters] = true;
 
         }//endif known branch
 
         else{
                 std::cout << "WARNING(MillRecoTree): ignoring unknown branch '"
                                   << branch->GetFullName() << "'" << std::endl;
         }//else
 
     }//for branches
 
     fIsVerified = true;
     return true;
 
 }//VerifyBranches()

Member Data Documentation

std::map<RecoBranch,bool> garana::MillRecoTree::branchToDrawOpt

private

Initial value:

= { 
                    {kEvent        , false},
                        {kTracks       , false},
                        {kVees         , false},
                        {kVertices     , false},
                    {kCalClusters  , false},
            {kVertTrackEnds, false},
            {kVeeTrackEnds , false}
                }

Definition at line 57 of file MillRecoTree.h.

StructuredRecoTree* garana::MillRecoTree::fRecoIn = nullptr

private

points to input structured recoTree we are flattening

Definition at line 33 of file MillRecoTree.h.

std::map<std::string,RecoBranch> garana::MillRecoTree::nameToRecoBranch

private

Initial value:

=       {
                        {"Event",         kEvent,       },
                        {"Tracks",        kTracks       },
                        {"Vees",          kVees         },
                        {"Vertices",      kVertices     },
                    {"CalClusters",   kCalClusters  },
            {"VertTrackEnds", kVertTrackEnds},
            {"VeeTrackEnds ", kVeeTrackEnds }
                }

Definition at line 47 of file MillRecoTree.h.

The documentation for this class was generated from the following files:

garana/garana/Processors/MillRecoTree.h
garana/garana/Processors/MillRecoTree.cxx

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation