diff --git a/macro/mcbm/mcbm_reco_event.C b/macro/mcbm/mcbm_reco_event.C index 9ca28f2fa10c0015c507ff533b23e6429f28a06e..bbb8a21af129c8e0ab6b35429138a91b47f0c8d3 100644 --- a/macro/mcbm/mcbm_reco_event.C +++ b/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); diff --git a/macro/mcbm/mcbm_transport_boxgen.C b/macro/mcbm/mcbm_transport_boxgen.C index e867d1202303e5576ffdb2eb4b96434732c5059d..b5b5ca88be025121ea45577b7d5356da4a413315 100644 --- a/macro/mcbm/mcbm_transport_boxgen.C +++ b/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); // ------------------------------------------------------------------------