#include <FlatG4Tree.h>

Inheritance diagram for garana::FlatG4Tree:

Public Member Functions
	FlatG4Tree ()

	FlatG4Tree (TTree *tree)

	FlatG4Tree (TTree *tree, char opt)

const UInt_t	NSim () const override
	number of particles More...

const UInt_t	NPoints (const UInt_t &iparticle) const override
	number of G4 steps (i.e. trajectory points) More...

const UInt_t	NRegions (const UInt_t &iparticle) const override
	number of regions traversed by particle More...

const Int_t	Region (const UInt_t &iparticle, const UInt_t &iregion) const override
	region number More...

const vector< const TLorentzVector * > *	SimMomEnter (const UInt_t &iparticle) const override
	particle 4-momentum at entry point, all regions More...

const vector< const TLorentzVector * > *	SimMomExit (const UInt_t &iparticle) const override
	particle 4-momentum at exit point, all regions More...

const vector< const TLorentzVector * > *	SimPosEnter (const UInt_t &iparticle) const override
	particle 4-position at entry point, all regions More...

const vector< const TLorentzVector * > *	SimPosExit (const UInt_t &iparticle) const override
	particle 4-position at exit point, all regions More...

const TLorentzVector *	SimMomEnter (const UInt_t &iparticle, const UInt_t &iregion) const override
	particle 4-momentum at entry point in region More...

const TLorentzVector *	SimMomExit (const UInt_t &iparticle, const UInt_t &iregion) const override
	particle 4-momentum at exit point in region More...

const TLorentzVector *	SimPosEnter (const UInt_t &iparticle, const UInt_t &iregion) const override
	particle 4-position at entry point in region More...

const TLorentzVector *	SimPosExit (const UInt_t &iparticle, const UInt_t &iregion) const override
	particle 4-position at exit point in region More...

const bool	IsPrimary (const UInt_t &iparticle) const override
	did particle come from generator? More...

const Int_t	PDG (const UInt_t &iparticle) const override
	particle PDG code More...

const int	ParentPDG (const UInt_t &iparticle) const override
	parent particle's PDG code More...

const int	ProgenitorPDG (const UInt_t &iparticle) const override
	PDG of primary that led this one. More...

const int	TrackID (const UInt_t &iparticle) const override
	G4 track ID (can be <0 if it fell below trking threshold) More...

const int	ParentTrackID (const UInt_t &iparticle) const override
	G4 track ID of parent particle. More...

const int	ProgenitorTrackID (const UInt_t &iparticle) const override
	G4 track ID of primary that led this one. More...

const Int_t	ProcessI (const UInt_t &iparticle) const override
	code for process that created this one More...

const Int_t	ProcessF (const UInt_t &iparticle) const override
	code for process that killed this one More...

void	GetEntry (const UInt_t &ientry) override

const UInt_t	NSubEntries () const

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

UInt_t const	GetTruthIndex (UInt_t iparticle) const
	index in gen tree subentry to truth match to this More...

const UInt_t	NPrimary () const

bool	HasPassedTPC (const UInt_t &iparticle) const
	did the G4Particle pass through any TPC drift volume(s)? More...

bool	HasPassedCalo (const UInt_t &iparticle) const
	did the G4Particle pass through any active ECal volume(s)? More...

bool	IsStoppedTPC (const UInt_t &iparticle) const
	did the G4Particle stop/decay in any TPC drift volume(s)? More...

bool	IsStoppedCalo (const UInt_t &iparticle) const
	did the G4Particle stop/decay in any active ECal volume(s)? More...

bool	IsContainedTPC (const UInt_t &iparticle) const
	if the G4Particle was produced in any TPC drift volume, does it remain in either drift volume? More...

bool	IsContainedCalo (const UInt_t &iparticle) const
	if the G4Particle was produced in any active ECal volume, does it remain there? More...

bool	IsCathodeCrosser (const UInt_t &iparticle) const
	did the G4Particle cross the TPC central cathode? More...

bool	IsContainedTPCEvent () const
	do all particles produced in any TPC drift volume in this event remain in either volume? More...

bool	IsContainedTPCPrimaries () const
	do all primaries produced in any TPC drift volume in this event remain in either volume? More...

bool	IsContainedCaloEvent () const
	do all particles produced in any active ECal volume in this event remain there? More...

bool	IsContainedCaloPrimaries () const
	do all primaries produced in any active ECal volume in this event remain there? More...

