#include <StructuredRecoTree.h>

Inheritance diagram for garana::StructuredRecoTree:

Public Member Functions
	StructuredRecoTree (TTree *tree=0)

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 4-position for vertex with index ivertex 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 CaloCluster *	GetCalCluster (const size_t &icluster) const
	================ ECal cluster ====================== More...

const TLorentzVector *	CalClustPosition (const size_t &icluster) const override

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 Member Functions
bool	SetBranchAddresses () override

Private Attributes
vector< Track > *	fTracks = nullptr

vector< Vee > *	fVees = nullptr

vector< Vertex > *	fVertices = nullptr

vector< CaloCluster > *	fCalClusters = nullptr

vector< vector< Int_t > > *	fVertTrackEnds = nullptr

vector< vector< Int_t > > *	fVeeTrackEnds = nullptr

TBranch *	b_Tracks = nullptr

TBranch *	b_Vees = nullptr

TBranch *	b_Vertices = nullptr

TBranch *	b_CalClusters = nullptr

TBranch *	b_VertTrackEnds = nullptr

TBranch *	b_VeeTrackEnds = nullptr

Additional Inherited Members
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 18 of file StructuredRecoTree.h.

Constructor & Destructor Documentation

StructuredRecoTree::StructuredRecoTree ( TTree * tree = 0 )

Definition at line 13 of file StructuredRecoTree.cxx.

                                                   {
         SetupRead(tree); //initialize tree pointer in TreeReader instance and set branch addresses
 }

Member Function Documentation

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

overridevirtual

Implements garana::RecoTree.

Definition at line 373 of file StructuredRecoTree.cxx.

                                                                                          {
         return fCalClusters->at(icluster).EigenVecs();
 }

const float StructuredRecoTree::CalClustEnergy ( const size_t & icluster ) const

overridevirtual

reconstructed ECal cluster energy

Implements garana::RecoTree.

Definition at line 331 of file StructuredRecoTree.cxx.

                                                                             {
         return GetCalCluster(icluster)->Energy();
 }

