#include <MillGenTree.h>

Inheritance diagram for garana::MillGenTree:

Public Member Functions
	MillGenTree ()

	MillGenTree (TTree treeIn, TTree treeOut)

void	MillTrees () override

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

bool	IsVerified () const

std::string	CharStarToString (const char *cstr)

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

	FlatGenTree (TTree *tree)

	FlatGenTree (TTree *tree, char opt)

const UInt_t	NGen () const override

const UInt_t	NFSParticles (const UInt_t &igen) const override

const Bool_t	IsGenie (const UInt_t &igen) const override

const int	NuRegion (const UInt_t &igen) const override
	region code where the vertex is located More...

const Bool_t	IsCC (const UInt_t &igen) const override
	whether interaction is CC or not More...

const Int_t	NuPDG (const UInt_t &igen) const override
	neutrino PDG code More...

const TLorentzVector *	NuP (const UInt_t &igen) override
	initial neutrino 4-momentum More...

const TLorentzVector *	NuVertex (const UInt_t &igen) override
	4-position of neutrino vertex More...

const int	ScatterCode (const UInt_t &igen) const override
	GENIE neutrino scattering code. More...

const int	InteractCode (const UInt_t &igen) const override
	GENIE interaction code. More...

const Int_t	TgtPDG (const UInt_t &igen) const override
	PDG of Target Nucleus, nucleon only if free. More...

const TLorentzVector *	TgtP4 (const UInt_t &igen) const override
	target 4-momentum More...

const int	TgtZ (const UInt_t &igen) const override
	target Z More...

const int	TgtA (const UInt_t &igen) const override
	target A More...

const int	HitNucPDG (const UInt_t &igen) const override
	hit nucleon PDG code More...

const TLorentzVector *	HitNucP4 (const UInt_t &igen) const override
	hit nucleon (initial) 4-momentum More...

const double	Q2 (const UInt_t &igen) const override

const double	q2 (const UInt_t &igen) const override

const double	W (const UInt_t &igen) const override

const double	T (const UInt_t &igen) const override

const double	X (const UInt_t &igen) const override

const double	Y (const UInt_t &igen) const override

const TLorentzVector *	FSLeptonP4 (const UInt_t &igen) const override
	final state primary lepton 4-momentum (LAB frame) More...

const TLorentzVector *	FSHadSystP4 (const UInt_t &igen) const override
	final state hadronic system 4-momentum(LAB frame) More...

const int	NumNuProton (const UInt_t &igen) const override
	number of protons after reaction, before FSI More...

const int	NumNuNeutron (const UInt_t &igen) const override
	number of neutrons after reaction, before FSI More...

const int	NumNuPi0 (const UInt_t &igen) const override
	number of pi0 after reaction, before FSI More...

const int	NumNuPiPlus (const UInt_t &igen) const override
	number of pi pluses after reaction, before FSI More...

const int	NumNuPiMinus (const UInt_t &igen) const override
	number of pi minuses after reaction, before FSI More...

const Float_t	FSTotEnergy (const UInt_t &igen) const override

const Float_t	FSEnergy (const UInt_t &igen, const UInt_t &ifsp) const override

const Float_t	FSPDG (const UInt_t &igen, const UInt_t &ifsp) const override

const Int_t	FSTrackId (const UInt_t &igen, const UInt_t &ifsp) const override

void	GetEntry (const UInt_t &ientry) override

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

const std::vector< Int_t > *	GetGIndex () const
	whether or not subentry is GENIE interaction, index of GENIE interaction More...

bool	NuInAV (const UInt_t &igen) const

bool	NuInIV (const UInt_t &igen) const

bool	NuInCalo (const UInt_t &igen) 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

Private Types
enum	GenBranch { kEvent, kGIndex, kGTruth, kFSParticle }

Private Member Functions
bool	VerifyBranches () override

void	FillGTruth (const GTruth *truth)

void	FillEmptyGTruth ()

void	FillFSParticle (const vector< FSParticle > *fsp)

Private Attributes
StructuredGenTree *	fGenIn = nullptr

