TCut and TCutG

infinidhi · April 24, 2009, 6:46am

Hello Root Users,

I need help with the following problem.

First I open the file that contains a tree (named “t”), two TCut (named “goodEPS” and “goodDPS”) and one TCutG (named “goodCvsA”)
root [0] TFile *f263=new TFile(“phase5/data/a3tp5_263_bche.root”,“update”)

Then, I set TCutG pointer
root [1] TCutG goodCvsA=(TCutG)goodCvsA

Then, I create a TCut named cuts
root [2] TCut cuts=Form("(%s) && (%s) && (goodCvsA)",goodEPS->GetTitle(),goodDPS->GetTitle())

*Note:
goodEPS contains
KEY: TCut goodEPS;1 (EPS>457956 && EPS<458062) || (EPS>458172 && EPS<458185) || (EPS>458998 && EPS<459021) || (EPS>459053 && EPS<459092) || (EPS>459101 && EPS<459595) || (EPS>459656 && EPS<459777) ||(EPS>459802 && EPS<460797)

and goodDPS contains
KEY: TCut goodDPS;1 DPS>0.021 && DPS<0.045

Then, I draw the tree
root [3] t->Draw(“C:A”,cuts)

This is when I get the following error:
Error in TObjArray::At: index 16 out of bounds (size: 16, this: 0x015ccbf8)
Error in TObjArray::At: index 16 out of bounds (size: 16, this: 0x015ccbf8)
Error in TObjArray::At: index 16 out of bounds (size: 16, this: 0x015ccbf8)
(Long64_t)1006

The histogram is drawn, but just curious what this error means and whether there are errors in the displayed histogram.

When I use any two combination of the cut or graphical cut:
root [4] TCut cuts=Form("(%s) && (%s)",goodEPS->GetTitle(),goodDPS->GetTitle())
root [5] t->Draw(“C:A”,cuts)
(Long64_t)201307
root [6] TCut cuts=Form("(%s) && (goodCvsA)",goodDPS->GetTitle())
root [7] t->Draw(“C:A”,cuts)
(Long64_t)2035
root [8] TCut cuts=Form("(%s) && (goodCvsA)",goodEPS->GetTitle())
root [9] t->Draw(“C:A”,cuts)
(Long64_t)1102

This error does not occur.

Thanks,
Nidhi

pcanal · April 24, 2009, 9:25am

Hi,

Which version of ROOT do you have this problem with?

Philippe.

infinidhi · April 24, 2009, 3:22pm

Sorry, I forgot to mention that…I’m using Root 5.22/00

pcanal · April 24, 2009, 4:17pm

Hi,

Could you send me your root file so that I can try to reproduce this problem?

Thanks,
Philippe

infinidhi · April 25, 2009, 1:34am

I forgot to mention that in t->Draw(“C:A”), “A” is an alias to “RT+RB+LD+LU”.

Also, call TCut goodEPS=(TCut)goodEPS and TCut goodDPS=(TCut)goodDPS before setting TCut cuts…

I’ve tried once again to reproduce the problem but now I get segmentation error shown below (the first line is repeated many times but I’ve only include it once for short)

Error in TObjArray::AddAt: index 10 out of bounds (size: 10, this: 0x014bd5f8)

