mcbm: exclude much from default tracking, improve box generator

macro/mcbm/mcbm_reco_event.C

@@ -258,8 +258,8 @@ void mcbm_reco_event(Int_t nEvents = 10, TString dataset = "data/test",
 
 
   // -----   Track reconstruction   ------------------------------------------
-  Double_t beamWidthX = 0.1;
-  Double_t beamWidthY = 0.1;
+  Double_t beamWidthX  = 0.1;
+  Double_t beamWidthY  = 0.1;
   Int_t iGenCor        = 1;
   Double_t dScalFac    = 1.;
   Double_t dChi2Lim2   = 3.5;
@@ -383,6 +383,10 @@ void mcbm_reco_event(Int_t nEvents = 10, TString dataset = "data/test",
   // L1 tracking
   auto l1 = (debugWithMC) ? new CbmL1("CA", 5, 1) : new CbmL1("CA");
   l1->SetMcbmMode();
+  l1->DisableTrackingStation(cbm::algo::ca::EDetectorID::kMuch, 0);
+  l1->DisableTrackingStation(cbm::algo::ca::EDetectorID::kMuch, 1);
+  l1->DisableTrackingStation(cbm::algo::ca::EDetectorID::kMuch, 2);
+
   // User configuration example for CA:
   //l1->SetConfigUser(srcDir + "/macro/L1/configs/ca_params_user_example.yaml");
   run->AddTask(l1);

macro/mcbm/mcbm_transport_boxgen.C

@@ -108,12 +108,38 @@ void mcbm_transport_boxgen(Int_t nEvents = 10,
 
   // ----- Box generator of a single -----------------------
   FairBoxGenerator* boxGen = new FairBoxGenerator(pdg, multiplicity);
-  boxGen->SetPRange(0.1, 3.);
+  boxGen->SetXYZ(targetPosX, targetPosY, targetPosZ);
+  boxGen->SetPRange(100., 100.);
   boxGen->SetPhiRange(0., 360.);
-  boxGen->SetThetaRange(0, 50.);
+  boxGen->SetThetaRange(0., 50.);
   boxGen->SetCosTheta();
   boxGen->Init();
+
+  run.SetEngine(kGeant3);
   run.AddInput(boxGen);
+
+  CbmGeant3Settings g3set;
+
+  g3set.SetProcessPairProduction(0);
+  g3set.SetProcessComptonScattering(0);
+  g3set.SetProcessPhotoEffect(0);
+  g3set.SetProcessPhotoFission(0);
+  g3set.SetProcessDeltaRay(0);
+  g3set.SetProcessAnnihilation(0);
+  g3set.SetProcessBremsstrahlung(0);
+  g3set.SetProcessHadronicInteraction(0);
+  g3set.SetProcessMuonNuclearInteraction(0);
+  g3set.SetProcessDecay(0);
+  //g3set.SetProcessEnergyLossModel(0); // energy loss is required by the STS digitization
+  g3set.SetProcessMultipleScattering(0);
+
+  g3set.SetProcessRayleighScattering(0);
+  g3set.SetProcessCherenkovProduction(0);
+  g3set.SetProcessEneryLossStraggling(0);
+
+  run.SetGeant3Settings(&g3set);
+
+
   // ----- Box generator of a single proton per event -----------------------
 
   run.SetOutFileName(outFile, overwrite);
@@ -128,6 +154,8 @@ void mcbm_transport_boxgen(Int_t nEvents = 10,
   if (nEvents <= 10)  // store only for small number of events
     run.StoreTrajectories();
   run.SetRandomSeed(randomSeed);
+
+
   run.Run(nEvents);
   // ------------------------------------------------------------------------