Message: Problem with absorption length Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problem with absorption length 

Forum: Processes Involving Optical Photons
Date: 22 Apr, 2013
From: Maxim Laletin <Maxim Laletin>

Hello,

I'm developing a simulation of light propagation in lead tungstate crystalls. I apply the following experimental data as "ABSLENGTH" (here goes an abstract of the code):

G4double x[M]= { 300,  310,  320,  330,  340,  350,  360,  370,  380,  390, 
    400,  410,  420,  430,  440,  450,  460,  470,  480,  490,  
    500,  510,  520,  530,  540,  550,  560,  570,  580,  590,  
    600,  610,  620,  630,  640,  650,  660,  670,  680,  690,  
    700,  710,  720,  730,  740,  750,  760,  770,  780,  790,  
    800,  810,  820,  830,  840,  850,  860,  870,  880,  890,  
    900,  910,  920,  930,  940,  950,  960,  970,  980,  990 };

  G4double y[M]= {  2.46,  2.44,  2.43,  2.43,  2.42,  2.42,  2.75,  4.11,  5.94,  8.48, 
    11.64, 15.09, 18.54, 21.98, 25.67, 29.89, 34.53, 39.01, 43.02, 46.70, 
    50.29, 54.14, 56.89, 59.09, 60.78, 62.23, 63.34, 64.83, 66.35, 67.48, 
    68.62, 69.53, 71.07, 72.56, 74.65, 76.77, 79.42, 82.17, 85.69, 89.88,
    94.29, 99.32,104.6, 111.4, 118.7, 127.2, 137.6, 147.1, 158.5, 166.4, 
    176.0, 176.4, 182.6, 185.9, 187.7, 203.9, 223.0, 246.0, 274.5, 310.5, 
    310.5, 310.5, 310.5, 310.5, 310.5, 311.0, 311.0, 311.0, 311.0, 311.0 };

    for(int i=0;i<M;i++){
      x[i]=((4.135666*3.*1.e+02)/x[i])*eV;
      y[i]=y[i]*cm;
    }

  G4MaterialPropertiesTable* myMPT1 = new G4MaterialPropertiesTable();

  myMPT1->AddProperty("RINDEX", PhotonEnergy, RefractiveIndex1,nEntries)
        ->SetSpline(true);
  myMPT1->AddProperty("ABSLENGTH", x, y,     M);
        //->SetSpline(true);

When I run the simulation the observed "path" of the photon with the appropriate energy doesn't correspond to the mentioned ABSLENGTH values. But if I set a single value for any energy everything works fine.

Where could have I make a mistake or am I doing something wrong?

Thank you for your answers!

Maxim

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

1 Feedback: Re: Problem with absorption length   (Gumplinger Peter - 23 Apr, 2013)
(_ Feedback: Re: Problem with absorption length   (Maxim Laletin - 06 May, 2013)
(_ Feedback: Re: Problem with absorption length   (Gumplinger Peter - 06 May, 2013)
 Add Message Add Message
to: "Problem with absorption length"

 Subscribe Subscribe

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