Functions
	LHCBeamSpotSize (lhcTimeString, lhcBeta1, lhcEmit1, lhcBeta2, lhcEmit2, betaGamma=3730. # Default value of betaGamma for 7 TeV)
	LHCEmittance (lhcTimeString, lhcEmit1, lhcEmit2, minEmittance=0., maxEmittance=8.)
	LHCBeamSpotSizeFill1022 (plot='sigx')
	LHCEmittanceFill1022 (plot='sigx')
	LHCBeamSpotSizeFill1058 (plot='sigx')

Variables
str	__author__ = 'Juerg Beringer'
str	__version__ = '$Id $'
list	fill1022TimeString = ['Apr 4 14:27:30 2010', 'Apr 4 22:27:30 2010', 'Apr 5 3:57:30 2010','Apr 5 7:57:30 2010' ]
	fill1022Emit1X = array([ 1.45, 1.8, 1.92, 2.0 ])
	fill1022Emit1Y = array([ 1.6, 2.1, 2.45, 2.55 ])
	fill1022Emit2X = array([ 1.8, 2.27, 2.27, 2.5 ])
	fill1022Emit2Y = array([ 2.6, 4.6, 5.4, 5.8 ])
float	fill1022Beta1X = 11.54
float	fill1022Beta1Y = 12.74
float	fill1022Beta2X = 9.89
float	fill1022Beta2Y = 12.70
list	fill1058TimeString = ['Apr 25 6:03:00 2010', 'Apr 24 6:57:00 2010']
	fill1058Emit1X = array([ 5.05, 2.6])
	fill1058Emit1Y = array([ 5.04, 2.6])
	fill1058Emit2X = array([ 3.88, 2.5])
	fill1058Emit2Y = array([10.64, 4.9])
float	fill1058Beta1X = 2.28
float	fill1058Beta1Y = 2.02
float	fill1058Beta2X = 1.92
float	fill1058Beta2Y = 2.10
	gr = LHCBeamSpotSizeFill1022()

Detailed Description

A library with miscellaneous beam spot related utilities and plots.

Function Documentation

◆ LHCBeamSpotSize()

python.PlotLibrary.LHCBeamSpotSize	(	lhcTimeString,
		lhcBeta1,
		lhcEmit1,
		lhcBeta2,
		lhcEmit2,
		betaGamma = 3730. # Default value of betaGamma for 7 TeV )

Definition at line 23 of file PlotLibrary.py.

                    ):
    lhcTime = array([ mktime(strptime(s,'%b %d %H:%M:%S %Y')) for s in lhcTimeString])
    print (lhcTime)
 
    # Normalized emittance quoted in um rad; lhcEps arrays are emittance in m rad
    lhcEps1 = lhcEmit1 * 1.0e-6 / betaGamma
    lhcEps2 = lhcEmit2 * 1.0e-6 / betaGamma
 
    # Beam sizes in mm
    lhcSigmaBeam1 = sqrt(lhcBeta1*lhcEps1) * 1.0e3
    print ('Size of beam 1 (mm): ',lhcSigmaBeam1)
    lhcSigmaBeam2 = sqrt(lhcBeta2*lhcEps2) * 1.0e3
    print ('Size of beam 2 (mm): ',lhcSigmaBeam2)
 
    # Beam spot size in mm
    lhcSigma = 1. / sqrt(1./lhcSigmaBeam1**2 + 1./lhcSigmaBeam2**2)
    print ('Beam spot size (mm): ',lhcSigma)
 
    # Draw on existing plot
    gr = ROOTUtils.protect( ROOT.TGraph(len(lhcTime),lhcTime,lhcSigma) )
    gr.SetMarkerColor(2)
    gr.SetLineColor(2)
    #gr.SetMarkerSize(1.0)
    #gr.SetMarkerStyle(25)
    gr.SetMarkerSize(1.5)
    gr.SetMarkerStyle(20)
    gr.Draw('PSAME')
    return [(gr,'Expected from  #sqrt{#epsilon #beta*}','P')]
 
 
#
# Plot normalized emittance with new axis on right hand side of plot
#

◆ LHCBeamSpotSizeFill1022()

python.PlotLibrary.LHCBeamSpotSizeFill1022 ( plot = 'sigx' )

Definition at line 128 of file PlotLibrary.py.

def LHCBeamSpotSizeFill1022(plot='sigx'):
    if plot=='sigx':
        return LHCBeamSpotSize(fill1022TimeString,fill1022Beta1X,fill1022Emit1X,fill1022Beta2X,fill1022Emit2X)
    if plot=='sigy':
        return LHCBeamSpotSize(fill1022TimeString,fill1022Beta1Y,fill1022Emit1Y,fill1022Beta2Y,fill1022Emit2Y)
 

◆ LHCBeamSpotSizeFill1058()

python.PlotLibrary.LHCBeamSpotSizeFill1058 ( plot = 'sigx' )

Definition at line 155 of file PlotLibrary.py.

def LHCBeamSpotSizeFill1058(plot='sigx'):
    if plot=='sigx':
        return LHCBeamSpotSize(fill1058TimeString,fill1058Beta1X,fill1058Emit1X,fill1058Beta2X,fill1058Emit2X)
    if plot=='sigy':
        return LHCBeamSpotSize(fill1058TimeString,fill1058Beta1Y,fill1058Emit1Y,fill1058Beta2Y,fill1058Emit2Y)
 
    
 
# Test code for modules

◆ LHCEmittance()

python.PlotLibrary.LHCEmittance	(	lhcTimeString,
		lhcEmit1,
		lhcEmit2,
		minEmittance = 0.,
		maxEmittance = 8. )

Definition at line 62 of file PlotLibrary.py.

                 ):
    lhcTime = array([ mktime(strptime(s,'%b %d %H:%M:%S %Y')) for s in lhcTimeString])
    ROOT.gPad.SetTicky(0)   # Remove tick marks on right side y axis
    ROOT.gPad.SetRightMargin(0.1)
    #maxEmittance = max(max(lhcEmit1),max(lhcEmit2))*1.7
    hmin = ROOT.gPad.GetUymin()
    hmax = ROOT.gPad.GetUymax()
 
    lhcEmit1Data = (lhcEmit1-minEmittance)/(maxEmittance-minEmittance)*(hmax-hmin)+hmin
    gr1 = ROOTUtils.protect( ROOT.TGraph(len(lhcTime),lhcTime,lhcEmit1Data) )
    gr1.SetMarkerColor(4)
    gr1.SetLineColor(4)
    gr1.SetMarkerStyle(25)
    gr1.Draw('PSAME')
    #gr1.Draw('PLSAME')
 
    lhcEmit2Data = (lhcEmit2-minEmittance)/(maxEmittance-minEmittance)*(hmax-hmin)+hmin
    gr2 = ROOTUtils.protect( ROOT.TGraph(len(lhcTime),lhcTime,lhcEmit2Data) )
    gr2.SetMarkerColor(4)
    gr2.SetLineColor(4)
    gr2.SetMarkerStyle(26)
    gr2.Draw('PSAME')
    #gr2.Draw('PLSAME')
 
    raxis = ROOTUtils.protect( ROOT.TGaxis(ROOT.gPad.GetUxmax(),hmin,ROOT.gPad.GetUxmax(),hmax,minEmittance,maxEmittance,510,'+L') )
    raxis.SetLineColor(4)
    raxis.SetTitle('Invariant Emittance (#mum rad)')
    raxis.SetTitleColor(4)
    raxis.SetTitleSize(0.05)
    raxis.SetTextFont(42)
    raxis.SetLabelFont(42)
    raxis.SetLabelSize(0.05)
    raxis.SetLabelColor(4)
    raxis.Draw()
    return [(gr1,'Emittance beam 1','P'),(gr2,'Emittance beam 2','P')]
 
 
