dch_noise_hist.C

Go to the documentation of this file.

// macro for getting 'noise channels' parameter value
// example runs for all run numbers from 12 to 688 (if run number is not present then error "table record wasn't found"
// wiil arise - it's ok)
void dch_noise_hist()
{
    TCanvas* c = new TCanvas("c", "Noise channels", 0, 0, 640, 480);
    // TGraph2D* gr2 = new TGraph2D();
    TH3I* histo = new TH3I("h3", "Noise channels", 20, 1, 20, 70, 1, 70, 689, 0, 688);
    histo->SetMarkerColor(kBlue);
    histo->GetXaxis()->SetTitle("Slot");
    histo->GetYaxis()->SetTitle("Channel");
    histo->GetZaxis()->SetTitle("Run №");
 
    UniqueRunNumber* run_numbers;
    int run_count = UniRun::GetRunNumbers(run_numbers, 1, 12, 3, 688);
    if (run_count <= 0)
        return;
 
    for (int i = 0; i < run_count; i++) {
        // get noise parameter values presented by IIStructure: Int+Int (slot:channel)
        UniDetectorParameter* pDetectorParameter = UniDetectorParameter::GetDetectorParameter(
            "DCH1", "noise", run_numbers[i].period_number,
            run_numbers[i].run_number);   //(detector_name, parameter_name, period_number, run_number)
        if (pDetectorParameter != nullptr) {
            vector<UniValue*> pValues;
            pDetectorParameter->GetValue(pValues);
 
            cout << "Element count: " << pValues.size() << " for run " << run_numbers[i].period_number << ":"
                 << run_numbers[i].run_number << " (period:run)" << endl;
            for (int j = 0; j < pValues.size(); j++)
                histo->Fill(((IIValue*)pValues[j])->value1, ((IIValue*)pValues[j])->value2, run_numbers[i].run_number);
 
            // clean memory after work
            for (int i = 0; i < pValues.size(); i++)
                delete pValues[i];
            delete pDetectorParameter;
        }
    }
 
    delete run_numbers;
 
    histo->Draw();
}