In This Package:

Opts.py

Go to the documentation of this file.

#!/usr/bin/env python

'''
Configure a DetSim from the command line.

This does not set up kinematics.  You can run it something like:

nuwa.py -n 10 -o output.root -m 'GenTools.Opts [opts]' -m 'DetSim.Opts [opts]'

Passing '--help' as an option will show what options the modules
provide.
'''

def configure(argv=[]):
    
    from optparse import OptionParser
    parser = OptionParser(usage=__doc__)
    parser.add_option("-p", "--physics-list",default="all",type="string",
                      help='Set physics list.  Can be "basic", "all" or list of physics constructor classes.')
    parser.add_option("-s", "--site",default="all",type="string",
                      help='Set what site geometry to load, can be "far", "dayabay", "lingao" or "all"')
    parser.add_option("-t", "--track-selection",default="",type="string",
                      help='Set the track selection rules for the particle history')
    parser.add_option("-v", "--vertex-selection",default="",type="string",
                      help='Set the vertex selection rules for the particle history')
    parser.add_option("-u", "--unobservable-statistics",default="[]",type="string",
                      help='Set the rules for the unobservable statistics, takes a list or tripplets.')
    parser.add_option("-V", "--verbose",default=False,
                      help='Turn up Geant4 verbosity.  Big log files!')

    (opts,args) = parser.parse_args(args=argv)

    from DetSim import Configure, physics_list_basic, physics_list_nuclear

    if opts.physics_list == "basic":
        opts.physics_list = physics_list_basic
    elif opts.physics_list == "all":
        opts.physics_list = physics_list_basic+physics_list_nuclear
    # convert string rep of list of strings to Python list of strings
    elif opts.physics_list: 
        opts.physics_list = eval(opts.physics_list)
    else:
        raise ValueError,'Unknown physics list: "%s"'%opts.physics_list
        
    # User should pass in a string rep of a list of tripplets of
    # strings.  Need to convert it to a Python list of tripplets of
    # strings.
    opts.unobservable_statistics = eval(opts.unobservable_statistics)

    if opts.site == 'All' or opts.site == 'all': opts.site='far,dayabay,lingao'

    print 'Site is',opts.site
    detsim = Configure(site=opts.site,physlist=opts.physics_list)

    if opts.track_selection or opts.vertex_selection:
        detsim.historian(opts.track_selection, opts.vertex_selection)
    if opts.unobservable_statistics:
        detsim.unobserver(opts.unobservable_statistics)

    if opts.verbose:
        # Turn up verbosity
        from GaussTools.GaussToolsConf import GiGaEventActionCommand
        geac = GiGaEventActionCommand("GiGa.GiGaEventActionCommand")
        geac.BeginOfEventCommands = [
            "/control/verbose 2",
            "/run/verbose 1",
            "/event/verbose 2",
            "/tracking/verbose 2",
            "/geometry/navigator/verbose 2"
            ]
        detsim.giga.EventAction = geac

    return

    
if '__main__' == __name__:
    #configure(["--help"])
    configure()

---