 << ist << "\t" << pdg << "\t" << eta << "\t" << pt << endl;
     //etaH->Fill(eta);
     if (pt > 0.) {
        //ptH->Fill(pt, 1./(2. * pt));
       //ep->Fill(eta,phi,pt);
       ep[iev]->Fill(eta,phi,pt);
       ep[iev]->SetMarkerStyle(20);

     }
  }
}

pythia8->PrintStatistics();

/*
TCanvas* c1 = new TCanvas("c1","Pythia8 test example",800,800);
c1->Divide(1, 2);
c1->cd(1);
etaH->Scale(5./Float_t(nev));
etaH->Draw();
etaH->SetXTitle("#eta");
etaH->SetYTitle("dN/d#eta");

c1->cd(2);
gPad->SetLogy();
ptH->Scale(5./Float_t(nev));
ptH->Draw();
ptH->SetXTitle("p_{t} [GeV/c]");
ptH->SetYTitle("dN/dp_{t}^{2} [GeV/c]^{-2}");
*/
//TCanvas* c2 = new TCanvas("c2","Pythia8 test example2",400,400);
//ep->Draw();
file->Write();
}