#!/usr/bin/perl
$rc = `fgrep "COMP = KELLY" ~/reacsource/NORED/*.dat`;
@rc=split/\n/,$rc;
print "$#rc\n";
$last="";
$REAC[1]="PI\\+ P --> PI\\+ P";
$REAC[2]="PI- P --> PI- P";
$REAC[3]="PI- P --> PI0 N";
$OBS[1]="POL";
$OBS[2]="SIG";
$OBS[3]="DSIG/DOMEGA";
foreach $r (1...3){
foreach $o (1...3){
$dir="output$r$o";
$html="/var/www/docs/kelly$r$o.html";
#$mk = `mkdir $dir`;
mkdir $dir;
open OUTPUT,">$dir/list";
open HTML,">$html";
print HTML "\n";
print HTML "\n";
print HTML "kelly11, \n";
print HTML "kelly12, \n";
print HTML "kelly13
\n";
print HTML "kelly21, \n";
print HTML "kelly22, \n";
print HTML "kelly23\n";
print HTML "kelly31, \n";
print HTML "kelly32, \n";
print HTML "kelly33\n";
print HTML "$REAC[$r]
\n";
print HTML "$OBS[$o]
\n";
print HTML "\n";
print HTML "| File | NAME | REFERENCE |
|---|
\n";
foreach $n (1...$#rc){
($file,$dummy)=split/:/,$rc[$n];
@name=split/\//,$file;
$name=$name[$#name];
($shortname,$dummy)=split/\./,$name;
open INPUT,"<$file";
$wantre="no";
$wantob="no";
while ($line=){
if($line=~m/^OBS/){
$line=~s/\(//g;
$line=~s/\)//g;
}
if($line=~m/RE = $REAC[$r]/) { $wantre="yes";}
if($line=~m/OBS = $OBS[$o]/) { $wantob="yes";}
}
close INPUT;
$wantre="yes";
$wantob="yes";
if($wantre eq "yes" && $wantob eq "yes"){
print OUTPUT "\n$n $name ";
open INPUT,"<$file";
$dataset=0;
while($line=){
chomp $line;
$oline=$line;
$line=~s/;//g;
if($line=~m/^SC =/){
$line=~s/SC = //;
print OUTPUT "$line";
$sc = $line;
print HTML "| ($n) $name | $line | \n";
}
elsif($line=~m/^R =/){
($ref,$type,$date)=split/\s*;\s*/,$oline;
$ref=~s/R = //;
$type=~s/TY = //;
$date=~s/D = //;
print OUTPUT "\t$ref\t$type\t$date";
print HTML "$ref | \n";
$ref =~ s/,\s+/,/g;
$db=`echo "select p.InspireId,p.RedId,p.HepdataId \
from Papers p, PaperRefs pr \
where p.PAPER_ID=pr.PAPER_ID \
and pr.Description like '$ref'" | mysql hepdata`;
undef @db;
@db=split/\n/,$db;
if($#db==1) {
print OUTPUT "\t$#db\t$db[1]";
print HTML "$db[1] | ";
}
else {
print HTML "- | ";
}
print OUTPUT "\n";
}
elsif($line=~m/^CR =/){
$comment=$line;
until ( $line eq "TABLE;" || $line eq "SYSTEMATICS;") {
$line=;
chomp $line;
if($line ne "TABLE;" && $line ne "SYSTEMATICS;") {
$comment=$comment.$line;
}
}
@comment=split/\/\/\//,$comment;
foreach $n (1...$#comment){
# print"($n) $comment[$n]\n";
}
}
elsif($line=~m/^RE =/){
($dummy,$reac)=split/\s*=\s*/,$line;
print OUTPUT "$reac";
}
elsif($line=~m/^OBS =/){
($dummy,$obs)=split/\s*=\s*/,$line;
print OUTPUT "\t$obs";
}
elsif($line=~m/^PLAB =/){
($dummy,$plab)=split/\s*=\s*/,$line;
($lower,$upper)=split/,/,$plab;
if($upper == $lower) {
$lower=1000*$lower;
printf OUTPUT "\t%0.1f MeV/c",$lower;
$dataset++;
# printf HTML "%0.1f | ",$lower;
}
else{ print OUTPUT "\t$plab"; }
}
elsif($line=~m/^PLAB. =/){
($dummy,$plab)=split/\s*=\s*/,$line;
($lower,$upper)=split/,/,$plab;
if($upper == $lower) {
$lower=1000*$lower;
printf OUTPUT "\t%0.1f",$lower;
}
else{ print OUTPUT "\t$plab"; }
}
elsif($line=~m/^N = YN/) {
@yname=split/\s*=\s*/,$line;
$yname=$yname[$#yname];
}
elsif($line=~m/^\*F/) {
# $dataset++;
($dummy,$xname) = split/\s*=\s*/,$last;
print OUTPUT "\t$xname v $yname";
$count=0;
open DATASET,">$dir/$shortname$dataset.dat";
print DATASET "file $name dataset $dataset\n";
print DATASET "$reac\n";
print DATASET "$obs\n";
print DATASET "$lower MeV/c\n";
print DATASET "$sc $ref $type $date\n";
print DATASET "$yname as a function of $xname\n";
until($line eq "*CF *"){
$line=;
chomp $line;
if($line ne "*CF *" ) { print DATASET "$line\n"; }
$count++;
}
$count--;
print OUTPUT "\t($count)";
print OUTPUT "\n";
}
# if($line=~m/^N =/){ print "$line\n";}
$last=$line;
}
}
}
close OUTPUT;
print HTML "
\n";
print HTML "\n";
print HTML "\n";
close HTML;
}
}
exit()