{ TStopwatch timer; timer.Start(); // Run over all events (or blocks) Int_t nev = -1; //Int_t nev = 2; // Create source with unpackers ---------------------------------------------- //FairRemoteSource* source = new FairRemoteSource("lxi047"); //FairLmdSource* source = new FairLmdSource(); FairbinSource* source = new FairbinSource(); source->AddFile("/home/mala/Echanges/R3B_Simulation/r3b_dataread_Marcello/Mixed-alpha-source/2014-03-06-Eth-b0.bin"); // char* strCalDir = "/home/work/example/"; // TODO: set proper values to the variable below when the binary file is from GSI MBS DAQ is used. // Short_t type = 0; Short_t subType = 0; Short_t procId = 0; Short_t subCrate = -1; // All sub-crates Short_t control = 0; source->AddUnpacker(new R3BStarTrackUnpack("", type, subType, procId, subCrate, control)); // --------------------------------------------------------------------------- // Create online run --------------------------------------------------------- FairRunOnline* run = new FairRunOnline(source); run->SetOutputFile("startrack_raw_data.root"); run->SetGenerateHtml(kTRUE); // --------------------------------------------------------------------------- // Create ALADIN field map --------------------------------------------------- R3BAladinFieldMap* magField = new R3BAladinFieldMap("AladinMaps"); Double_t fMeasCurrent = 2500.;// I_current [A] magField->SetCurrent(fMeasCurrent); magField->SetScale(1.); run->SetField(magField); // --------------------------------------------------------------------------- // Add analysis task --------------------------------------------------------- R3BStarTrackRawAna* ana = new R3BStarTrackRawAna(); run->AddTask(ana); // --------------------------------------------------------------------------- // Initialize ---------------------------------------------------------------- run->Init(); // --------------------------------------------------------------------------- // Runtime data base --------------------------------------------------------- FairRuntimeDb* rtdb = run->GetRuntimeDb(); R3BFieldPar* fieldPar = (R3BFieldPar*) rtdb->getContainer("R3BFieldPar"); fieldPar->SetParameters(magField); fieldPar->setChanged(); Bool_t kParameterMerged = kTRUE; FairParRootFileIo* parOut = new FairParRootFileIo(kParameterMerged); parOut->open("r3bpar.root"); rtdb->setOutput(parOut); rtdb->saveOutput(); rtdb->print(); // --------------------------------------------------------------------------- // Run ----------------------------------------------------------------------- run->Run(nev, 0); // --------------------------------------------------------------------------- timer.Stop(); Double_t rtime = timer.RealTime(); Double_t ctime = timer.CpuTime(); cout << endl << endl; cout << "Macro finished succesfully." << endl; cout << "Real time " << rtime << " s, CPU time " << ctime << "s" << endl << endl; }