const TLorentzVector *	SimMomBegin (const UInt_t &iparticle) const

const TLorentzVector *	SimMomEnd (const UInt_t &iparticle) const

const TLorentzVector *	SimPosBegin (const UInt_t &iparticle) const

const TLorentzVector *	SimPosEnd (const UInt_t &iparticle) const

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

void	SetupRead (TTree *tree)

TTree *	GetInputTree ()

size_t	NEntries () const

Int_t	Event () const

const TObjArray *	GetBranchList () const

void	Fill ()

void	Write ()

void	CheckOpt (char opt)

bool	BlockWrite () const

Protected Member Functions
bool	SetBranchAddresses () override

void	SetVecs ()

void	ClearVecs ()

void	SetLimits ()

void	FillIndexMap ()

const UInt_t	LocalToGlobalIndex (const UInt_t &iparticle) const

Protected Attributes
std::map< UInt_t, std::pair< UInt_t, UInt_t > >	fLimits
	for each array for a given particle More...

std::map< UInt_t, UInt_t >	fLocalToGlobalIndex

UInt_t	fNSim = UINT_MAX
	number of G4 particles per event More...

vector< UInt_t > *	fNPts = nullptr
	number of 4-vector "snapshots" (G4 steps) More...

vector< UInt_t > *	fNRegions = nullptr
	number of regions of interest cross by the particle More...

vector< Int_t > *	fRegions = nullptr
	region IDs More...

vector< Int_t > *	fTrkID = nullptr
	particle's G4 trackID More...

vector< Int_t > *	fPDG = nullptr
	particle's PDG code More...

vector< Int_t > *	fParentPdg = nullptr
	particle parent's PDG code More...

vector< Int_t > *	fProgenitorPdg = nullptr
	FS particle from gen stage that led to this one. More...

vector< Int_t > *	fParentTrackId = nullptr
	particle's parent's trackID More...

vector< Int_t > *	fProgenitorTrackId = nullptr
	FS particle from gen stage that led to this one. More...

vector< Int_t > *	fProcessI = nullptr
	process that produced the particle More...

vector< Int_t > *	fProcessF = nullptr
	process that killed the particle More...

vector< Float_t > *	fX = nullptr
	particle's x-position in lab frame [cm] More...

vector< Float_t > *	fY = nullptr
	particle's y-position in lab frame [cm] More...

vector< Float_t > *	fZ = nullptr
	particle's z-position in lab frame [cm] More...

vector< Float_t > *	fT = nullptr
	particle's time in lab frame [ns] More...

vector< Float_t > *	fPx = nullptr
	particle's x-momentum in lab frame [GeV/c] More...

vector< Float_t > *	fPy = nullptr
	particle's y-momentum in lab frame [GeV/c] More...

vector< Float_t > *	fPz = nullptr
	particle's z-momentum in lab frame [GeV/c] More...

vector< Float_t > *	fE = nullptr
	particle's total energy in lab frame [GeV] More...

TBranch *	b_NSim = nullptr

TBranch *	b_NPts = nullptr

TBranch *	b_NRegions = nullptr

TBranch *	b_Regions = nullptr

TBranch *	b_TrkID = nullptr

TBranch *	b_PDG = nullptr

TBranch *	b_ParentPdg = nullptr

TBranch *	b_ProgenitorPdg = nullptr

TBranch *	b_ParentTrackId = nullptr

TBranch *	b_ProgenitorTrackId = nullptr

TBranch *	b_ProcessI = nullptr

TBranch *	b_ProcessF = nullptr

TBranch *	b_X = nullptr

TBranch *	b_Y = nullptr

TBranch *	b_Z = nullptr

TBranch *	b_T = nullptr

TBranch *	b_Px = nullptr

TBranch *	b_Py = nullptr

TBranch *	b_Pz = nullptr

TBranch *	b_E = nullptr

Protected Attributes inherited from garana::G4Tree
vector< UInt_t > *	fG4TruthIndex = nullptr

vector< UInt_t > *	fG4FSIndex = nullptr

TBranch *	b_G4TruthIndex = nullptr

TBranch *	b_G4FSIndex = 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 19 of file FlatG4Tree.h.

Constructor & Destructor Documentation

garana::FlatG4Tree::FlatG4Tree ( )

inline

Definition at line 23 of file FlatG4Tree.h.

