Index: /trunk/MagicSoft/Mars/macros/multidimdist2.C =================================================================== --- /trunk/MagicSoft/Mars/macros/multidimdist2.C (revision 1608) +++ /trunk/MagicSoft/Mars/macros/multidimdist2.C (revision 1608) @@ -0,0 +1,201 @@ +/* ======================================================================== *\ +! +! * +! * 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/2002 +! Author(s): Rudy Bock, 11/2002 +! +! Copyright: MAGIC Software Development, 2000-2002 +! +! +\* ======================================================================== */ + +void multidimdist2() +{ + // + // This is an example program which reads image parameters for gammas + // and hadrons, builds a 'matrix' of look-alike events, + // and then finds quality numbers and acceptances for the same sample + + // Create an empty Parameter List and an empty Task List + // The tasklist is identified in the eventloop by its name + // + MParList plist; + MTaskList tlistg; + plist.AddToList(&tlistg); + + // + // ------------- user attention ----------------- + // + + // cut used both for the matrices and the sample + // for more information see the documentation of MF and MDataChain + MF filterenergy("MMcEvt.fEnergy > 0"); + + // matrix limitation for look-alike events (approximate number) + MFEventSelector selector; + selector.SetNumSelectEvts(2000); + + // setup an AND-Filterlist from the two filters to be used + // in the event selection for the filling of the matrices + MFilterList flist; + flist.AddToList(&filterenergy); + flist.AddToList(&selector); + + // + // --------------------------------------------------------------- + // Now set up the tasks and tasklist (first event loop, gammas) + // --------------------------------------------------------------- + // + + // --------------- user change ----------------- + // Give the names of the star-files to be read + // Here you give the trainings sample(s) for + // the hadrons + // --------------------------------------------- + MReadMarsFile readg("Events", "star_gammas.root"); + readg.DisableAutoScheme(); + tlistg.AddToList(&readg); + + MHMatrix matrix("MatrixGammas"); + matrix.AddColumn("MHillas.fWidth"); + matrix.AddColumn("MHillas.fLength"); + matrix.AddColumn("MHillas.fWidth*MHillas.fLength/MHillas.fSize"); + matrix.AddColumn("abs(MHillas.fAsym)"); + matrix.AddColumn("abs(MHillas.fM3Long)"); + matrix.AddColumn("abs(MHillas.fM3Trans)"); + matrix.AddColumn("abs(MHillasSrc.fHeadTail)"); + matrix.AddColumn("MHillas.fConc"); + matrix.AddColumn("MHillas.fConc1"); + matrix.AddColumn("MHillasSrc.fDist"); + matrix.AddColumn("log10(MHillas.fSize)"); + plist.AddToList(&matrix); + + MFillH fillmatg("MatrixGammas"); + fillmatg.SetFilter(&flist); + tlistg.AddToList(&flist); + tlistg.AddToList(&fillmatg); + + // + // --- Create and set up the eventloop (gammas) --- + // + MEvtLoop evtloop; + evtloop.SetParList(&plist); + + // + // --- Execute matrix buildup (gammas) --- + // + if (!evtloop.Eventloop()) + return; + + tlistg.PrintStatistics(); + + // + // ------------------------------------------------------------------ + // prepare second event loop, hadrons + // ------------------------------------------------------------------ + // + MTaskList tlisth; + plist.Replace(&tlisth); + + // --------------- user change ----------------- + // Give the names of the star-files to be read + // Here you give the trainings sample(s) for + // the hadrons + // --------------------------------------------- + MReadMarsFile readh("Events", "star_protons.root"); + readh.DisableAutoScheme(); + tlisth.AddToList(&readh); + + MHMatrix matrixh("MatrixHadrons"); + matrixh.AddColumns(matrix.GetColumns()); + plist.AddToList(&matrixh); + + MFillH fillmath("MatrixHadrons"); + fillmath.SetFilter(&flist); // filter list + tlisth.AddToList(&flist); + tlisth.AddToList(&fillmath); + + // + // Create and set up the eventloop (protons) + // + MEvtLoop evtloop; + evtloop.SetParList(&plist); + + // + // Execute matrix buildup (hadrons) + // + if (!evtloop.Eventloop()) + return; + + // sum up in log + tlistg.PrintStatistics(); + + matrix.Print("size"); + matrixh.Print("size"); + + // + // ---------------------------------------------------------- + // Go through full sample again, now for getting hadronness + // (third event loop) + // ---------------------------------------------------------- + // + + MTaskList tlist2; + plist.Replace(&tlist2); + + // ------------------- user change -------------------- + // Give the names of the star-files to be read as + // test samples you should at least have one hadron + // and one gamma file available + // ---------------------------------------------------- + MReadMarsFile read2("Events", "star_protons.root"); + read2.AddFile("star_gammas.root"); + read2.DisableAutoScheme(); + tlist2.AddToList(&read2); + + // ---------------- user attention ----------------- + // Here you may change the algorithm used + // You can switch from Kernel to Next Neighbor + // or change the number of used shortest distances + // For kernel you should always set this to 0 (all) + // ------------------------------------------------- + MMultiDimDistCalc calc; + calc.SetUseNumRows(25); + calc.SetUseKernelMethod(kFALSE); + tlist2.AddToList(&calc); + + MFillH fillh("MHHadronness"); + fillh.SetFilter(&filterenergy); + tlist2.AddToList(&filterenergy); + tlist2.AddToList(&fillh); + + // + // Execute analysis of gammas and hadrons + // + MProgressBar bar; + MEvtLoop evtloop; + evtloop.SetProgressBar(&bar); + evtloop.SetParList(&plist); + + if (!evtloop.Eventloop()) + return; + + tlist2.PrintStatistics(); + + plist.FindObject("MHHadronness")->DrawClone(); + plist.FindObject("MHHadronness")->Print(); +} Index: /trunk/MagicSoft/Mars/macros/star.C =================================================================== --- /trunk/MagicSoft/Mars/macros/star.C (revision 1607) +++ /trunk/MagicSoft/Mars/macros/star.C (revision 1608) @@ -35,6 +35,6 @@ { // - // This is a demonstration program which calculates the Hillas - // parameter out of a Magic root file. + // This is a demonstration program which calculates the image + // parameters from a Magic raw data root file. // @@ -47,4 +47,7 @@ plist.AddToList(&tlist); + MSrcPosCam src; + plist.AddToList(&src); + // // Uncomment this two line if you want to use MHillasExt instead @@ -56,22 +59,8 @@ // // The geometry container must be created by yourself to make sure - // that you don't choos a wrong geometry by chance + // that you don't choose a wrong geometry by mistake // MGeomCamMagic geomcam; plist.AddToList(&geomcam); - - // - // Craete the object which hlods the source positions in the camera - // plain in respect to which the image parameters will be calculated. - // For real data the containers will be filled by a task. - // - MSrcPosCam source("Source") - source.SetXY(0, 0); - - MSrcPosCam antisrc("AntiSrc"); - antisrc.SetXY(240, 0); - - plist.AddToList(&source); - plist.AddToList(&antisrc); // @@ -80,6 +69,9 @@ // MReadMarsFile read("Events"); - read.AddFile("Gammas*.root"); read.DisableAutoScheme(); + + // ------------- user change ----------------- + read.AddFile("Pro*.root"); + //read.AddFile("Gam*.root"); MMcPedestalCopy pcopy; @@ -92,16 +84,14 @@ MImgCleanStd clean; MHillasCalc hcalc; - MHillasSrcCalc csrc1("Source", "HillasSource"); - MHillasSrcCalc csrc2("AntiSrc", "HillasAntiSrc"); + MHillasSrcCalc scalc; // !!Preliminary!! Will be removed later! - MWriteRootFile write("star.root"); + // ------------- user change ----------------- + MWriteRootFile write("data/star_protons.root"); write.AddContainer("MHillas", "Events"); write.AddContainer("HillasSource", "Events"); - write.AddContainer("HillasAntiSrc", "Events"); write.AddContainer("MMcEvt", "Events"); + write.AddContainer("MHillasSrc", "Events"); write.AddContainer("MRawRunHeader", "RunHeaders"); write.AddContainer("MMcRunHeader", "RunHeaders"); - write.AddContainer("Source", "RunHeaders"); - write.AddContainer("AntiSource", "RunHeaders"); tlist.AddToList(&read); @@ -112,10 +102,9 @@ tlist.AddToList(&clean); tlist.AddToList(&hcalc); - tlist.AddToList(&csrc1); - tlist.AddToList(&csrc2); + tlist.AddToList(&scalc); tlist.AddToList(&write); // - // Create and setup the eventloop + // Create and set up the eventloop // MProgressBar bar;