Message: Scintillation process not invoked. Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Scintillation process not invoked. 

Forum: Processes Involving Optical Photons
Date: 11 Oct, 2012
From: Dat Vu <Dat Vu>

I am simulating the following process: Boron 10 + neutron -> Li 7 + alpha. Then: alpha + europium -> optical photon. I'm currently using Modular physics list. The first process works fine. The problem is at the second process. No optical photon/photon shows up at all although it is compiled and run very well. Is there anyway to register the optical physics process in G4VModularPhysicsList? Here is my code:

Physics list:

PhysicsList::PhysicsList()
:G4VModularPhysicsList()
{
  G4int verb = 1;  
  SetVerboseLevel(verb);
  // Hadron Elastic scattering
  RegisterPhysics( new G4HadronElasticPhysics(verb) );
  // Hadron Physics
  RegisterPhysics( new HadronPhysicsQGSP_BERT_HP(verb));
  // Ion Physics
  RegisterPhysics( new G4IonPhysics(verb));
  // Gamma Physics
  RegisterPhysics( new GammaPhysics("gamma"));
  // Optical Physics
  G4OpticalPhysics* optPhys = new G4OpticalPhysics();
  RegisterPhysics(optPhys);
  optPhys->SetTrackSecondariesFirst(kScintillation,true);
  optPhys->SetScintillationYieldFactor(10.0);
  optPhys->SetMaxNumPhotonsPerStep(1000);
  optPhys->SetScintillationByParticleType(false);

}

Optical properties of materials used:

  const G4int nEntries = 32;

  G4double PhotonEnergy[nEntries] =
            { 2.034*eV, 2.068*eV, 2.103*eV, 2.139*eV,
              2.177*eV, 2.216*eV, 2.256*eV, 2.298*eV,
              2.341*eV, 2.386*eV, 2.433*eV, 2.481*eV,
              2.532*eV, 2.585*eV, 2.640*eV, 2.697*eV,
              2.757*eV, 2.820*eV, 2.885*eV, 2.954*eV,
              3.026*eV, 3.102*eV, 3.181*eV, 3.265*eV,
              3.353*eV, 3.446*eV, 3.545*eV, 3.649*eV,
              3.760*eV, 3.877*eV, 4.002*eV, 4.136*eV };

  G4double RefractiveIndex1[nEntries] =
            { 1.3435, 1.344,  1.3445, 1.345,  1.3455,
              1.346,  1.3465, 1.347,  1.3475, 1.348,
              1.3485, 1.3492, 1.35,   1.3505, 1.351,
              1.3518, 1.3522, 1.3530, 1.3535, 1.354,
              1.3545, 1.355,  1.3555, 1.356,  1.3568,
              1.3572, 1.358,  1.3585, 1.359,  1.3595,
              1.36,   1.3608};

  G4double Absorption1[nEntries] =
           {3.448*m,  4.082*m,  6.329*m,  9.174*m, 12.346*m, 13.889*m,
           15.152*m, 17.241*m, 18.868*m, 20.000*m, 26.316*m, 35.714*m,
           45.455*m, 47.619*m, 52.632*m, 52.632*m, 55.556*m, 52.632*m,
           52.632*m, 47.619*m, 45.455*m, 41.667*m, 37.037*m, 33.333*m,
           30.000*m, 28.500*m, 27.000*m, 24.500*m, 22.000*m, 19.500*m,
           17.500*m, 14.500*m };

  G4MaterialPropertiesTable* myMPT1 = new G4MaterialPropertiesTable();

  myMPT1->AddProperty("RINDEX",       PhotonEnergy, RefractiveIndex1,nEntries)
        ->SetSpline(true);
  myMPT1->AddProperty("ABSLENGTH",    PhotonEnergy, Absorption1,     nEntries)
        ->SetSpline(true);

  myMPT1->AddConstProperty("SCINTILLATIONYIELD",50./eV);
  myMPT1->AddConstProperty("RESOLUTIONSCALE",1.0);
  myMPT1->AddConstProperty("FASTTIMECONSTANT", 1.*ns);
  myMPT1->AddConstProperty("SLOWTIMECONSTANT",10.*ns);
  myMPT1->AddConstProperty("YIELDRATIO",0.8);

  GlassB->SetMaterialPropertiesTable(myMPT1);

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

1 Feedback: Re: Scintillation process not invoked.   (Gumplinger Peter - 12 Oct, 2012)
 Add Message Add Message
to: "Scintillation process not invoked."

 Subscribe Subscribe

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