FdCherenkovSubtracter.cc

Go to the documentation of this file.

#include "FdCherenkovSubtracter.h"

#include <utl/ErrorLogger.h>
#include <utl/MathConstants.h>
#include <utl/Vector.h>
#include <utl/PhysicalConstants.h>
#include <utl/TabulatedFunctionErrors.h>

#include <evt/Event.h>
#include <evt/ShowerRecData.h>
#include <evt/ShowerFRecData.h>

#include <fevt/FEvent.h>
#include <fevt/Eye.h>
#include <fevt/Telescope.h>
#include <fevt/EyeRecData.h>

#include <det/Detector.h>
#include <atm/Atmosphere.h>
#include <atm/AttenuationResult.h>

#include <fdet/Eye.h>
#include <fdet/Pixel.h>
#include <fdet/Channel.h>

#include <utl/Vector.h>

using namespace FdCherenkovSubtracterOG;
using namespace fwk;
using namespace utl;
using namespace evt;
using namespace fevt;
using namespace det;
using namespace fdet;
using namespace std;
using namespace atm;


FdCherenkovSubtracter::FdCherenkovSubtracter(){}

fwk::VModule::ResultFlag FdCherenkovSubtracter::Init(){
  return eSuccess;
}

fwk::VModule::ResultFlag FdCherenkovSubtracter::Finish(){
  return eSuccess;
}

fwk::VModule::ResultFlag FdCherenkovSubtracter::Run(evt::Event& event){
  
  if (!event.HasFEvent()) return eSuccess;

  FEvent::EyeIterator eyeIter;
  
  for (eyeIter= event.GetFEvent().EyesBegin();
       eyeIter != event.GetFEvent().EyesEnd();
       ++eyeIter){                    
    
    fevt::Eye& eye = *eyeIter;
    
    if(!eye.GetRecData().HasLightFlux(fevt::FdConstants::eCherDirect)) continue;
    SubtractCherenkov(eye); 

    
  } // end loop over eyes

  return eSuccess;
}

void FdCherenkovSubtracter::SubtractCherenkov(fevt::Eye& eye){
  
  TabulatedFunctionErrors& totalFlux = 
    eye.GetRecData().GetLightFlux(fevt::FdConstants::eTotal);

   TabulatedFunctionErrors& fluoFlux = 
     eye.GetRecData().GetLightFlux(fevt::FdConstants::eFluorTotal);
  
  TabulatedFunctionErrors& directCherenkovFlux = 
    eye.GetRecData().GetLightFlux(fevt::FdConstants::eCherDirect);
  
  TabulatedFunctionErrors& rayscattCherenkovFlux = 
    eye.GetRecData().GetLightFlux(fevt::FdConstants::eCherRayleighScattered);
  
  TabulatedFunctionErrors& miescattCherenkovFlux = 
    eye.GetRecData().GetLightFlux(fevt::FdConstants::eCherMieScattered); 
  
  
  for (unsigned int i=0; i<totalFlux.GetNPoints(); ++i) {
    
    fluoFlux.GetY(i)  = totalFlux.GetY(i)
                      - directCherenkovFlux.GetY(i)
                      - rayscattCherenkovFlux.GetY(i)
                      - miescattCherenkovFlux.GetY(i);
    
  }     
  
} // end of SubtractCherenkov


// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: