DeadROCCounter_Phase1 Namespace Reference

Functions
def	countBadROCBarrel
	barrel######################################################## More...
def	countBadROCForward
	endacp######################################### More...
def	getRunNumber

Variables
int	bpix_tot_deadROC = 0
int	bpix_tot_ineffROC = 0
int	bpix_tot_Nrocspopulated = 0
int	bpix_tot_totalentries = 0
string	commonPath = "DQMData/Run "
tuple	fin = TFile(fname)
list	fname = sys.argv[1]
	main####################################### More...
int	fpix_tot_deadROC = 0
int	fpix_tot_ineffROC = 0
int	fpix_tot_Nrocspopulated = 0
int	fpix_tot_totalentries = 0
tuple	hnpixclus_bpix = fin.Get(commonPath + "charge_PXBarrel")
tuple	hnpixclus_fpix = fin.Get(commonPath + "charge_PXForward")
tuple	out_file = open(outname, "w")
string	outname = "PixZeroOccROCs_run"

Function Documentation

def DeadROCCounter_Phase1.countBadROCBarrel	(		fin,
			layerNo,
			os
	)

barrel########################################################

Definition at line 14 of file DeadROCCounter_Phase1.py.

References if(), print(), sistrip::SpyUtilities.range(), and str.

def countBadROCBarrel(fin, layerNo, os):
    global bpix_tot_deadROC
    global bpix_tot_ineffROC
    global bpix_tot_Nrocspopulated
    global bpix_tot_totalentries

    barrelPath  = commonPath + "PXBarrel/";
    histoname = ["digi_occupancy_per_SignedModuleCoord_per_SignedLadderCoord_PXLayer_1", "digi_occupancy_per_SignedModuleCoord_per_SignedLadderCoord_PXLayer_2",  
                 "digi_occupancy_per_SignedModuleCoord_per_SignedLadderCoord_PXLayer_3", "digi_occupancy_per_SignedModuleCoord_per_SignedLadderCoord_PXLayer_4"]
    digi2D = fin.Get(barrelPath + histoname[layerNo-1])
    #return status flag is histogram is empty!
    if digi2D.GetEntries() == 0 :  
        return 1;
    Nrocspopulated = 0
    totalEntries = 0
    NexpectedROC = [1536, 3584, 5632, 8192]
    nLadders_bpx = [6, 14, 22, 32]
    nx = digi2D.GetNbinsX()
    ny = digi2D.GetNbinsY()
    for xbin in range(1,nx+1):
        if xbin >= 33 and xbin <= 40:    continue;#region of cross on x-axis
        for ybin in range(1,ny+1):
            if (ybin == 2*nLadders_bpx[layerNo-1] + 1) or (ybin == 2*nLadders_bpx[layerNo-1] + 2):    continue;#region of cross on y-axis
            bentries = digi2D.GetBinContent(xbin,ybin)
            if(bentries > 0):
                Nrocspopulated+=1
                totalEntries += bentries
    meanEntries = float(totalEntries)/Nrocspopulated
    NineffROC = 0
    #Loop to chek inefficient ROC per layer
    for xbin in range(1,nx+1):
        if xbin >= 33 and xbin <= 40:
            continue;#region of cross on x-axis
        for ybin in range(1,ny+1):
            if (ybin == 2*nLadders_bpx[layerNo-1] + 1) or (ybin == 2*nLadders_bpx[layerNo-1] + 2):    continue;#region of cross on y-axis
            bentries = digi2D.GetBinContent(xbin,ybin);
            if(bentries > 0 and bentries < meanEntries/4. ):#Assume < 25% of MEAN = inefficient
                NineffROC+=1;
    ##Printing Layer no., #dead ROC, #inefficienct ROC, #mean occupancy of Non-zer roc
    tmpstr = "BPix L" + str(layerNo)
    print(tmpstr, '{0:4d} {1:4d} {2:4.1f}'.format(NexpectedROC[layerNo-1] - Nrocspopulated, NineffROC, round(meanEntries,1)), file=os)
    bpix_tot_deadROC += NexpectedROC[layerNo-1] - Nrocspopulated
    bpix_tot_ineffROC += NineffROC
    bpix_tot_Nrocspopulated += Nrocspopulated
    bpix_tot_totalentries += float(totalEntries)