[Bio] / Sprout / SproutLoad.pm Repository:
ViewVC logotype

Diff of /Sprout/SproutLoad.pm

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.13, Fri Sep 16 02:15:02 2005 UTC revision 1.14, Fri Sep 16 02:38:49 2005 UTC
# Line 629  Line 629 
629      # Get the genome hash. We'll use it to filter the genomes in each      # Get the genome hash. We'll use it to filter the genomes in each
630      # spreadsheet.      # spreadsheet.
631      my $genomeHash = $self->{genomes};      my $genomeHash = $self->{genomes};
632        my @genomes = sort keys %{$genomeHash};
633      # Get the subsystem hash. This lists the subsystems we'll process.      # Get the subsystem hash. This lists the subsystems we'll process.
634      my $subsysHash = $self->{subsystems};      my $subsysHash = $self->{subsystems};
635      my @subsysIDs = sort keys %{$subsysHash};      my @subsysIDs = sort keys %{$subsysHash};
# Line 669  Line 670 
670              }              }
671          }          }
672          # Now all roles for this subsystem have been filled in. We create the          # Now all roles for this subsystem have been filled in. We create the
673          # spreadsheet by matches roles to genomes. To do this, we need to          # spreadsheet by matching roles to genomes.
         # get the genomes on the sheet.  
674          Trace("Creating subsystem $subsysID spreadsheet.") if T(3);          Trace("Creating subsystem $subsysID spreadsheet.") if T(3);
675          my @genomes = map { $_->[0] } @{$fig->subsystem_genomes($subsysID)};          for my $genomeID (@genomes) {
         for my $genomeID (sort @genomes) {  
             # Only process this genome if it's one of ours.  
             if (exists $genomeHash->{$genomeID}) {  
                 Trace("Processing genome $genomeID for subsystem $subsysID.");  
                 # Connect the genome to the subsystem.  
                 $loadParticipatesIn->Put($genomeID, $subsysID);  
676                  # Count the PEGs and cells found for verification purposes.                  # Count the PEGs and cells found for verification purposes.
677                  my $pegCount = 0;                  my $pegCount = 0;
678                  my $cellCount = 0;                  my $cellCount = 0;
# Line 704  Line 698 
698                          }                          }
699                      }                      }
700                  }                  }
701                  Trace("$pegCount PEGs found in $cellCount cells.") if T(3);              # If we found some cells for this genome, denote it participates in the
702                # subsystem.
703                if ($pegCount > 0) {
704                    Trace("$pegCount PEGs in $cellCount cells for $genomeID.") if T(3);
705                    $loadParticipatesIn->Put($genomeID, $subsysID);
706              }              }
707          }          }
708      }      }

Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3