Ignore:
Timestamp:
12/09/03 18:25:21 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r2603 r2627  
    2323\* ======================================================================== */
    2424
    25 void calibration(TString pedname="../../Mars-0.8.2/20031102_02399_P_Unavailable_E.root",
    26                  TString calname="../../Mars-0.8.2/20031102_02400_D_Flip500Hz_E.root")
     25void calibration(TString pedname="/data/MAGIC/rootdata/2003_12_01/20031130_03340_P_CrabNebula_E.root",
     26                 TString calname="/data/MAGIC/rootdata/2003_12_01/20031130_03341_C_CrabNebula_E.root")
    2727{
    2828
    … …  
    6262    tlist.AddToList(&fill);
    6363
    64     //    MStatusDisplay *d1 = new MStatusDisplay;
     64    MStatusDisplay *d1 = new MStatusDisplay;
    6565 
    6666    // Set update time to 3s
    67     //    d1->SetUpdateTime(3000);
     67    d1->SetUpdateTime(3000);
    6868 
    6969    //
    … …  
    7272    MEvtLoop evtloop;
    7373    evtloop.SetParList(&plist);
    74     //    evtloop.SetDisplay(d1);       
     74    evtloop.SetDisplay(d1);     
    7575
    7676    //
    … …  
    114114    // we have to set the color of the pulser LED by hand
    115115    //
    116     calcalc.SetPulserColor(MCalibrationCalc::kEBlue);
     116    calcalc.SetPulserColor(MCalibrationCalc::kECT1);
    117117
    118118    tlist2.AddToList(&read2);
    … …  
    145145    cam.Print();
    146146
    147     Int_t pixnr;
    148 
    149     while (1)
    150       {
    151        
    152         cout << "Which pixel number do you want to display? (Press 0 to exit)" << endl;
    153         cin >> pixnr;
    154         if (pixnr == 0)
    155           break;
    156        
    157         if (pixnr >= 577)
    158           break;
    159        
    160         MCalibrationPix *pix = cam->GetCalibrationPix(pixnr);
    161         pix->Draw();
    162    
    163       }
    164    
    165147    //
    166148    // Here we are confronted to a serious bug in ROOT:
    … …  
    174156    MHCamera &disp0  = *h.GetHistByName();
    175157    MHCamera disp1   (geomcam, "MCalibrationCam;q", "Fitted Mean Charges");
    176     MHCamera disp2   (geomcam, "MCalibrationCam;errq", "Error of Fitted Mean Charges");
     158    //    MHCamera disp2   (geomcam, "MCalibrationCam;errq", "Error of Fitted Mean Charges");
    177159    MHCamera disp3   (geomcam, "MCalibrationCam;sigmaq", "Sigma of Fitted Mean Charges");
    178     MHCamera disp4   (geomcam, "MCalibrationCam;errsigmaq", "Error of Sigma of Fitted Mean Charges");
     160    //    MHCamera disp4   (geomcam, "MCalibrationCam;errsigmaq", "Error of Sigma of Fitted Mean Charges");
    179161    MHCamera disp5   (geomcam, "MCalibrationCam;probq", "Probability of Fit");
    180162    MHCamera disp6   (geomcam, "MCalibrationCam;t", "Arrival Times");
    … …  
    183165    MHCamera disp9   (geomcam, "MCalibrationCam;ped", "Pedestals");
    184166    MHCamera disp10  (geomcam, "MCalibrationCam;pedrms", "Pedestal RMS");
    185     MHCamera disp11  (geomcam, "MCalibrationCam;rq", "Reduced Charges");
    186     MHCamera disp12  (geomcam, "MCalibrationCam;rsigma", "Reduced Sigmas");
    187     MHCamera disp13  (geomcam, "MCalibrationCam;phe", "Nr. of Phe's (F-Factor Method)");
    188     MHCamera disp14  (geomcam, "MCalibrationCam;convphe", "Conversion Factor (F-Factor Method)");
    189     MHCamera disp15  (geomcam, "MCalibrationCam;photons", "Nr. of Photons (Blind Pixel Method)");
    190     MHCamera disp16  (geomcam, "MCalibrationCam;convphot", "Conversion Factor (Blind Pixel Method)");
     167    MHCamera disp11  (geomcam, "MCalibrationCam;rsigma", "Reduced Sigmas");
     168    MHCamera disp12  (geomcam, "MCalibrationCam;phe", "Nr. of Phe's (F-Factor Method)");
     169    MHCamera disp13  (geomcam, "MCalibrationCam;convphe", "Conversion Factor (F-Factor Method)");
     170    MHCamera disp14  (geomcam, "MCalibrationCam;photons", "Nr. of Photons (Blind Pixel Method)");
     171    MHCamera disp15  (geomcam, "MCalibrationCam;convphot", "Conversion Factor (Blind Pixel Method)");
    191172
    192173    disp1.SetCamContent(*cam, 0);
    193     disp2.SetCamContent(*cam, 1);
     174    disp1.SetCamError(*cam,1);
     175    //    disp2.SetCamContent(*cam, 1);
     176
    194177    disp3.SetCamContent(*cam, 2);
    195     disp4.SetCamContent(*cam, 3);
     178    disp3.SetCamError(*cam,3);
     179    //    disp4.SetCamContent(*cam, 3);
    196180    disp5.SetCamContent(*cam, 4);
     181
    197182    disp6.SetCamContent(*cam, 5);
     183    disp6.SetCamError(*cam,6);
    198184    disp7.SetCamContent(*cam, 6);
    199185    disp8.SetCamContent(*cam, 7);
     186
    200187    disp9.SetCamContent(*cam, 8);
     188    disp9.SetCamError(*cam,9);
    201189    disp10.SetCamContent(*cam, 9);
     190
    202191    disp11.SetCamContent(*cam, 10);
     192
    203193    disp12.SetCamContent(*cam, 11);
    204194    disp13.SetCamContent(*cam, 12);
    205195    disp14.SetCamContent(*cam, 13);
    206196    disp15.SetCamContent(*cam, 14);
    207     disp16.SetCamContent(*cam, 15);
     197    //    disp16.SetCamError(*cam, 16);
    208198
    209199
    210200    disp1.SetYTitle("Q [FADC counts]");
    211     disp2.SetYTitle("\\Delta Q [FADC counts]");
     201    //    disp2.SetYTitle("\\Delta Q [FADC counts]");
    212202    disp3.SetYTitle("\\sigma_{Q} [FADC counts]");
    213     disp4.SetYTitle("\\Delta {\\sigma_{Q}} [FADC counts]");
     203    //    disp4.SetYTitle("\\Delta {\\sigma_{Q}} [FADC counts]");
    214204    disp5.SetYTitle("P [au]");
    215205    disp6.SetYTitle("T [FADC slices]");
    … …  
    218208    disp9.SetYTitle("P [FADC counts/ slice ]");
    219209    disp10.SetYTitle("RMS_{P} [FADC counts / slice ]");
    220     disp11.SetYTitle("Q [FADC counts]");
    221     disp12.SetYTitle("\\sigma^2_{Q} - RMS^2_{P} [FADC counts^2]");
    222     disp13.SetYTitle("Nr Phe's");
    223     disp14.SetYTitle("Conversion Factor [Phe/FADC count]");
    224     disp15.SetYTitle("Nr Photons");
    225     disp16.SetYTitle("Conversion Factor [Ph/FADC count]");
     210    disp11.SetYTitle("\\sigma^2_{Q} - RMS^2_{P} [FADC counts^2]");
     211    disp12.SetYTitle("Nr Phe's");
     212    disp13.SetYTitle("Conversion Factor [Phe/FADC count]");
     213    disp14.SetYTitle("Nr Photons");
     214    disp15.SetYTitle("Conversion Factor [Ph/FADC count]");
    226215
    227216    MStatusDisplay *d2 = new MStatusDisplay;
    … …  
    230219    d2->SetUpdateTime(1000);
    231220
    232     TCanvas *c1 = &d2->AddTab("Charges Mean");
     221    TCanvas *c1 = &d2->AddTab("Fitted Charges");
    233222    c1->Divide(2, 2);
    234223
    … …  
    238227    gStyle->SetOptStat(1111);
    239228    obj=disp1.DrawCopy("hist");
     229    ((MHCamera*)obj)->AddNotify(*cam);
    240230
    241231    c1->cd(3);
    … …  
    245235    c1->cd(2);
    246236    gStyle->SetOptStat(1101);
    247     obj=disp2.DrawCopy("hist");
     237    obj=disp3.DrawCopy("hist");
     238    ((MHCamera*)obj)->AddNotify(*cam);
    248239
    249240    c1->cd(4);
    … …  
    251242    obj->Draw();
    252243
    253     TCanvas *c11 = &d2->AddTab("Charges Sigma");
    254     c11->Divide(2, 2);
    255 
    256     c11->cd(1);
    257     gStyle->SetOptStat(1101);
    258     obj=disp3.DrawCopy("hist");
    259 
    260     c11->cd(3);
    261     gPad->SetBorderMode(0);
    262     obj->Draw();
    263 
    264     c11->cd(2);
    265     gStyle->SetOptStat(1101);
    266     obj=disp4.DrawCopy("hist");
    267 
    268     c11->cd(4);
    269     gPad->SetBorderMode(0);
    270     obj->Draw();
    271 
    272244
    273245    TCanvas *c12 = &d2->AddTab("Fit Prob.");
    … …  
    277249    gStyle->SetOptStat(1101);
    278250    obj=disp5.DrawCopy("hist");
     251    ((MHCamera*)obj)->AddNotify(*cam);
    279252
    280253    c12->cd(2);
    … …  
    288261    gStyle->SetOptStat(1111);
    289262    obj=disp6.DrawCopy("hist");
     263    ((MHCamera*)obj)->AddNotify(*cam);
    290264
    291265    c2->cd(4);
    … …  
    295269    gStyle->SetOptStat(1101);
    296270    obj=disp7.DrawCopy("hist");
     271    ((MHCamera*)obj)->AddNotify(*cam);
    297272
    298273    c2->cd(5);
    … …  
    302277    gStyle->SetOptStat(1101);
    303278    obj=disp8.DrawCopy("hist");
     279    ((MHCamera*)obj)->AddNotify(*cam);
    304280
    305281    c2->cd(6);
    … …  
    312288    gStyle->SetOptStat(1111);
    313289    obj=disp9.DrawCopy("hist");
     290    ((MHCamera*)obj)->AddNotify(*cam);
    314291
    315292    c3->cd(3);
    … …  
    319296    gStyle->SetOptStat(1111);
    320297    obj=disp10.DrawCopy("hist");
     298    ((MHCamera*)obj)->AddNotify(*cam);
    321299
    322300    c3->cd(4);
    … …  
    324302
    325303    TCanvas *c4 = &d2->AddTab("Reduced Charges");
    326     c4->Divide(2, 2);
     304    c4->Divide(2,1);
    327305
    328306    c4->cd(1);
    329307    gStyle->SetOptStat(1111);
    330308    obj=disp11.DrawCopy("hist");
     309    ((MHCamera*)obj)->AddNotify(*cam);
    331310
    332311    c4->cd(3);
    333     obj->Draw();
    334 
    335     c4->cd(2);
    336     gStyle->SetOptStat(1101);
    337     obj=disp12.DrawCopy("hist");
    338 
    339     c4->cd(4);
    340312    obj->Draw();
    341313
    … …  
    345317    c5->cd(1);
    346318    gStyle->SetOptStat(1111);
     319    obj=disp12.DrawCopy("hist");
     320    ((MHCamera*)obj)->AddNotify(*cam);
     321
     322    c5->cd(3);
     323    obj->Draw();
     324
     325    c5->cd(2);
     326    gStyle->SetOptStat(1101);
    347327    obj=disp13.DrawCopy("hist");
    348 
    349     c5->cd(3);
    350     obj->Draw();
    351 
    352     c5->cd(2);
    353     gStyle->SetOptStat(1101);
    354     obj=disp14.DrawCopy("hist");
     328    ((MHCamera*)obj)->AddNotify(*cam);
    355329
    356330    c5->cd(4);
    … …  
    362336    c6->cd(1);
    363337    gStyle->SetOptStat(1111);
     338    obj=disp14.DrawCopy("hist");
     339    ((MHCamera*)obj)->AddNotify(*cam);
     340
     341    c6->cd(3);
     342    obj->Draw();
     343
     344    c6->cd(2);
     345    gStyle->SetOptStat(1101);
    364346    obj=disp15.DrawCopy("hist");
    365 
    366     c6->cd(3);
    367     obj->Draw();
    368 
    369     c6->cd(2);
    370     gStyle->SetOptStat(1101);
    371     obj=disp16.DrawCopy("hist");
     347    ((MHCamera*)obj)->AddNotify(*cam);
    372348
    373349    c6->cd(4);
    374350    obj->Draw();
    375351
    376 
    377 #endif
    378 
    379352}
    380353
Note: See TracChangeset for help on using the changeset viewer.