Message: GB02 Example Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question GB02 Example 

Forum: Biasing and Scoring
Date: Mar 26, 07:54
From: <modest24>

Dear Experts,

I am implementing the force collision scheme in GB02 to force collision in a volume. I have essential followed the example and copied relevant files so as to force neutron collision in a particular volume. I get an error about the implementation which I need your advice to resolve. My code is based on Hadr07 and I am using steppingAction and TrackingAction.

The code runs well if the number of events is less than 100, but nonetheless still produces an error message. The error is below:

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : BIAS.GEN.19
      issued by :  G4BOptrForceCollisionTrackData::~G4BOptrForceCollisionTrackData()
Track deleted while under G4BOptrForceCollision biasing scheme of operator `ForceCollisionForalpha'. Will result in inconsistencies.
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : BIAS.GEN.03
      issued by :  G4BiasingProcessInterface::PostStepDoIt(...)
 Negative interaction weight : w_I = 0 XS_I(phys) = inf XS_I(bias) = 2.02451e-25 step length = -0 Interaction law = `0x561d946c1b70'

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

Here is my code in detector construction for implementing the force interaction:

 for (G4int k=1; k<=fNbOfAbsor; k++) {
    G4Material* material = fAbsorMaterial[k];
    G4String matname = material->GetName();
     fSolidAbsor[k] =
      new G4Box(matname,fAbsorThickness[k]/2,fAbsorSizeYZ/2,fAbsorSizeYZ/2);

     fLogicAbsor[k] =
      new G4LogicalVolume(fSolidAbsor[k],           // solid
                          material,             // material
                          matname);             // name   

// force interaction in Volume 3 if(k==3){

  GB02BOptrMultiParticleForceCollision* testMany =  new GB02BOptrMultiParticleForceCollision();
  // testMany->AddParticle("gamma");
  testMany->AddParticle("neutron");
  // testMany->AddParticle("alpha");
  testMany->AttachTo(fLogicAbsor[3]);
  G4cout << " Attaching biasing operator " << testMany->GetName()
         << " to logical volume " << fLogicAbsor[3]->GetName()
         << G4endl;
}
    //

fXfront[k] = fXfront[k-1] + (fAbsorThickness[k-1]);    
    G4double xcenter = fXfront[k]+0.5*fAbsorThickness[k];
   G4double top= xcenter + 0.5*fAbsorThickness[k]; 
   G4double bottom =xcenter-0.5*fAbsorThickness[k];
    G4ThreeVector position = G4ThreeVector(xcenter,0.,0.);
    G4double thickness= fAbsorThickness[k];
//

     fPhysiAbsor[k]= new G4PVPlacement(0,                     //no rotation
                        position,             //position
                        fLogicAbsor[k],       //logical volume        
                        matname,               //name
                        logicWorld,            //mother
                        false,                 //no boulean operat
                        k);                    //copy number

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

1 None: Re: GB02 Example   (modest24 - Mar 27, 08:40)
 Add Message Add Message
to: "GB02 Example"

 Subscribe Subscribe

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