mc_tel.h

Go to the documentation of this file.

/* ============================================================================

   Copyright (C) 1997, 2001, 2009, 2010, 2019  Konrad Bernloehr

   This file is part of the eventio/hessio library.

   The eventio/hessio library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.

   This library is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with this library. If not, see <http://www.gnu.org/licenses/>.

============================================================================ */

/* ========================================================= */

#ifndef _MC_TEL_LOADED

#define _MC_TEL_LOADED 2

#ifndef _EVENTIO_BASIC_LOADED
# include "io_basic.h"
#endif

#include "mc_atmprof.h"

#ifdef __cplusplus
extern "C" {
#endif

/* Data types: */

#ifdef DOUBLEPREC
typedef double real;
#else
typedef float real;
#endif

typedef short INT16;
typedef unsigned short UINT16;
#ifdef OS_MSDOS
typedef long INT32;
typedef unsigned long UINT32;
#else
typedef int INT32;
typedef unsigned int UINT32;
#endif

/* Data structures: */
struct bunch
{
   float photons; 
   float x, y;    
   float cx, cy;  
   float ctime;   
   float zem;     
   float lambda;  
};

struct compact_bunch
{
   short photons; 
   short x, y;    
   short cx, cy;  
   short ctime;   
   short log_zem; 
   short lambda;  
};

struct photo_electron
{
   int pixel;     
   int lambda;    
   double atime;  
   /* double amplitude; */ 
};

struct linked_string
{
   char *text;
   struct linked_string *next;
};

struct shower_extra_parameters
{
   long id;       
   int is_set;    
   double weight; 
   size_t niparam;
   int *iparam;   
   size_t nfparam;
   float *fparam;   
};
typedef struct shower_extra_parameters ShowerExtraParam;

/* I/O item types: */

/* Never change the following numbers after MC data is created: */
#define IO_TYPE_MC_BASE     1200
#define IO_TYPE_MC_RUNH     (IO_TYPE_MC_BASE+0)
#define IO_TYPE_MC_TELPOS   (IO_TYPE_MC_BASE+1)
#define IO_TYPE_MC_EVTH     (IO_TYPE_MC_BASE+2)
#define IO_TYPE_MC_TELOFF   (IO_TYPE_MC_BASE+3)
#define IO_TYPE_MC_TELARRAY (IO_TYPE_MC_BASE+4)
#define IO_TYPE_MC_PHOTONS  (IO_TYPE_MC_BASE+5)
#define IO_TYPE_MC_LAYOUT   (IO_TYPE_MC_BASE+6)
#define IO_TYPE_MC_TRIGTIME (IO_TYPE_MC_BASE+7)
#define IO_TYPE_MC_PE       (IO_TYPE_MC_BASE+8)
#define IO_TYPE_MC_EVTE     (IO_TYPE_MC_BASE+9)
#define IO_TYPE_MC_RUNE     (IO_TYPE_MC_BASE+10)
#define IO_TYPE_MC_LONGI    (IO_TYPE_MC_BASE+11)
#define IO_TYPE_MC_INPUTCFG (IO_TYPE_MC_BASE+12)
#define IO_TYPE_MC_TELARRAY_HEAD (IO_TYPE_MC_BASE+13)
#define IO_TYPE_MC_TELARRAY_END (IO_TYPE_MC_BASE+14)
#define IO_TYPE_MC_EXTRA_PARAM (IO_TYPE_MC_BASE+15)
#define IO_TYPE_MC_ATMPROF  (IO_TYPE_MC_BASE+16)

/* Function prototypes: */

/* io_telescope.c */
int write_tel_block (IO_BUFFER *iobuf, int type, int num,
      real *data, int len);
int read_tel_block (IO_BUFFER *iobuf, int type, real *data,
      int maxlen);
int print_tel_block (IO_BUFFER *iobuf);

int write_input_lines (IO_BUFFER *iobuf, 
      struct linked_string *list);
int read_input_lines (IO_BUFFER *iobuf, 
      struct linked_string *list);

int write_tel_pos (IO_BUFFER *iobuf, int ntel, double *x,
      double *y, double *z, double *r);
int read_tel_pos (IO_BUFFER *iobuf, int max_tel, int *ntel,
      double *x, double *y, double *z, double *r);
int print_tel_pos (IO_BUFFER *iobuf);

int write_tel_offset (IO_BUFFER *iobuf, int narray, double toff,
      double *xoff, double *yoff);
int write_tel_offset_w (IO_BUFFER *iobuf, int narray, double toff,
      double *xoff, double *yoff, double *weight);
int read_tel_offset (IO_BUFFER *iobuf, int max_array, int *narray,
      double *toff, double *xoff, double *yoff);
int read_tel_offset_w (IO_BUFFER *iobuf, int max_array, int *narray,
      double *toff, double *xoff, double *yoff, double *weight);
int print_tel_offset (IO_BUFFER *iobuf);

int begin_write_tel_array (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih, int array);
int end_write_tel_array (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih);
int begin_read_tel_array (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih, int *array);
int end_read_tel_array (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih);

int write_tel_array_head (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih, int array);
int write_tel_array_end (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih, int array);
int read_tel_array_head (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih, int *array);
int read_tel_array_end (IO_BUFFER *iobuf, IO_ITEM_HEADER *ih, int *array);

int write_tel_photons (IO_BUFFER *iobuf, int array, int tel,
      double photons, struct bunch *bunches, int nbunches,
      int ext_bunches, char *ext_fname);
int write_tel_compact_photons (IO_BUFFER *iobuf, int array, int tel,
      double photons, struct compact_bunch *cbunches, int nbunches,
      int ext_bunches, char *ext_fname);
int read_tel_photons (IO_BUFFER *iobuf, int max_bunches, int *array,
      int *tel, double *photons, struct bunch *bunches, int *nbunches);
int print_tel_photons (IO_BUFFER *iobuf);

int write_shower_longitudinal (IO_BUFFER *iobuf, int event, int type, 
      double *data, int ndim, int np, int nthick, double thickstep);
int read_shower_longitudinal (IO_BUFFER *iobuf, int *event, 
      int *type, double *data, int ndim, int *np, int *nthick, 
      double *thickstep, int max_np);
int print_shower_longitudinal (IO_BUFFER *iobuf);

int write_camera_layout (IO_BUFFER *iobuf, int itel, int type, 
      int pixels, double *xp, double *yp);
int read_camera_layout (IO_BUFFER *iobuf, int max_pixels, int *itel,
      int *type, int *pixels, double *xp, double *yp);
int print_camera_layout (IO_BUFFER *iobuf);

int write_photo_electrons (IO_BUFFER *iobuf, int array, int tel, 
      int npe, int pixels, int flags, int *pe_counts, int *tstart, 
      double *t, double *a, int *photon_counts);
int read_photo_electrons (IO_BUFFER *iobuf, int max_pixel,
      int max_pe, int *array, int *tel, int *npe, int *pixels, int *flags, 
      int *pe_counts, int *tstart, double *t, double *a, int *photon_counts);
int print_photo_electrons (IO_BUFFER *iobuf);

int write_shower_extra_parameters (IO_BUFFER *iobuf, ShowerExtraParam *ep);
int read_shower_extra_parameters (IO_BUFFER *iobuf, ShowerExtraParam *ep);
int print_shower_extra_parameters (IO_BUFFER *iobuf);
int init_shower_extra_parameters(ShowerExtraParam *ep, size_t ni_max, size_t nf_max);
int clear_shower_extra_parameters (ShowerExtraParam *ep);
ShowerExtraParam *get_shower_extra_parameters(void);

void remember_corsika_atm_params(float *hlay, float *aatm, float *batm, float *catm);
int get_corsika_atm_params(double *hlay, double *aatm, double *batm, double *catm, double *datm);

int write_atmprof (IO_BUFFER *iobuf, AtmProf *atmprof);
int read_atmprof (IO_BUFFER *iobuf, AtmProf *atmprof);
int print_atmprof (IO_BUFFER *iobuf);

#ifdef __cplusplus
}
#endif

#endif