Re: [BLAST_ANAWARE] commands for adding ntuples (fwd)

From: Chris & Kim Crawford (chris2@MIT.EDU)
Date: Sat Nov 02 2002 - 16:24:01 EST


hi all,
  here is an quick way to chain the 'ntp' ntuples together
using the init.C script:

bash> root -l ntp.C 1479 1480
root[0] chain[0]->Draw("whatever")
root[1] Draw16("ttl#t:ttl#b") //draw's 16 histograms automatically
root[2] CutG() // after hitting enter, click on the bounds of your
                // graphical cut; double click on last point
root[3] Draw("ttl2t:ttl3b") //it will now draw all points in black
                // and superimpose points passing the cut in red
                // (root still has problems w/ TCutG & TChain together
                // so this may not work until we upgrade root)

  for the 'eep' ntuple, use init.C instead of ntp.C.

chaining is done internally using the 'TChain' class, which
works as an ntuple across multiple files.
  TChain ntp("ntp")
  ntp->Add("new_datafile")

--hope this helps, chris

  

Tancredi Botto wrote:
>
> Sometimes you just need to string ntuples together.
> Read on
>
> -- t
> ________________________________________________________________________________
> Tancredi Botto, phone: +1-617-253-9204 mobile: +1-978-490-4124
> research scientist MIT/Bates, 21 Manning Av Middleton MA, 01949
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>
> ---------- Forwarded message ----------
> Date: Thu, 31 Oct 2002 23:24:04 -0500
> From: vitaliy ziskin <vziskin@yahoo.com>
> To: Tancredi Botto <tancredi@lns.mit.edu>
> Subject: commands for adding ntuples
>
> TFile run1("filename1")
> TFile run1("filename2")
> TFile run1("filename3")
> TFile run1("filename4")
>
> TNtuple *ntp1
> TNtuple *ntp2
> TNtuple *ntp3
> TNtuple *ntp4
>
> ntp1 = (TNtuple*)run1.Get("eep")
> ntp2 = (TNtuple*)run2.Get("eep")
> ntp3 = (TNtuple*)run3.Get("eep")
> ntp4 = (TNtuple*)run4.Get("eep")
>
> ntp1->Draw("ncr>>htemp")
> ntp2->Draw("ncr>>+htemp")
> ntp13->Draw("ncr>>+htemp")
> ntp4->Draw("ncr>>+htemp")
>
> htemp->Draw()
>
> // you don't need to declaire htemp



This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:28 EST