ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/RateMonShiftTool_dev/AddTableInfo_db.py
(Generate patch)

Comparing UserCode/RateMonShiftTool_dev/AddTableInfo_db.py (file contents):
Revision 1.17 by grchrist, Mon Mar 26 15:50:19 2012 UTC vs.
Revision 1.28 by amott, Thu May 10 14:27:53 2012 UTC

# Line 9 | Line 9 | NHighStreamA=0
9  
10   def MoreTableInfo(parser,LumiRange,config,isCol=True):
11      print "Monitoring Run %d" % (parser.RunNumber,)
12 <    [AvInstLumi, AvLiveLumi, AvDeliveredLumi, AvDeadTime,PSCols] = parser.GetAvLumiInfo(LumiRange)
13 <    deadtimebeamactive=parser.GetDeadTimeBeamActive(LumiRange)
12 >    print "len=",len(LumiRange)
13 >    print "LSRange=", LumiRange
14 >    if len(LumiRange)>0:
15 >        
16 >        [AvInstLumi, AvLiveLumi, AvDeliveredLumi, AvDeadTime,PSCols] = parser.GetAvLumiInfo(LumiRange)
17 >        deadtimebeamactive=parser.GetDeadTimeBeamActive(LumiRange)*100
18 >        ##print "dtba=",deadtimebeamactive
19 >    else:
20 >        print "no lumisections to monitor"
21 >        return
22      try:
23          LastPSCol = PSCols[-1]
24      except:
25          LastPSCol = -1
26 +    if isCol:
27          
28 <    aRates = parser.GetTriggerRatesByLS("AOutput")
28 >        aRates = parser.GetTriggerRatesByLS("AOutput")
29 >        bRates = parser.GetTriggerRatesByLS("BOutput")
30 >        if len(bRates) == 0:
31 >            realARates = aRates
32 >        else:
33 >            realARates={}
34 >            for k,v in bRates.iteritems():
35 >                realARates[k]=aRates[k]-bRates[k]
36 >                #realARates = aRates - bRates*20;
37 >    else:
38 >        if len(parser.GetTriggerRatesByLS("AOutput"))>0:
39 >            aRates = parser.GetTriggerRatesByLS("AOutput")
40 >        else:
41 >            aRates = parser.GetTriggerRatesByLS("AForPPOutput")
42 >            
43 >    
44      expressRates = {}
45      if isCol:
46          expressRates = parser.GetTriggerRatesByLS("ExpressOutput")
47      else:
48 <        expressRates = parser.GetTriggerRatesByLS("ExpressCosmicsOutput")
48 >        if len(parser.GetTriggerRatesByLS("ExpressOutput"))>0:
49 >            expressRates=parser.GetTriggerRatesByLS("ExpressOutput")
50 >        else:
51 >            expressRates = parser.GetTriggerRatesByLS("ExpressForCosmicsOutput")
52      ExpRate=0
53      PeakRate=0
54      AvgExpRate=0
# Line 30 | Line 57 | def MoreTableInfo(parser,LumiRange,confi
57      PeakRateA=0
58      AvgRateA=0
59      
60 +    realARate=0
61 +    realPeakRateA=0
62 +    realAvgRateA=0
63 +    
64      if len(expressRates.values()) > 20:
65          AvgExpRate = sum(expressRates.values())/len(expressRates.values())
66  
# Line 43 | Line 74 | def MoreTableInfo(parser,LumiRange,confi
74          ARate+=thisRateA
75          if thisRateA>PeakRateA:
76              PeakRateA=thisRateA
77 <        
77 >
78 >        thisRealRateA = aRates.get(ls,0) - bRates.get(ls,0)
79 >        realARate+=thisRealRateA
80 >        if thisRealRateA > realPeakRateA:
81 >            realReakRateA = thisRealRateA
82          #ARate+=aRates.get(ls,0)
83      ## Print Stream A Rate --moved see below
84      ##print "Current Steam A Rate is: %0.1f Hz" % (ARate/len(LumiRange),)
# Line 65 | Line 100 | def MoreTableInfo(parser,LumiRange,confi
100          
101      write(text)
102      if badExpress:
103 <        if (ExpRate-PeakRate)/(len(LumiRange)-1) <=config.MaxExpressRate: ## one lumisection causes this
104 <            write("  <<  This appears to be due to a 1 lumisection spike, please monitor\n")
105 <        else:
106 <            if NHighExpress > 1:  # big problem, call HLT DOC
107 <                write(colored("  <<  WARNING: Current Express rate is too high!",'red',attrs=['reverse']) )
108 <                Warn = True
103 >        if len(LumiRange)>1:
104 >            if (ExpRate-PeakRate)/(len(LumiRange)-1) <=config.MaxExpressRate: ## one lumisection causes this
105 >                write("  <<  This appears to be due to a 1 lumisection spike, please monitor\n")
106 >            else:
107 >                if NHighExpress > 1:  # big problem, call HLT DOC
108 >                    write(colored("  <<  WARNING: Current Express rate is too high!",'red',attrs=['reverse']) )
109 >                    Warn = True
110  
111                  #    if AvgExpRate > config.MaxExpressRate:
112                  #        write( colored("\n\nWARNING: Average Express Stream Rate is too high (%0.1f Hz)  << CALL HLT DOC" % AvgExpRate,'red',attrs=['reverse']) )
113                  #        Warn = True
114 <    
114 >        
115 >            
116  
117  
118      #########################################
119      ##Check if Stream A is too high
120      #########################################
121      global NHighStreamA
122 <    badStreamA =ARate/len(LumiRange) > config.MaxStreamARate ##Cosmics Express Rate 300 Hz max
122 >    badStreamA =realARate/len(LumiRange) > config.MaxStreamARate ##Cosmics Express Rate 300 Hz max
123      baseTextA= "\nCurrent Steam A Rate is: %0.1f Hz" % (ARate/len(LumiRange),)
124 +    baseTextRealA= "\nCurrent PROMPT Steam A Rate is: %0.1f Hz" % (realARate/len(LumiRange),)
125      if badStreamA:
126          textA=colored(baseTextA,'red',attrs=['reverse'])  ## bad, make the text white on red
127 +        textRealA=colored(baseTextRealA,'red',attrs=['reverse'])  ## bad, make the text white on red
128          NHighStreamA+=1
129      else:
130          textA=baseTextA
131 +        textRealA=baseTextRealA
132          NHighStreamA=0
133  
134      write(textA)
135 +    write(textRealA)
136      if badStreamA:
137 <        if (ARate-PeakRateA)/(len(LumiRange)-1) <=config.MaxStreamARate: ## one lumisection causes this
138 <            write("  <<  This appears to be due to a 1 lumisection spike, please monitor\n")
139 <        else:
140 <            if NHighStreamA >1: ##Call HLT doc!
141 <                write(colored("  <<  WARNING: Current Stream A rate is too high!",'red',attrs=['reverse']) )
142 <                Warn = True
137 >        if len(LumiRange)>1:
138 >            if (realARate-realPeakRateA)/(len(LumiRange)-1) <=config.MaxStreamARate: ## one lumisection causes this
139 >                write("  <<  This appears to be due to a 1 lumisection spike, please monitor\n")
140 >            else:
141 >                if NHighStreamA >1: ##Call HLT doc!
142 >                    write(colored("  <<  WARNING: Current Stream A rate is too high!",'red',attrs=['reverse']) )
143 >                    Warn = True
144      write("\n\n")
145              
146      ######################################
# Line 135 | Line 177 | def MoreTableInfo(parser,LumiRange,confi
177              write(bcolors.OKBLUE)
178          write(str(round(deadtimebeamactive,2))+"%")
179          write(bcolors.ENDC+"\n")
180 <
181 <    print "Used prescale column(s): "+str(PrescaleColumnString)    
182 <    write("Lumisections: ")
180 >    write("Used prescale column(s): %s  " % (str(PrescaleColumnString),) )
181 >    if LastPSCol in config.ForbiddenCols and isCol:
182 >        write( colored("<< Using column %d! Please check in the documentation that this is the correct column" % (LastPSCol),'red',attrs=['reverse']) )
183 >    write("\nLumisections: ")
184      if not isSequential(LumiRange):
185          write(str(LumiRange)+"   Lumisections are not sequential (bad LS skipped)\n")
186      else:
187          write("%d - %d\n" % (min(LumiRange),max(LumiRange),))
188 <    print "\nLast Lumisection of the run is:        "+str(parser.GetLastLS())
189 <    write(  "Last Lumisection good for physics is:  "+str(parser.GetLastLS(True)) )
188 >    ##print "\nLast Lumisection of the run is:        "+str(parser.GetLastLS())
189 >    write(  "\nLast Lumisection good where DAQ is active is:  "+str(parser.GetLastLS(isCol)) )
190 >    ##write(  "Last Lumisection where DAQ is active is:  "+str(parser.GetLastLS(True)) )
191      write("\n\n\n")
192  
193 <    if isCol:
194 <        L1RatePredictions = config.GetExpectedL1Rates(AvInstLumi)
195 <        if len(L1RatePredictions):
196 <            print "Expected Level 1 Rates:"
197 <        for key,val in L1RatePredictions.iteritems():
198 <            write("Prescale Column "+str(key)+":  "+str(round(val/1000,1))+" kHz")
199 <            if key == LastPSCol:
200 <                write(' << taking data in this column')
201 <            write('\n')
193 >    ## if isCol:
194 > ##         L1RatePredictions = config.GetExpectedL1Rates(AvInstLumi)
195 > ##         if len(L1RatePredictions):
196 > ##             print "Expected Level 1 Rates:"
197 > ##         for key,val in L1RatePredictions.iteritems():
198 > ##             write("Prescale Column "+str(key)+":  "+str(round(val/1000,1))+" kHz")
199 > ##             if key == LastPSCol:
200 > ##                 write(' << taking data in this column')
201 > ##             write('\n')
202          
203      
204  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines