Index: trunk/MagicSoft/Mars/Changelog =================================================================== --- trunk/MagicSoft/Mars/Changelog (revision 2981) +++ trunk/MagicSoft/Mars/Changelog (revision 2982) @@ -9,4 +9,12 @@ - changed default FADC integration range (more adequate for coming camera version). + + * macros/mccalibrate.C, starmc2.C + - added. First macro transforms input files from camera simulation + into files with pixel data already calibrated in photons. The + second reads in an already calibrated data file and produces an + output file with Hillas parameters. Both together do the same as + starmc.C, but it is convenient to have files containing + calibrated pixel information. 2004/01/30: Thomas Bretz Index: trunk/MagicSoft/Mars/macros/mccalibrate.C =================================================================== --- trunk/MagicSoft/Mars/macros/mccalibrate.C (revision 2982) +++ trunk/MagicSoft/Mars/macros/mccalibrate.C (revision 2982) @@ -0,0 +1,195 @@ +/* ======================================================================== *\ + ! + ! * + ! * 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 1/2004 + ! Thomas Bretz 5/2002 + ! + ! Copyright: MAGIC Software Development, 2000-2004 + ! + ! + \* ======================================================================== */ + +///////////////////////////////////////////////////////////////////////////// +// +// MMCALIBRATE - Calibration of MC data +// +// This macro is a version of the standard procedure to convert raw MC data +// into calibrated data (photons per pixel) +// +///////////////////////////////////////////////////////////////////////////// + +#include "MImgCleanStd.h" + +void mccalibrate() +{ + // + // This is a demonstration program which reads in MC camera files + // and produces and output with calibrated events (signal in photons + // for all pixels, in MCerPhotEvt containers). + // + + + // ------------- user change ----------------- + TString* CalibrationFilename; + CalibrationFilename = new TString("nonoise/Gamma_zbin0_0*.root"); + // File to be used for the calibration (must be a camera file without added noise) + + Char_t* AnalysisFilename = "Gamma_zbin0*.root"; // File to be analyzed + + Char_t* OutFilename = "calibrated_data.root"; // Output file name + + + Int_t BinsHigh[2] = {5, 9}; // First and last FADC bin of the range to be integrated, + Int_t BinsLow[2] = {5, 9}; // for high and low gain respectively. + + // ------------------------------------------- + + // + // 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); + + MSrcPosCam src; + src.SetReadyToSave(); + + plist.AddToList(&src); + + // + // Now setup the tasks and tasklist: + // --------------------------------- + // + MReadMarsFile read("Events"); + + if (CalibrationFilename) + read.AddFile(CalibrationFilename->Data()); + + read.DisableAutoScheme(); + + MGeomApply geom; // Reads in geometry from MC file and sets the right sizes for + // several parameter containers. + + MMcPedestalCopy pcopy; + // Copies pedestal data from the MC file run fadc header to the MPedestalCam container. + + MExtractSignal sigextract; + // Define ADC slices to be integrated in high and low gain: + sigextract.SetRange(BinsHigh[0], BinsHigh[1], BinsLow[0], BinsLow[1]); + + MMcCalibrationUpdate mccalibupdate; + + MCalibrate calib; // Transforms signals from ADC counts into photons. + + // MBlindPixelCalc blind; + // blind.SetUseInterpolation(); + + // + // Applies tail cuts to image. Since the calibration is performed on + // noiseless camera files, the precise values of the cleaning levels + // are unimportant (in any case, only pixels without any C-photon will + // be rejected). + // + MImgCleanStd clean; + + + MHillasCalc hcalc; // Calculates Hillas parameters not dependent on source position. + + MMcCalibrationCalc mccalibcalc; + + tlist.AddToList(&read); + tlist.AddToList(&geom); + tlist.AddToList(&pcopy); + + tlist.AddToList(&sigextract); + tlist.AddToList(&mccalibupdate); + tlist.AddToList(&calib); + tlist.AddToList(&clean); + // tlist.AddToList(&blind); + tlist.AddToList(&hcalc); + + tlist.AddToList(&mccalibcalc); + + // + // Open output file: + // + MWriteRootFile write(OutFilename); // Writes output + write.AddContainer("MGeomCam", "RunHeaders"); + write.AddContainer("MMcConfigRunHeader", "RunHeaders"); + write.AddContainer("MMcCorsikaRunHeader", "RunHeaders"); + write.AddContainer("MMcFadcHeader", "RunHeaders"); + write.AddContainer("MMcRunHeader", "RunHeaders"); + write.AddContainer("MMcTrigHeader", "RunHeaders"); + write.AddContainer("MRawRunHeader", "RunHeaders"); + write.AddContainer("MSrcPosCam", "RunHeaders"); + + + write.AddContainer("MMcEvt", "Events"); + write.AddContainer("MMcTrig", "Events"); + write.AddContainer("MRawEvtHeader", "Events"); + write.AddContainer("MCerPhotEvt", "Events"); + write.AddContainer("MPedPhotCam", "Events"); + + // + // First loop: Calibration loop + // + + MProgressBar bar; + bar.SetWindowName("Calibrating..."); + + MEvtLoop evtloop; + evtloop.SetProgressBar(&bar); + evtloop.SetParList(&plist); + + if (CalibrationFilename) + { + if (!evtloop.Eventloop()) + return; + mccalibcalc->GetHist()->Write(); + } + + // + // Second loop: analysis loop + // + + // + // Change the read task by another one which reads the file we want to analyze: + // + + MReadMarsFile read2("Events"); + read2.AddFile(AnalysisFilename); + read2.DisableAutoScheme(); + tlist.AddToListBefore(&read2, &read, "All"); + tlist.RemoveFromList(&read); + + bar.SetWindowName("Writing..."); + + tlist.RemoveFromList(&clean); + tlist.RemoveFromList(&hcalc); + tlist.RemoveFromList(&mccalibcalc); + + tlist.AddToList(&write); // Add task to write output. + + if (!evtloop.Eventloop()) + return; + + + tlist.PrintStatistics(); +} Index: trunk/MagicSoft/Mars/macros/starmc2.C =================================================================== --- trunk/MagicSoft/Mars/macros/starmc2.C (revision 2982) +++ trunk/MagicSoft/Mars/macros/starmc2.C (revision 2982) @@ -0,0 +1,117 @@ +/* ======================================================================== *\ + ! + ! * + ! * 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 1/2004 + ! Thomas Bretz 5/2002 + ! + ! Copyright: MAGIC Software Development, 2000-2004 + ! + ! + \* ======================================================================== */ + +///////////////////////////////////////////////////////////////////////////// +// +// STARMC2 - STandard Analysis and Reconstruction (MC example) +// +// This macro converts into image parameters an input file of MC data +// previously calibrated (see mccalibrate.C). +// +// +///////////////////////////////////////////////////////////////////////////// + +#include "MImgCleanStd.h" + +void starmc2() +{ + Char_t* AnalysisFilename = "calibrated_data.root"; // File to be analyzed + Char_t* OutFilename = "star_mc.root"; // Output file name + + Float_t CleanLev[2] = {4., 3.}; // Tail cuts for image analysis + + Int_t BinsHigh[2] = {5, 9}; // First and last FADC bin of the range to be integrated, + Int_t BinsLow[2] = {5, 9}; // for high and low gain respectively. + + // ------------------------------------------------------------------ + + // + // 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); + + MSrcPosCam src; + src.SetReadyToSave(); + + plist.AddToList(&src); + + // + // Now setup the tasks and tasklist: + // --------------------------------- + // + MReadMarsFile read("Events"); + + read.AddFile(AnalysisFilename); + + read.DisableAutoScheme(); + + MImgCleanStd clean(CleanLev[0], CleanLev[1]); // Applies tail cuts to image. + + MHillasCalc hcalc; // Calculates Hillas parameters not dependent on source position. + MHillasSrcCalc scalc; // Calculates source-dependent Hillas parameters + + tlist.AddToList(&read); + tlist.AddToList(&clean); + // tlist.AddToList(&blind); + tlist.AddToList(&hcalc); + tlist.AddToList(&scalc); // Calculates Source-dependent Hillas parameters + + // + // Open output file: + // + MWriteRootFile write(OutFilename); // Writes output + write.AddContainer("MRawRunHeader", "RunHeaders"); + write.AddContainer("MMcRunHeader", "RunHeaders"); + write.AddContainer("MSrcPosCam", "RunHeaders"); + write.AddContainer("MMcEvt", "Events"); + write.AddContainer("MHillas", "Events"); + write.AddContainer("MHillasExt", "Events"); + write.AddContainer("MHillasSrc", "Events"); + write.AddContainer("MNewImagePar", "Events"); + + tlist.AddToList(&write); // Add task to write output. + + // + // analysis loop + // + + MEvtLoop evtloop; + MProgressBar bar; + bar.SetWindowName("Analyzing..."); + evtloop.SetProgressBar(&bar); + evtloop.SetParList(&plist); + + if (!evtloop.Eventloop()) + return; + + tlist.PrintStatistics(); + + return; +}