Hi,
In your previous email you said that your ROOTSYS is:
export ROOTSYS=/afs/cern.ch/sw/root/v3.04.02/rh72_gcc2953/root;
Try to use 3.05/05 instead. I have recently modified that "contour
business" quite a lot. There was a lot of problems (with none equidistant
contours in particular). Note that the non equidistant level are also
implemented with the COLZ option now.
Olivier
On Wed, 18 Jun 2003, Matthieu Lechowski wrote:
> Hello,
>
> Thank you for your answer.
>
> I think it does not work for me.
> When I retrieve the array contours (after the plot), I have
> 0.99
> 63.9933
> 126.997
> I have max=190, so the levels are not 0.1*max, 0.2*max and 0.7*max, they
> are (almost) equidistant !
> Very strange...
> So it would be because of the ROOT version or the machine.
>
> Matthieu.
>
> On Wed, 18 Jun 2003, Olivier Couet wrote:
>
> >
> >
> > Can you try this one:
> >
> > {
> > TH2F *h = new TH2F("h","test contours",60,-3,3,60,-3,3);
> > for (Int_t i=0;i<100000;i++) {
> > double x,y;
> > gRandom->Rannor(x,y);
> > h->Fill(x,y);
> > }
> > h->SetMinimum(0.99);
> > Double_t max = h->GetMaximum();
> >
> > double contours[3];
> > contours[0] = 0.1*max;
> > contours[1] = 0.2*max;
> > contours[2] = 0.7*max;
> >
> > h->SetContour(3,contours);
> > h->Draw("cont1");
> >
> > }
> >
> > it should produce 3 non equidistant contours. At least it does on my Linux
> > RedHat PC using root 3.05/05
> >
> > Cheers,
> >
> > O.Couet
> >
> >
> >
> >
> > On Wed, 18 Jun 2003, Matthieu Lechowski wrote:
> >
> > > Hello,
> > >
> > > I do not succeed specifying my own levels of contour for 2D histograms.
> > >
> > > In a ROOT script (I load a program .c then I call a function), I do:
> > > TH2F* contour=new TH2F("contour","contour",nbinsX,xbins,nbinsY,ybins);
> > > //...filling "contour"
> > > Double_t levels[17];
> > > levels[0]=1.864
> > > levels[1]=5.
> > > //...filling "levels" non-equidistant
> > > contour->SetContour(17,levels); // <<----- ######
> > > // NB: definition of this function in TH1 class
> > > contour->Draw("cont1");
> > >
> > > My plot has 17 levels, OK, but they are equidistant !
> > >
> > > When I want be sure that ROOT has understood "levels", I do (in my
> > > function too):
> > > Double_t levels2[50];
> > > int ncontours=contour->GetContour(levels2);
> > > for(int i=0;i<ncontours;++i)
> > > cout<<levels2[i]<<endl;
> > > and I have exactly the same values that "levels", so it seems OK.
> > >
> > > The funny thing is that, after having done all that in the function, I
> > > redo in interactive the check with "levels2", and I have not the
> > > values of "levels", but the values of 17 equidistant levels (so as it
> > > is drawn)!
> > >
> > > So what happens ??
> > > Have anybody an idea?
> > > (I hope that I gave a not too dark explanation)
> > > Thank you.
> > >
> > > Matthieu Lechowski.
> > >
> > > PS: I'm using ROOT on lxplus at cern, with environment:
> > > export ROOTSYS=/afs/cern.ch/sw/root/v3.04.02/rh72_gcc2953/root;
> > > export PATH=$ROOTSYS/bin:$PATH;
> > > export LD_LIBRARY_PATH=$ROOTSYS/lib:$LD_LIBRARY_PATH
> > >
> > >
> >
> > --
> > Org: CERN - European Laboratory for Particle Physics.
> > Mail: 1211 Geneve 23 - Switzerland Mailbox: J25910
> > E-Mail: Olivier.Couet@cern.ch Phone: +41 22 7676522
> > WWW: http://cern.ch/Olivier.Couet/ Fax: +41 22 7677155
> >
> >
>
--
Org: CERN - European Laboratory for Particle Physics.
Mail: 1211 Geneve 23 - Switzerland Mailbox: J25910
E-Mail: Olivier.Couet@cern.ch Phone: +41 22 7676522
WWW: http://cern.ch/Olivier.Couet/ Fax: +41 22 7677155
This archive was generated by hypermail 2b29 : Thu Jan 01 2004 - 17:50:12 MET