const float StructuredRecoTree::CalClustEnergyError ( const size_t & icluster ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 335 of file StructuredRecoTree.cxx.

                                                                                  {
         return GetCalCluster(icluster)->EnergyError();
 }

const float StructuredRecoTree::CalClustMaxDeposit ( const size_t & icluster ) const

overridevirtual

maximum true deposited energy from a single MCParticle

Implements garana::RecoTree.

Definition at line 348 of file StructuredRecoTree.cxx.

                                                                               {
         if(GetCalCluster(icluster))
                 return GetCalCluster(icluster)->MaxDeposit();
         else
                 return 0.;
 }

const size_t StructuredRecoTree::CalClustNTrueTrack ( const size_t & icluster ) const

overridevirtual

number of MCParticles associated with ith cluster

Implements garana::RecoTree.

Definition at line 342 of file StructuredRecoTree.cxx.

                                                                                    {
         return GetCalCluster(icluster)->NIdes();
 }

const float StructuredRecoTree::CalClustPhi ( const size_t & icluster ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 369 of file StructuredRecoTree.cxx.

                                                                          {
         return fCalClusters->at(icluster).Phi();
 }

const TLorentzVector * StructuredRecoTree::CalClustPosition ( const size_t & icluster ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 327 of file StructuredRecoTree.cxx.

                                                                                          {
         return GetCalCluster(icluster)->Position();
 }

const float * StructuredRecoTree::CalClustShape ( const size_t & icluster ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 361 of file StructuredRecoTree.cxx.

                                                                            {
         return fCalClusters->at(icluster).Shape();
 }

const float StructuredRecoTree::CalClustTheta ( const size_t & icluster ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 365 of file StructuredRecoTree.cxx.

                                                                            {
         return fCalClusters->at(icluster).Theta();
 }

const float StructuredRecoTree::CalClustTimeDifference ( const size_t & icluster ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 357 of file StructuredRecoTree.cxx.

                                                                                     {
         return fCalClusters->at(icluster).TimeDifference();
 }

const int StructuredRecoTree::CalClustTrkIdMaxDeposit ( const size_t & icluster ) const

overridevirtual

trackID of the MCParticle depositing the most energy in ith cluster

Implements garana::RecoTree.

Definition at line 345 of file StructuredRecoTree.cxx.

                                                                                    {
         return GetCalCluster(icluster)->TrackIdMaxDep();
 }

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

overridevirtual

Implements garana::RecoTree.

Definition at line 354 of file StructuredRecoTree.cxx.

                                                                                                                     {
         return GetCalCluster(icluster)->GetTrackIdEdep(itrack);
 }

const float StructuredRecoTree::CalClustTrueEnergy ( const size_t & icluster ) const

overridevirtual

total associated true energy deposited with ith cluster

Implements garana::RecoTree.

Definition at line 339 of file StructuredRecoTree.cxx.

                                                                                    {
         return GetCalCluster(icluster)->TotalTrueEnergy();
 }

const CaloCluster * StructuredRecoTree::GetCalCluster ( const size_t & icluster ) const

================ ECal cluster ======================

Definition at line 321 of file StructuredRecoTree.cxx.

                                                                                  {
         if(fCalClusters->size()==0)
                 return nullptr;
 
         return &(fCalClusters->at(icluster));
 }

const TLorentzVector * StructuredRecoTree::GetVertex ( const size_t & ivertex ) const

overridevirtual

vertex 4-position for vertex with index ivertex

Implements garana::RecoTree.

Definition at line 294 of file StructuredRecoTree.cxx.

                                                                                {
         return fVertices->at(ivertex).GetVertex();
 }

const size_t StructuredRecoTree::NCalCluster ( ) const

overridevirtual

number of ECal clusters in this event

Implements garana::RecoTree.

Definition at line 69 of file StructuredRecoTree.cxx.

                                                    {
         return fCalClusters->size();
 }

const size_t StructuredRecoTree::NTrack ( ) const

overridevirtual

number of tracks in this event

Implements garana::RecoTree.

Definition at line 60 of file StructuredRecoTree.cxx.

                                                    {
         return fTracks->size();
 }

const size_t StructuredRecoTree::NTrackHit ( const size_t & itrack ) const

overridevirtual

number of reconstructed hits used in track fit

Implements garana::RecoTree.

Definition at line 81 of file StructuredRecoTree.cxx.

                                                                      {
         return fTracks->at(itrack).fNHits;
 }

const size_t StructuredRecoTree::NVee ( ) const

overridevirtual

number of vees in this event

Implements garana::RecoTree.

Definition at line 66 of file StructuredRecoTree.cxx.

                                                    {
         return fVees->size();
 }

const size_t StructuredRecoTree::NVertex ( ) const

overridevirtual

number of vertices in this event

Implements garana::RecoTree.

Definition at line 63 of file StructuredRecoTree.cxx.

                                                    {
         return fVertices->size();
 }

bool StructuredRecoTree::SetBranchAddresses ( )

overrideprivatevirtual

Implements garana::TreeReader.

Definition at line 17 of file StructuredRecoTree.cxx.

                                             {
 
         cout << "setting StructuredRecoTree branch addresses" << endl;
 
     fTreeIn->SetBranchAddress("Event",            &fEvent,                   &b_Event                  );
     cout << "set Event" << endl;
     fTreeIn->SetBranchAddress("Tracks",           &fTracks          ,        &b_Tracks                 );
     cout << "set Track" << endl;
     fTreeIn->SetBranchAddress("Vees",             &fVees            ,        &b_Vees                   );
     cout << "set Vees" << endl;
     fTreeIn->SetBranchAddress("Vertices",         &fVertices        ,        &b_Vertices               );
     cout << "set Vertices" << endl;
     fTreeIn->SetBranchAddress("CalClusters",      &fCalClusters     ,        &b_CalClusters            );
     cout << "set CalClusters" << endl;
     fTreeIn->SetBranchAddress("TrackG4Indices",   &fTrackG4PIndices ,        &b_TrackG4PIndices        );
     cout << "set TrackG4Indices" << endl;
     fTreeIn->SetBranchAddress("VertTrackIndices", &fVertTrackIndices,        &b_VertTrackIndices       );
     fTreeIn->SetBranchAddress("VertTrackEnds",    &fVertTrackEnds   ,        &b_VertTrackEnds          );
     fTreeIn->SetBranchAddress("VeeTrackIndices",  &fVeeTrackIndices ,        &b_VeeTrackIndices        );
     fTreeIn->SetBranchAddress("VeeTrackEnds",     &fVeeTrackEnds    ,        &b_VeeTrackEnds           );
     fTreeIn->SetBranchAddress("CalTrackIndices",  &fCalClusterTrackIndices , &b_CalClusterTrackIndices );
     //fTreeIn->SetBranchAddress("CalTrackEnds",     &fCalTrackEnds    , &b_CalTrackEnds    );
     fTreeIn->SetBranchAddress("CalG4Indices",     &fCalClusterG4Indices ,    &b_CalClusterG4Indices    );
     /*if(fGeo->HasMuonDetector()){
             fTreeIn->SetBranchAddress("MuIDClusters", "vector<gar::rec::Cluster>",   &fMuClusters);
     }*/
 
     //TODO: add Assns for reco r&d products
     /*if(fAnaMode == "reco"){
         fRecoTree->Branch("TPCHits",       "vector<gar::rec::Hit>",        &fTPCHits);
         fRecoTree->Branch("TPCClusters",   "vector<gar::rec::TPCCluster>", &fTPCClusters);
         fRecoTree->Branch("CalHits",       "vector<gar::rec::CaloHit>",    &fCalHits);
         if(fGeo->HasMuonDetector()){
             fRecoTree->Branch("MuIDHits",       "vector<gar::rec::CaloHit>", &fMuHits);
         }
     }*/
 
     cout << "done." << endl;
 
     return true;
 }// end const'or

const int StructuredRecoTree::TrackChgBkd ( const size_t & itrack ) const

overridevirtual

charge sign of track if current hypothesis wrong

Implements garana::RecoTree.

Definition at line 148 of file StructuredRecoTree.cxx.

                                                                     {
         return fTracks->at(itrack).fChgBac;
 }

const int StructuredRecoTree::TrackChgFwd ( const size_t & itrack ) const

overridevirtual

charge sign of track if current hypothesis correct

Implements garana::RecoTree.

Definition at line 144 of file StructuredRecoTree.cxx.

                                                                     {
         return fTracks->at(itrack).fChgFwd;
 }

const int StructuredRecoTree::TrackChiSqrBkd ( const size_t & itrack ) const

overridevirtual

chi-squared of track fit in backward direction

Implements garana::RecoTree.

Definition at line 140 of file StructuredRecoTree.cxx.

                                                                        {
         return fTracks->at(itrack).fChiBac;
 }

const int StructuredRecoTree::TrackChiSqrFwd ( const size_t & itrack ) const

overridevirtual

chi-squared of track fit in forward direction

Implements garana::RecoTree.

Definition at line 136 of file StructuredRecoTree.cxx.

                                                                        {
         return fTracks->at(itrack).fChiFwd;
 }

void StructuredRecoTree::TrackCovarBeg	(	const size_t &	itrack,
		float	pars[15]
	)		const

overridevirtual

track fit covariance matrix (assume symmetry) at track's assumed start

Implements garana::RecoTree.

Definition at line 162 of file StructuredRecoTree.cxx.

                                                                                   {
         for(size_t i=0; i<15; i++)
                 covar[i] = fTracks->at(itrack).fCovMatBeg[i];
 }

void StructuredRecoTree::TrackCovarEnd	(	const size_t &	itrack,
		float	pars[15]
	)		const

overridevirtual

track fit covariance matrix (assume symmetry) at track's assumed end

Implements garana::RecoTree.

Definition at line 167 of file StructuredRecoTree.cxx.

                                                                                   {
         for(size_t i=0; i<15; i++)
                 covar[i] = fTracks->at(itrack).fCovMatEnd[i];
 }

const TLorentzVector * StructuredRecoTree::TrackEnd ( const size_t & itrack ) const

overridevirtual

4-position of track's assumed end point

Implements garana::RecoTree.

Definition at line 78 of file StructuredRecoTree.cxx.

                                                                              {
         return &(fTracks->at(itrack).fEnd);
 }

const float StructuredRecoTree::TrackEndDirectionX ( const size_t & itrack ) const

overridevirtual

x-direction cosine at track end

Implements garana::RecoTree.

Definition at line 108 of file StructuredRecoTree.cxx.

                                                                              {
         return fTracks->at(itrack).fEndDir.X();
 }

const float StructuredRecoTree::TrackEndDirectionY ( const size_t & itrack ) const

overridevirtual

y-direction cosine at track end

Implements garana::RecoTree.

Definition at line 112 of file StructuredRecoTree.cxx.

                                                                              {
         return fTracks->at(itrack).fEndDir.Y();
 }

const float StructuredRecoTree::TrackEndDirectionZ ( const size_t & itrack ) const

overridevirtual

z-direction cosine at track end

Implements garana::RecoTree.

Definition at line 116 of file StructuredRecoTree.cxx.

                                                                              {
         return fTracks->at(itrack).fEndDir.Z();
 }

const float StructuredRecoTree::TrackIonizBkd ( const size_t & itrack ) const

overridevirtual

track average ionization rate from backward fit

Implements garana::RecoTree.

Definition at line 132 of file StructuredRecoTree.cxx.

                                                                         {
         return fTracks->at(itrack).fLenBac;
 }

const float StructuredRecoTree::TrackIonizFwd ( const size_t & itrack ) const

overridevirtual

track average ionization rate from forward fit

Implements garana::RecoTree.

Definition at line 128 of file StructuredRecoTree.cxx.

                                                                         {
         return fTracks->at(itrack).fIonFwd;
 }

const float StructuredRecoTree::TrackLenBkd ( const size_t & itrack ) const

overridevirtual

track length from backward fit

Implements garana::RecoTree.

Definition at line 124 of file StructuredRecoTree.cxx.

                                                                       {
         return fTracks->at(itrack).fLenBac;
 }

const float StructuredRecoTree::TrackLenFwd ( const size_t & itrack ) const

overridevirtual

track length from forward fit

Implements garana::RecoTree.

Definition at line 120 of file StructuredRecoTree.cxx.

                                                                       {
         return fTracks->at(itrack).fLenFwd;
 }

const float StructuredRecoTree::TrackMaxDeposit ( const size_t & itrack ) const

overridevirtual

maximum true deposited energy from a single MCParticle

Implements garana::RecoTree.

Definition at line 268 of file StructuredRecoTree.cxx.

                                                                                {
         //TODO check if we actually want total EDep rather than just from the leading contributor
         auto it = std::max_element(fTracks->at(itrack).fTrueEnergy.begin(),fTracks->at(itrack).fTrueEnergy.end(),
                                       [](const std::pair<int,float>& lhs,const std::pair<int,float>& rhs) -> bool {return lhs.second < rhs.second; }
                                   );
         if(it==fTracks->at(itrack).fTrueEnergy.end()) {
                 cout << "max element not found for itrack = " << itrack << "! (this is just a warning; is it all noise hits?)" << endl;
                 //cout << "fTracks->at(itrack).fTrueEnergy.size() = " << fTracks->at(itrack).fTrueEnergy.size() << endl;
                 return 0.;
         }
         return (*it).second;
         //return fTracks->at(itrack).fTrueEnergy.at(TrackTrkIdMaxDeposit(itrack)).second;
 }

const TVector3 * StructuredRecoTree::TrackMomBeg ( const size_t & itrack ) const

overridevirtual

momentum 3-vector as track's assumed start point

Implements garana::RecoTree.

Definition at line 84 of file StructuredRecoTree.cxx.

                                                                           {
         TVector3* v = new TVector3(fTracks->at(itrack).fVtxDir);
         (*v) *= fTracks->at(itrack).fMomBeg;
     return v;
 }

const TVector3 * StructuredRecoTree::TrackMomEnd ( const size_t & itrack ) const

overridevirtual

momentum 3-vector as track's assumed end point

Implements garana::RecoTree.

Definition at line 90 of file StructuredRecoTree.cxx.

                                                                           {
         TVector3* v = new TVector3(fTracks->at(itrack).fEndDir);
         (*v) *= fTracks->at(itrack).fMomEnd;
     return v;
 }

const size_t StructuredRecoTree::TrackNTrueTrack ( const size_t & itrack ) const

overridevirtual

number of MCParticles associated with ith cluster

Implements garana::RecoTree.

Definition at line 239 of file StructuredRecoTree.cxx.

                                                                                {
         return fTracks->at(itrack).fTrueEnergy.size();
 }

void StructuredRecoTree::TrackParBeg	(	const size_t &	itrack,
		float	pars[5]
	)		const

overridevirtual

track fit parameters at the track's assumed start

Implements garana::RecoTree.

Definition at line 152 of file StructuredRecoTree.cxx.

                                                                               {
         for(size_t i=0; i<5; i++)
                 pars[i] = fTracks->at(itrack).fTrackParBeg[i];
 }

void StructuredRecoTree::TrackParEnd	(	const size_t &	itrack,
		float	pars[5]
	)		const

overridevirtual

track fit parameters at the track's assumed end

Implements garana::RecoTree.

Definition at line 157 of file StructuredRecoTree.cxx.

                                                                               {
         for(size_t i=0; i<5; i++)
                 pars[i] = fTracks->at(itrack).fTrackParEnd[i];
 }

const int StructuredRecoTree::TrackTrkIdMaxDeposit ( const size_t & itrack ) const

overridevirtual

trackID of the MCParticle depositing the most energy in ith cluster

Implements garana::RecoTree.

Definition at line 243 of file StructuredRecoTree.cxx.

                                                                                 {
         /*int imax=-1;
         float emax = -1.;
         for(size_t i=0; i<fTracks->at(itrack).fTrueEnergy.size(); i++) {
 
                 if(fTracks->at(itrack).fTrueEnergy.at(i).second > emax){
                         emax = fTracks->at(itrack).fTrueEnergy.at(i).second;
                         imax = i;
                 }
         }
 
         return fTracks->at(itrack).fTrueEnergy.at(imax).first;*/
 
         auto it = std::max_element(fTracks->at(itrack).fTrueEnergy.begin(),fTracks->at(itrack).fTrueEnergy.end(),
                               [](const std::pair<int,float>& lhs,const std::pair<int,float>& rhs) -> bool {return lhs.second < rhs.second; }
                                   );
         if(it==fTracks->at(itrack).fTrueEnergy.end()) {
                 cout << "max element not found for itrack = " << itrack << "! (this is just a warning; is it all noise hits?)" << endl;
                 //cout << "fTracks->at(itrack).fTrueEnergy.size() = " << fTracks->at(itrack).fTrueEnergy.size() << endl;
                 return -INT_MAX;
         }
         return (*it).first;
 
 }

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

overridevirtual

Implements garana::RecoTree.

Definition at line 282 of file StructuredRecoTree.cxx.

                                                                                                           {
 
         return &(fTracks->at(itrack).fTrueEnergy.at(itrue));
 
 }

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

overridevirtual

Implements garana::RecoTree.

Definition at line 288 of file StructuredRecoTree.cxx.

                                                                                                           {
 
         return &(fTracks->at(itrack).fTrueEnergy);
 }

const float StructuredRecoTree::TrackTrueEnergy ( const size_t & itrack ) const

overridevirtual

total associated true energy deposited with ith cluster

Implements garana::RecoTree.

Definition at line 230 of file StructuredRecoTree.cxx.

                                                                            {
         float etrue = 0.;
         for(auto const& e : fTracks->at(itrack).fTrueEnergy)
             etrue += e.second;
 
         //return FLT_MAX-itrack*0;
         return etrue;
 }

const TLorentzVector * StructuredRecoTree::TrackTrueMomBeg ( const size_t & itrack ) const

overridevirtual

true 4-momentum at track vertex [GeV/c,GeV]

Implements garana::RecoTree.

Definition at line 202 of file StructuredRecoTree.cxx.

                                                                                     {
         TLorentzVector outvec(0,0,0,0);
         float ptot = 0.;
         for(size_t imctrk=0; imctrk<fTracks->at(itrack).fTrueMomVtx.size(); imctrk++) {
                 TLorentzVector tmpvec = fTracks->at(itrack).fTrueMomVtx[imctrk].second;
                 tmpvec *= fTracks->at(itrack).fTrueMomVtx[imctrk].second.P();
                 outvec += tmpvec;
                 ptot += fTracks->at(itrack).fTrueMomVtx[imctrk].second.P();
         }
 
         outvec *= 1.0/ptot;
         return new TLorentzVector(outvec);
 }

const TLorentzVector * StructuredRecoTree::TrackTrueMomEnd ( const size_t & itrack ) const

overridevirtual

true 4-momentum at track end [GeV/c,GeV]

Implements garana::RecoTree.

Definition at line 216 of file StructuredRecoTree.cxx.

                                                                                     {
         TLorentzVector outvec(0,0,0,0);
         float ptot = 0.;
         for(size_t imctrk=0; imctrk<fTracks->at(itrack).fTrueMomEnd.size(); imctrk++) {
                 TLorentzVector tmpvec = fTracks->at(itrack).fTrueMomEnd[imctrk].second;
                 tmpvec *= fTracks->at(itrack).fTrueMomEnd[imctrk].second.P();
                 outvec += tmpvec;
                 ptot += fTracks->at(itrack).fTrueMomEnd[imctrk].second.P();
         }
 
         outvec *= 1.0/ptot;
         return new TLorentzVector(outvec);
 }

const TLorentzVector * StructuredRecoTree::TrackTruePosBeg ( const size_t & itrack ) const

overridevirtual

true 4-position at track vertex [cm,ns]

Implements garana::RecoTree.

Definition at line 172 of file StructuredRecoTree.cxx.

                                                                                     {
         TLorentzVector outvec(0,0,0,0);
         float ptot = 0.;
         for(size_t imctrk=0; imctrk<fTracks->at(itrack).fTruePosVtx.size(); imctrk++ ) {
                 TLorentzVector tmpvec = fTracks->at(itrack).fTruePosVtx[imctrk].second;
                 tmpvec *= fTracks->at(itrack).fTrueMomVtx[imctrk].second.P();
                 outvec += tmpvec;
                 ptot += fTracks->at(itrack).fTrueMomVtx[imctrk].second.P();
         }
 
         outvec *= 1.0/ptot;
         return new TLorentzVector(outvec);
 }

const TLorentzVector * StructuredRecoTree::TrackTruePosEnd ( const size_t & itrack ) const

overridevirtual

true 4-position at track end [cm,ns]

Implements garana::RecoTree.

Definition at line 186 of file StructuredRecoTree.cxx.

                                                                                     {
         TLorentzVector outvec(0,0,0,0);
         float ptot = 0.;
         std::cout << "loop over -- " << fTracks->at(itrack).fTruePosEnd.size() << " -- true positions" << std::endl;
         for(size_t imctrk=0; imctrk<fTracks->at(itrack).fTruePosEnd.size(); imctrk++) {
                 TLorentzVector tmpvec = fTracks->at(itrack).fTruePosEnd[imctrk].second;
                 std::cout << "  submomentum: " << fTracks->at(itrack).fTrueMomEnd[imctrk].second.P() << std::endl;
                 tmpvec *= fTracks->at(itrack).fTrueMomEnd[imctrk].second.P();
                 outvec += tmpvec;
                 ptot += fTracks->at(itrack).fTrueMomEnd[imctrk].second.P();
         }
     std::cout << "total true momentum: " << std::to_string(ptot) << std::endl;
         outvec *= 1.0/ptot;
         return new TLorentzVector(outvec);
 }

const TLorentzVector * StructuredRecoTree::TrackVertex ( const size_t & itrack ) const

overridevirtual

4-position of track's assumed start point

Implements garana::RecoTree.

Definition at line 74 of file StructuredRecoTree.cxx.

                                                                                 {
         return &(fTracks->at(itrack).fVtx);
 }

const float StructuredRecoTree::TrackVtxDirectionX ( const size_t & itrack ) const

overridevirtual

x-direction cosine at track vertex

Implements garana::RecoTree.

Definition at line 96 of file StructuredRecoTree.cxx.

                                                                              {
         return fTracks->at(itrack).fVtxDir.X();
 }

const float StructuredRecoTree::TrackVtxDirectionY ( const size_t & itrack ) const

overridevirtual

y-direction cosine at track vertex

Implements garana::RecoTree.

Definition at line 100 of file StructuredRecoTree.cxx.

                                                                              {
         return fTracks->at(itrack).fVtxDir.Y();
 }

const float StructuredRecoTree::TrackVtxDirectionZ ( const size_t & itrack ) const

overridevirtual

z-direction cosine at track vertex

Implements garana::RecoTree.

Definition at line 104 of file StructuredRecoTree.cxx.

                                                                              {
         return fTracks->at(itrack).fVtxDir.Z();
 }

const float StructuredRecoTree::VeeChiSquared ( const size_t & ivee ) const

overridevirtual

Implements garana::RecoTree.

Definition at line 316 of file StructuredRecoTree.cxx.

                                                                       {
         return fVees->at(ivee).GetChiSqr();
 }

void StructuredRecoTree::VeeCovariance	(	const size_t &	ivee,
		float	covar[][3]
	)		const

overridevirtual

given a vee index, fill given position covariance matrix

Implements garana::RecoTree.

Definition at line 308 of file StructuredRecoTree.cxx.

                                                                                  {
         return fVees->at(ivee).GetCovar(covar);
 }

const vector< TLorentzVector > * StructuredRecoTree::VeeMomentumPerHypothesis ( const size_t & ivee ) const

overridevirtual

4-momentum for vee if it matches {Kshort, Lambda1, Lambda2}

Implements garana::RecoTree.

Definition at line 312 of file StructuredRecoTree.cxx.

                                                                                                    {
         return fVees->at(ivee).GetMomentaPerHypothesis();
 }

const TLorentzVector * StructuredRecoTree::VeeVertex ( const size_t & ivee ) const

overridevirtual

=============== Vee =======================

Implements garana::RecoTree.

Definition at line 304 of file StructuredRecoTree.cxx.

                                                                             {
         return fVees->at(ivee).GetVertex();
 }

void StructuredRecoTree::VertexCovariance	(	const size_t &	ivertex,
		float	covar[][3]
	)		const

overridevirtual

given a vertex index, fill given position covariance matrix

Implements garana::RecoTree.

Definition at line 298 of file StructuredRecoTree.cxx.

                                                                                        {
         fVertices->at(ivertex).GetCovar(covar);
 }

Member Data Documentation

TBranch* garana::StructuredRecoTree::b_CalClusters = nullptr

private

Definition at line 109 of file StructuredRecoTree.h.

TBranch* garana::StructuredRecoTree::b_Tracks = nullptr

private

Definition at line 106 of file StructuredRecoTree.h.

TBranch* garana::StructuredRecoTree::b_Vees = nullptr

private

Definition at line 107 of file StructuredRecoTree.h.

TBranch* garana::StructuredRecoTree::b_VeeTrackEnds = nullptr

private

Definition at line 112 of file StructuredRecoTree.h.

TBranch* garana::StructuredRecoTree::b_Vertices = nullptr

private

Definition at line 108 of file StructuredRecoTree.h.

TBranch* garana::StructuredRecoTree::b_VertTrackEnds = nullptr

private

Definition at line 111 of file StructuredRecoTree.h.

vector<CaloCluster>* garana::StructuredRecoTree::fCalClusters = nullptr

private

Definition at line 101 of file StructuredRecoTree.h.

vector<Track>* garana::StructuredRecoTree::fTracks = nullptr

private

Definition at line 98 of file StructuredRecoTree.h.

vector<Vee>* garana::StructuredRecoTree::fVees = nullptr

private

Definition at line 99 of file StructuredRecoTree.h.

vector<vector<Int_t> >* garana::StructuredRecoTree::fVeeTrackEnds = nullptr

private

Definition at line 104 of file StructuredRecoTree.h.

vector<Vertex>* garana::StructuredRecoTree::fVertices = nullptr

private

Definition at line 100 of file StructuredRecoTree.h.

vector<vector<Int_t> >* garana::StructuredRecoTree::fVertTrackEnds = nullptr

private

Definition at line 103 of file StructuredRecoTree.h.

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

garana/garana/Accessors/StructuredRecoTree.h
garana/garana/Accessors/StructuredRecoTree.cxx

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation