Ignore:
Timestamp:
06/10/08 16:50:49 (16 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/datacenter/macros/fillstar.C

    r8931 r8934  
    183183        return 2;
    184184    }
    185 
    186     Double_t max = TMath::MaxElement(g->GetN(), g->GetY());
    187     TString maxhum = Form("%5.1f", max);
     185    TString maxhum = Form("%5.1f", TMath::MaxElement(g->GetN(), g->GetY()));
     186    TString avghum = Form("%5.1f", g->GetMean(2));
    188187
    189188    g = (TGraph*)arr.FindObjectInCanvas("Temperature", "TGraph", "MHWeather");
     
    193192        return 2;
    194193    }
    195 
    196     Double_t temp = g->GetMean(2);
    197     TString avgtemp = Form("%5.1f", temp);
     194    TString avgtemp = Form("%5.1f", g->GetMean(2));
    198195
    199196    g = (TGraph*)arr.FindObjectInCanvas("WindSpeed", "TGraph", "MHWeather");
     
    203200        return 2;
    204201    }
    205 
    206     Double_t wind = g->GetMean(2);
    207     TString avgwind = Form("%5.1f", wind);
     202    TString avgwind = Form("%5.1f", g->GetMean(2));
     203
     204    g = (TGraph*)arr.FindObjectInCanvas("Cloudiness", "TGraph", "MHWeather");
     205    if (!g)
     206        cout << "WARNING - Reading of Cloudiness failed." << endl;
     207    TString avgclouds = g ? Form("%5.1f", g->GetMean(2)) : 0;
     208    TString rmsclouds = g ? Form("%5.1f", g->GetRMS(2))  : 0;
     209
     210    g = (TGraph*)arr.FindObjectInCanvas("TempSky", "TGraph", "MHWeather");
     211    if (!g)
     212        cout << "WARNING - Reading of TempSky failed." << endl;
     213    TString avgsky = g ? Form("%5.1f", g->GetMean(2)+200) : 0;
    208214
    209215
    210216    g = (TGraph*)arr.FindObjectInCanvas("NumStars", "TGraph", "MHPointing");
    211217    if (!g)
    212     {
    213218        cout << "WARNING - Reading of NumStars failed." << endl;
    214 //        return 2;
    215     }
    216219
    217220    Double_t numstarmed = g ? TMath::Median(g->GetN(), g->GetY()) : -1;
     
    227230    }
    228231
    229     Double_t numcormed = g ? TMath::Median(g->GetN(), g->GetY()) : -1;
    230     TString numcorsmed = Form("%5.1f", numcormed);
    231     Double_t numcorrms = g ? g->GetRMS(2) : -1;
    232     TString numcorsrms = Form("%5.1f", numcorrms);
     232    TString numcorsmed = Form("%5.1f", g ? TMath::Median(g->GetN(), g->GetY()) : -1);
     233    TString numcorsrms = Form("%5.1f", g ? g->GetRMS(2) : -1);
    233234
    234235    g = (TGraph*)arr.FindObjectInCanvas("Brightness", "TGraph", "MHPointing");
     
    239240    }
    240241
    241     Double_t brightnessmed = g ? TMath::Median(g->GetN(), g->GetY()) : -1;
    242     TString skybrightnessmed = Form("%5.1f", brightnessmed);
    243     Double_t brightnessrms = g ? g->GetRMS(2) : -1;
    244     TString skybrightnessrms = Form("%5.1f", brightnessrms);
     242    TString skybrightnessmed = Form("%5.1f", g ? TMath::Median(g->GetN(), g->GetY()) : -1);
     243    TString skybrightnessrms = Form("%5.1f", g ? g->GetRMS(2) : -1);
    245244
    246245
     
    267266    cout << "  Rate of sparks [Hz]      " << sparkrate << endl;
    268267    cout << "  Maximum Humidity [%]     " << maxhum    << endl;
     268    cout << "  Average Humidity [%]     " << avghum    << endl;
     269    cout << "  Average WindSpeed [km/h] " << avgwind   << endl;
    269270    cout << "  Average Temp [°C]        " << avgtemp   << endl;
    270     cout << "  Average WindSpeed [km/h] " << avgwind   << endl;
    271     cout << "  Number of Stars          " << numstarsmed      << " +/- " << numstarsrms       << endl;
    272     cout << "  Number of cor. Stars     " << numcormed        << " +/- " << numcorrms         << endl;
    273     cout << "  Skybrightness            " << skybrightnessmed << " +/- " << skybrightnessrms  << endl;
     271    cout << "  Average Sky Temp [K]     " << avgsky    << endl;
     272    cout << "  Cloundiness [%]          " << avgclouds        << " +/- " << rmsclouds        << endl;
     273    cout << "  Number of Stars          " << numstarsmed      << " +/- " << numstarsrms      << endl;
     274    cout << "  Number of cor. Stars     " << numcorsmed       << " +/- " << numcorsrms       << endl;
     275    cout << "  Skybrightness            " << skybrightnessmed << " +/- " << skybrightnessrms << endl;
    274276
    275277    TString vars = Form(" fMeanNumberIslands=%s,"
     
    282284                        " fSparkRate=%s,"
    283285                        " fMaxHumidity=%s,"
     286                        " fAvgHumidity=%s, "
    284287                        " fAvgTemperature=%s,"
    285288                        " fAvgWindSpeed=%s,"
     289                        " fAvgTempSky=%s,"
     290                        " fAvgCloudiness=%s, "
     291                        " fRmsCloudiness=%s, "
    286292                        " fNumStarsMed=%s,"
    287293                        " fNumStarsRMS=%s,"
     
    295301                        muonrate.Data(), PSF.Data(),
    296302                        datarate, sparkrate.Data(), maxhum.Data(),
    297                         avgtemp.Data(), avgwind.Data(),
     303                        avghum.Data(), avgtemp.Data(), avgwind.Data(),
     304                        avgsky.Data(), avgclouds.Data(), rmsclouds.Data(),
    298305                        numstarsmed.Data(), numstarsrms.Data(),
    299306                        numcorsmed.Data(), numcorsrms.Data(),
Note: See TracChangeset for help on using the changeset viewer.