gmane.comp.lang.c++.root

Aux files for TSelector in PROOF

Christian Holm Christensen — 2013-05-22T09:26:53

Hi all, 

First off, sorry if this question has been asked before. 

I have a TSelector (wrapped by ALICE's ANALYSIS framework) which needs
an auxiliary file - that is, a file that's not part of the data-set
being analysed - the file contains correction objects (stored in a
TTree) needed by the TSelector to properly process the data-set. 

So my question is, how do I push this file to the workers along with the
job?  Preferably, the file should go into the working directory of the
individual jobs.  

Of course, one can put the file on some shared location e.g., NFS,
rootd://, or http(s):// server, but since the code should work uniformly
across environments and clusters, this solution seems more troublesome
than it's worth. 

The code must work in all three environments supported by ALICE: AliEn
(Grid), Proof, and local.  On AliEn, one can simply upload the auxiliary
file with the job.  For local jobs, the file can be a symlink or the
like.  But for Proof it is not so obvious to me. 

Note also, that I cannot really change the TSelector code - all I can
change is a limited set of functionality called by the TSelector. 

The protected member function TProof::SendFile or perhaps the private
TProof::BroadcastFile seems to do what I want, but they are obviously
not accessible.  In principle I could roll my own `SendFile': 


  Int_t SendFile(const TString& fileName) 
  {
    Int_t    bufSize = 32768;
    Char_t   buf[bufSize];
    Long64_t size = 0;
    Long_t   id = 0, flags = 0, modtime = 0;
    if (gSystem->GetPathInfo(fileName.Data(), &id, &size, &flags, &modtime)==1 
|| size <= 0) {
      Error("SendFile", "Cannot stat %s", fileName.Data());
      return -1;
    }
    TString fn(gSystem->BaseName(fileName.Data()));
    TList*  slaves = gProof->GetListOfActiveSlaves();
    TIter   next(slaves);
    TSlave* sl   = 0;
    Int_t   ret  = 0;
    Int_t   fd = open(fileName.Data(), O_RDONLY);
    while ((sl = static_cast<TSlave*>(next()))) {
      if (!sl->IsValid()) continue;
      if (sl->GetSlaveType() != TSlave::kSlave) continue;
      snprintf(buf,busSize,"%s %d %lld %d", fn.Data(), 1, size, 0);
      if (sl->GetSocket()->Send(buf, kPROOF_SENDFILE) == -1) {
Warning("SendFile", "Could not send kPROOF_SENDFILE request");
continue;
      }
      lseek(fd, 0, SEEK_SET);  // Begining of file 
      Int_t len = 0;
      do { 
while ((len = read(fd, buf, bufSize)) < 0 && TSystem::Errno() == EINTR)
  TSystem::ResetErrno();
if (len < 0) { 
  Error("SendFile", "error reading input");
  close(fd);
  return -1;
}
if (len > 0 && sl->GetSocket()->SendRaw(buf, len) == -1) {
  Error("SendFile", "error writing to slave");
  sl = 0;
  break;
}
      } while (len > 0);
      ret ++;

      if (sl) Collect(sl, gEnv->GetValue("Proof.CollectTimeout", -1));
    }
    close(fd);
  }
          
Except, of course, that TProof::GetListOfActiveSlaves and
TProof::Collect are protected/private and I cannot find an equivalent
public member functions. 

Any input would be welcome.  

Thank you. 
          
Yours,

Get histograms with `/' in the name

Federico Ferri — 2013-05-21T13:54:51

Dear experts,

I've been given a root file containing histograms with `/' in their
names. How can I get a hold on them from the TDirectory or gDirectory?

The TBrowser works, so I see them and I can plot them, but as soon as I
try to get a pointer to the histograms I get a null one.
Few lines to reproduce the behaviour are below.

From the code of TDirectory::Get(const char *namecycle) it is clear why
I experience this behaviour. Is there another method that can be used?

Thank you in advance,
Federico


--------------------
root [0] TFile * fout = new TFile("test.root", "recreate")
root [1] TH1F * h = new TH1F("ciao/pippo", "ciao/pippo", 100, 0., 1.)
root [2] h->Fill(0.5)
root [3] fout->Write()
root [4] fout->Close()

root -l test.root
root [0] 
Attaching file test.root as _file0...
root [1] _file0->Get("ciao/pippo")
(class TObject*)0x0
root [2]

SetTitleOffset for TEfficiency

Marek Sirendi — 2013-05-16T15:46:07

Bug in TMatrixTBase

Omar Andrés Zapata Mesa — 2013-05-13T16:35:39

Hello dear ROOTers,
I found a bug in TMatrixTBase using Option "F" for fortran representation
in the method GetMatrix2Array(Element *data,Option_t *option).

When I get tha matrix's data in an array with fortran representation and I
try to create a matrix again using TMatrixT(nrows,ncols,*data,"F"), the
matrix is corrupted.

I attach the macro to test the bug and the patch.

Thanks.

X axis title in a stack histogram

Pilar Casado — 2013-05-09T07:58:23

Hello,

How can I set the X axis title of a stack histogram.
I have tried the usual SetTitle on the stack histogram,
have also set two of the histogram in the stack, but
I don't manage.

Many thanks, Pilar

TTimer time resolution

Irakli Keshelashvili — 2013-05-07T11:36:44

Dear ROOTers,

I'm using TTimer for VMEbus readout. As I learned, the minimum time resolution is 10ms as it is defined in TSystem.h -

kItimerResolution = 10      // interval-timer resolution in ms

Question: is it O.K. if one changes it to 1ms? Or is there some other method to decrease time interval?

Best regards
Irakli

* * * * * * * * * * * * * * *
 Dr. Irakli Keshelashvili
 Department of Physics
 University of Basel
 Klingelbergstrasse 82,
 CH-4056 Basel, Switzerland
 Tel. +41 6126 73775 (Office)
      +41 7887 06657 (Mobile)
 Fax: +41 6126 73784
 * * * * * * * * * * * * * * *

TH1F* hists froma THStack plot

Maria Pilar Casado — 2013-05-03T13:29:14

Dear root expert,

I'm using root version 5.34/05  14 February 2013

I would need to retrieve the list of histograms from a THStack histogram
and add them into another histogram. I'm trying to do it in this way by
I get the complain you can see at the end. How can I do it?

Many thanks, Pilar

   TH1F * h_signalf = new TH1F("h_signalf","h_signalf",50,0.,500.);
   TList *thisList = h_stack->GetHists();

   TObjLink *lnk = thisList->FirstLink();
   while (lnk) {
    //h_signalf->Add( (TH1F*) lnk->GetObject());
    h_signalf->Add( (TH1F*) lnk);
    lnk = lnk->Next();
   }
--------------------------------
Error:
Error: Can't call TObjLink::Next() in current scope KolmoTest.C:45:
Possible candidates are...
(in TObjLink)

extracting systematics from profile

Nabil Ghodbane — 2013-05-02T16:02:21

dear experts,
I am trying to submit again this question to which i would be grateful 
if someone could confirm the bottom approach.
thanks in advance.
I discussed this with several people, but I would need a confirmation 
that what i am doing is consistent enough.
Basically, I am trying to extract some single source systematic 
uncertainty associated to one measurement. In order to extract the 
different systematics, I am proceeding as follow:

1) start with the histogram templates for signal (nominal and 
systematics) and data driven background and produce the Roofit workspace 
using the Histfactory machinery.

