| 1 | #include <TApplication.h> | 
|---|
| 2 |  | 
|---|
| 3 | #include <TThread.h> | 
|---|
| 4 | #include <TCanvas.h> | 
|---|
| 5 | #include <TMethod.h>          // GetMenuItems | 
|---|
| 6 |  | 
|---|
| 7 | #include <TGButton.h> | 
|---|
| 8 |  | 
|---|
| 9 | #include "MParList.h" | 
|---|
| 10 | #include "MTaskList.h" | 
|---|
| 11 | #include "MEvtLoop.h" | 
|---|
| 12 |  | 
|---|
| 13 | #include "MLogManip.h" | 
|---|
| 14 |  | 
|---|
| 15 | #include "MRawSocketRead.h" | 
|---|
| 16 | #include "MGeomApply.h" | 
|---|
| 17 | #include "MCerPhotAnal2.h" | 
|---|
| 18 | #include "MFillH.h" | 
|---|
| 19 | #include "MFRealTimePeriod.h" | 
|---|
| 20 | #include "MHTriggerLvl0.h" | 
|---|
| 21 | #include "MHCamera.h"         // ::Class(), SetFreezed | 
|---|
| 22 | #include "MHCamEvent.h"       // MHCamEvent | 
|---|
| 23 | #include "MHEvent.h"          // MHEvent::GetHist() | 
|---|
| 24 |  | 
|---|
| 25 | #include "MOnlineDump.h" | 
|---|
| 26 | #include "MOnlineDisplay.h" | 
|---|
| 27 |  | 
|---|
| 28 | #include "MImgCleanStd.h" | 
|---|
| 29 | #include "MHillasCalc.h" | 
|---|
| 30 | #include "MHillasSrcCalc.h" | 
|---|
| 31 | #include "MEventRateCalc.h" | 
|---|
| 32 |  | 
|---|
| 33 | #include "MArgs.h" | 
|---|
| 34 |  | 
|---|
| 35 | // -------------------------------------------------------------------- | 
|---|
| 36 |  | 
|---|
| 37 | #include "MTime.h" | 
|---|
| 38 | #include "MRawEvtHeader.h"       // MRawEvtHeader | 
|---|
| 39 | #include "MRawRunHeader.h"       // MRawRunHeader | 
|---|
| 40 |  | 
|---|
| 41 | using namespace std; | 
|---|
| 42 |  | 
|---|
| 43 | void StartUpMessage() | 
|---|
| 44 | { | 
|---|
| 45 | gLog << all << endl; | 
|---|
| 46 |  | 
|---|
| 47 | //                1         2         3         4         5 | 
|---|
| 48 | //       12345678901234567890123456789012345678901234567890 | 
|---|
| 49 | gLog << "==================================================" << endl; | 
|---|
| 50 | gLog << "                    MONA V" << MARSVER << "                      " << endl; | 
|---|
| 51 | gLog << "              Magic Online Analysis               " << endl; | 
|---|
| 52 | gLog << "            Compiled on <" << __DATE__ << ">"       << endl; | 
|---|
| 53 | gLog << "               Using ROOT v" << ROOTVER             << endl; | 
|---|
| 54 | gLog << "==================================================" << endl; | 
|---|
| 55 | gLog << endl; | 
|---|
| 56 | } | 
|---|
| 57 |  | 
|---|
| 58 | void Remove(TMethod *m, const char *name) | 
|---|
| 59 | { | 
|---|
| 60 | if (TString(m->GetName()).BeginsWith(name)) | 
|---|
| 61 | m->SetMenuItem(kMenuNoMenu); | 
|---|
| 62 | } | 
|---|
| 63 |  | 
|---|
| 64 | void ChangeContextMenus() | 
|---|
| 65 | { | 
|---|
| 66 | TList l; | 
|---|
| 67 | MHCamera::Class()->GetMenuItems(&l); | 
|---|
| 68 | TIter Next(&l); | 
|---|
| 69 | TMethod *m=NULL; | 
|---|
| 70 | while ((m=(TMethod*)Next())) | 
|---|
| 71 | { | 
|---|
| 72 | Remove(m, "DrawClone"); | 
|---|
| 73 | Remove(m, "SetName"); | 
|---|
| 74 | Remove(m, "Delete"); | 
|---|
| 75 | Remove(m, "DrawClass"); | 
|---|
| 76 | Remove(m, "Dump"); | 
|---|
| 77 | Remove(m, "Inspect"); | 
|---|
| 78 | Remove(m, "Smooth"); | 
|---|
| 79 | Remove(m, "Fit"); | 
|---|
| 80 | Remove(m, "Divide"); | 
|---|
| 81 | Remove(m, "Add"); | 
|---|
| 82 | Remove(m, "Multiply"); | 
|---|
| 83 | Remove(m, "Delete"); | 
|---|
| 84 | Remove(m, "SetLineAttributes"); | 
|---|
| 85 | Remove(m, "SetFillAttributes"); | 
|---|
| 86 | Remove(m, "SetMarkerAttributes"); | 
|---|
| 87 | Remove(m, "SetDrawOption"); | 
|---|
| 88 | } | 
|---|
| 89 | } | 
|---|
| 90 |  | 
|---|
| 91 | #include "MPedPhotCam.h" | 
|---|
| 92 | Bool_t Loop(MOnlineDisplay *display, Int_t port) | 
|---|
| 93 | { | 
|---|
| 94 | display->Reset(); | 
|---|
| 95 |  | 
|---|
| 96 | // | 
|---|
| 97 | // create the tasks which should be executed and add them to the list | 
|---|
| 98 | // in the case you don't need parameter containers, all of them can | 
|---|
| 99 | // be created by MRawFileRead::PreProcess | 
|---|
| 100 | // | 
|---|
| 101 | MRawSocketRead reader; | 
|---|
| 102 | reader.SetPort(port); | 
|---|
| 103 |  | 
|---|
| 104 | // | 
|---|
| 105 | // create a (empty) list of parameters which can be used by the tasks | 
|---|
| 106 | // and an (empty) list of tasks which should be executed | 
|---|
| 107 | // | 
|---|
| 108 | MParList plist; | 
|---|
| 109 |  | 
|---|
| 110 | MTaskList tasks; | 
|---|
| 111 | plist.AddToList(&tasks); | 
|---|
| 112 | tasks.AddToList(&reader); | 
|---|
| 113 |  | 
|---|
| 114 | MPedPhotCam pcam; | 
|---|
| 115 | plist.AddToList(&pcam); | 
|---|
| 116 |  | 
|---|
| 117 | MGeomApply geomapl; | 
|---|
| 118 | MCerPhotAnal2 ncalc; | 
|---|
| 119 | tasks.AddToList(&geomapl); | 
|---|
| 120 | tasks.AddToList(&ncalc); | 
|---|
| 121 |  | 
|---|
| 122 | MEventRateCalc tcalc; | 
|---|
| 123 | tcalc.SetNumEvents(100); | 
|---|
| 124 | tasks.AddToList(&tcalc); | 
|---|
| 125 |  | 
|---|
| 126 | MFillH fill1("MHEvent",    "MCerPhotEvt", "MFillEvent"); | 
|---|
| 127 | MFillH fill2("MHCamEvent", "MCerPhotEvt", "MFillCamEvent"); | 
|---|
| 128 |  | 
|---|
| 129 | MFRealTimePeriod filter; | 
|---|
| 130 | fill1.SetFilter(&filter); | 
|---|
| 131 |  | 
|---|
| 132 | MHTriggerLvl0 trigmap(128, "Above 128"); | 
|---|
| 133 | MFillH fill3(&trigmap, "MRawEvtData", "MFillTriggerLvl0"); | 
|---|
| 134 |  | 
|---|
| 135 | tasks.AddToList(&filter); | 
|---|
| 136 |  | 
|---|
| 137 | tasks.AddToList(&fill1); | 
|---|
| 138 | tasks.AddToList(&fill2); | 
|---|
| 139 | tasks.AddToList(&fill3); | 
|---|
| 140 |  | 
|---|
| 141 | MHCamEvent hist("PedestalRms"); | 
|---|
| 142 | hist.SetType(1); | 
|---|
| 143 | plist.AddToList(&hist); | 
|---|
| 144 |  | 
|---|
| 145 | // ------------------------------------------- | 
|---|
| 146 |  | 
|---|
| 147 | MHCamEvent maxhi("MaxIdxHi", "Index of slice with maximum content (hi-gain)"); | 
|---|
| 148 | MHCamEvent maxlo("MaxIdxLo", "Index of slice with maximum content (lo-gain)"); | 
|---|
| 149 | maxhi.SetType(3); | 
|---|
| 150 | maxlo.SetType(4); | 
|---|
| 151 | plist.AddToList(&maxhi); | 
|---|
| 152 | plist.AddToList(&maxlo); | 
|---|
| 153 |  | 
|---|
| 154 | // ------------------------------------------- | 
|---|
| 155 |  | 
|---|
| 156 | MFillH hfilla("MaxIdxHi", "MRawEvtData", "FillMaxIdxHi"); | 
|---|
| 157 | MFillH hfillb("MaxIdxLo", "MRawEvtData", "FillMaxIdxLo"); | 
|---|
| 158 | MFillH hfillc("Pedestals [MHCamEvent]", "MPedPhotCam", "FillPedestal"); | 
|---|
| 159 | MFillH hfilld("PedestalRms", "MPedPhotCam", "FillPedestalRms"); | 
|---|
| 160 | tasks.AddToList(&hfilla); | 
|---|
| 161 | tasks.AddToList(&hfillb); | 
|---|
| 162 | tasks.AddToList(&hfillc); | 
|---|
| 163 | tasks.AddToList(&hfilld); | 
|---|
| 164 |  | 
|---|
| 165 | MOnlineDump dump; | 
|---|
| 166 | dump.SetFilter(&filter); | 
|---|
| 167 | display->SetTask(&dump); | 
|---|
| 168 | tasks.AddToList(&dump); | 
|---|
| 169 |  | 
|---|
| 170 |  | 
|---|
| 171 | MImgCleanStd      clean; | 
|---|
| 172 | MHillasCalc       hcalc; | 
|---|
| 173 | MHillasSrcCalc    scalc; // !!Preliminary!! Will be removed later! | 
|---|
| 174 | MFillH hfill3("MHHillas",   "MHillas",    "MFillHillas"); | 
|---|
| 175 | MFillH hfill4("MHHillasSrc","MHillasSrc", "MFillHillasSrc"); | 
|---|
| 176 | tasks.AddToList(&clean); | 
|---|
| 177 | tasks.AddToList(&hcalc); | 
|---|
| 178 | tasks.AddToList(&scalc); | 
|---|
| 179 | tasks.AddToList(&hfill3); | 
|---|
| 180 | tasks.AddToList(&hfill4); | 
|---|
| 181 |  | 
|---|
| 182 | MEvtLoop magic; | 
|---|
| 183 | magic.SetParList(&plist); | 
|---|
| 184 |  | 
|---|
| 185 | magic.SetDisplay(display); | 
|---|
| 186 | magic.SetProgressBar((TGProgressBar*)NULL); | 
|---|
| 187 |  | 
|---|
| 188 | if (!magic.Eventloop()) | 
|---|
| 189 | { | 
|---|
| 190 | gLog << err << "Mona Eventloop error..." << endl; | 
|---|
| 191 | return kFALSE; | 
|---|
| 192 | } | 
|---|
| 193 |  | 
|---|
| 194 | tasks.PrintStatistics(); | 
|---|
| 195 | return kTRUE; | 
|---|
| 196 | } | 
|---|
| 197 |  | 
|---|
| 198 | // | 
|---|
| 199 | // By defining the function return type 'void' instead of | 
|---|
| 200 | // 'void*' we tell TThread to create a detached thread, | 
|---|
| 201 | // but calling Exit() in a detached thread results in a crash | 
|---|
| 202 | // when the TGApplication is terminated. | 
|---|
| 203 | // | 
|---|
| 204 | void *thread(void *ptr) | 
|---|
| 205 | { | 
|---|
| 206 | const Int_t port = *((Int_t**)ptr)[0]; | 
|---|
| 207 | MOnlineDisplay &d = *((MOnlineDisplay**)ptr)[1]; | 
|---|
| 208 |  | 
|---|
| 209 | // | 
|---|
| 210 | // If a loop is stopped reinitialize a new loop until the | 
|---|
| 211 | // user requested to exit the program. | 
|---|
| 212 | // | 
|---|
| 213 | while (d.CheckStatus()!=MStatusDisplay::kFileExit) | 
|---|
| 214 | Loop(&d, port); | 
|---|
| 215 |  | 
|---|
| 216 |  | 
|---|
| 217 | gLog << dbg << "Loop stopped... exit thread()" << endl; | 
|---|
| 218 |  | 
|---|
| 219 | return 0; | 
|---|
| 220 | } | 
|---|
| 221 |  | 
|---|
| 222 | /* | 
|---|
| 223 | #if ROOT_VERSION_CODE < ROOT_VERSION(3,10,02) | 
|---|
| 224 | #include <TRootGuiFactory.h> | 
|---|
| 225 | #include <TPluginManager.h> | 
|---|
| 226 | void InitGuiFactory() | 
|---|
| 227 | { | 
|---|
| 228 | if (gROOT->IsBatch()) | 
|---|
| 229 | gROOT->SetBatch(kFALSE); | 
|---|
| 230 |  | 
|---|
| 231 | // | 
|---|
| 232 | // Must be loaded by hand, because it is not loaded by TGApplication. | 
|---|
| 233 | // We could also use TApplication instead, but TApplication doesn't | 
|---|
| 234 | // deal with the DISPLAY variable in a convient way. | 
|---|
| 235 | // | 
|---|
| 236 | TPluginHandler *h; | 
|---|
| 237 | if ((h = gROOT->GetPluginManager()->FindHandler("TGuiFactory", "root"))) | 
|---|
| 238 | { | 
|---|
| 239 | if (h->LoadPlugin() == -1) | 
|---|
| 240 | return; | 
|---|
| 241 | gGuiFactory = (TGuiFactory*)h->ExecPlugin(0); | 
|---|
| 242 | } | 
|---|
| 243 | } | 
|---|
| 244 | #endif | 
|---|
| 245 | */ | 
|---|
| 246 |  | 
|---|
| 247 | static void Usage() | 
|---|
| 248 | { | 
|---|
| 249 | gLog << all << endl; | 
|---|
| 250 | gLog << "Sorry the usage is:" << endl; | 
|---|
| 251 | gLog << "   mona [-?] [-h] [-a0] [-h] [-pn] [-vn]" << endl << endl; | 
|---|
| 252 | gLog << "     -a0: Do not use Ansii codes." << endl; | 
|---|
| 253 | gLog << "     -pn: TCP/IP port n [default=7000]" << endl; | 
|---|
| 254 | gLog << "     -vn: Verbosity level n [default=2]" << endl << endl; | 
|---|
| 255 | gLog << "     -?/-h: This help" << endl << endl; | 
|---|
| 256 | } | 
|---|
| 257 |  | 
|---|
| 258 |  | 
|---|
| 259 | int main(int argc, char **argv) | 
|---|
| 260 | { | 
|---|
| 261 | StartUpMessage(); | 
|---|
| 262 |  | 
|---|
| 263 | // | 
|---|
| 264 | // Evaluate arguments | 
|---|
| 265 | // | 
|---|
| 266 | MArgs arg(argc, argv); | 
|---|
| 267 |  | 
|---|
| 268 | if (arg.HasOption("-?") || arg.HasOption("-h")) | 
|---|
| 269 | { | 
|---|
| 270 | Usage(); | 
|---|
| 271 | return -1; | 
|---|
| 272 | } | 
|---|
| 273 |  | 
|---|
| 274 | // | 
|---|
| 275 | // Set verbosity to highest level. | 
|---|
| 276 | // | 
|---|
| 277 | gLog.SetDebugLevel(arg.HasOption("-v") ? arg.GetIntAndRemove("-v") : 2); | 
|---|
| 278 |  | 
|---|
| 279 | if (arg.HasOption("-a") && arg.GetIntAndRemove("-a")==0) | 
|---|
| 280 | gLog.SetNoColors(); | 
|---|
| 281 |  | 
|---|
| 282 | const Int_t port = arg.HasOption("-p") ? arg.GetIntAndRemove("-p") : 7000; | 
|---|
| 283 |  | 
|---|
| 284 | // | 
|---|
| 285 | // check for the right usage of the program | 
|---|
| 286 | // | 
|---|
| 287 | if (arg.GetNumArguments()>0) | 
|---|
| 288 | { | 
|---|
| 289 | Usage(); | 
|---|
| 290 | return -1; | 
|---|
| 291 | } | 
|---|
| 292 |  | 
|---|
| 293 | TApplication app("Mona", &argc, argv); | 
|---|
| 294 | if (gROOT->IsBatch() || !gClient) | 
|---|
| 295 | { | 
|---|
| 296 | gLog << "Bombing... maybe your DISPLAY variable is not set correctly!" << endl; | 
|---|
| 297 | return 1; | 
|---|
| 298 | } | 
|---|
| 299 |  | 
|---|
| 300 | /* | 
|---|
| 301 | TGApplication app("Mona", &argc, argv); | 
|---|
| 302 |  | 
|---|
| 303 | #if ROOT_VERSION_CODE < ROOT_VERSION(3,10,02) | 
|---|
| 304 | InitGuiFactory(); | 
|---|
| 305 | #endif | 
|---|
| 306 | */ | 
|---|
| 307 |  | 
|---|
| 308 | ChangeContextMenus(); | 
|---|
| 309 |  | 
|---|
| 310 | // | 
|---|
| 311 | // Starting MStatusDisplay in the thread results in hangs | 
|---|
| 312 | // if the thread is terminated (by return) | 
|---|
| 313 | // | 
|---|
| 314 | gLog << dbg << "Starting Display..." << flush; | 
|---|
| 315 | MOnlineDisplay d; | 
|---|
| 316 | d.SetWindowName("Magic ONline Analysis"); | 
|---|
| 317 | d.SetIconName("Mona"); | 
|---|
| 318 |  | 
|---|
| 319 | d.SetBit(MStatusDisplay::kExitLoopOnExit); | 
|---|
| 320 | gLog << "done." << endl; | 
|---|
| 321 |  | 
|---|
| 322 | // gDebug=1; | 
|---|
| 323 |  | 
|---|
| 324 | gLog << dbg << "Starting Thread..." << flush; | 
|---|
| 325 | const void *ptr[2] = { &port, &d }; | 
|---|
| 326 | TThread t(thread, ptr); | 
|---|
| 327 | t.Run(); | 
|---|
| 328 | gLog << "done." << endl; | 
|---|
| 329 |  | 
|---|
| 330 | gLog << dbg << "Starting System loop... " << endl; | 
|---|
| 331 | app.Run(kTRUE); | 
|---|
| 332 | gLog << dbg << "System loop stopped." << endl; | 
|---|
| 333 |  | 
|---|
| 334 | d.UnmapWindow(); | 
|---|
| 335 |  | 
|---|
| 336 | gLog << dbg << "Waiting for termination of thread... (be patient)" << endl; | 
|---|
| 337 | while (t.GetState()!=TThread::kCanceledState) | 
|---|
| 338 | gSystem->ProcessEvents(); | 
|---|
| 339 |  | 
|---|
| 340 | gLog << dbg << "Thread terminated... joining." << endl; | 
|---|
| 341 |  | 
|---|
| 342 | TThread::Join(t.GetId()); //- seems that it is implicitly done... why? | 
|---|
| 343 |  | 
|---|
| 344 | gLog << dbg << "Exit MONA." << endl; | 
|---|
| 345 |  | 
|---|
| 346 | return 0; | 
|---|
| 347 | } | 
|---|