Ignore:
Timestamp:
11/26/03 02:21:16 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r2545 r2569  
    2424
    2525void calibration(TString pedname="./20031102_02399_P_Unavailable_E.root",
    26              TString calname="./20031102_02400_D_Flip500Hz_E.root")
     26                 TString calname="./20031102_02400_D_Flip500Hz_E.root")
    2727{
    2828
    … …  
    148148    MCalibrationCam *cam = plist2.FindObject("MCalibrationCam");
    149149    MCalibrationPix *pix = cam->GetCalibrationPix(523);
    150     pix->Draw();
    151 
    152     TCanvas *c1 = MH::MakeDefCanvas("Calibration1","Plots of FADC Charges",700,900);
    153     c1->Divide(5, 2);
     150    //pix->Draw();
     151
     152    gROOT->GetListOfCanvases()->Delete();
    154153
    155154    MHCamEvent &h    = *(MHCamEvent*)plist2->FindObject("MHCamEvent");
    156     MHCamera *disp0  = h.GetHistByName();
    157     MHCamera *disp1  = new MHCamera(geomcam, "MCalibrationCam;q", "Fitted Mean Charges");
    158     MHCamera *disp2  = new MHCamera(geomcam, "MCalibrationCam;errq", "Error of Fitted Mean Charges");
    159     MHCamera *disp3  = new MHCamera(geomcam, "MCalibrationCam;sigmaq", "Sigma of Fitted Mean Charges");
    160     MHCamera *disp4  = new MHCamera(geomcam, "MCalibrationCam;errsigmaq", "Error of Sigma of Fitted Mean Charges");
    161     MHCamera *disp5  = new MHCamera(geomcam, "MCalibrationCam;probq", "Probability of Fit");
    162     MHCamera *disp6  = new MHCamera(geomcam, "MCalibrationCam;t", "Arrival Times");
    163     MHCamera *disp7  = new MHCamera(geomcam, "MCalibrationCam;sigmat", "Sigma of Arrival Times");
    164     MHCamera *disp8  = new MHCamera(geomcam, "MCalibrationCam;probt", "Probability of Time Fit");
    165     MHCamera *disp9  = new MHCamera(geomcam, "MCalibrationCam;ped", "Pedestals");
    166     MHCamera *disp10 = new MHCamera(geomcam, "MCalibrationCam;pedrms", "Pedestal RMS");
    167     MHCamera *disp11 = new MHCamera(geomcam, "MCalibrationCam;rq", "Reduced Charges");
    168     MHCamera *disp12 = new MHCamera(geomcam, "MCalibrationCam;errrq", "Error of Reduced Charges");
    169 
    170     disp1->SetBit(kCanDelete);
    171     disp2->SetBit(kCanDelete);
    172     disp3->SetBit(kCanDelete);
    173     disp4->SetBit(kCanDelete);
    174     disp5->SetBit(kCanDelete);
    175     disp6->SetBit(kCanDelete);
    176     disp7->SetBit(kCanDelete);
    177     disp8->SetBit(kCanDelete);
    178     disp9->SetBit(kCanDelete);
    179     disp10->SetBit(kCanDelete);
    180     disp11->SetBit(kCanDelete);
    181     disp12->SetBit(kCanDelete);
    182 
    183     disp1->SetCamContent(*cam, 0);
    184     disp2->SetCamContent(*cam, 1);
    185     disp3->SetCamContent(*cam, 2);
    186     disp4->SetCamContent(*cam, 3);
    187     disp5->SetCamContent(*cam, 4);
    188     disp6->SetCamContent(*cam, 5);
    189     disp7->SetCamContent(*cam, 6);
    190     disp8->SetCamContent(*cam, 7);
    191     disp9->SetCamContent(*cam, 8);
    192     disp10->SetCamContent(*cam, 9);
    193     disp11->SetCamContent(*cam, 10);
    194     disp12->SetCamContent(*cam, 11);
    195 
    196     disp1->SetYTitle("Q [FADC counts]");
    197     disp2->SetYTitle("\\Delta_{Q} [FADC counts]");
    198     disp3->SetYTitle("\\sigma_{Q} [FADC counts]");
    199     disp4->SetYTitle("\\Delta_{\\sigma_{Q}} [FADC counts]");
    200     disp5->SetYTitle("P [au]");
    201     disp6->SetYTitle("T [FADC slices]");
    202     disp7->SetYTitle("\\Delta_{T} [FADC slices]");
    203     disp8->SetYTitle("P [au]");
    204     disp9->SetYTitle("P [FADC counts/ slice ]");
    205     disp10->SetYTitle("RMS_{P} [FADC counts / slice ]");
    206     disp11->SetYTitle("Q [FADC counts]");
    207     disp12->SetYTitle("\\Delta_{Q} [FADC counts]");
     155    MHCamera &disp0  = *h.GetHistByName();
     156    MHCamera disp1   (geomcam, "MCalibrationCam;q", "Fitted Mean Charges");
     157    MHCamera disp2   (geomcam, "MCalibrationCam;errq", "Error of Fitted Mean Charges");
     158    MHCamera disp3   (geomcam, "MCalibrationCam;sigmaq", "Sigma of Fitted Mean Charges");
     159    MHCamera disp4   (geomcam, "MCalibrationCam;errsigmaq", "Error of Sigma of Fitted Mean Charges");
     160    MHCamera disp5   (geomcam, "MCalibrationCam;probq", "Probability of Fit");
     161    MHCamera disp6   (geomcam, "MCalibrationCam;t", "Arrival Times");
     162    MHCamera disp7   (geomcam, "MCalibrationCam;sigmat", "Sigma of Arrival Times");
     163    MHCamera disp8   (geomcam, "MCalibrationCam;probt", "Probability of Time Fit");
     164    MHCamera disp9   (geomcam, "MCalibrationCam;ped", "Pedestals");
     165    MHCamera disp10  (geomcam, "MCalibrationCam;pedrms", "Pedestal RMS");
     166    MHCamera disp11  (geomcam, "MCalibrationCam;rq", "Reduced Charges");
     167    MHCamera disp12  (geomcam, "MCalibrationCam;errrq", "Error of Reduced Charges");
     168
     169    disp1.SetCamContent(*cam, 0);
     170    disp2.SetCamContent(*cam, 1);
     171    disp3.SetCamContent(*cam, 2);
     172    disp4.SetCamContent(*cam, 3);
     173    disp5.SetCamContent(*cam, 4);
     174    disp6.SetCamContent(*cam, 5);
     175    disp7.SetCamContent(*cam, 6);
     176    disp8.SetCamContent(*cam, 7);
     177    disp9.SetCamContent(*cam, 8);
     178    disp10.SetCamContent(*cam, 9);
     179    disp11.SetCamContent(*cam, 10);
     180    disp12.SetCamContent(*cam, 11);
     181
     182    disp1.SetYTitle("Q [FADC counts]");
     183    disp2.SetYTitle("\\Delta_{Q} [FADC counts]");
     184    disp3.SetYTitle("\\sigma_{Q} [FADC counts]");
     185    disp4.SetYTitle("\\Delta_{\\sigma_{Q}} [FADC counts]");
     186    disp5.SetYTitle("P [au]");
     187    disp6.SetYTitle("T [FADC slices]");
     188    disp7.SetYTitle("\\Delta_{T} [FADC slices]");
     189    disp8.SetYTitle("P [au]");
     190    disp9.SetYTitle("P [FADC counts/ slice ]");
     191    disp10.SetYTitle("RMS_{P} [FADC counts / slice ]");
     192    disp11.SetYTitle("Q [FADC counts]");
     193    disp12.SetYTitle("\\Delta_{Q} [FADC counts]");
    208194
    209195//    TText text(0.1, 0.5, &fname[fname.Last('/')+1]);
    210196//    text.SetTextSize(0.015);
    211197//    text.DrawClone();
     198    TCanvas *c1 = MH::MakeDefCanvas("Calibration1","Plots of FADC Charges",700,900);
     199    c1->Divide(5, 2);
     200
     201    TObject *obj;
    212202
    213203    c1->cd(1);
    214204    gStyle->SetOptStat(1111);
    215     disp1->Draw("hist");
    216     gPad->Update();
     205    obj=disp1.DrawCopy("hist");
     206
     207    c1->cd(6);
     208    gPad->SetBorderMode(0);
     209    obj->Draw();
    217210
    218211    c1->cd(2);
    219212    gStyle->SetOptStat(1101);
    220     disp2->Draw("hist");
    221     gPad->Update();
     213    obj=disp2.DrawCopy("hist");
     214
     215    c1->cd(7);
     216    gPad->SetBorderMode(0);
     217    obj->Draw();
    222218
    223219    c1->cd(3);
    224220    gStyle->SetOptStat(1101);
    225     disp3->Draw("hist");
    226     gPad->Update();
     221    obj=disp3.DrawCopy("hist");
     222
     223    c1->cd(8);
     224    gPad->SetBorderMode(0);
     225    obj->Draw();
    227226
    228227    c1->cd(4);
    229228    gStyle->SetOptStat(1101);
    230     disp4->Draw("hist");
    231     gPad->Update();
     229    obj=disp4.DrawCopy("hist");
     230
     231    c1->cd(9);
     232    gPad->SetBorderMode(0);
     233    obj->Draw();
    232234
    233235    c1->cd(5);
    234236    gStyle->SetOptStat(1101);
    235     disp5->Draw("hist");
    236     gPad->Update();
    237 
    238     c1->cd(6);
    239     gPad->SetBorderMode(0);
    240     gPad->Divide(1,1);
    241     gPad->cd(1);
    242     disp1->Draw();
    243 
    244     c1->cd(7);
    245     gPad->SetBorderMode(0);
    246     gPad->Divide(1,1);
    247     gPad->cd(1);
    248     disp2->Draw();
    249 
    250     c1->cd(8);
    251     gPad->SetBorderMode(0);
    252     gPad->Divide(1,1);
    253     gPad->cd(1);
    254     disp3->Draw();
    255 
    256     c1->cd(9);
    257     gPad->SetBorderMode(0);
    258     gPad->Divide(1,1);
    259     gPad->cd(1);
    260     disp4->Draw();
     237    obj=disp5.DrawCopy("hist");
    261238
    262239    c1->cd(10);
    263240    gPad->SetBorderMode(0);
    264     gPad->Divide(1,1);
    265     gPad->cd(1);
    266     disp5->Draw();
     241    obj->Draw();
    267242
    268243    TCanvas *c2 = MH::MakeDefCanvas("Calibration2","Plots of Arrival Times",700,900);
    … …  
    271246    c2->cd(1);
    272247    gStyle->SetOptStat(1111);
    273     disp6->Draw("hist");
    274     gPad->Update();
     248    obj=disp6.DrawCopy("hist");
     249
     250    c2->cd(4);
     251    obj->Draw();
    275252
    276253    c2->cd(2);
    277254    gStyle->SetOptStat(1101);
    278     disp7->Draw("hist");
    279     gPad->Update();
     255    obj=disp7.DrawCopy("hist");
     256
     257    c2->cd(5);
     258    obj->Draw();
    280259
    281260    c2->cd(3);
    282261    gStyle->SetOptStat(1101);
    283     disp8->Draw("hist");
    284     gPad->Update();
    285 
    286     c2->cd(4);
    287     gPad->SetBorderMode(0);
    288     gPad->Divide(1,1);
    289     gPad->cd(1);
    290     disp6->Draw();
    291 
    292     c2->cd(5);
    293     gPad->SetBorderMode(0);
    294     gPad->Divide(1,1);
    295     gPad->cd(1);
    296     disp7->Draw();
     262    obj=disp8.DrawCopy("hist");
    297263
    298264    c2->cd(6);
    299     gPad->SetBorderMode(0);
    300     gPad->Divide(1,1);
    301     gPad->cd(1);
    302     disp8->Draw();
    303 
     265    obj->Draw();
    304266
    305267    TCanvas *c3 = MH::MakeDefCanvas("Calibration3","Plots of Pedestals",700,900);
    … …  
    308270    c3->cd(1);
    309271    gStyle->SetOptStat(1111);
    310     disp9->Draw("hist");
    311     gPad->Update();
     272    obj=disp9.DrawCopy("hist");
     273
     274    c3->cd(3);
     275    obj->Draw();
    312276
    313277    c3->cd(2);
    314278    gStyle->SetOptStat(1101);
    315     disp10->Draw("hist");
    316     gPad->Update();
    317 
    318     c3->cd(3);
    319     gPad->SetBorderMode(0);
    320     gPad->Divide(1,1);
    321     gPad->cd(1);
    322     disp9->Draw();
     279    obj=disp10.DrawCopy("hist");
    323280
    324281    c3->cd(4);
    325     gPad->SetBorderMode(0);
    326     gPad->Divide(1,1);
    327     gPad->cd(1);
    328     disp10->Draw();
     282    obj->Draw();
    329283
    330284    TCanvas *c4 = MH::MakeDefCanvas("Calibration4","Plots of Reduced Charges",700,900);
    331     c3->Divide(2, 2);
     285    c4->Divide(2, 2);
    332286
    333287    c4->cd(1);
    334288    gStyle->SetOptStat(1111);
    335     disp11->Draw("hist");
    336     gPad->Update();
     289    obj=disp11.DrawCopy("hist");
     290
     291    c4->cd(3);
     292    obj->Draw();
    337293
    338294    c4->cd(2);
    339295    gStyle->SetOptStat(1101);
    340     disp12->Draw("hist");
    341     gPad->Update();
    342 
    343     c4->cd(3);
    344     gPad->SetBorderMode(0);
    345     gPad->Divide(1,1);
    346     gPad->cd(1);
    347     disp11->Draw();
     296    obj=disp12.DrawCopy("hist");
    348297
    349298    c4->cd(4);
    350     gPad->SetBorderMode(0);
    351     gPad->Divide(1,1);
    352     gPad->cd(1);
    353     disp12->Draw();
    354 
    355 
     299    obj->Draw();
    356300
    357301//    c->SaveAs(fname(0, fname.Last('.')+1) + "ps");
Note: See TracChangeset for help on using the changeset viewer.