2) fix all nuisance parameters to zero, fit and extract the 68% 
confidence interval.
-> This is the nominal signal strength factor, that i label: mu_nom_c, 
mu_nom_u, mu_nom_d (central, upper, lower)

3) leave all nuisance parameters free within [-2,2], perform the fit and 
extract again the associated 68% confidence interval.
-> I label these: mu_sys_c, mu_sys_u, mu_sys_d (central, upper, lower)
-> I have the set of associated values for the nuisance parameters.

4) loop over all nuisance parameters, for each of them, fix its value to 
the value found in step 3) and perform the fit w.r.t. all other ones 
left free and extract the associated 68% CL interval. I call this:
-> I label these: mu_alphaI_c, mu_alphaI_u, mu_alphaI_d (central, upper, 
lower)

5) fix all nuisance parameters to their value found in step 3) and 
perform the fit -> this allows me to extract the stat only uncertainty: 
mu_stat_c, mu_stat_u, mu_stat_d

then i compute the systematics (total and single sources ) as follow:
a) total systematics:
--------------------------
tot_sys_u = ( |(mu_sys_c - mu_sys_u)^2 - (mu_stat_c - mu_stat_u)^2 |)^(1/2)
tot_sys_d = ( |(mu_sys_c - mu_sys_d)^2 - (mu_stat_c - mu_stat_d)^2 |)^(1/2)

