Index: unk/MagicSoft/Mars/macros/tutorials/sumcurrents.C =================================================================== --- /trunk/MagicSoft/Mars/macros/tutorials/sumcurrents.C (revision 7156) +++ (revision ) @@ -1,162 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz, 6/2003 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -/////////////////////////////////////////////////////////////////////////// -// -// sumcurrents.C -// ============= -// -// This is a demonstration macro to display mean DC currents for all pixels. -// The input is cc report file. The output are histograms and plots. -// Using the MDirIter functionality you can process more than one file -// in one or more directories. For more information see MDirIter. -// -/////////////////////////////////////////////////////////////////////////// - -void ProcessFile(TString fname) -{ - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MParList plist; - - MTaskList tlist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - - // Create the magic geometry - MGeomCamMagic geom; - plist.AddToList(&geom); - - // First Task: Read file with image parameters - // (created with the star.C macro) - MReportFileRead read(fname); - read.SetHasNoHeader(); - read.AddToList("MReportCurrents"); - tlist.AddToList(&read); - - // create a task to fill a histogram - MFillH fill("MHCamEvent", "MCameraDC"); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - // - // Now display the result of the loop - // - MHCamEvent &h2 = *(MHCamEvent*)plist->FindObject("MHCamEvent"); - MHCamera &h = *(MHCamera*)h2.GetHistByName("sum"); -; - - TCanvas *c = MH::MakeDefCanvas(); - c->Divide(3, 2); - - MHCamera *disp1=h.Clone(); - MHCamera *disp2=h.Clone(); - MHCamera *disp3=h.Clone(); - disp2->SetCamContent(h, 1); - disp3->SetCamContent(h, 2); - - disp1->SetYTitle("I [nA]"); - disp2->SetYTitle("\\sigma_{I} [\\mu A]"); - disp3->SetYTitle("\\sigma_{I} [%]"); - disp1->SetName("Currents;avg"); - disp2->SetName("Currents;err"); - disp3->SetName("Currents;rel"); - disp1->SetTitle("Currents Average"); - disp2->SetTitle("Currents error"); - disp3->SetTitle("Currents relative error"); - - c->cd(1); - TText text(0.1, 0.95, &fname[fname.Last('/')+1]); - text.SetTextSize(0.03); - text.DrawClone(); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - gPad->SetLogy(); - disp1->Draw(); - disp1->SetBit(kCanDelete); - c->cd(2); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - gPad->SetLogy(); - disp2->Draw(); - disp2->SetBit(kCanDelete); - c->cd(3); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - gPad->SetLogy(); - disp3->Draw(); - disp3->SetBit(kCanDelete); - c->cd(4); - gPad->SetBorderMode(0); - disp1->Draw("EPhist"); - c->cd(5); - gPad->SetBorderMode(0); - gPad->SetLogy(); - disp2->Draw("Phist"); - c->cd(6); - gPad->SetBorderMode(0); - gPad->SetLogy(); - disp3->Draw("Phist"); - - c->SaveAs(fname(0, fname.Last('.')+1) + "ps"); - c->SaveAs(fname(0, fname.Last('.')+1) + "root"); -} - -void sumcurrents(const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "dc_*.txt", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname); - } -} Index: unk/MagicSoft/Mars/macros/tutorials/sumevents.C =================================================================== --- /trunk/MagicSoft/Mars/macros/tutorials/sumevents.C (revision 7156) +++ (revision ) @@ -1,155 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz, 6/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -void ProcessFile(TString fname) -{ - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MParList plist; - - MTaskList tlist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - - // First Task: Read file with image parameters - // (created with the star.C macro) - - MReadMarsFile read("Events", fname); - read.DisableAutoScheme(); - - MGeomApply geomapl; - MCerPhotAnal2 ncalc; - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&ncalc); - - MFillH fill("MHCamEvent", "MCerPhotEvt"); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - TCanvas *c = MH::MakeDefCanvas(); - c->Divide(3, 2); - - MHCamEvent &h = *(MHCamEvent*)plist->FindObject("MHCamEvent"); - MHCamera *disp0 = h.GetHistByName(); - MHCamera *disp1 = new MHCamera(geom, "MCerPhotEvt;avg", "Cerenkov Photons Avarage"); - MHCamera *disp2 = new MHCamera(geom, "MCerPhotEvt;err", "Cerenkov Photons Error"); - MHCamera *disp3 = new MHCamera(geom, "MCerPhotEvt;rel", "Cerenkov Photons ERR/VAL"); - disp1->SetBit(kCanDelete); - disp2->SetBit(kCanDelete); - disp3->SetBit(kCanDelete); - - disp1->SetCamContent(*disp0, 0); - disp2->SetCamContent(*disp0, 1); - disp3->SetCamContent(*disp0, 2); - - disp1->SetYTitle("S [au]"); - disp2->SetYTitle("\\sigma_{S} [au]"); - disp3->SetYTitle("\\sigma_{S} [%]"); - - TText text(0.1, 0.5, &fname[fname.Last('/')+1]); - text.SetTextSize(0.015); - text.DrawClone(); - - c->cd(1); - gStyle->SetOptStat(1111); - disp1->Draw("hist"); - gPad->Update(); - - c->cd(2); - gPad->SetLogy(); - gStyle->SetOptStat(1101); - disp2->Draw("hist"); - gPad->Update(); - - c->cd(3); - gPad->SetLogy(); - gStyle->SetOptStat(1101); - disp3->Draw("hist"); - gPad->Update(); - - c->cd(4); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp1->Draw(); - - c->cd(5); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp2->Draw(); - - c->cd(6); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp3->Draw(); - - c->SaveAs(fname(0, fname.Last('.')+1) + "ps"); - //c->SaveAs(fname(0, fname.Last('.')+1) + "root"); -} - -// ------------------------------------------------------------------------- -// -// plot.C -// -// This macro shows how to fill and display a histogram using Mars -// -void sumevents(const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "*1947*.root", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname); - return; - } -} Index: unk/MagicSoft/Mars/macros/tutorials/sumeventserr.C =================================================================== --- /trunk/MagicSoft/Mars/macros/tutorials/sumeventserr.C (revision 7156) +++ (revision ) @@ -1,158 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz, 6/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -void ProcessFile(TString fname) -{ - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MParList plist; - - MTaskList tlist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - - // First Task: Read file with image parameters - // (created with the star.C macro) - - MReadMarsFile read("Events", fname); - read.DisableAutoScheme(); - - //MPedestalSum ncalc; - //MCerPhotCalc ncalc; - //MCerPhotAnal ncalc; - MGeomApply geomapl; - MCerPhotAnal2 ncalc; - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&ncalc); - - MFillH fill("MHCamEvent", "MCerPhotEvt"); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - TCanvas *c = MH::MakeDefCanvas(); - c->Divide(3, 2); - - MHCamEvent &h = *(MHCamEvent*)plist->FindObject("MHCamEvent"); - MHCamera *disp0 = h.GetHistByName(); - MHCamera *disp1 = new MHCamera(geom, "MCerPhotEvt;avg", "Cerenkov Photons Avarage"); - MHCamera *disp2 = new MHCamera(geom, "MCerPhotEvt;err", "Cerenkov Photons Error"); - MHCamera *disp3 = new MHCamera(geom, "MCerPhotEvt;rel", "Cerenkov Photons ERR/VAL"); - disp1->SetBit(kCanDelete); - disp2->SetBit(kCanDelete); - disp3->SetBit(kCanDelete); - - disp1->SetCamContent(*disp0, 0); - disp2->SetCamContent(*disp0, 1); - disp3->SetCamContent(*disp0, 2); - - disp1->SetYTitle("S_{err} [au]"); - disp2->SetYTitle("\\sigma_{S_{err}} [au]"); - disp3->SetYTitle("\\sigma_{S_{err}} [%]"); - - TText text(0.1, 0.5, &fname[fname.Last('/')+1]); - text.SetTextSize(0.015); - text.DrawClone(); - - c->cd(1); - gStyle->SetOptStat(1111); - disp1->Draw("hist"); - gPad->Update(); - - c->cd(2); - gPad->SetLogy(); - gStyle->SetOptStat(1101); - disp2->Draw("hist"); - gPad->Update(); - - c->cd(3); - gPad->SetLogy(); - gStyle->SetOptStat(1101); - disp3->Draw("hist"); - gPad->Update(); - - c->cd(4); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp1->Draw(); - - c->cd(5); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp2->Draw(); - - c->cd(6); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp3->Draw(); - - c->SaveAs(fname(0, fname.Last('.')+1) + "ps"); - //c->SaveAs(fname(0, fname.Last('.')+1) + "root"); -} - -// ------------------------------------------------------------------------- -// -// plot.C -// -// This macro shows how to fill and display a histogram using Mars -// -void sumeventserr(const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "*1947*.root", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname); - return; - } -} Index: unk/MagicSoft/Mars/macros/tutorials/sumeventsrms.C =================================================================== --- /trunk/MagicSoft/Mars/macros/tutorials/sumeventsrms.C (revision 7156) +++ (revision ) @@ -1,159 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz, 6/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -void ProcessFile(TString fname) -{ - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MParList plist; - - MTaskList tlist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - - // First Task: Read file with image parameters - // (created with the star.C macro) - - MReadMarsFile read("Events", fname); - read.DisableAutoScheme(); - - MGeomApply geomapl; - MPedestalSum ncalc; - //MCerPhotCalc ncalc; - //MCerPhotAnal ncalc; - //MCerPhotAnal2 ncalc; - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&ncalc); - - MHCamEvent hist; - hist.SetType(3); - plist.AddToList(&hist); - - MFillH fill("MHCamEvent", "MCerPhotEvt"); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - TCanvas *c = MH::MakeDefCanvas(); - c->Divide(3, 2); - - MHCamEvent &h = *(MHCamEvent*)plist->FindObject("MHCamEvent"); - MHCamera *disp0 = h.GetHistByName(); - MHCamera *disp1 = new MHCamera(geom, "MCerPhotEvt;avg", "Cerenkov Photons RMS Avarage"); - MHCamera *disp2 = new MHCamera(geom, "MCerPhotEvt;err", "Cerenkov Photons RMS Error"); - MHCamera *disp3 = new MHCamera(geom, "MCerPhotEvt;rel", "Cerenkov Photons RMS ERR/VAL"); - disp1->SetBit(kCanDelete); - disp2->SetBit(kCanDelete); - disp3->SetBit(kCanDelete); - - disp1->SetCamContent(*disp0, 0); - disp2->SetCamContent(*disp0, 1); - disp3->SetCamContent(*disp0, 2); - - disp2->SetStats(kFALSE); - disp3->SetStats(kFALSE); - - disp1->SetYTitle("a.u."); - disp2->SetYTitle("err"); - disp3->SetYTitle("rel.err [%]"); - - TText text(0.1, 0.5, &fname[fname.Last('/')+1]); - text.SetTextSize(0.015); - text.DrawClone(); - - c->cd(1); - disp1->Draw("hist"); - - c->cd(2); - gPad->SetLogy(); - disp2->Draw("hist"); - - c->cd(3); - gPad->SetLogy(); - disp3->Draw("hist"); - - c->cd(4); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp1->Draw(); - - c->cd(5); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp2->Draw(); - - c->cd(6); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp3->Draw(); - - c->SaveAs(fname(0, fname.Last('.')+1) + "ps"); - //c->SaveAs(fname(0, fname.Last('.')+1) + "root"); -} - -// ------------------------------------------------------------------------- -// -// plot.C -// -// This macro shows how to fill and display a histogram using Mars -// -void sumeventsrms(const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "*ped*.root", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname); - return; - } -} Index: unk/MagicSoft/Mars/macros/tutorials/sumpedestalrms.C =================================================================== --- /trunk/MagicSoft/Mars/macros/tutorials/sumpedestalrms.C (revision 7156) +++ (revision ) @@ -1,157 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz, 6/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -void ProcessFile(TString fname) -{ - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MParList plist; - - MTaskList tlist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - - // First Task: Read file with image parameters - // (created with the star.C macro) - - MReadMarsFile read("Events", fname); - read.DisableAutoScheme(); - - MGeomApply geomapl; - MCerPhotAnal2 ncalc; - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&ncalc); - - MHCamEvent hist; - hist.SetType(1); - plist.AddToList(&hist); - - MFillH fill("MHCamEvent", "MPedestalCam"); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - TCanvas *c = MH::MakeDefCanvas(); - c->Divide(3, 2); - - MHCamEvent &h = *(MHCamEvent*)plist->FindObject("MHCamEvent"); - MHCamera *disp0 = h.GetHistByName(); - MHCamera *disp1 = new MHCamera(geom, "MPedestalCam;avg", "Pedestal-RMS Avarage"); - MHCamera *disp2 = new MHCamera(geom, "MPedestalCam;err", "Pedestal-RMS Avarage Error"); - MHCamera *disp3 = new MHCamera(geom, "MPedestalCam;rel", "Pedestal-RMS Avarage ERR/VAL"); - disp1->SetBit(kCanDelete); - disp2->SetBit(kCanDelete); - disp3->SetBit(kCanDelete); - - disp1->SetCamContent(*disp0, 0); - disp2->SetCamContent(*disp0, 1); - disp3->SetCamContent(*disp0, 2); - - disp1->SetYTitle("P_{rms} [au]"); - disp2->SetYTitle("\\sigma_{P_{rms}} [au]"); - disp3->SetYTitle("\\sigma_{P_{rms}} [%]"); - - TText text(0.1, 0.5, &fname[fname.Last('/')+1]); - text.SetTextSize(0.015); - text.DrawClone(); - - c->cd(1); - gStyle->SetOptStat(1111); - disp1->Draw("hist"); - gPad->Update(); - - c->cd(2); - gStyle->SetOptStat(1101); - disp2->Draw("hist"); - gPad->Update(); - - c->cd(3); - gStyle->SetOptStat(1101); - disp3->Draw("hist"); - gPad->Update(); - - c->cd(4); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp1->Draw(); - - c->cd(5); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp2->Draw(); - - c->cd(6); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp3->Draw(); - - c->SaveAs(fname(0, fname.Last('.')+1) + "ps"); - //c->SaveAs(fname(0, fname.Last('.')+1) + "root"); -} - -// ------------------------------------------------------------------------- -// -// plot.C -// -// This macro shows how to fill and display a histogram using Mars -// -void sumpedestalrms(const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "*1947*.root", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname); - return; - } -} Index: unk/MagicSoft/Mars/macros/tutorials/sumpedestals.C =================================================================== --- /trunk/MagicSoft/Mars/macros/tutorials/sumpedestals.C (revision 7156) +++ (revision ) @@ -1,153 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz, 6/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -void ProcessFile(TString fname) -{ - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MParList plist; - - MTaskList tlist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - - // First Task: Read file with image parameters - // (created with the star.C macro) - - MReadMarsFile read("Events", fname); - read.DisableAutoScheme(); - - MGeomApply geomapl; - MCerPhotAnal2 ncalc; - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&ncalc); - - MFillH fill("MHCamEvent", "MPedestalCam"); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - TCanvas *c = MH::MakeDefCanvas(); - c->Divide(3, 2); - - MHCamEvent &h = *(MHCamEvent*)plist->FindObject("MHCamEvent"); - MHCamera *disp0 = h.GetHistByName(); - MHCamera *disp1 = new MHCamera(geom, "MPedestalCam;avg", "Pedestals Avarage"); - MHCamera *disp2 = new MHCamera(geom, "MPedestalCam;rms", "Pedestals Avarage Error"); - MHCamera *disp3 = new MHCamera(geom, "MPedestalCam;rel", "Pedestals Avarage ERR/VAL"); - disp1->SetBit(kCanDelete); - disp2->SetBit(kCanDelete); - disp3->SetBit(kCanDelete); - - disp1->SetCamContent(*disp0, 0); - disp2->SetCamContent(*disp0, 1); - disp3->SetCamContent(*disp0, 2); - - disp1->SetYTitle("P [au]"); - disp2->SetYTitle("\\sigma_{P} [au]"); - disp3->SetYTitle("\\sigma_{P} [%]"); - - TText text(0.1, 0.5, &fname[fname.Last('/')+1]); - text.SetTextSize(0.015); - text.DrawClone(); - - c->cd(1); - gStyle->SetOptStat(1111); - disp1->Draw("hist"); - gPad->Update(); - - c->cd(2); - gStyle->SetOptStat(1101); - disp2->Draw("hist"); - gPad->Update(); - - c->cd(3); - gStyle->SetOptStat(1101); - disp3->Draw("hist"); - gPad->Update(); - - c->cd(4); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp1->Draw(); - - c->cd(5); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp2->Draw(); - - c->cd(6); - gPad->SetBorderMode(0); - gPad->Divide(1,1); - gPad->cd(1); - disp3->Draw(); - - c->SaveAs(fname(0, fname.Last('.')+1) + "ps"); - //c->SaveAs(fname(0, fname.Last('.')+1) + "root"); -} - -// ------------------------------------------------------------------------- -// -// plot.C -// -// This macro shows how to fill and display a histogram using Mars -// -void sumpedestals(const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "*1947*.root", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname); - return; - } -}