Changeset 9585 for trunk/MagicSoft/Mars/datacenter/scripts/runcorsika

Timestamp:

06/11/10 17:36:11 (15 years ago)

Author:

Daniela Dorner

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/datacenter/scripts/runcorsika (modified) (8 diffs)

trunk/MagicSoft/Mars/datacenter/scripts/runcorsika

@@ -18 +18 @@
 #
 #
-#   Author(s): Daniela Dorner  05/2007 <mailto:dorner@astro.uni-wuerzburg.de>
-#
-#   Copyright: MAGIC Software Development, 2000-2008
+#   Author(s): Daniela Dorner  05/2007 <mailto:daniela.dorner@unige.ch>
+#
+#   Copyright: MAGIC Software Development, 2000-2010
 #
 #
@@ -31 +31 @@
 printprocesslog "INFO starting $0"
 program=corsika
-column=fCorsikaFileAvail
+step=Corsika
 
 set -C
@@ -41 +41 @@
 {
    # Run number
-   echo "RUNNR "$corsikarunno
+   echo "RUNNR "$file
    # Starting event number
    echo "EVTNR 1"
@@ -60 +60 @@
    echo "FIXCHI "$startalt
    # Observation level, magnetic field, atmosphere
-   echo "OBSLEV 2200.E2"
+   echo "OBSLEV "$obslev".E2"
    # From www.noaa.gov/geomagmodels roughly for January 2009. (was 29.5 and 23.0)
    # MAGNET BX BZ (BX is the horizontal component (muT) to the x-direction of
    # the detector, BZ is the vertical component (muT) downwards)
-   echo "MAGNET 30.3 24.1"
+   echo "MAGNET "$magnet
    # Rotation angle between array x-direction and magnetic north direction
    # positive if array x-direction points to the west
-   echo "ARRANG -7.0"
+   echo "ARRANG "$magnetrot
    # Atmosphere (model, true/false for atmospheric refraction of cer-photons)
    # Atmosphere read from atmprof11.dat
    # My current understanding is that only refraction (bending of light
    # but no absorption of photons is taken into account) setting T.
-   echo "ATMOSPHERE 11 T"
+   echo "ATMOSPHERE "$atm" T"
+
+   # Current boundaries in reflector program
+   # echo "ATMLAY 400000 1000000 4000000 10000000"
+
+   # Boundaries as determined by CORSIKA
+   # Until we can read the result from the corsika fit we fix the
+   # values which come out of the fit. This is not ideal but forces
+   # things to be consistent. Be careful is using a different model -
+   # the fit might not be ideal.
+   echo "ATMLAY 775000 1650000 5000000 10500000"
+
    # Other parameters
    echo "RADNKG 200.E2"
@@ -82 +93 @@
    # Enable muon multiple scattering
    echo "MUMULT T"
-   # Cherenkov telescope options
-   echo "CWAVLG 290. 600."
-   echo "CSCAT "$reuse" 0. "$impactmax
+   # Wavelength range (was originally 600)
+   echo "CWAVLG "$wavelengths
+   # MMCS
+   echo "CSCAT 1 0. "$impactmax
+   # CORSIKA
+   # echo "CSCAT "$reuse" "$impactmax" 0."
    echo "CERSIZ 1."
+   # Write Cherenkov photons to cer-file
    echo "CERFIL T"
+   # Number of telescopes
    echo "CERTEL 1"
-   echo " 0. 0. 0. 0. 0. "$diameter" 1700."
+   # MMCS
+   echo " 0. 0. 0. 0. 0. "$diameter" "$diameter
+   # CORSIKA (x, y, z, r)
+   # echo "TELESCOPE 0. 0. 0. "$diameter
+   # This keyword has no relevance but allows to store the 
+   # maximum simulated impact parameter. The "magic" is that
+   # ceres interpretes the third and fourth argument as
+   # maximum simulated impact if they are identical and
+   # the other four arguments are 1
+   # echo "CERARY 1 1 "$diameter" "$diameter" 1 1"
    # How the longitudinal shower development is sampled
    # echo "LONGI T 10. T F"
@@ -103 +128 @@
    # Output path
    echo "DIRECT "$outpath
-   # User and host (for convinience)
+   # DIRECT /dev/null
+   # TELFIL filename
+   # TELFIL | gzip
+   # TELESCOPE X Y Z R
+   # User and host (for convenience)
    echo "USER "`whoami`
    echo "HOST "`hostname`
@@ -114 +143 @@
 gettodo "1" 
 run=${primaries[0]}
-# get corsika runno
-query="SELECT fCorsikaRunNumber FROM MCRunData where fMCRunNumber="$run
-# to be changed later again, when additional tables are available
-#corsikarunno=`sendquery`
-corsikarunno=$run
+file=${primaries[1]}
 
 # lock sequ
-lockfile=$lockpath/lock-$table-$column-$corsikarunno.txt
+lockfile=$lockpath/lock-$step-$run.$file.txt
 checklock 
 
-printprocesslog "INFO starting $program for run $run"
+printprocesslog "INFO starting $program for run $run and file $file"
 
 setstatus "start" 
 
-printprocesslog "INFO run $program for run $run " 
-printprocesslog "INFO create input card... " 
+printprocesslog "INFO run $program for run $run  and file $file" 
+printprocesslog "INFO create input card for run $run  and file $file" 
 
 # get values for inputcard
-query="SELECT fParticleTypeKEY FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fParticleTypeKEY FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 primparticle=`sendquery`
-query="SELECT fEnergyMin, fEnergyMax FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fEnergyMin, fEnergyMax FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 erange=`sendquery`
-query="SELECT fNumEvents FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fNumEvents FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 numevts=`sendquery`
-query="SELECT fEnergySlope FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fEnergySlope FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 slope=`sendquery`
-query="SELECT fZenithDistanceMin, fZenithDistanceMax FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fZenithDistanceMin, fZenithDistanceMax FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 zdrange=`sendquery`
 # Change from the MARS definition (Az=90 -> East) to Corsika (Az=90 -> West)
-query="SELECT 180-fAzimuthMin, 180-fAzimuthMax FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT 180-fAzimuthMin, 180-fAzimuthMax FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 azrange=`sendquery`
-query="SELECT DATE_FORMAT(fRunStart, '%Y/%m/%d') FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
-date=`sendquery`
-outpath=$mcpath"/"$program"/"$date
-query="SELECT fCorsikaSeed FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fCorsikaSeed1 FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 seed1=`sendquery`
-query="SELECT fCorsikaSeed+1 FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fCorsikaSeed2 FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 seed2=`sendquery`
-query="SELECT fCorsikaSeed+2 FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fCorsikaSeed3 FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 seed3=`sendquery`
-query="SELECT fImpactMax FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fImpactMax FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 impactmax=`sendquery`
-query="SELECT fViewConeMax FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fViewConeMax FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 viewconemax=`sendquery`
-query="SELECT fNumReUseShower FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
-reuse=`sendquery`
-query="SELECT fStartingAltitude FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fStartingAltitude FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 startalt=`sendquery`
-query="SELECT fMirrorDiameter FROM MCCorsikaRunData where fCorsikaRunNumber="$corsikarunno
+query="SELECT fMirrorDiameter FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
 diameter=`sendquery`
-
-logfile=$outpath/$program"-"`printf %06d $corsikarunno`".log"
+query="SELECT fAtmosphericModelKEY FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
+atm=`sendquery`
+query="SELECT fObsLevel FROM CorsikaInfo LEFT JOIN Observatory USING(fObservatoryKEY) WHERE fRunNumber=$run AND fFileNumber=$file"
+obslev=`sendquery`
+query="SELECT fMagnetBX, fMagnetBZ FROM CorsikaInfo LEFT JOIN Observatory USING(fObservatoryKEY) WHERE fRunNumber=$run AND fFileNumber=$file"
+magnet=`sendquery`
+query="SELECT fMagnetRotation FROM CorsikaInfo LEFT JOIN Observatory USING(fObservatoryKEY) WHERE fRunNumber=$run AND fFileNumber=$file"
+magnetrot=`sendquery`
+query="SELECT fWavelengthMin, fWavelengthMax FROM CorsikaInfo WHERE fRunNumber=$run AND fFileNumber=$file"
+wavelengths=`sendquery`
+
+outpath=$mcpath"/"$program"/"`printf %08d $run | cut -c 0-4`/`printf %08d $run | cut -c 5-8`
 makedir $outpath 
+logfile=$outpath/cer000`printf %06d $file | cut -c 4-6`".log"
 
 cd $corsikapath
 
-echo "INPUTCARD:" >| $logfile 2>&1
-printinputcard >> $logfile 2>&1
-echo "" >> $logfile 2>&1
-echo "" >> $logfile 2>&1
-printinputcard | ./cc6501p-linux >> $logfile 2>&1
-
-rm -f $outpath/dat`printf %06d $corsikarunno`
+printinputcard | ./cc6501p-linux >| $logfile 2>&1
 
 check1=$?
-
 case $check1 in
-   0)   printprocesslog "INFO $program finished successfully for corsika run number $corsikarunno (check1=$check1)"
+   0)   printprocesslog "INFO $program finished successfully for run $run file $file (check1=$check1)"
         ;;
-   *)   printprocesslog "ERROR $program failed for corsika run number $corsikarunno (check1=$check1)"
+   *)   printprocesslog "ERROR $program failed for run $run file $file (check1=$check1)"
         com=$Fcorsika
         check=$check1
@@ -189 +213 @@
 esac
 
+# removing empty datnnnnnn file
+rm -f $outpath/dat`printf %06d $file`
+
 setstatus "stop"