GarfieldPhysicsList Class Reference

#include <GarfieldPhysicsList.hh>

Inheritance diagram for GarfieldPhysicsList:

Public Member Functions
	GarfieldPhysicsList ()
virtual	~GarfieldPhysicsList ()
virtual void	SetCuts ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()

Protected Member Functions
void	AddParameterisation ()

Detailed Description

Definition at line 39 of file GarfieldPhysicsList.hh.

Constructor & Destructor Documentation

GarfieldPhysicsList()

GarfieldPhysicsList::GarfieldPhysicsList

(

)

Definition at line 49 of file GarfieldPhysicsList.cc.

                                         :
        G4VModularPhysicsList() {
    G4int verb = 0;
    SetVerboseLevel(verb);
    defaultCutValue = 1 * CLHEP::mm;
    QGSP_BERT_HP *physicsList = new QGSP_BERT_HP;
    for (G4int i = 0;; ++i) {
        G4VPhysicsConstructor* elem =
                const_cast<G4VPhysicsConstructor*>(physicsList->GetPhysics(i));
        if (elem == NULL)
            break;
        G4cout << "RegisterPhysics: " << elem->GetPhysicsName() << G4endl;
        RegisterPhysics(elem);
    }
 
}

~GarfieldPhysicsList()

GarfieldPhysicsList::~GarfieldPhysicsList ( )

virtual

Definition at line 68 of file GarfieldPhysicsList.cc.

                                          {
 
}

Member Function Documentation

AddParameterisation()

void GarfieldPhysicsList::AddParameterisation ( )

protected

Definition at line 72 of file GarfieldPhysicsList.cc.

                                              {
 
    GarfieldPhysics* garfieldPhysics = GarfieldPhysics::GetInstance();
 
    std::string ionizationModel = garfieldPhysics->GetIonizationModel();
 
    G4FastSimulationManagerProcess* fastSimProcess_garfield =
            new G4FastSimulationManagerProcess("G4FSMP_garfield");
 
        auto theParticleIterator = GetParticleIterator();
    theParticleIterator->reset();
    while ((*theParticleIterator)()) {
 
        G4ParticleDefinition* particle = theParticleIterator->value();
        G4ProcessManager* pmanager = particle->GetProcessManager();
        G4EmConfigurator* config =
                G4LossTableManager::Instance()->EmConfigurator();
        G4LossTableManager::Instance()->SetVerbose(1);
 
        if (garfieldPhysics->FindParticleName(particle->GetParticleName(),
                "garfield")) {
            pmanager->AddDiscreteProcess(fastSimProcess_garfield);
        }
 
        if (garfieldPhysics->FindParticleName(particle->GetParticleName(),
                "geant4")) {
            double minEnergy_MeV = garfieldPhysics->GetMinEnergyMeVParticle(
                    particle->GetParticleName(), "geant4");
            double maxEnergy_MeV = garfieldPhysics->GetMaxEnergyMeVParticle(
                    particle->GetParticleName(), "geant4");
            if (ionizationModel == "PAI") {
 
                G4PAIModel* pai = new G4PAIModel(particle, "G4PAIModel");
                if (particle->GetParticleName() == "e-"
                        || particle->GetParticleName() == "e+") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "eIoni", pai, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, pai);
 
                } else if (particle->GetParticleName() == "mu-"
                        || particle->GetParticleName() == "mu+") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "muIoni", pai, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, pai);
 
                } else if (particle->GetParticleName() == "proton"
                        || particle->GetParticleName() == "pi+"
                        || particle->GetParticleName() == "pi-") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "hIoni", pai, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, pai);
 
                } else if (particle->GetParticleName() == "alpha"
                        || particle->GetParticleName() == "He3"
                        || particle->GetParticleName() == "GenericIon") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "ionIoni", pai, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, pai);
 
                }
 
            } else if (ionizationModel == "PAIPhot") {
 
                G4PAIPhotModel* paiPhot = new G4PAIPhotModel(particle,
                        "G4PAIModel");
                if (particle->GetParticleName() == "e-"
                        || particle->GetParticleName() == "e+") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "eIoni", paiPhot, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, paiPhot);
 
                } else if (particle->GetParticleName() == "mu-"
                        || particle->GetParticleName() == "mu+") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "muIoni", paiPhot, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, paiPhot);
 
                } else if (particle->GetParticleName() == "proton"
                        || particle->GetParticleName() == "pi+"
                        || particle->GetParticleName() == "pi-") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "hIoni", paiPhot, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, paiPhot);
 
                } else if (particle->GetParticleName() == "alpha"
                        || particle->GetParticleName() == "He3"
                        || particle->GetParticleName() == "GenericIon") {
 
                    config->SetExtraEmModel(particle->GetParticleName(),
                            "ionIoni", paiPhot, "RegionGarfield", minEnergy_MeV * MeV,
                            maxEnergy_MeV * MeV, paiPhot);
 
                }
            }
        }
    }
}

Referenced by ConstructProcess().

ConstructParticle()

void GarfieldPhysicsList::ConstructParticle ( )

virtual

Definition at line 197 of file GarfieldPhysicsList.cc.

                                            {
    G4VModularPhysicsList::ConstructParticle();
}

ConstructProcess()

void GarfieldPhysicsList::ConstructProcess ( )

virtual

Definition at line 201 of file GarfieldPhysicsList.cc.

                                           {
    G4VModularPhysicsList::ConstructProcess();
    AddParameterisation();
}

SetCuts()

void GarfieldPhysicsList::SetCuts ( )

virtual

Definition at line 178 of file GarfieldPhysicsList.cc.

                                  {
    G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(100. * eV,
            100. * TeV);
 
    SetCutsWithDefault();
 
    G4Region *region = G4RegionStore::GetInstance()->GetRegion(
            "RegionGarfield");
    G4ProductionCuts * cuts = new G4ProductionCuts();
    cuts->SetProductionCut(1 * um, G4ProductionCuts::GetIndex("gamma"));
    cuts->SetProductionCut(1 * um, G4ProductionCuts::GetIndex("e-"));
    cuts->SetProductionCut(1 * um, G4ProductionCuts::GetIndex("e+"));
    if (region) {
        region->SetProductionCuts(cuts);
    }
 
    DumpCutValuesTable();
}

---

The documentation for this class was generated from the following files:

• GarfieldPhysicsList.hh
• GarfieldPhysicsList.cc