std::map< std::string, GenBranch >	nameToGenBranch

std::map< GenBranch, bool >	branchToDrawOpt

Additional Inherited Members
Protected Types inherited from garana::FlatGenTree
typedef UInt_t	fsLimit

typedef fsLimit	fsBegin

typedef fsLimit	fsEnd

Protected Member Functions inherited from garana::FlatGenTree
bool	SetBranchAddresses () override

void	FindFSLimits ()

std::pair< UInt_t, UInt_t >	FSParticlesFromGTruth (size_t igen) const

void	SetVecs ()

void	ClearVecs ()

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

TTree *	fTreeOut = nullptr

TFile *	fOutfile = nullptr

Protected Attributes inherited from garana::FlatGenTree
vector< std::pair< fsBegin, fsEnd > >	fGenToFSLimits

TLorentzVector	fNuP

TLorentzVector	fNuVertex

Int_t	fNGen = -1

const vector< Int_t > *	fGIndex = nullptr

vector< UInt_t >	fNFSObj
	number of FSParticles for igen^th GTruth len=NGen More...

vector< UInt_t >	fFSIndexObj

vector< UInt_t >	fFSPdgObj
	particle PDG code More...

vector< Float_t >	fFSPosXObj
	x-coordinate [cm] More...

vector< Float_t >	fFSPosYObj
	y-coordinate [cm] More...

vector< Float_t >	fFSPosZObj
	z-coordinate [cm] More...

vector< Float_t >	fFSTObj
	time [ns] More...

vector< Float_t >	fFSMomXObj
	momentum, x-component [GeV/c] More...

vector< Float_t >	fFSMomYObj
	momentum, y-component [GeV/c] More...

vector< Float_t >	fFSMomZObj
	momentum, z-component [GeV/c] More...

vector< Float_t >	fFSEObj
	total energy [GeV] More...

vector< double >	fNuXObj
	GTruth (one or more per genTree entry) //. More...

vector< double >	fNuYObj

vector< double >	fNuZObj

vector< double >	fNuTObj

vector< double >	fWeightObj

vector< double >	fProbabilityObj

vector< double >	fXsecObj

vector< double >	fDiffXsecObj

vector< int >	fGPhaseSpaceObj

vector< int >	fProbePDGObj

vector< double >	fProbePxObj

vector< double >	fProbePyObj

vector< double >	fProbePzObj

vector< double >	fProbeEObj

vector< double >	fTgtPxObj

vector< double >	fTgtPyObj

vector< double >	fTgtPzObj

vector< double >	fTgtEObj

vector< int >	fTgtZObj

vector< int >	fTgtAObj

vector< int >	fTgtPDGObj

vector< int >	fHitNucPDGObj

vector< int >	fHitQrkPDGObj

vector< bool >	fIsSeaQuarkObj

vector< double >	fHitNucPxObj

vector< double >	fHitNucPyObj

vector< double >	fHitNucPzObj

vector< double >	fHitNucEObj

vector< double >	fHitNucPosObj

vector< int >	fGscatterObj

vector< int >	fGintObj

vector< double >	fGQ2Obj

vector< double >	fGq2Obj

vector< double >	fGWObj

vector< double >	fGTObj

vector< double >	fGXObj

vector< double >	fGYObj

vector< double >	fFSleptonPxObj

vector< double >	fFSleptonPyObj

vector< double >	fFSleptonPzObj

vector< double >	fFSleptonEObj

vector< double >	fFShadSystPxObj

vector< double >	fFShadSystPyObj

vector< double >	fFShadSystPzObj

vector< double >	fFShadSystEObj

vector< bool >	fIsCharmObj

vector< int >	fCharmHadronPDGObj

vector< bool >	fIsStrangeObj

vector< int >	fStrangeHadronPDGObj

vector< int >	fNumProtonObj

vector< int >	fNumNeutronObj

vector< int >	fNumPi0Obj

vector< int >	fNumPiPlusObj

vector< int >	fNumPiMinusObj

vector< int >	fResNumObj

