Index: trunk/MagicSoft/Mars/macros/sumcurrents.C =================================================================== --- trunk/MagicSoft/Mars/macros/sumcurrents.C (revision 7157) +++ (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: trunk/MagicSoft/Mars/macros/sumevents.C =================================================================== --- trunk/MagicSoft/Mars/macros/sumevents.C (revision 7157) +++ (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: trunk/MagicSoft/Mars/macros/sumeventserr.C =================================================================== --- trunk/MagicSoft/Mars/macros/sumeventserr.C (revision 7157) +++ (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: trunk/MagicSoft/Mars/macros/sumeventsrms.C =================================================================== --- trunk/MagicSoft/Mars/macros/sumeventsrms.C (revision 7157) +++ (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: trunk/MagicSoft/Mars/macros/sumpedestalrms.C =================================================================== --- trunk/MagicSoft/Mars/macros/sumpedestalrms.C (revision 7157) +++ (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: trunk/MagicSoft/Mars/macros/sumpedestals.C =================================================================== --- trunk/MagicSoft/Mars/macros/sumpedestals.C (revision 7157) +++ (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; - } -} Index: trunk/MagicSoft/Mars/macros/tutorials/derotatedc.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/derotatedc.C (revision 7157) +++ (revision ) @@ -1,104 +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, 4/2004 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// derotatedc.C -// ============ -// -// Derotate a MCamEvent and fill a histogram with derotated data -// (sky-plot). -// -// As an input you need a merpped root file conataining DC information -// from a camera report file. To be able to derotate you also need -// aproproitate time-stamps and the corresponding pointing information. -// -// All events are filled into a 2D histograms - derotated. -// -// The example shows the usage of MHCamEventRot. The output is the derotated -// sky-plot. -// -///////////////////////////////////////////////////////////////////////////// - -void derotatedc() -{ - // - // 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); - - // Define Observatory location (for derotation) - MObservatory obs; - plist.AddToList(&obs); - - // Set the camera geometry (for histogram size) - MGeomCamMagic cam; - plist.AddToList(&cam); - - // setup pointing position - MPointingPos ppos; - ppos.SetSkyPosition(MAstro::Hms2Hor(5, 34, 31.9), MAstro::Dms2Deg(22, 0, 52.0)); - plist.AddToList(&ppos); - - // Define which file to read - MReadTree read("Currents", "../dc.root"); - read.DisableAutoScheme(); - - // Derotated histogram to fill - MHCamEventRot hist; - - // Set name of time container corresponding to your data - hist.SetNameTime("MTimeCurrents"); - - // Setup fill task - MFillH fill(&hist, "MCameraDC"); - - // Set a draw option for your 2D histogram - //fill.SetDrawOption("colz"); - - // Setup tasklist - tlist.AddToList(&read); - tlist.AddToList(&fill); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - MStatusDisplay *d = new MStatusDisplay; - evtloop.SetDisplay(d); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); -} Index: trunk/MagicSoft/Mars/macros/tutorials/evtrate.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/evtrate.C (revision 7157) +++ (revision ) @@ -1,78 +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, 12/2003 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -/////////////////////////////////////////////////////////////////////////// -// -// evtrate.C -// ========= -// -// Example macro how to calulate the eventrate (per event) and display -// the result versus time. -// -// As an input you need a merpped raw-data file containing correct times. -// The output is the plot: Eventrate vs. time. -// -/////////////////////////////////////////////////////////////////////////// - -void evtrate() -{ - // Setup parameter- and tasklist - MParList plist; - MTaskList tlist; - plist.AddToList(&tlist); - - // Setup reading task - MReadMarsFile read("Events"); - read.DisableAutoScheme(); - read.AddFile("test-time.root"); - - // Setup event rate calculator - MEventRateCalc calc; - // Setup number of events to be averaged - calc.SetNumEvents(200); - - // Setup histogram to be filles with rate - MHVsTime rate("MEventRate.fRate"); - - // Setup task to fill the histogram - MFillH fill(&rate, "MTime"); - - // Setup tasklist - tlist.AddToList(&read); - tlist.AddToList(&calc); - tlist.AddToList(&fill); - - // Execute your eventloop - MEvtLoop loop; - loop.SetParList(&plist); - - if (!loop.Eventloop()) - return; - - // print some execution statistics - tlist.PrintStatistics(); - - // Draw result - rate.DrawClone(); -} Index: trunk/MagicSoft/Mars/macros/tutorials/extendcam.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/extendcam.C (revision 7157) +++ (revision ) @@ -1,235 +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, 12/2004 -! -! Copyright: MAGIC Software Development, 2004-2005 -! -! -\* ======================================================================== */ - -////////////////////////////////////////////////////////////////////////////// -// -// extendcam.C -// -// This macro demonstrates how to extend the magic camera and performes -// a fourier transformation on it -// -///////////////////////////////////////////////////////////////////////////////// - -Bool_t HandleInput() -{ - // This must be there to get accesss to the GUI while the macro - // is still running! - - TTimer timer("gSystem->ProcessEvents();", 50, kFALSE); - while (1) - { - // - // While reading the input process gui events asynchronously - // - timer.TurnOn(); - TString input = Getline("Type 'q' to exit, to go on: "); - timer.TurnOff(); - - if (input=="q\n") - return kFALSE; - - if (input=="\n") - return kTRUE; - }; - - return kFALSE; -} - -TCanvas *c; -MHCamera *d1=new MHCamera(MGeomCamMagic()); -MHCamera *d2=new MHCamera; -MHCamera *d3=new MHCamera; -MHCamera *d4=new MHCamera; -MHCamera *d5=new MHCamera; -MHCamera *d6=new MHCamera(MGeomCamMagicXT()); - -MCerPhotEvt *evt=0; -MMcEvt *mc; - -MHexagonFreqSpace *fFreq1; -MHexagonFreqSpace *fFreq2; - -Int_t PreProcess(MParList *plist) -{ - fFreq1 = (MHexagonFreqSpace*)plist->FindObject("MHexagonFreqSpace1"); - if (!fFreq1) - return kFALSE; - fFreq2 = (MHexagonFreqSpace*)plist->FindObject("MHexagonFreqSpace2"); - if (!fFreq2) - return kFALSE; - - //Für Real-Raum - evt = (MCerPhotEvt*)plist->FindObject("MCerPhotEvt"); - if (!evt) - { - cout << "Fehler: MCerPhotEvt" << endl; - return kFALSE; - } - - mc = (MMcEvt*)plist->FindObject("MMcEvt"); - if (!mc) - { - cout << "Fehler: MMcEvt" << endl; - return kFALSE; - } - - c = new TCanvas("Events", "Real Events", 900, 600); - c->SetBorderMode(0); - c->Divide(3,2); - - MGeomCam *geomfreq=fFreq1->NewGeomCam(); - geomfreq->SetName("FreqSpace"); - - c->cd(1); - d1->SetPrettyPalette(); - d1->Draw(); - c->cd(2); - d2->SetGeometry(MGeomCamMagicXT()); - d2->SetPrettyPalette(); - d2->Draw(); - c->cd(4); - d3->SetGeometry(*geomfreq); - d3->SetPrettyPalette(); - d3->Draw(); - c->cd(5); - d4->SetGeometry(MGeomCamMagicXT()); - d4->SetPrettyPalette(); - d4->Draw(); - c->cd(3); - d5->SetGeometry(MGeomCamMagicXT()); - d5->SetPrettyPalette(); - d5->Draw(); - c->cd(6); - //d6->SetGeometry(MGeomCamMagicXT()); - d6->SetPrettyPalette(); - d6->Draw(); - return kTRUE; -} - -MClone *cl1; -MClone *cl2; -MClone *cl3; -MClone *cl4; - -Int_t Process() -{ - d1->SetCamContent(*(MCerPhotEvt*)cl1->GetClone()); - d2->SetCamContent(*(MCerPhotEvt*)cl2->GetClone()); - d3->SetCamContent(*fFreq2); - d4->SetCamContent(*(MCerPhotEvt*)cl4->GetClone()); - d5->SetCamContent(*(MCerPhotEvt*)cl3->GetClone()); - d6->SetCamContent(*evt); - - d2->SetMaximum(-1111); - d4->SetMaximum(-1111); - d2->SetMinimum(0); - d4->SetMinimum(0); - d2->SetMaximum(TMath::Max(d2->GetMaximum(), d4->GetMaximum())); - d4->SetMaximum(d2->GetMaximum()); - - for (int i=1; i<7; i++) - { - c->GetPad(i)->GetPad(1)->Modified(); - c->GetPad(i)->GetPad(1)->Update(); - } - - c->Modified(); - c->Update(); - - - return HandleInput(); -} - -void extendcam() -{ - MParList plist; - MTaskList tlist; - - MGeomCamMagicXT geom("MGeomCamMagicXT"); - plist.AddToList(&geom); - - plist.AddToList(&tlist); - - MReadMarsFile read("Events"); - read.DisableAutoScheme(); - read.AddFile("/home/kolodziejski/Diplomarbeit/marsdata/mcdata/*.root"); - //read.AddFile("/home/hoehne/data/mcdata/gammas/calib/calib_gamma_zbin0.root"); - //read.AddFile("/home/hoehne/data/mcdata/hadrons/calib/calib_proton_zbin0.root"); - - - // Setup intercative task calling the functions defined above - MGeomApply apply; //necessary for evt - MGeomCamMagicEnhance enhance; - // enhance.SetNameCerPhotEvtOut("MCerPhotEvt2"); - - MClone clone1("MCerPhotEvt", "Clone1"); - MClone clone2("MCerPhotEvt", "Clone2"); - MClone clone3("MCerPhotEvt", "Clone3"); - MClone clone4("MCerPhotEvt", "Clone4"); - - cl1 = &clone1; - cl2 = &clone2; - cl3 = &clone3; - cl4 = &clone4; - - MHexagonalFTCalc ftcalc; - ftcalc.SetNameGeomCam("MGeomCamMagicXT"); - // ftcalc.SetNameCerPhotEvt("MCerPhotEvt2"); - - MTaskInteractive mytask; - mytask.SetPreProcess(PreProcess); - mytask.SetProcess(Process); - - MImgCleanStd clean1(26, 18); - clean1.SetMethod(MImgCleanStd::kAbsolute); - clean1.SetNameGeomCam("MGeomCamMagicXT"); - - MImgCleanStd clean2(26, 18); - clean2.SetMethod(MImgCleanStd::kAbsolute); - clean2.SetNameGeomCam("MGeomCamMagicXT"); - - // Setup your tasklist - tlist.AddToList(&read); - tlist.AddToList(&apply); - tlist.AddToList(&clone1); - tlist.AddToList(&enhance); - tlist.AddToList(&clone2); - tlist.AddToList(&clean1); - tlist.AddToList(&clone3); - tlist.AddToList(&ftcalc); - tlist.AddToList(&clone4); - tlist.AddToList(&clean2); - tlist.AddToList(&mytask); - - // Run your analysis - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - if (!evtloop.Eventloop()) - return; - - // Print statistics information about your loop - tlist.PrintStatistics(); -} - Index: trunk/MagicSoft/Mars/macros/tutorials/hft.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/hft.C (revision 7157) +++ (revision ) @@ -1,236 +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 12/2004 -! -! Copyright: MAGIC Software Development, 2004-2005 -! -! -\* ======================================================================== */ - -/////////////////////////////////////////////////////////////////////////// -// -// hft.C -// -// example of using MHexagonal* -// -/////////////////////////////////////////////////////////////////////////// - -Bool_t HandleInput() -{ - // This must be there to get accesss to the GUI while the macro - // is still running! - - TTimer timer("gSystem->ProcessEvents();", 50, kFALSE); - while (1) - { - // - // While reading the input process gui events asynchronously - // - timer.TurnOn(); - TString input = Getline("Type 'q' to exit, to go on: "); - timer.TurnOff(); - - if (input=="q\n") - return kFALSE; - - if (input=="\n") - return kTRUE; - }; - - return kFALSE; -} - -// -// Setup the data-members of your 'virtual' class -// -MHCamera display[5]; - -TCanvas *c; -MParList *fParList; -MTaskList *fTaskList; - -MHexagonalFTCalc hft; -MHexagonFreqSpace *fFreq1; -MHexagonFreqSpace *fFreq2; - -MGeomCam *geomfreq; - -TH1F histo("", "", 35, -0.5, 35); - -// -// Called like all PreProcess functions of tasks. Get the access to -// the containers necessary for you. -// -Int_t PreProcess(MParList *plist) -{ - fFreq1 = (MHexagonFreqSpace*)plist->FindObject("MHexagonFreqSpace1"); - if (!fFreq1) - return kFALSE; - fFreq2 = (MHexagonFreqSpace*)plist->FindObject("MHexagonFreqSpace2"); - if (!fFreq2) - return kFALSE; - - geomfreq=fFreq1->NewGeomCam(); - geomfreq->SetName("FreqSpace"); - plist->AddToList(geomfreq); - - // return kTRUE; - - // Get parameter and tasklist, see Process - fParList = plist; - fTaskList = (MTaskList*)plist->FindObject("MTaskList"); - - // Get camera geoemtry - MGeomCam *geomcam = (MGeomCam*)plist->FindObject("MGeomCam"); - - // setup canvas and camera-histograms - c = new TCanvas("Events", "Real Events", 900, 600); - c->SetBorderMode(0); - c->Divide(3,2); - for (int i=0; i<3; i++) - { - display[i].SetGeometry(*geomfreq); - display[i].SetPrettyPalette(); - c->cd(i+1); - display[i].Draw(); - } - - display[1].SetName("In"); - display[2].SetName("Out"); - display[0].SetName("Freq"); - - return kTRUE; -} - -// -// Called like all Process functions of tasks. Process a single -// event - here display it. -// -Int_t Process() -{ - // For simplicity we search in the Process function for the - // objects. This is deprectaed! Store the pointers to the objects - // as data member and get the pointers in PreProcess. - MClone *clone1 = (MClone*)fTaskList->FindObject("Clone1"); - MClone *clone2 = (MClone*)fTaskList->FindObject("Clone2"); - MGeomCam *geom = (MGeomCam*)fParList->FindObject("MGeomCam"); - - // Fill the data into your camera-histograms - display[1].SetCamContent(*fFreq1, 1); - display[2].SetCamContent(*fFreq1, 2); - display[0].SetCamContent(*fFreq2, 0); - - display[1].SetMaximum(-1111); - display[2].SetMaximum(-1111); - display[1].SetMinimum(0); - display[2].SetMinimum(0); - display[1].SetMaximum(TMath::Max(display[1].GetMaximum(), - display[2].GetMaximum())); - display[2].SetMaximum(display[1].GetMaximum()); - - MHexagonalFT *ft = &hft.GetHFT(); - - histo.Reset(); - for (int i=0; iGetNumKnots(); i++) - { - Double_t val=0; - fFreq2->GetPixelContent(val, i, *display[1].GetGeometry()); - histo.Fill(ft->GetRow(i), val/(i+1)); - } - - c->cd(4); - histo.Draw(); - - c->cd(5); - - TH1D *obj1 = (TH1D*)display[1].Projection("Proj1"); - obj1->SetLineColor(kBlue); - obj1->Draw(); - obj1 = (TH1D*)display[2].Projection("Proj2"); - obj1->Draw("same"); - - c->cd(6); - display[0].DrawProjection(); - - // Update the display - for (int i=1; i<=3; i++) - { - c->GetPad(i)->GetPad(1)->Modified(); - c->GetPad(i)->GetPad(1)->Update(); - } - - return HandleInput(); -} - -// -// Called like all PostProcess functions of tasks. Delete -// instanciated objects. -// -Int_t PostProcess() -{ - delete c; - delete geomfreq; -} - -void hft(const char *fname="/home/hoehne/data/mcdata/hadrons/calib/calib_proton_zbin0.root") -{ - // Setup parameter- and tasklist - MParList plist; - MTaskList tlist; - - plist.AddToList(&tlist); - - // setup reading task - MReadMarsFile read("Events", fname); - read.DisableAutoScheme(); - - // Clone MCerPhotEvt befor eimage cleaning - MClone clone1("MCerPhotEvt", "Clone1"); - - // Setup image cleaning - MImgCleanStd clean(24, 12); - clean.SetMethod(MImgCleanStd::kAbsolute); - - // Clone MCerPhotEvt befor eimage cleaning - MClone clone2("MCerPhotEvt", "Clone2"); - - // Setup intercative task calling the functions defined above - MTaskInteractive mytask; - - mytask.SetPreProcess(PreProcess); - mytask.SetProcess(Process); - - // Setup your tasklist - tlist.AddToList(&read); - tlist.AddToList(&clone1); - tlist.AddToList(&hft); - tlist.AddToList(&clone2); - tlist.AddToList(&clean); - tlist.AddToList(&mytask); - - // Run your analysis - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - if (!evtloop.Eventloop()) - return; - - // Print statistics information about your loop - tlist.PrintStatistics(); -} - Index: trunk/MagicSoft/Mars/macros/tutorials/pedestalvstime.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/pedestalvstime.C (revision 7157) +++ (revision ) @@ -1,138 +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, Int_t idx) -{ - // - // 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); - - MHPixVsTime hist1(idx, "Pedestal"); - MHPixVsTime hist2(idx, "PedestalRMS"); - hist2.SetType(1); - plist.AddToList(&hist1); - plist.AddToList(&hist2); - - MFillH fill1("Pedestal", "MPedestalCam"); - MFillH fill2("PedestalRMS", "MPedestalCam"); - tlist.AddToList(&fill1); - tlist.AddToList(&fill2); - - // - // 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(2,2); - - c->cd(1); - TGraph &g1 = hist1.GetGraph(); - TH1 *h1 = g1.GetHistogram(); - h1->SetXTitle("Time [au]"); - h1->SetYTitle("P [au]"); - g1.DrawClone("A*"); - - c->cd(2); - TH1F h1f("Pedestal", "Pedestals", 21, 46, 56); - h1f.SetXTitle("P"); - h1f.SetYTitle("Counts"); - for (int i=0;iFit("gaus"); - - c->cd(3); - TGraph &g2 = hist2.GetGraph(); - TH1 *h2 = g2.GetHistogram(); - h2->SetXTitle("Time [au]"); - h2->SetYTitle("P_{rms} [au]"); - g2.DrawClone("A*"); - - c->cd(4); - TH1F h2f("PedestalRMS", "Pedestals RMS", 26, 0, 3.5); - h2f.SetXTitle("P_{rms}"); - h2f.SetYTitle("Counts"); - for (int i=0;iFit("gaus"); - - 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 pedestalvstime(Int_t idx=0, const char *dirname=".") -{ - MDirIter Next; - Next.AddDirectory(dirname, "raw*.root", -1); - - TString fname; - while (1) - { - fname = Next(); - if (fname.IsNull()) - break; - - ProcessFile(fname, idx); - return; - } -} Index: trunk/MagicSoft/Mars/macros/tutorials/pedvsevent.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/pedvsevent.C (revision 7157) +++ (revision ) @@ -1,214 +1,0 @@ - -///////////////////////////////////////////////////////////////////////////// -// // -// Author(s): S.C. Commichau, Javier Rico, 12/2003 // -// // -// Macro to generate pedestal vs time (event) plot for a single pixel // -// // -///////////////////////////////////////////////////////////////////////////// - - -const Int_t default_pixel = 1 ; - -void pedvsevent(Int_t pixel = default_pixel, - TString pname = "/disc02/Data/rootdata/Miscellaneous/2003_11_29/20031128_03118_P_Park-camera-closed_E.root") -{ - - TH1F* Pedestal = new TH1F("Pedestal","Pedestals",100,0,20); - TH1F* PedestalRMS = new TH1F("PedestalRMS","Pedestal RMS",100,0,10); - - //Create an empty parameter list and an empty task list - //the tasklist is identified in the eventloop by ist name - MParList plist; - MTaskList tlist; - - //Creates a MPedestalPix object for each pixel, i.e. it is a - //storage container for all Pedestal information in the camera - MPedestalCam cam; - - plist.AddToList(&cam); - - //MHCamEvent hist; - //hist.SetType(1); - //plist.AddToList(&hist); - - plist.AddToList(&tlist); - - //Setup task and tasklist for the pedestals - MReadMarsFile read("Events", pname); - read.DisableAutoScheme(); - - //Apply the geometry to geometry dependant containers. - //MGeomApply changes the size of the arrays in the containers to a size - //matching the number of pixels, eg: MPedestalCam, MBlindPixels - //Default geometry is MGeomCamMagic - MGeomApply geomapl; - - //Task to calculate pedestals - MPedCalcPedRun ped; - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&ped); - - //ped.SetPixel(pixel); - //ped.Draw(); - - //Create and setup the 1st Eventloop - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - //Execute first analysis, pedestals... - if (!tlist.PreProcess(&plist)) - return; - - const Int_t nevents = read.GetEntries(); - - Float_t x[nevents], rms[nevents], mean[nevents]; - - - Int_t i = 0; - while (tlist.Process()) - { - mean[i] = cam[pixel].GetPedestal(); - rms[i] = cam[pixel].GetPedestalRms(); - x[i] = i; - i++; - } - - TGraphErrors* pedgraph = new TGraphErrors(nevents,x,mean,NULL,NULL); - TGraphErrors* rmsgraph = new TGraphErrors(nevents,x,rms,NULL,NULL); - - //plist.FindObject("MPedestalCam")->Print(); - - tlist.PostProcess(); - //if (!evtloop.Eventloop()) - // return; - - Float_t finalmean = cam[pixel].GetPedestal(); - Float_t finalrms = cam[pixel].GetPedestalRms(); - - TLine* pedline = new TLine(0,finalmean,nevents,finalmean); - TLine* rmsline = new TLine(0,finalrms,nevents,finalrms); - - //The draw area - - - - gROOT->Reset(); - gStyle->SetOptStat(0); - // Set statistics options, 1111111 - // ||||||| - // ||||||histogram name - // |||||number of entries - // ||||mean value - // |||RMS - // ||underflows - // |overflows - // sum of bins - - // Set gPad options - gStyle->SetFrameBorderMode(0); - gStyle->SetPalette(1); - // Delete Frames of subCanvas' does not work, hook it Mr. gPad! - gStyle->SetFrameBorderSize(0); - gStyle->SetCanvasColor(0); - gStyle->SetFrameFillColor(0); - gStyle->SetTitleFont(102); - gStyle->SetTitleFillColor(0); - gStyle->SetTitleBorderSize(0); - gStyle->SetStatColor(0); - gStyle->SetStatBorderSize(0); - - // Set Canvas options - TCanvas *MyC1 = new TCanvas("MyC","Pedestal vs Event", 0, 100, 900, 500); - MyC->SetBorderMode(0); // Delete the Canvas' border line - - MyC->cd(); - gPad->SetBorderMode(0); - - // TLine* pedline = new TLine(0,finalmean,nevents,finalmean); - //TLine* rmsline = new TLine(0,finalrms,nevents,finalrms); - - tlist.PrintStatistics(); - //plist.FindObject("MPedestalCam")->Print(); - Size_t pos = pname.Last('/')+10; - TString iRun = TString(pname(pos,5)); - - char str[64]; - - sprintf(str,"Run %s Pixel %d",iRun.Data(),pixel); - - pedgraph->SetMaximum(30); - pedgraph->SetMinimum(0); - pedgraph->SetMarkerStyle(24); - pedgraph->SetMarkerSize(.5); - pedgraph->GetXaxis()->SetTitleFont(102); - pedgraph->GetYaxis()->SetTitleFont(102); - pedgraph->GetXaxis()->SetLabelFont(102); - pedgraph->GetYaxis()->SetLabelFont(102); - pedgraph->SetTitle(str); -// pedgraph->SetTitleFont(102); - pedgraph->GetYaxis()->SetTitleFont(102); - pedgraph->GetXaxis()->SetTitle("Event"); - pedgraph->GetYaxis()->SetTitle("[FADC Counts]"); - pedgraph->GetXaxis()->SetLimits(0,nevents-1); - - rmsgraph->SetMarkerStyle(25); - rmsgraph->SetMarkerSize(.5); - rmsgraph->SetMarkerColor(8); - - pedline->SetLineColor(2); - rmsline->SetLineColor(4); - pedline->SetLineWidth(2); - rmsline->SetLineWidth(2); - - pedgraph->Draw("AP"); - rmsgraph->Draw("P"); - - pedline->Draw("same"); - rmsline->Draw("same"); - - TLegend* leg = new TLegend(.14,.68,.39,.88); - leg->SetHeader(""); - leg->AddEntry(pedgraph,"Event based Pedestal","P"); - leg->AddEntry(pedline,"Run based Pedestal","L"); - leg->AddEntry(rmsgraph,"Event based RMS","P"); - leg->AddEntry(rmsline,"Run based RMS","L"); - leg->SetFillColor(0); - leg->SetLineColor(1); - leg->SetBorderSize(1); - - - leg->Draw("same"); -} - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Index: trunk/MagicSoft/Mars/macros/tutorials/pixfirerate.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/pixfirerate.C (revision 7157) +++ (revision ) @@ -1,133 +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): Abelardo Moralejo -! Author(s): Thomas Bretz -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// pixfirerate.C -// ============= -// -// This macro can help to find "hot" pixels firing too often -// or pixels firing too rarely. It plots camera displays showing how many -// times the FADC integral of each pixel has been found to be above pedestal -// by 10, 20, 50, 100 or 200 ADC counts. As "pedestal" we now use simply -// the signal in the first ADC slice, which seemed reasonable from a first -// look at the available test data. -// -///////////////////////////////////////////////////////////////////////////// - -void pixfirerate(TString filename="rawfile.root") -{ - // - // Update frequency by default = 1Hz - // - MStatusDisplay *d = new MStatusDisplay; - - // Set update time to 5s - // d->SetUpdateTime(5000); - - // Disable online update - // d->SetUpdateTime(-1); - - d->SetLogStream(&gLog, kTRUE); // Disables output to stdout - gLog.SetOutputFile("status.log", kTRUE); // Enable output to file - //gLog.EnableOutputDevice(MLog::eStdout); // Enable output to stdout again - - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MTaskList tlist; - MParList plist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - MReadMarsFile read("Events"); - read.DisableAutoScheme(); - read.AddFile(filename); - tlist.AddToList(&read); - - MGeomApply geomapl; - tlist.AddToList(&geomapl); - - // A list of threshold which should be displayed. The last entry - // MUST be -1. - Double_t threshold[] = { 10, 20, 100, 200, -1 }; - - Int_t cnt = 0; - while (threshold[cnt]>0) cnt++; - - // Create the corresponding fill tasks and containers - for (int i=0; iSetBit(kCanDelete); - fillh->SetBit(kCanDelete); - plist.AddToList(trigmap); - tlist.AddToList(fillh); - } - - // create the eventloop - MEvtLoop evtloop; - evtloop.SetParList(&plist); - evtloop.SetDisplay(d); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - // - // Make sure the display hasn't been deleted by the user while the - // eventloop was running. - // - if ((d = evtloop.GetDisplay())) - { - // Save data in a postscriptfile (status.ps) - d->SaveAsPS(); - /* - * ----------- Write status to a root file ------------ - * - TFile file("status.root", "RECREATE"); - d->Write(); - file.Close(); - delete d; - */ - } - -} Index: trunk/MagicSoft/Mars/macros/tutorials/pixsatrate.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/pixsatrate.C (revision 7157) +++ (revision ) @@ -1,120 +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): Abelardo Moralejo -! Author(s): Thomas Bretz -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// pixfirerate.C -// ============= -// -// This macro can help to find "hot" pixels firing too often -// or pixels firing too rarely. It plots camera displays showing how many -// times the FADC integral of each pixel has been found to be above pedestal -// by 10, 20, 50, 100 or 200 ADC counts. As "pedestal" we now use simply -// the signal in the first ADC slice, which seemed reasonable from a first -// look at the available test data. -// -///////////////////////////////////////////////////////////////////////////// - -void pixsatrate(TString filename="rawfile.root") -{ - // - // Update frequency by default = 1Hz - // - MStatusDisplay *d = new MStatusDisplay; - - // Set update time to 5s - // d->SetUpdateTime(5000); - - // Disable online update - // d->SetUpdateTime(-1); - - d->SetLogStream(&gLog, kTRUE); // Disables output to stdout - gLog.SetOutputFile("status.log", kTRUE); // Enable output to file - //gLog.EnableOutputDevice(MLog::eStdout); // Enable output to stdout again - - // - // Create a empty Parameter List and an empty Task List - // The tasklist is identified in the eventloop by its name - // - MTaskList tlist; - MParList plist; - plist.AddToList(&tlist); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - MReadMarsFile read("Events"); - read.DisableAutoScheme(); - read.AddFile(filename); - tlist.AddToList(&read); - - MGeomApply geomapl; - tlist.AddToList(&geomapl); - - // Create histograms with saturation limits at 254 - MHTriggerLvl0 trighi(254, "SaturationHi", "Saturation Rate of Hi Gains"); - MHTriggerLvl0 triglo(254, "SaturationLo", "Saturation Rate of Lo Gains"); - trighi.SetType(1); - triglo.SetType(2); - - // craete fill tasks to fill the histogarms - MFillH fillhi(&trighi, "MRawEvtData"); - MFillH filllo(&triglo, "MRawEvtData"); - tlist.AddToList(&fillhi); - tlist.AddToList(&filllo); - - // create eventloop - MEvtLoop evtloop; - evtloop.SetParList(&plist); - evtloop.SetDisplay(d); - - // - // Execute your analysis - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - // - // Make sure the display hasn't been deleted by the user while the - // eventloop was running. - // - if ((d = evtloop.GetDisplay())) - { - // Save data in a postscriptfile (status.ps) - d->SaveAsPS(); - /* - * ----------- Write status to a root file ------------ - * - TFile file("status.root", "RECREATE"); - d->Write(); - file.Close(); - delete d; - */ - } - -} Index: trunk/MagicSoft/Mars/macros/tutorials/readpix.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/readpix.C (revision 7157) +++ (revision ) @@ -1,97 +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 12/2000 (tbretz@uni-sw.gwdg.de) -! -! Copyright: MAGIC Software Development, 2000-2001 -! -! -\* ======================================================================== */ - - -#include - -void readpix() -{ - // - // open the file - // - TFile input("delme.root", "READ"); - - // - // open the Run Header tree - // - TTree *runtree = (TTree*) input.Get("RunHeaders") ; - - if (!runtree) - return; - - cout << " Entries in Tree RunHeaders: " << runtree->GetEntries() << endl ; - - // - // create an instance of MRawRunHeader, enable the branch and - // read the header, print it - // - MRawRunHeader *runheader = new MRawRunHeader(); - runtree->GetBranch("MRawRunHeader")->SetAddress(&runheader); - runtree->GetEvent(0); - runheader->Print(); - - // - // open the Data Tree - // - TTree *evttree = (TTree*) input.Get("Data") ; - - // - // create the instances of the data to read in - // - MRawEvtData *evtdata = new MRawEvtData(); - - // - // enable the corresponding branches - // - evttree->GetBranch("MRawEvtData")->SetAddress(&evtdata); - - evttree->GetEvent(0); - - MRawEvtPixelIter pixel(evtdata); - - while (pixel.Next()) - { - - cout << "Pixel Nr: " << dec << pixel.GetPixelId() << " "; - cout << setfill('0') << endl; - - cout << "Hi Gains: "; - for (int i=0; iGetNumHiGainSamples(); i++) - cout << setw(3) << (UInt_t)pixel.GetHiGainFadcSamples()[i]; - cout << endl; - - cout << "Lo Gains: "; - if (pixel.IsLoGain()) - for (int i=0; iGetNumLoGainSamples(); i++) - cout << setw(3) << (UInt_t)pixel.GetLoGainFadcSamples()[i]; - cout << endl << endl; - - if (pixel.GetPixelId() == 11) - pixel.Draw(); - }; -} - - - - Index: trunk/MagicSoft/Mars/macros/tutorials/readraw.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/readraw.C (revision 7157) +++ (revision ) @@ -1,95 +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 12/2000 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - - -void readraw(const char *fname="/data/MAGIC/Period016/mcdata/spot_1cm/standard/gamma/Gamma_zbin9_90_7_1740to1749_w0.root") -{ - // - // open the file - // - TFile input(fname, "READ"); - - // - // open the Run Header tree - // - TTree *runtree = (TTree*) input.Get("RunHeaders") ; - - if (!runtree) - return; - - cout << " Entries in Tree RunHeaders: " << runtree->GetEntries() << endl ; - - // - // create an instance of MRawRunHeader, enable the branch and - // read the header, print it - // - MRawRunHeader *runheader = new MRawRunHeader(); - runtree->GetBranch("MRawRunHeader.")->SetAddress(&runheader); - runtree->GetEvent(0); - runheader->Print(); - - // - // open the Data Tree - // - TTree *evttree = (TTree*) input.Get("Events"); - - // - // create the instances of the data to read in - // - MRawEvtHeader *evtheader = 0; - MTime *evttime = 0; - MRawEvtData *evtdata = 0; - MRawCrateArray *evtcrate = 0; - - // - // enable the corresponding branches - // - evttree->GetBranch("MRawEvtHeader.")->SetAddress(&evtheader); - evttree->GetBranch("MRawEvtData.")->SetAddress(&evtdata); - - // Use this for real data only - //evttree->GetBranch("MTime.")->SetAddress(&evttime); - //evttree->GetBranch("MRawCrateArray.")->SetAddress(&evtcrate); - - // - // loop over all events and print the information - // - Int_t iEnt = evttree->GetEntries(); - cout << " Entries in Tree Data: " << iEnt << endl; - - for (Int_t i = 0; iGetEvent(i); - - evtheader->Print(); - evtdata->Print(); - - // Use this for real data only - //evttime->Print(); - //evtcrate->Print(); - } -} - - Index: trunk/MagicSoft/Mars/macros/tutorials/readrep.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/readrep.C (revision 7157) +++ (revision ) @@ -1,125 +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, 11/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -////////////////////////////////////////////////////////////////////////////// -// -// This macro demonstrates how to read a central control report file. -// (rootification, see merpp, too) -// -// In a second eventloop it gives an example on how to read such a root file. -// -////////////////////////////////////////////////////////////////////////////// - -void readrep(const char *fname="CC_2003_11_04_23_53_18.rep") -{ - // - // Read a report file and write containers into a root file - // - MParList plist; - MTaskList tlist; - plist.AddToList(&tlist); - - MReportFileRead read(fname); - tlist.AddToList(&read); - - read.AddToList("DAQ"); - read.AddToList("Drive"); - read.AddToList("Camera"); - read.AddToList("Trigger"); - - MWriteRootFile write("test.root"); - write.AddContainer("MReportCamera", "Camera"); - write.AddContainer("MTimeCamera", "Camera"); - write.AddContainer("MCameraAUX", "Camera"); - write.AddContainer("MCameraCalibration", "Camera"); - write.AddContainer("MCameraCooling", "Camera"); - write.AddContainer("MCameraHV", "Camera"); - write.AddContainer("MCameraLV", "Camera"); - write.AddContainer("MCameraLids", "Camera"); - write.AddContainer("MReportTrigger", "Trigger"); - write.AddContainer("MTimeTrigger", "Trigger"); - write.AddContainer("MReportDrive", "Drive"); - write.AddContainer("MTimeDrive", "Drive"); - tlist.AddToList(&write); - - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - // ---------------------------------------------------------------- - // - // Read a report file and write containers into a root file - // - MTaskList tlist2; - plist.Replace(&tlist2); - - // Create a tasklist to process the read events from the Camera tree - MTaskList list1("ProcessCamera"); - MPrint print1("MTimeCamera"); - list1.AddToList(&print1); - - // Create a tasklist to process the read events from the Drive tree - MTaskList list2("ProcessDrive"); - MPrint print2("MTimeDrive"); - list2.AddToList(&print2); - - // Tell the reader to read the trees Drive, Trigger and Camera - MReadReports read; - read.AddTree("Drive"); - read.AddTree("Trigger"); - read.AddTree("Camera"); - //read.AddTree("Events", "MTime."); // for later use! - - // Now (not earlier!) set the file to read! - read.AddFile("test.root"); - - // First read the events - tlist.AddToList(&read); - // Process the events from the Camera tree with the task list list1 - tlist.AddToList(&list1, "Camera"); - // Process the events from the Camera tree with the task list list2 - tlist.AddToList(&list2, "Drive"); - - // The task lists list1 and list2 are only executed (according to - // their stream id - the second argument of AddToList) if a - // corresponding event was read and MReadReports has set the stream - // id accoringly. MReadReports always sets the stream id to the name - // of the tree from which the event was read - - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute the eventloop which should print the time-stamps of the subsystem - // events from Drive and Camera in the correct order... - // - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); -} Index: trunk/MagicSoft/Mars/macros/tutorials/readrfl.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/readrfl.C (revision 7157) +++ (revision ) @@ -1,106 +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, 5/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -Bool_t HandleInput() -{ - TTimer timer("gSystem->ProcessEvents();", 50, kFALSE); - while (1) - { - // - // While reading the input process gui events asynchronously - // - timer.TurnOn(); - TString input = Getline("Type 'q' to exit, to go on: "); - timer.TurnOff(); - - if (input=="q\n") - return kFALSE; - - if (input=="\n") - return kTRUE; - }; - - return kFALSE; -} - -Bool_t JumpTo(MTaskList &tlist, MReadRflFile &read, int runno, int evtno) -{ - if (runno<0 || evtno<0) - return tlist.Process(); - - return read.SearchFor(runno, evtno); -} - -void readrfl(int runno=-1, int evtno=-1, const char *fname="Gamma_zbin0_90_7_50700to50704_w0.rfl") -{ - MParList plist; - - MGeomCamMagic geomcam; - MRflEvtData event; - MRflEvtHeader evthead; - MRflRunHeader runhead; - MTaskList tlist; - - plist.AddToList(&geomcam); - plist.AddToList(&event); - plist.AddToList(&evthead); - plist.AddToList(&runhead); - plist.AddToList(&tlist); - - MReadRflFile read(fname); - tlist.AddToList(&read); - - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - if (!evtloop.PreProcess()) - return; - - MHCamera display(geomcam); - display.Draw(); - gPad->cd(1); - event.Draw(); - - - - cout << "Runno: " << runno << " Eventno: " << evtno << endl; - - while (JumpTo(tlist, read, runno, evtno)) - { - runno = -1; - - cout << "Run #" << runhead.GetRunNumber() << " "; - cout << "Event #" << evthead.GetEvtNumber() << endl; - - display.SetCamContent(event); - - gPad->Modified(); - gPad->Update(); - - if (!HandleInput()) - break; - } - - evtloop.PostProcess(); -} Index: trunk/MagicSoft/Mars/macros/tutorials/runbook.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/runbook.C (revision 7157) +++ (revision ) @@ -1,150 +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 12/2000 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -/////////////////////////////////////////////////////////////////////////// -// -// runbook.C -// ========= -// -// This is an example macro how to extract the runbook information of -// a given date from the runbook files. -// -// The argument is a date in regexp format, eg. -// .x runbook.C("2004-05-1[0-9]") -// -// With an MDirIter you can tell the macro which directories and files -// should be searched. -// -// The output are the runbook entries matching your query. -// -/////////////////////////////////////////////////////////////////////////// - -void Wrap(TString &str) -{ - // Wrap lines to 80 chars - while (str.Length()>80) - { - TString ss = str(0, 80); - Int_t l = ss.Last(' '); - if (l<0) - break; - cout << str(0, l) << endl; - str.Remove(0, l); - } - cout << str << endl; -} - -void ProcessFile(TString fname, TString date) -{ - // File processed - cout << fname << endl; - - // Open file - ifstream fin(fname); - - TString line; - - int i=81; - - TRegexp r0("20[0-2][0-9]-[0-1][0-9]-[0-3][0-9] [0-2][0-9]:[0-5][0-9]:[0-5][0-9]"); - TRegexp r1(date+" [0-2][0-9]:[0-5][0-9]:[0-5][0-9]"); - - int lcnt=0; - - int print = 0; - - while (!fin.eof()) - { - // Read file line by line - char *txt=new char[i]; - fin.getline(txt, i-1); - - if (lcnt++<4) - continue; - - if (!fin) - txt[i-1]=0; - - line += txt; - delete txt; - - if (line.Length()==0) - { - if (print) - cout << endl; - continue; - } - - if (!fin) - { - i *= 2; - fin.clear(); - continue; - } - - // Check whether line matches regexp - if (!line(r0).IsNull()) - print = !line(r1).IsNull(); - - // Wrap lines to 80 chars - if (print) - Wrap(line); - line = ""; - i=81; - } - cout << "DONE." << endl; -} - -void runbook(const char *d="20[0-2][0-9]-[0-1][0-9]-[0-3][0-9]") -{ - // Regexp to check for valid date - TString regexp = "20[0-2][0-9]-[0-1][0-9]-[0-3][0-9]"; - - if (regexp!=d) - { - TRegexp r(regexp); - - TString date=d; - if (date(r).IsNull()) - { - cout << "ERROR - Date has wrong format (2003-05-06)" << endl; - return; - } - } - - // Tell which dierctories and files to search - MDirIter Next; - Next.AddDirectory("/home/MAGIC/online_data/runbook", "CC_*.rbk"); - - // Loop over files - TString name; - while (1) - { - name=Next(); - if (name.IsNull()) - break; - - ProcessFile(name, d); - } -} Index: trunk/MagicSoft/Mars/macros/tutorials/sectorvstime.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sectorvstime.C (revision 7157) +++ (revision ) @@ -1,100 +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 12/2000 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// sectorvstime.C -// ============== -// -// In this example we plot the mean content of the right and left half of -// camera. As an input we use a class derived from MCamEvent. Here this -// are dc currents read directly from a camera control report file. -// -// The output are two histograms one for each half. -// -///////////////////////////////////////////////////////////////////////////// - -void sectorvstime() -{ - // Initialize Mars environment - MParList plist; - MTaskList tlist; - plist.AddToList(&tlist); - - // Create Magic camera geometry - MGeomCamMagic cam; - plist.AddToList(&cam); - - // Which DC file to read? - MReportFileRead read("/data/MAGIC/Period013/cacodata/2004_01_26/dc_2004_01_26_05_35_10_12117_OffMrk421-1.txt"); - read.SetHasNoHeader(); - read.AddToList("MReportCurrents"); - - // Initialize histogram - MHSectorVsTime hist1; - hist1.SetNameTime("MTimeCurrents"); - - // Define sectors you want to display the mean from - TArrayI s0(3); - s0[0] = 6; - s0[1] = 1; - s0[2] = 2; - - // Define area index [0=inner, 1=outer] - TArrayI inner(1); - inner[0] = 0; - - // Don't call this if you want to have all sectors - hist1.SetSectors(s0); - - // Don't call this if you want to have all area indices - hist1.SetAreaIndex(inner); - - // Task to fill the histogram - MFillH fill1(&hist1, "MCameraDC"); - - // Also fill a histogram with the mean of all pixels - MHCamEvent hist2; - MFillH fill2(&hist2, "MCameraDC"); - - // Setup Tasklist - tlist.AddToList(&read); - tlist.AddToList(&fill1); - tlist.AddToList(&fill2); - - // Setup Eventloop - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // Run Eventloop - if (!evtloop.Eventloop()) - return; - - // Print some statistics - tlist.PrintStatistics(); - - // Draw clones of the histograms - hist1.DrawClone(); - hist2.DrawClone(); -} Index: trunk/MagicSoft/Mars/macros/tutorials/spline.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/spline.C (revision 7157) +++ (revision ) @@ -1,109 +1,0 @@ -/* This macro is defined as a class for debugging (with CInt) reasons - -To use it at the root prompt type: - -root [0] .L spline.C -root [1] TestSpline::sp() -*/ -/* Example of a spline. You can use it as Test. If you think there are some - bugs in the MCubicSpline class please mail to: raducci@fisica.uniud.it */ - -class TestSpline -{ -public: - void sp(); -}; - -void TestSpline::sp() -{ -gROOT -> Reset(); - -//Here are defined the points. X goes from 0 to 14 (as the fadc slices...) -//Y are arbitrary values - -/* User Change */ -const Byte_t y[]={0x0F,0x10,0x2F,0x7F,0xAA,0x6C,0x14,0x13,0x15,0x18,0x21,0x12,0x11,0x14,0x13}; -/* End user Change */ -const Byte_t x[]={0x00,0x01,0x02,0x03,0x04,0x05,0x06,0x07,0x08,0x09,0x0A,0x0B,0x0C,0x0D,0x0E}; - -/*This cast is needed only to show graphically the output. Don' t needed if you - use the spline to calc the arrival times */ - -Int_t *newX = new Int_t[15]; -Int_t *newY = new Int_t[15]; - -for (Int_t i = 0; i < 15; i++) -{ - newX[i] = (Int_t) x[i]; - newY[i] = (Int_t) y[i]; -} - -//Canvas to display output -TCanvas *c = new TCanvas ("c1","c1",800,600); - -//Graph containting only the points (displayed as stars) -TGraph *g1 = new TGraph(15,newX,newY); - -g1 -> Draw("A*"); - -//Spline constructor(specialized for 15 slices using Bytes as values. There exist another constructor. -MCubicSpline *s = new MCubicSpline(y); - -//*spline and *ab are two arrays containing some values evaluated from the spline -Double_t *spline = new Double_t[139]; -Double_t *ab = new Double_t[139]; -Double_t step = 0.0; - -for (Int_t i = 0; i < 139; i++) -{ - spline[i] = s->Eval(step); - ab[i] = step; - step += 0.1; -} - -//Graph of the sline. The points calculated are joined with a red line. If the stars lie -//on the red line, then the Spline class is working properly -TGraph *g2 = new TGraph(139,ab,spline); - -g2 -> SetLineColor(2); -g2 -> Draw("C"); - -//Maximum and minimum evaluation -Double_t *mm = new Double_t[2]; -Double_t *abmm = new Double_t[2]; - -mm[0] = s->EvalMin(); -mm[1] = s->EvalMax(); -abmm[0] = s->EvalAbMin(); -abmm[1] = s->EvalAbMax(); - -//Display the max and the min using two black squares. If they lie on the max and min -//of the red line, then the Spline class is working properly - -TGraph *g3 = new TGraph(2,abmm,mm); - -g3 -> SetMarkerStyle(21); -g3 -> Draw("P"); - -//Test of the Cardan formula. Find the point(abval) where the Spline value is val -Double_t val = 82.5; -Double_t abval = s->FindVal(val, abmm[1], 'l'); - -//Display this point. Again, if it lies on the red line, then we are right. -//It's a black triangle - -TGraph *g4 = new TGraph(1,&abval,&val); - -g4 -> SetMarkerStyle(22); -g4 -> Draw("P"); - -//Free memory -s->~MCubicSpline(); -delete [] newX; -delete [] newY; -delete [] spline; -delete [] ab; -delete [] mm; -delete [] abmm; - -} Index: trunk/MagicSoft/Mars/macros/tutorials/starplot.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/starplot.C (revision 7157) +++ (revision ) @@ -1,115 +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 et al, 08/2002 -! -! Copyright: MAGIC Software Development, 2000-2002 -! -! -\* ======================================================================== */ - - -void starplot(const char *filename="Gamma_*.root") -{ - // - // This is a demonstration program which plots the Hillas - // parameter from a file created with star.C - // - - // - // 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); - - // - // Use this if you want to change the binning of one of - // the histograms. You can use: - // BinningConc, BinningConc1, BinningAsym, BinningM3Long, - // BinningM3Trans, BinningWidth, BinningLength, BinningDist, - // BinningHeadTail, BinningAlpha, BinningSize, BinningDelta, - // BinningPixels and BinningCamera - // - // For more information see MBinning and the corresponding - // histograms - // - // MBinning binsalpha("BinningAlpha"); - // binsalpha.SetEdges(90, 0, 90); // 90 bins from 0 to 90 deg - // plist.AddToList(&binsalpha); - - // MBinning binssize("BinningSize"); - // binssize.SetEdgesLog(50, 1, 1e7); - // plist.AddToList(&binssize); - - // - // Now setup the tasks and tasklist: - // --------------------------------- - // - // The first argument is the tree you want to read. - // Events: Cosmic ray events - // PedEvents: Pedestal Events - // CalEvents: Calibration Events - // - MReadMarsFile read("Events", filename); - read.DisableAutoScheme(); - - MGeomApply geomapl; - - MFillH hfill1("MHHillas", "MHillas"); - MFillH hfill2("MHHillasExt"); - MFillH hfill3("MHStarMap", "MHillas"); - MFillH hfill4("HistExtSource [MHHillasExt]", "MHillasSrc"); - MFillH hfill5("HistSource [MHHillasSrc]", "MHillasSrc"); - - tlist.AddToList(&read); - tlist.AddToList(&geomapl); - tlist.AddToList(&hfill1); - tlist.AddToList(&hfill2); - tlist.AddToList(&hfill3); - tlist.AddToList(&hfill4); - tlist.AddToList(&hfill5); - - // - // Create and setup the eventloop - // - MEvtLoop evtloop; - evtloop.SetParList(&plist); - - // - // Execute your analysis - // - MProgressBar bar; - evtloop.SetProgressBar(&bar); - if (!evtloop.Eventloop()) - return; - - tlist.PrintStatistics(); - - // - // After the analysis is finished we can display the histograms - // - plist.FindObject("MHHillas")->DrawClone(); - plist.FindObject("MHHillasExt")->DrawClone(); - plist.FindObject("MHStarMap")->DrawClone(); - plist.FindObject("HistSource")->DrawClone(); - plist.FindObject("HistExtSource")->DrawClone(); -} - Index: trunk/MagicSoft/Mars/macros/tutorials/sumcurrents.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sumcurrents.C (revision 7158) +++ trunk/MagicSoft/Mars/macros/tutorials/sumcurrents.C (revision 7158) @@ -0,0 +1,162 @@ +/* ======================================================================== *\ +! +! * +! * 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: trunk/MagicSoft/Mars/macros/tutorials/sumevents.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sumevents.C (revision 7158) +++ trunk/MagicSoft/Mars/macros/tutorials/sumevents.C (revision 7158) @@ -0,0 +1,155 @@ +/* ======================================================================== *\ +! +! * +! * 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: trunk/MagicSoft/Mars/macros/tutorials/sumeventserr.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sumeventserr.C (revision 7158) +++ trunk/MagicSoft/Mars/macros/tutorials/sumeventserr.C (revision 7158) @@ -0,0 +1,158 @@ +/* ======================================================================== *\ +! +! * +! * 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: trunk/MagicSoft/Mars/macros/tutorials/sumeventsrms.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sumeventsrms.C (revision 7158) +++ trunk/MagicSoft/Mars/macros/tutorials/sumeventsrms.C (revision 7158) @@ -0,0 +1,159 @@ +/* ======================================================================== *\ +! +! * +! * 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: trunk/MagicSoft/Mars/macros/tutorials/sumpedestalrms.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sumpedestalrms.C (revision 7158) +++ trunk/MagicSoft/Mars/macros/tutorials/sumpedestalrms.C (revision 7158) @@ -0,0 +1,157 @@ +/* ======================================================================== *\ +! +! * +! * 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: trunk/MagicSoft/Mars/macros/tutorials/sumpedestals.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/sumpedestals.C (revision 7158) +++ trunk/MagicSoft/Mars/macros/tutorials/sumpedestals.C (revision 7158) @@ -0,0 +1,153 @@ +/* ======================================================================== *\ +! +! * +! * 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; + } +} Index: trunk/MagicSoft/Mars/macros/tutorials/testenv.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/testenv.C (revision 7157) +++ (revision ) @@ -1,120 +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, 12/2000 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -/////////////////////////////////////////////////////////////////////////// -// -// testenv.C -// ========= -// -// This example reads an config-file (steering card, input card, ...) -// The contents of the file a forwarded to the apropriate eventloop -// setup in your file. -// -// All tasks and containers in an eventloop should implement the -// ReadEnv/WriteEnv function (for an example see the tasks used below). -// -// The functions gets the corresponding setup data from the file as an -// argument and can change their behaviour and setup on this information. -// -/////////////////////////////////////////////////////////////////////////// - -void testenv() -{ - // Setup for all MHMatrix objects is done by: - // MHMatrix.Column0: ... - // - // This can be overwritten for a MHMatrix called MatrixGammas by: - // MatrixGammas.Column0; - // - // This can be overwritten for all MHMatrix in one Job by: - // Job1.MHMatrix.Column0; - // - // This can be overwritten for a MHMatrix called MatrixGammas in one Job by: - // Job1.MatrixGammas.Column0; - // - TEnv env(".marsrc"); - - // - // For developers: Set this to kTRUE to see how the TEnv file - // entries are checked. - // - Bool_t print = kFALSE; - - // ------------ Job 1 ------------------- - if (env.GetValue("Job1", kFALSE)) - { - cout << "++++++++++++++++++ Job 1 +++++++++++++++++++" << endl; - MParList plist1; - MTaskList tlist1; - - plist1.AddToList(&tlist1); - - MReadMarsFile read1("Events"); - - MHMatrix matrix1("MatrixGammas"); - - MFillH fillm1(&matrix1); - plist1.AddToList(&matrix1); - - tlist1.AddToList(&read1); - tlist1.AddToList(&fillm1); - - MEvtLoop evtloop1("Job1"); - evtloop1.SetParList(&plist1); - cout << "--------------------------------------------" << endl; - evtloop1.ReadEnv(env, "", print); - cout << "--------------------------------------------" << endl; - evtloop1.Eventloop(); - cout << endl; - } - - // ------------ Job 2 ------------------- - - if (env.GetValue("Job2", kTRUE)) - { - cout << "++++++++++++++++++ Job 2 +++++++++++++++++++" << endl; - MParList plist2; - MTaskList tlist2; - - plist2.AddToList(&tlist2); - - MReadMarsFile read2("Events"); - - MHMatrix matrix2("MatrixGammas"); - - MFillH fillm2(&matrix2); - plist2.AddToList(&matrix2); - - tlist2.AddToList(&read2); - tlist2.AddToList(&fillm2); - - MEvtLoop evtloop2("Job2"); - evtloop2.SetParList(&plist2); - cout << "--------------------------------------------" << endl; - evtloop2.ReadEnv(env, "", print); - cout << "--------------------------------------------" << endl; - evtloop2.Eventloop(); - cout << endl; - } -} Index: trunk/MagicSoft/Mars/macros/tutorials/weights.C =================================================================== --- trunk/MagicSoft/Mars/macros/tutorials/weights.C (revision 7157) +++ (revision ) @@ -1,138 +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): Marcos Lopez, 10/2003 -! -! Copyright: MAGIC Software Development, 2000-2004 -! -! -\* ======================================================================== */ - -////////////////////////////////////////////////////////////////////////////// -// -// weights.C -// ========= -// -// This macro shows how to use the class MMcWeightEnergySpecCalc -// to convert the energy spectrum of the MC showers generated with Corsika, -// to a different one. -// -////////////////////////////////////////////////////////////////////////////// - -Double_t myfunction(Double_t *x, Double_t *par) -{ - Double_t xx = x[0]; - - return pow(xx,-2)*exp(-xx/20); -} - -void weights(TString filename="/up1/data/Magic-MC/CameraAll/Gammas/zbin0/Gamma_zbin0_0_7_1000to1009_w0-4:4:2.root") -{ - - // - // PartList - // - MParList parlist; - MTaskList tasklist; - - MHMcEnergyImpact h1("h1"); - MHMcEnergyImpact h2("h2"); - parlist.AddToList(&h1); - parlist.AddToList(&h2); - - MBinning binsenergy("BinningEnergy"); - binsenergy.SetEdgesLog(100, 1, 1e5); - parlist.AddToList(&binsenergy); - - MBinning binsimpact("BinningImpact"); - binsimpact.SetEdges(100, 0, 450); - parlist.AddToList(&binsimpact); - - parlist.AddToList(&tasklist); - - - // - // TaskList - // - MReadMarsFile reader("Events", filename); - reader.EnableBranch("fEnergy"); - reader.EnableBranch("fImpact"); - - - // ------------------------------------------------------------- - // - // Option 1. Just change the slope of the MC power law spectrum - // - //MMcWeightEnergySpecCalc wcalc(-2.0); //<-- Uncomment me - - // - // Option 2. A completely differente specturm pass as a TF1 function - // e.g. spectrum with exponential cutoff - // - //TF1 spec("spectrum","pow(x,[0])*exp(-x/[1])"); //<-- Uncomment me - //spec->SetParameter(0,-2.0); //<-- Uncomment me - //spec->SetParameter(1,50); //<-- Uncomment me - //MMcWeightEnergySpecCalc wcalc(spec); //<-- Uncomment me - - // - // Option 3. A completely differente specturm pass as a cahr* - // - //char* func = "pow(x,-2)"; //<-- Uncomment me - //MMcWeightEnergySpecCalc wcalc(func); //<-- Uncomment me - - // - // Option 4. A completely differente specturm pass as a c++ function - // - MMcWeightEnergySpecCalc wcalc((void*)myfunction); //<-- Uncomment me - // - //------------------------------------------------------------- - - MFillH hfill(&h1,"MMcEvt"); - MFillH hfill2(&h2,"MMcEvt"); - hfill2.SetWeight("MWeight"); - - tasklist.AddToList(&reader); - tasklist.AddToList(&wcalc); - tasklist.AddToList(&hfill); - tasklist.AddToList(&hfill2); - - // - // EventLoop - // - MEvtLoop magic; - magic.SetParList(&parlist); - - if (!magic.Eventloop()) - return; - - tasklist.PrintStatistics(); - parlist.Print(); - - // - // Draw the Results - // - TCanvas *c = new TCanvas(); - c->SetLogy(); - c->SetLogx(); - - TH1D* hist1 = (h1->GetHist())->ProjectionX(); - TH1D* hist2 = (h2->GetHist())->ProjectionX(); - hist2->SetLineColor(2); - - hist1->DrawClone(); - hist2->DrawClone("same"); -}