diff --git a/macro/mcbm/mcbm_reco.C b/macro/mcbm/mcbm_reco.C
index 0714f2a1781a53c44a081fd65a42f7c197ae4059..70c64ee793af1c98cf438a3239af869dbc03b2dc 100644
--- a/macro/mcbm/mcbm_reco.C
+++ b/macro/mcbm/mcbm_reco.C
@@ -13,11 +13,12 @@
// Binned tracker for track reconstruction
//
// V. Friese 11.06.2018
-// S. Roy 18.11.2021 // time-based mode event building and reconstruction. The "Real" event builder is used.
+// S. Roy 18.11.2021 // time-based mode event building and reconstruction. The "Real" event builder is used.
// --------------------------------------------------------------------------
-void mcbm_reco(Int_t nEvents = 10, TString dataset = "data/test", TString sEvBuildRaw = "Real", const char* setupName = "mcbm_beam_2022_05")
+void mcbm_reco(Int_t nEvents = 10, TString dataset = "data/test", TString sEvBuildRaw = "Real",
+ const char* setupName = "mcbm_beam_2022_05")
{
// ========================================================================
// Adjust this part according to your requirements
@@ -55,7 +56,6 @@ void mcbm_reco(Int_t nEvents = 10, TString dataset = "data/test", TString sEvBui
// setup->RemoveModule(ECbmModuleId::kTrd);
// ------------------------------------------------------------------------
-
// ----- Some global switches -----------------------------------------
Bool_t eventBased = !sEvBuildRaw.IsNull();
@@ -129,55 +129,54 @@ void mcbm_reco(Int_t nEvents = 10, TString dataset = "data/test", TString sEvBui
FairLogger::GetLogger()->SetLogVerbosityLevel(logVerbosity.Data());
// ------------------------------------------------------------------------
-
// ----- Raw event building from digis --------------------------------
if (sEvBuildRaw.EqualTo("Real", TString::ECaseCompare::kIgnoreCase)) {
- CbmTaskBuildRawEvents* evBuildRaw = new CbmTaskBuildRawEvents();
-
- //Choose between NoOverlap, MergeOverlap, AllowOverlap
- evBuildRaw->SetEventOverlapMode(EOverlapModeRaw::AllowOverlap);
-
- // Remove detectors where digis not found
- if (!useRich) evBuildRaw->RemoveDetector(kRawEventBuilderDetRich);
- if (!useMuch) evBuildRaw->RemoveDetector(kRawEventBuilderDetMuch);
- if (!usePsd) evBuildRaw->RemoveDetector(kRawEventBuilderDetPsd);
- if (!useTof) evBuildRaw->RemoveDetector(kRawEventBuilderDetTof);
- if (!useTrd) evBuildRaw->RemoveDetector(kRawEventBuilderDetTrd);
-
- if (!useSts) {
- std::cerr << "-E- " << myName << ": Sts must be present for raw event "
- << "building using ``Real2019'' option. Terminating macro." << std::endl;
- return;
- }
- // Set STS as reference detector
- evBuildRaw->SetReferenceDetector(kRawEventBuilderDetSts);
-
- // Use sliding window seed builder with STS
- // evBuildRaw->SetReferenceDetector(kRawEventBuilderDetUndef);
- // evBuildRaw->AddSeedTimeFillerToList(kRawEventBuilderDetSts);
- // evBuildRaw->SetSlidingWindowSeedFinder(10, 40, 100);
- // evBuildRaw->SetSeedFinderQa(true); // optional QA information for seed finder
+ CbmTaskBuildRawEvents* evBuildRaw = new CbmTaskBuildRawEvents();
- evBuildRaw->SetTsParameters(0.0, 1.e7, 0.0);
+ //Choose between NoOverlap, MergeOverlap, AllowOverlap
+ evBuildRaw->SetEventOverlapMode(EOverlapModeRaw::AllowOverlap);
- // Use CbmMuchDigi instead of CbmMuchBeamtimeDigi
- evBuildRaw->ChangeMuchBeamtimeDigiFlag(kFALSE);
+ // Remove detectors where digis not found
+ if (!useRich) evBuildRaw->RemoveDetector(kRawEventBuilderDetRich);
+ if (!useMuch) evBuildRaw->RemoveDetector(kRawEventBuilderDetMuch);
+ if (!usePsd) evBuildRaw->RemoveDetector(kRawEventBuilderDetPsd);
+ if (!useTof) evBuildRaw->RemoveDetector(kRawEventBuilderDetTof);
+ if (!useTrd) evBuildRaw->RemoveDetector(kRawEventBuilderDetTrd);
- evBuildRaw->SetTriggerMinNumber(ECbmModuleId::kSts, 100);
- evBuildRaw->SetTriggerMaxNumber(ECbmModuleId::kSts, -1);
- evBuildRaw->SetTriggerWindow(ECbmModuleId::kSts, -10, 65);
-
- run->AddTask(evBuildRaw);
- std::cout << "-I- " << myName << ": Added task " << evBuildRaw->GetName() << std::endl;
- eventBased = kTRUE;
- } //? Real raw event building
- else {
- std::cerr << "-E- " << myName << ": Unknown option " << sEvBuildRaw
- << " for raw event building! Terminating macro execution." << std::endl;
+ if (!useSts) {
+ std::cerr << "-E- " << myName << ": Sts must be present for raw event "
+ << "building using ``Real2019'' option. Terminating macro." << std::endl;
return;
}
- //? event-based reco
+ // Set STS as reference detector
+ evBuildRaw->SetReferenceDetector(kRawEventBuilderDetSts);
+
+ // Use sliding window seed builder with STS
+ // evBuildRaw->SetReferenceDetector(kRawEventBuilderDetUndef);
+ // evBuildRaw->AddSeedTimeFillerToList(kRawEventBuilderDetSts);
+ // evBuildRaw->SetSlidingWindowSeedFinder(10, 40, 100);
+ // evBuildRaw->SetSeedFinderQa(true); // optional QA information for seed finder
+
+ evBuildRaw->SetTsParameters(0.0, 1.e7, 0.0);
+
+ // Use CbmMuchDigi instead of CbmMuchBeamtimeDigi
+ evBuildRaw->ChangeMuchBeamtimeDigiFlag(kFALSE);
+
+ evBuildRaw->SetTriggerMinNumber(ECbmModuleId::kSts, 100);
+ evBuildRaw->SetTriggerMaxNumber(ECbmModuleId::kSts, -1);
+ evBuildRaw->SetTriggerWindow(ECbmModuleId::kSts, -10, 65);
+
+ run->AddTask(evBuildRaw);
+ std::cout << "-I- " << myName << ": Added task " << evBuildRaw->GetName() << std::endl;
+ eventBased = kTRUE;
+ } //? Real raw event building
+ else {
+ std::cerr << "-E- " << myName << ": Unknown option " << sEvBuildRaw
+ << " for raw event building! Terminating macro execution." << std::endl;
+ return;
+ }
+ //? event-based reco
// ------------------------------------------------------------------------
@@ -520,9 +519,9 @@ void mcbm_reco(Int_t nEvents = 10, TString dataset = "data/test", TString sEvBui
// ----- Database update ----------------------------------------------
-// rtdb->setOutput(parIo1);
-// rtdb->saveOutput();
-// rtdb->print();
+ // rtdb->setOutput(parIo1);
+ // rtdb->saveOutput();
+ // rtdb->print();
// ------------------------------------------------------------------------
@@ -547,17 +546,14 @@ void mcbm_reco(Int_t nEvents = 10, TString dataset = "data/test", TString sEvBui
std::cout << " All ok " << std::endl;
// ------------------------------------------------------------------------
// save all historgrams
-
+
gROOT->LoadMacro("save_hst.C");
- TString FSave = Form("save_hst(\"CluStatus%d_%d_Cal_%s.hst.root\")",
- iCalSet,
- iSel2in,
- cCalId.Data());
+ TString FSave = Form("save_hst(\"CluStatus%d_%d_Cal_%s.hst.root\")", iCalSet, iSel2in, cCalId.Data());
gInterpreter->ProcessLine(FSave.Data());
-
+
// ----- Resource monitoring ------------------------------------------
-if (hasFairMonitor) { // FairRoot Version >= 15.11
+ if (hasFairMonitor) { // FairRoot Version >= 15.11
// Extract the maximal used memory an add is as Dart measurement
// This line is filtered by CTest and the value send to CDash
FairSystemInfo sysInfo;
@@ -575,5 +571,5 @@ if (hasFairMonitor) { // FairRoot Version >= 15.11
tempMon->Print();
}
-// RemoveGeoManager();
+ // RemoveGeoManager();
}