Message: Re: program stops when scintillation occurs :( Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Disagree Re: program stops when scintillation occurs :( 

Forum: Processes Involving Optical Photons
Re: Question program stops when scintillation occurs :( (Yen-Yung Chang)
Re: Feedback Re: program stops when scintillation occurs :( (Gumplinger Peter)
Re: None Re: program stops when scintillation occurs :( (Yen-Yung Chang)
Date: 15 Oct, 2012
From: Gumplinger Peter <Gumplinger Peter>

> SetProcessOrderingToLast(CerenkovPross, idxAtRest);
> 
> for Cerenkov process because exampleN06 has this line...

But sorry, exampleN06 has:

    if (theCerenkovProcess->IsApplicable(*particle)) {
      pmanager->AddProcess(theCerenkovProcess);
      pmanager->SetProcessOrdering(theCerenkovProcess,idxPostStep);
    }
    if (theScintillationProcess->IsApplicable(*particle)) {
      pmanager->AddProcess(theScintillationProcess);
      pmanager->SetProcessOrderingToLast(theScintillationProcess, idxAtRest);
      pmanager->SetProcessOrderingToLast(theScintillationProcess, idxPostStep);
    }

I see part of the problem:

You attach ScintXtal_MPT to ScintXtal and then:

  YSO->AddMaterial(ScintXtal, 99.8*perCent);

  crystalCell_log = new G4LogicalVolume(crystalCell_box,YSO,....

This doesn't work. The MaterialPropertiesTable has to be directly attached to the material that's making up the G4LogicalVolume. If the 0.2perCent of inert material in the scintillator matters, then you'll have to reduce the yield by that amount to account for it before filling the table.

Also, the arrays of properties are in energy of the optical photon (eV not keV).

 Add Message Add Message
to: "Re: program stops when scintillation occurs :("

 Subscribe Subscribe

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