Changeset 3169 for trunk


Ignore:
Timestamp:
02/15/04 22:16:40 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r3159 r3169  
    1818   * manalysis/MPedestalPix.h
    1919     - removed include of MHPedestalPixel
     20
     21   * macros/calibration.C
     22     - changed pedestal displays accordingly
     23
    2024
    2125   * mjobs/MJPedestal.h
    … …  
    3337   * mcalib/MHGausEvent.[h,cc], mcalib/Makefile, mcalib/CalibLinkDef.h:
    3438     - replaced by the improved version: MHGausEvents.[h,cc]
     39
    3540
    3641

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

    r3124 r3169  
    3939{
    4040
     41    gStyle->SetOptStat(1111);
     42    gStyle->SetOptFit();
     43
    4144    MStatusDisplay *display = new MStatusDisplay;
    4245    display->SetUpdateTime(3000);
    … …  
    6366    MExtractSignal  sigcalc;
    6467    //
    65     // We saw that the signal jumps between slices,
    66     // thus set the extraction range as high as possible
     68    // Set the extraction range higher:
    6769    //         
    68     sigcalc.SetRange(1,14,1,14);
     70    //sigcalc.SetRange(1,14,1,14);
    6971
    7072    MPedCalcPedRun pedcalc;
    71 //    pedcalc.SetNumHiGainSamples(14);
    72     MFillH fill("MPedestalCam", "MExtractedSignalCam");
     73
     74    //
     75    // Additionally to calculating the pedestals,
     76    // you can fill histograms and look at them
     77    //
     78    MFillH fill("MHPedestalCam", "MExtractedSignalCam");
    7379
    7480    tlist.AddToList(&read);
    … …  
    8086    MGeomCamMagic  geomcam;
    8187    MPedestalCam   pedcam;
     88    MHPedestalCam  hpedcam;
    8289    plist.AddToList(&geomcam);
    8390    plist.AddToList(&pedcam);
     91    plist.AddToList(&hpedcam);
    8492
    8593    //
    … …  
    98106    tlist.PrintStatistics();
    99107
    100     //    pedcam(17).DrawClone();
     108    //
     109    // Look at one specific pixel, after all the histogram manipulations:
     110    //
     111    // hpedcam[17].DrawClone("fourierevents");
    101112
    102113    MHCamera dispped0  (geomcam, "Ped;Pedestal",               "Mean per Slice");
    … …  
    104115    MHCamera dispped2  (geomcam, "Ped;PedestalRms",            "RMS per Slice");
    105116    MHCamera dispped3  (geomcam, "Ped;PedestalRmsErr",         "RMS Error per Slice");
     117
    106118    MHCamera dispped4  (geomcam, "Ped;Mean",                   "Fitted Mean per Slice");
    107119    MHCamera dispped5  (geomcam, "Ped;MeanErr",                "Fitted Error of Mean per Slice");
    … …  
    110122    MHCamera dispped8  (geomcam, "Ped;Prob",                   "Probability of Fit");
    111123    MHCamera dispped9  (geomcam, "Ped;DeltaPedestalMean",      "Rel. Diff. Mean per Slice (Calc.-Fitte)");
    112     MHCamera dispped11 (geomcam, "Ped;DeltaPedestalMeanError", "Rel. Diff. Mean Error per Slice (Calc.-Fitted)");
    113     MHCamera dispped12  (geomcam, "Ped;DeltaRmsSigma", "Rel. Diff. RMS per Slice (Calc.-Fitted)");
    114     MHCamera dispped14  (geomcam, "Ped;DeltaRmsSigmaError", "Rel. Diff. RMS Error per Slice (Calc.-Fitted)");
    115 
    116     dispped0.SetCamContent(pedcam, 0);
    117     dispped0.SetCamError(pedcam, 1);
    118     dispped1.SetCamContent(pedcam, 1);
    119     dispped2.SetCamContent(pedcam, 2);
    120     dispped2.SetCamError(pedcam,3);
    121     dispped3.SetCamContent(pedcam, 3);
    122     dispped4.SetCamContent(pedcam, 4);
    123     dispped4.SetCamError(pedcam, 5);
    124     dispped5.SetCamContent(pedcam, 5);
    125     dispped6.SetCamContent(pedcam, 6);
    126     dispped6.SetCamError(pedcam, 7);
    127     dispped7.SetCamContent(pedcam, 7);
    128     dispped8.SetCamContent(pedcam, 8);
    129     dispped9.SetCamContent(pedcam, 9);
    130     dispped9.SetCamError(pedcam, 10);
    131     dispped11.SetCamContent(pedcam, 11);
    132     dispped12.SetCamContent(pedcam, 12);
    133     dispped12.SetCamError(pedcam, 13);
    134     dispped14.SetCamContent(pedcam, 14);
     124    MHCamera dispped10 (geomcam, "Ped;DeltaPedestalMeanError", "Rel. Diff. Mean Error per Slice (Calc.-Fitted)");
     125    MHCamera dispped11  (geomcam, "Ped;DeltaRmsSigma",         "Rel. Diff. RMS per Slice (Calc.-Fitted)");
     126    MHCamera dispped12  (geomcam, "Ped;DeltaRmsSigmaError",    "Rel. Diff. RMS Error per Slice (Calc.-Fitted)");
     127    MHCamera dispped13  (geomcam, "Ped;FitOK",                 "Gaus Fit not OK");
     128    MHCamera dispped14  (geomcam, "Ped;FourierOK",             "Fourier Analysis not OK");
     129
     130    dispped0.SetCamContent(  pedcam, 0);
     131    dispped0.SetCamError(    pedcam, 1);
     132    dispped1.SetCamContent(  pedcam, 1);
     133    dispped2.SetCamContent(  pedcam, 2);
     134    dispped2.SetCamError(    pedcam, 3);
     135    dispped3.SetCamContent(  pedcam, 3);
     136
     137    dispped4.SetCamContent( hpedcam, 0);
     138    dispped4.SetCamError(   hpedcam, 1);
     139    dispped5.SetCamContent( hpedcam, 1);
     140    dispped6.SetCamContent( hpedcam, 2);
     141    dispped6.SetCamError(   hpedcam, 3);
     142    dispped7.SetCamContent( hpedcam, 3);
     143    dispped8.SetCamContent( hpedcam, 4);
     144    dispped9.SetCamContent( hpedcam, 5);
     145    dispped9.SetCamError(   hpedcam, 6);
     146    dispped10.SetCamContent(hpedcam, 7);
     147    dispped11.SetCamContent(hpedcam, 8);
     148    dispped11.SetCamError(  hpedcam, 9);
     149    dispped12.SetCamContent(hpedcam, 10);
     150    dispped13.SetCamContent(hpedcam, 11);
     151    dispped14.SetCamContent(hpedcam, 12);
    135152
    136153    dispped0.SetYTitle("Calc. Pedestal per slice [FADC counts]");
    … …  
    144161    dispped8.SetYTitle("Fit Probability [1]");
    145162    dispped9.SetYTitle("Rel. Diff. Pedestal Calc.-Fitted per slice [1]");
    146     dispped11.SetYTitle("Rel. Diff. Pedestal Error Calc.-Fitted per slice [1]");
    147     dispped12.SetYTitle("Rel. Diff. Pedestal RMS Calc.-Fitted per slice [1]");
    148     dispped14.SetYTitle("Rel. Diff. Pedestal RMS Error Calc.-Fitted per slice [1]");
    149 
    150     gStyle->SetOptStat(1111);
    151     gStyle->SetOptFit();
     163    dispped10.SetYTitle("Rel. Diff. Pedestal Error Calc.-Fitted per slice [1]");
     164    dispped11.SetYTitle("Rel. Diff. Pedestal RMS Calc.-Fitted per slice [1]");
     165    dispped12.SetYTitle("Rel. Diff. Pedestal RMS Error Calc.-Fitted per slice [1]");
     166    dispped13.SetYTitle("[1]");
     167    dispped14.SetYTitle("[1]");
    152168   
    153169    // Histogram values
    … …  
    155171    b1.Divide(4,3);
    156172
    157     CamDraw(b1,dispped0,pedcam,1,4,0);
    158     CamDraw(b1,dispped1,pedcam,2,4,1);
    159     CamDraw(b1,dispped2,pedcam,3,4,1);
    160     CamDraw(b1,dispped3,pedcam,4,4,1);
     173    CamDraw(b1,dispped0,pedcam,1,4,1);
     174    CamDraw(b1,dispped1,pedcam,2,4,2);
     175    CamDraw(b1,dispped2,pedcam,3,4,2);
     176    CamDraw(b1,dispped3,pedcam,4,4,2);
    161177
    162178    // Fitted values
    … …  
    164180    b2.Divide(4,3);
    165181
    166     CamDraw(b2,dispped4,pedcam,1,4,0);
    167     CamDraw(b2,dispped5,pedcam,2,4,1);
    168     CamDraw(b2,dispped6,pedcam,3,4,1);
    169     CamDraw(b2,dispped7,pedcam,4,4,1);
     182    CamDraw(b2,dispped4,hpedcam,1,4,1);
     183    CamDraw(b2,dispped5,hpedcam,2,4,2);
     184    CamDraw(b2,dispped6,hpedcam,3,4,2);
     185    CamDraw(b2,dispped7,hpedcam,4,4,2);
    170186
    171187
    … …  
    174190    b3.Divide(1,3);
    175191
    176     CamDraw(b3,dispped8,pedcam,1,1,3);
     192    CamDraw(b3,dispped8,hpedcam,1,1,3);
    177193
    178194    // Differences
    … …  
    180196    c4.Divide(4,3);
    181197
    182     CamDraw(c4,dispped9,pedcam,1,4,1);
    183     CamDraw(c4,dispped11,pedcam,2,4,1);
    184     CamDraw(c4,dispped12,pedcam,3,4,1);
    185     CamDraw(c4,dispped14,pedcam,4,4,1);
    186 
     198    CamDraw(c4,dispped9,hpedcam,1,4,1);
     199    CamDraw(c4,dispped10,hpedcam,2,4,1);
     200    CamDraw(c4,dispped11,hpedcam,3,4,1);
     201    CamDraw(c4,dispped12,hpedcam,4,4,1);
     202
     203    // Defects
     204    TCanvas &c5 = display->AddTab("Defects");
     205    c5.Divide(2,2);
     206
     207    CamDraw(c5,dispped13,hpedcam,1,2,0);
     208    CamDraw(c5,dispped14,hpedcam,2,2,0);
    187209
    188210    //
    … …  
    269291    //
    270292    MCalibrationBlindPix *bp = calcam.GetBlindPixel();
    271 //    bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPolya);
    272     bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson5);
     293    //  bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson5);
    273294
    274295    //
    … …  
    285306    // of the arrival times using a spline, uncomment the next line
    286307    //
    287     //    tlist2.AddToList(&timecalc);
     308    // tlist2.AddToList(&timecalc);
    288309    tlist2.AddToList(&cont);
    289310    tlist2.AddToList(&calcalc);
    … …  
    312333    // just one example how to get the plots of individual pixels
    313334    //
    314     //    calcam[563].DrawClone();
    315     //    calcam[564].DrawClone();
     335    //  calcam[563].DrawClone();
     336    //  calcam[564].DrawClone();
    316337
    317338    // Create histograms to display
    … …  
    353374    disp2.SetCamContent(calcam,  2);
    354375    disp2.SetCamError(  calcam,  3);
     376
    355377    // Fit probabilities
    356378    disp3.SetCamContent(calcam,  4);
    … …  
    547569  gPad->SetBorderMode(0);
    548570  MHCamera *obj1=(MHCamera*)cam.DrawCopy("hist");
    549   obj1->AddNotify(evt);
     571  //  obj1->AddNotify(evt);
    550572 
    551573  c.cd(i+j);
Note: See TracChangeset for help on using the changeset viewer.