QuirksPhysicsTool Class Reference

#include <QuirkPhysicsTool.h>

Inheritance diagram for QuirksPhysicsTool:

Collaboration diagram for QuirksPhysicsTool:

Public Member Functions
	QuirksPhysicsTool (const std::string &type, const std::string &name, const IInterface *parent)
	Standard constructor. More...
virtual	~QuirksPhysicsTool ()
	Destructor. More...
virtual StatusCode	initialize () override final
	Initialize method. More...
virtual void	ConstructParticle () override final
virtual void	ConstructProcess () override final
virtual QuirksPhysicsTool *	GetPhysicsOption () override final

Private Attributes
G4double	m_mass
G4double	m_charge
G4int	m_pdgid
G4double	m_stringForce
G4double	m_firstStringLength
G4double	m_maxBoost
G4double	m_maxMergeT
G4double	m_maxMergeMag
G4int	m_enableDebug
G4double	m_debugStep
G4int	m_numDebugSteps

Detailed Description

Tool for the concrete implementation of a Physics List selection class

Author

Edoardo Farina

Date

15-05-2015

Definition at line 20 of file QuirkPhysicsTool.h.

Constructor & Destructor Documentation

QuirksPhysicsTool()

QuirksPhysicsTool::QuirksPhysicsTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Standard constructor.

Definition at line 49 of file QuirkPhysicsTool.cxx.

  : base_class ( type, nam , parent )
{
  m_physicsOptionType = G4AtlasPhysicsOption::Type::BSMPhysics;
 std::ifstream f;
    f.open("quirks_setup.txt");
    f >> m_mass;                m_mass *= CLHEP::MeV;
    f >> m_charge;              m_charge *= CLHEP::eplus;
    f >> m_pdgid;
    f >> m_stringForce;         m_stringForce *= CLHEP::MeV/CLHEP::mm;
    f >> m_firstStringLength;   m_firstStringLength *= CLHEP::mm;
    f >> m_maxBoost;
    f >> m_maxMergeT;           m_maxMergeT *= CLHEP::mm;
    f >> m_maxMergeMag;         m_maxMergeMag *= CLHEP::mm;
    f >> m_enableDebug;
    if (m_enableDebug) {
        f >> m_debugStep;           m_debugStep *= CLHEP::mm;
        f >> m_numDebugSteps;
    }
    f.close();
}

~QuirksPhysicsTool()

QuirksPhysicsTool::~QuirksPhysicsTool ( )

virtual

Destructor.

Definition at line 76 of file QuirkPhysicsTool.cxx.

{
 
}

Member Function Documentation

ConstructParticle()

void QuirksPhysicsTool::ConstructParticle ( )

finaloverridevirtual

Definition at line 97 of file QuirkPhysicsTool.cxx.

{
  ATH_MSG_DEBUG("ConstructParticle for the Quirks being run");
 
    Quirk* q1 = new Quirk(
        /*name*/ "quirk+",    m_mass,               /*width*/ 0.0,       +m_charge,
        /*2*spin*/ 1,         /*parity*/ 0,         /*C-conjugation*/ 0,
        /*2*Isospin*/ 0,      /*2*Isospin3*/ 0,     /*G-parity*/ 0,
        /*type*/ "lepton",    /*lepton number*/ +1, /*baryon number*/ 0, +m_pdgid,
        /*stable*/ true,      /*lifetime*/ -1.0,    /*decay table*/ 0,
        /*shortlived*/ false
    );
    new Quirk(
        /*name*/ "quirk-",    m_mass,               /*width*/ 0.0,       -m_charge,
        /*2*spin*/ 1,         /*parity*/ 0,         /*C-conjugation*/ 0,
        /*2*Isospin*/ 0,      /*2*Isospin3*/ 0,     /*G-parity*/ 0,
        /*type*/ "lepton",    /*lepton number*/ -1, /*baryon number*/ 0, -m_pdgid,
        /*stable*/ true,      /*lifetime*/ -1.0,    /*decay table*/ 0,
        /*shortlived*/ false
    );
 
    InfracolorForce& s = q1->GetStringIn();
    s.SetStringForce(m_stringForce);
    s.SetFirstStringLength(m_firstStringLength);
    s.SetMaxBoost(m_maxBoost);
    s.SetMaxMergeT(m_maxMergeT);
    s.SetMaxMergeMag(m_maxMergeMag);
}

ConstructProcess()

void QuirksPhysicsTool::ConstructProcess ( )

finaloverridevirtual

Definition at line 126 of file QuirkPhysicsTool.cxx.

