Commit f7643dab authored by Valentina Fioretti's avatar Valentina Fioretti
Browse files

EXACRAD and AREMBESPhys update

parent b89bf1b1
Loading
Loading
Loading
Loading
+27 −0
Original line number Diff line number Diff line
@@ -84,10 +84,37 @@ G4VPhysicalVolume* EXACRAD::Construct() {
    World_log = gm.World_log;
    World_phys = ConstructGeometry(World_phys);
    
    ConfigureRegionCuts();

    return World_phys;
    
}

void EXACRAD::ConfigureRegionCuts()
{

    G4String regName;
    G4Region* region = nullptr;
    G4ProductionCuts* cuts = nullptr;

    // Target region
    G4double targetCut = 0.001 * mm;
    gm.config->readInto(targetCut, "REGION.TARGET.CUT");
    G4double targetCutProton = 0.001 * mm;
    gm.config->readInto(targetCutProton, "REGION.TARGET.PROTON.CUT");

    regName = "EXACRAD_target";
    region = G4RegionStore::GetInstance()->GetRegion(regName);
    if (region) {
        cuts = new G4ProductionCuts;
        cuts->SetProductionCut(targetCut);
        cuts->SetProductionCut(targetCutProton, G4ProductionCuts::GetIndex("proton"));
        region->SetProductionCuts(cuts);
    }
    
    G4cout << "[EXACRAD] Applying custom production cuts to region" << G4endl;

}

G4VPhysicalVolume* EXACRAD::ConstructGeometry(G4VPhysicalVolume* World_phys) {
    
+1 −0
Original line number Diff line number Diff line
@@ -61,6 +61,7 @@ public:
        G4VPhysicalVolume* Construct();
        void SetPhysicalWorld(G4VPhysicalVolume* World_phys);
        void DefineMaterials();
        void ConfigureRegionCuts();
        G4VPhysicalVolume* ConstructGeometry(G4VPhysicalVolume* World_phys);
    

+34 −487

File changed.

Preview size limit exceeded, changes collapsed.

+3 −32
Original line number Diff line number Diff line
//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
// ***************************************************************************
//                          AREMBESPhys.hh  -  description
//                             -------------------
//    Author		: Valentina Fioretti, Andrea Bulgarelli
//    creation date	: 10/02/2015
//    email		: fioretti@iasfbo.inaf.it, bulgarelli@iasfbo.inaf.it
// ***************************************************************************/
// Copyright 2025 Valentina Fioretti
//
@@ -92,6 +60,9 @@ private:

    std::vector<G4VPhysicsConstructor*> hadronAREMBESPhys;
    
    G4bool stepLimiterEnabled = false;  // default false
    G4bool SS_physics = false;
    
    G4double cutForGamma;
    G4double cutForElectron;
    G4double cutForPositron;
+231 −441

File changed.

Preview size limit exceeded, changes collapsed.

Loading