In This Package:

mixing.py

Go to the documentation of this file.

#!/usr/bin/env python
"""
nuwa.py -A none -l 4 -n 1000 mixing.py

"""

inputpath = ""
outputpath = ""

import Mixing
mixing=Mixing.Configure()

#Set the samples being mixed
mixing.MixSvc.SampleInputs={
            "IBD1":[inputpath+"IBD_DayaBayAD1.root"],
            "IBD2":[inputpath+"IBD_DayaBayAD1.root"],
            "RAD1":[inputpath+"Co60_PMT_DayaBayAD1.root"], 
            "RAD2":[inputpath+"Co60_PMT_DayaBayAD2.root"],
            "RAD3":[inputpath+"Co60_sst_DayaBayAD1.root"],
            "RAD4":[inputpath+"Co60_sst_DayaBayAD2.root"],
            "RAD5":[inputpath+"K40_gds_DayaBayAD1.root"], 
            "RAD6":[inputpath+"K40_gds_DayaBayAD2.root"],
            "RAD7":[inputpath+"K40_lso_DayaBayAD1.root"],
            "RAD8":[inputpath+"K40_lso_DayaBayAD2.root"],
            "RAD9":[inputpath+"K40_PMT_DayaBayAD1.root"],
            "RAD10":[inputpath+"K40_PMT_DayaBayAD2.root"],
            "RAD11":[inputpath+"K40_sst_DayaBayAD1.root"],
            "RAD12":[inputpath+"K40_sst_DayaBayAD2.root"],
            "RAD13":[inputpath+"Th232_gds_DayaBayAD1.root"], 
            "RAD14":[inputpath+"Th232_gds_DayaBayAD2.root"],
            "RAD15":[inputpath+"Th232_lso_DayaBayAD1.root"],
            "RAD16":[inputpath+"Th232_lso_DayaBayAD2.root"],
            "RAD17":[inputpath+"Th232_PMT_DayaBayAD1.root"],
            "RAD18":[inputpath+"Th232_PMT_DayaBayAD2.root"],
            "RAD19":[inputpath+"Th232_sst_DayaBayAD1.root"],
            "RAD20":[inputpath+"Th232_sst_DayaBayAD2.root"],
            "RAD21":[inputpath+"U238_gds_DayaBayAD1.root"], 
            "RAD22":[inputpath+"U238_gds_DayaBayAD2.root"],
            "RAD23":[inputpath+"U238_lso_DayaBayAD1.root"],
            "RAD24":[inputpath+"U238_lso_DayaBayAD2.root"],
            "RAD25":[inputpath+"U238_PMT_DayaBayAD1.root"],
            "RAD26":[inputpath+"U238_PMT_DayaBayAD2.root"],
            "RAD27":[inputpath+"U238_sst_DayaBayAD1.root"],
            "RAD28":[inputpath+"U238_sst_DayaBayAD2.root"]
         }
if mixing.MixSvc.SampleInputs:
    print 'input_samples=', str(mixing.MixSvc.SampleInputs)

#Set the event rates of the samples    
mixing.MixSvc.SampleTaus={"IBD1":86400./10000,
                        "IBD2": 86400./2000,
                        "RAD1": 86400./2000,
                        "RAD2": 86400./2000,
                        "RAD3": 86400./2000,
                        "RAD4": 86400./2000,
                        "RAD5": 86400./2000,
                        "RAD6": 86400./2000,
                        "RAD7": 86400./2000,
                        "RAD8": 86400./2000,
                        "RAD9": 86400./2000,
                        "RAD10": 86400./2000,
                        "RAD11": 86400./2000,
                        "RAD12": 86400./2000,
                        "RAD13": 86400./2000,
                        "RAD14": 86400./2000,
                        "RAD15": 86400./2000,
                        "RAD16": 86400./2000,
                        "RAD17": 86400./2000,
                        "RAD18": 86400./2000,
                        "RAD19": 86400./2000,
                        "RAD20": 86400./2000,
                        "RAD21": 86400./2000,
                        "RAD22": 86400./2000,
                        "RAD23": 86400./2000,
                        "RAD24": 86400./2000,
                        "RAD25": 86400./2000,
                        "RAD26": 86400./2000,
                        "RAD27": 86400./2000,
                        "RAD28": 86400./2000}
if mixing.MixSvc.SampleTaus:
    print 'sample_deltaT=', str(mixing.MixSvc.SampleTaus)

#Set the time window of mixing
mixing.MixSvc.MixingTime= 10000
if mixing.MixSvc.MixingTime:
    print 'Mixing Time=', mixing.MixSvc.MixingTime

#Set the begin time of this time
mixing.MixSvc.BeginTime=0.0
print 'Begin Time = ', mixing.MixSvc.BeginTime

mixing.MixSvc.DefaultOutput=outputpath+"mixing.root"

---