Message: Re: Empty plot??? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Idea Re: Empty plot??? 

Forum: Documentation and Examples
Re: Question Empty plot??? (Muneerah)
Date: 22 Jan, 2016
From: Muneerah <Muneerah>

Hi

I need help please. I am trying to simulate straightforward proton decay in silicon detector and find the best thickness (either 300 micrometer or 1 mm). I've started by this example (rdecay01) cuz ifound it the close one to my case. Please correct me if I'm in the wrong way.

I made some changes in the source files and the related include files. However, I got all the plots empty???? These changes are: In the DetectorConstruction.cc file:

 //added diodLog1, diodPhys1, physiWorld
DetectorConstruction::DetectorConstruction()
  :diodeLog1(0),diodePhys1(0), physiWorld(0)
{
  fWorldSize = 0.5*m;
}

and

//changed Air to nist
  G4NistManager* nist = G4NistManager::Instance();
  //add Air
  G4Material* Air = nist ->FindOrBuildMaterial("G4_AIR"); 

  //add the silicon detector
  G4Material* silicon = nist->FindOrBuildMaterial("G4_Si");

And

  //Define detector width (hy) and height(hz) and thickness (hx)
G4double diode_hy = 8*mm;
G4double diode_hz = 8*mm;
G4double diode_hx = 300* micrometer;

And Finally

//Silicon detector "Diode"

diodeBox = new G4Box("Diode", diode_hx, diode_hy, diode_hz);

diodeLog1 = new G4LogicalVolume(diodeBox, silicon, "Diode1");

//I think but not sure this is the position of the detectors G4double pos_x = 0*mm; G4double pos_y = 0*mm; G4double pos_z = 0*mm;

diodePhys1 = new G4PVPlacement(0, G4ThreeVector(pos_x, pos_y, pos_z), diodeLog1, "Diode1",logicWorld , false, 0);

In the HistoManager.cc file:

//add proton
    {	    "energy spectrum (%): proton",              //0

In PhysicsList.cc file:

 //changed to true
  radioactiveDecay->SetARM(true);               //Atomic Rearangement

In PrimeryGeneratorAction.cc file:

 //add proton
 G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4String particleName;
  G4ParticleDefinition* particle  = particleTable->FindParticle(particleName="proton");
  fParticleGun->SetParticleDefinition(particle);

  fParticleGun->SetParticleEnergy(2.05*MeV);

And commant the geantino part

In the RunAction.cc file:
 //add o
 analysisManager->ScaleH1(0,factor);

In TrackingAction.cc file:

  //add proton
  if (particle == G4Proton::Proton())  ih = 0;

These are all the changes that I've made so far and the result is empty plots!! Anyone can help?

Thanks in advance Muneerah

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

1 None: Re: Empty plot???   (michel maire - 22 Jan, 2016)
(_ None: Re: Empty plot???   (Muneerah - 22 Jan, 2016)
(_ None: Re: Empty plot???   (michel maire - 23 Jan, 2016)
(_ None: Re: Empty plot???   (Muneerah - 25 Jan, 2016)
(_ None: Re: Empty plot???   (Muneerah - 25 Jan, 2016)
 Add Message Add Message
to: "Re: Empty plot???"

 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 ]