00001 // $Id: Example_readfz_plot.C,v 1.4 2000/04/18 20:37:25 kathy Exp $ 00002 // $Log: Example_readfz_plot.C,v $ 00003 // Revision 1.4 2000/04/18 20:37:25 kathy 00004 // St_DataSet,St_DataSetIter,St_Table classes are nowchanged to TDataSet,TDataSetIter,TTable 00005 // 00006 // Revision 1.3 1999/05/21 15:45:02 kathy 00007 // fixed owner on Example_readfz_plot.C 00008 // 00009 // Revision 1.2 1999/05/21 15:33:49 kathy 00010 // made sure Log & Id are in each file and also put in standard comment line with name of owner 00011 // 00012 //======================================================================= 00013 // owner: Pavel Nevski 00014 // what it does: 00015 //======================================================================= 00016 00017 // readfz_plot.C 00018 // 00019 // reads in geant fz file, gives a browser and then you can make a plot 00020 00021 { 00022 gROOT->LoadMacro("geant.C"); 00023 Load(); 00024 geant(1,"/disk1/star/test/psc0049_08_40evts.fzd"); 00025 //this will run 1 event and give you a browser 00026 // if you want to plot something from geant/event/g2t_track then 00027 00028 TDataSet *yp = chain->DataSet("geant"); 00029 //pointer to geant set 00030 00031 TDataSetIter geantIter(yp); 00032 // iterator 00033 00034 St_g2t_track *pt = (St_g2t_track *) geantIter.Find("g2t_track"); 00035 // pointer to g2t_track 00036 00037 g2t_track_st *t = pt->GetTable(); 00038 // pointer to table 00039 00040 TH1F *h1 = new TH1F("g2t_rap","rapidity",100,-7.0,7.0); 00041 //define histogram 00042 00043 for (Int_t i = 0; i < pt->GetNRows(); i++,t++){ h1->Fill(t->rapidity);} 00044 //fill histogram 00045 00046 TCanvas *c1 = new TCanvas("c1","c1"); 00047 // make canvas 00048 00049 h1->Draw(); 00050 // on screen 00051 } 00052
1.5.9