vector< int >	fDecayModeObj

vector< UInt_t > *	fNFS = nullptr
	number of FSParticles for igen^th GTruth len=NGen More...

vector< UInt_t > *	fFSIndex = nullptr

vector< UInt_t > *	fFSTrackId = nullptr
	particle generator track ID More...

vector< UInt_t > *	fFSPdg = nullptr
	particle PDG code More...

vector< Float_t > *	fFSPosX = nullptr
	x-coordinate [cm] More...

vector< Float_t > *	fFSPosY = nullptr
	y-coordinate [cm] More...

vector< Float_t > *	fFSPosZ = nullptr
	z-coordinate [cm] More...

vector< Float_t > *	fFST = nullptr
	time [ns] More...

vector< Float_t > *	fFSMomX = nullptr
	momentum, x-component [GeV/c] More...

vector< Float_t > *	fFSMomY = nullptr
	momentum, y-component [GeV/c] More...

vector< Float_t > *	fFSMomZ = nullptr
	momentum, z-component [GeV/c] More...

vector< Float_t > *	fFSE = nullptr
	total energy [GeV] More...

vector< double > *	fNuX = nullptr
	GTruth (one or more per genTree entry) //. More...

vector< double > *	fNuY = nullptr

vector< double > *	fNuZ = nullptr

vector< double > *	fNuT = nullptr

vector< int > *	fNuRegion = nullptr

vector< double > *	fWeight = nullptr

vector< double > *	fProbability = nullptr

vector< double > *	fXsec = nullptr

vector< double > *	fDiffXsec = nullptr

vector< int > *	fGPhaseSpace = nullptr

vector< int > *	fProbePDG = nullptr

vector< double > *	fProbePx = nullptr

vector< double > *	fProbePy = nullptr

vector< double > *	fProbePz = nullptr

vector< double > *	fProbeE = nullptr

vector< double > *	fTgtPx = nullptr

vector< double > *	fTgtPy = nullptr

vector< double > *	fTgtPz = nullptr

vector< double > *	fTgtE = nullptr

vector< int > *	fTgtZ = nullptr

vector< int > *	fTgtA = nullptr

vector< int > *	fTgtPDG = nullptr

vector< int > *	fHitNucPDG = nullptr

vector< int > *	fHitQrkPDG = nullptr

vector< bool > *	fIsSeaQuark = nullptr

vector< double > *	fHitNucPx = nullptr

vector< double > *	fHitNucPy = nullptr

vector< double > *	fHitNucPz = nullptr

vector< double > *	fHitNucE = nullptr

vector< double > *	fHitNucPos = nullptr

vector< int > *	fGscatter = nullptr

vector< int > *	fGint = nullptr

vector< double > *	fGQ2 = nullptr

vector< double > *	fGq2 = nullptr

vector< double > *	fGW = nullptr

vector< double > *	fGT = nullptr

vector< double > *	fGX = nullptr

vector< double > *	fGY = nullptr

vector< double > *	fFSleptonPx = nullptr

vector< double > *	fFSleptonPy = nullptr

vector< double > *	fFSleptonPz = nullptr

vector< double > *	fFSleptonE = nullptr

vector< double > *	fFShadSystPx = nullptr

vector< double > *	fFShadSystPy = nullptr

vector< double > *	fFShadSystPz = nullptr

vector< double > *	fFShadSystE = nullptr

vector< bool > *	fIsCharm = nullptr

vector< int > *	fCharmHadronPDG = nullptr

vector< bool > *	fIsStrange = nullptr

vector< int > *	fStrangeHadronPDG = nullptr

vector< int > *	fNumProton = nullptr

vector< int > *	fNumNeutron = nullptr

vector< int > *	fNumPi0 = nullptr

vector< int > *	fNumPiPlus = nullptr

vector< int > *	fNumPiMinus = nullptr

vector< int > *	fResNum = nullptr

vector< int > *	fDecayMode = nullptr

TBranch *	b_NGen = nullptr

TBranch *	b_GIndex = nullptr

TBranch *	b_NFS = nullptr

