run_prod Namespace Reference

Functions
None	makeProd (prefix, DY, bool y=False, bool phiRandom=False, X=None)
None	copy2EOS ()
	makeProd("muon812",10,False,True) # –< 831 copied back, done 16.3.2015 makeProd("muon822",10,True,True) makeProd("muon821",10,True,True) makeProd("muon822",10,True,True)
None	copyFromEOS ()

Function Documentation

copy2EOS()

None run_prod.copy2EOS

(

)

makeProd("muon812",10,False,True) # –< 831 copied back, done 16.3.2015 makeProd("muon822",10,True,True) makeProd("muon821",10,True,True) makeProd("muon822",10,True,True)

makeProd("muon813",10,False,True) # makeProd("muon823",10,True,True)

Definition at line 228 of file run_prod.py.

def copy2EOS() -> None:
    import os
 
    eos = "/afs/cern.ch/project/eos/installation/0.3.15/bin/eos.select"
    for prod in [610, 620]:
        for run in range(0, 10):
            prefix = "muon" + str(prod + run)
            if prod in [610, 620] and run == 0:
                prefix = "muon" + str(int(prod / 100))
            for i in range(1, 10):
                # requires full path
                cmd = (
                    eos
                    + " cp -r "
                    + os.path.abspath(".")
                    + "/"
                    + prefix
                    + str(i)
                    + "/ /eos/experiment/ship/data/muonBackground/"
                    + prefix
                    + str(i)
                    + "/"
                )
                print(cmd)
                os.system(cmd)
 
 

copyFromEOS()

None run_prod.copyFromEOS

(

)

Definition at line 255 of file run_prod.py.

def copyFromEOS() -> None:
    import os
 
    eos = "/afs/cern.ch/project/eos/installation/0.3.15/bin/eos.select"
    for prod in [610, 620]:
        for run in range(0, 10):
            prefix = "muon" + str(prod + run)
            if prod in [610, 620] and run == 0:
                prefix = "muon" + str(int(prod / 100))
            for i in range(1, 10):
                # requires full path
                cmd = (
                    eos
                    + " cp -r  /eos/experiment/ship/data/muonBackground/"
                    + prefix
                    + str(i)
                    + "/ "
                    + os.path.abspath(".")
                    + "/"
                    + prefix
                    + str(i)
                    + "/"
                )
                print(cmd)
                os.system(cmd)

makeProd()

None run_prod.makeProd	(		prefix,
			DY,
		bool	y = False,
		bool	phiRandom = False,
			X = None
	)

Definition at line 10 of file run_prod.py.

def makeProd(prefix, DY, y: bool = False, phiRandom: bool = False, X=None) -> None:
    ncpu = 9
    shipsoft = os.environ["SHIPSOFT"].replace("/dev", "")
    if not y:
        f = shipsoft + "/data/pythia8_Geant4_onlyMuons.root"
        cmd = "python $FAIRSHIP/macro/run_simScript.py --MuonBack -f " + f + " -Y " + str(float(DY))  # --display"
    elif y == "Jpsi":
        f = shipsoft + "/data/pythia8_Geant4_Jpsi_onlyMuons.root"
        cmd = "python $FAIRSHIP/macro/run_simScript.py --MuonBack -f " + f + " -Y " + str(float(DY))  # --display"
    else:
        f = shipsoft + "/data/pythia8_Geant4_Yandex_onlyMuons.root"
        cmd = "python $FAIRSHIP/macro/run_simScript.py --MuonBack -f " + f + " -Y " + str(float(DY))  # --display"
    if phiRandom:
        cmd = cmd + " --phiRandom"
    fn = ROOT.TFile(f)
    sTree = fn.FindObjectAny("pythia8-Geant4")
    ntot = sTree.GetEntries()
    fn.Close()
    ns = 0
    n3 = int(ntot / ncpu)
    for i in range(1, ncpu + 1):
        d = prefix + str(i)
        if d not in os.listdir("."):
            os.system("mkdir " + d)
    os.chdir("./" + prefix + "1")
    for i in range(1, ncpu + 1):
        if i == ncpu:
            n3 = ntot - (i - 1) * n3
        if X:
            if i == X:
                os.system("cp $FAIRSHIP/macro/run_simScript.py .")
                os.system(cmd + " -n " + str(n3) + " -i " + str(ns) + " > log &")
        else:
            os.system("cp $FAIRSHIP/macro/run_simScript.py .")
            os.system(cmd + " -n " + str(n3) + " -i " + str(ns) + " > log &")
            time.sleep(5)
        # print " -n "+str(n3)+" -i "+str(ns)
        ns += n3
        if i == ncpu:
            os.chdir("../")
            break
        os.chdir("../" + prefix + str(i + 1))
 
 
