GSI Forum: General » Reading a tree HELP!

Home » PANDA » PandaRoot » General » Reading a tree HELP!

Show: Today's Messages :: Polls :: Message Navigator

Re: Reading a tree HELP! [message #17784 is a reply to message #17783]

Thu, 22 January 2015 15:48

asanchez
Messages: 350
Registered: March 2006

first-grade participant

From: *kph.uni-mainz.de

Hi Mamaen,
here you are some lines as example how to get read of
the TParticle Info in a TTree Smile

.

TFile *tf1 = new TFile("the file");
TTree *data = (TTree*)tf1->Get("data");
TClonesArray* hit_array=new TClonesArray("TParticle",100);
data->SetBranchAddress("Particles",&hit_array);//Branch names

for (Int_t j=0; j<data->GetEntries(); j++)
{
// if ((j)%100000==0)cout <<"evt: "<<j<<endl;

data->GetEntry(j); // Loop over particles in TClonesArray

for (Int_t iPart=0; iPart < nParts; iPart++) {

TParticle* part = (TParticle*) hit_array->At(iPart);
Int_t pdgType = part->GetPdgCode();

//you can use the member function of TPArticle to get the kinematic information

TLorentzVector lv1;

part->Momentum(lv1);

TVector3 v1;

v1= lv1.Vect();

// Total momentum

Double_t P;

P = part->P();

// theta

Double_t theta;

theta = part->Theta();

}

}

Report message to a moderator

[Message index]

		Reading a tree HELP! By: Mamen on Thu, 22 January 2015 15:31
		Re: Reading a tree HELP! By: asanchez on Thu, 22 January 2015 15:48
		Re: Reading a tree HELP! By: Mamen on Thu, 22 January 2015 18:08
		Re: Reading a tree HELP! By: asanchez on Fri, 23 January 2015 00:13
		Re: Reading a tree HELP! By: Ralf Kliemt on Thu, 22 January 2015 15:54
		Re: Reading a tree HELP! By: Mamen on Thu, 22 January 2015 16:10
		Re: Reading a tree HELP! By: Ralf Kliemt on Fri, 23 January 2015 09:40

Previous Topic:	Howto for eventdisplay
Next Topic:	simulation output; sge scripts in macro/prod

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Apr 26 22:05:15 CEST 2024

Total time taken to generate the page: 0.01212 seconds