TBranch *	b_FSIndex = nullptr

TBranch *	b_FSPdg = nullptr

TBranch *	b_FSPosX = nullptr

TBranch *	b_FSPosY = nullptr

TBranch *	b_FSPosZ = nullptr

TBranch *	b_FST = nullptr

TBranch *	b_FSMomX = nullptr

TBranch *	b_FSMomY = nullptr

TBranch *	b_FSMomZ = nullptr

TBranch *	b_FSE = nullptr

TBranch *	b_NuX = nullptr

TBranch *	b_NuY = nullptr

TBranch *	b_NuZ = nullptr

TBranch *	b_NuT = nullptr

TBranch *	b_NuRegion = nullptr

TBranch *	b_Weight = nullptr

TBranch *	b_Probability = nullptr

TBranch *	b_Xsec = nullptr

TBranch *	b_DiffXsec = nullptr

TBranch *	b_GPhaseSpace = nullptr

TBranch *	b_ProbePDG = nullptr

TBranch *	b_ProbePx = nullptr

TBranch *	b_ProbePy = nullptr

TBranch *	b_ProbePz = nullptr

TBranch *	b_ProbeE = nullptr

TBranch *	b_TgtPx = nullptr

TBranch *	b_TgtPy = nullptr

TBranch *	b_TgtPz = nullptr

TBranch *	b_TgtE = nullptr

TBranch *	b_TgtZ = nullptr

TBranch *	b_TgtA = nullptr

TBranch *	b_TgtPDG = nullptr

TBranch *	b_HitNucPDG = nullptr

TBranch *	b_HitQrkPDG = nullptr

TBranch *	b_IsSeaQuark = nullptr

TBranch *	b_HitNucPx = nullptr

TBranch *	b_HitNucPy = nullptr

TBranch *	b_HitNucPz = nullptr

TBranch *	b_HitNucE = nullptr

TBranch *	b_HitNucPos = nullptr

TBranch *	b_Gscatter = nullptr

TBranch *	b_Gint = nullptr

TBranch *	b_GQ2 = nullptr

TBranch *	b_Gq2 = nullptr

TBranch *	b_GW = nullptr

TBranch *	b_GT = nullptr

TBranch *	b_GX = nullptr

TBranch *	b_GY = nullptr

TBranch *	b_FSleptonPx = nullptr

TBranch *	b_FSleptonPy = nullptr

TBranch *	b_FSleptonPz = nullptr

TBranch *	b_FSleptonE = nullptr

TBranch *	b_FShadSystPx = nullptr

TBranch *	b_FShadSystPy = nullptr

TBranch *	b_FShadSystPz = nullptr

TBranch *	b_FShadSystE = nullptr

TBranch *	b_IsCharm = nullptr

TBranch *	b_CharmHadronPDG = nullptr

TBranch *	b_IsStrange = nullptr

TBranch *	b_StrangeHadronPDG = nullptr

TBranch *	b_NumProton = nullptr

TBranch *	b_NumNeutron = nullptr

TBranch *	b_NumPi0 = nullptr

TBranch *	b_NumPiPlus = nullptr

TBranch *	b_NumPiMinus = nullptr

TBranch *	b_ResNum = nullptr

TBranch *	b_DecayMode = nullptr

Protected Attributes inherited from garana::GenTree
const std::vector< Int_t > *	fGIndex = nullptr

TBranch *	b_GIndex = 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 MillGenTree.h.

Member Enumeration Documentation

enum garana::MillGenTree::GenBranch

private

Enumerator
kEvent
kGIndex
kGTruth
kFSParticle

Definition at line 37 of file MillGenTree.h.

                                {
                         kEvent,
                         kGIndex,
                         kGTruth,
                         kFSParticle
                 };//

Constructor & Destructor Documentation

garana::MillGenTree::MillGenTree ( )

inline

Definition at line 20 of file MillGenTree.h.

20 {};

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

