Changeset 9496


Ignore:
Timestamp:
08/17/09 12:46:45 (15 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Cosy/tpoint
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Cosy/tpoint/plot_m1.C

    r9479 r9496  
    5858    ifstream fin(fname);
    5959
    60     cout << "Reading " << setw(23) << fname << "..." << flush;
     60    cout << "Reading " << setw(23) << gSystem->BaseName(fname) << "..." << flush;
    6161
    6262    while (1)
    … …  
    101101};
    102102
    103 const Int_t counts = 29+10+18+1+13+11;
     103const Int_t counts = 29+10+18+1+13+11+31;
    104104Description_t desc[counts] =
    105105{
    … …  
    229229    {"090517", "TPoints 09/05/17" ,"tpoint/m1/2009_05_17/tpoint_20090516_234825.txt"},
    230230    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_05_17/tpoint_20090517_023340.txt"},
    231 
    232231    // New period
    233232    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_01/tpoint_20090531_215148.txt"},
    234233    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_01/tpoint_20090531_222549.txt"},
    235 
    236234    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_02/tpoint_20090601_223009.txt"},
    237 
    238235    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_03/tpoint_20090602_213509.txt"},
    239236    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_03/tpoint_20090603_011936.txt"},
    240 
    241237    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_04/tpoint_20090603_215840.txt"},
    242238    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_04/tpoint_20090603_230510.txt"},
    243 
    244239    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_05/tpoint_20090604_215943.txt"},
    245240    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_05/tpoint_20090604_232320.txt"},
    246 
     241    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_10/tpoint_20090609_231948.txt"},
     242    {"+090517","TPoints 09/05/17", "tpoint/m1/2009_06_11/tpoint_20090611_011155.txt"},
     243
     244    // ------- Something happened (reason unknown) --------
     245
     246    {"090613","TPoints 09/05/17", "tpoint/m1/2009_06_11/tpoint_20090611_023625.txt"},
     247
     248    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_06_13/tpoint_20090613_033838.txt"},
     249    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_06_14/tpoint_20090614_021257.txt"},
     250    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_06_18/tpoint_20090618_041433.txt"},
     251
     252    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_02/tpoint_20090701_215304.txt"},
     253    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_02/tpoint_20090701_222059.txt"},
     254    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_02/tpoint_20090701_224051.txt"},
     255    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_02/tpoint_20090701_225615.txt"},
     256    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_02/tpoint_20090701_230946.txt"},
     257
     258    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_03/tpoint_20090702_225940.txt"},
     259    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_04/tpoint_20090703_224721.txt"},
     260    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_05/tpoint_20090705_012638.txt"},
     261
     262    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_12/tpoint_20090712_014300.txt"},
     263    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_12/tpoint_20090712_024710.txt"},
     264    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_13/tpoint_20090713_004241.txt"},
     265    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_13/tpoint_20090713_025934.txt"},
     266
     267    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_14/tpoint_20090714_024729.txt"},
     268    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_15/tpoint_20090715_021320.txt"},
     269    {"+090613","TPoints 09/05/17", "tpoint/m1/2009_07_15/tpoint_20090715_025237.txt"},
     270
     271    // ----- AMC adjust (23.7.) -----
     272
     273    {"090723","TPoints 09/05/17", "tpoint/m1/2009_07_24/tpoint_20090724_051142.txt"},
     274    {"+090723", "TPoints 08/2009" , "tpoint/m1/2009_08_02/tpoint_20090801_223848.txt"},
     275    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_03/tpoint_20090802_224434.txt"},
     276    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_04/tpoint_20090803_230214.txt"},
     277    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_10/tpoint_20090809_233524.txt"},
     278    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_10/tpoint_20090810_014642.txt"},
     279    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_10/tpoint_20090810_020616.txt"},
     280    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_11/tpoint_20090810_232543.txt"},
     281    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_12/tpoint_20090812_001846.txt"},
     282    {"+090723","TPoints 08/2009" , "tpoint/m1/2009_08_13/tpoint_20090813_030807.txt"},
     283
     284    // ------ 09/08/17 new pointing model -------
    247285};
    248286
    … …  
    270308        cout << setw(2) << num << ": " << flush;
    271309        n[num] = fill(desc[i].fFile, g, &hx[num]);
    272         cout << " Mean: " << setw(5) << setprecision(2) << hx[num].GetMean() << " deg +/- " <<  hx[num].GetRMS()<< endl;
     310        cout << " Mean: " << setw(5) << setprecision(3) << hx[num].GetMean()*3600 << " sec +/- " <<  hx[num].GetRMS()*3600 << endl;
    273311    }
    274312
    … …  
    359397        MTime(2009,  5, 12).GetAxisTime(),
    360398        MTime(2009,  5, 13).GetAxisTime(),
    361         MTime(2009, 12, 30).GetAxisTime(),
     399        MTime(2009,  6, 11).GetAxisTime(),
     400        MTime(2009,  7, 23).GetAxisTime(),
     401        MTime(2009,  8, 17).GetAxisTime(),
     402        MTime(2009, 12, 31).GetAxisTime(),
    362403    };
    363404
    … …  
    365406
    366407    MBinning bins;
    367     bins.SetEdges(TArrayD(17, time));
     408    bins.SetEdges(TArrayD(20, time));
    368409
    369410    bins.Apply(histres[0]);
    … …  
    376417    for (int i=0; i<n.GetSize(); i++)
    377418    {
    378         cout << i+1 << " - Mean:  " << Form("%.4f +- %.4f", hx[i].GetMean(), hx[i].GetRMS());
     419        cout << i+1 << " - Mean:  " << Form("%5.1f +- %5.1f", hx[i].GetMean()*3600, hx[i].GetRMS()*3600);
    379420        cout << "   (Overflows=" << hx[i].GetBinContent(hx[i].GetNbinsX()+1)*hx[i].GetEntries() << ")   " << (int)hx[i].GetEntries() << endl;
    380421
    … …  
    473514
    474515    line.DrawLine(time[0], 0, time[0], histres[2].GetMaximum());
    475     for (int i=0; i<11; i++)
     516    for (int i=0; i<bins.GetNumBins()-1; i++)
    476517        line.DrawLine(time[i+1], 0, time[i+1], histres[2].GetBinContent(i+1));
    477518

  • trunk/MagicSoft/Cosy/tpoint/plot_m2.C

    r9479 r9496  
    101101};
    102102
    103 const Int_t counts = 27;
     103const Int_t counts = 46;
    104104Description_t desc[counts] =
    105105{
    … …  
    134134    {"+090512",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_04/tpoint_20090603_215841.txt"},
    135135    {"+090512",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_05/tpoint_20090604_215941.txt"},
     136
     137    {"+090512",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_10/tpoint_20090609_232320.txt"},
     138    {"+090512",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_11/tpoint_20090611_011148.txt"},
     139
     140    // Something happened (reason unknown)
     141
     142    {"090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_11/tpoint_20090611_023623.txt"},
     143    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_13/tpoint_20090613_033839.txt"},
     144    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_14/tpoint_20090614_021304.txt"},
     145    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_06_18/tpoint_20090618_041436.txt"},
     146    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_07_12/tpoint_20090712_014259.txt"},
     147    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_07_12/tpoint_20090712_024706.txt"},
     148    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_07_13/tpoint_20090713_004244.txt"},
     149    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_07_13/tpoint_20090713_025925.txt"},
     150    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_07_14/tpoint_20090714_022322.txt"},
     151    {"+090611",  "TPoints Residuals 11/2004" , "tpoint/m2/2009_07_15/tpoint_20090715_021316.txt"},
     152    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_02/tpoint_20090801_225305.txt"},
     153    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_03/tpoint_20090802_224437.txt"},
     154    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_04/tpoint_20090803_231853.txt"},
     155    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_11/tpoint_20090810_234604.txt"},
     156    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_11/tpoint_20090811_020844.txt"},
     157    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_12/tpoint_20090812_001850.txt"},
     158    {"+090611",  "TPoints Residuals 08/2009" , "tpoint/m2/2009_08_13/tpoint_20090813_030805.txt"},
     159
     160    // ------ 09/08/17 new pointing model -------
     161
    136162};
    137163
    … …  
    232258        MTime(2009,  4, 11).GetAxisTime(),
    233259        MTime(2009,  5, 11).GetAxisTime(),
    234         MTime(2009,  5, 31).GetAxisTime(),
    235         MTime(2009,  6,  6).GetAxisTime(),
     260        MTime(2009,  6, 11).GetAxisTime(),
     261        MTime(2009,  8, 17).GetAxisTime(),
     262        MTime(2009, 12, 31).GetAxisTime(),
    236263    };
    237264
    … …  
    239266
    240267    MBinning bins;
    241     bins.SetEdges(TArrayD(5, time));
     268    bins.SetEdges(TArrayD(6, time));
    242269
    243270    bins.Apply(histres[0]);
    … …  
    347374
    348375    line.DrawLine(time[0], 0, time[0], histres[2].GetMaximum());
    349     for (int i=0; i<bins.GetNumBins(); i++)
     376    for (int i=0; i<bins.GetNumBins()-1; i++)
    350377        line.DrawLine(time[i+1], 0, time[i+1], histres[2].GetBinContent(i+1));
    351378
    … …  
    379406        line.DrawLine(monl, -0.12, monl, 0.12);
    380407    }*/
    381     for (int m=1; m<13; m++)
     408    for (int m=4; m<13; m++)
    382409    {
    383410        Double_t monl = MTime(2009, m,   1,0).GetAxisTime();
Note: See TracChangeset for help on using the changeset viewer.