Message: Re: Problems for transportation of high energy optical photons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Problems for transportation of high energy optical photons 

Keywords: optical photon transportation OpAbsorption
Forum: Processes Involving Optical Photons
Date: 17 Mar, 2008
From: Keito Horie <kate@rcnp.osaka-u.ac.jp>

Dear Peter Gumplinger,

Thank you very much for your prompt response. I agree that something is wrong in the setup, as you guess. I am sending a typical example which is the file for PbW0 crystals. I am looking forward to your answer. Thank you very much in advance.

Best regards, Keito Horie

------------
 G4Element* O = new G4Element("Oxygen"  , "O", z=8 , a=16.00*g/mole);
 G4Element* Pb = new G4Element("Lead", "Pb", z=82 , a=207.20*g/mole);
 G4Element* W = new G4Element("Tungsten", "W", z=74 , a=183.84*g/mole);
 G4Material* PWO =
   new G4Material("PWO",density= 8.28*g/cm3, nelements=3);
 PWO->AddElement(Pb, natoms=1);
 PWO->AddElement(W, natoms=1);
 PWO->AddElement(O, natoms=4);

//
// PWO property
//
 const G4int nEntries = 32;
 G4double PhotonEnergy[nEntries] =
           { 2.034*eV, 2.068*eV, 2.103*eV, 2.139*eV,
             2.177*eV, 2.216*eV, 2.256*eV, 2.298*eV,
             2.341*eV, 2.386*eV, 2.433*eV, 2.481*eV,
             2.532*eV, 2.585*eV, 2.640*eV, 2.697*eV,
             2.757*eV, 2.820*eV, 2.885*eV, 2.954*eV,
             3.026*eV, 3.102*eV, 3.181*eV, 3.265*eV,
             3.353*eV, 3.446*eV, 3.545*eV, 3.649*eV,
             3.760*eV, 3.877*eV, 4.002*eV, 4.136*eV };
 G4double RefractiveIndex5[nEntries] =
           { 2.16, 2.16, 2.16, 2.16, 2.16, 2.16, 2.16,
             2.16, 2.16, 2.16, 2.16, 2.16, 2.16, 2.16,
             2.16, 2.16, 2.16, 2.16, 2.16, 2.16, 2.16,
             2.16, 2.16, 2.16, 2.16, 2.16, 2.16, 2.16,
             2.16, 2.16, 2.16, 2.16 };
 G4double Absorption4[nEntries] =
          {4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,
           4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,
           4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,
           4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,  4.20*m,
           4.20*m,  4.20*m };
 G4double ScintilFast[nEntries] =
            { 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00 };

 G4MaterialPropertiesTable* myMPT5 = new G4MaterialPropertiesTable();
 myMPT5->AddProperty("RINDEX", PhotonEnergy, RefractiveIndex5, nEntries);
 myMPT5->AddProperty("ABSLENGTH",    PhotonEnergy, Absorption4, nEntries);
 myMPT5->AddProperty("FASTCOMPONENT",PhotonEnergy, ScintilFast,     nEntries);
 myMPT5->AddProperty("SLOWCOMPONENT",PhotonEnergy, ScintilFast,     nEntries);

 myMPT5->AddConstProperty("SCINTILLATIONYIELD",500./MeV);
 myMPT5->AddConstProperty("RESOLUTIONSCALE",1.0);
 myMPT5->AddConstProperty("FASTTIMECONSTANT", 10.*ns);
 myMPT5->AddConstProperty("SLOWTIMECONSTANT",1000.*ns);
 myMPT5->AddConstProperty("YIELDRATIO",1.);

PWO->SetMaterialPropertiesTable(myMPT5);

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

1 Disagree: Re: Problems for transportation of high energy optical photons   (Peter Gumplinger - 17 Mar, 2008)
(_ None: Re: Problems for transportation of high energy optical photons (solved.)   (Keito Horie - 18 Mar, 2008)
 Add Message Add Message
to: "Re: Problems for transportation of high energy optical photons"

 Subscribe Subscribe

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