+ Float_t cut = 1.e-3; // 1MeV cut by default
+ Float_t tofmax = 1.e10;
+
+ gMC->SetCut("CUTGAM", cut);
+ gMC->SetCut("CUTELE", cut);
+ gMC->SetCut("CUTNEU", cut);
+ gMC->SetCut("CUTHAD", cut);
+ gMC->SetCut("CUTMUO", cut);
+ gMC->SetCut("BCUTE", cut);
+ gMC->SetCut("BCUTM", cut);
+ gMC->SetCut("DCUTE", cut);
+ gMC->SetCut("DCUTM", cut);
+ gMC->SetCut("PPCUTM", cut);
+ gMC->SetCut("TOFMAX", tofmax);
+
+ // Set apply cuts
+ geant4->ProcessGeantCommand("/run/particle/applyCuts");
+ geant4->ProcessGeantCommand("/mcVerbose/geometryManager 2");
+ geant4->ProcessGeantCommand("/mcDet/volNameSeparator !");
+ geant4->ProcessGeantCommand("/mcControl/g3Defaults");
+ geant4->ProcessGeantCommand("/mcPhysics/setStackPopperSelection e+ e- pi+ pi- kaon+ kaon- gamma");
+// geant4->ProcessGeantCommand("/tracking/verbose 1");