b) Luminosity:
------------------
lumi_sys_u = ( |(mu_lumi_c - mu_lumi_u)^2 - (mu_stat_c - mu_stat_u)^2 
|)^(1/2)
lumi_sys_d = ( |(mu_lumi_c - mu_lumi_d)^2 - (mu_stat_c - mu_stat_d)^2 
|)^(1/2)

c) Nuisance parameter associated systematics:
---------------------------------------------------------------
alpha_sys_u = ( | (mu_sys_c - mu_sys_u)^2 - (mu_alpha_c - mu_alpha_u)^2 
|)^(1/2);
alpha_sys_d = ( | (mu_sys_c - mu_sys_d)^2 - (mu_alpha_c - mu_alpha_d)^2 
|)^(1/2);

To do this on the technical aspect, I am using Minuit2 like:
mc->GetPdf()->fitTo(*data,Save(kTRUE),InitialHesse(true), 
Hesse(true),Minos(kTRUE), RooFit::Minimizer("Minuit2","minos"), 
Strategy(1), PrintLevel(3), Save(true) );

and to extract the confidence interval, i use:
ProfileLikelihoodCalculator* PL = new 
ProfileLikelihoodCalculator(*data,*mc);
PL->SetConfidenceLevel(0.68);
LikelihoodInterval* interval = PL->GetInterval();
double intervalValue = POI->getVal();
double intervalError = POI->getError();
double intervalLow = interval->LowerLimit(*POI);
double intervalUp = interval->UpperLimit(*POI);

thanks in advance.

Unsubscribing ...

Doug Schouten — 2013-04-26T16:04:34

Hello,

I found an article describing how to unsubscribe to roottalk:

http://root.cern.ch/phpBB3/viewtopic.php?f=7&t=11954

But the majordomo< at >root.cern.ch system replies with "unknown list 'roottalk'"

Can some please tell me how to unsubscribe from roottalk? I have a 
number of email aliases that I can no longer send from (but which are 
being forwarded to my current address). I need to remove all of these 
from the roottalk list.

thanks,

ROOT has moved to Git...

Fons Rademakers — 2013-04-26T15:23:59

Today ROOT has moved to Git for version control.

We have made the change because Git offers many benefits over Subversion. 
Basically all large Open Source projects have preceded us in this move. For 
more details about the why, how and where see this presentation: 
https://indico.cern.ch/conferenceDisplay.py?confId=246803

