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 37 of file GarfieldPhysicsList.hh.

Constructor & Destructor Documentation

GarfieldPhysicsList()

GarfieldPhysicsList::GarfieldPhysicsList

(

)

Definition at line 50 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 66 of file GarfieldPhysicsList.cc.

{}

Member Function Documentation

AddParameterisation()

void GarfieldPhysicsList::AddParameterisation ( )

protected

Definition at line 68 of file GarfieldPhysicsList.cc.

                                              {
  GarfieldPhysics* garfieldPhysics = GarfieldPhysics::GetInstance();
 
  std::string ionizationModel = garfieldPhysics->GetIonizationModel();
 
  auto 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);
 
    auto particleName = particle->GetParticleName();
    if (garfieldPhysics->FindParticleName(particleName, "garfield")) {
      pmanager->AddDiscreteProcess(fastSimProcess_garfield);
    }
 
    if (garfieldPhysics->FindParticleName(particleName, "geant4")) {
      double eMin = MeV * garfieldPhysics->GetMinEnergyMeVParticle(
          particleName, "geant4");
      double eMax = MeV * garfieldPhysics->GetMaxEnergyMeVParticle(
          particleName, "geant4");
      if (ionizationModel == "PAI") {
        G4PAIModel* pai = new G4PAIModel(particle, "G4PAIModel");
        if (particleName == "e-" || particleName == "e+") {
          config->SetExtraEmModel(particleName, "eIoni", pai,
                                  "RegionGarfield", eMin, eMax, pai);
        } else if (particleName == "mu-" || particleName == "mu+") {
          config->SetExtraEmModel(particleName, "muIoni", pai,
                                  "RegionGarfield", eMin, eMax, pai);
        } else if (particleName == "proton" ||
                   particleName == "pi+" || particleName == "pi-") {
          config->SetExtraEmModel(particleName, "hIoni", pai,
                                  "RegionGarfield", eMin, eMax, pai);
        } else if (particleName == "alpha" || particleName == "He3" ||
                   particleName == "GenericIon") {
          config->SetExtraEmModel(particleName, "ionIoni", pai,
                                  "RegionGarfield", eMin, eMax, pai);
        }
      } else if (ionizationModel == "PAIPhot") {
        G4PAIPhotModel* paiPhot = new G4PAIPhotModel(particle, "G4PAIModel");
        if (particleName == "e-" || particleName == "e+") {
          config->SetExtraEmModel(particleName, "eIoni", paiPhot,
                                  "RegionGarfield", eMin, eMax, paiPhot);
        } else if (particleName == "mu-" || particleName == "mu+") {
          config->SetExtraEmModel(particleName, "muIoni", paiPhot, 
                                  "RegionGarfield", eMin, eMax, paiPhot);
        } else if (particleName == "proton" ||
                   particleName == "pi+" || particleName == "pi-") {
          config->SetExtraEmModel(particleName, "hIoni", paiPhot,
                                  "RegionGarfield", eMin, eMax, paiPhot);
        } else if (particleName == "alpha" || particleName == "He3" ||
                   particleName == "GenericIon") {
          config->SetExtraEmModel(particleName, "ionIoni", paiPhot, 
                                  "RegionGarfield", eMin, eMax, paiPhot);
        }
      }
    }
  }
}

Referenced by ConstructProcess().

ConstructParticle()

void GarfieldPhysicsList::ConstructParticle ( )

virtual

Definition at line 150 of file GarfieldPhysicsList.cc.

                                            {
  G4VModularPhysicsList::ConstructParticle();
}

ConstructProcess()

void GarfieldPhysicsList::ConstructProcess ( )

virtual

Definition at line 154 of file GarfieldPhysicsList.cc.

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

SetCuts()

void GarfieldPhysicsList::SetCuts ( )

virtual

Definition at line 132 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