Message: Issues with WLS initialization Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Issues with WLS initialization 

Forum: Processes Involving Optical Photons
Date: 05 Feb, 2015
From: Chris Ferraioli <Chris Ferraioli>

Hello,

I'm a new user of geant trying to get some WLS fibers set up. So far I have built the main part of my detector and have created Cerenkov showers. WLS fibers are in position to collect the light and direct it towards a PMT, but I've been having trouble getting the fibers to direct the optical photons. I've been testing my system only visually, running interactively and checking to see if any of the showering photons are redirected by the fibers.

I followed the geant documentation WLSOpticalPhysics as an example. I have included the relevant code below. Could someone please let me know if I've missed something?

Thanks in advance, Chris

Added to detector construction:

 const G4int nEntries = 9;
 G4double PhotonEnergy[nEntries] = { 6.6*eV, 6.7*eV, 6.8*eV, 6.9*eV,
                                   7.0*eV, 7.1*eV, 7.2*eV, 7.3*eV, 7.4*eV };
 G4double RIndexFiber[nEntries] =
           { 1.60, 1.60, 1.60, 1.60, 1.60, 1.60, 1.60, 1.60, 1.60 };
 G4double AbsFiber[nEntries] =
           {0.1*mm,0.2*mm,0.3*mm,0.4*cm,1.0*cm,10*cm,1.0*m,10.0*m,10.0*m};
 G4double EmissionFiber[nEntries] =
           {0.0, 0.0, 0.0, 0.1, 0.5, 1.0, 5.0, 10.0, 10.0 };

  G4Material* WLSFiber;
  G4MaterialPropertiesTable* MPTFiber = new G4MaterialPropertiesTable();

  MPTFiber->AddProperty("RINDEX",PhotonEnergy,RIndexFiber,nEntries);
  MPTFiber->AddProperty("WLSABSLENGTH",PhotonEnergy,AbsFiber,nEntries);
  MPTFiber->AddProperty("WLSCOMPONENT",PhotonEnergy,EmissionFiber,nEntries);
  MPTFiber->AddConstProperty("WLSTIMECONSTANT", 0.5*ns);

  Pethylene->SetMaterialPropertiesTable(MPTFiber);

Added to my physics list:

theWLSProcess = new G4OpWLS("WLS");
  G4ProcessManager * pManager = 0;
  pManager = G4OpticalPhoton::OpticalPhoton()->GetProcessManager();
  theWLSProcess->UseTimeProfile("delta");
  pManager->AddDiscreteProcess(theWLSProcess);

Physics reference list info:

  QGSP_BERT_EMV * MyQGSP = new QGSP_BERT_EMV;
   MyQGSP->RegisterPhysics(new MyAdditionalPhysics("WLS") );
 runManager-> SetUserInitialization(MyQGSP);

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

1 Feedback: Re: Issues with WLS initialization   (Gumplinger Peter - 12 Mar, 2015)
 Add Message Add Message
to: "Issues with WLS initialization"

 Subscribe Subscribe

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