move_field_2_run.C

Go to the documentation of this file.

// macro for moving magnetic field voltage from parameter values to table RUN_
class UniqueRunNumber;
void move_field_2_run(int period)
{
    UniqueRunNumber* run_numbers;
    int run_count = UniRun::GetRunNumbers(run_numbers, period, 1, period, 100000);
    if (run_count <= 0)
        return;
 
    const char* detector_name = "magnet";
    const char* parameter_name = "average_field_mv";
 
    for (int i = 0; i < run_count; i++) {
        // get the average magnetic field from Runs of the Unified Condition Database
        UniDetectorParameter* pDetectorParameter = UniDetectorParameter::GetDetectorParameter(
            detector_name, parameter_name, run_numbers[i].period_number, run_numbers[i].run_number);
        if (pDetectorParameter == nullptr) {
            cout << "No field voltage was found for run (" << run_numbers[i].period_number << ":"
                 << run_numbers[i].run_number << ")" << endl;
            continue;
        }
        double field_voltage = ((DoubleValue*)pDetectorParameter->GetValue())->value;
        delete pDetectorParameter;
 
        // get run
        UniRun* pRun = UniRun::GetRun(run_numbers[i].period_number, run_numbers[i].run_number);
        if (pRun == nullptr) {
            cout << "The function encountered with errors: no experimental run was found ("
                 << run_numbers[i].period_number << ":" << run_numbers[i].run_number << "). This run will be skipped!"
                 << endl;
            continue;
        }
        pRun->SetFieldVoltage(&field_voltage);
        delete pRun;
 
        cout << "Average magnetic field for run " << run_numbers[i].period_number << "-" << run_numbers[i].run_number
             << ": " << field_voltage << " mV" << endl;
    }
 
    delete run_numbers;
    cout << "Macro finished successfully" << endl;
}
 
void convert_current_2_voltage(int period)
{
    UniqueRunNumber* run_numbers;
    int run_count = UniRun::GetRunNumbers(run_numbers, period, 1, period, 100000);
    if (run_count <= 0)
        return;
 
    const char* detector_name = "magnet";
    const char* parameter_name = "field_current";
 
    for (int i = 0; i < run_count; i++) {
        // get the field current from parameter values of the Unified Condition Database
        UniDetectorParameter* pDetectorParameter = UniDetectorParameter::GetDetectorParameter(
            detector_name, parameter_name, run_numbers[i].period_number, run_numbers[i].run_number);
        if (pDetectorParameter == nullptr) {
            cout << "No field current was found for run (" << run_numbers[i].period_number << ":"
                 << run_numbers[i].run_number << ")" << endl;
            continue;
        }
        int field_current = ((IntValue*)pDetectorParameter->GetValue())->value;
        delete pDetectorParameter;
 
        // convert current to voltage approximately
        double field_voltage = field_current / 900.0 * 55.87;
 
        // get run
        UniRun* pRun = UniRun::GetRun(run_numbers[i].period_number, run_numbers[i].run_number);
        if (pRun == nullptr) {
            cout << "The function encountered with errors: no experimental run was found ("
                 << run_numbers[i].period_number << ":" << run_numbers[i].run_number << "). This run will be skipped!"
                 << endl;
            continue;
        }
 
        double* current_voltage = pRun->GetFieldVoltage();
        if (current_voltage == nullptr) {
            pRun->SetFieldVoltage(&field_voltage);
            cout << "Field voltage was written for run " << run_numbers[i].period_number << "-"
                 << run_numbers[i].run_number << ": " << field_voltage << " mV" << endl;
        } else
            cout << "Field voltage exists for run " << run_numbers[i].period_number << "-" << run_numbers[i].run_number
                 << ": " << *current_voltage << " mV" << endl;
 
        delete pRun;
    }
 
    delete run_numbers;
    cout << "Macro finished successfully" << endl;
}