23 {}

FlatG4Tree::FlatG4Tree ( TTree * tree )

Definition at line 14 of file FlatG4Tree.cxx.

 {
     SetupRead(tree);
 
 }

FlatG4Tree::FlatG4Tree	(	TTree *	tree,
		char	opt
	)

Definition at line 20 of file FlatG4Tree.cxx.

 {
         CheckOpt(opt);//sets TreeReader::fOpt
 
     if(fOpt=='r'){
         std::cout << "constructed FlatG4Tree object in read-only mode" << std::endl;
         SetupRead(tree);
     }
     else {
         std::cout << "constructed FlatG4Tree object in write mode" << std::endl;
         fTreeIn = tree;
         SetBranchAddresses();
         SetVecs();
     }
 
 }

Member Function Documentation

void FlatG4Tree::ClearVecs ( )

protected

Definition at line 130 of file FlatG4Tree.cxx.

                             {
 
          fNPts->clear();
          fTrkID->clear();
          fPDG->clear();
          fParentPdg->clear();
          fProgenitorPdg->clear();
          fParentTrackId->clear();
          fProgenitorTrackId->clear();
          fProcessI->clear();
          fProcessF->clear();
          fX->clear();
          fY->clear();
          fZ->clear();
          fT->clear();
          fPx->clear();
          fPy->clear();
          fPz->clear();
          fE->clear();
 
  }

void FlatG4Tree::FillIndexMap ( )

protected

Definition at line 300 of file FlatG4Tree.cxx.

                               {
         fLocalToGlobalIndex.clear();
         UInt_t iparticle=0;
         int trkid = -1; fTrkID->at(0);
         for(UInt_t i=0; i<NSubEntries(); i++){
                 if(trkid!=fTrkID->at(i)){
                         fLocalToGlobalIndex[iparticle] = i;
                         trkid = fTrkID->at(i);
                         iparticle++;
                 }
         }
 }

void FlatG4Tree::GetEntry ( const UInt_t & ientry )

overridevirtual

Reimplemented from garana::TreeReader.

Definition at line 38 of file FlatG4Tree.cxx.

                                               {
 
         FillIndexMap();
         if(ientry != fCurrentEntry) {
                 fCurrentEntry = ientry;
                 fTreeIn->GetEntry(ientry);
         }
 }

const bool FlatG4Tree::IsPrimary ( const UInt_t & iparticle ) const

overridevirtual

did particle come from generator?

Implements garana::G4Tree.

Definition at line 160 of file FlatG4Tree.cxx.

                                                                  {
         if(fParentPdg->at(iparticle)==INT_MAX) return true;
         else return false;
  }

const UInt_t FlatG4Tree::LocalToGlobalIndex ( const UInt_t & iparticle ) const

protected

Definition at line 313 of file FlatG4Tree.cxx.

                                                                          {
 
         return fLocalToGlobalIndex.at(iparticle);
 }

const UInt_t FlatG4Tree::NPoints ( const UInt_t & iparticle ) const

overridevirtual

number of G4 steps (i.e. trajectory points)

Implements garana::G4Tree.

Definition at line 156 of file FlatG4Tree.cxx.

                                                                {
         return fNPts->at(iparticle);
  }//

const UInt_t FlatG4Tree::NRegions ( const UInt_t & iparticle ) const

overridevirtual

number of regions traversed by particle

Implements garana::G4Tree.

Definition at line 169 of file FlatG4Tree.cxx.

                                                                {
         return fNRegions->at(iparticle);
 }

const UInt_t FlatG4Tree::NSim ( ) const

overridevirtual

number of particles

Implements garana::G4Tree.

Definition at line 152 of file FlatG4Tree.cxx.

                                      {
          return fNSim;
  }//

const UInt_t FlatG4Tree::NSubEntries ( ) const

Definition at line 296 of file FlatG4Tree.cxx.

                                            {
         return fTrkID->size();
 }

const int FlatG4Tree::ParentPDG ( const UInt_t & iparticle ) const

overridevirtual

parent particle's PDG code

Implements garana::G4Tree.

Definition at line 268 of file FlatG4Tree.cxx.

                                                                      {
         return fParentPdg->at(LocalToGlobalIndex(iparticle));
 }

const int FlatG4Tree::ParentTrackID ( const UInt_t & iparticle ) const

overridevirtual

G4 track ID of parent particle.

