GSI Forum: Bugs, Fixes, Releases » Number of entries in TTree

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » Number of entries in TTree

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

Number of entries in TTree [message #12198]

Tue, 05 July 2011 13:19

Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...

first-grade participant

From: *fuw.edu.pl

Dear colleagues,

I have observed a strangeness with number of entries in TTree for simulation, digitization, reconstruction and pid, i.e. these number are different and I expect them to be equal.

I do simulation with eta_c macros from /macro/run/tdrct/eta_c.
With 2000 events simulation with trunk rev.12565

Running simple macro on produced data

void test_data()
{
	TString inSimFileName = "evt_points_stt.root";
	TString inDigiFileName = "evt_digi_stt.root";
	TString inRecoFileName = "evt_reco_stt.root";
	TString inPidFileName  = "evt_pid_stt.root";
	
	TFile *inSimFile = TFile::Open(inSimFileName,"READ");
	TTree *treeSim=(TTree *) inSimFile->Get("cbmsim");

	TFile *inDigiFile = TFile::Open(inDigiFileName,"READ");
	TTree *treeDigi=(TTree *) inDigiFile->Get("cbmsim");

	TFile *inRecoFile = TFile::Open(inRecoFileName,"READ");
	TTree *treeReco=(TTree *) inRecoFile->Get("cbmsim");

	TFile *inPidFile = TFile::Open(inPidFileName,"READ");
	TTree *treePid=(TTree *) inPidFile->Get("cbmsim");

	std::cout<<"treeSim->GetEntriesFast()="<<treeSim->GetEntriesFast()<<std::endl;
	std::cout<<"treeDigi->GetEntriesFast()="<<treeDigi->GetEntriesFast()<<std::endl;
	std::cout<<"treeReco->GetEntriesFast()="<<treeReco->GetEntriesFast()<<std::endl;
	std::cout<<"treePid->GetEntriesFast()="<<treePid->GetEntriesFast()<<std::endl;
}

I have the following output

treeSim->GetEntriesFast()=2000
treeDigi->GetEntriesFast()=2001
treeReco->GetEntriesFast()=2002
treePid->GetEntriesFast()=2001

And those number are the same for stt and tpc simulation.

Should those numbers be equal?

Dima

Report message to a moderator

Re: Number of entries in TTree [message #12199 is a reply to message #12198]

Tue, 05 July 2011 13:32

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi,
you are correct, they must be equal.
I sent a mail to Mohammad and Florian last week because of this bug. Meanwhile, if you use the latest may11, this problem is not present, we moved to an older version of base w/o such feature.

Report message to a moderator

Re: Number of entries in TTree [message #12202 is a reply to message #12198]

Tue, 05 July 2011 13:44

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi,

the difference in the number of entries is due to modifications for the time based simulation.

They should not cause any problems for standard simulations (the additional entries are just empty) but are necessary to store the content of some buffers after the time based simulation.

Cheers,

Tobias

Report message to a moderator

Re: Number of entries in TTree [message #12203 is a reply to message #12202]

Tue, 05 July 2011 13:53

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

But the AddFriend returns a warning message, saying that the numbers of entries are different.
Is there no other place where to put that information? Maybe the parameter file?

Report message to a moderator

Re: Number of entries in TTree [message #12204 is a reply to message #12203]

Tue, 05 July 2011 14:27

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Hi Stefano,

I do not see a different place where to put it. This is ordinary data which has to be processed by the subdetector tasks. Maybe we can introduce another variable which switches of this feature for ordinary simulation or we swith off the warning if the difference between two trees is just one entry.

Cheers,

Tobias

Report message to a moderator

Re: Number of entries in TTree [message #12206 is a reply to message #12204]

Tue, 05 July 2011 14:56

Bernhard Ketzer
Messages: 3
Registered: March 2006

occasional visitor

From: *dfn.mwn.de

Dear Tobias,

I think we should find a different way to store this data. The time-based simulation will (or should) become the standard for MC simulations in PANDA. If we continue to carry on this error, we will constantly run into troubles or questions associating the events.

Best regards,

Bernhard

Report message to a moderator

Re: Number of entries in TTree [message #12208 is a reply to message #12206]

Tue, 05 July 2011 15:19

Tobias Stockmanns
Messages: 489
Registered: May 2007

first-grade participant

From: *ikp.kfa-juelich.de

Dear Bernhard and Stefano,

I do not think this is an error.

In terms of time based simulation it does not make sense to think about number of events despite in the MC simulation at the beginning. In the digitization stage the digis will be randomized and you access them not any longer by taking the same event but asking for time windows. The FairRootManager cares about providing you with the correct data.
The additional entry at the end of the tree is necessary to store all those data which is still sitting in the data buffers in the digitization tasks. This data does not differ in any kind from the rest of the data and is processed later on in the same way. To store this data at a different place does not make sense in my point of view.

We sure can discuss this on one of the next EVO meetings or during the collaboration meeting in September.

Cheers,

Tobias

Report message to a moderator

Re: Number of entries in TTree [message #12296 is a reply to message #12208]

Tue, 19 July 2011 13:19

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi,
another problem of the additional events appears when merging more than one file.

Let us assume that we have 2 mc files of 1000 events; after recosntrution we will have 2 pid files of 1001 events, the last is the dummy one.

When I want to add these 2 files, I will have 2000 mc files and 2002 pid files, but in this case the mc event #1001 will correspond to pid event #1002, breaking the friend mechanism.
This is quite nasty.

Report message to a moderator

Previous Topic:	Crash in PndTpcRiemannTrackingTask
Next Topic:	D+ and D- vertex resolution from Psi analysis

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Mon Dec 02 20:47:53 CET 2024

Total time taken to generate the page: 0.00871 seconds