code: def _identify_heterogeneity_blocks_shared(in_file, segment_fn, params, work_dir, somatic_info):
      
      out_file = os.path.join(work_dir, "%s-hetblocks.bed" % utils.splitext_plus(os.path.basename(in_file))[0])
      if not utils.file_uptodate(out_file, in_file):
          with file_transaction(somatic_info.tumor_data, out_file) as tx_out_file:
              with open(tx_out_file, "w") as out_handle:
                  for chrom, freqs, coords in _freqs_by_chromosome(in_file, params, somatic_info):
                      for start, end in segment_fn(chrom, freqs, coords):
                          out_handle.write("%s\t%s\t%s\n" % (chrom, start, end))
      return out_file
docs: Identify heterogeneity blocks corresponding to segmentation from CNV input file.