Message: Re: PhysicsLists and another physics process Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: PhysicsLists and another physics process 

Keywords: a general scheme to do this
Forum: Physics List
Re: Question PhysicsLists and another physics process (Taylan Yetkin)
Re: None Re: PhysicsLists and another physics process (Vladimir Ivanchenko)
Re: None Re: PhysicsLists and another physics process (Taylan Yetkin)
Date: 23 Mar, 2004
From: Wellisch, J.P. <Wellisch, J.P.>

Hi,

  as discussed, but here again for more general consumption:

the simplest way is to write a physics constructor, inheriting from G4VPhysicsConstructor (let's call it MyAdditionalPhysics for the sake of this discussion), that constructs the Cerencov process, and registers it with the process managers in the particles for which you want Cerencov radiation.

Then, in your main, you can use the RegisterPhysics interface of ModularPhysicsList to regiter this.

Many greetings,

Hans-Peter.

P.S.: A few code snippets:

#ifndef MyAdditionalPhysics_hh
#define MyAdditionalPhysics_hh

class MyAdditionalPhysics : public G4VPhysicsConstructor
{
  public:
    MyAdditionalPhysics(G4String aS) : G4VPhysicsConstructor(aS) {}
    virtual void ConstructParticle() 
    {
     // here we do nothing - or may be optical photons?
    }
    virtual void ConstructProcess() 
    {
     // here we can construct and register Cerencov and other stuff.
     // Make sure nothing is there twice, though !!!
    }

  private:

};

#endif

and from the main:

...etc...

#include "MyAdditionalPhysics.hh"

...etc...

  QGSP * theQGSP = new QGSP;
  theQGSP->RegisterPhysics(new MyAdditionalPhysics("Cerencov") );
  runManager->SetUserInitialization(theQGSP);

...etc...

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "Re: PhysicsLists and another physics process"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]