{
  double m = 1, cm = 0.01, nanometer = 0.000000001, GeV = 1000000000;
  double LambdaE = 2.0 * 3.14159265358979323846 * 1.973269602e-16 * m * GeV;
  const int num = 36;
  double WaveLength[num];
  double Absorption[num]; // default value for absorption
  double AirAbsorption[num]; // absorption value for air
  double AirRefractiveIndex[num]; // air refractive index
  double PhotonEnergy[num]; // energy of photons which correspond to the given 

  {
    // absorption of KamLandOil per 50 cm - from jjv
    double KamLandOilAbsorption[num]=
      {0.97469022,0.976603956,0.978511548,0.980400538,0.982258449,0.984072792,
       0.985831062,0.987520743,0.989129303,0.990644203,0.992052894,
       0.993342822,0.994501428,0.995516151,0.996374433,0.997063719,
       0.997571464,0.997885132,0.997992205,0.997880183,0.997536591,
       0.99,0.98,0.97,0.96,0.94,0.93,0.924507,0.89982,0.883299,
       0.85657,0.842637,0.77020213,0.65727,0.324022,0.019192};

    // absorption of quartz per 1 m - from jjv
    double QuartzAbsorption[num] = 
      {0.999572036,0.999544661,0.999515062,0.999483019,0.999448285,
       0.999410586,0.999369611,0.999325013,0.999276402,0.999223336,
       0.999165317,0.999101778,0.999032079,0.998955488,0.998871172,
       0.998778177,0.99867541,0.998561611,0.998435332,0.998294892,0.998138345,
       0.997963425,0.997767484,0.997547418,
       0.99729958,0.99701966,0.99670255,0.996342167,0.995931242,0.995461041,
       0.994921022,0.994298396,0.993577567,0.992739402,0.991760297,0.990610945};
  
    // absorption of epotek per one layer - thicknes 0.001'' - from jjv
    double EpotekAbsorption[num] = 
      {0.99999999,0.99999999,0.99999999,0.99999999,0.99999999,
       0.99999999,0.99999999,0.99999999,0.99999999,0.99999999,
       0.99999999,0.99999999,0.99999999,0.99999999,0.99999999,
       0.99999999,0.99999999,0.99999999,0.99999999,0.99999999,
       0.99999999,0.99999999,0.99999999,0.99999999,0.99999999,
       0.9999,0.9998,0.9995,0.999,0.998,0.997,0.996,0.9955,0.993,
       0.9871,0.9745};

        double Nlak33aAbsorption2[76]={371813,352095,331021,310814,291458,272937,255238,238342,222234,206897,192313,178463,165331,152896,141140,130043,119585,109747,100507,91846.3,83743.1,76176.7,69126.1,62570.2,56488,50858.3,45660.1,40872.4,36474.6,32445.8,28765.9,25414.6,22372.2,19619.3,17136.9,14906.5,12910.2,11130.3,9550.13,8153.3,6924.25,5848.04,4910.46,4098.04,3398.06,2798.54,2288.32,1856.99,1494.92,1193.28,943.973,739.657,573.715,440.228,333.94,250.229,185.064,134.967,96.9664,68.5529,47.6343,32.4882,21.7174,14.2056,9.07612,5.65267,3.4241,2.01226,1.14403,0.62722,0.330414,0.166558,0.0799649,0.0363677,0.0155708,0.00623089};

    double Nlak33aAbsorption[753]={371813,371813,369586,367369,365160,362961,360770,358588,356415,354250,352095,349948,347810,345681,343560,341449,339346,337252,335166,333089,331021,328962,326911,324869,322835,320810,318794,316786,314787,312796,310814,308840,306875,304918,302970,301030,299099,297176,295261,293355,291458,289568,287687,285814,283950,282094,280246,278407,276575,274752,272937,271131,269332,267542,265760,263986,262220,260462,258713,256971,255238,253512,251795,250085,248384,246690,245005,243327,241658,239996,238342,236696,235058,233428,231806,230191,228584,226985,225394,223810,222234,220666,219106,217553,216008,214471,212941,211418,209904,208397,206897,205405,203921,202444,200974,199512,198057,196610,195170,193738,192313,190895,189485,188082,186686,185297,183916,182542,181175,179816,178463,177118,175780,174449,173125,171809,170499,169196,167901,166612,165331,164056,162789,161528,160274,159027,157787,156554,155328,154108,152896,151690,150491,149299,148113,146934,145762,144596,143437,142285,141140,140001,138868,137742,136623,135510,134404,133304,132210,131123,130043,128968,127901,126839,125784,124735,123693,122656,121626,120602,119585,118574,117568,116569,115576,114590,113609,112634,111666,110703,109747,108796,107851,106913,105980,105053,104132,103217,102308,101405,100507,99615.5,98729.4,97849.2,96974.6,96105.8,95242.6,94385.1,93533.2,92687,91846.3,91011.3,90181.7,89357.8,88539.3,87726.3,86918.8,86116.8,85320.1,84528.9,83743.1,82962.7,82187.6,81417.8,80653.3,79894.1,79140.2,78391.5,77648,76909.8,76176.7,75448.8,74726,74008.3,73295.7,72588.2,71885.8,71188.4,70496,69808.5,69126.1,68448.6,67776,67108.3,66445.6,65787.6,65134.5,64486.3,63842.8,63204.1,62570.2,61941,61316.6,60696.8,60081.7,59471.2,58865.4,58264.2,57667.6,57075.5,56488,55905,55326.5,54752.5,54183,53617.9,53057.2,52500.9,51949,51401.5,50858.3,50319.4,49784.8,49254.5,48728.4,48206.6,47688.9,47175.5,46666.2,46161.1,45660.1,45163.2,44670.4,44181.6,43696.9,43216.3,42739.6,42266.9,41798.1,41333.3,40872.4,40415.5,39962.4,39513.1,39067.7,38626.1,38188.3,37754.2,37324,36897.4,36474.6,36055.4,35639.9,35228.1,34819.9,34415.3,34014.3,33616.9,33223,32832.7,32445.8,32062.5,31682.6,31306.1,30933.1,30563.5,30197.3,29834.4,29474.9,29118.7,28765.9,28416.3,28069.9,27726.9,27387,27050.4,26716.9,26386.6,26059.5,25735.5,25414.6,25096.8,24782,24470.3,24161.6,23856,23553.3,23253.6,22956.9,22663.1,22372.2,22084.2,21799,21516.8,21237.3,20960.7,20686.9,20415.9,20147.6,19882.1,19619.3,19359.2,19101.8,18847,18594.9,18345.4,18098.6,17854.3,17612.6,17373.5,17136.9,16902.8,16671.3,16442.2,16215.5,15991.4,15769.6,15550.3,15333.3,15118.8,14906.5,14696.7,14489.1,14283.9,14080.9,13880.2,13681.8,13485.6,13291.6,13099.8,12910.2,12722.7,12537.4,12354.3,12173.2,11994.3,11817.4,11642.6,11469.8,11299.1,11130.3,10963.6,10798.9,10636.1,10475.2,10316.3,10159.3,10004.2,9850.99,9699.64,9550.13,9402.46,9256.6,9112.55,8970.28,8829.79,8691.04,8554.04,8418.76,8285.19,8153.3,8023.1,7894.56,7767.66,7642.4,7518.75,7396.71,7276.25,7157.36,7040.04,6924.25,6809.99,6697.24,6586,6476.23,6367.94,6261.11,6155.71,6051.75,5949.19,5848.04,5748.27,5649.87,5552.84,5457.14,5362.78,5269.74,5178,5087.54,4998.37,4910.46,4823.8,4738.38,4654.19,4571.21,4489.42,4408.82,4329.4,4251.14,4174.02,4098.04,4023.19,3949.45,3876.8,3805.24,3734.76,3665.34,3596.97,3529.64,3463.34,3398.06,3333.77,3270.48,3208.18,3146.84,3086.46,3027.02,2968.53,2910.95,2854.3,2798.54,2743.68,2689.7,2636.59,2584.34,2532.94,2482.38,2432.64,2383.73,2335.63,2288.32,2241.8,2196.06,2151.09,2106.88,2063.41,2020.69,1978.69,1937.42,1896.85,1856.99,1817.82,1779.33,1741.51,1704.36,1667.87,1632.02,1596.81,1562.23,1528.27,1494.92,1462.18,1430.03,1398.47,1367.49,1337.07,1307.22,1277.93,1249.18,1220.96,1193.28,1166.12,1139.48,1113.34,1087.71,1062.56,1037.9,1013.72,990.011,966.764,943.973,921.632,899.733,878.269,857.234,836.621,816.424,796.635,777.249,758.258,739.657,721.439,703.598,686.128,669.022,652.275,635.881,619.834,604.128,588.757,573.715,558.998,544.599,530.513,516.735,503.258,490.079,477.192,464.591,452.271,440.228,428.456,416.951,405.708,394.721,383.987,373.501,363.257,353.252,343.481,333.94,324.624,315.529,306.651,297.985,289.528,281.275,273.223,265.367,257.703,250.229,242.94,235.831,228.901,222.144,215.558,209.139,202.884,196.788,190.85,185.064,179.429,173.942,168.598,163.395,158.329,153.399,148.6,143.931,139.387,134.967,130.668,126.486,122.42,118.466,114.622,110.886,107.254,103.726,100.297,96.9664,93.7312,90.5894,87.5385,84.5765,81.7012,78.9105,76.2023,73.5747,71.0255,68.5529,66.1549,63.8297,61.5755,59.3903,57.2725,55.2203,53.2321,51.306,49.4406,47.6343,45.8853,44.1924,42.5538,40.9682,39.4341,37.9501,36.5149,35.127,33.7852,32.4882,31.2347,30.0235,28.8534,27.7232,26.6318,25.5779,24.5607,23.5788,22.6314,21.7174,20.8358,19.9856,19.1658,18.3757,17.6141,16.8803,16.1734,15.4925,14.8368,14.2056,13.5981,13.0134,12.451,11.9099,11.3896,10.8894,10.4086,9.9465,9.50255,9.07612,8.66662,8.27347,7.8961,7.53398,7.18656,6.85335,6.53383,6.22752,5.93396,5.65267,5.38323,5.1252,4.87815,4.64169,4.41543,4.19898,3.99197,3.79405,3.60487,3.4241,3.2514,3.08648,2.92902,2.77872,2.63532,2.49852,2.36807,2.2437,2.12518,2.01226,1.90471,1.8023,1.70483,1.61207,1.52385,1.43995,1.36019,1.2844,1.2124,1.14403,1.07912,1.01752,0.959082,0.903663,0.851125,0.801337,0.754172,0.709508,0.667228,0.62722,0.589375,0.553591,0.519768,0.487811,0.457629,0.429134,0.402244,0.376878,0.352959,0.330414,0.309173,0.289169,0.270338,0.252619,0.235953,0.220284,0.20556,0.191729,0.178744,0.166558,0.155127,0.14441,0.134366,0.124958,0.11615,0.107907,0.100197,0.0929896,0.0862548,0.0799649,0.0740936,0.0686159,0.0635081,0.0587477,0.0543134,0.0501852,0.046344,0.0427718,0.0394517,0.0363677,0.0335045,0.0308479,0.0283844,0.0261015,0.023987,0.0220297,0.0202191,0.0185452,0.0169987,0.0155708,0.0142532,0.0130381,0.0119184,0.0108873,0.00993835,0.00906562,0.00826355,0.00752692,0.00685088,0.00623089,0.0056627,0.01
    };
    

    double EpotekThickness = 0.001*2.54*cm;
    for(int i=0;i<num;i++){
      WaveLength[i]= (300 +i*10)*nanometer;
      Absorption[i]= 100*m;
      AirAbsorption[i] = 4.*cm; 
      AirRefractiveIndex[i] = 1.; 
      PhotonEnergy[num-(i+1)]= LambdaE/WaveLength[i];
      
      EpotekAbsorption[i] = (-1)/log(EpotekAbsorption[i])*EpotekThickness;
      QuartzAbsorption[i] = (-1)/log(QuartzAbsorption[i])*100*cm;
      KamLandOilAbsorption[i] = (-1)/log(KamLandOilAbsorption[i])*50*cm;
    }

  }

  {
    double QuartzRefractiveIndex[num]=
      {1.456535,1.456812,1.4571,1.457399,1.457712,1.458038,1.458378,
       1.458735,1.459108,1.4595,1.459911,1.460344,1.460799,1.46128,
       1.461789,1.462326,1.462897,1.463502,1.464146,1.464833,
       1.465566,1.46635,1.46719,1.468094,1.469066,1.470116,1.471252,1.472485,
       1.473826,1.475289,1.476891,1.478651,1.480592,1.482739,1.485127,1.487793};

    double EpotekRefractiveIndex[num]=
      {1.554034,1.555575,1.55698,1.558266,1.559454,1.56056,1.561604,
       1.562604,1.563579,1.564547,1.565526,1.566536,1.567595,
       1.568721,1.569933,1.57125,1.57269,1.574271,1.576012,
       1.577932,1.580049,1.582381,1.584948,1.587768,1.590859,
       1.59424,1.597929,1.601946,1.606307,1.611033,1.616141,1.621651,1.62758,
       1.633947,1.640771,1.64807};

  

    double KamLandOilRefractiveIndex[num]={
      1.433055,1.433369,1.433698,1.434045,1.434409,1.434793,1.435198,
      1.435626,1.436077,1.436555,1.4371,1.4376,1.4382,1.4388,1.4395,
      1.4402,1.4409,1.4415,1.4425,1.4434,1.4444,1.4455,1.4464,1.4479,1.4501,
      1.450428,1.451976,1.453666,1.455513,1.45754,1.45977,1.462231,1.464958,
      1.467991,1.471377,1.475174};


    double enNlak33a2[76] = {1,1.2511,1.26386,1.27687,1.29016,1.30372,1.31758,1.33173,1.34619,1.36097,1.37607,1.39152,1.40731,1.42347,1.44,1.45692,1.47425,1.49199,1.51016,1.52878,1.54787,1.56744,1.58751,1.6081,1.62923,1.65092,1.6732,1.69609,1.71961,1.7438,1.76868,1.79427,1.82062,1.84775,1.87571,1.90452,1.93423,1.96488,1.99652,2.0292,2.06296,2.09787,2.13398,2.17135,2.21006,2.25017,2.29176,2.33492,2.37973,2.42631,2.47473,2.52514,2.57763,2.63236,2.68946,2.7491,2.81143,2.87666,2.94499,3.01665,3.09187,3.17095,3.25418,3.34189,3.43446,3.53231,3.6359,3.74575,3.86244,3.98663,4.11908,4.26062,4.41225,4.57506,4.75035,4.93961};

    double enNlak33a[753]={1,1.23984,1.24108,1.24233,1.24357,1.24482,1.24607,1.24733,1.24858,1.24984,1.2511,1.25237,1.25363,1.2549,1.25617,1.25745,1.25872,1.26,1.26128,1.26257,1.26386,1.26514,1.26644,1.26773,1.26903,1.27033,1.27163,1.27294,1.27425,1.27556,1.27687,1.27819,1.27951,1.28083,1.28215,1.28348,1.28481,1.28614,1.28748,1.28882,1.29016,1.2915,1.29285,1.2942,1.29555,1.29691,1.29826,1.29962,1.30099,1.30235,1.30372,1.3051,1.30647,1.30785,1.30923,1.31062,1.312,1.31339,1.31478,1.31618,1.31758,1.31898,1.32039,1.32179,1.3232,1.32462,1.32603,1.32745,1.32888,1.3303,1.33173,1.33316,1.3346,1.33604,1.33748,1.33892,1.34037,1.34182,1.34327,1.34473,1.34619,1.34765,1.34912,1.35059,1.35206,1.35354,1.35502,1.3565,1.35799,1.35948,1.36097,1.36246,1.36396,1.36546,1.36697,1.36848,1.36999,1.37151,1.37303,1.37455,1.37607,1.3776,1.37913,1.38067,1.38221,1.38375,1.3853,1.38685,1.3884,1.38996,1.39152,1.39308,1.39465,1.39622,1.39779,1.39937,1.40095,1.40254,1.40412,1.40572,1.40731,1.40891,1.41051,1.41212,1.41373,1.41534,1.41696,1.41858,1.42021,1.42184,1.42347,1.42511,1.42675,1.42839,1.43004,1.43169,1.43334,1.435,1.43666,1.43833,1.44,1.44168,1.44335,1.44504,1.44672,1.44841,1.45011,1.45181,1.45351,1.45521,1.45692,1.45864,1.46036,1.46208,1.4638,1.46553,1.46727,1.46901,1.47075,1.4725,1.47425,1.476,1.47776,1.47952,1.48129,1.48306,1.48484,1.48662,1.48841,1.49019,1.49199,1.49379,1.49559,1.49739,1.4992,1.50102,1.50284,1.50466,1.50649,1.50832,1.51016,1.512,1.51385,1.5157,1.51755,1.51941,1.52128,1.52315,1.52502,1.5269,1.52878,1.53067,1.53256,1.53446,1.53636,1.53827,1.54018,1.54209,1.54401,1.54594,1.54787,1.5498,1.55174,1.55369,1.55564,1.55759,1.55955,1.56151,1.56348,1.56546,1.56744,1.56942,1.57141,1.5734,1.5754,1.57741,1.57942,1.58143,1.58345,1.58548,1.58751,1.58954,1.59158,1.59363,1.59568,1.59773,1.5998,1.60186,1.60394,1.60601,1.6081,1.61018,1.61228,1.61438,1.61648,1.61859,1.62071,1.62283,1.62496,1.62709,1.62923,1.63137,1.63352,1.63568,1.63784,1.64,1.64217,1.64435,1.64654,1.64873,1.65092,1.65312,1.65533,1.65754,1.65976,1.66199,1.66422,1.66645,1.6687,1.67095,1.6732,1.67546,1.67773,1.68,1.68228,1.68457,1.68686,1.68916,1.69146,1.69377,1.69609,1.69841,1.70074,1.70308,1.70542,1.70777,1.71013,1.71249,1.71486,1.71723,1.71961,1.722,1.7244,1.7268,1.72921,1.73162,1.73404,1.73647,1.73891,1.74135,1.7438,1.74626,1.74872,1.75119,1.75367,1.75615,1.75864,1.76114,1.76364,1.76616,1.76868,1.7712,1.77374,1.77628,1.77883,1.78138,1.78395,1.78652,1.78909,1.79168,1.79427,1.79687,1.79948,1.8021,1.80472,1.80735,1.80999,1.81263,1.81529,1.81795,1.82062,1.8233,1.82598,1.82868,1.83138,1.83409,1.8368,1.83953,1.84226,1.845,1.84775,1.85051,1.85328,1.85605,1.85883,1.86162,1.86442,1.86723,1.87005,1.87287,1.87571,1.87855,1.8814,1.88426,1.88713,1.89,1.89289,1.89578,1.89869,1.9016,1.90452,1.90745,1.91039,1.91334,1.91629,1.91926,1.92224,1.92522,1.92821,1.93122,1.93423,1.93725,1.94028,1.94333,1.94638,1.94944,1.95251,1.95559,1.95868,1.96178,1.96488,1.968,1.97113,1.97427,1.97742,1.98058,1.98375,1.98693,1.99012,1.99331,1.99652,1.99974,2.00298,2.00622,2.00947,2.01273,2.016,2.01929,2.02258,2.02589,2.0292,2.03253,2.03587,2.03921,2.04257,2.04594,2.04933,2.05272,2.05612,2.05954,2.06296,2.0664,2.06985,2.07331,2.07679,2.08027,2.08377,2.08728,2.0908,2.09433,2.09787,2.10143,2.10499,2.10857,2.11217,2.11577,2.11939,2.12302,2.12666,2.13031,2.13398,2.13766,2.14135,2.14506,2.14877,2.1525,2.15625,2.16,2.16377,2.16756,2.17135,2.17516,2.17898,2.18282,2.18667,2.19053,2.19441,2.1983,2.20221,2.20612,2.21006,2.214,2.21796,2.22194,2.22593,2.22993,2.23395,2.23798,2.24203,2.24609,2.25017,2.25426,2.25836,2.26249,2.26662,2.27077,2.27494,2.27912,2.28332,2.28753,2.29176,2.296,2.30026,2.30454,2.30883,2.31314,2.31746,2.3218,2.32616,2.33053,2.33492,2.33932,2.34375,2.34819,2.35264,2.35711,2.3616,2.36611,2.37063,2.37518,2.37973,2.38431,2.38891,2.39352,2.39815,2.40279,2.40746,2.41214,2.41685,2.42157,2.42631,2.43106,2.43584,2.44063,2.44545,2.45028,2.45513,2.46,2.46489,2.4698,2.47473,2.47968,2.48465,2.48964,2.49465,2.49968,2.50473,2.5098,2.51489,2.52,2.52514,2.53029,2.53546,2.54066,2.54588,2.55111,2.55637,2.56166,2.56696,2.57229,2.57763,2.583,2.5884,2.59381,2.59925,2.60471,2.61019,2.6157,2.62123,2.62678,2.63236,2.63796,2.64359,2.64923,2.65491,2.6606,2.66633,2.67207,2.67784,2.68364,2.68946,2.69531,2.70118,2.70708,2.713,2.71895,2.72493,2.73093,2.73696,2.74301,2.7491,2.7552,2.76134,2.7675,2.7737,2.77991,2.78616,2.79244,2.79874,2.80507,2.81143,2.81782,2.82424,2.83069,2.83717,2.84367,2.85021,2.85678,2.86338,2.87,2.87666,2.88335,2.89007,2.89683,2.90361,2.91043,2.91727,2.92416,2.93107,2.93801,2.94499,2.952,2.95905,2.96613,2.97324,2.98039,2.98757,2.99479,3.00204,3.00932,3.01665,3.024,3.0314,3.03883,3.04629,3.0538,3.06134,3.06892,3.07653,3.08418,3.09187,3.0996,3.10737,3.11518,3.12303,3.13091,3.13884,3.14681,3.15481,3.16286,3.17095,3.17908,3.18725,3.19547,3.20373,3.21203,3.22037,3.22875,3.23719,3.24566,3.25418,3.26274,3.27135,3.28,3.28871,3.29745,3.30624,3.31509,3.32397,3.33291,3.34189,3.35092,3.36001,3.36914,3.37832,3.38755,3.39683,3.40616,3.41554,3.42498,3.43446,3.44401,3.4536,3.46325,3.47295,3.4827,3.49251,3.50238,3.5123,3.52228,3.53231,3.54241,3.55256,3.56276,3.57303,3.58336,3.59374,3.60419,3.6147,3.62527,3.6359,3.64659,3.65735,3.66817,3.67906,3.69001,3.70102,3.7121,3.72325,3.73446,3.74575,3.7571,3.76852,3.78001,3.79157,3.8032,3.8149,3.82667,3.83852,3.85044,3.86244,3.87451,3.88665,3.89887,3.91117,3.92355,3.93601,3.94854,3.96116,3.97385,3.98663,3.99949,4.01243,4.02546,4.03857,4.05177,4.06506,4.07843,4.09189,4.10544,4.11908,4.13281,4.14663,4.16054,4.17455,4.18865,4.20285,4.21715,4.23154,4.24603,4.26062,4.27532,4.29011,4.30501,4.32001,4.33511,4.35032,4.36564,4.38107,4.3966,4.41225,4.42801,4.44388,4.45986,4.47596,4.49218,4.50852,4.52497,4.54155,4.55824,4.57506,4.59201,4.60908,4.62628,4.6436,4.66106,4.67865,4.69637,4.71423,4.73222,4.75035,4.76862,4.78703,4.80559,4.82429,4.84313,4.86212,4.88127,4.90056,4.92001,4.93961,4.95937,10};

    double Nlak33aRefractiveIndex2[76]={1.73816,1.73836,1.73858,1.73881,1.73904,1.73928,1.73952,1.73976,1.74001,1.74026,1.74052,1.74078,1.74105,1.74132,1.7416,1.74189,1.74218,1.74249,1.74279,1.74311,1.74344,1.74378,1.74412,1.74448,1.74485,1.74522,1.74562,1.74602,1.74644,1.74687,1.74732,1.74779,1.74827,1.74878,1.7493,1.74985,1.75042,1.75101,1.75163,1.75228,1.75296,1.75368,1.75443,1.75521,1.75604,1.75692,1.75784,1.75882,1.75985,1.76095,1.76211,1.76335,1.76467,1.76608,1.76758,1.7692,1.77093,1.77279,1.7748,1.77698,1.77934,1.7819,1.7847,1.78775,1.79111,1.79481,1.79889,1.80343,1.8085,1.81419,1.82061,1.8279,1.83625,1.84589,1.85713,1.87039};

    double Nlak33aRefractiveIndex[753]={1.73816,1.73816,1.73819,1.73821,1.73823,1.73825,1.73827,1.7383,1.73832,1.73834,1.73836,1.73838,1.73841,1.73843,1.73845,1.73847,1.7385,1.73852,1.73854,1.73856,1.73858,1.73861,1.73863,1.73865,1.73868,1.7387,1.73872,1.73874,1.73877,1.73879,1.73881,1.73883,1.73886,1.73888,1.7389,1.73893,1.73895,1.73897,1.739,1.73902,1.73904,1.73907,1.73909,1.73911,1.73914,1.73916,1.73918,1.73921,1.73923,1.73925,1.73928,1.7393,1.73932,1.73935,1.73937,1.7394,1.73942,1.73944,1.73947,1.73949,1.73952,1.73954,1.73956,1.73959,1.73961,1.73964,1.73966,1.73968,1.73971,1.73973,1.73976,1.73978,1.73981,1.73983,1.73986,1.73988,1.73991,1.73993,1.73996,1.73998,1.74001,1.74003,1.74006,1.74008,1.74011,1.74013,1.74016,1.74018,1.74021,1.74023,1.74026,1.74028,1.74031,1.74033,1.74036,1.74039,1.74041,1.74044,1.74046,1.74049,1.74052,1.74054,1.74057,1.74059,1.74062,1.74065,1.74067,1.7407,1.74072,1.74075,1.74078,1.7408,1.74083,1.74086,1.74088,1.74091,1.74094,1.74097,1.74099,1.74102,1.74105,1.74107,1.7411,1.74113,1.74116,1.74118,1.74121,1.74124,1.74127,1.74129,1.74132,1.74135,1.74138,1.7414,1.74143,1.74146,1.74149,1.74152,1.74154,1.74157,1.7416,1.74163,1.74166,1.74169,1.74172,1.74174,1.74177,1.7418,1.74183,1.74186,1.74189,1.74192,1.74195,1.74198,1.74201,1.74204,1.74206,1.74209,1.74212,1.74215,1.74218,1.74221,1.74224,1.74227,1.7423,1.74233,1.74236,1.74239,1.74242,1.74245,1.74249,1.74252,1.74255,1.74258,1.74261,1.74264,1.74267,1.7427,1.74273,1.74276,1.74279,1.74283,1.74286,1.74289,1.74292,1.74295,1.74298,1.74302,1.74305,1.74308,1.74311,1.74315,1.74318,1.74321,1.74324,1.74328,1.74331,1.74334,1.74337,1.74341,1.74344,1.74347,1.74351,1.74354,1.74357,1.74361,1.74364,1.74367,1.74371,1.74374,1.74378,1.74381,1.74384,1.74388,1.74391,1.74395,1.74398,1.74402,1.74405,1.74409,1.74412,1.74416,1.74419,1.74423,1.74426,1.7443,1.74433,1.74437,1.74441,1.74444,1.74448,1.74451,1.74455,1.74459,1.74462,1.74466,1.7447,1.74473,1.74477,1.74481,1.74485,1.74488,1.74492,1.74496,1.745,1.74503,1.74507,1.74511,1.74515,1.74519,1.74522,1.74526,1.7453,1.74534,1.74538,1.74542,1.74546,1.7455,1.74554,1.74558,1.74562,1.74566,1.7457,1.74574,1.74578,1.74582,1.74586,1.7459,1.74594,1.74598,1.74602,1.74606,1.7461,1.74614,1.74619,1.74623,1.74627,1.74631,1.74635,1.7464,1.74644,1.74648,1.74653,1.74657,1.74661,1.74665,1.7467,1.74674,1.74679,1.74683,1.74687,1.74692,1.74696,1.74701,1.74705,1.7471,1.74714,1.74719,1.74723,1.74728,1.74732,1.74737,1.74741,1.74746,1.74751,1.74755,1.7476,1.74765,1.74769,1.74774,1.74779,1.74784,1.74788,1.74793,1.74798,1.74803,1.74808,1.74813,1.74818,1.74822,1.74827,1.74832,1.74837,1.74842,1.74847,1.74852,1.74857,1.74862,1.74867,1.74873,1.74878,1.74883,1.74888,1.74893,1.74898,1.74904,1.74909,1.74914,1.74919,1.74925,1.7493,1.74935,1.74941,1.74946,1.74952,1.74957,1.74963,1.74968,1.74974,1.74979,1.74985,1.7499,1.74996,1.75002,1.75007,1.75013,1.75019,1.75024,1.7503,1.75036,1.75042,1.75047,1.75053,1.75059,1.75065,1.75071,1.75077,1.75083,1.75089,1.75095,1.75101,1.75107,1.75113,1.75119,1.75126,1.75132,1.75138,1.75144,1.75151,1.75157,1.75163,1.7517,1.75176,1.75182,1.75189,1.75195,1.75202,1.75208,1.75215,1.75222,1.75228,1.75235,1.75242,1.75248,1.75255,1.75262,1.75269,1.75275,1.75282,1.75289,1.75296,1.75303,1.7531,1.75317,1.75324,1.75332,1.75339,1.75346,1.75353,1.7536,1.75368,1.75375,1.75382,1.7539,1.75397,1.75405,1.75412,1.7542,1.75427,1.75435,1.75443,1.7545,1.75458,1.75466,1.75474,1.75482,1.75489,1.75497,1.75505,1.75513,1.75521,1.7553,1.75538,1.75546,1.75554,1.75562,1.75571,1.75579,1.75587,1.75596,1.75604,1.75613,1.75622,1.7563,1.75639,1.75648,1.75656,1.75665,1.75674,1.75683,1.75692,1.75701,1.7571,1.75719,1.75728,1.75737,1.75747,1.75756,1.75765,1.75775,1.75784,1.75794,1.75803,1.75813,1.75823,1.75832,1.75842,1.75852,1.75862,1.75872,1.75882,1.75892,1.75902,1.75912,1.75922,1.75933,1.75943,1.75954,1.75964,1.75975,1.75985,1.75996,1.76007,1.76017,1.76028,1.76039,1.7605,1.76061,1.76072,1.76083,1.76095,1.76106,1.76117,1.76129,1.7614,1.76152,1.76164,1.76175,1.76187,1.76199,1.76211,1.76223,1.76235,1.76247,1.7626,1.76272,1.76284,1.76297,1.76309,1.76322,1.76335,1.76348,1.76361,1.76374,1.76387,1.764,1.76413,1.76426,1.7644,1.76453,1.76467,1.7648,1.76494,1.76508,1.76522,1.76536,1.7655,1.76564,1.76579,1.76593,1.76608,1.76622,1.76637,1.76652,1.76667,1.76682,1.76697,1.76712,1.76727,1.76743,1.76758,1.76774,1.7679,1.76805,1.76821,1.76837,1.76854,1.7687,1.76886,1.76903,1.7692,1.76936,1.76953,1.7697,1.76987,1.77005,1.77022,1.77039,1.77057,1.77075,1.77093,1.77111,1.77129,1.77147,1.77166,1.77184,1.77203,1.77222,1.77241,1.7726,1.77279,1.77299,1.77318,1.77338,1.77358,1.77378,1.77398,1.77418,1.77439,1.7746,1.7748,1.77501,1.77523,1.77544,1.77565,1.77587,1.77609,1.77631,1.77653,1.77675,1.77698,1.77721,1.77744,1.77767,1.7779,1.77813,1.77837,1.77861,1.77885,1.77909,1.77934,1.77958,1.77983,1.78008,1.78034,1.78059,1.78085,1.78111,1.78137,1.78164,1.7819,1.78217,1.78244,1.78271,1.78299,1.78327,1.78355,1.78383,1.78412,1.78441,1.7847,1.78499,1.78529,1.78558,1.78589,1.78619,1.7865,1.78681,1.78712,1.78744,1.78775,1.78808,1.7884,1.78873,1.78906,1.78939,1.78973,1.79007,1.79041,1.79076,1.79111,1.79146,1.79182,1.79218,1.79254,1.79291,1.79328,1.79366,1.79404,1.79442,1.79481,1.7952,1.79559,1.79599,1.79639,1.7968,1.79721,1.79762,1.79804,1.79847,1.79889,1.79933,1.79976,1.80021,1.80065,1.8011,1.80156,1.80202,1.80249,1.80296,1.80343,1.80392,1.8044,1.8049,1.80539,1.8059,1.80641,1.80692,1.80744,1.80797,1.8085,1.80904,1.80959,1.81014,1.8107,1.81126,1.81183,1.81241,1.813,1.81359,1.81419,1.8148,1.81541,1.81603,1.81666,1.8173,1.81795,1.8186,1.81926,1.81993,1.82061,1.8213,1.82199,1.8227,1.82341,1.82414,1.82487,1.82561,1.82637,1.82713,1.8279,1.82869,1.82948,1.83029,1.8311,1.83193,1.83277,1.83362,1.83449,1.83536,1.83625,1.83715,1.83807,1.839,1.83994,1.84089,1.84186,1.84285,1.84385,1.84486,1.84589,1.84694,1.848,1.84908,1.85017,1.85129,1.85242,1.85357,1.85474,1.85593,1.85713,1.85836,1.85961,1.86088,1.86217,1.86348,1.86481,1.86617,1.86756,1.86896,1.87039,1.87185,3.5};

  }

  TCanvas *c1 = new TCanvas("c1","hists with different scales",600,400);

  TGraph *gQuartz = new TGraph(num,PhotonEnergy,QuartzRefractiveIndex);
  TGraph *gEpotek = new TGraph(num,PhotonEnergy,EpotekRefractiveIndex);
  TGraph *gKoil = new TGraph(num,PhotonEnergy,KamLandOilRefractiveIndex);
  TGraph *gNlak33a = new TGraph(76,enNlak33a2,Nlak33aRefractiveIndex2);
  
  gQuartz->SetMarkerStyle(7);
  gEpotek->SetMarkerStyle(7);
  gKoil->SetMarkerStyle(7);
  gNlak33a->SetMarkerStyle(7);
  gQuartz->GetYaxis()->SetRangeUser(1,2);
  gQuartz->GetYaxis()->SetTitle("refractive index");
  gQuartz->GetXaxis()->SetTitle("photon energy, [eV]");

  gQuartz->SetLineColor(2);
  gEpotek->SetLineColor(3);
  gKoil->SetLineColor(1);
  gNlak33a->SetLineColor(4);

  gQuartz->SetMarkerColor(2);
  gEpotek->SetMarkerColor(3);
  gKoil->SetMarkerColor(1);
  gNlak33a->SetMarkerColor(4);

  gQuartz->Draw("ACP");
  gEpotek->Draw("same CP");
  gKoil->Draw("same CP");
  gNlak33a->Draw("same CP");

  gStyle->SetOptTitle(0);
  gQuartz->GetXaxis()->SetRangeUser(2,4);
  c1->Update();  
  TF1 *f1=new TF1("f1","1.24/x",620,310); //2-4 eV

  TGaxis *A1 = new TGaxis(gPad->GetUxmin()+0.02,gPad->GetUymax(), gPad->GetUxmax(), gPad->GetUymax(),"f1",510,"-");
  A1->SetTitle("photon wavelength, [nm]");
  A1->SetTitleFont(42);
  A1->SetTitleSize(0.035);
  A1->SetLabelFont(42);
  A1->SetLabelSize(0.035);
  A1->Draw();

  Double_t xl1=.5, yl1=0.75, xl2=xl1+.3, yl2=yl1+.125;
  TLegend *leg = new TLegend(xl1,yl1,xl2,yl2);
  leg->SetFillColor(0);
  leg->SetLineColor(0);
  leg->AddEntry(gQuartz,"quartz","l");
  leg->AddEntry(gEpotek,"epotek","l");
  leg->AddEntry(gKoil,"kamland oil","l");
  leg->AddEntry(gNlak33a,"n-lak33a","l");
  leg->Draw("same");

  TCanvas *c2 = new TCanvas("c2","hists with different scales",600,400);

  TGraph *g2Quartz = new TGraph(num,PhotonEnergy,QuartzAbsorption);
  TGraph *g2Epotek = new TGraph(num,PhotonEnergy,EpotekAbsorption);
  TGraph *g2Koil = new TGraph(num,PhotonEnergy,KamLandOilAbsorption);
  TGraph *g2Nlak33a = new TGraph(76,enNlak33a2,Nlak33aAbsorption2);
  
  g2Quartz->SetMarkerStyle(7);
  g2Epotek->SetMarkerStyle(7);
  g2Koil->SetMarkerStyle(7);
  g2Nlak33a->SetMarkerStyle(7);
  // g2Quartz->GetYaxis()->SetRangeUser(1,2);
  g2Quartz->GetYaxis()->SetTitleOffset(1.42);
  g2Quartz->GetYaxis()->SetTitle("mean free path length, [cm]");
  g2Quartz->GetXaxis()->SetTitle("photon energy, [eV]");

  g2Quartz->SetLineColor(2);
  g2Epotek->SetLineColor(3);
  g2Koil->SetLineColor(1);
  g2Nlak33a->SetLineColor(4);

  g2Quartz->SetMarkerColor(2);
  g2Epotek->SetMarkerColor(3);
  g2Koil->SetMarkerColor(1);
  g2Nlak33a->SetMarkerColor(4);

  g2Quartz->Draw("ALP");
  g2Epotek->Draw("same LP");
  g2Koil->Draw("same LP");
  g2Nlak33a->Draw("same LP");

  gStyle->SetOptTitle(0);
  g2Quartz->GetXaxis()->SetRangeUser(2,4);
  c2->Update();  
  TF1 *f2=new TF1("f2","1.24/x",620,310); //2-4 eV

  TGaxis *A2 = new TGaxis(gPad->GetUxmin()+0.02,gPad->GetUymax(), gPad->GetUxmax(), gPad->GetUymax(),"f2",510,"-");
  A2->SetTitle("photon wavelength, [nm]");
  A2->SetTitleFont(42);
  A2->SetTitleSize(0.035);
  A2->SetLabelFont(42);
  A2->SetLabelSize(0.035);
  A2->Draw();

  TLegend *leg2 = new TLegend(xl1,yl1,xl2,yl2);
  leg2->SetFillColor(0);
  leg2->SetLineColor(0);
  leg2->AddEntry(g2Quartz,"quartz","l");
  leg2->AddEntry(g2Epotek,"epotek","l");
  leg2->AddEntry(g2Koil,"kamland oil","l");
  leg2->AddEntry(g2Nlak33a,"n-lak33a","l");
  leg2->Draw("same");
}