Ignore:
Timestamp:
02/16/04 12:32:04 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r3169 r3185  
    6565    MGeomApply      geomapl;
    6666    MExtractSignal  sigcalc;
     67
     68
     69    MExtractPINDiode pidcalc;
     70    MExtractedSignalPINDiode pindiode;
     71    //    MHCalibrationChargePINDiode pinhist;
    6772    //
    6873    // Set the extraction range higher:
    … …  
    7782    //
    7883    MFillH fill("MHPedestalCam", "MExtractedSignalCam");
     84    //    MFillH fillp("MHCalibrationChargePINDiode", "MExtractedSignalPINDiode");
    7985
    8086    tlist.AddToList(&read);
    8187    tlist.AddToList(&geomapl);
    8288    tlist.AddToList(&sigcalc);
     89    tlist.AddToList(&pidcalc);
    8390    tlist.AddToList(&pedcalc);
    8491    tlist.AddToList(&fill);
     92    //    tlist.AddToList(&fillp);
    8593
    8694    MGeomCamMagic  geomcam;
    … …  
    9098    plist.AddToList(&pedcam);
    9199    plist.AddToList(&hpedcam);
     100    plist.AddToList(&pindiode);
     101    //    plist.AddToList(&pinhist);
    92102
    93103    //
    … …  
    109119    // Look at one specific pixel, after all the histogram manipulations:
    110120    //
    111     // hpedcam[17].DrawClone("fourierevents");
     121    hpedcam[9].DrawClone("fourier");
     122
    112123
    113124    MHCamera dispped0  (geomcam, "Ped;Pedestal",               "Mean per Slice");
    … …  
    215226    plist2.AddToList(&tlist2);
    216227
    217     MExtractedSignalCam   sigcam;
    218     MCalibrationCam       calcam;
     228    MExtractedSignalCam     sigcam;
     229    MCalibrationCam         calcam;
     230    MHCalibrationRelTimeCam timecam;
    219231    //
    220232    // Get the previously created MPedestalCam into the new Parameter List
    … …  
    224236    plist2.AddToList(&sigcam);
    225237    plist2.AddToList(&calcam);
    226 
     238    plist2.AddToList(&timecam);
    227239    //
    228240    // Get the MAGIC geometry
    … …  
    242254    //         
    243255    MExtractSignal2      sigcalc2;
    244     MArrivalTimeCalc     timecalc;
     256    MArrivalTimeCalc2    timecalc;
    245257    MCalibrationCalc     calcalc;
    246    
     258
     259    MFillH filltime("MHCalibrationRelTimeCam", "MArrivalTime");
    247260    //
    248261    // Set the range (other than default)
    … …  
    291304    //
    292305    MCalibrationBlindPix *bp = calcam.GetBlindPixel();
    293     //  bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson5);
     306    bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson5);
    294307
    295308    //
    … …  
    303316    tlist2.AddToList(&sigcalc2);
    304317    //
     318    // In case, you want to skip the cosmics rejection,
     319    // uncomment the next line
     320    //
     321    tlist2.AddToList(&cont);
     322    //
    305323    // In case, you want to skip the somewhat lengthy calculation
    306     // of the arrival times using a spline, uncomment the next line
    307     //
    308     // tlist2.AddToList(&timecalc);
    309     tlist2.AddToList(&cont);
     324    // of the arrival times using a spline, uncomment the next two lines
     325    //
     326    tlist2.AddToList(&timecalc);
     327    tlist2.AddToList(&filltime);
     328    //
    310329    tlist2.AddToList(&calcalc);
    311 
    312330    //
    313331    // Create and setup the eventloop
    … …  
    357375    MHCamera disp19 (geomcam, "Cal;Saturation",     "Pixels with saturated Hi Gain");
    358376
    359 
    360     MHCamera disp20 (geomcam, "cal;Ped",         "Pedestals");
    361     MHCamera disp21 (geomcam, "cal;PedRms",      "Pedestal RMS");
    362 
    363     MHCamera disp22 (geomcam, "cal;Time",        "Rel. Arrival Times");
    364     MHCamera disp23 (geomcam, "cal;SigmaTime",   "Sigma of Rel. Arrival Times");
    365     MHCamera disp24 (geomcam, "cal;TimeProb",    "Probability of Time Fit");
    366 
    367     MHCamera disp25 (geomcam, "cal;AbsTimeMean", "Abs. Arrival Times");
    368     MHCamera disp26 (geomcam, "cal;AbsTimeRms",  "RMS of Arrival Times");
    369 
     377    MHCamera disp20 (geomcam, "Cal;Ped",         "Pedestals");
     378    MHCamera disp21 (geomcam, "Cal;PedRms",      "Pedestal RMS");
     379
     380    MHCamera disp22 (geomcam, "time;Time",        "Rel. Arrival Times");
     381    MHCamera disp23 (geomcam, "time;SigmaTime",   "Sigma of Rel. Arrival Times");
     382    MHCamera disp24 (geomcam, "time;TimeProb",    "Probability of Time Fit");
     383    MHCamera disp25 (geomcam, "time;NotFitValid", "Pixels with not valid fit results");
     384    MHCamera disp26 (geomcam, "time;Oscillating", "Oscillating Pixels");
     385
     386    MHCamera disp27 (geomcam, "Cal;AbsTimeMean", "Abs. Arrival Times");
     387    MHCamera disp28 (geomcam, "Cal;AbsTimeRms",  "RMS of Arrival Times");
    370388
    371389    // Fitted charge means and sigmas
    … …  
    426444
    427445    // Relative Times
    428     disp20.SetCamContent(calcam,39);
    429     disp20.SetCamError(  calcam,40);
    430     disp21.SetCamContent(calcam,41);
    431     disp21.SetCamError(  calcam,42);
    432     disp22.SetCamContent(calcam,43);
     446    disp22.SetCamContent(timecam,0);
     447    disp22.SetCamError(  timecam,1);
     448    disp23.SetCamContent(timecam,2);
     449    disp23.SetCamError(  timecam,3);
     450    disp24.SetCamContent(timecam,4);
     451    disp25.SetCamContent(timecam,5);
     452    disp26.SetCamContent(timecam,6);
    433453
    434454    // Absolute Times
    435     disp23.SetCamContent(calcam,44);
    436     disp23.SetCamError(  calcam,45);
    437     disp24.SetCamContent(calcam,46);
    438     disp24.SetCamError(  calcam,47);
     455    disp27.SetCamContent(calcam,39);
     456    disp27.SetCamError(  calcam,40);
     457    disp28.SetCamContent(calcam,41);
     458    disp28.SetCamError(  calcam,42);
    439459
    440460
    … …  
    462482    disp17.SetYTitle("[1]");
    463483    disp18.SetYTitle("[1]");
    464 
    465     disp19.SetYTitle("Ped [FADC Counts ]");
    466     disp20.SetYTitle("RMS_{Ped} [FADC Counts ]");
    467 
    468     disp21.SetYTitle("Rel. Arr. Time [ns]");
    469     disp22.SetYTitle("\\sigma_{Time} [ns]");
    470     disp23.SetYTitle("P_{Time} [1]");
    471 
    472     disp24.SetYTitle("Mean Abs. Time [FADC slice]");
    473     disp25.SetYTitle("RMS Abs. Time [FADC slices]");
    474 
     484    disp19.SetYTitle("[1]");
     485
     486    disp20.SetYTitle("Ped [FADC Counts ]");
     487    disp21.SetYTitle("RMS_{Ped} [FADC Counts ]");
     488
     489    disp22.SetYTitle("Time Offset [ns]");
     490    disp23.SetYTitle("Timing resolution [ns]");
     491    disp24.SetYTitle("P_{Time} [1]");
     492
     493    disp25.SetYTitle("[1]");
     494    disp26.SetYTitle("[1]");
     495
     496    disp27.SetYTitle("Mean Abs. Time [FADC slice]");
     497    disp28.SetYTitle("RMS Abs. Time [FADC slices]");
    475498
    476499    gStyle->SetOptStat(1111);
    … …  
    503526    CamDraw(c4, disp6,calcam,1, 2 , 2);
    504527    CamDraw(c4, disp7,calcam,2, 2 , 2);
    505     //    CamDraw(c4, disp8,calcam,3, 3 , 2);
     528    CamDraw(c4, disp8,calcam,3, 3 , 2);
    506529
    507530    // Blind Pixel Method
    … …  
    541564    c9.Divide(2,3);
    542565
    543     CamDraw(c9,disp20,calcam,1,3,1);
    544     CamDraw(c9,disp21,calcam,2,3,2);
     566    CamDraw(c9,disp20,calcam,1,2,1);
     567    CamDraw(c9,disp21,calcam,2,2,2);
    545568
    546569
    … …  
    553576    CamDraw(c10,disp24,calcam,3,3,4);
    554577
     578    // Time Defects
     579    TCanvas &c11 = display->AddTab("Time Def.");
     580    c11.Divide(2,2);
     581
     582    CamDraw(c11, disp25,calcam,1,2, 0);
     583    CamDraw(c11, disp26,calcam,2,2, 0);
    555584
    556585    // Abs. Times
    557     TCanvas &c11 = display->AddTab("Abs. Times");
    558     c11.Divide(2,3);
    559 
    560     CamDraw(c11,disp25,calcam,1,2,2);
    561     CamDraw(c11,disp26,calcam,2,2,2);
     586    TCanvas &c12 = display->AddTab("Abs. Times");
     587    c12.Divide(2,3);
     588
     589    CamDraw(c12,disp27,calcam,1,2,2);
     590    CamDraw(c12,disp28,calcam,2,2,2);
    562591
    563592}
Note: See TracChangeset for help on using the changeset viewer.