コード例 #1
0
ファイル: cvs2local.py プロジェクト: khotilov/cmssw
     #print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass
     continue
 if channel == "vhbb" or channel == "hmm" or channel == "hbb":
     for category in categories[channel]:
         if options.verbose:
             print "copying datacards for:", period, channel, category, mass
         if options.analysis == "mssm":
             os.system(
                 "cp {INPUT}/{CHN}/{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root {OUTPUT}/common/{PRE}{CHN}.input_{PERIOD}-{MASSCAT}.root"
                 .format(INPUT=input,
                         ANA=options.analysis,
                         CHN=channel,
                         OUTPUT=options.out,
                         PRE=prefix,
                         PERIOD=period,
                         MASSCAT=mass_category(
                             mass, category, channel)))
             os.system(
                 "cp {INPUT}/{CHN}/{CHN}_{CAT}_{PERIOD}-{MASS}.txt {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt"
                 .format(INPUT=input,
                         CHN=channel,
                         CAT=category,
                         PERIOD=period,
                         MASS=mass,
                         MASSDIR=massdir(mass),
                         OUTPUT=options.out,
                         PRE=prefix))
             os.system(
                 "perl -pi -e 's/{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root/..\/common\/{PRE}{CHN}.input_{PERIOD}-{MASSCAT}.root/g' {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt"
                 .format(CHN=channel,
                         ANA=options.analysis,
                         PRE=prefix,
コード例 #2
0
            ## check validity of run period
            if not period in valid_periods[channel] :
                #print "drop due to failing period: ", channel, valid_periods[channel], period
                continue
            for category in categories[channel] :
                if options.verbose :
                    print "copying datacards for:", period, channel, category, mass
                ## check validity of run period
                if not period in valid_periods[channel] :
                    #print "drop due to failing period: ", channel, valid_periods[channel], period
                    continue
                if options.analysis == "mssm" :
                    add_mass("htt_{CHN}_{CAT}_{PERIOD}".format(CHN=channel, CAT=category, PERIOD=period), mass)
                    if options.no_update:
                        files = ' '.join(os.listdir("{OUTPUT}/common/".format(OUTPUT=options.out)))
                        if not "htt_"+channel+'.inputs-mssm-'+period in files :
                            os.system("cp {INPUT}/htt_{CHN}/htt_{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root* {OUTPUT}/common/".format(
                                INPUT=input, CHN=channel, ANA=options.analysis, PERIOD=period, OUTPUT=options.out, MASSCAT=mass_category(mass,category, channel)))
                        else :
                            pass
                        #print "no update of files needed."
                    else :
                        os.system("cp {INPUT}/htt_{CHN}/htt_{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root* {OUTPUT}/common/".format(
                            INPUT=input, CHN=channel, ANA=options.analysis, PERIOD=period, OUTPUT=options.out, MASSCAT=mass_category(mass,category, channel)))
                    os.system("cp {INPUT}/htt_{CHN}/htt_{CHN}_{CAT}_{PERIOD}-100.txt {OUTPUT}/{MASSDIR}/htt_{CHN}_{CAT}_{PERIOD}.txt".format( #mA=110 always for lowmH!
                        INPUT=input, CHN=channel, CAT=category, PERIOD=period, MASSDIR=massdir(mass), OUTPUT=options.out))
                    os.system("perl -pi -e 's/htt_{CHN}.inputs/..\/common\/htt_{CHN}.inputs/g' {OUTPUT}/{MASSDIR}/htt_{CHN}_{CAT}_{PERIOD}.txt".format(
                            CHN=channel, OUTPUT=options.out, MASS=mass, MASSDIR=massdir(mass), CAT=category, PERIOD=period))
                else :
                    print "options.analysis has to be 'mssm' for lowmH"
コード例 #3
0
for channel in channels:
    for period in periods:
        for cat in ['0']:
            for mass in parseArgs(args):
                for process in ['ggH', 'bbH']:
                    exe = "{CMSSW_BASE}/src/HiggsAnalysis/HiggsToTauTau/macros/rescaleSignal.C+".format(
                        CMSSW_BASE=os.environ.get("CMSSW_BASE"))
                    os.system(
                        r"root -l -b -q {EXE}\(true,{SCALE},\"{PATH}/{CHN}/htt_{CHN}.inputs-mssm-{PER}-{MASSCAT}.root\",\"{PROCESS}\",0\)"
                        .format(EXE=exe,
                                PATH=source_path,
                                SCALE=acceptance_correction(
                                    process, mass, period),
                                CHN=channel,
                                PER=period,
                                MASSCAT=mass_category(mass, cat,
                                                      'htt_' + channel),
                                PROCESS=process + str(mass)))
                    if channel == "mm":
                        os.system(
                            r"root -l -b -q {EXE}\(true,{SCALE},\"{PATH}/{CHN}/htt_{CHN}.inputs-mssm-{PER}-{MASSCAT}-msv.root\",\"{PROCESS}\",0\)"
                            .format(EXE=exe,
                                    PATH=source_path,
                                    SCALE=acceptance_correction(
                                        process, mass, period),
                                    CHN=channel,
                                    PER=period,
                                    MASSCAT=mass_category(
                                        mass, cat, 'htt_' + channel),
                                    PROCESS=process + str(mass)))
コード例 #4
0
     if options.analysis == "sm" and "mm" in channel :
         ## this has been removed since there is no problem anymore with morphed templates in SM mm / ee
         ##mass = closest_simulated_masspoint(mass)
         ##fudge_mm_datacards = float(fudge_mass)-float(mass)!=0
         pass
 ## check validity of mass
 if (float(mass)< valid_masses[channel][0] or float(mass)> valid_masses[channel][1]) :
     #print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass
     continue
 print "creating datacard for:", options.analysis, period, channel, cat, fudge_mass
 if options.analysis == "mssm" :
     os.system("create-datacard.py -i {CHN}.inputs-{ANA}-{PER}-{MASSCAT}.root -o {CHN}_{CAT}_{PER}-{LABEL}.txt {MASS}".format(
         CHN=prefix+channel,
         ANA=options.analysis,
         PER=period,
         MASSCAT=mass_category(mass,cat,prefix+channel),
         CAT=cat,
         MASS='' if options.ignore_mass_argument else mass,
         LABEL='125' if options.ignore_mass_argument else mass
         ))
 if options.analysis == "sm" :
     if options.SM4 :
         print "rescaling signal cross sections accoring to SM4 cross sections"
         os.system(r"root -q -l -b {CMSSW_BASE}/src/HiggsAnalysis/HiggsToTauTau/macros/rescale2SM4.C+\(true,\"{CHN}.inputs-{per}.root\"\)".format(
             CMSSW_BASE=os.environ['CMSSW_BASE'],
             CHN=channel,
             per=period.lower()
             ))
     os.system("create-datacard.py -i {CHN}.inputs-{ANA}-{per}.root -o {CHN}_{CAT}_{PER}-{LABEL}.txt {MASS} {SM_HIGGS_BKG}".format(
         CHN=prefix+channel,
         ANA=options.analysis,
コード例 #5
0
     # print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass
     continue
 if channel == "vhbb" or channel == "hmm" or channel == "hbb":
     for category in categories[channel]:
         if options.verbose:
             print "copying datacards for:", period, channel, category, mass
         if options.analysis == "mssm":
             os.system(
                 "cp {INPUT}/{CHN}/{PRE}{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root* {OUTPUT}/common".format(
                     INPUT=input,
                     ANA=options.analysis,
                     CHN=channel,
                     OUTPUT=options.out,
                     PRE=prefix,
                     PERIOD=period,
                     MASSCAT=mass_category(mass, category, channel),
                 )
             )
             os.system(
                 "cp {INPUT}/{CHN}/{CHN}_{CAT}_{PERIOD}-{MASS}.txt {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     INPUT=input,
                     CHN=channel,
                     CAT=category,
                     PERIOD=period,
                     MASS=mass,
                     MASSDIR=massdir(mass),
                     OUTPUT=options.out,
                     PRE=prefix,
                 )
             )
             os.system(
コード例 #6
0
channels = options.channels.split()
for idx in range(len(channels)) : channels[idx] = channels[idx].rstrip(',')

for channel in channels :
    for period in periods :
        print period
        for cat in ['0'] :
            for mass in parseArgs(args) :
                for process in ['ggH', 'bbH'] :
                    exe = "{CMSSW_BASE}/src/HiggsAnalysis/HiggsToTauTau/macros/rescaleSignal.C+".format(CMSSW_BASE=os.environ.get("CMSSW_BASE"))
                    os.system(r"root -l -b -q {EXE}\(true,{SCALE},\"{PATH}/{CHN}/htt_{CHN}.inputs-mssm-{PER}-{MASSCAT}.root\",\"{PROCESS}\",0\)".format(
                        EXE=exe,
                        PATH=source_path,
                        SCALE=acceptance_correction(process, mass, period),
                        CHN=channel,
                        PER=period,
                        MASSCAT=mass_category(mass,cat,'htt_'+channel),
                        PROCESS=process+str(mass)
                        ))
                    if channel=="mm" :
                        os.system(r"root -l -b -q {EXE}\(true,{SCALE},\"{PATH}/{CHN}/htt_{CHN}.inputs-mssm-{PER}-{MASSCAT}-msv.root\",\"{PROCESS}\",0\)".format(
                            EXE=exe,
                            PATH=source_path,
                            SCALE=acceptance_correction(process, mass, period),
                            CHN=channel,
                            PER=period,
                            MASSCAT=mass_category(mass,cat,'htt_'+channel),
                            PROCESS=process+str(mass)
                            ))
                    
コード例 #7
0
         ##fudge_mm_datacards = float(fudge_mass)-float(mass)!=0
         pass
 ## check validity of mass
 if (float(mass) < valid_masses[channel][0]
         or float(mass) > valid_masses[channel][1]):
     #print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass
     continue
 print "creating datacard for:", options.analysis, period, channel, cat, fudge_mass
 if options.analysis == "mssm":
     os.system(
         "create-datacard.py -i {CHN}.inputs-{ANA}-{PER}-{MASSCAT}.root -o {CHN}_{CAT}_{PER}-{LABEL}.txt {MASS}"
         .format(
             CHN=prefix + channel,
             ANA=options.analysis,
             PER=period,
             MASSCAT=mass_category(mass, cat, prefix + channel),
             CAT=cat,
             MASS='' if options.ignore_mass_argument else mass,
             LABEL='125'
             if options.ignore_mass_argument else mass))
 if options.analysis == "sm":
     if options.SM4:
         print "rescaling signal cross sections accoring to SM4 cross sections"
         os.system(
             r"root -q -l -b {CMSSW_BASE}/src/HiggsAnalysis/HiggsToTauTau/macros/rescale2SM4.C+\(true,\"{CHN}.inputs-{per}.root\"\)"
             .format(CMSSW_BASE=os.environ['CMSSW_BASE'],
                     CHN=channel,
                     per=period.lower()))
     os.system(
         "create-datacard.py -i {CHN}.inputs-{ANA}-{per}.root -o {CHN}_{CAT}_{PER}-{LABEL}.txt {MASS} {SM_HIGGS_BKG}"
         .format(
コード例 #8
0
 for mass in parseArgs(args) :
     ## check validity of run period
     if not period in valid_periods[channel] :
         #print "drop due to failing period: ", channel, valid_periods[channel], period
         continue
     ## check validity of mass
     if (float(mass)< valid_masses[channel][0] or float(mass)> valid_masses[channel][1]) :
         #print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass
         continue
     if channel == "vhbb" or channel == "hmm" or channel == "hbb" :
         for category in categories[channel] :
             if options.verbose :
                 print "copying datacards for:", period, channel, category, mass
             if options.analysis == "mssm" :
                 os.system("cp {INPUT}/{CHN}/{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root {OUTPUT}/common/{PRE}{CHN}.input_{PERIOD}-{MASSCAT}.root".format(
                     INPUT=input, ANA=options.analysis, CHN=channel, OUTPUT=options.out, PRE=prefix, PERIOD=period, MASSCAT=mass_category(mass,category, channel)))
                 os.system("cp {INPUT}/{CHN}/{CHN}_{CAT}_{PERIOD}-{MASS}.txt {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     INPUT=input, CHN=channel, CAT=category, PERIOD=period, MASS=mass, MASSDIR=massdir(mass), OUTPUT=options.out, PRE=prefix))
                 os.system("perl -pi -e 's/{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root/..\/common\/{PRE}{CHN}.input_{PERIOD}-{MASSCAT}.root/g' {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     CHN=channel, ANA=options.analysis, PRE=prefix, OUTPUT=options.out, MASS=mass, MASSDIR=massdir(mass), CAT=category, PERIOD=period, MASSCAT=mass_category(mass,category,channel)))
             else :
                 os.system("cp {INPUT}/{CHN}/{CHN}.inputs-{ANA}-{PERIOD}.root {OUTPUT}/common/{PRE}{CHN}.input_{PERIOD}.root".format(
                     INPUT=input, ANA=options.analysis, CHN=channel, OUTPUT=options.out, PRE=prefix, PERIOD=period))
                 os.system("cp {INPUT}/{CHN}/{CHN}_{CAT}_{PERIOD}-{MASS}.txt {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     INPUT=input, CHN=channel, CAT=category, PERIOD=period, MASS=mass, MASSDIR=massdir(mass), OUTPUT=options.out, PRE=prefix))
                 os.system("perl -pi -e 's/{CHN}.inputs-{ANA}-{PERIOD}.root/..\/common\/{PRE}{CHN}.input_{PERIOD}.root/g' {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     CHN=channel, ANA=options.analysis, PRE=prefix, OUTPUT=options.out, MASS=mass, MASSDIR=massdir(mass), CAT=category, PERIOD=period))
             if options.analysis == "mssm" :
                 add_mass("{CHN}_{CAT}_{PERIOD}".format(CHN=channel, CAT=category, PERIOD=period), mass)
     elif channel == "vhtt":
         for category in categories[channel] :
コード例 #9
0
        for cat in categories[channel] :
            ## here the normal workflow continues
            prefix = "" if (channel == "vhtt" or channel == "vhbb") else "htt_"
            os.chdir("{PWD}/{CHN}".format(CHN=prefix+channel, PWD=base))
            ## check validity of run period
            if not period in valid_periods[channel] :
                #print "drop due to failing period: ", channel, valid_periods[channel], period 
                continue
            os.system("datacard-project.py -i {PATH} -c {CHN} -e {ANA}-{PER}-0{CAT} {PER}-0{CAT}".format(PATH=os.path.abspath(options.input), CHN=channel, ANA=options.analysis, PER=period, CAT=cat))
            os.chdir("{PWD}/{CHN}/{PER}-0{CAT}".format(CHN=prefix+channel, PER=period, PWD=base, CAT=cat))
            for mass in parseArgs(args) :
                ## check validity of mass
                if (float(mass)< valid_masses[channel][0] or float(mass)> valid_masses[channel][1]) :
                    #print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass 
                    continue
                print "creating datacard for:", options.analysis, period, channel, cat, mass
                if options.analysis == "mssm" :
                    os.system("create-datacard.py -i {CHN}.inputs-{ana}-{PER}-{MASSCAT}.root -o {CHN}_{CAT}_{PER}-{MASS}.txt {MASS}".format(
                        CHN=prefix+channel, ana=options.analysis, PER=period, MASSCAT=mass_category(mass), CAT=cat, MASS=mass))
                if options.analysis == "sm" :
                    if options.SM4 :
                        print "rescaling signal cross sections accoring to SM4 cross sections"
                        os.system(r"root -q -l -b {CMSSW_BASE}/src/HiggsAnalysis/HiggsToTauTau/macros/rescale2SM4.C+\(true,\"{CHN}.inputs-{per}.root\"\)".format(
                            CMSSW_BASE=os.environ['CMSSW_BASE'], CHN=channel, per=period.lowe()))
                    os.system("create-datacard.py -i {CHN}.inputs-{ANA}-{per}.root -o {CHN}_{CAT}_{PER}-{MASS}.txt {MASS}".format(
                        CHN=prefix+channel, ANA=options.analysis, per=period, CAT=cat, PER=period, MASS=mass))
            os.system("mv *.* ../")
            os.chdir("{PWD}/{CHN}".format(CHN=prefix+channel, PWD=base))
            os.system("rm -r {PER}-0{CAT}".format(PER=period, CAT=cat))
            os.system("rm -r cgs.* unc.*")
コード例 #10
0
 for mass in parseArgs(args) :
     ## check validity of run period
     if not period in valid_periods[channel] :
         #print "drop due to failing period: ", channel, valid_periods[channel], period
         continue
     ## check validity of mass
     if (float(mass)< valid_masses[channel][0] or float(mass)> valid_masses[channel][1]) :
         #print "drop due to failing mass:" , channel, valid_masses[channel][0], valid_masses[channel][1], ":", mass
         continue
     if channel == "vhbb" or channel == "hmm" or channel == "hbb" :
         for category in categories[channel] :
             if options.verbose :
                 print "copying datacards for:", period, channel, category, mass
             if options.analysis == "mssm" :
                 os.system("cp {INPUT}/{CHN}/{PRE}{CHN}.inputs-{ANA}-{PERIOD}-{MASSCAT}.root* {OUTPUT}/common".format(
                     INPUT=input, ANA=options.analysis, CHN=channel, OUTPUT=options.out, PRE=prefix, PERIOD=period, MASSCAT=mass_category(mass,category, channel)))
                 os.system("cp {INPUT}/{CHN}/{CHN}_{CAT}_{PERIOD}-{MASS}.txt {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     INPUT=input, CHN=channel, CAT=category, PERIOD=period, MASS='100' if options.model=='lowmH' or options.model=='2HDM' else mass, MASSDIR=massdir(mass), OUTPUT=options.out, PRE=prefix))
                 os.system("perl -pi -e 's/{CHN}.inputs/..\/common\/{PRE}{CHN}.inputs/g' {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     CHN=channel, ANA=options.analysis, PRE=prefix, OUTPUT=options.out, MASS=mass, MASSDIR=massdir(mass), CAT=category, PERIOD=period, MASSCAT=mass_category(mass,category,channel)))
             else :
                 os.system("cp {INPUT}/{CHN}/{CHN}.inputs-{ANA}-{PERIOD}.root {OUTPUT}/common/{PRE}{CHN}.input_{PERIOD}.root".format(
                     INPUT=input, ANA=options.analysis, CHN=channel, OUTPUT=options.out, PRE=prefix, PERIOD=period))
                 os.system("cp {INPUT}/{CHN}/{CHN}_{CAT}_{PERIOD}-{MASS}.txt {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     INPUT=input, CHN=channel, CAT=category, PERIOD=period, MASS='300' if options.model=='lowmH' or options.model=='2HDM' else mass, MASSDIR=massdir(mass), OUTPUT=options.out, PRE=prefix))
                 os.system("perl -pi -e 's/{CHN}.inputs-{ANA}-{PERIOD}.root/..\/common\/{PRE}{CHN}.input_{PERIOD}.root/g' {OUTPUT}/{MASSDIR}/{PRE}{CHN}_{CAT}_{PERIOD}.txt".format(
                     CHN=channel, ANA=options.analysis, PRE=prefix, OUTPUT=options.out, MASS=mass, MASSDIR=massdir(mass), CAT=category, PERIOD=period))
     elif channel == "vhtt":
         for category in categories[channel] :
             if options.verbose :
                 print "copying datacards for:", period, channel, category, mass
コード例 #11
0
channels = options.channels.split()
for idx in range(len(channels)) : channels[idx] = channels[idx].rstrip(',')

for channel in channels :
    for period in periods :
        for cat in ['0'] :
            for mass in parseArgs(args) :
                for process in ['ggH', 'bbH'] :
                    exe = "{CMSSW_BASE}/src/HiggsAnalysis/HiggsToTauTau/macros/rescaleSignal.C+".format(CMSSW_BASE=os.environ.get("CMSSW_BASE"))
                    os.system(r"root -l -b -q {EXE}\(true,{SCALE},\"{PATH}/{CHN}/htt_{CHN}.inputs-{MSSM}-{PER}{MASSCAT}.root\",\"{PROCESS}\",0\)".format(
                        EXE=exe,
                        PATH=source_path,
                        SCALE=acceptance_correction(process, mass, period),
                        CHN=channel,
                        MSSM="Hhh" if options.Hhh=="Hhh" else "mssm",
                        PER=period,
                        MASSCAT="-" + str(mass_category(mass,cat,'htt_'+channel)) if options.Hhh=="" else "",
                        PROCESS=process+str(mass)
                        ))
                    if channel=="mm" :
                        os.system(r"root -l -b -q {EXE}\(true,{SCALE},\"{PATH}/{CHN}/htt_{CHN}.inputs-mssm-{PER}-{MASSCAT}-msv.root\",\"{PROCESS}\",0\)".format(
                            EXE=exe,
                            PATH=source_path,
                            SCALE=acceptance_correction(process, mass, period),
                            CHN=channel,
                            PER=period,
                            MASSCAT=mass_category(mass,cat,'htt_'+channel),
                            PROCESS=process+str(mass)
                            ))