Message: Re: how2 access to G4PSCylinderSurfaceFlux Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: how2 access to G4PSCylinderSurfaceFlux 

Forum: Biasing and Scoring
Re: Question how2 access to G4PSCylinderSurfaceFlux (Aliya Mukanova)
Re: None Re: how2 access to G4PSCylinderSurfaceFlux (Makoto Asai)
Date: 29 Apr, 2014
From: Aliya Mukanova <Aliya Mukanova>

Hi, Makoto

Thank you for answer

I have written MyRun.cc that way

     G4double nSf = 0.;
     G4THitsMap<G4double>* nsfl
	        =  static_cast<G4THitsMap<G4double>*>(HCE->GetHC(nSurfID));
     for (itr = nsfl->GetMap()->begin(); itr != nsfl->GetMap()->end(); itr++) 
    {
     nSf = *(itr->second);
    }
    nSurf += nSf;

and then in Merge()

     nSurf += localRun -> nSurf;

Also I added GetNeutronSurfFlux in MyRun.hh

     G4double GetNeutrSurfFlux()   const { return nSurf; }

and called it from MyRunAction.cc as follows

    G4double neutSurf = run->GetNeutrSurfFlux();

However, it gives me Zero. This is my File.ouput:

Primary Neutrons 1000;

Neutron SurfFlux 0 /cm2;

Tritons in Li 972;

Alphas in Li 973;

Triton Track Length in Li 117.183 um.

Besides as you can see in results I have implemented this way with others scorers and they calculate well.

Could you point at the possible causes why it happens

Thanks Best regards, Aliya

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

1 None: Re: how2 access to G4PSCylinderSurfaceFlux   (Makoto Asai - 01 May, 2014)
(_ None: Re: how2 access to G4PSCylinderSurfaceFlux   (Aliya Mukanova - 05 May, 2014)
(_ Sad: Re: how2 access to G4PSCylinderSurfaceFlux   (Aliya Mukanova - 12 May, 2014)
 Add Message Add Message
to: "Re: how2 access to G4PSCylinderSurfaceFlux"

 Subscribe Subscribe

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