*** Break *** segmentation violation
Attaching to program: /proc/10290/exe, process 10290
[Thread debugging using libthread_db enabled]
0x00007fbd8dbcf5d5 in waitpid () from /lib64/libc.so.6
#1 0x00007fbd8db6b861 in ?? () from /lib64/libc.so.6
#2 0x00007fbd8febbdff in TUnixSystem::StackTrace() () from /usr/local/lib64/root/libCore.so
#3 0x00007fbd8febc9b5 in TUnixSystem::DispatchSignals(ESignals) () from /usr/local/lib64/root/libCore.so
#4
#5 0x00007fbd8bd8e606 in TBranch::GetEntry(long long, int) () from /usr/local/lib64/root/libTree.so
#6 0x00007fbd87497d69 in TTreeFormula::EvalInstance(int, char const**) () from /usr/local/lib64/root/libTreePlayer.so
#7 0x00007fbd874828c4 in TSelectorDraw::ProcessFill(long long) () from /usr/local/lib64/root/libTreePlayer.so
#8 0x00007fbd874a27b0 in TTreePlayer::Process(TSelector*, char const*, long long, long long) ()
from /usr/local/lib64/root/libTreePlayer.so
#9 0x00007fbd874a4059 in TTreePlayer::DrawSelect(char const*, char const*, char const*, long long, long long) ()
from /usr/local/lib64/root/libTreePlayer.so
#10 0x00007fbd8be55db4 in G__G__Tree_109_0_40(G__value*, char const*, G__param*, int) () from /usr/local/lib64/root/libTree.so
#11 0x00007fbd8f07ab1f in Cint::G__ExceptionWrapper(int ()(G__value, char const*, G__param*, int), G__value*, char*, G__param*, int)
() from /usr/local/lib64/root/libCint.so
#12 0x00007fbd8f128b30 in G__execute_call () from /usr/local/lib64/root/libCint.so
#13 0x00007fbd8f129a33 in G__call_cppfunc () from /usr/local/lib64/root/libCint.so
#14 0x00007fbd8f105bbc in G__interpret_func () from /usr/local/lib64/root/libCint.so
#15 0x00007fbd8f0f2662 in G__getfunction () from /usr/local/lib64/root/libCint.so
#16 0x00007fbd8f1d7f34 in G__getstructmem () from /usr/local/lib64/root/libCint.so
#17 0x00007fbd8f1d1c14 in G__getvariable () from /usr/local/lib64/root/libCint.so
#18 0x00007fbd8f0c538f in G__getitem () from /usr/local/lib64/root/libCint.so
#19 0x00007fbd8f0c5c8b in G__getitem () from /usr/local/lib64/root/libCint.so
#20 0x00007fbd8f0cb16f in G__getexpr () from /usr/local/lib64/root/libCint.so
#21 0x00007fbd8f158be8 in G__exec_statement () from /usr/local/lib64/root/libCint.so
#22 0x00007fbd8f0b1fd9 in G__exec_tempfile_core () from /usr/local/lib64/root/libCint.so
#23 0x00007fbd8f0b227e in G__exec_tempfile_fp () from /usr/local/lib64/root/libCint.so
#24 0x00007fbd8f15e0bb in G__process_cmd () from /usr/local/lib64/root/libCint.so
#25 0x00007fbd8fea8175 in TCint::ProcessLine(char const*, TInterpreter::EErrorCode*) () from /usr/local/lib64/root/libCore.so
#26 0x00007fbd8fde1660 in TApplication::ProcessLine(char const*, bool, int*) () from /usr/local/lib64/root/libCore.so
#27 0x00007fbd8ea3267f in TRint::HandleTermInput() () from /usr/local/lib64/root/libRint.so
#28 0x00007fbd8febab99 in TUnixSystem::CheckDescriptors() () from /usr/local/lib64/root/libCore.so
#29 0x00007fbd8febb1c1 in TUnixSystem::DispatchOneEvent(bool) () from /usr/local/lib64/root/libCore.so
#30 0x00007fbd8fe3b8e6 in TSystem::InnerLoop() () from /usr/local/lib64/root/libCore.so
#31 0x00007fbd8fe3da1c in TSystem::Run() () from /usr/local/lib64/root/libCore.so
#32 0x00007fbd8fddf5af in TApplication::Run(bool) () from /usr/local/lib64/root/libCore.so
#33 0x00007fbd8ea33f08 in TRint::Run(bool) () from /usr/local/lib64/root/libRint.so
#34 0x00000000004011cd in main ()
The program is running. Quit anyway (and detach it)? (y or n) [answered Y; input not from terminal]
Detaching from program: /proc/10290/exe, process 10290

Also, I’m running ROOT 5.22/00 on openSUSE 11.1
2.6.27.21-0.1-default #1 SMP 2009-03-31 14:50:44 +0200 x86_64 x86_64 x86_64 GNU/Linux

with gcc version
gcc (SUSE Linux) 4.3.2 [gcc-4_3-branch revision 141291]

incase you need this info.

The root file size is 13.2MB to add as an attachment. Any suggestions on how to send it to you?

Thanks for helping.
Nidhi

pcanal · April 25, 2009, 1:03pm

Hi,

You can send it to mail via email or better yet post it on a website, ftp site or shared files system (afs for example).

Cheers,
Philippe.

infinidhi · April 26, 2009, 7:23pm

Hi,

I tried the following,

TCut cuts=Form("(%s) && (goodCvsA) && (%s)",goodEPS->GetTitle(),goodDPS->GetTitle())

and

TCut cuts=Form("(goodCvsA) && (%s) && (%s)",goodEPS->GetTitle(),goodDPS->GetTitle())

I don’t get errors as I mentioned before. I don’t understand why the order of the graphical cut matters.

Do you still want me to email you the root file?

Thanks for all your quick replies.

Nidhi

pcanal · April 27, 2009, 2:30pm

Hi,

If you need any further assistance (i.e. fixing the original problem), I would still need the root file.

Cheers,
Philippe.

infinidhi · April 27, 2009, 8:47pm

Hi,

I’ve modified the original tree with fewer entries and removed non-essential branch to shrink the file size. With this new tree, I’ve been able to reproduce the problem. Only, two things have changed, the tree is now called “t1” and you don’t need to set the alias for “A”.

Let me know what you find out.

Nidhi
a3tp5_263_bche.root (1.19 MB)

pcanal · April 28, 2009, 7:23pm

Hi,

Thank you for reporting this problem. It is fixed in the trunk by the revision #28379.

Cheers,
Philippe.

infinidhi · April 28, 2009, 8:16pm

Hey thanks for fixing it. I’m glad I can help.

ROOT & ROOT support rocks! $\:D/$

Nidhi