Implements garana::G4Tree.

Definition at line 280 of file FlatG4Tree.cxx.

                                                                      {
         return fParentTrackId->at(LocalToGlobalIndex(iparticle));
 }

const Int_t FlatG4Tree::PDG ( const UInt_t & iparticle ) const

overridevirtual

particle PDG code

Implements garana::G4Tree.

Definition at line 165 of file FlatG4Tree.cxx.

                                                                  {
          return fPDG->at(iparticle);
  }

const Int_t FlatG4Tree::ProcessF ( const UInt_t & iparticle ) const

overridevirtual

code for process that killed this one

Implements garana::G4Tree.

Definition at line 292 of file FlatG4Tree.cxx.

                                                               {
         return fProcessF->at(LocalToGlobalIndex(iparticle));
 }

const Int_t FlatG4Tree::ProcessI ( const UInt_t & iparticle ) const

overridevirtual

code for process that created this one

Implements garana::G4Tree.

Definition at line 288 of file FlatG4Tree.cxx.

                                                               {
         return fProcessI->at(LocalToGlobalIndex(iparticle));
 }

const int FlatG4Tree::ProgenitorPDG ( const UInt_t & iparticle ) const

overridevirtual

PDG of primary that led this one.

Implements garana::G4Tree.

Definition at line 272 of file FlatG4Tree.cxx.

                                                                      {
         return fProgenitorPdg->at(LocalToGlobalIndex(iparticle));
 }

const int FlatG4Tree::ProgenitorTrackID ( const UInt_t & iparticle ) const

overridevirtual

G4 track ID of primary that led this one.

Implements garana::G4Tree.

Definition at line 284 of file FlatG4Tree.cxx.

                                                                      {
         return fProgenitorTrackId->at(LocalToGlobalIndex(iparticle));
 }

const Int_t FlatG4Tree::Region	(	const UInt_t &	iparticle,
		const UInt_t &	iregion
	)		const

overridevirtual

region number

Implements garana::G4Tree.

Definition at line 173 of file FlatG4Tree.cxx.

                                                                                    {
         return fRegions->at(LocalToGlobalIndex(iparticle)+iregion);
 }

bool FlatG4Tree::SetBranchAddresses ( )

overrideprotectedvirtual

Implements garana::TreeReader.

