In This Package:

Pool.py

Go to the documentation of this file.

#!/usr/bin/env python
'''
Configuration helpers for pool validatation simulations.
'''

class PoolMuons:
    ''' Default configuration for sending muons through a pool.  This
    makes use of an external HepEVT muon generator, turns off
    scintilation to increase processing speed and sets up for the
    validation ntuple.  It is based on the original ran*Muon.py
    scripts.
    '''

    def __init__(self,site="far",histogram_filename = '../mc_ntuple/%sMuonPool.root'):
        ''' Construct the default configuration for muon pool
        validation.  The site can be set to "far" or "near" (dayabay
        near site, near sites non-identicality was not considered).  A
        non-default histogram_filename can be specified to hold the
        resulting histograms.  If a "%s" is included the site will be
        substituted.
        
        Warning: this hard codes generation of 1000 muons.  Simulating
        less is okay, but more is bad.
        '''

        # key off site.
        volume = "/dd/Structure/Pool/db-ows"
        musicsite = 'DYB'
        if 'far' == site:
            volume = "/dd/Structure/Pool/far-ows"
            musicsite = 'Far'
        else:
            site = 'dayabay'    # normalize

        # let any %s be dereferenced.
        try:
            histogram_filename = histogram_filename%site
        except TypeError:
            pass

        import GaudiKernel.SystemOfUnits as units

        # Try to find MUSIC's data directory 
        import os
        muonroot = os.getenv('MUONROOT',None)
        if not muonroot:
            print "Can not find $MUONROOT needed to locate MUSIC directory."
            import sys
            sys.exit(1)
            pass
        
        source = "Muon.exe -n 1000 -s %s -music_dir %s/data|"%(musicsite,
                                                               muonroot)

        from GenTools.Helpers import HepEVT
        hepevt = HepEVT(source)
        hepevt.timerator.LifeTime = 1*units.second
        hepevt.transformer.Volume = volume
        # Geometry impedeance mismatch between Muon.exe and XmlDetDesc
        hepevt.transformer.Offset = [0., 0., (0.042)*units.meter]

        import GenTools
        gt = GenTools.Configure(helper=hepevt)


        # Configure DetSim

        import DetSim
        detsim = DetSim.Configure(site)
        
        from DetSim.DetSimConf import DsPhysConsOptical
        optical = DsPhysConsOptical()
        optical.CerenPhotonScaleWeight = 3.0
        optical.ScintPhotonScaleWeight = 3.0

        # use temps to keep lines succinct
        ismuon = "(pdg==13 or pdg==-13) and "
        material = "MaterialName == '/dd/Materials/%s"
        stats = [
            ["MuonTrkLengthInOws","dx",ismuon+material%"OwsWater'"],
            ["MuonTrkLengthInIws", "dx", ismuon+material%"IwsWater'"],
            ["MuonTrkLengthInLS", "dx", ismuon+material%"LiquidScintillator'"],
            ["MuonTrkLengthInGdLS","dx", ismuon+material%"GdDopedLS'"],
            ["MuonStop", "dx", ismuon+" (IsStopping>0)"],    
            ]
        detsim.unobserver(stats=stats)
        detsim.historian(trackSelection = "pdg==13 or pdg==-13")

        # Finally, DetSimVali itself
        from DetSimValidation.DetSimValidationConf import DetSimVali
        dsv = DetSimVali()
        self.detsimvali = dsv
        
        from Gaudi.Configuration import ApplicationMgr
        theApp = ApplicationMgr()
        theApp.TopAlg.append(dsv)

        from GaudiSvc.GaudiSvcConf import THistSvc
        histsvc = THistSvc()
        histsvc.Output =["file1 DATAFILE='%s' OPT='RECREATE' TYP='ROOT' "%histogram_filename]
        return

if '__main__' == __name__:
    from Gaudi.Configuration import ApplicationMgr
    theApp = ApplicationMgr()

    from DybPython.Control import main
    nuwa = main()

    import XmlDetDesc
    XmlDetDesc.Configure()
    
    from DetSimValidation import PoolMuons
    pool = PoolMuons("far")   # "far" by default otherwise will do far

    print theApp.TopAlg
    nuwa.run()

---