#
# Plots and data for fill 1022
#

◆ LHCEmittanceFill1022()

python.PlotLibrary.LHCEmittanceFill1022 ( plot = 'sigx' )

Definition at line 134 of file PlotLibrary.py.

def LHCEmittanceFill1022(plot='sigx'):
    if plot=='sigx':
        return LHCEmittance(fill1022TimeString,fill1022Emit1X,fill1022Emit2X,0.1,8.5)
    if plot=='sigy':
        return LHCEmittance(fill1022TimeString,fill1022Emit1Y,fill1022Emit2Y,0.1,8.5)
 
 
 
#
# Plots and data for fill 1058
#

Variable Documentation

◆ author

str python.PlotLibrary.__author__ = 'Juerg Beringer'

private

Definition at line 8 of file PlotLibrary.py.

◆ version

str python.PlotLibrary.__version__ = '$Id $'

private

Definition at line 9 of file PlotLibrary.py.

◆ fill1022Beta1X

float python.PlotLibrary.fill1022Beta1X = 11.54

Definition at line 112 of file PlotLibrary.py.

◆ fill1022Beta1Y

float python.PlotLibrary.fill1022Beta1Y = 12.74

Definition at line 113 of file PlotLibrary.py.

◆ fill1022Beta2X

float python.PlotLibrary.fill1022Beta2X = 9.89