Definition at line 47 of file FlatG4Tree.cxx.

                                    {
 
 std::cout << "FlatG4Tree SetBranchAddresses()" << std::endl;
 
         if(fOpt=='r'){
                 fTreeIn->SetBranchAddress("Event",             &fEvent,             &b_Event            );
             fTreeIn->SetBranchAddress("TruthIndex",        &fG4TruthIndex,      &b_G4TruthIndex     );
             //fTreeIn->SetBranchAddress("FSIndex",      &fG4FSIndex,    &b_G4FSIndex);
                 fTreeIn->SetBranchAddress("NSim",              &fNSim,              &b_NSim             );
                 fTreeIn->SetBranchAddress("NPts",              &fNPts,              &b_NPts             );
                 fTreeIn->SetBranchAddress("NRegions",          &fNRegions,          &b_NRegions         );
                 fTreeIn->SetBranchAddress("Region",            &fRegions,           &b_Regions          );
                 fTreeIn->SetBranchAddress("TrkID",             &fTrkID,             &b_TrkID            );
                 fTreeIn->SetBranchAddress("Pdg",               &fPDG,               &b_PDG              );
                 fTreeIn->SetBranchAddress("ParentPdg",         &fParentPdg,         &b_ParentPdg        );
                 fTreeIn->SetBranchAddress("ProgenitorPdg",     &fProgenitorPdg,     &b_ProgenitorPdg    );
                 fTreeIn->SetBranchAddress("ParentTrackId",     &fParentTrackId,     &b_ParentTrackId    );
                 fTreeIn->SetBranchAddress("ProgenitorTrackId", &fProgenitorTrackId, &b_ProgenitorTrackId);
                 fTreeIn->SetBranchAddress("ProcessI",          &fProcessI,          &b_ProcessI         );
                 fTreeIn->SetBranchAddress("ProcessF",          &fProcessF,          &b_ProcessF         );
                 fTreeIn->SetBranchAddress("X",                 &fX,                 &b_X                );
                 fTreeIn->SetBranchAddress("Y",                 &fY,                 &b_Y                );
                 fTreeIn->SetBranchAddress("Z",                 &fZ,                 &b_Z                );
                 fTreeIn->SetBranchAddress("T",                 &fT,                 &b_T                );
                 fTreeIn->SetBranchAddress("Px",                &fPx,                &b_Px               );
                 fTreeIn->SetBranchAddress("Py",                &fPy,                &b_Py               );
                 fTreeIn->SetBranchAddress("Pz",                &fPz,                &b_Pz               );
                 fTreeIn->SetBranchAddress("E",                 &fE,                 &b_E                );
 
         }
 
 
         if(fOpt=='w'){
 
                 fTreeIn->Branch("Event",             &fEvent,            "Event/I");
                 fTreeIn->Branch("NSim",              &fNSim,             "NSim/I" );
                 fTreeIn->Branch("NPts",              &fNPts                       );
                 fTreeIn->Branch("NRegions",          &fNRegions                   );
                 fTreeIn->Branch("Region",            &fRegions                    );
                 fTreeIn->Branch("TrkID",             &fTrkID                      );
                 fTreeIn->Branch("Pdg",               &fPDG                        );
                 fTreeIn->Branch("ParentPdg",         &fParentPdg                  );
                 fTreeIn->Branch("ProgenitorPdg",     &fProgenitorPdg              );
                 fTreeIn->Branch("ParentTrackId",     &fParentTrackId              );
                 fTreeIn->Branch("ProgenitorTrackId", &fProgenitorTrackId          );
                 fTreeIn->Branch("ProcessI",          &fProcessI                   );
                 fTreeIn->Branch("ProcessF",          &fProcessF                   );
                 fTreeIn->Branch("X",                 &fX                          );
                 fTreeIn->Branch("Y",                 &fY                          );
                 fTreeIn->Branch("Z",                 &fZ                          );
                 fTreeIn->Branch("T",                 &fT                          );
                 fTreeIn->Branch("Px",                &fPx                         );
                 fTreeIn->Branch("Py",                &fPy                         );
                 fTreeIn->Branch("Pz",                &fPz                         );
                 fTreeIn->Branch("E",                 &fE                          );
         }
 
         return true;
 }

void garana::FlatG4Tree::SetLimits ( )

protected

void FlatG4Tree::SetVecs ( )

protected

Definition at line 108 of file FlatG4Tree.cxx.

                           {
 
          fNPts                 = new vector<UInt_t>();
          fTrkID                = new vector<Int_t>();
          fPDG                  = new vector<Int_t>();
          fParentPdg            = new vector<Int_t>();
          fProgenitorPdg        = new vector<Int_t>();
          fParentTrackId        = new vector<Int_t>();
          fProgenitorTrackId    = new vector<Int_t>();
          fProcessI             = new vector<Int_t>();
          fProcessF             = new vector<Int_t>();
          fX                    = new vector<Float_t>();
          fY                    = new vector<Float_t>();
          fZ                    = new vector<Float_t>();
          fT                    = new vector<Float_t>();
          fPx                   = new vector<Float_t>();
          fPy                   = new vector<Float_t>();
          fPz                   = new vector<Float_t>();
          fE                    = new vector<Float_t>();
 
  }

const vector< const TLorentzVector * > * FlatG4Tree::SimMomEnter ( const UInt_t & iparticle ) const

overridevirtual

particle 4-momentum at entry point, all regions

Implements garana::G4Tree.

Definition at line 177 of file FlatG4Tree.cxx.

                                                                                            {
 
                 //output vector
                 auto v = new vector<const TLorentzVector*>();
 
                 //get index range
                 size_t ireg = LocalToGlobalIndex(iparticle);
                 size_t ireg_max = ireg+2*fNRegions->at(iparticle);
 
                 //fill for all regions
                 //entry points are every other value and precede exit points
                 for(; ireg<ireg_max; ireg+=2){
                         v->push_back(new TLorentzVector(fPx->at(ireg), fPy->at(ireg),
                                                                   fPz->at(ireg), fE->at(ireg)));
                 }
                 return v;
 }

const TLorentzVector * FlatG4Tree::SimMomEnter	(	const UInt_t &	iparticle,
		const UInt_t &	iregion
	)		const

overridevirtual

particle 4-momentum at entry point in region

Implements garana::G4Tree.

