In This Package:

DsPhysConsGeneral.cc

Go to the documentation of this file.

#include "DsPhysConsGeneral.h"

#include "GaudiKernel/DeclareFactoryEntries.h" 

DECLARE_TOOL_FACTORY(DsPhysConsGeneral);

DsPhysConsGeneral::DsPhysConsGeneral(const std::string& type,
                                     const std::string& name,
                                     const IInterface* parent)
    : GiGaPhysConstructorBase(type,name,parent)
{
    declareProperty("LowProdCut", m_lowProdCut=250*eV,
                    "Lower bound on production energy");
    declareProperty("HighProdCut", m_highProdCut=100*GeV,
                    "Upper bound on production energy");

    declareProperty("GammaCut",m_gammaCut=1.0*mm, 
                    "Cut for gammas");
    declareProperty("ElectronCut",m_electronCut=100.0*micrometer,
                    "Cut for electrons");
    declareProperty("PositronCut", m_positronCut=100.0*micrometer,
                    "Cut for positrons");
}

DsPhysConsGeneral::~DsPhysConsGeneral()
{
}

#include "G4ParticleTypes.hh"
#include "G4BosonConstructor.hh"
#include "G4LeptonConstructor.hh"
#include "G4MesonConstructor.hh"
#include "G4BaryonConstructor.hh"
#include "G4IonConstructor.hh"
#include "G4ShortLivedConstructor.hh"
#include "G4Decay.hh"
#include "G4RadioactiveDecay.hh"
#include "DsG4RadioactiveDecay.hh"
#include "G4ProcessManager.hh"
#include "G4ProductionCutsTable.hh"
#include "G4ParticlePropertyTable.hh"

void DsPhysConsGeneral::ConstructParticle()
{
    {   // All bosons
        G4BosonConstructor con;
        con.ConstructParticle();
    }
    {   // All leptons
        G4LeptonConstructor con;
        con.ConstructParticle();
    }
    { /* Adopted a subset of mesons to avoid segv at the very end of
       * the run (seems to be related with B mesons and hLowEIoni???).
       * Carbon-copied the meson list from the example in
       * advanced/underground/physics. Jianglai 05/08/2006
       */
        G4PionPlus::PionPlusDefinition();
        G4PionMinus::PionMinusDefinition();
        G4PionZero::PionZeroDefinition();
        G4KaonPlus::KaonPlusDefinition();
        G4KaonMinus::KaonMinusDefinition();
        G4Eta::EtaDefinition();
        G4EtaPrime::EtaPrimeDefinition();
        G4KaonZero::KaonZeroDefinition();
        G4AntiKaonZero::AntiKaonZeroDefinition();
        G4KaonZeroLong::KaonZeroLongDefinition();
        G4KaonZeroShort::KaonZeroShortDefinition();
    }
    {   // All baryons
        G4BaryonConstructor con;
        con.ConstructParticle();
    }
    {   // All ions
        G4IonConstructor con;
        con.ConstructParticle();  
    }

    {   // all short lived (resonaces and quarks)
        G4ShortLivedConstructor con;
        con.ConstructParticle();  
    }
    
    G4ParticlePropertyTable* propTable = G4ParticlePropertyTable::GetParticlePropertyTable();
    assert(propTable);
    G4ParticlePropertyData* photonData = propTable->GetParticleProperty(G4OpticalPhoton::Definition());
    assert(photonData);
    photonData->SetPDGEncoding(20022);
    photonData->SetAntiPDGEncoding(20022);
    if(propTable->SetParticleProperty(*photonData))
      info() << "Set PDG code for opticalphoton to " << G4OpticalPhoton::Definition()->GetPDGEncoding() << endreq;
    else 
      warning() << "Failed to reset PDG code on opticalphoton.. it's still set to " 
                << G4OpticalPhoton::Definition()->GetPDGEncoding() << endreq;
} 

void DsPhysConsGeneral::ConstructProcess()
{
    // can't call this from a GiGaPhysConstructorBase, but
    // G4VModularPhysicsList will do it for us.
    // AddTransportation();

    G4Decay* theDecayProcess = new G4Decay();
    theParticleIterator->reset();
    while( (*theParticleIterator)() ) {
        G4ParticleDefinition* particle = theParticleIterator->value();
        G4ProcessManager* pmanager = particle->GetProcessManager();
        if (theDecayProcess->IsApplicable(*particle)) {
            pmanager->AddDiscreteProcess(theDecayProcess);
            pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
            pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
        }
    }

    const G4IonTable* theIonTable 
        = G4ParticleTable::GetParticleTable()->GetIonTable();
    //G4RadioactiveDecay* theRadioactiveDecay = new G4RadioactiveDecay();
    DsG4RadioactiveDecay* theRadioactiveDecay = new DsG4RadioactiveDecay();
    for (G4int i=0; i<theIonTable->Entries(); i++) {
        G4String particleName = theIonTable->GetParticle(i)->GetParticleName();
        if (particleName == "GenericIon") {
            G4ProcessManager* pmanager 
                = theIonTable->GetParticle(i)->GetProcessManager();
            pmanager ->AddProcess(theRadioactiveDecay);
            pmanager ->SetProcessOrdering(theRadioactiveDecay, idxPostStep);
            pmanager ->SetProcessOrdering(theRadioactiveDecay, idxAtRest);
        }
    }

}

#if 0
void DsPhysConsGeneral::SetCuts()
{
    // special for low energy physics
    G4ProductionCutsTable::GetProductionCutsTable()->
        SetEnergyRange(m_lowProdCut,m_highProdCut);

    // set cut values for gamma at first and for e- second and next for e+,
    // because some processes for e+/e- need cut values for gamma 
    SetCutValue(m_gammaCut, "gamma");
    SetCutValue(m_electronCut, "e-");
    SetCutValue(m_positronCut, "e+");
}
#endif

---