Line data    Source code

       1             : /**
       2             :    \file path_measure.cpp
       3             :    
       4             :    Renamed path_measure.cc 2023
       5             : 
       6             :    Bojan Nikolic <b.nikolic@mrao.cam.ac.uk>, <bojan@bnikolic.co.uk>
       7             :    
       8             : */   
       9             : 
      10             : #include "path_measure.h"
      11             : #include "model_iface.h"
      12             : #include "numalgo.h"
      13             : 
      14             : namespace LibAIR2 {
      15             : 
      16           0 :   PathMeasure::PathMeasure(WVRAtmoQuantModel &model) :
      17             :     ALMAMeasure(model),
      18           0 :     Tb_obs(4),
      19           0 :     Tb_sigma(4),
      20           0 :     dT_dL(4),
      21           0 :     dT_dL_sigma(4)
      22             :   {
      23           0 :   }
      24             : 
      25           0 :   void PathMeasure::modelObs(void)
      26             :   {
      27           0 :     model.eval(Tb_obs);
      28           0 :     model.dTdL_ND(dT_dL);    
      29           0 :   }
      30             : 
      31           0 :   double PathMeasure::lLikely (void) const
      32             :   {
      33           0 :     double ltot=0;
      34             : 
      35           0 :     std::vector<double> res;
      36           0 :     model.eval(res);
      37             : 
      38           0 :     ltot+=GaussError( Tb_obs.begin(), res, Tb_sigma.thermNoise);
      39             :     
      40           0 :     for (size_t i =0 ; i < 4; ++i)
      41           0 :       res[i]=model.dTdL_ND(i);
      42             : 
      43           0 :     ltot+=GaussError( dT_dL.begin(), res, dT_dL_sigma.thermNoise);    
      44             : 
      45           0 :     return ltot;
      46           0 :   }
      47             : 
      48             : }
      49             :