Definition at line 249 of file FlatG4Tree.cxx.

                                                                                                   {
         UInt_t index = LocalToGlobalIndex(iparticle) + iregion;
         return new TLorentzVector(fPx->at(index), fPy->at(index), fPz->at(index), fE->at(index));
 }

const vector< const TLorentzVector * > * FlatG4Tree::SimMomExit ( const UInt_t & iparticle ) const

overridevirtual

particle 4-momentum at exit point, all regions

Implements garana::G4Tree.

Definition at line 195 of file FlatG4Tree.cxx.

                                                                                            {
 
         //output vector
         auto v = new vector<const TLorentzVector*>();
 
         //get index range
         size_t ireg = LocalToGlobalIndex(iparticle)+1;
         const size_t ireg_max = ireg-1 + 2*fNRegions->at(iparticle);
 
         //fill for all regions
         //exit points are every other value and follow entry points
         for(; ireg<ireg_max; ireg+=2){
                 v->push_back(new TLorentzVector(fPx->at(ireg), fPy->at(ireg),
                                                           fPz->at(ireg), fE->at(ireg)));
         }
         return v;
 }

const TLorentzVector * FlatG4Tree::SimMomExit	(	const UInt_t &	iparticle,
		const UInt_t &	iregion
	)		const

overridevirtual

particle 4-momentum at exit point in region

Implements garana::G4Tree.

Definition at line 254 of file FlatG4Tree.cxx.

                                                                                                  {
         UInt_t index = LocalToGlobalIndex(iparticle) + 1 + iregion;
         return new TLorentzVector(fPx->at(index), fPy->at(index), fPz->at(index), fE->at(index));
 }

const vector< const TLorentzVector * > * FlatG4Tree::SimPosEnter ( const UInt_t & iparticle ) const

overridevirtual

particle 4-position at entry point, all regions

Implements garana::G4Tree.

Definition at line 213 of file FlatG4Tree.cxx.

                                                                                            {
 
         //output vector
         auto v = new vector<const TLorentzVector*>();
 
         //get index range
         size_t ireg = LocalToGlobalIndex(iparticle);
         size_t ireg_max = ireg+2*fNRegions->at(iparticle);
 
         //fill for all regions
         //exit points are every other value and follow entry points
         for(; ireg<ireg_max; ireg+=2){
                 v->push_back(new TLorentzVector(fX->at(ireg), fY->at(ireg),
                                                           fZ->at(ireg), fT->at(ireg)));
         }
         return v;
 }

const TLorentzVector * FlatG4Tree::SimPosEnter	(	const UInt_t &	iparticle,
		const UInt_t &	iregion
	)		const

overridevirtual

particle 4-position at entry point in region

Implements garana::G4Tree.

Definition at line 259 of file FlatG4Tree.cxx.

                                                                                                   {
         UInt_t index = LocalToGlobalIndex(iparticle) + iregion;
         return new TLorentzVector(fX->at(index), fY->at(index), fZ->at(index), fT->at(index));
 }

const vector< const TLorentzVector * > * FlatG4Tree::SimPosExit ( const UInt_t & iparticle ) const

overridevirtual

particle 4-position at exit point, all regions

Implements garana::G4Tree.

Definition at line 231 of file FlatG4Tree.cxx.

                                                                                            {
 
         //output vector
         auto v = new vector<const TLorentzVector*>();
 
         //get index range
         size_t ireg = LocalToGlobalIndex(iparticle)+1;
         size_t ireg_max = ireg-1 + 2*fNRegions->at(iparticle);
 
         //fill for all regions
         //exit points are every other value and follow entry points
         for(; ireg<ireg_max; ireg+=2){
                 v->push_back(new TLorentzVector(fX->at(ireg), fY->at(ireg),
                                                           fZ->at(ireg), fT->at(ireg)));
         }
         return v;
 }

const TLorentzVector * FlatG4Tree::SimPosExit	(	const UInt_t &	iparticle,
		const UInt_t &	iregion
	)		const

overridevirtual

particle 4-position at exit point in region

Implements garana::G4Tree.

Definition at line 263 of file FlatG4Tree.cxx.

                                                                                                  {
         UInt_t index = LocalToGlobalIndex(iparticle) + 1 + iregion;
         return new TLorentzVector(fX->at(index), fY->at(index), fZ->at(index), fT->at(index));
 }

const int FlatG4Tree::TrackID ( const UInt_t & iparticle ) const

overridevirtual

