Message: Re: No hits on shield Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: No hits on shield 

Forum: Hits, Digitization and Pileup
Re: Question No hits on shield (Kazutaka Nakahara)
Date: 22 Feb, 2007
From: Vladimir IVANTCHENKO <vnivanch@mail.cern.ch>

On Thu, 22 Feb 2007, Kazutaka Nakahara wrote:

> *** Discussion title: Hits, Digitization and Pileup
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/hitsdigits/143"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> Hi,
> 
> I'm having some problems with registering hits for 4.26MeV positive
> muons incident on ~1.3 meters of lead. In principle, it should stop near
> the surface, but the muons pass right through.
> 
> To debug, I've set the stepsize to 1mm, and put statements in
> UserSteppingAction() to get the stepsize, Energy deposit, position, etc
> for each step. A sample below:
> 
> ----------------------------- Stepsize = 1 mm Edep = 0 eV Position =
> 157.369 -251.233 497.341 Momentum = 8.51397 -13.0988 25.967 K.E = 4.26
> MeV ----------------------------- At every step throughout the track,
> the energy deposit is 0. The physic processes I've included are (in my
> PhysicsList.cc)
> 
> -------------------------------
>  } else if( particleName == "mu+" || 
>                particleName == "mu-"    ) {
>            G4cout << "Constructing Mu+ and Mu- processes" << G4endl;
>       //muon  
>       pmanager->AddProcess(new G4MultipleScattering);
>       pmanager->AddProcess(new G4MuIonisation);
>       pmanager->AddDiscreteProcess(new G4MuBremsstrahlung);
>       pmanager->AddDiscreteProcess(new G4MuPairProduction);        
>       pmanager = G4MuonPlus::MuonPlus()->GetProcessManager();
>       pmanager->AddProcess(StepLimiter);
> ----------------------------------------
> 
> The  definition of the lead shielding is as follows:
> -------------------------------------
>     G4RotationMatrix rmlead;
>     rmlead.rotateY(0.);
>         G4ThreeVector positionshield1 = G4ThreeVector(0,0,184.8*cm);
>         solidLeadShield1 = new G4Tubs("lead_shield1", 32*cm,100*cm,133.5*cm,0,2*3.1415);
>         logicLeadShield1 = new G4LogicalVolume(solidLeadShield1,LeadShieldMater
>                         ,"lead_shield1",0,0,0);
>         physiLeadShield1 = new G4PVPlacement(0,positionshield1,
>                                                         logicLeadShield1,
>                                                         "lead_shield1",
>                                                         logicWorld,
>                                                         0,
>                                                         false,
>                                                         0);     
> ---------------------------------------
> 
> I've based everything on the N02 example, but I'm pretty sure I'm
> missing something VERY basic. Any help would be appreciated.
> 
> Thanks in advance, Kaz
> 

Hello,

Please, consult PhysicsLists in 
$G4INSTALL/source/physics_lists/builders/src/G4EmStandardPhysics.cc

or in $G4INSTALL/example/extended/electromagnetic/TestEm3 (any other 
example as well)

VI

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

 Add Message Add Message
to: "Re: No hits on shield"

 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 ]