{

	Double_t particle_energy;
	Double_t particle_xvertex;
	Double_t particle_yvertex;
	Double_t particle_zvertex;
	Int_t event;
	
	TreePGunInfo  = new TTree("TreePGunInfo","ParticleGunInfo");
	TreePGunInfo->Branch("particleenergy",&particle_energy,"particle_energy/D");
	TreePGunInfo->Branch("particlexvertex",&particle_xvertex,"particle_xvertex/D");
	TreePGunInfo->Branch("particleyvertex",&particle_yvertex,"particle_yvertex/D");
	TreePGunInfo->Branch("particlezvertex",&particle_zvertex,"particle_zvertex/D");
	TreePGunInfo->Branch("event",&event,"event/I");
	
	TFile* f= new TFile("f.root","RECREATE");
	f->cd();
	
	for (Int_t i =0; i< 300001; i++) {
		
		event=i;
		particle_energy=gRandom->Uniform(i,i+3);
		particle_xvertex=gRandom->Uniform();
		particle_yvertex=gRandom->Uniform();
		particle_zvertex=gRandom->Uniform(i,i+1);

		TreePGunInfo->Fill();

		if (i%100000==0) { 
			TreePGunInfo->Print();
			TreePGunInfo->Write();
			TreePGunInfo->Reset();
		}
	}

f->Close();
}