| 1 | /* ======================================================================== *\
|
|---|
| 2 | !
|
|---|
| 3 | ! *
|
|---|
| 4 | ! * This file is part of MARS, the MAGIC Analysis and Reconstruction
|
|---|
| 5 | ! * Software. It is distributed to you in the hope that it can be a useful
|
|---|
| 6 | ! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
|
|---|
| 7 | ! * It is distributed WITHOUT ANY WARRANTY.
|
|---|
| 8 | ! *
|
|---|
| 9 | ! * Permission to use, copy, modify and distribute this software and its
|
|---|
| 10 | ! * documentation for any purpose is hereby granted without fee,
|
|---|
| 11 | ! * provided that the above copyright notice appear in all copies and
|
|---|
| 12 | ! * that both that copyright notice and this permission notice appear
|
|---|
| 13 | ! * in supporting documentation. It is provided "as is" without express
|
|---|
| 14 | ! * or implied warranty.
|
|---|
| 15 | ! *
|
|---|
| 16 | !
|
|---|
| 17 | !
|
|---|
| 18 | ! Author(s): Thomas Bretz, 8/2004 <mailto:tbretz@astro.uni-wuerzburg.de>
|
|---|
| 19 | !
|
|---|
| 20 | ! Copyright: MAGIC Software Development, 2004-2007
|
|---|
| 21 | !
|
|---|
| 22 | !
|
|---|
| 23 | \* ======================================================================== */
|
|---|
| 24 |
|
|---|
| 25 | /////////////////////////////////////////////////////////////////////////////
|
|---|
| 26 | //
|
|---|
| 27 | // MSequence
|
|---|
| 28 | //
|
|---|
| 29 | // This class describes a sequence. For sequences see:
|
|---|
| 30 | // http://db.astro.uni-wuerzburg.de
|
|---|
| 31 | //
|
|---|
| 32 | // A sequence is a collection of runs which should be used together.
|
|---|
| 33 | // Any run can be contained only once.
|
|---|
| 34 | //
|
|---|
| 35 | // Here is an example how a file describing a sequence could look like.
|
|---|
| 36 | //
|
|---|
| 37 | // Reading the file is based on TEnv. For more details see also
|
|---|
| 38 | // the class reference of TEnv.
|
|---|
| 39 | //
|
|---|
| 40 | // ===========================================================================
|
|---|
| 41 | //
|
|---|
| 42 | // sequence.txt
|
|---|
| 43 | // ------------
|
|---|
| 44 | //
|
|---|
| 45 | // # Sequence number (identifier) - necessary if the path
|
|---|
| 46 | // # contains the sequence number, e.g. image files
|
|---|
| 47 | // Sequence: 31015
|
|---|
| 48 | //
|
|---|
| 49 | // # Observation Period (not needed)
|
|---|
| 50 | // Period: 18
|
|---|
| 51 | //
|
|---|
| 52 | // # Date of sunrise of the observation night - necessary if the path
|
|---|
| 53 | // # contains the date, e.g. raw data
|
|---|
| 54 | // Night: 2004-06-24
|
|---|
| 55 | //
|
|---|
| 56 | // # Start time of the sequence (first data run, not needed)
|
|---|
| 57 | // Start: 2004-06-24 03:12:42
|
|---|
| 58 | //
|
|---|
| 59 | // # Run number of last data run in sequence (not needed)
|
|---|
| 60 | // LastRun: 31032
|
|---|
| 61 | //
|
|---|
| 62 | // # Project name of data-runs of sequence (not needed)
|
|---|
| 63 | // Project: 3EG2033+41
|
|---|
| 64 | //
|
|---|
| 65 | // # Source name of all runs of sequence (not needed)
|
|---|
| 66 | // Source: 3EG2033+41
|
|---|
| 67 | //
|
|---|
| 68 | // # Trigger table of data-runs of sequence (not needed)
|
|---|
| 69 | // TriggerTable: L1_4NN:L2_DEFAULT
|
|---|
| 70 | //
|
|---|
| 71 | // # HV Setting table of data-runs of sequence (not needed)
|
|---|
| 72 | // HvSettings: HVSettings_FF36q
|
|---|
| 73 | //
|
|---|
| 74 | // # Total number of data-events in sequence (not needed)
|
|---|
| 75 | // NumEvents: 250914
|
|---|
| 76 | //
|
|---|
| 77 | // # Whether this is MC data or not (necessary in case of MCs if
|
|---|
| 78 | // # default paths should be used)
|
|---|
| 79 | // MonteCarlo: Yes
|
|---|
| 80 | //
|
|---|
| 81 | // # List of all runs of this sequence (not needed)
|
|---|
| 82 | // Runs: 31015 31016 31017 31018 31019 31020 31021 31022 31023 31024 31025 31026 31027 31028 31029 31030 31031 31032
|
|---|
| 83 | //
|
|---|
| 84 | // # List of all calibration runs of this sequence (necessary if accessed)
|
|---|
| 85 | // CalRuns: 31015 31016 31017
|
|---|
| 86 | // # List of pedestal runs belonging to the calibration runs of this sequence (necessary if accessed)
|
|---|
| 87 | // PedRuns: 31018
|
|---|
| 88 | // # List of all data runs belonging to this sequence (necessary)
|
|---|
| 89 | // DatRuns: 31019 31020 31022 31023 31024 31025 31027 31028 31030 31032
|
|---|
| 90 | //
|
|---|
| 91 | // # Just for fun ;-) (not needed, but helpful)
|
|---|
| 92 | // Comment: This is a template for a sequence file
|
|---|
| 93 | //
|
|---|
| 94 | // ===========================================================================
|
|---|
| 95 | //
|
|---|
| 96 | // For special cases you can also setup a sequence directly from a macro,
|
|---|
| 97 | // for example:
|
|---|
| 98 | //
|
|---|
| 99 | // MDirIter pediter, datiter, caliter;
|
|---|
| 100 | //
|
|---|
| 101 | // MSequence seq;
|
|---|
| 102 | // seq.SetNight("2004-07-06");
|
|---|
| 103 | // seq.AddPedRuns(31751);
|
|---|
| 104 | // seq.AddCalRuns(31752);
|
|---|
| 105 | // seq.AddDatRuns(31753, 31764);
|
|---|
| 106 | // seq.SetupPedRuns(pediter);
|
|---|
| 107 | // seq.SetupCalRuns(caliter);
|
|---|
| 108 | // seq.SetupDatRuns(datiter);
|
|---|
| 109 | //
|
|---|
| 110 | // or
|
|---|
| 111 | //
|
|---|
| 112 | // MDirIter iter;
|
|---|
| 113 | //
|
|---|
| 114 | // MSequence seq;
|
|---|
| 115 | // seq.SetNight("2004-07-06");
|
|---|
| 116 | // seq.AddRuns(31753, 31764);
|
|---|
| 117 | // seq.SetupRuns(iter);
|
|---|
| 118 | // seq.SetupPedRuns(iter, "/mypath", "[DPC]");
|
|---|
| 119 | //
|
|---|
| 120 | // ===========================================================================
|
|---|
| 121 | //
|
|---|
| 122 | // Class Version 2:
|
|---|
| 123 | // + fMonteCarlo
|
|---|
| 124 | //
|
|---|
| 125 | // Class Version 3:
|
|---|
| 126 | // + fComment
|
|---|
| 127 | //
|
|---|
| 128 | // Class Version 4:
|
|---|
| 129 | // + fExclRuns
|
|---|
| 130 | //
|
|---|
| 131 | /////////////////////////////////////////////////////////////////////////////
|
|---|
| 132 | #include "MSequence.h"
|
|---|
| 133 |
|
|---|
| 134 | #include <stdlib.h>
|
|---|
| 135 |
|
|---|
| 136 | #include <TEnv.h>
|
|---|
| 137 | #include <TRegexp.h>
|
|---|
| 138 | #include <TSystem.h> // TSystem::ExpandPath
|
|---|
| 139 |
|
|---|
| 140 | #include "MLog.h"
|
|---|
| 141 | #include "MLogManip.h"
|
|---|
| 142 |
|
|---|
| 143 | #include "MEnv.h"
|
|---|
| 144 | #include "MJob.h"
|
|---|
| 145 | #include "MAstro.h"
|
|---|
| 146 | #include "MString.h"
|
|---|
| 147 | #include "MDirIter.h"
|
|---|
| 148 |
|
|---|
| 149 | ClassImp(MSequence);
|
|---|
| 150 |
|
|---|
| 151 | using namespace std;
|
|---|
| 152 |
|
|---|
| 153 | MSequence::~MSequence()
|
|---|
| 154 | {
|
|---|
| 155 | /*
|
|---|
| 156 | TExMapIter iter(&fFileNames);
|
|---|
| 157 |
|
|---|
| 158 | Long_t key, val;
|
|---|
| 159 |
|
|---|
| 160 | while (iter.Next(key, val))
|
|---|
| 161 | delete (TString*)val;
|
|---|
| 162 | */
|
|---|
| 163 | }
|
|---|
| 164 |
|
|---|
| 165 |
|
|---|
| 166 | // --------------------------------------------------------------------------
|
|---|
| 167 | //
|
|---|
| 168 | // Copy the run numbers from the TString runs into the TArrayI data.
|
|---|
| 169 | // Runs which are twice in the list are only added once. In this case
|
|---|
| 170 | // a warning is emitted.
|
|---|
| 171 | //
|
|---|
| 172 | void MSequence::Split(TString &runs, TArrayI &data) const
|
|---|
| 173 | {
|
|---|
| 174 | const TRegexp regexp("[0-9]+");
|
|---|
| 175 |
|
|---|
| 176 | data.Set(0);
|
|---|
| 177 |
|
|---|
| 178 | runs.ReplaceAll("\t", " ");
|
|---|
| 179 | runs = runs.Strip(TString::kBoth);
|
|---|
| 180 |
|
|---|
| 181 | while (!runs.IsNull())
|
|---|
| 182 | {
|
|---|
| 183 | const TString num = runs(regexp);
|
|---|
| 184 |
|
|---|
| 185 | if (num.IsNull())
|
|---|
| 186 | {
|
|---|
| 187 | *fLog << warn << "WARNING - Run is NaN (not a number): '" << runs << "'" << endl;
|
|---|
| 188 | break;
|
|---|
| 189 | }
|
|---|
| 190 |
|
|---|
| 191 | const Int_t run = atoi(num.Data());
|
|---|
| 192 | const Int_t n = data.GetSize();
|
|---|
| 193 |
|
|---|
| 194 | // skip already existing entries
|
|---|
| 195 | int i;
|
|---|
| 196 | for (i=0; i<n; i++)
|
|---|
| 197 | if (data[i] == run)
|
|---|
| 198 | break;
|
|---|
| 199 |
|
|---|
| 200 | if (i<n)
|
|---|
| 201 | *fLog << warn << "WARNING - Run #" << run << " already in list... skipped." << endl;
|
|---|
| 202 | else
|
|---|
| 203 | {
|
|---|
| 204 | // set new entry
|
|---|
| 205 | data.Set(n+1);
|
|---|
| 206 | data[n] = run;
|
|---|
| 207 | }
|
|---|
| 208 |
|
|---|
| 209 | runs.Remove(0, runs.First(num)+num.Length());
|
|---|
| 210 | }
|
|---|
| 211 |
|
|---|
| 212 | MJob::SortArray(data);
|
|---|
| 213 | }
|
|---|
| 214 |
|
|---|
| 215 | UInt_t MSequence::SetupRuns(MDirIter &iter, const TArrayI &arr, FileType_t type, const char *path) const
|
|---|
| 216 | {
|
|---|
| 217 | TString d(path);
|
|---|
| 218 | if (d.IsNull())
|
|---|
| 219 | d = fDataPath;
|
|---|
| 220 |
|
|---|
| 221 | const Bool_t def = d.IsNull();
|
|---|
| 222 |
|
|---|
| 223 | // For this particular case we assume that the files are added one by
|
|---|
| 224 | // one without wildcards.
|
|---|
| 225 | const Int_t n0 = iter.GetNumEntries();
|
|---|
| 226 |
|
|---|
| 227 | // Setup path
|
|---|
| 228 | if (def)
|
|---|
| 229 | {
|
|---|
| 230 | d = GetStandardPath();
|
|---|
| 231 | switch (type)
|
|---|
| 232 | {
|
|---|
| 233 | case kRawDat: // rawdata
|
|---|
| 234 | case kRawPed:
|
|---|
| 235 | case kRawCal:
|
|---|
| 236 | case kRawAll:
|
|---|
| 237 | case kRootDat: // mcdata
|
|---|
| 238 | case kRootPed:
|
|---|
| 239 | case kRootCal:
|
|---|
| 240 | case kRootAll:
|
|---|
| 241 | d += "rawfiles/";
|
|---|
| 242 | d += fNight.GetStringFmt("%Y/%m/%d");
|
|---|
| 243 | break;
|
|---|
| 244 | case kCalibrated:
|
|---|
| 245 | d += Form("callisto/%04d/%08d", fSequence/10000, fSequence);
|
|---|
| 246 | break;
|
|---|
| 247 | case kImages:
|
|---|
| 248 | d += Form("star/%04d/%08d", fSequence/10000, fSequence);
|
|---|
| 249 | break;
|
|---|
| 250 | }
|
|---|
| 251 | }
|
|---|
| 252 | else
|
|---|
| 253 | gSystem->ExpandPathName(d);
|
|---|
| 254 |
|
|---|
| 255 | if (!d.EndsWith("/"))
|
|---|
| 256 | d += '/';
|
|---|
| 257 |
|
|---|
| 258 | Int_t excluded = 0;
|
|---|
| 259 | for (int i=0; i<arr.GetSize(); i++)
|
|---|
| 260 | {
|
|---|
| 261 | if (IsExcluded(arr[i]))
|
|---|
| 262 | {
|
|---|
| 263 | excluded++;
|
|---|
| 264 | continue;
|
|---|
| 265 | }
|
|---|
| 266 |
|
|---|
| 267 | // R. DeLosReyes and T. Bretz
|
|---|
| 268 | // Changes to read the DAQ numbering format. Changes takes place
|
|---|
| 269 | // between runs 35487 and 00035488 (2004_08_30)
|
|---|
| 270 | const char *fmt = arr[i]>35487 || fMonteCarlo ? "%08d_%s_*_E" : "%05d_%s_*_E";
|
|---|
| 271 |
|
|---|
| 272 | TString n;
|
|---|
| 273 | const char *id="_";
|
|---|
| 274 | switch (type)
|
|---|
| 275 | {
|
|---|
| 276 | case kRawDat:
|
|---|
| 277 | case kRootDat:
|
|---|
| 278 | id = "D";
|
|---|
| 279 | break;
|
|---|
| 280 | case kRawPed:
|
|---|
| 281 | case kRootPed:
|
|---|
| 282 | id = "P";
|
|---|
| 283 | break;
|
|---|
| 284 | case kRawCal:
|
|---|
| 285 | case kRootCal:
|
|---|
| 286 | id = "C";
|
|---|
| 287 | break;
|
|---|
| 288 | case kRawAll:
|
|---|
| 289 | case kRootAll:
|
|---|
| 290 | id = "[PCD]";
|
|---|
| 291 | break;
|
|---|
| 292 | case kCalibrated:
|
|---|
| 293 | id = "Y";
|
|---|
| 294 | break;
|
|---|
| 295 | case kImages:
|
|---|
| 296 | id = "I";
|
|---|
| 297 | break;
|
|---|
| 298 | }
|
|---|
| 299 |
|
|---|
| 300 | // Create file name
|
|---|
| 301 | n = fNight.GetStringFmt("%Y%m%d_");
|
|---|
| 302 | n += Form(fmt, arr[i], id);
|
|---|
| 303 |
|
|---|
| 304 | switch (type)
|
|---|
| 305 | {
|
|---|
| 306 | case kRawDat:
|
|---|
| 307 | case kRawPed:
|
|---|
| 308 | case kRawCal:
|
|---|
| 309 | case kRawAll:
|
|---|
| 310 | n += ".raw.?g?z?";
|
|---|
| 311 | break;
|
|---|
| 312 | default:
|
|---|
| 313 | n += ".root";
|
|---|
| 314 | }
|
|---|
| 315 |
|
|---|
| 316 | // Check existance and accessibility of file
|
|---|
| 317 | MDirIter file(d, n, 0);
|
|---|
| 318 | TString name = file();
|
|---|
| 319 | gSystem->ExpandPathName(name);
|
|---|
| 320 | if (gSystem->AccessPathName(name, kFileExists))
|
|---|
| 321 | {
|
|---|
| 322 | *fLog << err;
|
|---|
| 323 | *fLog << "ERROR - File " << d << n << " not accessible!" << endl;
|
|---|
| 324 | return 0;
|
|---|
| 325 | }
|
|---|
| 326 | if (!file().IsNull())
|
|---|
| 327 | {
|
|---|
| 328 | *fLog << err;
|
|---|
| 329 | *fLog << "ERROR - Searching for file " << d << n << " gave more than one result!" << endl;
|
|---|
| 330 | return 0;
|
|---|
| 331 | }
|
|---|
| 332 |
|
|---|
| 333 | // Add Path/File to TIter
|
|---|
| 334 | iter.AddDirectory(d, n, 0);
|
|---|
| 335 | }
|
|---|
| 336 |
|
|---|
| 337 | // n0: The previous contents of the iter
|
|---|
| 338 | // n1: The number of files which have been added to the iter
|
|---|
| 339 | // n2: The number of files which should have been added from the array
|
|---|
| 340 | const Int_t n1 = iter.GetNumEntries()-n0;
|
|---|
| 341 | const Int_t n2 = arr.GetSize()-excluded;
|
|---|
| 342 | if (n1==0)
|
|---|
| 343 | {
|
|---|
| 344 | *fLog << err;
|
|---|
| 345 | *fLog << "ERROR - No input files for sequence #" << GetSequence() << endl;
|
|---|
| 346 | *fLog << " read from " << GetName() << endl;
|
|---|
| 347 | *fLog << " found in" << (def?" default-path ":" ") << d << endl;
|
|---|
| 348 | return 0;
|
|---|
| 349 | }
|
|---|
| 350 |
|
|---|
| 351 | if (n1==n2)
|
|---|
| 352 | return n1;
|
|---|
| 353 |
|
|---|
| 354 | *fLog << err;
|
|---|
| 355 | *fLog << "ERROR - " << n1 << " input files for sequence #" << GetSequence() << " found in" << endl;
|
|---|
| 356 | *fLog << " " << (def?" default-path ":" ") << d << endl;
|
|---|
| 357 | *fLog << " but " << n2 << " files were defined in sequence file" << endl;
|
|---|
| 358 | *fLog << " " << GetName() << endl;
|
|---|
| 359 | if (fLog->GetDebugLevel()<=4)
|
|---|
| 360 | return 0;
|
|---|
| 361 |
|
|---|
| 362 | *fLog << inf << "Files which are searched for this sequence:" << endl;
|
|---|
| 363 | iter.Print();
|
|---|
| 364 | return 0;
|
|---|
| 365 | }
|
|---|
| 366 |
|
|---|
| 367 | // --------------------------------------------------------------------------
|
|---|
| 368 | //
|
|---|
| 369 | // Read the file fname as setup file for the sequence.
|
|---|
| 370 | //
|
|---|
| 371 | //void MSequence::GetFileNames(TEnv &env, const TArrayI &arr)
|
|---|
| 372 | //{
|
|---|
| 373 | /*
|
|---|
| 374 | for (int i=0; i<arr.GetSize(); i++)
|
|---|
| 375 | {
|
|---|
| 376 | // Get run number
|
|---|
| 377 | const Int_t num = arr[i];
|
|---|
| 378 |
|
|---|
| 379 | // Check if name already set
|
|---|
| 380 | if (fFileNames.GetValue(num))
|
|---|
| 381 | continue;
|
|---|
| 382 |
|
|---|
| 383 | TString *str = new TString(env.GetValue(Form("%d", num), ""));
|
|---|
| 384 | fFileNames.Add(num, (Long_t)str);
|
|---|
| 385 | }
|
|---|
| 386 | */
|
|---|
| 387 | //}
|
|---|
| 388 |
|
|---|
| 389 | // --------------------------------------------------------------------------
|
|---|
| 390 | //
|
|---|
| 391 | // Get a file name corresponding to the run-number num, returns 0 if n/a
|
|---|
| 392 | //
|
|---|
| 393 | //const char *MSequence::GetFileName(UInt_t num)
|
|---|
| 394 | //{
|
|---|
| 395 | // return 0;
|
|---|
| 396 | /*
|
|---|
| 397 | TString *str = (TString*)fFileNames.GetValue(num);
|
|---|
| 398 | return str ? str->Data() : 0;*/
|
|---|
| 399 | //}
|
|---|
| 400 |
|
|---|
| 401 | MSequence::LightCondition_t MSequence::ReadLightCondition(TEnv &env, const char *prefix) const
|
|---|
| 402 | {
|
|---|
| 403 | TString str = GetEnvValue2(env, prefix, "LightConditions", "n/a");
|
|---|
| 404 | if (!str.CompareTo("n/a", TString::kIgnoreCase))
|
|---|
| 405 | return kNA;
|
|---|
| 406 | if (!str.CompareTo("No_Moon", TString::kIgnoreCase))
|
|---|
| 407 | return kNoMoon;
|
|---|
| 408 | if (!str.CompareTo("NoMoon", TString::kIgnoreCase))
|
|---|
| 409 | return kNoMoon;
|
|---|
| 410 | if (!str.CompareTo("Twilight", TString::kIgnoreCase))
|
|---|
| 411 | return kTwilight;
|
|---|
| 412 | if (!str.CompareTo("Moon", TString::kIgnoreCase))
|
|---|
| 413 | return kMoon;
|
|---|
| 414 | if (!str.CompareTo("Day", TString::kIgnoreCase))
|
|---|
| 415 | return kDay;
|
|---|
| 416 |
|
|---|
| 417 | gLog << warn;
|
|---|
| 418 | gLog << "WARNING - in " << fFileName << ":" << endl;
|
|---|
| 419 | gLog << " LightCondition-tag is '" << str << "' but must be n/a, no_moon, twilight, moon or day." << endl;
|
|---|
| 420 | return kNA;
|
|---|
| 421 | }
|
|---|
| 422 |
|
|---|
| 423 | // --------------------------------------------------------------------------
|
|---|
| 424 | //
|
|---|
| 425 | // Read the file fname as setup file for the sequence.
|
|---|
| 426 | //
|
|---|
| 427 | MSequence::MSequence(const char *fname, const char *path, UInt_t seq)
|
|---|
| 428 | {
|
|---|
| 429 | fName = fname;
|
|---|
| 430 | fTitle = path;
|
|---|
| 431 |
|
|---|
| 432 | fFileName = fname;
|
|---|
| 433 | fDataPath = path;
|
|---|
| 434 |
|
|---|
| 435 | gSystem->ExpandPathName(fName);
|
|---|
| 436 | gSystem->ExpandPathName(fTitle);
|
|---|
| 437 |
|
|---|
| 438 | const Bool_t rc1 = gSystem->AccessPathName(fName, kFileExists);
|
|---|
| 439 | const Bool_t rc2 = !fTitle.IsNull() && gSystem->AccessPathName(fTitle, kFileExists);
|
|---|
| 440 |
|
|---|
| 441 | if (rc1)
|
|---|
| 442 | gLog << err << "ERROR - Sequence file '" << fName << "' doesn't exist." << endl;
|
|---|
| 443 | if (rc2)
|
|---|
| 444 | gLog << err << "ERROR - Directory '" << fTitle << "' doesn't exist." << endl;
|
|---|
| 445 |
|
|---|
| 446 | MEnv env(fName);
|
|---|
| 447 |
|
|---|
| 448 | fSequence = (UInt_t)env.GetValue("Sequence", (Int_t)seq);
|
|---|
| 449 | if (rc1 || rc2)
|
|---|
| 450 | fSequence = (UInt_t)-1;
|
|---|
| 451 |
|
|---|
| 452 | const TString prefix = Form("Sequence%08d", fSequence);
|
|---|
| 453 |
|
|---|
| 454 | fLastRun = GetEnvValue2(env, prefix, "LastRun", -1);
|
|---|
| 455 | fNumEvents = GetEnvValue2(env, prefix, "NumEvents", -1);
|
|---|
| 456 | fPeriod = GetEnvValue2(env, prefix, "Period", -1);
|
|---|
| 457 |
|
|---|
| 458 | fLightCondition = ReadLightCondition(env, prefix);
|
|---|
| 459 |
|
|---|
| 460 | TString str;
|
|---|
| 461 | str = GetEnvValue2(env, prefix, "Start", "");
|
|---|
| 462 | fStart.SetSqlDateTime(str);
|
|---|
| 463 | str = GetEnvValue2(env, prefix, "Night", "");
|
|---|
| 464 | str += " 00:00:00";
|
|---|
| 465 | fNight.SetSqlDateTime(str);
|
|---|
| 466 |
|
|---|
| 467 | fProject = GetEnvValue2(env, prefix, "Project", "");
|
|---|
| 468 | fSource = GetEnvValue2(env, prefix, "Source", "");
|
|---|
| 469 | fTriggerTable = GetEnvValue2(env, prefix, "TriggerTable", "");
|
|---|
| 470 | fHvSettings = GetEnvValue2(env, prefix, "HvSettings", "");
|
|---|
| 471 | fMonteCarlo = GetEnvValue2(env, prefix, "MonteCarlo", kFALSE);
|
|---|
| 472 | fComment = GetEnvValue2(env, prefix, "Comment", "");
|
|---|
| 473 |
|
|---|
| 474 | str = GetEnvValue2(env, prefix, "Runs", "");
|
|---|
| 475 | Split(str, fRuns);
|
|---|
| 476 | str = GetEnvValue2(env, prefix, "CalRuns", "");
|
|---|
| 477 | Split(str, fCalRuns);
|
|---|
| 478 | str = GetEnvValue2(env, prefix, "PedRuns", "");
|
|---|
| 479 | Split(str, fPedRuns);
|
|---|
| 480 | str = GetEnvValue2(env, prefix, "DatRuns", "");
|
|---|
| 481 | Split(str, fDatRuns);
|
|---|
| 482 | str = GetEnvValue2(env, prefix, "Exclude", "");
|
|---|
| 483 | Split(str, fExclRuns);
|
|---|
| 484 |
|
|---|
| 485 | // GetFileNames(env, fRuns);
|
|---|
| 486 | // GetFileNames(env, fCalRuns);
|
|---|
| 487 | // GetFileNames(env, fPedRuns);
|
|---|
| 488 | // GetFileNames(env, fDatRuns);
|
|---|
| 489 |
|
|---|
| 490 | // Dummies:
|
|---|
| 491 | env.Touch("ZdMin");
|
|---|
| 492 | env.Touch("ZdMax");
|
|---|
| 493 | env.Touch("L1TriggerTable");
|
|---|
| 494 | env.Touch("L2TriggerTable");
|
|---|
| 495 |
|
|---|
| 496 | if (seq<0 && env.GetNumUntouched()>0)
|
|---|
| 497 | {
|
|---|
| 498 | gLog << warn << "WARNING - At least one resource in the sequence-file has not been touched!" << endl;
|
|---|
| 499 | env.PrintUntouched();
|
|---|
| 500 | }
|
|---|
| 501 | }
|
|---|
| 502 |
|
|---|
| 503 | //---------------------------------------------------------------------------
|
|---|
| 504 | //
|
|---|
| 505 | // Make sure that the name used for writing doesn't contain a full path
|
|---|
| 506 | //
|
|---|
| 507 | const char *MSequence::GetName() const
|
|---|
| 508 | {
|
|---|
| 509 | const char *pos = strrchr(GetRcName(), '/');
|
|---|
| 510 | return pos>0 ? pos+1 : GetRcName();
|
|---|
| 511 | }
|
|---|
| 512 |
|
|---|
| 513 | // --------------------------------------------------------------------------
|
|---|
| 514 | //
|
|---|
| 515 | // Print the contents of the sequence
|
|---|
| 516 | //
|
|---|
| 517 | void MSequence::Print(Option_t *o) const
|
|---|
| 518 | {
|
|---|
| 519 | const TString pre = TString(o).Contains("prefixed") ? Form("Sequence%08d.", fSequence) : "";
|
|---|
| 520 |
|
|---|
| 521 | gLog << all;
|
|---|
| 522 | if (!IsValid())
|
|---|
| 523 | {
|
|---|
| 524 | gLog << pre << "Sequence: " << fFileName << " <invalid>" << endl;
|
|---|
| 525 | return;
|
|---|
| 526 | }
|
|---|
| 527 | gLog << "# Path: " << GetRcName() << endl;
|
|---|
| 528 | gLog << "# Name: " << GetName() << endl;
|
|---|
| 529 | gLog << endl;
|
|---|
| 530 | if (pre.IsNull())
|
|---|
| 531 | gLog << "Sequence: " << fSequence << endl;
|
|---|
| 532 | if (fMonteCarlo)
|
|---|
| 533 | gLog << pre << "MonteCarlo: Yes" << endl;
|
|---|
| 534 | gLog << pre <<"Period: " << fPeriod << endl;
|
|---|
| 535 | gLog << pre << "Night: " << fNight << endl << endl;
|
|---|
| 536 | gLog << pre << "LightCondition: ";
|
|---|
| 537 | switch (fLightCondition)
|
|---|
| 538 | {
|
|---|
| 539 | case kNA: gLog << "n/a" << endl; break;
|
|---|
| 540 | case kNoMoon: gLog << "NoMoon" << endl; break;
|
|---|
| 541 | case kTwilight: gLog << "Twilight" << endl; break;
|
|---|
| 542 | case kMoon: gLog << "Moon" << endl; break;
|
|---|
| 543 | case kDay: gLog << "Day" << endl; break;
|
|---|
| 544 | }
|
|---|
| 545 | gLog << pre << "Start: " << fStart << endl;
|
|---|
| 546 | gLog << pre << "LastRun: " << fLastRun << endl;
|
|---|
| 547 | gLog << pre << "NumEvents: " << fNumEvents << endl;
|
|---|
| 548 | gLog << pre << "Project: " << fProject << endl;
|
|---|
| 549 | gLog << pre << "Source: " << fSource << endl;
|
|---|
| 550 | gLog << pre << "TriggerTable: " << fTriggerTable << endl;
|
|---|
| 551 | gLog << pre << "HvSettings: " << fHvSettings << endl << endl;
|
|---|
| 552 | if (fRuns.GetSize()>0)
|
|---|
| 553 | {
|
|---|
| 554 | gLog << pre << "Runs:";
|
|---|
| 555 | for (int i=0; i<fRuns.GetSize(); i++)
|
|---|
| 556 | gLog << " " << fRuns[i];
|
|---|
| 557 | gLog << endl;
|
|---|
| 558 | }
|
|---|
| 559 | if (fCalRuns.GetSize()>0)
|
|---|
| 560 | {
|
|---|
| 561 | gLog << pre << "CalRuns:";
|
|---|
| 562 | for (int i=0; i<fCalRuns.GetSize(); i++)
|
|---|
| 563 | gLog << " " << fCalRuns[i];
|
|---|
| 564 | gLog << endl;
|
|---|
| 565 | }
|
|---|
| 566 | if (fPedRuns.GetSize()>0)
|
|---|
| 567 | {
|
|---|
| 568 | gLog << pre << "PedRuns:";
|
|---|
| 569 | for (int i=0; i<fPedRuns.GetSize(); i++)
|
|---|
| 570 | gLog << " " << fPedRuns[i];
|
|---|
| 571 | gLog << endl;
|
|---|
| 572 | }
|
|---|
| 573 | if (fDatRuns.GetSize()>0)
|
|---|
| 574 | {
|
|---|
| 575 | gLog << pre << "DatRuns:";
|
|---|
| 576 | for (int i=0; i<fDatRuns.GetSize(); i++)
|
|---|
| 577 | gLog << " " << fDatRuns[i];
|
|---|
| 578 | gLog << endl;
|
|---|
| 579 | }
|
|---|
| 580 | if (fExclRuns.GetSize()>0)
|
|---|
| 581 | {
|
|---|
| 582 | gLog << pre << "Exclude:";
|
|---|
| 583 | for (int i=0; i<fExclRuns.GetSize(); i++)
|
|---|
| 584 | gLog << " " << fExclRuns[i];
|
|---|
| 585 | gLog << endl;
|
|---|
| 586 | }
|
|---|
| 587 | if (!fDataPath.IsNull())
|
|---|
| 588 | gLog << endl << pre << "DataPath: " << fDataPath << endl;
|
|---|
| 589 |
|
|---|
| 590 | gLog << endl << pre << "Comment: " << fComment << endl;
|
|---|
| 591 | }
|
|---|
| 592 |
|
|---|
| 593 | // --------------------------------------------------------------------------
|
|---|
| 594 | //
|
|---|
| 595 | // Add all ped runs from the sequence to MDirIter.
|
|---|
| 596 | // If path==0 fDataPath is used instead. If it is also empty
|
|---|
| 597 | // the standard path of the data-center is assumed.
|
|---|
| 598 | // If you have the runs locally use path="."
|
|---|
| 599 | // Using raw=kTRUE you get correspodning raw-files setup.
|
|---|
| 600 | // Return the number of files added.
|
|---|
| 601 | //
|
|---|
| 602 | // Runs which are in fExlRuns are ignored.
|
|---|
| 603 | //
|
|---|
| 604 | UInt_t MSequence::SetupPedRuns(MDirIter &iter, const char *path, Bool_t raw) const
|
|---|
| 605 | {
|
|---|
| 606 | return SetupRuns(iter, fPedRuns, raw?kRawPed:kRootPed, path);
|
|---|
| 607 | }
|
|---|
| 608 |
|
|---|
| 609 | // --------------------------------------------------------------------------
|
|---|
| 610 | //
|
|---|
| 611 | // Add all data runs from the sequence to MDirIter.
|
|---|
| 612 | // If path==0 fDataPath is used instead. If it is also empty
|
|---|
| 613 | // the standard path of the data-center is assumed.
|
|---|
| 614 | // If you have the runs locally use path="."
|
|---|
| 615 | // Using raw=kTRUE you get correspodning raw-files setup.
|
|---|
| 616 | // Return the number of files added.
|
|---|
| 617 | //
|
|---|
| 618 | // Runs which are in fExlRuns are ignored.
|
|---|
| 619 | //
|
|---|
| 620 | UInt_t MSequence::SetupDatRuns(MDirIter &iter, const char *path, Bool_t raw) const
|
|---|
| 621 | {
|
|---|
| 622 | return SetupRuns(iter, fDatRuns, raw?kRawDat:kRootDat, path);
|
|---|
| 623 | }
|
|---|
| 624 |
|
|---|
| 625 | // --------------------------------------------------------------------------
|
|---|
| 626 | //
|
|---|
| 627 | // Add all runs from the sequence to MDirIter.
|
|---|
| 628 | // If path==0 fDataPath is used instead. If it is also empty
|
|---|
| 629 | // the standard path of the data-center is assumed.
|
|---|
| 630 | // If you have the runs locally use path="."
|
|---|
| 631 | // Using raw=kTRUE you get correspodning raw-files setup.
|
|---|
| 632 | // Return the number of files added.
|
|---|
| 633 | //
|
|---|
| 634 | // Runs which are in fExlRuns are ignored.
|
|---|
| 635 | //
|
|---|
| 636 | UInt_t MSequence::SetupAllRuns(MDirIter &iter, const char *path, Bool_t raw) const
|
|---|
| 637 | {
|
|---|
| 638 | return SetupRuns(iter, fRuns, raw?kRawAll:kRootAll, path);
|
|---|
| 639 | }
|
|---|
| 640 |
|
|---|
| 641 | // --------------------------------------------------------------------------
|
|---|
| 642 | //
|
|---|
| 643 | // Add all calibration runs from the sequence to MDirIter.
|
|---|
| 644 | // If path==0 fDataPath is used instead. If it is also empty
|
|---|
| 645 | // the standard path of the data-center is assumed.
|
|---|
| 646 | // If you have the runs locally use path="."
|
|---|
| 647 | // Using raw=kTRUE you get correspodning raw-files setup.
|
|---|
| 648 | // Return the number of files added.
|
|---|
| 649 | //
|
|---|
| 650 | // Runs which are in fExlRuns are ignored.
|
|---|
| 651 | //
|
|---|
| 652 | UInt_t MSequence::SetupCalRuns(MDirIter &iter, const char *path, Bool_t raw) const
|
|---|
| 653 | {
|
|---|
| 654 | return SetupRuns(iter, fCalRuns, raw?kRawCal:kRootCal, path);
|
|---|
| 655 | }
|
|---|
| 656 |
|
|---|
| 657 | // --------------------------------------------------------------------------
|
|---|
| 658 | //
|
|---|
| 659 | // Add all data runs from the sequence to MDirIter.
|
|---|
| 660 | // If path==0 fDataPath is used instead. If it is also empty
|
|---|
| 661 | // the standard path of the data-center is assumed.
|
|---|
| 662 | // If you have the runs locally use path="."
|
|---|
| 663 | // Using raw=kTRUE you get correspodning raw-files setup.
|
|---|
| 664 | // Return the number of files added.
|
|---|
| 665 | //
|
|---|
| 666 | // Runs which are in fExlRuns are ignored.
|
|---|
| 667 | //
|
|---|
| 668 | UInt_t MSequence::SetupDatRuns(MDirIter &iter, FileType_t type, const char *path) const
|
|---|
| 669 | {
|
|---|
| 670 | return SetupRuns(iter, fDatRuns, type, path);
|
|---|
| 671 | }
|
|---|
| 672 |
|
|---|
| 673 | // --------------------------------------------------------------------------
|
|---|
| 674 | //
|
|---|
| 675 | // If you want to add runs manually, use this function.
|
|---|
| 676 | //
|
|---|
| 677 | UInt_t MSequence::AddRuns(UInt_t first, UInt_t last, TArrayI *runs)
|
|---|
| 678 | {
|
|---|
| 679 | if (last<first)
|
|---|
| 680 | {
|
|---|
| 681 | *fLog << warn << "MSequence::AddRuns - WARNING: Last runnumber " << last;
|
|---|
| 682 | *fLog << " smaller than first " << first << "... ignored." << endl;
|
|---|
| 683 | return 0;
|
|---|
| 684 | }
|
|---|
| 685 | if (!IsValid())
|
|---|
| 686 | {
|
|---|
| 687 | *fLog << inf << "Setting Sequence number to #" << first << endl;
|
|---|
| 688 | fSequence = first;
|
|---|
| 689 | }
|
|---|
| 690 |
|
|---|
| 691 | const UInt_t nall = fRuns.GetSize();
|
|---|
| 692 | const UInt_t nrun = runs ? runs->GetSize() : 0;
|
|---|
| 693 | const UInt_t add = last-first+1;
|
|---|
| 694 |
|
|---|
| 695 | fRuns.Set(nall+add);
|
|---|
| 696 | if (runs)
|
|---|
| 697 | runs->Set(nrun+add);
|
|---|
| 698 |
|
|---|
| 699 | for (UInt_t i=0; i<add; i++)
|
|---|
| 700 | {
|
|---|
| 701 | fRuns[nall+i] = first+i;
|
|---|
| 702 | if (runs)
|
|---|
| 703 | (*runs)[nrun+i] = first+i;
|
|---|
| 704 | }
|
|---|
| 705 | return add;
|
|---|
| 706 | }
|
|---|
| 707 |
|
|---|
| 708 | Bool_t MSequence::IsContained(const TArrayI &arr, UInt_t run) const
|
|---|
| 709 | {
|
|---|
| 710 | for (int i=0; i<arr.GetSize(); i++)
|
|---|
| 711 | if (run==(UInt_t)arr[i])
|
|---|
| 712 | return kTRUE;
|
|---|
| 713 | return kFALSE;
|
|---|
| 714 | }
|
|---|
| 715 |
|
|---|
| 716 | /*
|
|---|
| 717 | // --------------------------------------------------------------------------
|
|---|
| 718 | //
|
|---|
| 719 | // Check if num is found in arr. If it is found remove it. Copy all
|
|---|
| 720 | // following entries one back and decrease the array size by 1.
|
|---|
| 721 | //
|
|---|
| 722 | void MSequence::ExcludeRun(TArrayI &arr, UInt_t num)
|
|---|
| 723 | {
|
|---|
| 724 | UInt_t *ptr = (UInt_t*)arr.GetArray();
|
|---|
| 725 |
|
|---|
| 726 | Int_t i = 0;
|
|---|
| 727 | for (i=0; i<arr.GetSize(); i++)
|
|---|
| 728 | if (ptr[i]==num)
|
|---|
| 729 | break;
|
|---|
| 730 |
|
|---|
| 731 | if (i==arr.GetSize())
|
|---|
| 732 | return;
|
|---|
| 733 |
|
|---|
| 734 | for (; i<arr.GetSize()-1; i++)
|
|---|
| 735 | ptr[i] = ptr[i+1];
|
|---|
| 736 |
|
|---|
| 737 | arr.Set(arr.GetSize()-1);
|
|---|
| 738 | }
|
|---|
| 739 | */
|
|---|
| 740 |
|
|---|
| 741 | // --------------------------------------------------------------------------
|
|---|
| 742 | //
|
|---|
| 743 | // Exclude this run (i.e. add it to fExclRuns)
|
|---|
| 744 | //
|
|---|
| 745 | void MSequence::ExcludeRun(UInt_t num)
|
|---|
| 746 | {
|
|---|
| 747 | /*
|
|---|
| 748 | ExcludeRun(fRuns, num);
|
|---|
| 749 | ExcludeRun(fCalRuns, num);
|
|---|
| 750 | ExcludeRun(fPedRuns, num);
|
|---|
| 751 | ExcludeRun(fDatRuns, num);
|
|---|
| 752 | */
|
|---|
| 753 |
|
|---|
| 754 | if (IsExcluded(num))
|
|---|
| 755 | return;
|
|---|
| 756 |
|
|---|
| 757 | // set new entry
|
|---|
| 758 | const Int_t n = fExclRuns.GetSize();
|
|---|
| 759 | fExclRuns.Set(n+1);
|
|---|
| 760 | fExclRuns[n] = num;
|
|---|
| 761 |
|
|---|
| 762 | MJob::SortArray(fExclRuns);
|
|---|
| 763 | }
|
|---|
| 764 |
|
|---|
| 765 | // --------------------------------------------------------------------------
|
|---|
| 766 | //
|
|---|
| 767 | // Exclude all runs which are found in the list, e.g. "100 102 105"
|
|---|
| 768 | //
|
|---|
| 769 | void MSequence::ExcludeRuns(TString runs)
|
|---|
| 770 | {
|
|---|
| 771 | TArrayI data;
|
|---|
| 772 | Split(runs, data);
|
|---|
| 773 |
|
|---|
| 774 | for (int i=0; i<data.GetSize(); i++)
|
|---|
| 775 | ExcludeRun(data[i]);
|
|---|
| 776 | }
|
|---|
| 777 |
|
|---|
| 778 | const TString MSequence::GetExcludedRuns() const
|
|---|
| 779 | {
|
|---|
| 780 | TString rc;
|
|---|
| 781 | for (int i=0; i<fExclRuns.GetSize(); i++)
|
|---|
| 782 | {
|
|---|
| 783 | rc += fExclRuns[i];
|
|---|
| 784 | rc += " ";
|
|---|
| 785 | }
|
|---|
| 786 | return rc(0, rc.Length()-1);
|
|---|
| 787 | }
|
|---|
| 788 |
|
|---|
| 789 | // --------------------------------------------------------------------------
|
|---|
| 790 | //
|
|---|
| 791 | // If you want to change or set the night manually.
|
|---|
| 792 | // The Format is
|
|---|
| 793 | // SetNight("yyyy-mm-dd");
|
|---|
| 794 | //
|
|---|
| 795 | void MSequence::SetNight(const char *txt)
|
|---|
| 796 | {
|
|---|
| 797 | TString night(txt);
|
|---|
| 798 | night += " 00:00:00";
|
|---|
| 799 | fNight.SetSqlDateTime(night);
|
|---|
| 800 |
|
|---|
| 801 | fPeriod = MAstro::GetMagicPeriod(fNight.GetMjd());
|
|---|
| 802 | }
|
|---|
| 803 |
|
|---|
| 804 | // --------------------------------------------------------------------------
|
|---|
| 805 | //
|
|---|
| 806 | // If the sequence name seq is just a digit it is inflated to a full
|
|---|
| 807 | // path following the datacenter standard.
|
|---|
| 808 | //
|
|---|
| 809 | // Returns if file accessible or not.
|
|---|
| 810 | //
|
|---|
| 811 | Bool_t MSequence::InflatePath(TString &seq, Bool_t ismc)
|
|---|
| 812 | {
|
|---|
| 813 | if (seq.IsDigit())
|
|---|
| 814 | {
|
|---|
| 815 | const Int_t numseq = seq.Atoi();
|
|---|
| 816 | seq = "/magic/";
|
|---|
| 817 | if (ismc)
|
|---|
| 818 | seq += "montecarlo/";
|
|---|
| 819 | seq += Form("sequences/%04d/sequence%08d.txt", numseq/10000, numseq);
|
|---|
| 820 | gLog << inf << "Inflated sequence file: " << seq << endl;
|
|---|
| 821 | }
|
|---|
| 822 |
|
|---|
| 823 | if (!gSystem->AccessPathName(seq, kFileExists))
|
|---|
| 824 | return kTRUE;
|
|---|
| 825 |
|
|---|
| 826 | gLog << err << "Sorry, sequence file '" << seq << "' doesn't exist." << endl;
|
|---|
| 827 | return kFALSE;
|
|---|
| 828 | }
|
|---|
| 829 |
|
|---|