@transform(copyCapseqReplicatedBedFiles, suffix(".bed"), ".h3k4me1.bed" )
def getH3K4Me1Overlap( infile, outfile ):
    to_cluster = True
    annotation_file = PARAMS["bed_h3k4me1"]
    if len(annotation_file) > 0:
        statement = """intersectBed -a %(infile)s -b %(annotation_file)s -u > %(outfile)s"""
        P.run()

############################################################
Update the terminal regions (=leaves) of the given tree and
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/getH3K4Me1Overlap--crossover with--dataset/ETHPy150Open scikit-learn/scikit-learn/sklearn/ensemble/gradient_boosting.py/LossFunction.update_terminal_regions