{
  ATH_MSG_DEBUG(" ConstructProcess for Quirks being run");
    PARTICLEITERATOR->reset();
    while( (*PARTICLEITERATOR)() ){
        G4ParticleDefinition* particle = PARTICLEITERATOR->value();
        if (dynamic_cast<Quirk*>(particle) != 0) {
            G4ProcessManager* pmanager = particle->GetProcessManager();
            //G4cout << "# processes = " << pmanager->GetProcessListLength() << G4endl;
            //for (int i = 0; i < pmanager->GetProcessList()->size(); i++) {
            //    G4cout << "process i: " << (*pmanager->GetProcessList())[i]->GetProcessName() << G4endl;
            //}
            while (pmanager->GetProcessListLength() != 0) pmanager->RemoveProcess(0);
            G4VProcess* aTransportation = new QuirkTransportation;
            pmanager->AddProcess(aTransportation);
            pmanager->SetProcessOrderingToFirst(aTransportation, idxAlongStep);
            pmanager->SetProcessOrderingToFirst(aTransportation, idxPostStep);
            pmanager->AddProcess(new G4hMultipleScattering,-1, 1, 1);
            pmanager->AddProcess(new G4hIonisation,        -1, 2, 2);
            pmanager->AddProcess(new G4hBremsstrahlung,    -1, 3, 3);
            pmanager->AddProcess(new G4hPairProduction,    -1, 4, 4);
            //pmanager->AddProcess(new TestProcess,           -1, 1, 1);
            pmanager->AddProcess(new QuirkWatcher,          -1,-1, 5);
            //G4cout << "# processes = " << pmanager->GetProcessListLength() << G4endl;
            //for (int i = 0; i < pmanager->GetProcessList()->size(); i++) {
            //    G4cout << "process i: " << (*pmanager->GetProcessList())[i]->GetProcessName() << G4endl;
            //}
         }
    }
    G4RunManager* runManager = G4RunManager::GetRunManager();
    runManager->SetNumberOfAdditionalWaitingStacks(1);
    G4UserStackingAction* defaultStackingAction = G4EventManager::GetEventManager()->GetUserStackingAction();
    ATH_MSG_WARNING("This tool is interfering with the normal G4Atlas workflow, by forcing a differentcustom StackingAction instead of the default one. This is acceptable only as a temporary solution, and should be fixed asap.");
    runManager->SetUserAction(new QuirkStackingAction(defaultStackingAction));
 
    // ADS this must be done on the python side
 //f (m_enableDebug) {
 //   DebugSteppingAction* debugger = new DebugSteppingAction(m_debugStep, m_numDebugSteps);
    //#//f QUIRKS_STANDALONE
 //   runManager->SetUserAction(debugger);
    //#//
 //   FADS::FadsSteppingAction::GetSteppingAction()->RegisterAction(debugger);
    //#//f
 //
}

GetPhysicsOption()

QuirksPhysicsTool * QuirksPhysicsTool::GetPhysicsOption ( )

finaloverridevirtual

Definition at line 92 of file QuirkPhysicsTool.cxx.

{
  return this;
}

initialize()

StatusCode QuirksPhysicsTool::initialize ( )

finaloverridevirtual

Initialize method.

Definition at line 84 of file QuirkPhysicsTool.cxx.

{
  ATH_MSG_DEBUG("QuirksPhysicsTool initialize( )");
  this->SetPhysicsName(name());
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_charge

G4double QuirksPhysicsTool::m_charge

private

Definition at line 38 of file QuirkPhysicsTool.h.

m_debugStep

G4double QuirksPhysicsTool::m_debugStep

private

Definition at line 46 of file QuirkPhysicsTool.h.

m_enableDebug

G4int QuirksPhysicsTool::m_enableDebug

private

Definition at line 45 of file QuirkPhysicsTool.h.

m_firstStringLength

G4double QuirksPhysicsTool::m_firstStringLength

private

Definition at line 41 of file QuirkPhysicsTool.h.

m_mass

G4double QuirksPhysicsTool::m_mass

private

Definition at line 37 of file QuirkPhysicsTool.h.

m_maxBoost

G4double QuirksPhysicsTool::m_maxBoost

private

Definition at line 42 of file QuirkPhysicsTool.h.

m_maxMergeMag

G4double QuirksPhysicsTool::m_maxMergeMag

private

Definition at line 44 of file QuirkPhysicsTool.h.

m_maxMergeT

G4double QuirksPhysicsTool::m_maxMergeT

private

Definition at line 43 of file QuirkPhysicsTool.h.

m_numDebugSteps

G4int QuirksPhysicsTool::m_numDebugSteps

private

Definition at line 47 of file QuirkPhysicsTool.h.

m_pdgid

G4int QuirksPhysicsTool::m_pdgid

private

Definition at line 39 of file QuirkPhysicsTool.h.

m_stringForce

G4double QuirksPhysicsTool::m_stringForce

private

Definition at line 40 of file QuirkPhysicsTool.h.

---

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

• QuirkPhysicsTool.h
• QuirkPhysicsTool.cxx