Telescope.h

Go to the documentation of this file.

//  Author : florent robinet (IJCLab - Orsay): florent.robinet@ijclab.in2p3.fr
#ifndef __Telescope__
#define __Telescope__
 
#include <TGraph.h>
#include <TH3.h>
#include "Camera.h"
#include "Source.h"
 
using namespace std;
 
class Telescope: public Camera {
  
 public:
  
 
  Telescope(const int aNbits=14, const int aNy=256, const int aNz=256);
  
 
  virtual ~Telescope(void);  
 
  void TakePicture(const ULong64_t aTimeStart);
 
 
  void StartObservation(TFile *aRootFile);
 
 
  void StopObservation(void);
 
 
  void SetName(const string aName);
  
 
  inline string GetName(void){ return tel_name; };
  
 
  void AddSource(Source *aSource, TGraph *aEffectiveArea, TH2F *aAcceptance);
 
 
  TGraph* GetEffectiveArea(const UInt_t aSourceIndex){
    return tel_effarea[aSourceIndex];
  };
  
 
  double GetEffectiveArea(const UInt_t aSourceIndex, const double aEnergy);
  
 
  double GetEffectiveArea(const UInt_t aSourceIndex, const double aEnergyMin, const double aEnergyMax);
  
 
  inline void SetFocalLength(const int aFocalLength=1000){ tel_foclen=aFocalLength; };
 
 
  inline int GetFocalLength(void){ return tel_foclen; };
 
 
  bool SetPsfFile(const string aFileName);
 
 
  bool SetPsf(TH3F* aPsf);
 
 
  bool SetPsf(TH2F* aPsf, const double aEnergyMin=0.0, const double aEnergyMax=1000000);
 
 
  double SetPsfAiry(const double aOpeningAngle,
                    const int aYN, const int aZN,
                    const int aEnergyN, const double aEnergyMin, const double aEnergyMax);
 
 
  inline TH3F* GetPsf(void){ return tel_psf; };
 
 
  TH2F* GetPsf(const double aEnergy);
 
 
  TH2F* GetPsf(const UInt_t aEnergyBinIndex);
 
 private:
 
  string tel_name;                   
  int tel_foclen;                    
  int tel_psf_nebins;                
  TArrayD *tel_psf_ebins;            
  TH3F *tel_psf;                     
  TH2F **tel_psfe;                   
  vector <TGraph*> tel_effarea;      
  vector <TH2F*> tel_acc;            
  vector <Source*> tel_src;          
  TFile *tel_file;                   
 
 
  void ApplyPsf(double& aYf, double& aZf, const double aY0, const double aZ0, const double aEnergy);
 
  ClassDef(Telescope,0)
};
 
#endif