i.e somefunction() returns the difference of its two arguments and I have a TTree called MyTree which contains two branches called ‘x’ and ‘y’. By executing the TTree::Draw() command as I did above, a histogram will be produced of the difference (or whatever I write somefunction() to do) of x and y for each entry. What are the limitations of this feature? I tried doing something similar where I pass a pointer to an array of floats to somefunction() but I receive a segmentation fault. A clear indication of what is and is not possible would be greatly appreciated. Thanks in advance.
The functions used by TTree::Draw have to take simple numerical arguments and return simple numerical arguments. I.e. only int, floats, double, etc. and no array or objects.
To use function works on arrays, you should use the new TTree Proxy mechanism released in ROOT 4.00/08 where you can run a function in context where the branchname can be used as variables.
But is it possible to use to MakeProxy mechanism in the selection variable of TTree::Draw() (as it was in PAW)? For example,
MyTree.Draw(“somefunc.C++”,“func1.C++<a && func2.C++>b”), where a and b are numbers.
If not is there any ROOT analog of such a mechanism in PAW which allows to use automatic histogram creating (like Draw(), in contrast with MakeClass) with user defined function of Ntuple variables (such as arrays) in varexp and selection?
You can use:MyTree.Draw("somefunc.C+","somesel.C+");(+ vs ++ to avoid unnecessary recompilation).
Currently you can not pass argument to the generated code from inside the Draw command itself. However an alternative is to use something like:[code]// file somefunc.C
double a;
double b;
double func1() {
…
}
double func2() {
}
void somefunc() {
if (func1()<a && func2()>b) {
double value = …;
htemp->Fill(value);
}
}[/code]
and use with:MyTree->MakeProxy("sel","somefunc.C","","nohist");
gROOT->ProcessLine(".L sel.h+");
sel *mysel = new sel();
sel->a = val1;
sel->b = val2;
MyTree->Process(mysel);
TTree::Draw can not pass object to a function, it can only pass numerical values. [I.e. all three of your example should have failed with a nice warnign message ]
root [1] tree->MakeProxy("proxytest01", "somefunc.C", "", "nohist");
Fatal in <ROOT::TBranchProxyClassDescriptor::>: strcmp(fInfo->GetName(), type)==0 violated at line 117 of `treeplayer/src/TBranchProxyClassDescriptor.cxx' Humm … this seems to indicate an issue with the ROOT type system, it is missing some information about your classes! Could you send me your ROOT file so that I can try to reproduce this problem?
]