/** @file FHodoLabSetup ** @author Florian Uhlig ** @date 20.06.2016 ** ** ROOT macro to read tsa files which have been produced with the new data transport ** Convert data into cbmroot format. ** Uses CbmFlibTestSource as source task. */ //void unpack_tsa(Int_t nEvt=100, TString FileId = "cosmic_2016110701_safe_4links_4") void unpack_tsa(Int_t nEvt=100, Double_t dDeltaT=50., Int_t iReqDet = 0x00001006, TString FileId = "r0001_20171215_1342") { TString srcDir = gSystem->Getenv("VMCWORKDIR"); TString inDir = "./input/star2018/" + FileId + "/"; TString inFile = "*.tsa"; // --- Specify number of events to be produced. // --- -1 means run until the end of the input file. Int_t bEvt=-1; Int_t nEvents = 10; if(nEvt>0) { bEvt=0; nEvents = nEvt; } // --- Specify output file name (this is just an example) TString outDir = "./data/"; TString Mode = Form("_DT%d_0x%08x",(Int_t)dDeltaT,iReqDet); TString outFile = outDir + FileId + Mode + ".root"; TString parFile = outDir + FileId + Mode + ".params.root"; // --- Set log output levels //FairLogger::GetLogger(); //FairLogger::GetLogger()->SetLogScreenLevel("WARNING"); FairLogger::GetLogger()->SetLogScreenLevel("INFO"); //FairLogger::GetLogger()->SetLogScreenLevel("DEBUG"); FairLogger::GetLogger()->SetLogVerbosityLevel("MEDIUM"); // --- Define parameter files TList *parFileList = new TList(); TString paramDir = "./"; TString paramFileTof = paramDir + "MapTofHD_v18b.par"; TObjString* parTofFileName = new TObjString(paramFileTof); parFileList->Add(parTofFileName); // --- Set debug level gDebug = 0; std::cout << std::endl; std::cout << ">>> unpack_tsa: output file is " << outFile << std::endl; // ======================================================================== // ======================================================================== std::cout << std::endl; std::cout << ">>> unpack_tsa: Initialising..." << std::endl; // Get4 Unpacker CbmUnpackTofStar2018* unpacker_tof = new CbmUnpackTofStar2018(); //argument is number of gDpb // unpacker_tof->SetEpochSuppressedMode(bEpSupp); //unpacker_tof->SetTShiftRef(1700000.); // Reference signal time shift in ns, run 008 <-> 0018 unpacker_tof->SetTShiftRef(21000.); // Reference signal time shift in ns, run 0026 // --- Source task CbmFlibTestSource* source = new CbmFlibTestSource(); source->SetMaxDeltaT(dDeltaT); source->SetReqDigiAddr(iReqDet); //0x00005006); // request diamond for output events source->SetReqMode(1); // Look for events with ANY detector in setup //source->SetFileName(inFile); source->AddPath(inDir,inFile); source->AddUnpacker(unpacker_tof, 0x60, 6);//gDPB A & B // --- Event header FairEventHeader* event = new CbmTbEvent(); event->SetRunId(1); // --- Run FairRunOnline *run = new FairRunOnline(source); run->SetOutputFile(outFile); run->SetEventHeader(event); FairRuntimeDb* rtdb = run->GetRuntimeDb(); // ----- Runtime database --------------------------------------------- Bool_t kParameterMerged = kTRUE; FairParRootFileIo* parOut = new FairParRootFileIo(kParameterMerged); FairParAsciiFileIo* parIn = new FairParAsciiFileIo(); parOut->open(parFile.Data()); parIn->open(parFileList, "in"); rtdb->setFirstInput(parIn); rtdb->setOutput(parOut); run->Init(); // --- Start run TStopwatch timer; timer.Start(); std::cout << ">>> unpack_tsa: Starting run..." << std::endl; if(bEvt == -1) run->Run(bEvt, 0); // run until end of input file else run->Run(0, nEvents); // process nEvents timer.Stop(); std::cout << "Processed " << std::dec << source->GetTsCount() << " timeslices" << std::endl; // --- End-of-run info Double_t rtime = timer.RealTime(); Double_t ctime = timer.CpuTime(); std::cout << std::endl << std::endl; std::cout << ">>> unpack_tsa: Macro finished successfully." << std::endl; std::cout << ">>> unpack_tsa: Output file is " << outFile << std::endl; std::cout << ">>> unpack_tsa: Real time " << rtime << " s, CPU time " << ctime << " s" << std::endl; std::cout << std::endl; /// --- Screen output for automatic tests std::cout << " Test passed" << std::endl; std::cout << " All ok " << std::endl; }