Index: trunk/MagicSoft/Simulation/Corsika/Checkmc/Makefile =================================================================== --- trunk/MagicSoft/Simulation/Corsika/Checkmc/Makefile (revision 328) +++ trunk/MagicSoft/Simulation/Corsika/Checkmc/Makefile (revision 328) @@ -0,0 +1,220 @@ +################################################################## +# +# makefile +# +# @file makefile +# @title Simulation of the camera and trigger logic +# @author J C Gonz\'alez +# @email gonzalez@mppmu.mpg.de +# @date Fri Mar 12 11:51:11 MET 1999 +# +#_______________________________________________________________ +# +# Created: Fri Mar 12 11:51:11 MET 1999 +# Author: Jose Carlos Gonzalez +# Purpose: Makefile for the compilation of the camera program +# Notes: +# +#--------------------------------------------------------------- +# +# $RCSfile: Makefile,v $ +# $Revision: 1.1.1.1 $ +# $Author: magiccvs $ +# $Date: 1999-11-25 20:34:58 $ +# +################################################################## +# @maintitle + +# @code + +INCLUDEMK = config.mk.${OSTYPE} +include ${INCLUDEMK} + +# @endcode + +# @code + +# includes + +INCLUDE = ../../Detector/include-GENERAL +INCLUDE_COR = ../../Detector/include-CORSIKA +INCLUDE_CHECK = . + + +# common flags + +INCLUDES = -I${INCLUDE} \ + -I${INCLUDE_COR} \ + -I${INCLUDE_CHECK} \ + -I${INCLUDE_CERN} \ + -I/usr/include + +CERNLIBDIR = ${CERNDIR}/pro/lib/ +CERNLIB = -L${CERNLIBDIR} -lgraflib -lgrafX11 -lpacklib \ + -lkernlib -lpawlib -lmathlib + +RANLIBDIR = ../../Detector/lib +RANLIB = -L${RANLIBDIR} -lranlib + +# special flags + +osf_FORLIBS = -lUfor -lfor -lutil -lots -lm +linux_FORLIBS = -L/usr/local/lib/gcc-lib/i686-pc-linux-gnu/egcs-2.91.57 \ +-L/usr/local/i686-pc-linux-gnu/lib \ +-L/usr/local/lib \ +/usr/X11R6/lib/libXt.a \ +/usr/X11R6/lib/libXext.a \ +/usr/X11R6/lib/libXp.a \ +/usr/X11R6/lib/libX11.a \ +-lcrypt -ldl -lg2c -lm -lgcc -lc -lgcc /usr/lib/crtn.o +generic_FORLIBS = -lm + +FORLIBS = ${${SYSTEM}_FORLIBS} + +# compilation and linking flags + +CXXFLAGS = -D__${SYSTEM}__ ${INCLUDES} ${OPTIM} ${DEBUG} +CFLAGS = ${CXXFLAGS} +FFLAGS = ${CXXFLAGS} +LIBS = ${RANLIB} ${CERNLIB} ${FORLIBS} ${LIBS_linux} +#LIBS = ${CERNLIB} ${RANLIB} ${FORLIBS} + +#------------------------------------------------------------------------------ + +#.SILENT: + +.SUFFIXES: .c .cxx .C .c++ .h .hxx .H .h++ .o .so .f + +SRCS = \ + ${INCLUDE_COR}/COREventHeader.cxx \ + ${INCLUDE_COR}/CORParticle.cxx \ + ${INCLUDE_COR}/CORStatfile.cxx \ + checkmc.cxx + +HEADERS = \ + COREventHeader.hxx \ + CORParticle.hxx \ + CORStatfile.hxx \ + checkmc.h + +OBJS = \ + ${INCLUDE_COR}/COREventHeader.o \ + ${INCLUDE_COR}/CORParticle.o \ + ${INCLUDE_COR}/CORStatfile.o \ + checkmc.o + + +PROGRAM=checkmc_${OSTYPE} + +############################################################ + +all: ${PROGRAM} + +depend: + @makedepend $(SRCS) -fMakefile 2> /dev/null + +doc: reflector-doc + +reflector-doc: + @echo "Generating documentation for camera . . . " + $(DOCUM) -latex -o reflector.tex \ + reflector.cxx reflector.h \ + readparam.cxx readparam.h \ + atm.cxx atm.h + latex "\nonstopmode\input{reflector.tex}" && \ + makeindex reflector && \ + latex "\nonstopmode\input{reflector.tex}" && \ + latex "\nonstopmode\input{reflector.tex}" + @echo "Files reflector.tex and reflector.dvi generated." + +rate: + @echo "Rating documentation inside code . . . " + $(RATE) \ + reflector.cxx reflector.h \ + readparam.cxx readparam.h \ + atm.cxx atm.h + +${PROGRAM}: $(OBJS) + @echo "Linking..." $@ + $(CXX) $(CXXFLAGS) $(OBJS) $(LIBS) -o $@ + @echo "done." + +.cxx.o: + @echo "Compiling " $< + $(CXX) $(CXXFLAGS) -c $< -o $@ + +.c.o: + @echo "Compiling " $< + $(CC) $(CFLAGS) -c $< -o $@ + +.f.o: + @echo "Compiling " $< + $(F77) $(FFLAGS) -c $< -o $@ + +lclean: + @echo "Cleanning..." + @rm -f *.o core + +clean: + @echo "Cleanning..." + @rm -f $(OBJS) core + +mrproper: clean + @echo "Mr.Proper in action . . ." + @rm -rf test + @rm -f $(PROGRAM) + +ctags: + @echo "Creating CTAGS file . . ." + @ctags -txw $(SRCS) $(HEADERS) > CTAGS + +etags: + @echo "Creating TAGS file . . ." + @etags -C $(SRCS) $(HEADERS) + +listsrc: + @ls -m $(SRCS) $(HEADERS) | sed 's/,//g' + +redo: clean all + +# @endcode + +# DO NOT DELETE THIS LINE -- make depend depends on it. + +../../Detector/include-CORSIKA/COREventHeader.o: ../../Detector/include-CORSIKA/COREventHeader.hxx +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/stdlib.h +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/standards.h +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/getopt.h +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/sys/types.h +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/mach/machine/vm_types.h +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/sys/select.h +../../Detector/include-CORSIKA/COREventHeader.o: /usr/include/math.h +../../Detector/include-CORSIKA/CORParticle.o: ../../Detector/include-CORSIKA/CORParticle.hxx +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/stdlib.h +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/standards.h +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/getopt.h +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/sys/types.h +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/mach/machine/vm_types.h +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/sys/select.h +../../Detector/include-CORSIKA/CORParticle.o: /usr/include/math.h +../../Detector/include-CORSIKA/CORStatfile.o: ../../Detector/include-CORSIKA/CORStatfile.hxx +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/stdlib.h +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/standards.h +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/getopt.h +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/sys/types.h +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/mach/machine/vm_types.h +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/sys/select.h +../../Detector/include-CORSIKA/CORStatfile.o: /usr/include/math.h +checkmc.o: checkmc.h /usr/include/stdlib.h /usr/include/standards.h +checkmc.o: /usr/include/getopt.h /usr/include/sys/types.h +checkmc.o: /usr/include/mach/machine/vm_types.h /usr/include/sys/select.h +checkmc.o: /usr/include/stdio.h /usr/include/sys/seek.h +checkmc.o: /usr/include/va_list.h /usr/include/sys/limits.h +checkmc.o: /usr/include/sys/machine/machlimits.h /usr/include/sys/syslimits.h +checkmc.o: /usr/include/sys/machine/machtime.h /usr/include/sys/rt_limits.h +checkmc.o: /usr/include/string.h /usr/include/strings.h /usr/include/stdarg.h +checkmc.o: /usr/include/math.h /usr/include/float.h /usr/include/fp_class.h +checkmc.o: /usr/include/dirent.h +checkmc.o: ../../Detector/include-CORSIKA/COREventHeader.hxx +checkmc.o: ../../Detector/include-CORSIKA/CORParticle.hxx +checkmc.o: ../../Detector/include-CORSIKA/CORStatfile.hxx Index: trunk/MagicSoft/Simulation/Corsika/Checkmc/checkmc =================================================================== --- trunk/MagicSoft/Simulation/Corsika/Checkmc/checkmc (revision 328) +++ trunk/MagicSoft/Simulation/Corsika/Checkmc/checkmc (revision 328) @@ -0,0 +1,32 @@ +#!/bin/csh -f +date + +if ( -e test) then + +else + mkdir ${PWD}/test +endif + +cp checkmc.kumac test/ +cp checkmc_reference.hbook test/ + + +rm -f test/checkmc.log +rm -f test/checkmc.ps +rm -f test/checkmc_paw.log + + +echo "checkmc_${OSTYPE} < input_example >& test/checkmc.log " +checkmc_${OSTYPE} < input_example >& test/checkmc.log + +cd test + +pawX11 >& checkmc_paw.log << EOF +1 +exec checkmc +exit +EOF +echo " Checkmc finished ... " +echo " Output in files checkmc.log" +echo " and checkmc.ps" +date Index: trunk/MagicSoft/Simulation/Corsika/Checkmc/checkmc.cxx =================================================================== --- trunk/MagicSoft/Simulation/Corsika/Checkmc/checkmc.cxx (revision 328) +++ trunk/MagicSoft/Simulation/Corsika/Checkmc/checkmc.cxx (revision 328) @@ -0,0 +1,1466 @@ +////////////////////////////////////////////////////////////////// +// +// checkmc +// +// @file checkmc.cxx +// @title Check of the MC generated data +// @author J C Gonz\'alez, P. Morawitz, D. Petry +// @email gonzalez@mppmu.mpg.de, morawitz@cern.ch, petry@ifae.es +// +// @maintitle +//_______________________________________________________________ +// +// Created: Thu May 7 16:24:22 1998 +// Author: Jose Carlos Gonzales et al. +// Purpose: Check of the MC generated data +// Notes: version 27-10-98 +// +// + +// Below you can find a sort of, incomplete for sure, documentation +// concerning this program. Please, do not hesiate to contact the +// author in case of problems. + +// @T \newpage + +// @section Source code of {\tt checkmc.cxx} + +/* @text +In this section we show the (commented) code of the program +for for the check of the MC generated data, in the current version. +@endtext */ + +// @subsection Includes and Global variables definition + +/* @text +All the defines are located in the file {\tt checkmc.h}. +@endtext */ + +// @code +#include "checkmc.h" +// @endcode + +// @subsubsection Definition of global variables. + +/* @text +Here we define the global variables where the values from +the parameters file are stored. See the section +\ref{paramdesc} in page \pageref{paramdesc} for a more +complete description of the available commands in the +parameters file. +@endtext */ + +// @code +static char **Paths_list; // list of paths to get data from +static int Num_of_paths = 1; // number of paths +static char Output_filename[PATH_MAX_LENGTH]; // output filename +static char CT_filename[PATH_MAX_LENGTH]; // name of the CT def. file +static VerboseLevel verbose = VERBOSE_DEFAULT; // verboseness +static float fixed_Theta; // fixed target theta and phi +static float fixed_Phi; +static float low_Ecut = 0.0, high_Ecut = 100000.0; +static int is_Fixed_Target = FALSE; // are we using fixed target? +static int max_Events = 100000; // maximum number of events to read + +static float gammas = 1.; +static float protons = 14.; +static float helium = 402.; + + +inline int max(int a, int b) +{if (a>b) return a; + else return b;} + +inline int min(int a, int b) +{if (a xmax/dichte == penetration depth +// depth defined as == 120km=0 +// (120km - penetration depth - observation level (2km)) == real height of interaction +// + "xcore", +// x position of core on ground in cm + "ycore", +// y position of core on ground in cm + "Nph", +// number of cerenkov photons + "RunNum", + "EvtNum", + "DateRun", + "Version", + "HeigLev", + "ESlope", + "ELowLim", + "EUppLim", +// + "Ecutofh", + "Ecutofm", + "Ecutofe", + "Ecutofg", +// + "EGS4yn", + "NKGyn", + "GHEISHA", + "VENUS", + "CERENKO", + "NEUTRIN", + "HORIZON", + "COMPUTE", + "ThetMin", + "ThetMax", + "PhiMin", + "PhiMax", +// + "StepLen", + "CWavLow", + "CWavUpp" + }; + // char chTags[nvar * 10]; + float rec[nvar]; + int record_size=4096; + // @endcode + + + // @text + // Variables for the cerenkov photon N-Tuple + // @endtext + + // @code + int hid2=2,istat2=0,icycle2=0; + const int nvar2=11; + char chTags2[nvar2][8]={ +// +// First info about primary + "energy", + "cored", + "theta", +// +// Now cerenkov photon info + "Wavel", +// wavelength in nm + "xph", + "yph", + "uph", + "vph", + "tph", + "hph", + "EvtNum" + }; + // char chTags2[nvar2 * 10]; + float rec2[nvar2]; + int record_size2=4096; + // @endcode + + + // @text + // We start in the main program. First we (could) make some presentation, + // and follows the reading of the parameters file (now from the {\tt stdin}), + // the reading of the CT parameters file, and the creation of the + // output file, where the processed data will be stored. + // @endtext + // @code + + //++ + // START + //-- + + // make unbuffered output + + cout.setf ( ios::stdio ); + + // HBOOK core memory + + HLIMIT(PAWC_SIZE); + + /* + strcpy( chTags, "primary:" ); + strcat( chTags, "energy:" ); + strcat( chTags, "cored:" ); + strcat( chTags, "theta:" ); + strcat( chTags, "phi:" ); + strcat( chTags, "thick:" ); + strcat( chTags, "height:" ); + strcat( chTags, "target1:" ); + strcat( chTags, "px:" ); + strcat( chTags, "py:" ); + strcat( chTags, "pz:" ); + strcat( chTags, "tfirst:" ); + strcat( chTags, "tlast:" ); + strcat( chTags, "xmax" ); + */ + + // make some sort of presentation + + present(); + + // read parameters from the stdin + + readparam(); + if ( verbose >= VERBOSE_MINIMAL ) + log( SIGNATURE, "... reading parameters OK\n"); + // get name of the output file + + strcpy( outname, get_output_filename() ); + + // open HBOOK file + + if ( verbose >= VERBOSE_MINIMAL ) + log( SIGNATURE, "Openning N-Tuple file %s\n", outname ); + + // TFile hfile(outname,"RECREATE","Check MC generated data"); + + + HROPEN(1, "CHECKMC", outname, "N", record_size, istat); + if (istat != 0) + {fprintf(stderr, "\n Output Histo couldn't be opened. Adios.\n"); + exit(999);} + + // profile (final) histograms + + // TH1F *htime = new TH1F("htime", "htime", 60, 0., 60.); + // TH1F *helec = new TH1F("helec", "helec", NTHSTEP, 0., 790.); + // TH1F *hlight = new TH1F("hlight", "hlight", NTHSTEP, 0., 790.); + + double fstime[18][60], fs2time[18][60]; + + double ftstime[18][60], fts2time[18][60]; + double ftselec[18][NTHSTEP], fts2elec[18][60]; + double ftslight[18][NTHSTEP], fts2light[18][60]; + + for (i=0; i<60; ++i) + for (k=0; k<18; ++k) + ftstime[k][i] = fts2time[k][i] = 0.0; + + for (i=0; i= VERBOSE_MINIMAL ) + log( SIGNATURE, "Openning directory %s\n", pathname ); + if ( (directory = opendir(pathname)) == NULL ) + error( SIGNATURE, + "Cannot open directory %s\n", pathname ); + // trace the number of events (cer* files) + num_cer_files = 0; + max_num_cer_files = get_max_events(); + // for each directory entry (files) + while ( ( (de = readdir( directory )) != NULL) && + (num_cer_files < max_num_cer_files) ) { + // skip removed files + if ( de->d_ino == 0) + continue; + // keep only cer* files + if ( strstr(de->d_name, "cer") != de->d_name ) + continue; + // it is a real cer* file + ++num_cer_files; + if ( verbose >= VERBOSE_NORMAL ) + log( SIGNATURE, "cerfile: %s\n", de->d_name ); + // get cer* file number (number of the shower) + // increments the pointer in 3 to calculate the number + ncer = atoi(de->d_name + 3); + // full names of the input files cer* and sta* + sprintf(cername, "%s/cer%06d", pathname, ncer); + sprintf(staname, "%s/sta%06d", pathname, ncer); + // try to open cer* file + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "Openning %s\n", cername ); + cerfile.open( cername ); + if ( cerfile.bad() ) + error( SIGNATURE, + "Cannot open input file: %s\n", cername ); + // try to open the statistics file + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "Openning %s\n", staname ); + stat.openfile ( staname ); + // reading data from this sta* file + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "Reading %s\n", staname ); + stat.read(); + // read the Event Header from the cer* file + evth.read( cerfile ); + if (evth.PrimaryID == gammas) + ++nevents[0]; + else if (evth.PrimaryID == protons) + ++nevents[1]; + else if (evth.PrimaryID == helium) + ++nevents[2]; + else + ++nevents[3]; + cerfile.close(); + stat.closefile();}} + + + cout<<"NMAXBINS"<0) +// Primary histos + {sprintf(htit, "Log(Energy) of primary %s;log(E/GeV)", primary_name[i]); + HBOOK1(100+i,htit,min(nevents[i]/25,NMAXBINS),0.,4.5,0.); + + sprintf(htit, "Core distance (horizontal plane) of primary %s;(m)", primary_name[i]); + HBOOK1(110+i,htit,min(nevents[i]/25,NMAXBINS),0.,500.,0.); + + sprintf(htit, "Theta of primary %s;(degrees)", primary_name[i]); + HBOOK1(120+i,htit,min(nevents[i]/15,NMAXBINS),0.,70.,0.); + + sprintf(htit, "Phi of primary %s;(degrees)", primary_name[i]); + HBOOK1(130+i,htit,min(nevents[i]/25,NMAXBINS),0.,360.,0.); + + sprintf(htit, "Core position for primary %s;x (m);y (m)", primary_name[i]); + HBOOK2(140+i,htit,50,-500.,500.,50,-500.,500.,0.); + + sprintf(htit, "Height of first interaction of primary %s;(km)", primary_name[i]); + HBOOK1(150+i,htit,min(nevents[i]/25,NMAXBINS),0.,100.,0.); + + sprintf(htit, "Time (last-first) of C-Photons for %s;(ns)", primary_name[i]); + HBOOK1(160+i,htit,min(nevents[i]/25,NMAXBINS),-1.,160.,0.); + + sprintf(htit, "Shower maximum XMAX for %s;(g cm^-2!)", primary_name[i]); + HBOOK1(170+i,htit,min(nevents[i]/10,NTHSTEP),0.,790.,0.); + + sprintf(htit, "Log(Num of C-photons produced) for primary %s below 300 GeV;log(N)", primary_name[i]); + HBOOK1(180+i,htit,min(nevents[i]/25,NMAXBINS),-.5,6.,0.); + + sprintf(htit, "Log(Num of C-photons produced) for primary %s above 300GeV;log(N)", primary_name[i]); + HBOOK1(190+i,htit,min(nevents[i]/25,NMAXBINS),-.5,6.,0.); + +// C-Photon histos + sprintf(htit, "C-Photon position in horiz. plane for prim. %s;x (m);y (m)", primary_name[i]); + HBOOK2(300+i,htit,50,-15.,15.,50,-15.,15.,0.); + + sprintf(htit, "C-Photon distance to (0,0) (horizontal plane) for %s;(m)", primary_name[i]); + HBOOK1(310+i,htit,(min(nevents[i]/5,NMAXBINS)),0.,15.,0.); + + sprintf(htit, "C-Photon spectrum for primary %s; [l] (nm)", primary_name[i]); + HBOOK1(320+i,htit,(min(nevents[i]/5,NMAXBINS)),200.,700.,0.); + + sprintf(htit, "Production height of C-Photon for primary %s;(km)", primary_name[i]); + HBOOK1(330+i,htit,(min(nevents[i]/5,NMAXBINS)),0.,30.,0.); + } + + ntotfiles = 0; + + // for each path of data files + + for (np=0; np= VERBOSE_MINIMAL ) + log( SIGNATURE, "Openning directory %s\n", pathname ); + if ( (directory = opendir(pathname)) == NULL ) + error( SIGNATURE, + "Cannot open directory %s\n", pathname ); + + // trace the number of events (cer* files) + + num_cer_files = 0; + + // book n-tuple for this directory + + // hid = np+1; + + // if ( verbose >= VERBOSE_MINIMAL ) + // log( SIGNATURE, "Create N-Tuple #%d for path %s . . .\n", + // hid, pathname ); + + // HBOOKN(hid, pathname, nvar, " ", 4096, chTags); + + max_num_cer_files = get_max_events(); + + // for each directory entry (files) + + firstevent=-100; + while ( ( (de = readdir( directory )) != NULL) && + (num_cer_files < max_num_cer_files) ) { + + // skip removed files + + if ( de->d_ino == 0) + continue; + + // keep only cer* files + + if ( strstr(de->d_name, "cer") != de->d_name ) + continue; + + // it is a real cer* file + + ++num_cer_files; + + if ( verbose >= VERBOSE_NORMAL ) + log( SIGNATURE, "cerfile: %s\n", de->d_name ); + + // get cer* file number (number of the shower) + + // increments the pointer in 3 to calculate the number + + ncer = atoi(de->d_name + 3); + + // full names of the input files cer* and sta* + + sprintf(cername, "%s/cer%06d", pathname, ncer); + sprintf(staname, "%s/sta%06d", pathname, ncer); + + // try to open cer* file + + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "Openning %s\n", cername ); + + cerfile.open( cername ); + + if ( cerfile.bad() ) + error( SIGNATURE, + "Cannot open input file: %s\n", cername ); + + // @endcode + // @text + // Each shower has associated three files: {\em Cherenkov-photons file\/} + // ({\tt cer*}), {\em Particles file} ({\tt dat*}) and + // {\em Statistics file} ({\tt sta*}). First we read the data from + // this last file, which is explained below, and then we read the + // {\em Event-Header} block from the {\tt cer*} file. + // @endtext + // @code + + // try to open the statistics file + + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "Openning %s\n", staname ); + + stat.openfile ( staname ); + + // reading data from this sta* file + + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "Reading %s\n", staname ); + + stat.read(); + + // read the Event Header from the cer* file + + evth.read( cerfile ); + + // Distinquish between gamma/proton/helium/+heavier nuclei + if (evth.PrimaryID == gammas) + hprimary=0; + else if (evth.PrimaryID == protons) + hprimary=1; + else if (evth.PrimaryID == helium) + hprimary=2; + else + hprimary=3; + + // set bin for theta/energy + + kbin = ((int)floor(evth.get_theta()/10.)*6 + + (int)floor(log10(evth.get_energy()))); + + // get core distance and position + + coreD = evth.get_core(&coreX, &coreY); + + t1 = 100000; + t2 = 0.0; + + ++ntotfiles; + + // initialize time-histogram + + for (i=0; i<60; ++i) + fstime[kbin][i] = 0.0; + + // read each and every single photon (particle) from cer* file + + rec2[0] = evth.Etotal; + rec2[1] = coreD; + rec2[2] = evth.Theta; + rec2[10] = evth.EvtNumber; + + if (firstevent == -100) + {firstevent = evth.EvtNumber; + // save random number seed of first event per run + for (i=0;i<3;++i) + for (j=0;j<3;++j) + rnd[i][j]=evth.RndData[i][j];} + + lastevent = evth.EvtNumber; + + nCphotons = 0; + while ( ! cerfile.eof() ) { + + // read photon data from the file + + photon.read( cerfile ); + + wl = photon.get_wl(); + + HF2(300+hprimary,(photon.x-coreX)/100.,(photon.y-coreY)/100.,1.); + HF1(310+hprimary,sqrt(pow(((photon.x-coreX)/100.),2.)+pow((photon.y-coreY)/100.,2.)),1.); + HF1(320+hprimary,wl,1.); + HF1(330+hprimary,photon.h/100000.,1.); + + rec2[3] = wl; + rec2[4] = photon.x; + rec2[5] = photon.y; + rec2[6] = photon.u; + rec2[7] = photon.v; + rec2[8] = photon.t; + rec2[9] = photon.h; +// HFN(hid2,rec2); + + + if ( wl < 1.0 ) + break; + + // fill times + + t = photon.get_t() - stat.get_tfirst(); + + if ( (t>0.0) && (t<60.0)) { + i = (int)t/1.; + ++fstime[kbin][i]; + } + + // hdmytime->Fill(t); + + // increase number of Cphotons written + + ++nCphotons; + + } // while still there are photons left + + // calculate errors for time distribution + + for (i=0;i<60;++i) + fs2time[kbin][i] = sqrt(fstime[kbin][i]); + + // save intermediate vectors in final vector + + for (i=0;i<60;++i) { + ftstime[kbin][i] += fstime[kbin][i]; + fts2time[kbin][i] += fs2time[kbin][i]; + } + + for (i=0;i vmax) { + vmax = v; + j = i; + } + } + + // to calculate a rough estimate of Xmax, we make a weigthed mean + // of the values at maximum, and the previous and next bins + xmax=-100.; + if (j > 0) { + x0 = j*10. + 5.; + xa = (j-1)*10. + 5.; + xb = (j+1)*10. + 5.; + vmaxa = (stat.plong[2][j-1] + stat.plong[1][j-1]); // e- & e+ + vmaxb = (stat.plong[2][j+1] + stat.plong[1][j+1]); // e- & e+ + + xmax = (xa*vmaxa + x0*vmax + xb*vmaxb) / (vmaxa + vmax + vmaxb);} + // else there was no e+e- shower in this event, skip ... + + // fill N-Tuple record + + if (evth.Etotal>0) + HF1(100+hprimary,log10(evth.Etotal),1.); + HF1(110+hprimary,coreD/100.,1.); + HF1(120+hprimary,evth.Theta/3.14*180.,1.); + HF1(130+hprimary,evth.Phi/3.14*180.,1.); + HF2(140+hprimary,coreX/100.,coreY/100.,1.); + HF1(150+hprimary,evth.zFirstInt/100000.,1.); + if (stat.get_tlast()-stat.get_tfirst()>0) + HF1(160+hprimary,stat.get_tlast()-stat.get_tfirst(),1.); + else + HF1(160+hprimary,-5.,1.); + + HF1(170+hprimary,xmax,1.); + if (evth.Etotal<300){ + if (nCphotons>0) + HF1(180+hprimary,log10(nCphotons),1.); + else + HF1(180+hprimary,-1.,1.); + } + else{ + if (nCphotons>0) + HF1(190+hprimary,log10(nCphotons),1.); + else + HF1(190+hprimary,-1.,1.); + } + rec[0] = evth.PrimaryID; + rec[1] = evth.Etotal; + rec[2] = coreD; + rec[3] = evth.Theta; + rec[4] = evth.Phi; + rec[5] = thick(evth.zFirstInt, evth.Theta); + rec[6] = evth.zFirstInt; + rec[7] = evth.FirstTarget; + rec[8] = evth.p[0]; + rec[9] = evth.p[1]; + rec[10] = evth.p[2]; + rec[11] = stat.get_tfirst(); + rec[12] = stat.get_tlast(); + rec[13] = xmax; + rec[14] = coreX; + rec[15] = coreY; + rec[16] = nCphotons; + + rec[17] = evth.RunNumber; + rec[18] = evth.EvtNumber; + rec[19] = evth.DateRun; + rec[20] = evth.VersionPGM; + rec[21] = evth.HeightLev[0]; + rec[22] = evth.SlopeSpec; + rec[23] = evth.ELowLim; + rec[24] = evth.EUppLim; + + + rec[25] = evth.Ecutoffh; + rec[26] = evth.Ecutoffm; + rec[27] = evth.Ecutoffe; + rec[28] = evth.Ecutoffg; + + rec[29] = evth.EGS4yn; + rec[30] = evth.NKGyn; + rec[31] = evth.GHEISHAyn; + rec[32] = evth.VENUSyn; + rec[33] = evth.CERENKOVyn; + rec[34] = evth.NEUTRINOyn; + rec[35] = evth.HORIZONTyn; + rec[36] = evth.COMPUTER; + rec[37] = evth.ThetaMin; + rec[38] = evth.ThetaMax; + rec[39] = evth.PhiMin; + rec[40] = evth.PhiMax; + + rec[41] = evth.StepLength; + rec[42] = evth.CWaveLower; + rec[43] = evth.CWaveUpper; + + // write record to the file + +// HFN(hid,rec); + + // ntuple->Fill( rec ); + + // close files + + cerfile.close(); + stat.closefile(); + + // show how many photons were written + + if ( verbose >= VERBOSE_MAXIMAL ) + log( SIGNATURE, "%d C-photons written.\n", nCphotons ); + + + } // while there are still directory entries left + + // check parameters for each RUN... + // Monitor run statistics: + + ++rndiag; + sprintf(ctmp, "%s\n", "-------------------------------------------------------------"); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + ++rndiag; + sprintf(ctmp,"%s%f%s%f%s%f\n", "Run: ",evth.RunNumber, + " Date: ",evth.DateRun, " CorsikaV: ",evth.VersionPGM); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + ++rndiag; + sprintf(ctmp, "%s%s\n", " Primary Particle: ", + primary_name[hprimary]); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + ++rndiag; + sprintf(ctmp, "%s%d%s%d\n", " First/Last Event Number monitored: ", + firstevent," / ",lastevent); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + ++rndiag; + sprintf(ctmp, "%s %f %f %f\n", " Random number seeds: ", + rnd[0][0],rnd[0][1],rnd[0][2]); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + ++rndiag; + sprintf(ctmp, "%s %f %f %f\n", " : ", + rnd[1][0],rnd[1][1],rnd[1][2]); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + ++rndiag; + sprintf(ctmp, "%s %f %f %f\n", " : ", + rnd[2][0],rnd[2][1],rnd[2][2]); + for (i=0;i<256;++i) rdiag[i][rndiag]=ctmp[i]; + + // check ranges of i/p params + // + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.HeightLev[0],mon_HeightLev, + ": Observation level is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.SlopeSpec,mon_SlopeSpec[hprimary], + ": Energy-spectrum slope is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.ELowLim,mon_ElowLim[hprimary], + ": Energy-low lim is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.EUppLim,mon_EUppLim[hprimary], + ": Energy-up lim is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.Ecutoffh,mon_Ecutoffh[hprimary], + ": Energy-cutoff hadrons is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.Ecutoffm,mon_Ecutoffm[hprimary], + ": Energy-cutoff muons is "); + + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.Ecutoffe,mon_Ecutoffe[hprimary], + ": Energy-cutoff electrons is "); + + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.Ecutoffg,mon_Ecutoffg[hprimary], + ": Energy-cutoff gammas is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.EGS4yn,mon_EGS4yn[hprimary], + ": EGS4yn flag is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.NKGyn,mon_NKGyn[hprimary], + ": NKGyn flag is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.GHEISHAyn,mon_GHEISHAyn[hprimary], + ": GHEISHAyn flag is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.VENUSyn,mon_VENUSyn[hprimary], + ": VENUSyn flag is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.CERENKOVyn,mon_CERENKOVyn[hprimary], + ": CERENKOVyn flag is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.ThetaMin,mon_ThetaMin[hprimary], + ": ThetaMin is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.ThetaMax,mon_ThetaMax[hprimary], + ": ThetaMax is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.PhiMin,mon_PhiMin[hprimary], + ": PhiMin is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.PhiMax,mon_PhiMax[hprimary], + ": PhiMax is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.CWaveLower,mon_CWaveLower[hprimary], + ": CWaveLower is "); + monitor_event(diag,&ndiag,primary_name[hprimary], + evth.CWaveUpper,mon_CWaveUpper[hprimary], + ": CWaveUpper is "); + + } // for each data directory + + // store contents of temporary 1D histograms to profiles + + /* + Int_t bin; + Stat_t value; + + for (kbin=0; kbin<18; ++kbin) { + + for (i=0; i<60; ++i) { + ftstime[kbin][i] /= ntotfiles; + fts2time[kbin][i] = sqrt(fts2time[kbin][i]/ntotfiles - + SQR(ftstime[kbin][i])); + bin=i+1, value=ftstime[kbin][i]; + cout << bin << ' ' << value << endl << flush; + htime->SetBinContent( bin, value ); + htime->SetBinError( i+1, fts2time[kbin][i] ); + } + + for (i=0; iSetBinContent( i+1, ftselec[kbin][i] ); + helec->SetBinError( i+1, fts2elec[kbin][i] ); + } + + for (i=0; iSetBinContent( i+1, ftslight[kbin][i] ); + hlight->SetBinError( i+1, fts2light[kbin][i] ); + } + + TH1F h1; + TH1F h2; + TH1F h3; + + htime->Copy( h1 ); + helec->Copy( h2 ); + hlight->Copy( h3 ); + + hfile.Write(); + + } + */ + + // close output file + + if ( verbose >= VERBOSE_MINIMAL ) + log( SIGNATURE, "Closing N-Tuple file %s\n", outname ); + + // hfile.Write(); + // hfile.Close(); + + HROUT(0,icycle," "); + + fprintf(stdout, "%s\n", "-------------------------------------------------------------"); + fprintf(stdout, "%s\n", "---------------------- RUN STATISTICS -----------------------"); + + for (i=0;i<=rndiag;++i) + {char tmp[256]; int j; + {for (j=0;j<256;++j) tmp[j]=rdiag[j][i];} + fprintf(stdout,"%s",tmp);} + + fprintf(stdout, "%s\n", "-------------------------------------------------------------"); + fprintf(stdout, "%s\n", "----------------------- DIAGNOSTICS -------------------------"); + fprintf(stdout, "%s\n", "-------------------------------------------------------------"); + + if (ndiag<0) + {fprintf(stdout, "%s\n", " None, everything's splendidly perfect ..."); + fprintf(stdout, "%s\n", " Nada, todo fue bien, hacemos fiesta ..."); + fprintf(stdout, "%s\n", " Keine Fehler gefunden ...");} + else + {for (i=0;i<=ndiag;++i) + {char tmp[256]; int j; + {for (j=0;j<256;++j) tmp[j]=diag[j][i];} + fprintf(stdout,"%s",tmp);}} + + + HREND("CHECKMC"); + + // program finished + + if ( verbose >= VERBOSE_MINIMAL ) + log( SIGNATURE, "Done.\n"); + +} +// @endcode + +// @T \newpage + +// @subsection Functions definition + +//------------------------------------------------------------ +// @name present +// +// @desc Make some presentation +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//------------------------------------------------------------ +// present +// +// make some presentation +//------------------------------------------------------------ + +void +present(void) +{ +cout << "##################################################\n" +<< SIGNATURE << '\n' << '\n' +<< "Check of the MC data generated by MMCS\n" +<< "J C Gonzalez, P Morawitz, D Petry, 27 Oct 1998\n" +<< "##################################################\n\n"; +} +// @endcode + +//------------------------------------------------------------ +// @name log +// +// @desc function to send log information +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//------------------------------------------------------------ +// log +// +// function to send log information +//------------------------------------------------------------ + +void +log(const char *funct, char *fmt, ...) +{ + va_list args; + + // Display the name of the function that called error + printf("[%s]: ", funct); + + // Display the remainder of the message + va_start(args, fmt); + vprintf(fmt, args); + va_end(args); +} +// @endcode + +//------------------------------------------------------------ +// @name error +// +// @desc function to send an error message, and abort the program +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//------------------------------------------------------------ +// error +// +// function to send an error message, and abort the program +//------------------------------------------------------------ + +void +error(const char *funct, char *fmt, ...) +{ + va_list args; + + // Display the name of the function that called error + fprintf(stderr, "ERROR in %s: ", funct); + + // Display the remainder of the message + va_start(args, fmt); + vfprintf(stderr, fmt, args); + va_end(args); + + fprintf(stderr, "\nTerminating program.\n"); + exit(999); +} +// @endcode + + +//------------------------------------------------------------ +// @name readparam +// +// @desc function to read the parameters of the check +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//------------------------------------------------------------ +// readparam +// +// function to read the parameters of the check +//------------------------------------------------------------ + +void +readparam(void) +{ + char sign[] = GLUE_postp( PROGRAM, VERSION ); // initialize sign + char line[LINE_MAX_LENGTH]; // line to get from the stdin + char token[ITEM_MAX_LENGTH]; // a single token + int i, j; // dummy counters + + // get signature + cout << "\nExpecting input from stdin (see example input file) ...\n"; + cin.getline(sign, strlen(SIGNATURE) + 1); + // cin.getline(line, LINE_MAX_LENGTH); + if (strcmp(sign, SIGNATURE) != 0) { + cerr << "ERROR: Signature of parameters file is not correct\n"; + cerr << '"' << sign << '"' << '\n'; + cerr << "should be: " << SIGNATURE << '\n'; + exit(1); + } + + // loop till the "end" directive is reached + int is_end = FALSE; + while (! is_end) { + + // get line from stdin + cin.getline(line, LINE_MAX_LENGTH); + + // look for each item at the beginning of the line + for (i=0; i<=end_file; i++) + if (strstr(line, PAR_ITEM_NAMES[i]) == line) + break; + + // if it is not a valid line, just ignore it + if (i == end_file+1) { + cerr << " Skipping unknown token in [" << line << "]\n"; + continue; + } + + // case block for each directive + switch ( i ) { + + case data_paths: // beginning of the list of valid paths + + //get number of paths to read + sscanf(line, "%s %d", token, &Num_of_paths); + + if ( verbose == TRUE ) + cout << " " << Num_of_paths << " path(s) will be used.\n"; + + // allocate memory for paths list + Paths_list = (char**)malloc(Num_of_paths * sizeof(char*)); + for (i=0; i thickness [g cm-2] +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//-------------------------------------------------- +// thick +// +// converts height [cm] -> thickness [g cm-2] +//-------------------------------------------------- + +float +thick( float height, float theta ) +{ + float h, thk; + + // Ibarra's modification + h = sqrt( SQR(REarth) + + SQR(height/cos(theta)) + + (2.0*REarth*height)) - REarth; + + if ( h < 4.e5 ) { + thk = AATM[0] + BATM[0] * exp ( -h * DATM[0] ); + } else if ( h < 1.e6 ) { + thk = AATM[1] + BATM[1] * exp ( -h * DATM[1] ); + } else if ( h < 4.e6 ) { + thk = AATM[2] + BATM[2] * exp ( -h * DATM[2] ); + } else if ( h < 1.e7 ) { + thk = AATM[3] + BATM[3] * exp ( -h * DATM[3] ); + } else { + thk = AATM[4] - h * CATM[4]; + } + + return ( thk ); + +} +// @endcode + +//------------------------------------------------------------ +// @name height +// +// @desc converts thickness [g cm-2] -> height [cm] +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//-------------------------------------------------- +// height +// +// converts thickness [g cm-2] -> height [cm] +//-------------------------------------------------- + +float +height( float thk, float theta ) +{ + float h, height; + + if ( thk > 631.1e0 ) { + h = CATM[0] * log ( BATM[0] / (thk - AATM[0]) ); + } else if ( thk > 271.7e0 ) { + h = CATM[1] * log ( BATM[1] / (thk - AATM[1]) ); + } else if ( thk > 3.0395e0 ) { + h = CATM[2] * log ( BATM[2] / (thk - AATM[2]) ); + } else if ( thk > 0.00128292e0 ) { + h = CATM[3] * log ( BATM[3] / (thk - AATM[3]) ); + } else { + h = (AATM[4] - thk) * DATM[4]; + } + + // Ibarra's modification + height = SQR(cos(theta)) * + (-REarth + sqrt(SQR(REarth) + + (SQR(h + (2.0*REarth*h))/SQR(cos(theta))))); + + return ( height ); + +} +// @endcode + +//------------------------------------------------------------ +// @name monitor_event +// +// @desc compare parameter with reference and report if not equal +// +// @date Sat Jun 27 05:58:56 MET DST 1998 +// @function @code +//-------------------------------------------------- +// monitor_event +// +// compare parameter with reference +//-------------------------------------------------- +void monitor_event(char diag[DIAGX][DIAGY], int *ndiag, const char* primary, + const float a, const float b, const char* message) +{char diagt[DIAGX]; int i; + if (a!=b){ + if (*ndiag<(DIAGY-1)) + {++*ndiag; + sprintf(diagt,"%s%s%s%s%f%s%f\n","Primary ",primary," ", + message,a," but should be ",b);} + else{ + sprintf(diagt,"%s", + "Two many errors encountered, won't print any further diagnostics\n"); + for (i=0;i +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include "jcmacros.h" +#include "jcdebug.h" + +#include "COREventHeader.hxx" +#include "CORParticle.hxx" +#include "CORStatfile.hxx" + +// @endcode + +// @subsection Macro-definitions, and constants + +// @code + + + +#ifndef VERSION + +#define PROGRAM checkmc +#define VERSION 1.0 + +#define GLUE_prep(x,y) #x" "#y +#define GLUE_postp(x,y) GLUE_prep(x,y) + +const char SIGNATURE[] = GLUE_postp( PROGRAM, VERSION ); + +#endif // ! VERSION + + + +#define NORM(x) (sqrt(SQR(x[0])+SQR(x[1])+SQR(x[2]))) /* Norm(vector) */ +#define RAD(x) ((x)*0.0174532925199433) +#define DEG(x) ((x)*57.2957795130823) + + + +// now we define the list PAR_ITEM_LIST of possible items in +// the CT definition file. + +#define PAR_ITEM_LIST /* LIST OF ITEMS IN THE CT DEFINITION FILE */ \ +T(data_paths), /* begin of the list of source data paths */ \ +T(output_file), /* output file */ \ +T(verbose_level), /* defines verbose level of the output */ \ +T(max_events), /* maximum number of event to read */ \ +T(reference_histo), /* reference histogram */ \ +T(end_file) /* end of the parameters file */ + +#define T(x) x // define T() as the name as it is + +enum PAR_ITEM_TYPE { + PAR_ITEM_LIST +}; + +#undef T + +#define T(x) #x // define T() as the string of x + +const char *const PAR_ITEM_NAMES[] = { + PAR_ITEM_LIST +}; + +#undef T + + +#define LINE_MAX_LENGTH 120 +#define ITEM_MAX_LENGTH 40 +#define PATH_MAX_LENGTH 120 + +#define NTHSTEP 79 // assumes thstep = 10 g/cm2 (790/10 = 79) + +// maximum number of bins for the 1D histograms +#define NMAXBINS 100 + + +// Verbose Levels +enum VerboseLevel { + VERBOSE_QUIET, + VERBOSE_MINIMAL, + VERBOSE_NORMAL, + VERBOSE_MAXIMAL +}; + +#define VERBOSE_DEFAULT VERBOSE_MINIMAL + +// @endcode + +/* @text +Definition of the dimensions of the diagnostics messages array in + monitor_event() +@endtext */ +// @code +#define DIAGX 256 +#define DIAGY 20 +// @endcode + + +/* @text +Data for the atmosphere: useful for conversion between height and thickness. +@endtext */ + +// @code + +float AATM[5] = { + -186.5562e0, + -94.919e0, + 0.61289e0, + 0.e0, + .01128292e0 }; +float BATM[5] = { + 1222.6562e0, + 1144.9069e0, + 1305.5948e0, + 540.1778e0, + 0.e0 }; +float CATM[5] = { + 994186.38e0, + 878153.55e0, + 636143.04e0, + 772170.16e0, + 1.e-9}; +float DATM[5] = { + 1.00584761581626e-06, + 1.13875301193054e-06, + 1.57197349828743e-06, + 1.29505133946124e-06, + 1.00000000000000e+09}; // = 1 / CATM + +// @endcode + + +// @subsection Prototypes of functions + +//++ +// prototypes +//-- + +//@code + +void present(void); +void log(const char *funct, char *fmt, ...); +void error(const char *funct, char *fmt, ...); +void readparam(void); +int get_num_of_paths(void); +char *get_path_name(int i); +char *get_output_filename(void); +int get_max_events(void); +float thick( float height, float theta ); +float height( float thk, float theta ); +void monitor_event(char diag[DIAGX][DIAGY], int *ndiag, const char* primary, + const float a, const float b, const char* message); + +// @endcode + +// @text +// CFortran routines will be used. +// @endtext + +// @code +#include +#include +// #include +#include "hbook.h" +//#include "minuit.h" +//#include "kuip.h" +//#include "zebra.h" + +#define PAWC_SIZE 1000000 + +float pawc_[PAWC_SIZE]; + +// extern "C" +//{ +// extern void hlimit_(int*); +// extern void hropen_(int*, char*, char*, char*, int*, int*, int, int, int); +// extern void hbookn_(int*, char*, int*, char*, int*, char*, int, int, int); +// extern void hrout_(int*, int*, char*, int); +// extern void hrend_(char*, int); +// extern void hfn_(int*, float*); +// } +// @endcode + + + char* strconcat(const char* s1, const char* s2, char* s3) + { + for (int i=0;s3[i]=s1[i]; ++i); + for (int i=strlen(s1);s3[i]=s2[i-strlen(s1)]; ++i); + return s3; + } + + + Index: trunk/MagicSoft/Simulation/Corsika/Checkmc/config.mk.linux-gnu =================================================================== --- trunk/MagicSoft/Simulation/Corsika/Checkmc/config.mk.linux-gnu (revision 328) +++ trunk/MagicSoft/Simulation/Corsika/Checkmc/config.mk.linux-gnu (revision 328) @@ -0,0 +1,54 @@ +################################################################## +# +# config.mk +# +# @file config.mk +# @title small configuration file for Makefile +# @author J C Gonz\'alez +# @email gonzalez@mppmu.mpg.de +# @date Fri Mar 12 11:51:11 MET 1999 +# +#_______________________________________________________________ +# +# Created: Fri Mar 12 11:51:11 MET 1999 +# Author: Jose Carlos Gonzalez +# Purpose: Makefile for the compilation of the reflector program +# Notes: +# +#--------------------------------------------------------------- +# $RCSfile: config.mk.linux-gnu,v $ +# $Revision: 1.1.1.1 $ +# $Author: magiccvs $ +# $Date: 1999-11-25 20:34:58 $ +################################################################## +# @maintitle + +# @code + +# program + +PROGRAM = camera + +# compilers & tools + +CC = gcc +CXX = g++ +F77 = g77 + +DOCUM = ../sus/sus +RATE = ../../comstat/comstat -g -p + +OPTIM = -Df2cFortran +DEBUG = -g + +# libraries + +CERNDIR = /cern +INCLUDE_CERN = ${CERNDIR}/pro/include/cfortran + +# system + +SYSTEM = linux + +# @endcode +##EOF Index: trunk/MagicSoft/Simulation/Corsika/Checkmc/config.mk.osf1 =================================================================== --- trunk/MagicSoft/Simulation/Corsika/Checkmc/config.mk.osf1 (revision 328) +++ trunk/MagicSoft/Simulation/Corsika/Checkmc/config.mk.osf1 (revision 328) @@ -0,0 +1,54 @@ +################################################################## +# +# config.mk +# +# @file config.mk +# @title small configuration file for Makefile +# @author J C Gonz\'alez +# @email gonzalez@mppmu.mpg.de +# @date Fri Mar 12 11:51:11 MET 1999 +# +#_______________________________________________________________ +# +# Created: Fri Mar 12 11:51:11 MET 1999 +# Author: Jose Carlos Gonzalez +# Purpose: Makefile for the compilation of the reflector program +# Notes: +# +#--------------------------------------------------------------- +# $RCSfile: config.mk.osf1,v $ +# $Revision: 1.1.1.1 $ +# $Author: magiccvs $ +# $Date: 1999-11-25 20:34:58 $ +################################################################## +# @maintitle + +# @code + +# program + +PROGRAM = camera + +# compilers & tools + +CC = cc +CXX = cxx +F77 = f77 + +DOCUM = ../sus/sus +RATE = ../../comstat/comstat -g -p + +OPTIM = -O +DEBUG = -g + +# libraries + +CERNDIR = /CERN +INCLUDE_CERN = ${CERNDIR}/cfortran + +# system + +SYSTEM = osf + +# @endcode +##EOF Index: trunk/MagicSoft/Simulation/Corsika/Checkmc/input_example =================================================================== --- trunk/MagicSoft/Simulation/Corsika/Checkmc/input_example (revision 328) +++ trunk/MagicSoft/Simulation/Corsika/Checkmc/input_example (revision 328) @@ -0,0 +1,8 @@ +checkmc 1.0 +data_paths 1 +/hd61/Maggi/Data/mmcs_spez_221001 +output_file test/checkmc.hbook +verbose_level 1 +max_events 10000 +end_file +