Changeset 7453 for trunk/MagicSoft/Mars/datacenter/scripts/linkmc

Timestamp:

12/12/05 16:30:50 (19 years ago)

Author:

Daniela Dorner

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/datacenter/scripts/linkmc (modified) (9 diffs)

trunk/MagicSoft/Mars/datacenter/scripts/linkmc

@@ -53 +53 @@
 mccampath=/montecarlo/camera
 mcpath=/montecarlo/rawfiles
-next=$mcpath/.next
-processed=$mcpath/.processed
-processeddir=$mcpath/.processed.dir
-readme=$mcpath/README.txt
+next=$mcpath/.next #in .next the next runno is stored
+processed=$mcpath/.processed #in .processed the linked files are stored
+readme=$mcpath/README.txt #file in which the information about the modes is redirected to have always an updated explanation
 
 if ! ls $next >> $scriptlog 2>&1
@@ -68 +67 @@
 
 
-modes=("" "Gammawobble+" "Gammanowobble0" "GammawobbleHE+" "GammanowobbleHE0" "Gammawobble0" "GammawobbleHE0" "Proton" "DiffuseGamma" )
+modes=("" "Gammawobble+" "Gammanowobble0" "GammawobbleHE+" "GammanowobbleHE0" "Gammawobble0" "GammawobbleHE0" "Gammadiffuse0" "Protonnowobble0" )
+
+#be carful: 
+# w- not yet foreseen in this script
 
 echo "" >> $readme 2>&1
@@ -101 +103 @@
 
 
+#get runnumber
 runno=`cat $next`
 
+#get files, which have to be linked
 camfiles=`find $mccampath -type f | grep -v Cal_and_Ped`
 
 for camfile in ${camfiles[@]}
 do 
+   #continue, if file is already linked 
    if grep $camfile $processed >> $scriptlog 2>&1
    then 
       continue
    fi
-   file=`basename $camfile`
-   no=`printf %08d $runno | cut -c 0-5`
-#   no2=`printf %08d $runno`
+   file=`basename $camfile` #filename
+   no=`printf %08d $runno | cut -c 0-5` #first 5 digits of a 8digit runno -> for path
+#   no2=`printf %08d $runno` #runno with 8 digits
 # workaround due to 5digit runnumber for data with runnumber < 35487
    no2=`printf %05d $runno`
-   zbin=`echo $file | cut -d_ -f2 | cut -c 5-6`
-   zbin=`printf %02d $zbin`
-   wobble=`echo $file | cut -d_ -f6 | cut -c 2`
-   particle=`echo $file | cut -d_ -f1`
-   psf=`echo $camfile | cut -d/ -f5 | cut -c 6,8`
-   
-   particledir=`echo $camfile | cut -d/ -f4`
-   wobbledir=`echo $camfile | cut -d/ -f6`
-   
+   zbin=`echo $file | cut -d_ -f2 | cut -c 5-6` #zbin from filename 
+   zbin=`printf %02d $zbin` #2digit
+   wobble=`echo $file | cut -d_ -f6 | cut -c 2` #mode from filename
+   particle=`echo $file | cut -d_ -f1` #particle type from filename
+   psf=`echo $camfile | cut -d/ -f5 | cut -c 6,8` #psf from path
+   
+   particledir=`echo $camfile | cut -d/ -f4` #particletype from path
+   wobbledir=`echo $camfile | cut -d/ -f6` #mode from path
+   
+   #build mode name
    testmode=$particle$wobbledir$wobble 
-#   echo "testmode:"$testmode >> $scriptlog 2>&1
-
+
+   #get mode extension for filename
    case $wobble in 
        0)   wobblemode="";;
@@ -136 +142 @@
    esac
    
+   #get no of mode from array $modes
    for (( i=1 ; i <= 12 ; i++ )) 
    do 
        if [ "${modes[$i]}" == "$testmode" ]
        then 
-#          echo "modes[$i]=${modes[$i]}" >> $scriptlog 2>&1
           totalmode=$i
           break
@@ -146 +152 @@
    done
 
-   totalmode=`printf %02d $totalmode`
-#   echo "totalmode: "$totalmode >> $scriptlog 2>&1
-
-   project=${particle}${zbin}${wobblemode}
+   totalmode=`printf %02d $totalmode` #2 digits
+
+   project=${particle}${zbin}${wobblemode} #build project name
+   #create new filename
    newfile="$mcpath/19$zbin/$totalmode/$psf/19${zbin}${totalmode}${psf}_${no2}_D_${project}_E.root"
    newdir=`dirname $newfile`
    makedir $newdir >> $scriptlog 2>&1
 
-   runno=`expr $runno + 1`
+   runno=`expr $runno + 1` #next runnumber
    echo $runno >| $next 
 
-   ln -sv $camfile $newfile >> $scriptlog 2>&1
+   #link file
+   ln -sv $camfile $newfile >> $scriptlog 2>&1 
+   #add filename to processed file
    echo $camfile >> $processed
 
@@ -163 +171 @@
 
 echo "linking cal and ped file" >> $scriptlog 2>&1
+#get files
 pedfile=`find $mccampath/Cal_and_Ped -name *_P_*.root`
 calfile=`find $mccampath/Cal_and_Ped -name *_C_*.root`
 echo "calfile"$calfile >> $scriptlog 2>&1
 echo "pedfile"$pedfile >> $scriptlog 2>&1
+#check number of files
 numfiles=`echo $pedfile $calfile | wc -w`
 if [ "$numfiles" != "2" ]
@@ -175 +185 @@
 fi
 
+#get all directories in the linked structure
 dirs=`find $mcpath -type d`
 
 for dir in ${dirs[@]}
 do 
-   if grep $dir $processeddir >> $scriptlog 2>&1
+   #continue, if directory has already linked C and P run
+   $cont=`ls $dir/*_0000[12]_[CP]_MonteCarlo_E.root | wc -w`
+   if [ "$cont" == "2" ] >> $scriptlog 2>&1
    then 
       continue
    fi
 
+   #continue, if directory is not at the lowest level of the structure
    cont=`echo $dir | cut -d/ -f6`
    if [ "$cont" == "" ]
@@ -190 +204 @@
    fi
 
+   #get date for filename from directory name
    date=`echo $dir | cut -c 22-25,27,28,30,31`
+   
+   #create new filenames and link files
 # workaround due to 5digit runnumber for data with runnumber < 35487
 #   newcalfile="${dir}/${date}_00000001_C_MonteCarlo_E.root"
@@ -199 +216 @@
    newpedfile="${dir}/${date}_00002_P_MonteCarlo_E.root"
    ln -sv $pedfile $newpedfile >> $scriptlog 2>&1
-   echo $dir >> $processeddir
 done