G4 track ID (can be <0 if it fell below trking threshold)

Implements garana::G4Tree.

Definition at line 276 of file FlatG4Tree.cxx.

                                                                      {
         return fTrkID->at(LocalToGlobalIndex(iparticle));
 }

Member Data Documentation

TBranch* garana::FlatG4Tree::b_E = nullptr

protected

Definition at line 113 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_NPts = nullptr

protected

Definition at line 95 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_NRegions = nullptr

protected

Definition at line 96 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_NSim = nullptr

protected

Definition at line 94 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_ParentPdg = nullptr

protected

Definition at line 100 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_ParentTrackId = nullptr

protected

Definition at line 102 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_PDG = nullptr

protected

Definition at line 99 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_ProcessF = nullptr

protected

Definition at line 105 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_ProcessI = nullptr

protected

Definition at line 104 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_ProgenitorPdg = nullptr

protected

Definition at line 101 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_ProgenitorTrackId = nullptr

protected

Definition at line 103 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_Px = nullptr

protected

Definition at line 110 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_Py = nullptr

protected

Definition at line 111 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_Pz = nullptr

protected

Definition at line 112 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_Regions = nullptr

protected

Definition at line 97 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_T = nullptr

protected

Definition at line 109 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_TrkID = nullptr

protected

Definition at line 98 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_X = nullptr

protected

Definition at line 106 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_Y = nullptr

protected

Definition at line 107 of file FlatG4Tree.h.

TBranch* garana::FlatG4Tree::b_Z = nullptr

protected

Definition at line 108 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fE = nullptr

protected

particle's total energy in lab frame [GeV]

Definition at line 91 of file FlatG4Tree.h.

std::map<UInt_t, std::pair<UInt_t,UInt_t> > garana::FlatG4Tree::fLimits

protected

for each array for a given particle

tree is completely flat so need first and last index

Definition at line 64 of file FlatG4Tree.h.

std::map<UInt_t,UInt_t> garana::FlatG4Tree::fLocalToGlobalIndex

protected

Definition at line 66 of file FlatG4Tree.h.

vector<UInt_t>* garana::FlatG4Tree::fNPts = nullptr

protected

number of 4-vector "snapshots" (G4 steps)

Definition at line 73 of file FlatG4Tree.h.

vector<UInt_t>* garana::FlatG4Tree::fNRegions = nullptr

protected

number of regions of interest cross by the particle

Definition at line 74 of file FlatG4Tree.h.

UInt_t garana::FlatG4Tree::fNSim = UINT_MAX

protected

number of G4 particles per event

Definition at line 72 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fParentPdg = nullptr

protected

particle parent's PDG code

Definition at line 78 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fParentTrackId = nullptr

protected

particle's parent's trackID

Definition at line 80 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fPDG = nullptr

protected

particle's PDG code

Definition at line 77 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fProcessF = nullptr

protected

process that killed the particle

Definition at line 83 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fProcessI = nullptr

protected

process that produced the particle

Definition at line 82 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fProgenitorPdg = nullptr

protected

FS particle from gen stage that led to this one.

Definition at line 79 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fProgenitorTrackId = nullptr

protected

FS particle from gen stage that led to this one.

Definition at line 81 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fPx = nullptr

protected

particle's x-momentum in lab frame [GeV/c]

Definition at line 88 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fPy = nullptr

protected

particle's y-momentum in lab frame [GeV/c]

Definition at line 89 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fPz = nullptr

protected

particle's z-momentum in lab frame [GeV/c]

Definition at line 90 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fRegions = nullptr

protected

region IDs

Definition at line 75 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fT = nullptr

protected

particle's time in lab frame [ns]

Definition at line 87 of file FlatG4Tree.h.

vector<Int_t>* garana::FlatG4Tree::fTrkID = nullptr

protected

particle's G4 trackID

Definition at line 76 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fX = nullptr

protected

particle's x-position in lab frame [cm]

Definition at line 84 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fY = nullptr

protected

particle's y-position in lab frame [cm]

Definition at line 85 of file FlatG4Tree.h.

vector<Float_t>* garana::FlatG4Tree::fZ = nullptr

protected

particle's z-position in lab frame [cm]

Definition at line 86 of file FlatG4Tree.h.

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

garana/garana/Accessors/FlatG4Tree.h
garana/garana/Accessors/FlatG4Tree.cxx

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation