Message: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD 

Forum: Processes Involving Optical Photons
Re: Question G4Exception-START: ELECTRONSCINTILLATIONYIELD (Sergio Lo Meo)
Re: Feedback Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD (Gumplinger Peter)
Re: Question Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD (Sergio Lo Meo)
Re: Feedback Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD (Gumplinger Peter)
Date: 13 May, 2016
From: Sergio Lo Meo <Sergio Lo Meo>

Dear Peter,

thank you for the answer. I modified my code (OpticalPhysics.cc) as written in Application Manual (pdf) pag 208, example 5.6.

Now, I have these lines in the ConstructProcess of my OpticalPhysics.cc

-----------------------------------------------------------------------
 G4Scintillation* ScintillationProcess = new G4Scintillation();
  ScintillationProcess->SetScintillationYieldFactor(fYieldFactor);
  ScintillationProcess->SetScintillationExcitationRatio(fExcitationRatio);
  ScintillationProcess->SetFiniteRiseTime(fFiniteRiseTime);
  ScintillationProcess->SetScintillationByParticleType(fScintillationByParticleType);
  ScintillationProcess->SetTrackSecondariesFirst(fProcessTrackSecondariesFirst[kScintillation]);
  G4EmSaturation* emSaturation = G4LossTableManager::Instance()->EmSaturation();
  ScintillationProcess->AddSaturation(emSaturation);
  UIhelpers::buildCommands(ScintillationProcess);
  OpProcesses[kScintillation] = ScintillationProcess;

  // new lines

  G4Scintillation* theProtonScintProcess = new G4Scintillation("Scintillation");

  theProtonScintProcess->SetTrackSecondariesFirst(true);
  theProtonScintProcess->SetScintillationYieldFactor(0.5);

  aParticleIterator->reset();
  while( (*aParticleIterator)() ){
    G4ParticleDefinition* particle = aParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();
    if (theProtonScintProcess->IsApplicable(*particle)) {
       if (particleName == "proton") {
          pmanager->AddProcess(theProtonScintProcess);
          pmanager->SetProcessOrderingToLast(theProtonScintProcess, idxAtRest);
          pmanager->SetProcessOrderingToLast(theProtonScintProcess, idxPostStep);
       }
    }
  }

  // end new lines
---------------------------------------------------------------------

cmake ok, make -j2 ok, but when I gun a 500 MeV proton, I see two warnings:

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : ProcMan114
      issued by : G4ProcessManager::SetProcessOrderingToLast()
Set Ordering Last is invoked twice for Scintillation to proton
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : ProcMan114
      issued by : G4ProcessManager::SetProcessOrderingToLast()
Set Ordering Last is invoked twice for Scintillation to proton
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

Are they a problem? I think so... see below

Where is my mistake?

Then I make another test (same macro and same seed):

The ScintillationYield of my Scintillator is 10000./MeV optical photons. fScintillationByParticleType(false)

Always 500 MeV proton.

First Test) no lines for proton in my OpticalPhysics.cc (/* G4Scintillation* theProtonScintProcess = new G4Scintillation("Scintillation"); ....   */)
 In my PMT (QE = 100%) I counted 3594
optical photons (  1.138 MeV deposited in my Scint)

Second Test) with lines for proton In my PMT (QE = 100%) I counted 5513 optical photons ( 1.27 MeV deposited in my Scint)

It is impossible... It seems that the line: theProtonScintProcess->SetScintillationYieldFactor(0.5); of the first test is completely ignored...

I wish that the number of optical photons in the first test was about two times that the second test..

Thanks again Best Regards

Sergio Lo Meo

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

1 Feedback: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Gumplinger Peter - 13 May, 2016)
(_ Ok: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Sergio Lo Meo - 22 May, 2016)
(_ Question: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Nicholas DePorzio - 09 Aug, 2016)
(_ Feedback: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Gumplinger Peter - 11 Aug, 2016)
 Add Message Add Message
to: "Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD"

 Subscribe Subscribe

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