Definition at line 114 of file PlotLibrary.py.

◆ fill1022Beta2Y

float python.PlotLibrary.fill1022Beta2Y = 12.70

Definition at line 115 of file PlotLibrary.py.

◆ fill1022Emit1X

python.PlotLibrary.fill1022Emit1X = array([ 1.45, 1.8, 1.92, 2.0 ])

Definition at line 108 of file PlotLibrary.py.

◆ fill1022Emit1Y

python.PlotLibrary.fill1022Emit1Y = array([ 1.6, 2.1, 2.45, 2.55 ])

Definition at line 109 of file PlotLibrary.py.

◆ fill1022Emit2X

python.PlotLibrary.fill1022Emit2X = array([ 1.8, 2.27, 2.27, 2.5 ])

Definition at line 110 of file PlotLibrary.py.

◆ fill1022Emit2Y

python.PlotLibrary.fill1022Emit2Y = array([ 2.6, 4.6, 5.4, 5.8 ])

Definition at line 111 of file PlotLibrary.py.

◆ fill1022TimeString

list python.PlotLibrary.fill1022TimeString = ['Apr 4 14:27:30 2010', 'Apr 4 22:27:30 2010', 'Apr 5 3:57:30 2010','Apr 5 7:57:30 2010' ]

Definition at line 107 of file PlotLibrary.py.

◆ fill1058Beta1X

float python.PlotLibrary.fill1058Beta1X = 2.28

Definition at line 150 of file PlotLibrary.py.

◆ fill1058Beta1Y

float python.PlotLibrary.fill1058Beta1Y = 2.02

Definition at line 151 of file PlotLibrary.py.

◆ fill1058Beta2X

float python.PlotLibrary.fill1058Beta2X = 1.92

Definition at line 152 of file PlotLibrary.py.

◆ fill1058Beta2Y

float python.PlotLibrary.fill1058Beta2Y = 2.10

Definition at line 153 of file PlotLibrary.py.

◆ fill1058Emit1X

python.PlotLibrary.fill1058Emit1X = array([ 5.05, 2.6])

Definition at line 146 of file PlotLibrary.py.

◆ fill1058Emit1Y

python.PlotLibrary.fill1058Emit1Y = array([ 5.04, 2.6])

Definition at line 147 of file PlotLibrary.py.

◆ fill1058Emit2X

python.PlotLibrary.fill1058Emit2X = array([ 3.88, 2.5])

Definition at line 148 of file PlotLibrary.py.

◆ fill1058Emit2Y

python.PlotLibrary.fill1058Emit2Y = array([10.64, 4.9])

Definition at line 149 of file PlotLibrary.py.

◆ fill1058TimeString

list python.PlotLibrary.fill1058TimeString = ['Apr 25 6:03:00 2010', 'Apr 24 6:57:00 2010']

Definition at line 145 of file PlotLibrary.py.

◆ gr

python.PlotLibrary.gr = LHCBeamSpotSizeFill1022()

Definition at line 165 of file PlotLibrary.py.

Functions

Variables