run_digi_task.C

Go to the documentation of this file.

// run_digi_task.C: runs analysis tasks on a simulated or experimental DIGI file
#include "BmnGemDigitQa.h"
 
void run_digi_task(TString digiFileName = "bmnsim.root", TString outFileName = "digi_qa.root", Int_t nEvents = 0)
{
    FairRunAna* fRunAna = new FairRunAna();
 
    // the same choice of the input source as in run_reco_bmn.C
    Bool_t isExp = (BmnFunctionSet::isSimulationFile(digiFileName) == 0);
    if (isExp) {
        Int_t run_period = -1, run_number = -1;   // filled from DigiRunHeader
        fRunAna->SetSource(new BmnFileSource(digiFileName, run_period, run_number));
    } else {
        fRunAna->SetSource(new FairFileSource(digiFileName));
        // parameter containers of the simulation file: geometry, magnetic field, run parameters
        FairParRootFileIo* parIo = new FairParRootFileIo();
        parIo->open(digiFileName.Data());
        fRunAna->GetRuntimeDb()->setFirstInput(parIo);
    }
    fRunAna->SetSink(new FairRootFileSink(outFileName));
 
    fRunAna->AddTask(new BmnGemDigitQa(isExp));
 
    fRunAna->Init();
    fRunAna->Run(0, nEvents);
}