A detailed Git HowTo (http://root.cern.ch/drupal/content/git-howto) has 
been added to the ROOT website.

To checkout the ROOT sources, clone in Git parlance, do:

    git clone http://root.cern.ch/git/root.git

For the rest, in respect to building ROOT, nothing has changed.



Cheers, Fons.

how to plot intervals

Dr. Harinder Singh Bawa — 2013-04-26T08:34:55

Hello,

I am trying to plot the intervals/range in root. I mean I have say 10
points and I got 10 intervals of that point like following

100   0.5-1.5
200    0.6-1.8
300   0.7-2.3
........ like this


Actually I am trying to see how much intervals goes on increasing when I
increase the points. Any suggestion is appreciated as what would be the
best way to plot which gives the idea of the the increase in intervals as
we go high and high in points.


Thanks

screen problem finding libCore.so

Israel Ariel Gonzalez Medina — 2013-04-24T19:35:32

hi
i am facing some troubles while running a script using the screen bash
command. It may nothing to do with these forum but i ll be grateful if
someone give me a hint.
i run the script in a direct manner, just that when i try to run the
program under the screen command it doesn't find the libCore.so  library
I need it because im using ssh conextion and i need to keep running the
program even if i close the ssh connection.
thanks

Weird behavior plotting TGraphAsymmErrors

Barbara Storaci — 2013-04-24T09:41:43

Dear Root Expert,

I'm facing a weird behavior when I'm trying to plot with the option "3" a TGraphAsymmErrors
object. In fact not all the points range are filled, but for example in the simple case I attached
no filled region is present before 0, while if I draw the same graph with the option "p"
all the points are displaced. Am I missing something in the way the drawing is working?

You can find in attachment a short script reproducing the "problem", and two pdf files
with the graph plotted with option "3" and "p".

Can you suggest me the correct way to obtain my band-plot?

Thanks a lot in advance,

      Barbara

=============================================================
Barbara Storaci
Universität Zürich

CERN, office 13-1-034
CH-1211 Geneve
Postbox: F26110

Tel: +41 (0) 227679942
Email: Barbara.Storaci< at >cern.ch<mailto:Barbara.Storaci< at >cern.ch>
=============================================================

TFile

Pilar Casado — 2013-04-21T07:13:14

Dear root experts,

I have the following piece of code:
---------------------------------------------------------
#include <string>
#include <iomanip>
#include "AtlasStyle.C"
#include "AtlasUtils.C"

void StackPlots(){
//#ifdef __CINT__
//   gROOT->LoadMacro("AtlasStyle.C");
//   gROOT->LoadMacro("AtlasUtils.C");
//#endif

   SetAtlasStyle();
   gStyle->SetPadRightMargin(0.15);

   gStyle->SetPalette(1);

   TFile *file_data=new TFile("output_data.root", "READ");
}
------------------------------------------------------------
I'm using root version 5.34/05 14 Feb. 2013

When I try to compile I get the error:
------------------------------------------------------------
/nfs/pic.es/user/c/casado/HpTriggerMar2013/./StackPlots.C:17: error: 
invalid use of incomplete type âstruct TFileâ
/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase/x86_64/root/5.34.05-x86_64-slc5-gcc4.3/include/TUUID.h:39: 
error: forward declaration of âstruct TFileâ
---------------------------------------------------------------
How should I do the declaration of file_data?

Many thanks, Pilar

Installing root in cygwin

Dr. Harinder Singh Bawa — 2013-04-19T15:10:21

Hello

I am installing root in cygwin (installed on Win7). I downloaded .tar file
and untarred in /usr/local   dir.

Setup ROOTSYS, PATH and LD_LIBRARY_PATH.


Everything seems to be well, but when I give command as root

It displays Welcome logo but Root prompt doesnt come. Initially it looks
like its taking time but I waiting for 5-10 min and looks like it doesnt
get through. Can anybody suggest what could be the problem? I can kill
session using CTRL C but again when I type root command, it gets hanged
after root logo.


Harinder

Fitting with power function

Anatoliy Khmelevskiy — 2013-04-18T14:06:08

Hi,

I have difficulty with the fitting the data by a power function.
I have

TF1 *myfit = new TF1("myfit","[0]*pow(x,[1])*exp(x*[2])", 0, 2.5);
myfit->SetParName(0,"alpha");
myfit->SetParName(1,"beta");
myfit->SetParName(2,"gamma");

...

gr = new TGraph(18, x, y);
gr->Fit("myfit", "w");

I get the result

****************************************
Minimizer is Minuit / Migrad
Chi2                      =  4.67788e+23
NDf                       =           15
Edm                       =         -nan
NCalls                    =          143
alpha                     =        30303   +/-   4729.9
beta                      =            0   +/-   0
gamma                     =      6.81021   +/-   0.0548613

beta always = 0

where I have a mistake?

thanks.

Using mailing list for other than ROOT activity....

Fons Rademakers — 2013-04-18T09:51:36

Dear All,

    of course Nick has not really been banned for life and his offense was 
not that severe and advocated a good cause. Nick has apologized about using 
the list, without prior approval, for non ROOT related matters and promised 
to not do it again. Case closed.

Cheers, Fons.

Reading and displaying root files in a browser

Pavlo Prykhodko — 2013-04-18T08:45:22

Dear ROOTtalkers,

I'm working on a project where I have to create a web-GUI program. In this GUI I'd like to display histograms, pedestals, graphs etc. that are encapsulated in .root files. Basically the program should be able to read .root files (using Python preferably) and display the necessary graphs in a browser.

I've seen that there are people who work on the development of JavaScript based solution for this purpose. And it's promising. However, I'd like to know if anyone has some other solutions?

Have a nice Thursday!

readding tree error

Israel Ariel Gonzalez Medina — 2013-04-17T19:19:02

hello community.
i am facing a reading tntuple problem.
the thing is that i need to read a tntuple several time from differents
branches and i can not do it in a proper manner. I need to read a
tntuple from a text file, so i create a class containing the tntuple
DataHelper_Hyp.hh. I check the tntuple by get it (line 4) an print it
and showing all the entries. the varlist is
"Cos:n1:n2:p1:p2:nn:np" (line 2)
Now i am about to read. Basically what i do is to read the 25 th entry
fist from the n1 branch and later from the n2 branch. 

1gROOT->ProcessLine(".L helpers/DataHelperH.cc+");
2DataHelper_Hyp *DhhCos = new
DataHelper_Hyp("data/hyperon/C12/C12_CosTheta_0.out","Cos:n1:n2:p1:p2:nn:np");
3DhhCos->ReadData();
4TNtuple *t = DhhCos->GetTNtuple();
5t->Print();
6const char* Branchx = "Cos";
7const char* Branchy1 = "n1";
8const char* Branchy2 = "n2";
9t->Show(25);
10Float_t x,y;
11x = new Float_t[1];
12y = new Float_t[1];
14t->SetBranchAddress(Branchx, &x);
15t->SetBranchAddress(Branchy1, &y);
16t->GetEntry(25);
17t->Show(25);
18std::cout << "the values for branch " << Branchy1  << " are: x = " <<
x << ", y  " << y << std::endl;

19x = new Float_t[1];
20y = new Float_t[1];
21t->SetBranchAddress(Branchx, &x);
22t->SetBranchAddress(Branchy2, &y);
23t->GetEntry(25);
24t->Show(25);
25std::cout << "the values for branch " << Branchy2  << " are: x = " <<
x << ", y  " << y << std::endl;
26t->Show(25);

The problem is that, if i prefer to read in first place the n2 banch
(note that the n2 branchs is define after the n1's in the varlist) by
change line 7 nd 8:
7 const char* Branchy1 = "n2";
8 const char* Branchy2 = "n1";

the reading for teh first value for n2 branch is ok, but when i read the
n1 branch valor, the reading data is wrong, even more, when i print the
all entry by the command t->Show(25); i realize that the tntuple entry
for the n1 branches has changed by the value of the previous read value
for the n2 branches.
can anybody explain me this?
thanks.

Pojects for kids.

Nick van Eijndhoven — 2013-04-17T11:22:35

Dear all,
Allow me for this one time to use this mailing list to forward an email
from my daughter Anouk who,
together with a colleague, has taken up an initiative to help traumatized
children under the umbrella
of the non-profit organisation Dream4kids.
In case this initiative appeals to you, please forward this email to
friends, colleagues etc.

Cheers,
Nick.


Dear all,



We would like to draw your attention to the following:

We, Nadine and Anouk, work as volunteers for the Netherlands foundation
Dream4kids. This foundation tries to create a positive change in the lives
of traumatized children. The purpose of this foundation is to break through
the negative spiral in the lives of these children by organizing a ‘dream
day’. By realizing one dream, we give these children courage and energy to
enjoy life again and to look forward to new opportunities. However, in
order to achieve this, we need your help!

A dream day obviously has a price tag and via the Dream4kids foundation we
have taken up the effort to gather the necessary funding, on a non-profit
basis, to provide several children the "day of their life".
By means of this email we ask you to spread this message among as many
people as possible and consider to make a donation as well, to enable us
to organize several "dream days" and make a lot of wishes come true.

After the dream day we will write reports of these special moments.
In addition to publishing these reports on the foundation's website, we
will also send our sponsors these reports, so you know to which beautiful
goal you have contributed.


So do not hesitate, and support us!

Further information (incl. reports of the "dream days" of the past few
years) can be obtained via http://www.dream4kids.nl
.

Donations can be made to the Dream4kids-account: Rabobank 1034.08.487 or
IBAN-code: NL63 RABO 0103 4084 87 or BIC-code: RABO NL 2U (see the website
mentioned above) *under mentioning our names* (Nadine Vos and Anouk van
Eijndhoven) to make
sure that your donation gets coupled to the right project.

Thanks in advance!

Regards,
Nadine Vos and Anouk van Eijndhoven
 www.dream4kids.nl



Nick van Eijndhoven
Professor of Physics
Vrije Universiteit Brussel (VUB-ELEM)
Inter-university Institute for High Energies (IIHE)
Pleinlaan 2, B-1050 Brussel, Belgium
===============================================
Office : 0G-122  tel: +32.(0)2.629.3212   gsm: +32.(0)488.074166
Secr. : 1G-006  tel: +32.(0)2.629.3203    fax : +32.(0)2.629.3816
Email : Nick.van.Eijndhoven< at >vub.ac.be  or : nick< at >icecube.wisc.edu
Web  : http://www.iihe.ac.be   http://sites.google.com/site/nickveweb
Skype: nickve.nl <nickve.nl< at >gmail.com>   Gsm-NL: +31.(0)6.40181589

*"In the beginning there was nothing, which exploded".*

TTree GetEntry problem

Egon Pavlica — 2013-04-17T07:15:37

Dear Rooters,

I have encountered a strange behavior in ubuntu 12.04.2, which might be 
the result of system upgrades during last month.

I work with MPI, each thread creates a TTree in separate file. Then at the 
end the zero rank thread copies tracks and join them in single file. It 
was working, but now stops at GetEntry(0). If i do it in my script, 
without MPI, then it still works.

If anyone has a hint, please help. Thanks for great code.

Best Regards,
Egon

ps.the script, which is the same as in the compiled code:

{
   gSystem->Load("arch/libMCdata.so");//library with LTrack
   TFile f("data0.root");//input file
   TFile b("data.root","UPDATE");//outputfile
   f.cd();
   TTree* tree=(TTree*)gDirectory->Get("tree");//get tree from f
   LTrack* tfrom;
   TTree* tnew=new TTree("tree","tree");
   LTrack* tto=new LTrack;
   tnew->Branch("track",&tto,8000,1);
   tree->SetBranchAddress("track",&tfrom);
   //copy
   for (Int_t i=0;i<tree->GetEntries();i++){
         tree->GetEntry(i);//<-here it stops when launched with mpirun
         tto=tfrom;//LTrack has this operator
         tnew->Fill();
   }
   tnew->Write("tree",TFile::kOverwrite);
   b.Close();
   f.Close();
}