Hi Thomas,
Your result below indicates that the weight you use in your selection
contains negative values. Could you run the same test with no
selections.
Let me know.
Rene Brun
Thomas Hadig wrote:
>
> 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