Definition at line 26 of file MillGenTree.cxx.

 {
     fGenIn = new StructuredGenTree(treeIn);
 
     fTreeIn = treeOut;
     //fTreeOut = new TTree("genTree", "generator level information");
     //fGenOut= new FlatGenTree(treeOut,'w');
 
     fOpt='w';
         this->VerifyBranches();
         this->SetVecs();
         this->SetBranchAddresses();
         this->MillTrees();
 
 }//constructor()

Member Function Documentation

void MillGenTree::FillEmptyGTruth ( )

private

Definition at line 228 of file MillGenTree.cxx.

                                  {
 
         FillEmpty(fNuX);
         FillEmpty(fNuY);
         FillEmpty(fNuZ);
         FillEmpty(fNuT);
         FillEmpty(fWeight);
         FillEmpty(fProbability);
         FillEmpty(fXsec);
         FillEmpty(fDiffXsec);
         FillEmpty(fGPhaseSpace);
         FillEmpty(fProbePDG);
         FillEmpty(fProbePx);
         FillEmpty(fProbePy);
         FillEmpty(fProbePz);
         FillEmpty(fProbeE);
         FillEmpty(fTgtPx);
         FillEmpty(fTgtPy);
         FillEmpty(fTgtPz);
         FillEmpty(fTgtE);
 
         FillEmpty(fTgtZ);
         FillEmpty(fTgtA);
         FillEmpty(fTgtPDG);
         FillEmpty(fHitNucPDG);
         FillEmpty(fHitQrkPDG);
         FillEmpty(fIsSeaQuark);
         FillEmpty(fHitNucPx);
         FillEmpty(fHitNucPy);
         FillEmpty(fHitNucPz);
         FillEmpty(fHitNucE);
         FillEmpty(fHitNucPos);
 
         FillEmpty(fGscatter);
         FillEmpty(fGint);
         FillEmpty(fGQ2);
         FillEmpty(fGq2);
         FillEmpty(fGW);
         FillEmpty(fGT);
         FillEmpty(fGX);
         FillEmpty(fGY);
         FillEmpty(fFSleptonPx);
         FillEmpty(fFSleptonPy);
         FillEmpty(fFSleptonPz);
         FillEmpty(fFSleptonE);
         FillEmpty(fFShadSystPx);
         FillEmpty(fFShadSystPy);
         FillEmpty(fFShadSystPz);
         FillEmpty(fFShadSystE);
 
         FillEmpty(fIsCharm);
         FillEmpty(fCharmHadronPDG);
         FillEmpty(fIsStrange);
         FillEmpty(fStrangeHadronPDG);
         FillEmpty(fNumProton);
         FillEmpty(fNumNeutron);
         FillEmpty(fNumPi0);
         FillEmpty(fNumPiPlus);
         FillEmpty(fNumPiMinus);
         FillEmpty(fResNum);
         FillEmpty(fDecayMode);
 }//FillEmptyGTruth

void MillGenTree::FillFSParticle ( const vector< FSParticle > * fsp )

private

Definition at line 291 of file MillGenTree.cxx.

                                                               {
 
         for( auto const& fsp : *fsps ) {
 
                 fFSPdg->push_back(fsp.PDG());
                 fFSPosX->push_back(fsp.X());
                 fFSPosY->push_back(fsp.Y());
                 fFSPosZ->push_back(fsp.Z());
                 fFST->push_back(fsp.T());
                 fFSMomX->push_back(fsp.Px());
                 fFSMomY->push_back(fsp.Py());
                 fFSMomZ->push_back(fsp.Pz());
                 fFSE->push_back(fsp.E());
 
         }//for FSParticles
 }//FillFSParticle

void MillGenTree::FillGTruth ( const GTruth * truth )

private

