Message: Forced-flight biasing: Cumulated weight is null ! Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Forced-flight biasing: Cumulated weight is null ! 

Forum: Biasing and Scoring
Date: 24 Jan, 2017
From: Jayson Vavrek <Jayson Vavrek>

I'm running geant4-10-00-patch-02 and attempting to use the G4BOptnForceFreeFlight() class for forced-flight biasing. While I can see in the visualization that my primary gammas are indeed flying freely through the volumes to which I've attached the biasing operator, it appears that the track weights are not being computed under the hood:

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : BIAS.GEN.06
      issued by :  G4BOptnForceFreeFlight::DenyProcessPostStepDoIt(...)
 Cumulated weight is null ! 
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

In my main.cc I have:

 physicsList *thePhysics = new physicsList();
G4GenericBiasingPhysics *theBiasingPhysics = new G4GenericBiasingPhysics();
theBiasingPhysics->Bias("gamma");
thePhysics->RegisterPhysics(theBiasingPhysics);
runManager->SetUserInitialization(thePhysics);

and I have written my own concrete operator class, forcedFlightOperator

file forcedFlightOperator.hh:

#ifndef forcedFlightOperator_hh
#define forcedFlightOperator_hh 1

#include "G4VBiasingOperator.hh"
#include "G4BOptnForceFreeFlight.hh"

class forcedFlightOperator : public G4VBiasingOperator
{
public:
	forcedFlightOperator();
	virtual ~forcedFlightOperator();

	G4BOptnForceFreeFlight *GetForceFreeFlightLaw() {return fForceFreeFlightInteractionLaw;};

private:
	G4BOptnForceFreeFlight *fForceFreeFlightInteractionLaw;

	// used for forced flight
	virtual G4VBiasingOperation *ProposeOccurenceBiasingOperation(const G4Track*, const G4BiasingProcessInterface* );

	// unused but required pure virtual functions
	virtual G4VBiasingOperation *ProposeNonPhysicsBiasingOperation(const G4Track* track, const G4BiasingProcessInterface* callingProcess) { return 0; };
        virtual G4VBiasingOperation *ProposeFinalStateBiasingOperation(const G4Track*, const G4BiasingProcessInterface* ) { return 0; };

};

#endif

file forcedFlightOperator.cc:

#include "forcedFlightOperator.hh"
forcedFlightOperator::forcedFlightOperator() : G4VBiasingOperator("forcedFlightOperator")
{
	fForceFreeFlightInteractionLaw = new G4BOptnForceFreeFlight("forcedFlight");
}

forcedFlightOperator::~forcedFlightOperator()
{
	delete fForceFreeFlightInteractionLaw;
}

G4VBiasingOperation *forcedFlightOperator::ProposeOccurenceBiasingOperation(const G4Track*, const G4BiasingProcessInterface*)
{
	return fForceFreeFlightInteractionLaw;
}

Am I missing something?

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

1 None: Re: Forced-flight biasing: Cumulated weight is null !   (Marc Verderi - 24 Jan, 2017)
(_ None: Re: Forced-flight biasing: Cumulated weight is null !   (Jayson Vavrek - 27 Jan, 2017)
(_ None: Re: Forced-flight biasing: Cumulated weight is null !   (Marc Verderi - 28 Jan, 2017)
(_ None: Re: Forced-flight biasing: Cumulated weight is null !   (Adam Geva - 12 Apr, 2017)
 Add Message Add Message
to: "Forced-flight biasing: Cumulated weight is null !"

 Subscribe Subscribe

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