Changeset 4756


Ignore:
Timestamp:
08/27/04 14:45:08 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 added
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r4755 r4756  
    2929       been set (e.g. by the macro).
    3030
    31    
    3231
    3332
    … …  
    3635   * Makefile:
    3736     - added comments how to link statically
     37     - added star again
    3838
    3939   * callisto.cc:
    … …  
    6969       instead of accesing the arrays many times inside the loops.
    7070       Somebody should do the same for logain.
     71
     72   * mraw/MRawEvtData.cc:
     73     - fixed a bug in InitArrays
     74
     75   * showlog.cc:
     76     - added some more functionality
     77
     78   * star.cc:
     79     - changed to use new MJStar
     80
     81   * mfileio/MChain.cc:
     82     - added namespace
     83
     84   * mfileio/MReadTree.cc:
     85     - added a call to ResetTree as first call in PreProcess. Using
     86       MReadReports with a tree missing in the file mysteriously
     87       crashed because TChain::fTreeNumber was set already by
     88       a root-miracle.
     89
     90   * mjobs/JobsLinkDef.h, mjobs/Makefile:
     91     - added MJStar
     92
     93   * mjobs/MJCalibrateSignal.cc:
     94     - replaced MReadMarsFile by MReadReports
     95
     96   * mmain/MEventDisplay.cc:
     97     - fixed reading of resource file
    7198
    7299

  • trunk/MagicSoft/Mars/Makefile

    r4752 r4756  
    2121
    2222#PROGRAMS = readraw merpp mars test mona status
    23 PROGRAMS = readdaq readraw merpp mars mona showlog callisto showplot
     23PROGRAMS = readdaq readraw merpp mars mona showlog callisto showplot star
    2424SOLIB    = libmars.so
    2525CINT     = M

  • trunk/MagicSoft/Mars/mfileio/MChain.cc

    r2173 r4756  
    1818!   Author(s): Thomas Bretz, 12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
    1919!
    20 !   Copyright: MAGIC Software Development, 2000-2003
     20!   Copyright: MAGIC Software Development, 2000-2004
    2121!
    2222!
    … …  
    2424
    2525/////////////////////////////////////////////////////////////////////////////
    26 //                                                                         //
    27 // MChain                                                                  //
    28 //                                                                         //
    29 // Helper class for MReadTree                                              //
    30 //                                                                         //
     26//
     27// MChain
     28//
     29// Helper class for MReadTree
     30//
    3131/////////////////////////////////////////////////////////////////////////////
    3232#include "MChain.h"
    3333
    3434ClassImp(MChain);
     35
     36using namespace std;
    3537
    3638Int_t MChain::LoadTree(Int_t entry)

  • trunk/MagicSoft/Mars/mfileio/MReadTree.cc

    r3682 r4756  
    609609
    610610    //
     611    // It seems, that TFile and TTree are not completely independant if
     612    // many times the same file is openes (MReadReports) and some of
     613    // the files in the chains don't have one tree. TChain::fTreeNumber
     614    // is already set before LoadTree from GetFile is called and
     615    // it crashes. ResetTree makes sure, that the tree number is -1
     616    //
     617    fChain->ResetTree();
     618
     619    //
    611620    // check for files and for the tree!
    612621    //

  • trunk/MagicSoft/Mars/mmain/MEventDisplay.cc

    r4722 r4756  
    3737#include <TStyle.h>              // gStyle->SetOptStat
    3838#include <TCanvas.h>             // TCanvas::cd
     39#include <TSystem.h>             // TSystem::BaseName
    3940
    4041//
    … …  
    245246        plist->AddToList(ccam);
    246247
    247     fEvtLoop = new MEvtLoop;
    248     fEvtLoop->SetOwner();
    249     fEvtLoop->SetParList(plist);
    250     fEvtLoop->ReadEnv("mars.rc");
    251 
    252248    MHEvent *evt01 = new MHEvent(MHEvent::kEvtSignalRaw);
    253249    MHEvent *evt02 = new MHEvent(MHEvent::kEvtSignalRaw);
    … …  
    391387    // Now distribute Display to all tasks
    392388    //
    393     tlist->SetDisplay(this);
     389    fEvtLoop = new MEvtLoop(Form("Mars %s", gSystem->BaseName(fname)));
     390    fEvtLoop->SetOwner();
     391    fEvtLoop->SetParList(plist);
     392    fEvtLoop->SetDisplay(this);
     393    fEvtLoop->ReadEnv("mars.rc");
     394
    394395}
    395396

  • trunk/MagicSoft/Mars/showlog.cc

    r4732 r4756  
    22#include <fstream>
    33
     4#include "MArgs.h"
     5
    46#include "MLog.h"
     7#include "MLogManip.h"
    58
    69using namespace std;
    … …  
    1013    gLog << endl;
    1114    gLog << "Sorry the usage is:" << endl;
    12     gLog << "   showlog filename" << endl;
     15    gLog << "   showlog [options] filename" << endl;
    1316    gLog << " or" << endl;
    14     gLog << "   showlog < filename" << endl << endl;
    15     gLog << " This program converts colored output made with ansi codes" << endl;
    16     gLog << " (like it is done by MLog) and redirected into a file back" << endl;
     17    gLog << "   showlog [options] < filename" << endl << endl;
     18    gLog << "   -?, -h, --help            This help" << endl << endl;
     19    gLog << endl;
     20    gLog << " This program converts colored output  made with ansi codes" << endl;
     21    gLog << " (like it is done by MLog)  and redirected into a file back" << endl;
    1722    gLog << " into colored output." << endl << endl;
     23    gLog << " It cannot be used to get rid of the Ansi codes. To display" << endl;
     24    gLog << " colored output with less use the option -R, eg." << endl;
     25    gLog << "   less -R logfile.log" << endl << endl;
    1826}
    1927
    … …  
    2230int main(int argc, char **argv)
    2331{
    24     if (argc>2)
     32    MArgs arg(argc, argv);
     33
     34    if (arg.HasOnly("-?") || arg.HasOnly("-h") || arg.HasOnly("--help"))
     35    {
     36        Usage();
     37        return -1;
     38    }
     39
     40    if (arg.GetNumOptions()>0)
     41    {
     42        gLog << warn << "WARNING - Unknown commandline options..." << endl;
     43        arg.Print("options");
     44        gLog << endl;
     45    }
     46
     47    //
     48    // check for the right usage of the program
     49    //
     50    if (arg.GetNumArguments()>1)
    2551    {
    2652        Usage();
    … …  
    2955
    3056    // casts necessary for gcc 2.95.3
    31     istream *in = argc==2 ? (istream*)new ifstream(argv[1]) : (istream*)&cin;
     57    istream *in = arg.GetNumArguments()==1 ? (istream*)new ifstream(arg.GetArgumentStr(0)) : (istream*)&cin;
    3258    if (!*in)
    3359    {
    34         gLog << "Cannot open file " << argv[1] << ": " << strerror(errno) << endl;
     60        gLog << "Cannot open file " << arg.GetArgumentStr(0) << ": " << strerror(errno) << endl;
    3561        return -1;
    3662    }

  • trunk/MagicSoft/Mars/star.cc

    r4459 r4756  
     1#include <TROOT.h>
     2#include <TClass.h>
    13#include <TSystem.h>
    2 
    3 #include "MParList.h"
    4 #include "MTaskList.h"
    5 #include "MEvtLoop.h"
    6 /*'
     4#include <TGClient.h>
     5#include <TApplication.h>
     6
    77#include "MLog.h"
    8  */
    98#include "MLogManip.h"
    109
    1110#include "MArgs.h"
    12 
    1311#include "MArray.h"
    14 #include "MReadMarsFile.h"
    15 #include "MGeomApply.h"
    16 #include "MMcPedestalCopy.h"
    17 #include "MMcPedestalNSBAdd.h"
    18 #include "MCerPhotCalc.h"
    19 #include "MCerPhotAnal2.h"
    20 #include "MSigmabarCalc.h"
    21 #include "MImgCleanStd.h"
    22 #include "MHillasCalc.h"
    23 #include "MHillasSrcCalc.h"
    24 #include "MWriteRootFile.h"
    25 #include "MFDataMember.h"
    26 #include "MFillH.h"
    27 #include "MReportDrive.h"    // PRELIMINARY:
    28 //#include "MPointingPos.h"    // PRELIMINARY:
    29 //#include "MPointingPosCalc.h"
     12#include "MDirIter.h"
     13
     14#include "MStatusDisplay.h"
     15
     16#include "MSequence.h"
     17#include "MJStar.h"
    3018
    3119using namespace std;
    32 
    33 //////////////////////////////////////////////////////////////////////////////
    34 //
    35 // This is an easy implementation of the Star process
    36 // (as compilable prog)
    37 //
    38 //////////////////////////////////////////////////////////////////////////////
    3920
    4021static void StartUpMessage()
    … …  
    4223    gLog << all << endl;
    4324
    44     //                1         2         3         4         5
    45     //       12345678901234567890123456789012345678901234567890
    46     gLog << "==================================================" << endl;
    47     gLog << "                 STAR - MARS V" << MARSVER          << endl;
    48     gLog << "   MARS - STandard Analysis and Reconstruction"     << endl;
    49     gLog << "            Compiled on <" << __DATE__ << ">"       << endl;
    50     gLog << "               Using ROOT v" << ROOTVER             << endl;
    51     gLog << "==================================================" << endl;
     25    //                1         2         3         4         5         6
     26    //       123456789012345678901234567890123456789012345678901234567890
     27    gLog << "========================================================" << endl;
     28    gLog << "                   Star - MARS V" << MARSVER              << endl;
     29    gLog << "      MARS -- STandard Analysis and Reconstruction"       << endl;
     30    gLog << "               Compiled on <" << __DATE__ << ">"          << endl;
     31    gLog << "                  Using ROOT v" << ROOTVER                << endl;
     32    gLog << "========================================================" << endl;
    5233    gLog << endl;
    5334}
    … …  
    5536static void Usage()
    5637{
     38    //                1         2         3         4         5         6         7         8
     39    //       12345678901234567890123456789012345678901234567890123456789012345678901234567890
    5740    gLog << all << endl;
    5841    gLog << "Sorry the usage is:" << endl;
    59     gLog << "   star [-a0] [-vn] [-cn] inputfile[.root] outputfile[.root]" << endl << endl;
    60     gLog << "     input file: Merpped or MC root file." << endl;
    61     gLog << "     ouput file: Star-file (image parameter file)" << endl;
    62     gLog << "     -a0: Do not use Ansii codes." << endl;
    63     gLog << "     -cn: Compression level n=1..9 [default=2]" << endl;
    64     gLog << "     -vn: Verbosity level n [default=2]" << endl;
    65     gLog << "     -u1: Update File (instead of Recreate)" << endl;
    66     gLog << "     -tn: Telescope Serial Number n [default=0]" << endl << endl;
    67     gLog << " -> Further setup is not possible at the moment, please use" << endl;
    68     gLog << "    the star.C root macro instead. Using an input card will" << endl;
    69     gLog << "    be implemented in the future." << endl << endl;
     42    gLog << " star [options] sequence.txt" << endl << endl;
     43    gLog << " Arguments:" << endl;
     44    gLog << "   sequence.txt:             An ascii file defining a sequence of runs" << endl;
     45    gLog << " Root Options:" << endl;
     46    gLog << "   -b                        Batch mode (no graphical output to screen)" << endl<<endl;
     47    gLog << " Options:" << endl;
     48    gLog.Usage();
     49    gLog << "   --debug-env               Debug setting resources from file" << endl << endl;
     50    gLog << endl;
     51    gLog << "   -f                        Force overwrite of existing files" << endl;
     52    gLog << "   -ff                       Force execution if not all files found" << endl;
     53    gLog << "   --ind=path                Path where to search for the calibrated data (Y)" << endl;
     54    gLog << "                             [default=standard path in datacenter]" << endl;
     55    gLog << "   --out=path                Path to write the all results to [def=local path]" << endl;
     56    gLog << "                             (overwrites --outc and --outy)" << endl;
     57    gLog << "   --print-seq               Print Sequence information" << endl;
     58    gLog << "   --print-files             Print Files taken from Sequence" << endl;
     59    gLog << "   --print-only              Do not excute anything except print" << endl;
     60    gLog << "   --config=star.rc          Resource file [default=star.rc]" << endl;
     61    gLog << endl;
     62    gLog << "   -?, -h, --help            This help" << endl << endl;
    7063}
    7164
    72 class MJStar : public MParContainer
    73 {
    74 private:
    75     TString fInputFile;
    76     TString fOutputFile;
    77 
    78     Bool_t fIsUpdate;
    79     Byte_t fComprLevel;
    80 
    81     Byte_t fSerialIdx;
    82 
    83     Bool_t CheckFiles()
    84     {
    85         if (fOutputFile.IsNull())
    86         {
    87             fOutputFile = fInputFile;
    88 
    89             if (fOutputFile.EndsWith(".raw"))
    90                 fOutputFile = fOutputFile(0, fOutputFile.Length()-4);
    91 
    92             if (fOutputFile.EndsWith(".rep"))
    93                 fOutputFile = fOutputFile(0, fOutputFile.Length()-4);
    94 
    95             if (fOutputFile.EndsWith(".txt"))
    96                 fOutputFile = fOutputFile(0, fOutputFile.Length()-4);
    97         }
    98 
    99         if (fInputFile.IsNull())
    100         {
    101             *fLog << err << "Sorry, no input file." << endl;
    102             return kFALSE;
    103         }
    104 
    105         if (fOutputFile.IsNull())
    106         {
    107             *fLog << err << "Sorry, no output file." << endl;
    108             return kFALSE;
    109         }
    110 
    111         if (!fInputFile.EndsWith(".root"))
    112             fInputFile += ".root";
    113 
    114         if (!fOutputFile.EndsWith(".root"))
    115             fOutputFile += ".root";
    116 
    117         //
    118         // check whether the given files are OK.
    119         //
    120         if (gSystem->AccessPathName(fInputFile, kFileExists))
    121         {
    122             *fLog << err << "Sorry, the input file '" << fInputFile << "' doesn't exist." << endl;
    123             return kFALSE;
    124         }
    125 
    126         if (!gSystem->AccessPathName(fOutputFile, kFileExists))
    127         {
    128             if (fIsUpdate)
    129                 gLog << warn << "Warning: File doesn't '" << fOutputFile << "' exist... recreating." << endl;
    130         }
    131 
    132         if (fIsUpdate || gSystem->AccessPathName(fOutputFile, kFileExists))
    133             if (!gSystem->AccessPathName(fOutputFile, kWritePermission))
    134         {
    135             gLog << err << "Sorry, you don't have write permission for '" << fOutputFile << "'." << endl;
    136             return kFALSE;
    137         }
    138         return kTRUE;
    139     }
    140 
    141 public:
    142     MJStar(const char *name=0, const char *title=0)
    143         : fIsUpdate(kFALSE), fComprLevel(1), fSerialIdx(0)
    144     {
    145         fName  = name;
    146         fTitle = title;
    147     }
    148 
    149     void SetInputFile(const char *f)  { fInputFile  = f; }
    150     void SetOutputFile(const char *f) { fOutputFile = f; }
    151     void SetSerialIdx(Byte_t i)       { fSerialIdx  = i; }
    152     void SetComprLevel(Byte_t l)      { fComprLevel = l; }
    153     void SetUpdate(Bool_t u=kTRUE)    { fIsUpdate   = u; }
    154 
    155     Bool_t Process()
    156     {
    157         if (!CheckFiles())
    158             return kFALSE;
    159 
    160         //
    161         // Create a empty Parameter List and an empty Task List
    162         // The tasklist is identified in the eventloop by its name
    163         //
    164         MParList plist;
    165 
    166         MTaskList tlist;
    167         plist.AddToList(&tlist);
    168 
    169         // PRELIMINARY:
    170         /*
    171          MReportDrive rep;
    172          plist.AddToList(&rep);
    173          MPointingPos pos;
    174          plist.AddToList(&pos);
    175          */
    176 
    177         //
    178         // Now setup the tasks and tasklist:
    179         // ---------------------------------
    180         //
    181         MReadMarsFile read("Events", fInputFile);
    182         read.DisableAutoScheme();
    183 
    184         MGeomApply        apply;
    185         MMcPedestalCopy   pcopy;
    186         MMcPedestalNSBAdd pnsb;
    187 
    188         //MPointingPosCalc  pcalc;
    189 
    190         MCerPhotCalc      ncalc;
    191         MCerPhotAnal2     nanal;
    192 
    193         MFDataMember f1("MRawRunHeader.fRunType", '>', 255.5);
    194         MFDataMember f2("MRawRunHeader.fRunType", '<', 255.5);
    195 
    196         ncalc.SetFilter(&f1);
    197         nanal.SetFilter(&f2);
    198 
    199         MSigmabarCalc     sgcal;
    200         //MFillH            fills("MHSigmaTheta",  "", "FillSigmaTheta");
    201         //MFillH            fillb("MHBlindPixels", "", "FillBlindPixels");
    202         MImgCleanStd      clean;
    203         MHillasCalc       hcalc;
    204         MHillasSrcCalc    scalc; // !!Preliminary!! Will be removed later!
    205         MWriteRootFile    write(fOutputFile, fIsUpdate?"UPDATE":"RECREATE", "Star output", fComprLevel);
    206 
    207         tlist.AddToList(&read);
    208         tlist.AddToList(&f1);
    209         tlist.AddToList(&f2);
    210         tlist.AddToList(&apply);
    211         tlist.AddToList(&pcopy);
    212         tlist.AddToList(&pnsb);
    213         //tlist.AddToList(&pcalc);
    214         tlist.AddToList(&ncalc);
    215         tlist.AddToList(&nanal);
    216         tlist.AddToList(&sgcal);
    217         //tlist.AddToList(&fills);
    218         //tlist.AddToList(&fillb);
    219         tlist.AddToList(&clean);
    220         tlist.AddToList(&hcalc);
    221         tlist.AddToList(&scalc);
    222         tlist.AddToList(&write);
    223 
    224         //
    225         // Set the serial number for all tasks in the current tasklist
    226         //
    227         tlist.SetSerialNumber(fSerialIdx);
    228 
    229         //
    230         // Setup tasks
    231         //
    232         write.AddContainer(write.AddSerialNumber("MMcEvt"),       "Events", kFALSE);
    233         write.AddContainer(write.AddSerialNumber("MSigmabar"),    "Events");
    234         write.AddContainer(write.AddSerialNumber("MHillas"),      "Events");
    235         write.AddContainer(write.AddSerialNumber("MHillasExt"),   "Events");
    236         write.AddContainer(write.AddSerialNumber("MHillasSrc"),   "Events");
    237         write.AddContainer(write.AddSerialNumber("MNewImagePar"), "Events");
    238         write.AddContainer(write.AddSerialNumber("MSrcPosCam"),   "RunHeaders");
    239         //write.AddContainer(write.AddSerialNumber("MHSigmaTheta"), "RunHeaders");
    240         if (!fIsUpdate)
    241         {
    242             write.AddContainer("MRawRunHeader", "RunHeaders");
    243             write.AddContainer("MMcRunHeader",  "RunHeaders", kFALSE);
    244         }
    245 
    246         //
    247         // Create and set up the eventloop
    248         //
    249         MEvtLoop evtloop;
    250         evtloop.SetParList(&plist);
    251 
    252         //
    253         // Execute your analysis
    254         //
    255         if (!evtloop.Eventloop())
    256         {
    257             gLog << err << "ERROR: Star eventloop failed!" << endl;
    258             return kFALSE;
    259         }
    260 
    261         tlist.PrintStatistics();
    262 
    263         //plist.FindObject("MHSigmaTheta")->Write();
    264 
    265         gLog << all << "Star finished successfull!" << endl;
    266 
    267         return kTRUE;
    268     }
    269 //    ClassDef(MJStar, 0)
    270 };
    271 
    272 //ClassImp(MJStar);
    273 
    274 
    275 int main(const int argc, char **argv)
     65int main(int argc, char **argv)
    27666{
    27767    StartUpMessage();
    … …  
    28272    MArgs arg(argc, argv);
    28373
    284     //
    285     // Set verbosity to highest level.
    286     //
    287     gLog.SetDebugLevel(arg.HasOption("-v") ? arg.GetIntAndRemove("-v") : 2);
    288 
    289     if (arg.HasOption("-a") && arg.GetIntAndRemove("-a")==0)
    290         gLog.SetNoColors();
    291 
    292     const int  kComprlvl = arg.HasOption("-c") ? arg.GetIntAndRemove("-c") : 1;
    293     const int  kTelIndex = arg.HasOption("-t") ? arg.GetIntAndRemove("-t") : 0;
    294     const bool kUpdate   = arg.HasOption("-u") && arg.GetIntAndRemove("-u")==1;
     74    if (arg.HasOnly("-?") || arg.HasOnly("-h") || arg.HasOnly("--help"))
     75    {
     76        Usage();
     77        return -1;
     78    }
     79
     80    gLog.Setup(arg);
     81
     82    const TString kConfig     = arg.GetStringAndRemove("--config=", "star.rc");
     83
     84    const Bool_t  kPrintSeq   = arg.HasOnlyAndRemove("--print-seq");
     85    const Bool_t  kPrintFiles = arg.HasOnlyAndRemove("--print-files");
     86    const Bool_t  kPrintOnly  = arg.HasOnlyAndRemove("--print-only");
     87    const Bool_t  kDebugEnv   = arg.HasOnlyAndRemove("--debug-env");
     88
     89    const Bool_t  kOverwrite  = arg.HasOnlyAndRemove("-f");
     90    const Bool_t  kForceExec  = arg.HasOnlyAndRemove("-ff");
     91
     92    const TString kInpath     = arg.GetStringAndRemove("--in=",  "");
     93    const TString kOutpath    = arg.GetStringAndRemove("--out=", ".");
     94
     95    if (arg.GetNumOptions()>0)
     96    {
     97        gLog << warn << "WARNING - Unknown commandline options..." << endl;
     98        arg.Print("options");
     99        gLog << endl;
     100    }
    295101
    296102    //
    297103    // check for the right usage of the program
    298104    //
    299     if (arg.GetNumArguments()!=2)
     105    if (arg.GetNumArguments()!=1)
    300106    {
    301107        Usage();
    … …  
    304110
    305111    //
    306     // Initialize Non-GUI (batch) mode
    307     //
    308     gROOT->SetBatch();
    309 
     112    // Setup sequence file and check for its existance
     113    //
     114    const TString kSequence = arg.GetArgumentStr(0);
     115
     116    if (gSystem->AccessPathName(kSequence, kFileExists))
     117    {
     118        gLog << err << "Sorry, sequence file '" << kSequence << "' doesn't exist." << endl;
     119        return -1;
     120    }
     121
     122    //
     123    // Setup sequence and check its validity
     124    //
     125    MSequence seq(kSequence);
     126    if (kPrintSeq)
     127    {
     128        gLog << all;
     129        gLog.Separator(kSequence);
     130        seq.Print();
     131        gLog << endl;
     132    }
     133    if (!seq.IsValid())
     134    {
     135        gLog << err << "Sequence read but not valid!" << endl << endl;
     136        return -1;
     137    }
     138
     139    //
     140    // Process print options
     141    //
     142    MDirIter iter;
     143    seq.SetupDatRuns(iter, kInpath);
     144
     145    const Int_t n0 = seq.SetupDatRuns(iter, kInpath);
     146    const Int_t n1 = seq.GetNumDatRuns();
     147
     148    if (kPrintFiles)
     149    {
     150        gLog << all;
     151        gLog.Separator("Data Files");
     152        iter.Print("all");
     153        gLog << endl;
     154    }
     155
     156    //
     157    // Check for existance of all files
     158    //
     159    if (n0 != n1)
     160    {
     161        if (kForceExec)
     162            gLog << warn << "WARNING";
     163        else
     164            gLog << err << "ERROR";
     165        gLog << " - " << n1 << " files in sequence defined, but " << n0 << " found in ";
     166        gLog << (kInpath.IsNull() ? "<defaultpath>" : kInpath.Data()) << endl;
     167        if (!kForceExec)
     168            return -1;
     169    }
     170
     171    if (kPrintOnly)
     172        return 0;
     173
     174    //
     175    // Initialize root
     176    //
    310177    MArray::Class()->IgnoreTObjectStreamer();
    311178    MParContainer::Class()->IgnoreTObjectStreamer();
    312179
    313     MJStar star;
    314     star.SetInputFile(arg.GetArgumentStr(0));
    315     star.SetOutputFile(arg.GetArgumentStr(1));
    316     star.SetComprLevel(kComprlvl);
    317     star.SetSerialIdx(kTelIndex);
    318     star.SetUpdate(kUpdate);
    319 
    320     if (!star.Process())
    321     {
    322         gLog << err << "Star failed!" << endl;
    323         return -1;
    324     }
    325 
    326     gLog << inf << "Star finished successfull!" << endl;
    327 
     180    TApplication app("Star", &argc, argv);
     181    if (gROOT->IsBatch() || !gClient)
     182    {
     183        gLog << err << "Bombing... maybe your DISPLAY variable is not set correctly!" << endl;
     184        return 1;
     185    }
     186
     187    //
     188    // Update frequency by default = 1Hz
     189    //
     190    MStatusDisplay *d = new MStatusDisplay;
     191
     192    // From now on each 'Exit' means: Terminate the application
     193    d->SetBit(MStatusDisplay::kExitLoopOnExit);
     194    d->SetTitle(kSequence);
     195
     196    //
     197    // Do calibration
     198    //
     199    MJStar job(Form("MJStar #%d", seq.GetSequence()));
     200    job.SetSequence(seq);
     201    job.SetEnv(kConfig);
     202    job.SetEnvDebug(kDebugEnv);
     203    job.SetDisplay(d);;
     204    job.SetOverwrite(kOverwrite);
     205    job.SetPathOut(kOutpath);
     206    job.SetPathData(kInpath);
     207    // job.SetPathIn(kInpath); // not yet needed
     208
     209    if (!job.ProcessFile())
     210    {
     211        gLog << err << "Calculation of image parameters failed." << endl << endl;
     212        return -1;
     213    }
     214
     215    if (!job.GetDisplay())
     216    {
     217        gLog << warn << "Display closed by user... execution aborted." << endl << endl;
     218        return 1;
     219    }
     220
     221    // From now on each 'Close' means: Terminate the application
     222    d->SetBit(MStatusDisplay::kExitLoopOnClose);
     223
     224    // Wait until the user decides to exit the application
     225    app.Run(kFALSE);
    328226    return 0;
    329 /*
    330     //
    331     // This is to make argv[i] more readable inside the code
    332     //
    333     TString kNamein  = arg.GetArgumentStr(0);
    334     TString kNameout = arg.GetArgumentStr(1);
    335 
    336     if (!kNamein.EndsWith(".root"))
    337         kNamein += ".root";
    338 
    339     if (!kNameout.EndsWith(".root"))
    340         kNameout += ".root";
    341 
    342     //
    343     // check whether the given files are OK.
    344     //
    345     if (gSystem->AccessPathName(kNamein, kFileExists))
    346     {
    347         gLog << err << "Sorry, the input file '" << kNamein << "' doesn't exist." << endl;
    348         return -1;
    349     }
    350 
    351     if (!gSystem->AccessPathName(kNameout, kFileExists))
    352     {
    353         if (kUpdate)
    354             gLog << warn << "Warning: File doesn't '" << kNameout << "' exist... recreating." << endl;
    355     }
    356 
    357     if (kUpdate || gSystem->AccessPathName(kNameout, kFileExists))
    358         if (!gSystem->AccessPathName(kNameout, kWritePermission))
    359         {
    360             gLog << err << "Sorry, you don't have write permission for '" << kNameout << "'." << endl;
    361             return -1;
    362         }
    363 
    364     //
    365     // Create a empty Parameter List and an empty Task List
    366     // The tasklist is identified in the eventloop by its name
    367     //
    368     MParList plist;
    369 
    370     MTaskList tlist;
    371     plist.AddToList(&tlist);
    372 
    373     // PRELIMINARY:
    374     // MReportDrive rep;
    375     // plist.AddToList(&rep);
    376     // MPointingPos pos;
    377     // plist.AddToList(&pos);
    378 
    379     //
    380     // Now setup the tasks and tasklist:
    381     // ---------------------------------
    382     //
    383     MReadMarsFile read("Events", kNamein);
    384     read.DisableAutoScheme();
    385 
    386     MGeomApply        apply;
    387     MMcPedestalCopy   pcopy;
    388     MMcPedestalNSBAdd pnsb;
    389 
    390     //MPointingPosCalc  pcalc;
    391 
    392     MCerPhotCalc      ncalc;
    393     MCerPhotAnal2     nanal;
    394 
    395     MFDataMember f1("MRawRunHeader.fRunType", '>', 255.5);
    396     MFDataMember f2("MRawRunHeader.fRunType", '<', 255.5);
    397 
    398     ncalc.SetFilter(&f1);
    399     nanal.SetFilter(&f2);
    400 
    401     MSigmabarCalc     sgcal;
    402     //MFillH            fills("MHSigmaTheta",  "", "FillSigmaTheta");
    403     MImgCleanStd      clean;
    404     MHillasCalc       hcalc;
    405     MHillasSrcCalc    scalc; // !!Preliminary!! Will be removed later!
    406     MWriteRootFile    write(kNameout, kUpdate?"UPDATE":"RECREATE", "Star output", kComprlvl);
    407 
    408     tlist.AddToList(&read);
    409     tlist.AddToList(&f1);
    410     tlist.AddToList(&f2);
    411     tlist.AddToList(&apply);
    412     tlist.AddToList(&pcopy);
    413     tlist.AddToList(&pnsb);
    414     //tlist.AddToList(&pcalc);
    415     tlist.AddToList(&ncalc);
    416     tlist.AddToList(&nanal);
    417     tlist.AddToList(&sgcal);
    418     //tlist.AddToList(&fills);
    419     //tlist.AddToList(&fillb);
    420     tlist.AddToList(&clean);
    421     tlist.AddToList(&hcalc);
    422     tlist.AddToList(&scalc);
    423     tlist.AddToList(&write);
    424 
    425     //
    426     // Set the serial number for all tasks in the current tasklist
    427     //
    428     tlist.SetSerialNumber(kTelIndex);
    429 
    430     //
    431     // Setup tasks
    432     //
    433 
    434     write.AddContainer(write.AddSerialNumber("MMcEvt"),       "Events", kFALSE);
    435     write.AddContainer(write.AddSerialNumber("MSigmabar"),    "Events");
    436     write.AddContainer(write.AddSerialNumber("MHillas"),      "Events");
    437     write.AddContainer(write.AddSerialNumber("MHillasExt"),   "Events");
    438     write.AddContainer(write.AddSerialNumber("MHillasSrc"),   "Events");
    439     write.AddContainer(write.AddSerialNumber("MNewImagePar"), "Events");
    440     write.AddContainer(write.AddSerialNumber("MSrcPosCam"),   "RunHeaders");
    441     //write.AddContainer(write.AddSerialNumber("MHSigmaTheta"), "RunHeaders");
    442     if (!kUpdate)
    443     {
    444         write.AddContainer("MRawRunHeader", "RunHeaders");
    445         write.AddContainer("MMcRunHeader",  "RunHeaders", kFALSE);
    446     }
    447 
    448     //
    449     // Create and set up the eventloop
    450     //
    451     MEvtLoop evtloop;
    452     evtloop.SetParList(&plist);
    453 
    454     //
    455     // Execute your analysis
    456     //
    457     if (!evtloop.Eventloop())
    458     {
    459         gLog << err << "ERROR: Star eventloop failed!" << endl;
    460         return -1;
    461     }
    462 
    463     tlist.PrintStatistics();
    464 
    465     //plist.FindObject("MHSigmaTheta")->Write();
    466 
    467     gLog << all << "Star finished successfull!" << endl;
    468     return 0;
    469     */
    470227}
Note: See TracChangeset for help on using the changeset viewer.