Definition at line 161 of file MillGenTree.cxx.

                                                {
 
         fNuX->push_back(truth->fVertex.X());
         fNuY->push_back(truth->fVertex.Y());
         fNuZ->push_back(truth->fVertex.Z());
         fNuT->push_back(truth->fVertex.T());
 
         fWeight->push_back(truth->fweight);
         fProbability->push_back(truth->fprobability);
         fXsec->push_back(truth->fXsec);
         fDiffXsec->push_back(truth->fDiffXsec);
         fGPhaseSpace->push_back(truth->fGPhaseSpace);
 
         fProbePDG->push_back(truth->fProbePDG);
         fProbePx->push_back(truth->fProbeP4.Px());
         fProbePy->push_back(truth->fProbeP4.Py());
         fProbePz->push_back(truth->fProbeP4.Pz());
         fProbeE->push_back(truth->fProbeP4.E());
 
         fTgtPx->push_back(truth->fProbeP4.Px());
         fTgtPy->push_back(truth->fProbeP4.Py());
         fTgtPz->push_back(truth->fProbeP4.Pz());
         fTgtE->push_back(truth->fProbeP4.E());
         fTgtZ->push_back(truth->ftgtZ);
         fTgtA->push_back(truth->ftgtA);
         fTgtPDG->push_back(truth->ftgtPDG);
         fHitNucPDG->push_back(truth->fHitNucPDG);
         fHitQrkPDG->push_back(truth->fHitQrkPDG);
         fIsSeaQuark->push_back(truth->fIsSeaQuark);
 
         fHitNucPx->push_back(truth->fHitNucP4.Px());
         fHitNucPy->push_back(truth->fHitNucP4.Py());
         fHitNucPz->push_back(truth->fHitNucP4.Pz());
         fHitNucE->push_back(truth->fHitNucP4.E());
         fHitNucPos->push_back(truth->fHitNucPos);
         fGscatter->push_back(truth->fGscatter);
         fGint->push_back(truth->fGint);
         fGQ2->push_back(truth->fgQ2);
         fGq2->push_back(truth->fgq2);
         fGW->push_back(truth->fgW);
         fGT->push_back(truth->fgT);
         fGX->push_back(truth->fgX);
         fGY->push_back(truth->fgY);
 
         fFSleptonPx->push_back(truth->fFSleptonP4.Px());
         fFSleptonPy->push_back(truth->fFSleptonP4.Py());
         fFSleptonPz->push_back(truth->fFSleptonP4.Pz());
         fFSleptonE->push_back(truth->fFSleptonP4.E());
 
         fFShadSystPx->push_back(truth->fFShadSystP4.Px());
         fFShadSystPy->push_back(truth->fFShadSystP4.Py());
         fFShadSystPz->push_back(truth->fFShadSystP4.Pz());
         fFShadSystE->push_back(truth->fFShadSystP4.E());
 
         fIsCharm->push_back(truth->fIsCharm);
         fCharmHadronPDG->push_back(truth->fCharmHadronPdg);
         fIsStrange->push_back(truth->fIsStrange);
         fStrangeHadronPDG->push_back(truth->fStrangeHadronPdg);
         fNumProton->push_back(truth->fNumProton);
         fNumNeutron->push_back(truth->fNumNeutron);
         fNumPi0->push_back(truth->fNumPi0);
         fNumPiPlus->push_back(truth->fNumPiPlus);
         fNumPiMinus->push_back(truth->fNumPiMinus);
         fResNum->push_back(truth->fResNum);
         fDecayMode->push_back(truth->fDecayMode);
 }//FillGTruth()

void MillGenTree::MillTrees ( )

overridevirtual

Implements garana::Mill.

