| 21 |
|
my $nBadJobs = 0; |
| 22 |
|
my $counting = 0; |
| 23 |
|
my %exitCodes; |
| 24 |
+ |
my %signals; |
| 25 |
|
my %crossSections; |
| 26 |
|
my $integratedLuminosity = 10000; |
| 27 |
|
$integratedLuminosity = $opt{"luminosity"} if $opt{"luminosity"}; |
| 42 |
|
close (FILE); |
| 43 |
|
my $fileContents = join ("", @fileContents); |
| 44 |
|
$fileContents =~ s/\n/ /g; |
| 45 |
< |
$fileContents =~ s/.*\(return value ([^)]*)\).*/$1/g; |
| 45 |
< |
$nGoodJobs++; |
| 46 |
< |
if ($fileContents != 0) |
| 45 |
> |
if ($fileContents =~ m/return value/) |
| 46 |
|
{ |
| 47 |
+ |
$fileContents =~ s/.*\(return value ([^)]*)\).*/$1/g; |
| 48 |
+ |
$nGoodJobs++; |
| 49 |
+ |
print "WARNING: Nonzero exit code for job $jobNumber. (return value $fileContents)\n" if $fileContents != 0; |
| 50 |
+ |
$exitCodes{$dir}{$jobNumber} = $fileContents; |
| 51 |
+ |
$counting = 1; |
| 52 |
+ |
} |
| 53 |
+ |
if ($fileContents =~ m/signal/) |
| 54 |
+ |
{ |
| 55 |
+ |
$fileContents =~ s/.*\(signal ([^)]*)\).*/$1/g; |
| 56 |
|
$nBadJobs++; |
| 57 |
< |
print "Skipping job $jobNumber. (return value $fileContents)\n"; |
| 57 |
> |
print "WARNING: Skipping job $jobNumber. (signal $fileContents)\n"; |
| 58 |
> |
$signals{$dir}{$jobNumber} = $fileContents; |
| 59 |
> |
$counting = 1; |
| 60 |
|
} |
| 51 |
– |
$exitCodes{$dir}{$jobNumber} = $fileContents; |
| 52 |
– |
$counting = 1; |
| 61 |
|
} |
| 62 |
|
if ($file =~ m/^.*\/crossSectionInPicobarn\.txt$/) |
| 63 |
|
{ |
| 78 |
|
{ |
| 79 |
|
my $jobNumber = $file; |
| 80 |
|
$jobNumber =~ s/^.*_([^_]*)\.root$/$1/; |
| 81 |
< |
$badJob = defined $exitCodes{$dir} && defined $exitCodes{$dir}{$jobNumber} && $exitCodes{$dir}{$jobNumber}; |
| 81 |
> |
$badJob = defined $signals{$dir} && defined $signals{$dir}{$jobNumber}; |
| 82 |
|
} |
| 83 |
|
next if $badJob; |
| 84 |
|
if ($file =~ m/^.*\.root$/) |
| 98 |
|
my $nTotalEvents = countEvents ("$opt{'prefix'}.root", $cutFlow); |
| 99 |
|
unlink ("$opt{'prefix'}.root"); |
| 100 |
|
@weights = (); |
| 101 |
+ |
my %weights = (); |
| 102 |
|
foreach my $file (@$files) |
| 103 |
|
{ |
| 104 |
|
next if $file eq "."; |
| 110 |
|
{ |
| 111 |
|
my $jobNumber = $file; |
| 112 |
|
$jobNumber =~ s/^.*_([^_]*)\.root$/$1/; |
| 113 |
< |
$badJob = defined $exitCodes{$dir} && defined $exitCodes{$dir}{$jobNumber} && $exitCodes{$dir}{$jobNumber}; |
| 113 |
> |
$badJob = defined $signals{$dir} && defined $signals{$dir}{$jobNumber}; |
| 114 |
|
} |
| 115 |
|
next if $badJob; |
| 116 |
|
if ($file =~ m/^.*\.root$/) |
| 119 |
|
push (@weights, ($opt{"xsection"} * $integratedLuminosity) / $nTotalEvents) if !$opt{"weight"} && $opt{"xsection"}; |
| 120 |
|
push (@weights, ($crossSections{$dir} * $integratedLuminosity) / $nTotalEvents) if !$opt{"weight"} && !$opt{"xsection"} && defined $crossSections{$dir}; |
| 121 |
|
push (@weights, 1.0) if !$opt{"weight"} && !$opt{"xsection"} && !(defined $crossSections{$dir}); |
| 122 |
+ |
foreach my $arg (@ARGV) |
| 123 |
+ |
{ |
| 124 |
+ |
if (substr ($file, 0, length ($arg)) eq $arg) |
| 125 |
+ |
{ |
| 126 |
+ |
$arg =~ s/^.*\/([^\/]*)$/$1/; |
| 127 |
+ |
$weights{$arg} = $weights[-1] |
| 128 |
+ |
} |
| 129 |
+ |
} |
| 130 |
|
} |
| 131 |
|
} |
| 132 |
|
my $weights = join (",", @weights); |
| 133 |
|
system ("mergeTFileServiceHistograms -i $rootFiles -o $opt{'prefix'}.root -w $weights"); |
| 134 |
< |
system ("cutFlowTable $opt{'prefix'}.root $cutFlow >& $opt{'prefix'}.tex"); |
| 134 |
> |
print "Weights:\n"; |
| 135 |
> |
foreach my $arg (keys %weights) |
| 136 |
> |
{ |
| 137 |
> |
printf " $arg: %.1f (%.1f pb)\n", $weights{$arg}, $weights{$arg} / $integratedLuminosity; |
| 138 |
> |
} |
| 139 |
|
my $goodEvents = countEvents ("$opt{'prefix'}.root", $cutFlow); |
| 140 |
|
print "$nGoodJobs jobs ran successfully over $nTotalEvents ($goodEvents weighted) events.\n" if $counting; |
| 141 |
|
print "$nBadJobs jobs failed to run.\n" if $counting; |
| 183 |
|
$exeName =~ s/^.*\/([^\/]*)$/$1/; |
| 184 |
|
|
| 185 |
|
print "Usage: $exeName [OPTION]... -p PREFIX DIRECTORIES_AND_FILES\n"; |
| 186 |
< |
print "Merges ROOT files containing histograms and produces a cutflow table in the\n"; |
| 187 |
< |
print "form of a LaTeX document. If there are Condor logs in the specified\n"; |
| 188 |
< |
print "directories, checks for nonzero return values. If the directories where created\n"; |
| 189 |
< |
print "by \"osusub\", uses the cross section from the database to weight all histograms.\n"; |
| 186 |
> |
print "Merges ROOT files containing histograms. If there are Condor logs in the\n"; |
| 187 |
> |
print "specified directories, checks for nonzero return values. If the directories\n"; |
| 188 |
> |
print "where created by \"osusub\", uses the cross section from the database to weight\n"; |
| 189 |
> |
print "all histograms.\n"; |
| 190 |
|
print "\n"; |
| 191 |
|
print "Mandatory arguments to long options are mandatory for short options too.\n"; |
| 192 |
|
printf "%-29s%s\n", " -c, --cutflow HISTOGRAM", "name of histogram to use for the cutflow (default:"; |
| 194 |
|
printf "%-29s%s\n", " -h, --help", "print this help message"; |
| 195 |
|
printf "%-29s%s\n", " -l, --luminosity", "integrated luminosity to which the histograms are"; |
| 196 |
|
printf "%-29s%s\n", " ", "weighted (default: 10000/pb)"; |
| 197 |
< |
printf "%-29s%s\n", " -p, --prefix PREFIX", "output files are named PREFIX.root and PREFIX.tex"; |
| 197 |
> |
printf "%-29s%s\n", " -p, --prefix PREFIX", "output is named PREFIX.root"; |
| 198 |
|
printf "%-29s%s\n", " -w, --weight WEIGHT", "scale the output by WEIGHT, overriding the"; |
| 199 |
|
printf "%-29s%s\n", " ", "automatic weighting using the cross section from"; |
| 200 |
|
printf "%-29s%s\n", " ", "the database"; |
