JMergeTuneHV.cc File Reference

Auxiliary program to merge multiple HV tuning calibrations. More...

#include <iostream>
#include <iomanip>
#include <vector>
#include <set>
#include "JSystem/JDate.hh"
#include "JLang/JUUID.hh"
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
#include "Jeep/JProperties.hh"
#include "JLang/JPredicate.hh"
#include "JSupport/JMeta.hh"
#include "JSupport/JFilenameSupportkit.hh"
#include "JDetector/JDetectorCalibration.hh"
#include "JDB/JDBSupportkit.hh"
#include "JSon/JSon.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Auxiliary program to merge multiple HV tuning calibrations.

Note that only the calibration in the first given file will be maintained
in case a PMT appears in multiple files.

Author

bjung

Definition in file JMergeTuneHV.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 36 of file JMergeTuneHV.cc.

{
  using namespace std;
  using namespace JPP;

  string         login;
  string         locationID;
  int            elapsedTime = 0;
  
  vector<string> inputFiles;
  string         outputFile;
  
  int            debug;

  try {

    JProperties properties;

    properties.insert(gmake_property(login));
    properties.insert(gmake_property(locationID));
    properties.insert(gmake_property(elapsedTime));
    
    JParser<> zap("Auxiliary program to merge high-voltage calibration results.");

    zap['f'] = make_field(inputFiles,   "input  files");
    zap['o'] = make_field(outputFile,   "output file");
    zap['#'] = make_field(properties,   "database information") = JPARSER::initialised();
    zap['d'] = make_field(debug,        "debug")                = 2;

    zap(argc, argv);
    
  } catch(const exception &error) {
    
    FATAL(error.what() << endl);
  }

  if (login.empty() || locationID.empty()) {
    FATAL("Missing user information (please specify via -#login and -#locationID)." << endl);
  }


  const JUUID& UUID = JUUID::rndm();

  JDateAndTime       timer;

  timer.sub(elapsedTime);
  
  
  // Merge high-voltage calibrations
  
  JDBAPIVersion      DBAPIVersion;  
  set<string>        testTypes;
  string             metaInfoStr = MAKE_STRING(JMeta(argc, argv));

  JHVCalibration     merged;
  
  for (vector<string>::const_iterator fileIt = inputFiles.cbegin(); fileIt != inputFiles.cend(); ++fileIt) {

    NOTICE("Merging " << (*fileIt) << endl);

    if (isJSONFile(fileIt->c_str())) {

      json     js;
      ifstream ifs(fileIt->c_str());
      
      ifs >> js;
      ifs.close();

      // Extract data
    
      json::const_iterator i0 = js.find(APIVersion_t);

      if (i0 != js.cend()) {
      
        istringstream iss(i0->get<string>());
      
        iss >> DBAPIVersion;
      }

      JHVCalibration HVcals;
      
      json::const_iterator i1 = js.find(Data_t);
    
      if ( (DBAPIVersion.getMajorVersion() == 2) &&
           (i1 != js.cend() && i1->size() > 0) ) {

        set<string>::const_iterator j = testTypes.insert(testTypes.cbegin(), (*i1)[0].at(Test_t + Type_t).get<string>());

        JHVCalibration_t::setVersion(getDBVersionTuneHV(*j));
        
        HVcals         = (*i1)[0].at(Tests_t).get<JHVCalibration>();
        metaInfoStr   += (*i1)[0].at(Provenance_t + Info_t).at(Configuration_t).get<string>();
      
      } else {
        
        set<string>::const_iterator j = testTypes.insert(testTypes.cbegin(), js.at(Test_t + Type_t).get<string>());

        JHVCalibration_t::setVersion(getDBVersionTuneHV(*j));
        
        HVcals         = js.at(Tests_t).get<JHVCalibration>();
      }

      // Merge calibrations

      for (JHVCalibration::const_iterator cal = HVcals.cbegin(); cal != HVcals.cend(); ++cal) {
        
        if (find_if(merged.cbegin(), merged.cend(), make_predicate(&JHVCalibration_t::getNumber, cal->getNumber())) == merged.cend()) {
          merged.push_back(*cal);
        }
      }

    } else {

      WARNING(*fileIt << " is not a JSON file; skip" << endl);
      continue;
    }
  }

  if (testTypes.size() > 1) {
   
    WARNING("Ambiguous database test types; Assuming " << *testTypes.crbegin() << endl);
    
  } else if (testTypes.size() == 0) {
    
    FATAL("No database test type specified.");
  }


  json js;

  if (DBAPIVersion.getMajorVersion() == 2) {

    json error           = { {Message_t,              ""                      },
                             {Code_t,                 OK_t                    },
                             {Arguments_t,            json::array()           } };
    
    json metaData        = { {Configuration_t,        metaInfoStr             },
                             {UUID_t,                 MAKE_STRING(UUID)       } };
    
    json data            = { {Provenance_t + Info_t,  json(metaData)          },
                             {User_t,                 login                   },
                             {Location_t,             locationID              },
                             {Start_t + Time_t,       timer.toString()        },
                             {End_t   + Time_t,       timer().toString()      },
                             {Test_t  + Type_t,     *(testTypes.crbegin())    },
                             {Tests_t,                json(merged)            } };
      
    js[APIVersion_t]     = MAKE_STRING(DBAPIVersion);
    js[Data_t + Type_t]  = MAKE_STRING("ProductTestSession");
    js[Encoding_t]       = MAKE_STRING("NativeJSON");
    js[Error_t]          = json(error);
    js[Start_t]          = timer.toString();
    js[End_t]            = timer().toString();
    js[Data_t][0]        = json(data);
    
  } else {

    js[User_t]           = login;
    js[Location_t]       = locationID;
    js[Test_t + Type_t]  = *(testTypes.crbegin());
    js[Start_t + Time_t] = timer.toString();
    js[End_t + Time_t]   = timer().toString();
    js[Tests_t]          = json(merged);
  }
  
  ofstream ofs(outputFile.c_str());

  ofs << setw(2) << setprecision(8);
  ofs << js;

  ofs.close();

  return 0;
}