msf_helper.py

Go to the documentation of this file.

#!/usr/bin/env python
#------------------------------------------------------------#
# warning! do not make *too many* in one go, or FTS will die
# usage:  python msf_helper.py
# prod5, execute as novapro from /nova/app/home/novapro/PROD5_SHIFTERS/mcgen/fcls 
# script to make make_sim_fcl less error-prone
# help:   lasquith@fnal.gov, scalvez@fnal.gov
#------------------------------------------------------------#
from __future__ import print_function
import sys, shlex, subprocess, math, os
from subprocess import Popen, PIPE
debug = True
#** ND : [fhc,rhc][nonswap] make 65k, use 15k
#** FD : [fhc,rhc][nonswap,fluxswap] make 12k, use 6k
# systs   = ["cherenkov","llup","lldown","detectoraging"]
systs   = []#["mcnp"]#["detectoraging"]
dets    = ["nd"]#,"fd"]
swaps   = ["nonswap"]#, "swap"]
horns   = ["0hc"]
pot = "2.5e16"
potper = "spill"
events = "1000"
lastrun = "1001999"
numsubruns = "64"
gainmode = "100"
generator = "genie"
genietune = "N1810j0211a"

iteration = "1"
batchnumber = "1"
nppfx = "0"
timespec = "rocksingles"
nogenierw = True
rockonly = True
singlemod = True
geantonly = True
fts = True
tags = "--nogenierw --fts"
o ="/nova/app/home/novapro/PROD5_SHIFTERS/mcgen/fcls/"

NOVAPRODUCTION_DIR = os.getenv('NOVAPRODUCTION_DIR')

goodruns_fd_fhc = NOVAPRODUCTION_DIR + "/goodruns/fd_grl_v0_period12359.txt"
goodruns_nd_fhc = NOVAPRODUCTION_DIR + "/goodruns/nd_grl_v0_period12359.txt"
goodruns_fd_rhc = NOVAPRODUCTION_DIR + "/goodruns/fd_grl_v0_period4_6-8.txt"
goodruns_nd_rhc = NOVAPRODUCTION_DIR + "/goodruns/nd_grl_v0_period4_6-8.txt"

mix_fd_fhc = "cosmic lasquith_fd_cosmics_p12359_goodruns_prod5"
mix_fd_rhc = "cosmic scalvez_artdaq_fd_cosmics_p4678_goodruns_prod5"
mix_nd_fhc = "rock prod5_nd_fhc_rock_g4"
mix_nd_rhc = "rock prod5_nd_rhc_rock_g4"

plist_fhc="12359"
plist_rhc="4678"

num_fd = "12000"
num_nd = "10000"
            
mk_outdir_cmd="null"

for syst in systs:
    for det in dets:
        for horn in horns:
            if det=="fd":
                swaps   = ["nonswap", "swap"]
                if horn=="fhc":
                        mix = mix_fd_fhc
                        num = num_fd
                        goodruns = goodruns_fd_fhc
                        plist = plist_fhc
                elif horn=="rhc":
                        mix = mix_fd_rhc
                        num = num_fd
                        goodruns = goodruns_fd_rhc
                        plist = plist_rhc
            elif det=="nd":
                swaps   = ["nonswap"]
                if horn=="fhc":
                        mix = mix_nd_fhc
                        num = num_nd
                        goodruns = goodruns_nd_fhc
                        plist = plist_fhc
                elif horn=="rhc":
                        mix = mix_nd_rhc
                        num = num_nd
                        goodruns = goodruns_nd_rhc
                        plist = plist_rhc

            for swap in swaps:
                log = "new_"+syst+"_"+det+"_"+horn+"_"+swap+".log"
                if not debug: log_file = open(log, 'w')
                outdir = o+"new_"+syst+"_"+det+"_"+horn+"_"+swap
                cmd = "make_sim_fcl --syst "+syst+" -det "+det+" -hp "+horn+" -fs "+swap+" -ts periods"+plist+" --mix "+mix+" -rsr "+num+" "+goodruns+" -gt "+genietune+" -bn "+batchnumber+" -i "+iteration+" -gen "+generator+" --nppfx "+nppfx+" "+tags+" --fts -o "+outdir
                if det is "nd":
                    cmd += " --potper"+potper+ " --pot " + pot
                else:
                    cmd += " --events " + events
                mk_outdir_cmd = "mkdir "+outdir

                print("mk_outdir_cmd: ",mk_outdir_cmd)
                print("cmd: ",cmd)
                if not debug:
                    if not os.path.isdir(outdir):
                        p1=Popen( shlex.split(mk_outdir_cmd),stdout=PIPE, stderr=PIPE)
                    p2 = Popen( shlex.split(cmd) ,stdout=log_file ,stderr=subprocess.STDOUT )
           # swaps
        # horns
    # detectors
#systematics
print("Fin.")