Changeset 1543 for trunk


Ignore:
Timestamp:
10/16/02 10:45:18 (22 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r1542 r1543  
    2727     - removed the second empty line after Print
    2828     - added Getter-function for the DAQEvtNumber
     29
     30   * macros/star.C, macros/starplot.C, macros/threshold.C, macros/trigrate.C,
     31     macros/plot.C, macros/plot2.C, macros/MagicHillas.C, macros/CT1Hillas.C,
     32     macros/collarea.C:
     33     - added ProgressBar
     34
     35   * macros/flux.C:
     36     - fixed some coding bugs (redeclaration of Draw, Draw starting with a
     37       capital)
     38
     39   * macros/MagicHillas.C:
     40     - fixed some typos
    2941
    3042

  • trunk/MagicSoft/Mars/NEWS

    r1540 r1543  
    7979   - Changed the image cleaning so that pixels with to many 'used'
    8080     neighbors are left used (to get rid of 'holes' in events)
     81
     82   - Introduced more output to the camera displays
    8183
    8284

  • trunk/MagicSoft/Mars/macros/CT1Hillas.C

    r1490 r1543  
    159159    // Execute your analysis
    160160    //
     161    MProgressBar bar;
     162    evtloop.SetProgressBar(&bar);
    161163    if (!evtloop.Eventloop())
    162164        return;

  • trunk/MagicSoft/Mars/macros/MagicHillas.C

    r1490 r1543  
    107107    //
    108108
    109     MCameraSmooth   smooth;
    110109    MImgCleanStd    clean;
    111110    MBlindPixelCalc blind;
    … …  
    121120
    122121    //
    123     // Uncomment this two line if you want to use MHillasExt instead
    124     // of MHillas
     122    // To use MHillasExt instead of MHillas
    125123    //
    126     //MHillasExt hext;
    127     //plist.AddToList(&hext);
     124    MHillasExt hext;
     125    plist.AddToList(&hext);
    128126
    129127    MFillH hfill1("MHHillas", "MHillas");
    … …  
    131129    MFillH hfill3("MHStarMap", "MHillas");
    132130    MFillH hfill4("HistExtSource [MHHillasExt]", "HillasSource");
    133     MFillH hfill5("HistExtAntiSource [MHHillasExt]", "HillasAntiSource");
     131    MFillH hfill5("HistExtAntiSrc [MHHillasExt]", "HillasAntiSrc");
    134132    MFillH hfill6("HistSource [MHHillasSrc]", "HillasSource");
    135133    MFillH hfill7("HistAntiSrc [MHHillasSrc]", "HillasAntiSrc");
    … …  
    141139    write.AddContainer("MMcEvt","Hillas");
    142140    write.AddContainer("Source","RunHeaders");
    143     write.AddContainer("AntiSource","RunHeaders");
     141    write.AddContainer("AntiSrc","RunHeaders");
    144142
    145143    tlist.AddToList(&read);
    … …  
    172170    // Execute your analysis
    173171    //
     172    MProgressBar bar;
     173    evtloop.SetProgressBar(&bar);
    174174    if (!evtloop.Eventloop())
    175175        return;
    … …  
    184184    plist.FindObject("MHStarMap")->DrawClone();
    185185    plist.FindObject("HistSource")->DrawClone();
    186     plist.FindObject("HistAntiSource")->DrawClone();
     186    plist.FindObject("HistAntiSrc")->DrawClone();
    187187    plist.FindObject("HistExtSource")->DrawClone();
    188     plist.FindObject("HistExtAntiSource")->DrawClone();
     188    plist.FindObject("HistExtAntiSrc")->DrawClone();
    189189}
    190190

  • trunk/MagicSoft/Mars/macros/collarea.C

    r1301 r1543  
    5353    // Start to loop over all events
    5454    //
     55    MProgressBar bar;
     56    magic.SetProgressBar(&bar);
    5557    if (!magic.Eventloop())
    5658        return;

  • trunk/MagicSoft/Mars/macros/flux.C

    r1417 r1543  
    420420     //                         and for different intervals in Theta     
    421421
    422      char *varnametime;
    423      varnametime = "Time";
    424      char *unittime;
    425      unittime    =  "[ s]";
    426      MHEffOnTime effontime(varnametime, unittime);
    427 
    428      char *varnametheta;
    429      varnametheta = "Theta";
    430      char *unittheta;
    431      unittheta    =  "[  \\circ]";
    432      MHEffOnTime effontheta(varnametheta, unittheta);
     422     MHEffOnTime effontime ("Time",  "[s]");
     423     MHEffOnTime effontheta("Theta", "[\\circ]");
    433424
    434425     effontime.SetupFill(&parlist);
    … …  
    438429     // Draw == 0  don't draw the individual distributions of time differences
    439430     //      != 0  draw them
    440      Bool_t Draw=kFALSE;
    441      effontime.Calc (dtimetime->GetHist(), Draw);
    442      effontheta.Calc(dtimetheta->GetHist(),Draw);
     431     Bool_t draw=kFALSE;
     432     effontime.Calc (dtimetime->GetHist(), draw);
     433     effontheta.Calc(dtimetheta->GetHist(),draw);
    443434
    444435
    … …  
    487478     MHGamma gamma;
    488479
    489      Bool_t Draw = kFALSE;
    490480     TH3D *hsubtime  = gamma.Subtract( evtsptime->GetHist(),
    491481                                       evtasptime->GetHist(),
    492                         "Al-En-time", "3D-plot of Alpha,E-est,time", Draw);   
     482                        "Al-En-time", "3D-plot of Alpha,E-est,time", draw);
    493483
    494484     TH3D *hsubtheta = gamma.Subtract( evtsptheta->GetHist(),
    495485                                       evtasptheta->GetHist(),
    496                         "Al-En-time", "3D-plot of Alpha,E-est,Theta", Draw);
     486                        "Al-En-time", "3D-plot of Alpha,E-est,Theta", draw);
    497487   
    498488
    … …  
    504494     Axis_t lo =   0;    // [deg]
    505495     Axis_t up =  10;    // [deg]
    506      Bool_t Draw = kFALSE;
    507496     const TH2D &evttime  = *(gamma.GetAlphaProjection(hsubtime,
    508                               lo, up, Draw));
     497                              lo, up, draw));
    509498
    510499     const TH2D &evttheta = *(gamma.GetAlphaProjection(hsubtheta,
    511                               lo, up, Draw));
     500                              lo, up, draw));
    512501     
    513502
    … …  
    547536     // Draw == kTRUE  draw the no.of photons vs. E-est
    548537     //                for the individual bins of the variable Var
    549      Bool_t Draw=kTRUE;
    550      MHFlux fluxtime(evttime, Draw, "Time", "  [  s]");
    551      fluxtime.Unfold(Draw);
     538     draw=kTRUE;
     539     MHFlux fluxtime(evttime, draw, "Time", "[s]");
     540     fluxtime.Unfold(draw);
    552541     fluxtime.CalcFlux(effontime.GetHist(), bartime.GetHist(),
    553                        &aeff, Draw);
     542                       &aeff, draw);
    554543     
    555544     fluxtime.DrawClone();
    … …  
    559548     //..............................................
    560549     // get flux spectrum for different bins in Theta
    561      Bool_t Draw=kTRUE;
    562      MHFlux fluxtheta(evttheta, Draw, "Theta", "  [  \\circ]");
    563      fluxtheta.Unfold(Draw);
     550     MHFlux fluxtheta(evttheta, draw, "Theta", "[\\circ]");
     551     fluxtheta.Unfold(draw);
    564552     fluxtheta.CalcFlux(effontheta.GetHist(), bartheta.GetHist(),
    565                        &aeff, Draw);
     553                       &aeff, draw);
    566554     
    567555     fluxtheta.DrawClone();

  • trunk/MagicSoft/Mars/macros/plot.C

    r1410 r1543  
    9999    // Execute your analysis
    100100    //
     101    MProgressBar bar;
     102    evtloop.SetProgressBar(&bar);
    101103    if (!evtloop.Eventloop())
    102104        return;

  • trunk/MagicSoft/Mars/macros/plot2.C

    r1410 r1543  
    104104    // Execute your analysis
    105105    //
     106    MProgressBar bar;
     107    evtloop.SetProgressBar(&bar);
    106108    if (!evtloop.Eventloop())
    107109        return;

  • trunk/MagicSoft/Mars/macros/star.C

    r1491 r1543  
    119119    // Create and setup the eventloop
    120120    //
     121    MProgressBar bar;
     122
    121123    MEvtLoop evtloop;
     124    evtloop.SetProgressBar(&bar);
    122125    evtloop.SetParList(&plist);
    123126

  • trunk/MagicSoft/Mars/macros/starplot.C

    r1490 r1543  
    105105    // Execute your analysis
    106106    //
     107    MProgressBar bar;
     108    evtloop.SetProgressBar(&bar);
    107109    if (!evtloop.Eventloop())
    108110        return;

  • trunk/MagicSoft/Mars/macros/threshold.C

    r1008 r1543  
    9090    // don't try to draw the results
    9191    //
     92    MProgressBar bar;
     93    evtloop.SetProgressBar(&bar);
    9294    if (!evtloop.Eventloop())
    9395        return;

  • trunk/MagicSoft/Mars/macros/trigrate.C

    r1380 r1543  
    2929!
    3030\* ======================================================================== */
     31
     32Int_t GetNSBEvents(TString name, Float_t *BgR)
     33{
     34    Int_t numnsbevents;
     35
     36    TChain events("Events");
     37    events.Add(name);
     38
     39    TH1F h("h","",5,.5,5.5);
     40
     41    for (Int_t i = from; i <= to; i++)
     42    {
     43        TString plot = "MMcTrig;";
     44        plot+=i;
     45        plot+=".fNumFirstLevel>>h";
     46
     47        if (i == 0)
     48            plot = "MMcTrig.fNumFirstLevel>>h";
     49
     50        events.Draw(plot,"", "goff");
     51
     52        // Get total number of L1 triggers from histogram:
     53        BgR[dim>0? i-1: 0] = h.Integral()*h.GetMean();
     54
     55        numnsbevents = (Float_t) h.GetEntries();
     56    }
     57
     58    return numnsbevents;
     59}
     60
    3161
    3262void trigrate(int dim=0, char *filename = "data/camera.root",
    … …  
    97127    tasklist.AddToList(&reader);
    98128
    99  
    100129    // Now we have to build the BgR array, containing the number
    101130    // of triggers (may be more than 1 trigger/event!) from the
    … …  
    103132    // If no nsbfile is supplied, we assume no triggers from NSB
    104133
     134    Float_t* BgR = new Float_t[num];
     135    memset(BgR, 0, num*sizeof(Float_t));
     136
    105137    Float_t numnsbevents = 5.e4; // some default value.
    106 
    107     Float_t* BgR = new float[num];
    108     for(Int_t i = 0; i < num; i++)
    109       BgR[i] = 0.;
    110 
    111138    if (nsbfile)
    112       {
    113         TChain nsb_events("Events");
    114         nsb_events.Add(nsbfile);
    115 
    116         TH1F h("h","",5,.5,5.5);
    117 
    118         for (Int_t i = from; i <= to; i++)
    119           {
    120             TString plot = "MMcTrig;";
    121             plot+=i;
    122             plot+=".fNumFirstLevel>>h";
    123 
    124             if (i == 0)
    125               plot = "MMcTrig.fNumFirstLevel>>h";
    126 
    127             nsb_events.Draw(plot,"", "goff");
    128             // Get total number of L1 triggers from histogram:
    129             BgR[dim>0? i-1: 0] = h.Integral()*h.GetMean();
    130 
    131             numnsbevents = (Float_t) h.GetEntries();
    132           }
    133       }
     139        numnsbevents = GetNSBEvents(nsbfile, BgR);
    134140
    135141    cout << "Number of Trigger conditions: " << num << endl;
    … …  
    147153    // Start to loop over all events
    148154    //
     155    MProgressBar bar;
     156    magic.SetProgressBar(&bar);
    149157    if (!magic.Eventloop())
    150158        return;
    151159
     160    delete BgR;
     161
    152162    hists.Print();
    153163}
Note: See TracChangeset for help on using the changeset viewer.