Changeset 2217 for trunk


Ignore:
Timestamp:
06/23/03 15:00:09 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r2216 r2217  
    132132     - added -I../mhist
    133133
     134   * macros/sumcurrents.C:
     135     - adapted to changes
     136     - added filename to canvas
     137
    134138
    135139
  • trunk/MagicSoft/Mars/macros/sumcurrents.C

    r2192 r2217  
    4747    tlist.AddToList(&read);
    4848
    49     MFillH fill("MHCurrents");
     49    MFillH fill("MHCamEvent", "MCurrents");
    5050    tlist.AddToList(&fill);
    5151
     
    6464    tlist.PrintStatistics();
    6565
    66     MHCurrents &h = *(MHCurrents*)plist->FindObject("MHCurrents");
     66    MHCamEvent &h2 = *(MHCamEvent*)plist->FindObject("MHCamEvent");
     67    MHCamera &h = *(MHCamera*)h2.GetHistByName();
     68;
    6769
    6870    TCanvas *c = MH::MakeDefCanvas();
    6971    c->Divide(3, 2);
    7072
    71     MCamDisplay *disp1=new MCamDisplay(&geom);
    72     MCamDisplay *disp2=new MCamDisplay(&geom);
    73     MCamDisplay *disp3=new MCamDisplay(&geom);
    74     disp1->Fill(h.GetSum());
    75     disp2->Fill(h.GetRms());
     73    h.SetXTitle("Pixel Index");
    7674
    77     TArrayF arr(577);
    78     for (int i=0;i<577;i++)
    79     {
    80         TArrayF &r = h.GetRms();
    81         TArrayF &v = h.GetSum();
    82         arr[i] = (v[i]==0 ? 0 : r[i]/v[i]);
    83     }
    84     disp3->Fill(arr);
     75    MHCamera *disp1=h.Clone();
     76    MHCamera *disp2=h.Clone();
     77    MHCamera *disp3=h.Clone();
     78    disp2->SetCamContent(h, 1);
     79    disp3->SetCamContent(h, 2);
     80
     81    disp1->SetYTitle("I [nA]");
     82    disp2->SetYTitle("\\sigma_{I} [nA]");
     83    disp3->SetYTitle("\\sigma_{I} [%]");
     84    disp1->SetName("Currents;avg");
     85    disp2->SetName("Currents;err");
     86    disp3->SetName("Currents;rel");
     87    disp1->SetTitle("Currents Average");
     88    disp2->SetTitle("Currents error");
     89    disp3->SetTitle("Currents relative error");
    8590
    8691    c->cd(1);
     92    TText text(0.1, 0.95, &fname[fname.Last('/')+1]);
     93    text.SetTextSize(0.03);
     94    text.DrawClone();
    8795    gPad->SetBorderMode(0);
    8896    gPad->Divide(1,1);
    8997    gPad->cd(1);
    90     gPad->SetLogz();
     98    gPad->SetLogy();
    9199    disp1->Draw();
    92100    disp1->SetBit(kCanDelete);
     
    95103    gPad->Divide(1,1);
    96104    gPad->cd(1);
    97     gPad->SetLogz();
     105    gPad->SetLogy();
    98106    disp2->Draw();
    99107    disp2->SetBit(kCanDelete);
     
    102110    gPad->Divide(1,1);
    103111    gPad->cd(1);
    104     gPad->SetLogz();
     112    gPad->SetLogy();
    105113    disp3->Draw();
    106114    disp3->SetBit(kCanDelete);
    107115    c->cd(4);
    108116    gPad->SetBorderMode(0);
    109     h.GetHist().SetStats(kFALSE);
    110     h.GetHist().DrawCopy();
     117    disp1->Draw("EPhist");
    111118    c->cd(5);
    112119    gPad->SetBorderMode(0);
    113120    gPad->SetLogy();
    114     TH1F h1("currents;rms", "Currents Rms", 577, -0.5, 576.5);
    115     h1.SetXTitle("Pixel Index");
    116     h1.SetStats(kFALSE);
    117     for (int i=1;i<=577; i++)
    118         h1.SetBinContent(i, h.GetHist().GetBinError(i));
    119     h1.DrawCopy();
     121    disp2->Draw("Phist");
    120122    c->cd(6);
    121123    gPad->SetBorderMode(0);
    122124    gPad->SetLogy();
    123     TH1F h2("currents;relerr", "Currents rel. Error [%]", 577, -0.5, 576.5);
    124     h2.SetXTitle("Pixel Index");
    125     h2.SetStats(kFALSE);
    126     for (int i=1;i<=577; i++)
    127         h2.SetBinContent(i, arr[i-1]*100);
    128     h2.DrawCopy();
     125    disp3->Draw("Phist");
    129126
    130127    c->SaveAs(fname(0, fname.Last('.')+1) + "ps");
Note: See TracChangeset for help on using the changeset viewer.