// -----------------------------------------------------------------------------
// ----- SetCuts.C                                                         -----
// -----                                                                   -----
// ----- created by C. Simon on 2018-06-08                                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------

void SetCuts()
{
  // NOTE: The processes "PAIR", "COMP", "PHOT", "PFIS", "DRAY", "ANNI", "BREM",
  //       "HADR", "MUNU", and "DCAY" can be enabled without resulting secondaries
  //       being stored/processed by Geant (option '2').

  gMC->SetProcess("PAIR", 1);    // pair production mechanism
  gMC->SetProcess("COMP", 1);    // Compton scattering
  gMC->SetProcess("PHOT", 1);    // photo electric effect
  gMC->SetProcess("PFIS", 0);    // photo fission mechanism
  gMC->SetProcess("DRAY", 1);    // delta rays mechanism
  gMC->SetProcess("ANNI", 1);    // positron annihilation
  gMC->SetProcess("BREM", 1);    // bremsstrahlung
  gMC->SetProcess("HADR", 1);    // hadronic interactions
  gMC->SetProcess("MUNU", 1);    // muon nuclear interactions
  gMC->SetProcess("DCAY", 1);    // Decay mechanism
  gMC->SetProcess("LOSS", 1);    // energy loss
  gMC->SetProcess("MULS", 1);    // multiple scattering
  gMC->SetProcess("RAYL", 1);    // Rayleigh scattering
  gMC->SetProcess("STRA", 0);    // Photo-Absorption Ionization model
//  gMC->SetProcess("SYNC", 0);    // Synchrotron radiation in magnetic field
  gMC->SetProcess("CKOV", 1);    // Cherenkov production

  Double_t cut1 = 1.0E-3; // [GeV]
  Double_t tofmax = 1.0;  // [s]

  gMC->SetCut("CUTGAM", cut1);   // gammas
  gMC->SetCut("CUTELE", cut1);   // electrons
  gMC->SetCut("CUTHAD", cut1);   // charged hadrons
  gMC->SetCut("CUTNEU", cut1);   // neutral hadrons
  gMC->SetCut("CUTMUO", cut1);   // muons
  gMC->SetCut("BCUTE", cut1);    // electron bremsstrahlung
  gMC->SetCut("BCUTM", cut1);    // muon and hadron bremsstrahlung
  gMC->SetCut("DCUTE", cut1);    // delta-rays by electrons
  gMC->SetCut("DCUTM", cut1);    // delta-rays by muons
  gMC->SetCut("PPCUTM", cut1);   // direct pair production by muons
  gMC->SetCut("TOFMAX", tofmax); // time of flight cut

}