mevsimtest.C

{
   //  Delete the "old" log file
   gSystem.Exec("rm *.log");

printf(" Creating the tables \n");

   St_particle  *pa = new St_particle(30000);

printf(" Calling the 'mevsim' module \n");

   mevsim(pa);

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);

   hfile = new TFile("mevsim.root","RECREATE","mevsim simulation");
   particle_st *particle = pa.GetTable();

// Create "ntuple"

   TNtuple *ntuple = new TNtuple("ntuple","Demo ntuple","idhep:p1:p2:p3:p4:p5");

// Fill "Ntuple"
printf(" Filling ntuple \n");

   for ( Int_t l=0; l<t1_h.nok; l++) 
   {
      ntuple->Fill(particle[l]->idhep
                  ,particle[l]->phep[0]
                  ,particle[l]->phep[1]
                  ,particle[l]->phep[2]
                  ,particle[l]->phep[3]
                  ,particle[l]->phep[4]
                  );
   }

// Show some results

printf(" some KUMACS like results to show \n");

//*-*
//*-* nt/plot 666.phep(3) \n");
//
     ntuple->Draw("p3");

//*-*
//*-*  nt/plot 666.sqrt(phep(1)**2+phep(2)**2) (0.ne.phep(1))
//
     ntuple->Draw("sqrt(p1**2+p2**2)", "p1 != 0");

//*-*
//*-* nt/plot 666.log(tan(0.5*atan(phep(3)/sqrt(phep(1)**2+phep(2)**2))+.7854)) (0.ne.phep(1))
//
     ntuple->Draw("log(tan(0.5*atan(p3/sqrt(p1**2+p2**2))+.7854))", "p1 != 0");

//*-*
//*-* 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");
//
     ntuple->Draw("log(tan(0.5*atan(p3/sqrt(p1**2+p2**2))+.7854))", "idhep == 211 || idhep==-211");
     
hfile->Write();

}