#
# 11-19 with QGSP_BERT_EMV instead of QGSP_BERT_HP_PEN
# 51-59 passive shielding
# 41-49 fixed field horizontal
# 61-69  chamber with vacuum
# 71-79 without field in the wings
# 81-89 add shielding
# 91-99 also 111-119 switch off  muBrems, muPair
# 101-109 vacuum tube as cone, Al->Steel
# 121-129 switch off only muIoni
# 151-159 vacuum tube as cone, Al->Steel, air in cave
# 161-169 new geometry
# 251-259 passive shielding with concrete walls
# 171-179 new increased Bdl,  2.5m clearance from start, 9m space for tau
# 181-189 even more field, narrow tunnel around first magnet, add 2m hadron absorber
# 311-314 change numbering of veto volumes
# 321-324 1 Teslas
# 331-339 added TauSensitive, fixed C1, using muons from Yandex
# 341-349 should be the same as 331-339
# 351-354 fixed RpcPoints, fixed field !  also - 359 made on t3
# 361-364 4m high magnets!
# 371-374 new strawpoints  with magC1 fixed
# 381-389 yandex file, overwrite with magC1 fixed
# 391-399 change to dY = 2m
# 401-409 change to dY = 2m Yandex
# 411-419 back to dY = 4m, add sensitive plance downstream det2
# 421-429 reduce size of vacuum tank before veto station, add pdgcode for muonPoint, move det2 closer
# 431-439 reduce size of vacuum tank before veto station, add pdgcode for muonPoint, move det2 closer Yandex
# 441-449 try with increased tau mu detector, and reduced height of active muon shield, dy=1m
# 451-459 try with increased tau mu detector, and reduced height of active muon shield, dy=1m  Yandex
# 461-469 dy=4m and then 1m, RPC fixed, display on
# 471-479 dy=4m and then 1m, RPC fixed, Yandex
# 481-489 dy=4m and then 1m, RPC fixed, display OFF
# 491-499 back to dy=1m  display on
# 501-509 dy=4m for last magnet part,  display on
# 511-519 dy=2,3,4m for last magnet parts,  display on
# 521-529 removal of almost overlaps
# 551-554 ultimate slim design
# 561-564 fixing overlaps hcal and almost overlaps in magnets
# 571-574 more undetected overlaps fixed
# 581-584 test with replacing tungsten core with iron
# reuse the following, minor mod det2 with shielding, different lengths, new fine tuned active shield, no hcal overlap
# new scaling of L vs Y
# 591-599 6m height
# 601-609 6m height, Yandex data
# 610-619 10m height
# 620-629 10m height, Yandex data
#
# 630-639 10m height, change RPC width 4.5->3.6m
# 640-649  6m height, change RPC width 4.5->3.6m
# makeProd("muon59",6,False)  #
# makeProd("muon60",6,True)    #
# makeProd("muon61",10,False)    #
# makeProd("muon62",10,True)     #
 
# makeProd("muon63",10,False)     #
# makeProd("muon64",10,True)     #
# makeProd("muon65",6,False)     #
# makeProd("muon66",6,True)     #
 
# run more statistics, 61, 62
# makeProd("muon611",10,False,True)
# makeProd("muon621",10,True,True)
# makeProd("muon612",10,False,True)
# makeProd("muon622",10,True,True)
# makeProd("muon613",10,False,True)
# makeProd("muon623",10,True,True)
# makeProd("muon614",10,False,True)
# makeProd("muon624",10,True,True)
# makeProd("muon615",10,False,True)
# makeProd("muon625",10,True,True)
# makeProd("muon616",10,False,True)
# makeProd("muon626",10,True,True)
# makeProd("muon617",10,False,True)
# makeProd("muon627",10,True,True)
# makeProd("muon618",10,False,True)
# makeProd("muon628",10,True,True)
# makeProd("muon619",10,False,True)
# makeProd("muon629",10,True,True)
# makeProd("muon700",10,False,False) # switch off field of active shielding    # prefix,DY,y=False,phiRandom=False,X=None
# makeProd("muon710",10,False,False) # start production with beam smeared on r=3cm disk
# makeProd("muon720",10,True,False)
# makeProd("muon711",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon721",10,True,True)
# makeProd("muon712",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon722",10,True,True)
# makeProd("muon713",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon723",10,True,True)
# makeProd("muon714",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon724",10,True,True)
# makeProd("muon715",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon725",10,True,True)
# makeProd("muon716",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon726",10,True,True)
# makeProd("muon717",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon777",10,False,True) # in case the other one does not work 717
# makeProd("muon727",10,True,True)  # ?
# makeProd("muon718",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon728",10,True,True)
# makeProd("muon719",10,False,True) # start production with beam smeared on r=3cm disk
# makeProd("muon729",10,True,True)     #
# makeProd("muon730",10,'Jpsi',False)  # made with E50
# makeProd("muon731",10,'Jpsi',True)   # made with E50
# makeProd("muon732",10,'Jpsi',True)   # made with E50
# makeProd("muon733",10,'Jpsi',True)   # back to pencil beam
# makeProd("muon630",10,False,True) # test with new muonShield code, 3cm smearing
# makeProd("muon631",10,False,True) # run with concrete wall enabled as sensitive
# makeProd("muon632",10,False,True) # run with concrete wall enabled as sensitive, active shielding polarity fixed
# but wrong geometry
# makeProd("muon810",10,False,False) # start production with latest geometry
# makeProd("muon820",10,True,False)
# makeProd("muon811",10,False,True) #
# makeProd("muon821",10,True,True)