test9.C

{
// #include "St_mevsim_Module.h"
// void mevsimtest()

//   Read XDF file
   St_XDFFile  xdf("/star/mds/data/SD98/auau200/evg/central/hijing/set0001/regular/auau_ce_b0-2_4801_5000.xdf");
//   St_XDFFile xdf = St_XDFFile("/afs/rhic/star/packages/dev/params/tpc/tpg_pars.xdf");
   St_DataSet *event = xdf.NextEventGet(); 
   event->ls("*");
   St_DataSet *set=0;
   St_DataSetIter root(event);
   set = root.Cd("/evgen/particle");
   if (set) { 
     set->ls("*");
     printf(" Getting the tables \n");

     St_particle *pa=set->GetTableObj();


     printf(" Checking the results \n");

     table_head_st *t1_h = pa->GetHeader();

     printf(" Name=%s type=%s t1_h.nok = %i \n", t1_h->name, t1_h->type, t1_h->nok);
     printf(" t1_h.rbytes = %i \n", t1_h.rbytes);

     particle_st *particle = pa.GetTable();

     // Create "histograms"

//*-*  nt/plot 666.phep(3)

     TH1F h1("h1","phep(3)",100);
     h1->SetFillColor(16);

//*-*  nt/plot 666.sqrt(phep(1)**2+phep(2)**2) (0.ne.phep(1))

     TH1F h1Sqrt("h1Sqrt","666.sqrt",100);


//*-* nt/plot 666.log(tan(0.5*atan(phep(3)/sqrt(phep(1)**2+phep(2)**2))+.7854)) (0.ne.phep(1))

     TH1F h2Log1("h2Log1","666.log1",100);
     h2Log1->SetFillColor(42);

//*-* nt/plot 666.log(tan(0.5*atan(phep(3)/sqrt(phep(1)**2+phep(2)**2))+.7854)) (idhep.eq.211.or.idhep.eq.-211)\n");

     TH1F h3Log2("h3Log2","666.log2",100);
     h3Log2->SetFillColor(46);

     TSlider *slider = 0;

     Int_t kUPDATE = t1_h.nok/20;

//*-* Create a canvas to show the result (based on  root/tutorials/hsum.C )

     c1 = new TCanvas("c1","The reading STAF table: \"particle.h\"",200,10,600,400);
     c1->SetGrid();
     gBenchmark->Start("hsum");

     Int_t l = 0;
     for (l=0; l < pa->GetNRows(); l++) 
     {
       particle_st *p = particle[l];

    //*-* Fill a'la nt/plot 666.phep(3)

        h1->Fill(p->phep[2]);
    //
    //
    // Fill a'la nt/plot 666.sqrt(phep(1)**2+phep(2)**2) (0.ne.phep(1))
    //                         and 
    //           nt/plot 666.log(tan(0.5*atan(phep(3)/sqrt(phep(1)**2+phep(2)**2))+.7854)) (0.ne.phep(1))
    //
    // in a single step because of the common condition:  (0.ne.phep(1))
    //

    // ******  be aware CINT feature -> Non standard expressions "**" follow

        if ( p->phep[0] != 0 ) 
        {
          h1Sqrt->Fill(sqrt(p->phep[0]**2+p->phep[1]**2)); 
          h2Log1->Fill(sqrt(tan(0.5*atan(p->phep[2]/sqrt(p->phep[0]**2+p->phep[1]**2))+.7854)));
        }

    //*-*
    //*-* Fill a'la nt/plot 666.log(tan(0.5*atan(phep(3)/sqrt(phep(1)**2+phep(2)**2))+.7854)) (idhep->eq.211.or.idhep->eq.-211)\n");
    //
        if ( p->idhep ==211 || p->idhep == -211) 
          h3Log2->Fill(log(tan(0.5*atan(p->phep[2]/sqrt(p->phep[0]**2+p->phep[1]**2))+.7854)));

    // Update the view of these histograms (just for fun)

        if (l && (l%kUPDATE) == 0) {
          if (l == kUPDATE) {
             h1Sqrt->Draw("e1p");
             h2Log1->Draw("same");
             h1->Draw("same");
             h3Log2->Draw("same");
             c1->Update();
             slider = new TSlider("slider","test",1.05,0,1.1,h1Sqrt->GetMaximum()*1.3,38);
             slider->SetFillColor(46);
          }
          if (slider) slider->SetRange(0,Float_t(l)/t1_h->nok);
          c1->Modified();
          c1->Update();
        }
     }

     slider->SetRange(0,1);
     c1->Modified();
     gBenchmark->Show("hsum");
  }
  printf(" This is a finish \n");
}