Revision 3028e1f6d0cfa2a14a4fcc3cfe4cf5de60eb7edf authored by Wouter Verkerke on 22 March 2005, 22:11:41 UTC, committed by Wouter Verkerke on 22 March 2005, 22:11:41 UTC
o RooDataHist - Import hidden plotOn() virtual functions from RooAbsData using 'using' git-svn-id: http://root.cern.ch/svn/root/trunk@11411 27541ba8-7e3a-0410-8455-c3a389f83636
1 parent e74327f
jets.C
// script illustrating the use of a Tree using the JetEvent class.
// The JetEvent class has several collections (TClonesArray)
// and other collections (TRefArray) referencing objects
// in the TClonesArrays.
// The JetEvent class is in $ROOTSYS/tutorials/JetEvent.h,cxx
// to execute the script, do
// .x jets.C
void write(Int_t nev=100) {
//write nev Jet events
TFile f("JetEvent.root","recreate");
TTree *T = new TTree("T","Event example with Jets");
JetEvent *event = new JetEvent;
T->Branch("event","JetEvent",&event,8000,2);
for (Int_t ev=0;ev<nev;ev++) {
event->Build();
T->Fill();
}
T->Print();
T->Write();
}
void read() {
//read the JetEvent file
TFile f("JetEvent.root");
TTree *T = (TTree*)f.Get("T");
JetEvent *event = 0;
T->SetBranchAddress("event", &event);
Int_t nentries = (Int_t)T->GetEntries();
for (Int_t ev=0;ev<nentries;ev++) {
T->GetEntry(ev);
if (ev) continue; //dump first event only
cout << " Event: "<< ev
<< " Jets: " << event->GetNjet()
<< " Tracks: " << event->GetNtrack()
<< " Hits A: " << event->GetNhitA()
<< " Hits B: " << event->GetNhitB() << endl;
}
}
void pileup(Int_t nev=200) {
//make nev pilepup events, each build with LOOPMAX events selected
//randomly among the nentries
TFile f("JetEvent.root");
TTree *T = (TTree*)f.Get("T");
Int_t nentries = (Int_t)T->GetEntries();
const Int_t LOOPMAX=10;
JetEvent *events[LOOPMAX];
Int_t loop;
for (loop=0;loop<LOOPMAX;loop++) events[loop] = 0;
for (Int_t ev=0;ev<nev;ev++) {
if (ev%10 == 0) printf("building pileup: %d\n",ev);
for (loop=0;loop<LOOPMAX;loop++) {
Int_t rev = gRandom->Uniform(LOOPMAX);
T->SetBranchAddress("event", &events[loop]);
T->GetEntry(rev);
}
}
}
void jets(Int_t nev=100, Int_t npileup=200) {
gSystem->Load("libPhysics");
gROOT->ProcessLine(".L JetEvent.cxx+");
write(nev);
read();
pileup(npileup);
}
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...