Message: Re: How to turn off scintillation and Cerenkov? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: How to turn off scintillation and Cerenkov? 

Forum: Processes Involving Optical Photons
Re: Question How to turn off scintillation and Cerenkov? (Kareem Kazkaz)
Re: Feedback Re: How to turn off scintillation and Cerenkov? (Gumplinger Peter)
Date: 29 Jul, 2009
From: Kareem Kazkaz <Kareem Kazkaz>

Hi Peter. Thanks for the response.

This code is indeed called after the run manager has been initialized. And unfortunately, as far as the particle construction goes, I'm using the methods from GEANT4:

G4BosonConstructor pBoConstructor;
pBoConstructor.ConstructParticle();
 
G4LeptonConstructor pLConstructor;
pLConstructor.ConstructParticle();
 
G4MesonConstructor pMConstructor;
pMConstructor.ConstructParticle();
 
G4BaryonConstructor pBaConstructor;
pBaConstructor.ConstructParticle();
 
G4IonConstructor pIConstructor;
pIConstructor.ConstructParticle();
 
G4ShortLivedConstructor pSConstructor;
pSConstructor.ConstructParticle();

If there's a problem with the particle construction, I assume it would be found deeper down it the GEANT4 code than in the simulation I have.

I set the "/process/verbose 3" in the macro file, and I modified my code by adding a "/run/physicsModified" messenger command:

// Turn the optical processes on or off
G4cout << "useOpticalProcesses = " << useOpticalProcesses << G4endl;
if( useOpticalProcesses ) {
    UI->ApplyCommand( "/process/activate Scintillation" );
    UI->ApplyCommand( "/process/activate Cerenkov" );
    UI->ApplyCommand( "/run/physicsModified" );
} else {
    UI->ApplyCommand( "/process/inactivate Scintillation" );
    UI->ApplyCommand( "/process/inactivate Cerenkov" );
    UI->ApplyCommand( "/run/physicsModified" );
}
 
// Finally, run the beamOn command
stringstream command;
command << "/run/beamOn " << numEvents;
UI->ApplyCommand( command.str() );

Here's the (beginning of) the output after the simulation has started:

### Run 0 start.
G4ProcessTable::Insert : arguments are 0 pointer
G4ProcessTable::Insert Process[RadioactiveDecay] Particle[GenericIon]
This Process Manager is registered !!
G4ProcessTable::Insert Process[RadioactiveDecay] Particle[GenericIon]
G4ProcessTable::Insert : arguments are 0 pointer
G4ProcessTable::Insert Process[hLowEIoni] Particle[GenericIon]
This Process Manager is registered !!
G4ProcessTable::Insert Process[msc] Particle[GenericIon]
This Process Manager is registered !!
G4ProcessTable::Insert Process[Transportation] Particle[GenericIon]
This Process Manager is registered !!
G4ProcessTable::Insert Process[Transportation] Particle[GenericIon]
G4ProcessTable::Insert Process[msc] Particle[GenericIon]
G4ProcessTable::Insert Process[hLowEIoni] Particle[GenericIon]
G4ProcessTable::Insert : arguments are 0 pointer
G4ProcessTable::Insert Process[RadioactiveDecay] Particle[GenericIon]
G4ProcessTable::Insert : arguments are 0 pointer
G4ProcessTable::Insert Process[LUXSimMinEkineCuts] Particle[GenericIon]
This Process Manager is registered !!

...and so on and so on and so on, for a bit over 3000 lines (see the attached file).

Do you have any other recommendations to hunt down what's going wrong? Is it starting to sound like a bug somewhere in the guts of GEANT4?

Kareem Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/07/29/12.55-88865-ProcessVerboseOutput.txt

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

1 None: Re: How to turn off scintillation and Cerenkov?   (Kareem Kazkaz - 29 Jul, 2009)
 Add Message Add Message
to: "Re: How to turn off scintillation and Cerenkov?"

 Subscribe Subscribe

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