TH1F *eppx = new TH1F ("eppx", "eppx (All)", 200, -2000, 5000.);
(.../...)
RhoTuple *ntp  = new RhoTuple("RecoTuple","Reco_analysis");
(.../...)
PndAnalysis* theAnalysis = new PndAnalysis();
if (nevts==0) nevts= theAnalysis->GetEntries();
// *** RhoCandLists for the analysis
RhoCandList eplus;
(.../...)
while (theAnalysis->GetEvent() && i++<nevts)
	{
		if ((i%100)==0) cout<<"evt " << i << endl;
		// *** Select with no PID info ('All'); type and mass are set 		
		theAnalysis->FillList(chrg,    "Charged");
		theAnalysis->FillList(eplus,  "ElectronAllPlus");
(.../...)

	for (j=0;j<eplus.GetLength();++j) 
		{
		 eppx->Fill(eplus[j]->Px());
		 ntp->Column("eppx",           (Float_t) eplus[j]->Px(),                 -999.0f);
		 ntp->DumpData();
                }
(.../...)
}

out->cd();
eppx->Write();
out->Save();
ntp->GetInternalTree()->Write();
out->Close();