Message: Re: Transmission of Optical Photons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Transmission of Optical Photons 

Forum: Processes Involving Optical Photons
Re: None Transmission of Optical Photons
Re: Feedback Re: Transmission of Optical Photons (Gumplinger Peter)
Re: Feedback Re: Transmission of Optical Photons
Re: Feedback Re: Transmission of Optical Photons (Gumplinger Peter)
Date: 24 Jul, 2009
From: <jdswado@gatech.edu>

Hello again,

So I found an interesting thing. The "No Refractive Index" problem comes up only when one of the materials involved is my material "Aluminum". Currently my materials property table for aluminum looks like this:

// Aluminum Material Properties Table Stuff

    G4double RealRefractive[nEntries] =
            { 0.444,  0.467,  0.49,   0.5156, 0.5412, 0.5668, 0.5924,
              0.618,  0.6482, 0.6784, 0.7086, 0.7388, 0.769,  0.8068,
              0.8446, 0.8824, 0.9202, 0.958,  1.0064, 1.0548, 1.1032,
              1.1516, 1.2,    1.263,  1.326,  1.389,  1.452,  1.515,
              1.578,  1.641,  1.704,  1.767,  1.83,   1.927,  2.024,
	          2.121,  2.218,  2.315};

    G4double ImaginaryRefractive[nEntries] =
	    {-4.612, -4.736, -4.86,  -4.982, -5.104,  -5.226, -5.348,
	     -5.47,  -5.592, -5.714, -5.836, -5.958,  -6.08,  -6.202,
	     -6.324, -6.446, -6.568, -6.69,  -6.804,  -6.918, -7.032,
	     -7.146, -7.26,  -7.365, -7.47,  -7.575,  -7.68,  -7.785,
	     -7.89,  -7.995, -8.1,   -8.205, -8.31,   -8.324, -8.338,
	     -8.352, -8.366, -8.38}; 

	G4double ALAbsorption[nEntries] =      
	        {0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm,
                 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm,
                 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm,
                 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm,
                 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm, 0.1*mm,
                 0.1*mm, 0.1*mm, 0.1*mm};

  G4MaterialPropertiesTable* myMPT3 = new G4MaterialPropertiesTable();
  myMPT3->AddProperty("REALRINDEX", PhotonEnergy, RealRefractive, nEntries);
myMPT3->AddProperty("IMAGINARYRINDEX",PhotonEnergy,ImaginaryRefractive,nEntries);
  myMPT3->AddProperty("ABSLENGTH", PhotonEnergy, ALAbsorption,nEntries);

////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////

As I said, using this I get a great deal of "No Refractive Index". However, if I remove the IMAGINARYRINDEX and replace REALRINDEX with RINDEX, the program runs with no instances of "No Refractive Index". Any idea why this could be?

Also note that I have a surface defined from the aluminum but it currently doesn't seem to be affecting anything at all.

Aluminum surface definition:

G4OpticalSurface* OpALSurface = new G4OpticalSurface("ALSurface");
  OpALSurface->SetType(dielectric_metal);// should be dielectric_metal
  OpALSurface->SetFinish(polished);
  OpALSurface->SetModel(glisur);
  G4double ALRIndex[nEntries]=            
            { 0.444,  0.467,  0.49,   0.5156, 0.5412, 0.5668, 0.5924,
              0.618,  0.6482, 0.6784, 0.7086, 0.7388, 0.769,  0.8068,
              0.8446, 0.8824, 0.9202, 0.958,  1.0064, 1.0548, 1.1032,
              1.1516, 1.2,    1.263,  1.326,  1.389,  1.452,  1.515,
              1.578,  1.641,  1.704,  1.767,  1.83,   1.927,  2.024,
	      2.121,  2.218,  2.315}; 

G4double ALImgRIndex[nEntries]= 	       
            {-4.612, -4.736, -4.86,  -4.982, -5.104,  -5.226, -5.348,
	     -5.47,  -5.592, -5.714, -5.836, -5.958,  -6.08,  -6.202,
	     -6.324, -6.446, -6.568, -6.69,  -6.804,  -6.918, -7.032,
	     -7.146, -7.26,  -7.365, -7.47,  -7.575,  -7.68,  -7.785,
	     -7.89,  -7.995, -8.1,   -8.205, -8.31,   -8.324, -8.338,
	     -8.352, -8.366, -8.38};  

 G4MaterialPropertiesTable* myST_AL = new G4MaterialPropertiesTable();

  myST3->AddProperty("REALRINDEX", PhotonEnergy, ALRIndex, nEntries);  
  myST3->AddProperty("IMAGINARYRINDEX", PhotonEnergy, ALImgRIndex, nEntries);  

  OpCrystalSurfaceTwo->SetMaterialPropertiesTable(myST3);
  G4LogicalBorderSurface* ALSurface =new G4LogicalBorderSurface
("TubeTwoSurface",TubeTwo_phys,experimentalHall_phys,OpCrystalSurfaceTwo);

Have I implemented the imaginary R Index incorrectly?

Thanks!

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

1 Feedback: Re: Transmission of Optical Photons   (Gumplinger Peter - 24 Jul, 2009)
 Add Message Add Message
to: "Re: Transmission of Optical Photons"

 Subscribe Subscribe

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