[ROOT] accessing infos of profile histos

From: Rosario Turrisi (rosario.turrisi@pd.infn.it)
Date: Mon Apr 10 2000 - 12:19:07 MEST


Dear Rooters,
I am building profile histograms from a 2D
via myhist->ProfileY method.
Then I need to access the information stored, 
i.e. the mean and RMS, what I am not able to do 
by the usual 1D-hist methods, as I supposed it was meant.
It seems that in any case ProfileY builds a histogram
with NbinsX as the original 2D, no matter if I do ProfileY or -X.
It seems also that, when NbinsX>NbinsY, it puts the points 
corresponding to neighboring (Y) bins every NbinsX/NbinsY bins...
For the sake of clarity, I include here the macro wishing 
to do this job:

void makermsplot() {
 TFile* p = new TFile("fitvertexs-200.root","UPDATE");
 TH2D* d0vsphi=(TH2D*)p->Get("d0vsphi");
 Int_t nbins=(Int_t)d0vsphi->GetNbinsY();
 Int_t modbins=(Int_t)d0vsphi->GetNbinsX();
 Float_t xlow = (Float_t)d0vsphi->GetYaxis()->GetXmin();
 Float_t xup = (Float_t)d0vsphi->GetYaxis()->GetXmax();
 TH1F *hres = new TH1F("hres","res vs phi",nbins,xlow,xup);
 d0vsphi->ProfileY("profy",1,1000,"s"); 
 for (Int_t i=1; i<=nbins; i++) {
   Int_t newbin =
(Int_t)((Float_t)i*(Float_t)modbins/(Float_t)nbins)-15;  // here I tried
empirically to get the filled bins...
   Float_t cont= profy->GetBinError(newbin);
   Float_t xb = (Float_t)i;
   //   hres->Fill(xb,cont);
 }
 // hres->Draw();
 //cin>>xb;
 p->Close();
 return;
}


Any suggestion?
Thanks in advance!
		Rosario

-- 
-----------------------------------------------------
Rosario Turrisi
Dipartimento di Fisica dell'Universita' e INFN 
via Marzolo, 8 - Padova
tel. +39-049-8277106
fax  +39-049-8762641

      **** ALICE & NA57 Collaboration ****

'Il faut decider si on veut etre Cesar a Rome
  ou un reveur au milieu du desert'
	  	Marguerite Yourcenar
-----------------------------------------------------



This archive was generated by hypermail 2b29 : Tue Jan 02 2001 - 11:50:23 MET