Message: Refraction has not occur Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Refraction has not occur  

Forum: Processes Involving Optical Photons
Date: 18 Jan, 2011
From: Davinder Siwal <Davinder Siwal>

Hiii Geant4 Users

I am trying to simulate LYSO crystaland i have set up the material properties and physics list for gammas and optical photons as

void PlanarPhysicsList::ConstructOpticalProcess()
{
  G4Scintillation* theScintillationProcess = new G4Scintillation("Scintillation");
  G4OpAbsorption* theAbsorptionProcess = new G4OpAbsorption();
  G4OpRayleigh* theRayleighScatteringProcess = new G4OpRayleigh();
   G4OpMieHG* theMieHGScatteringProcess    = new G4OpMieHG();
  G4OpBoundaryProcess* theBoundaryProcess  = new G4OpBoundaryProcess();

  //-------- Set Verbosity 
  theScintillationProcess->SetVerboseLevel(1);
  theAbsorptionProcess->SetVerboseLevel(1);
  theRayleighScatteringProcess->SetVerboseLevel(1);
  theBoundaryProcess->SetVerboseLevel(1);
  //-----------------------

  // G4ProcessManager* pManager = 0;
  //pManager = G4OpticalPhoton::OpticalPhoton()->GetProcessManager();
  //pManager->AddDiscreteProcess(theAbsorptionProcess);
  //pManager->AddDiscreteProcess(theRayleighScatteringProcess);
  //pManager->AddDiscreteProcess(theBoundaryProcess);

  theScintillationProcess->SetScintillationYieldFactor(1);
  theScintillationProcess->SetTrackSecondariesFirst(true);

  G4OpticalSurfaceModel themodel = unified;
  theBoundaryProcess->SetModel(themodel);

  theParticleIterator->reset();
  while( (*theParticleIterator)())
    {
      G4ParticleDefinition* particle = theParticleIterator->value();
      G4ProcessManager* pmanager = particle->GetProcessManager();
      G4String particleName = particle->GetParticleName();

      if (theScintillationProcess->IsApplicable(*particle)) 
	{
	  pmanager->AddProcess(theScintillationProcess);
	  pmanager->SetProcessOrderingToLast(theScintillationProcess, idxAtRest);
	  pmanager->SetProcessOrderingToLast(theScintillationProcess, idxPostStep);
	}
      if (particleName == "opticalphoton") 
	{
	  G4cout << " AddDiscreteProcess to OpticalPhoton " << G4endl;
	  pmanager->AddDiscreteProcess(theAbsorptionProcess);
	  pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
	  pmanager->AddDiscreteProcess(theBoundaryProcess);
	  pmanager->AddDiscreteProcess(theMieHGScatteringProcess);
	  pmanager->AddDiscreteProcess(theBoundaryProcess);
	}

    }
}

and my material properties are

//-----------------Detector Scintillation properties----//

 const G4int nEntries = 2;
 G4double PhotonEnergy[nEntries] ={2.85,2.90}; // actual value 2,89@428nm

 G4double RefractiveIndex[nEntries] = {1.82,1.82};

 G4double AbsLength[nEntries] = {1.10*cm,1.11*cm};

 G4double DecayTime[nEntries] = {40*ns};

 G4MaterialPropertiesTable* myDet = new G4MaterialPropertiesTable();
 myDet->AddProperty("FASTCOMPONENT", PhotonEnergy, DecayTime, nEntries);
 myDet->AddProperty("RINDEX", PhotonEnergy, RefractiveIndex, nEntries);
 myDet->AddProperty("ABSLENGTH", PhotonEnergy, AbsLength, nEntries);

 myDet->AddConstProperty("SCINTILLATIONYIELD",26000./MeV);
 myDet->AddConstProperty("RESOLUTIONSCALE",1.0);
 myDet->AddConstProperty("FASTTIMECONSTANT",40*ns);

 LYSO->SetMaterialPropertiesTable(myDet);

 //------------------------------------------------------------------//

 //
 // Air
 //

 G4double RefractiveIndex2[nEntries] = {1.00,1.00};
 G4MaterialPropertiesTable* myMPT1 = new G4MaterialPropertiesTable();
 myMPT1->AddProperty("RINDX",PhotonEnergy,RefractiveIndex2, nEntries);
 Air->SetMaterialPropertiesTable(myMPT1);

 //------------------------------------------------------------------//

 //------------- Detector Surface --------//

 G4OpticalSurface* DetOptSurface = new G4OpticalSurface("DetectorSurface");
 DetOptSurface->SetType(dielectric_dielectric);
 DetOptSurface->SetFinish(ground);
 DetOptSurface->SetModel(unified);

  new G4LogicalBorderSurface("Surface",hpge_phys,experimentalHall_phys,DetOptSurface);

my problem is that there is no reflection occurs at the bounday of the the detector can anybody tell me what is going to be wrong.

thanks in advance

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

1 Feedback: Re: Refraction has not occur   (Gumplinger Peter - 18 Jan, 2011)
 Add Message Add Message
to: "Refraction has not occur "

 Subscribe Subscribe

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