Dear ROOTers,
i have a problem with TChain:
I created several PAW NTupels, converted them to ROOT TFiles and made a
plot after chaining the file. Obviously, not all events are put in the
file.
I am using ROOT * Version 2.22/10 28 July 1999 *
and tested it on SUN and SGI (IRIX 6.2) here at DESY.
Below the macro is given. It first plots the quantity for each file
seperatly and in the end gives the chained version of all files, which
should be the sum of all values before.
The output is:
root [1] plottest()
TH1.Print Name= q21, Total sum= 37063.6
fSumw[0]=13706.5
fSumw[1]=15168
fSumw[2]=13924.7
fSumw[3]=7970.94
fSumw[4]=180
TH1.Print Name= q22, Total sum= 52277.5
fSumw[0]=19210.4
fSumw[1]=21285.9
fSumw[2]=19620
fSumw[3]=11371.5
fSumw[4]=270.534
TH1.Print Name= q22a, Total sum= 16461
fSumw[0]=5469.4
fSumw[1]=6605.63
fSumw[2]=6189.12
fSumw[3]=3666.27
fSumw[4]=87.6057
TH1.Print Name= q22b, Total sum= 16392.7
fSumw[0]=5515.03
fSumw[1]=6469.1
fSumw[2]=6360.5
fSumw[3]=3563.09
fSumw[4]=100.022
TH1.Print Name= q22c, Total sum= 2.5015e+06
fSumw[0]=13522.1
fSumw[1]=2.48533e+06
fSumw[2]=10452.4
fSumw[3]=5721.22
fSumw[4]=2.52638e+08
TH1.Print Name= q23, Total sum= 51065.3
fSumw[0]=16219.9
fSumw[1]=21262.8
fSumw[2]=18916
fSumw[3]=10886.6
fSumw[4]=254.686
TH1.Print Name= q2a, Total sum= 89340.3
fSumw[0]=32917
fSumw[1]=36453.9
fSumw[2]=33545
fSumw[3]=19341.5
fSumw[4]=450.534
Obviously only the first to files are checked. Juggling with the sequence
of the files changes the number of files checked, but i could not see
a "broken" file of a maximal number of events (changing the TCut string
can reduce the number of events, but does not change number of files checked).
Chaining the original paw files in PAW works fine, so the paw files are ok.
HERE is the code :
plottest()
{
Float_t bins[4] = { 200, 300, 600, 5000};
Text_t *Cmc = "Zwgt1*Zwgt2*Weight*(Ieltyp>0 && Rele>11 && Relth>0.087266463 && Relth<2.6703538)";
TChain *c1 = new TChain("h100");
c1->Add("cteq1.root");
TChain *c2 = new TChain("h100");
c2->Add("cteq1a.root");
TChain *c2a = new TChain("h100");
c2a->Add("cteq1b2.root");
TChain *c2b = new TChain("h100");
c2b->Add("cteq1b3.root");
TChain *c2c = new TChain("h100");
c2c->Add("cteq1b.root");
TChain *c3 = new TChain("h100");
c3->Add("cteq1c.root");
TChain *ca = new TChain("h100");
ca->Add("cteq1.root");
ca->Add("cteq1a.root");
ca->Add("cteq1b2.root");
ca->Add("cteq1b3.root");
ca->Add("cteq1c.root");
ca->Add("cteq1b.root");
TH1F *h1 = new TH1F("q21","",3,bins);
c1->Draw("Q2jo>>q21",Cmc,"goff");
h1->Print("all");
TH1F *h2 = new TH1F("q22","",3,bins);
c2->Draw("Q2jo>>q22",Cmc,"goff");
h2->Print("all");
TH1F *h2a = new TH1F("q22a","",3,bins);
c2a->Draw("Q2jo>>q22a",Cmc,"goff");
h2a->Print("all");
TH1F *h2b = new TH1F("q22b","",3,bins);
c2b->Draw("Q2jo>>q22b",Cmc,"goff");
h2b->Print("all");
TH1F *h2c = new TH1F("q22c","",3,bins);
c2c->Draw("Q2jo>>q22c",Cmc,"goff");
h2c->Print("all");
TH1F *h3 = new TH1F("q23","",3,bins);
c3->Draw("Q2jo>>q23",Cmc,"goff");
h3->Print("all");
TH1F *ha = new TH1F("q2a","",3,bins);
ca->Draw("Q2jo>>q2a",Cmc,"goff");
ha->Print("all");
}
Ciao
Thomas
----------------------------------------------------------------------------
Thomas Hadig, Ottensen, Donnerstr. 20, D-22763 Hamburg (+49 40 39901186)
office : I. Phys. Institut AAC, Room 1c354, Notkestr. 85, D-22607 Hamburg
hadig@rwth-aachen.de Tel. +49 40 8998 2312
http://www-h1.desy.de/~hadig/ FAX +49 40 8998 4385
----------------------------------------------------------------------------
Windows - the user friendly interface! Brought to you by the makers of EDLIN
----------------------------------------------------------------------------
This archive was generated by hypermail 2b29 : Tue Jan 04 2000 - 00:43:39 MET