Definition at line 43 of file MillGenTree.cxx.

                             {
 
     if(!IsVerified()){
         std::cerr << "WARNING(MillGenTree::MillTrees): "
                           << "trying to mill trees that have not been verified"
                           << std::endl;
         return;
     }
 
     std::cout << "MillGenTree: loop over " << fGenIn->NEntries() << " gen entries" << std::endl;
 
     //loop over input genTree events (entries)
     for(size_t ientry=0; ientry<fGenIn->NEntries(); ientry++){
 
         //fGenOut->ClearVecs(); //clear previous entry's data (if any)
         ClearVecs();
         fGenIn->GetEntry(ientry);
 
         if(branchToDrawOpt[kGIndex])
                 fGIndex = fGenIn->GetGIndex();
 
         if(branchToDrawOpt[kFSParticle]) {
             fNGen = fGenIn->NGen();
                         std::cout << "MillGen: output NGen = " << fNGen << std::endl;
         }
 
                 //loop over gen entries (usually GTruths with 1 per event)
         std::cout << "MillGenTree: loop over " << fGenIn->NGen() << " igen" <<std::endl;
                 for(size_t igen=0; igen<fGenIn->NGen(); igen++){
 
                         if(branchToDrawOpt[kFSParticle]) {
 
                                 //loop over FSparticles associated with this GTruth
                                 auto const& particles = fGenIn->GetParticles(igen);
                                 fNFS->push_back(particles->size());
                                 std::cout << "MillGenTree: igen->" << igen << " NFSP->"
                                                   << fNFS->back() << std::endl;
                                 FillFSParticle(particles);
 
                         }//if filling FSParticle branches
 
                         if(branchToDrawOpt[kGTruth]) {
 
                                 const GTruth* truth = fGenIn->GetGTruth(igen);
 
                                 if(!truth)
                                         FillEmptyGTruth();
                                 else
                                         FillGTruth(truth);
 
                         }//if filling GTruth branches
                 }//for gen entries
 
         fTreeIn->Fill(); //actually our fTreeOut, but it's called fTreeIn in FlatGenTree.h, fill once per event
     }//for genTree entries
 
     fTreeIn->Write();
 
 }//Mill()

bool MillGenTree::VerifyBranches ( )

overrideprivatevirtual

Implements garana::Mill.

Definition at line 103 of file MillGenTree.cxx.

                                  {
 
     // get list of branches and check it matches what we expect
     const TObjArray* branches = fGenIn->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(MillGenTree::VerifyBranches): no branches found in passed input tree"
                           << std::endl;
         return false;
     }
 
     if(branches->GetEntriesFast() != (Int_t)nameToGenBranch.size())
         std::cout << "WARNING(MillGenTree::VerifyBranches): Mismatch in number of branches (expected "
                   << nameToGenBranch.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(nameToGenBranch.find(CharStarToString(branch->GetFullName()))!=nameToGenBranch.end()) {
 
                 std::cout << "   chopping branch '" << branch->GetFullName() << "'" << std::endl;
 
                 if(nameToGenBranch[CharStarToString(branch->GetFullName())] == kEvent)
                         branchToDrawOpt[kEvent] = true;
 
                 if(nameToGenBranch[CharStarToString(branch->GetFullName())] == kGIndex)
                         branchToDrawOpt[kGIndex] = true;
 
                 if(nameToGenBranch[CharStarToString(branch->GetFullName())] == kGTruth)
                         branchToDrawOpt[kGTruth] = true;
 
                 if(nameToGenBranch[CharStarToString(branch->GetFullName())] == kFSParticle)
                         branchToDrawOpt[kFSParticle] = true;
 
         }//endif known branch
 
         else{
                 std::cout << "WARNING(MillGenTree): ignoring unknown branch '"
                                   << branch->GetFullName() << "'" << std::endl;
         }//else
 
     }//for branches
 
     fIsVerified = true;
     return true;
 
 }//VerifyBranches()

Member Data Documentation

std::map<GenBranch,bool> garana::MillGenTree::branchToDrawOpt

private

Initial value:

= {
                    {kEvent,      false},
                        {kGIndex,     false},
                        {kGTruth,     false},
                        {kFSParticle, false}
                }

Definition at line 51 of file MillGenTree.h.

StructuredGenTree* garana::MillGenTree::fGenIn = nullptr

private

Definition at line 32 of file MillGenTree.h.

std::map<std::string,GenBranch> garana::MillGenTree::nameToGenBranch

private

Initial value:

=       {
                        {"Event",      kEvent},
                        {"GIndex",     kGIndex},
                        {"GTruth",     kGTruth},
                        {"FSParticles", kFSParticle}
                }

Definition at line 44 of file MillGenTree.h.

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

garana/garana/Processors/MillGenTree.h
garana/garana/Processors/MillGenTree.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