Changeset 9848 for trunk/Mars/mjtrain


Ignore:
Timestamp:
08/11/10 17:14:28 (14 years ago)
Author:
tbretz
Message:
Some small changes to MJTrainDisp.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/mjtrain/MJTrainDisp.cc

    r9824 r9848  
    132132    //gPad->SetRightMargin(0.12);
    133133
    134     const Float_t cutval = hist.GetYaxis()->GetBinLowEdge(4);
     134    const Int_t bin = 3;
     135    const Float_t cutval = hist.GetYaxis()->GetBinUpEdge(bin);
    135136
    136137    TH1D heff;
     
    143144    MH::SetBinning(&heff, hist.GetXaxis());
    144145
    145 
    146146    for (int x=0; x<=hist.GetNbinsX()+1; x++)
    147147    {
    148148        const Double_t n0 = hist.Integral(x, x, -1, -1);
    149149        if (n0>0)
    150             heff.SetBinContent(x, hist.Integral(x, x, -1, 3)/n0);
     150            heff.SetBinContent(x, hist.Integral(x, x, -1, bin)/n0);
    151151    }
    152152
     
    155155    heff.DrawCopy();
    156156
    157     line.DrawLine(10, 0.5, 31623, 0.5);
     157    line.DrawLine(hist.GetXaxis()->GetXmin(), 0.5,
     158                  hist.GetXaxis()->GetXmax(), 0.5);
    158159
    159160    c.cd(i+0);
     
    169170    for (int x=0; x<=hist.GetNbinsX(); x++)
    170171    {
    171         Float_t n = 0;
    172         for (int y=1; y<=2; y++)
    173             n += hist.GetBinContent(x,y);
    174 
     172        const Float_t n = hist.GetBinContent(1); // + hist.GetBinContent(2)
    175173        if (n==0)
    176174            continue;
    177175
    178176        for (int y=0; y<=hist.GetNbinsY(); y++)
    179             hist.SetBinContent(x, y, 200*hist.GetBinContent(x,y)/n);
     177            hist.SetBinContent(x, y, 100*hist.GetBinContent(x,y)/n);
    180178    }
    181179
     
    183181    hist.DrawCopy("colz");
    184182
    185     line.DrawLine(10, 0.04, 31623, 0.04);
     183    line.DrawLine(hist.GetXaxis()->GetXmin(), cutval,
     184                  hist.GetXaxis()->GetXmax(), cutval);
    186185
    187186    c.cd(i+2);
     
    197196    for (int x=0; x<=hist.GetNbinsX(); x++)
    198197    {
    199         Float_t n = 0;
    200         for (int y=0; y<=hist.GetNbinsY(); y++)
    201             n += hist.GetBinContent(x,y);
    202 
     198        const Float_t n = hist.Integral(x, x, -1, 9999);
    203199        if (n==0)
    204200            continue;
     
    211207    hist.DrawCopy("colz");
    212208
    213     line.DrawLine(10, 0.04, 31623, 0.04);
     209    line.DrawLine(hist.GetXaxis()->GetXmin(), cutval,
     210                  hist.GetXaxis()->GetXmax(), cutval);
    214211}
    215212
Note: See TracChangeset for help on using the changeset viewer.