batch_write_ntuple.C

Go to the documentation of this file.

//TITLE Writing TNtuple(s) and ASCII files
 
{
 
    TFile *f = new TFile("ntuple.root", "RECREATE");
    TNtuple *ntuple = new TNtuple("ntuple", "data from Pluto events", "eta_px:eta_py:eta_pz:eta_m");
    
    TFile *f2 = new TFile("ntuple2.root", "RECREATE");
    TNtuple *ntuple2 = new TNtuple("ntuple2", "data from Pluto events", "cos_theta_eta");
 
    //Define the reaction
    PReaction my_reaction("3.5", "p", "p", "p p eta [pi+ pi- pi0]");
   
    //Writing to the ntuple is done by defining the variables of the branches:
    my_reaction.Output(ntuple,"eta_px = [eta]->Px() ; eta_py = [eta]->Py() ; eta_pz = [eta]->Pz(); eta_m = [eta]->M()");
    my_reaction.Output(ntuple2,"myeta = [eta]; myeta->Boost([p + p]); cos_theta_eta = myeta->CosTheta()");
 
    //Writing to ASCII output files works with the same syntax as the "echo" command:
    my_reaction.Output("eta_sample.txt", "echo $eta_px $eta_py $eta_pz $eta_m");
 
    my_reaction.Print();
 
    cout << my_reaction.Loop(100000) << " events recorded" << endl;
 
    f->Write();
    f2->Write();
 
}