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

Question No hits on shield 

Keywords: hits shielding energy loss
Forum: Hits, Digitization and Pileup
Date: 22 Feb, 2007
From: Kazutaka Nakahara <nakahara@post.kek.jp>

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

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

1 None: Re: No hits on shield   (Vladimir IVANTCHENKO - 22 Feb, 2007)
2 None: Re: No hits on shield   (michel maire - 22 Feb, 2007)
 Add Message Add Message
to: "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 ]