File indexing completed on 2023-03-17 10:56:40
0001
0002
0003 from __future__ import print_function
0004 import sys
0005 import string
0006 from ROOT import *
0007 from array import array
0008
0009 def getFileInPath(rfile):
0010 import os
0011 for dir in os.environ['CMSSW_SEARCH_PATH'].split(":"):
0012 if os.path.exists(os.path.join(dir,rfile)): return os.path.join(dir,rfile)
0013 return None
0014
0015
0016 detIDsFileName = getFileInPath('DQM/SiStripMonitorClient/data/detids.dat')
0017
0018 filename_online=sys.argv[1]
0019 filename_offline=sys.argv[2]
0020
0021 runNum=filename_offline[19:25]
0022
0023 dir="DQMData/Run " + runNum + "/PixelPhase1/Run summary/Pahse1_MechanicalView/"
0024
0025
0026 dirBPix=dir + "PXBarrel/"
0027 dirFPix=dir + "PXForward/"
0028
0029 hnameB="digi_occupancy_per_SignedModuleCoord_per_SignedLadderCoord_PXLayer_"
0030
0031 minlad=[-6,-14,-22,-32]
0032
0033 shell=""
0034
0035
0036 def BPIX_list(inputFile):
0037 DQMfile=TFile(inputFile)
0038 BPIXCounter = []
0039 BPIXCounter_v0 = []
0040 for l in range(1,5):
0041 hname=hnameB + str(l)
0042
0043 RocMap=DQMfile.FindObjectAny(hname)
0044
0045 for j in range(1,RocMap.GetNbinsY()+1):
0046
0047 lad=minlad[l-1] + int(j-1)/2
0048 alad=abs(lad)
0049
0050 for i in range(1,RocMap.GetNbinsX()+1):
0051
0052 roc=0
0053
0054 bin=RocMap.GetBin(i,j)
0055 digi=RocMap.GetBinContent(bin)
0056
0057 if (digi!=0): continue
0058
0059 mod=-4 + int((i-1)/8)
0060
0061 if (lad==0 or mod==0): continue
0062
0063 if (lad < 0 and mod < 0) :
0064
0065 shell="BmO"
0066
0067 if (alad%2==1):
0068 if (j%2==1): roc=(i-1)%8
0069 if (j%2==0): roc= 15-(i-1)%8
0070
0071 elif (alad%2==0):
0072 if (j%2 ==0): roc=(i-1)%8
0073 if (j%2 ==1): roc= 15-(i-1)%8
0074
0075 if (lad > 0 and mod < 0):
0076
0077 shell= "BmI";
0078
0079 if (lad%2==1):
0080 if (j%2 ==0): roc=(i-1)%8
0081 if (j%2 ==1): roc= 15-(i-1)%8
0082
0083 if (lad%2==0):
0084 if (j%2 ==1): roc=(i-1)%8
0085 if (j%2 ==0): roc= 15-(i-1)%8
0086
0087 if (lad > 0 and mod > 0):
0088
0089 shell= "BpI"
0090
0091 if (lad%2==1):
0092 if (j%2 ==1): roc=7-(i-1)%8
0093 if (j%2 ==0): roc=8+(i-1)%8
0094
0095 if (lad%2==0):
0096 if (j%2 ==0): roc=7-(i-1)%8
0097 if (j%2 ==1): roc=8+(i-1)%8
0098
0099 if (lad < 0 and mod > 0):
0100 shell= "BpO"
0101
0102 if (alad%2==1):
0103 if (j%2 ==0): roc=7-(i-1)%8
0104 if (j%2 ==1): roc=8+(i-1)%8
0105
0106 if (alad%2==0):
0107 if (j%2 ==1): roc=7-(i-1)%8
0108 if (j%2 ==0): roc=8+(i-1)%8
0109
0110 f1=open(getFileInPath('DQM/SiStripMonitorClient/data/detids.dat'))
0111 modwritten=False
0112 Mod_check = "LYR"+str(l) + "_LDR" + str(abs(lad)) + "F_MOD" +str(abs(mod))
0113 shell_check = "BPix_" + str(shell)
0114
0115 for line in f1:
0116 refName=line.split(" ")[1]
0117 if modwritten: break
0118 shell_ref = str(refName[:8]).strip()
0119 module_ref = str(refName[14:]).strip()
0120
0121 if (Mod_check == module_ref) and (shell_check == shell_ref):
0122
0123 ModuleName_BPIX = refName.strip()+"_ROC "
0124 BmLYR1_check = ModuleName_BPIX.split('_')
0125 if ((BmLYR1_check[1] == "BmI" or BmLYR1_check[1] == "BmO") and (BmLYR1_check[3] == "LYR1")):
0126 if int(roc) <= 7:
0127 roc = str(int(roc)+8)
0128 elif int(roc) >= 8:
0129 roc =str(int(roc)-8)
0130
0131 BPix_Name = ModuleName_BPIX + str(roc)
0132 BPIXCounter_v0.append(BPix_Name)
0133 BPIXCounter = list(set(BPIXCounter_v0))
0134 modwritten=True
0135 return BPIXCounter
0136
0137
0138
0139
0140
0141 hnameF="digi_occupancy_per_SignedDiskCoord_per_SignedBladePanelCoord_PXRing_"
0142 minbld=[-11,-17]
0143
0144 def FPIX_list(inputFile):
0145 FPIXCounter = []
0146 DQMfile=TFile(inputFile)
0147 for r in range (1,3):
0148
0149 hname=hnameF + str(r)
0150 RocMap=DQMfile.FindObjectAny(hname)
0151
0152 for j in range(1,RocMap.GetNbinsY()+1):
0153
0154 bld=minbld[r-1] + int(j-1)/4
0155 abld=abs(bld)
0156
0157 for i in range(1,RocMap.GetNbinsX()+1):
0158
0159 roc=0
0160
0161 bin=RocMap.GetBin(i,j)
0162 digi=RocMap.GetBinContent(bin)
0163
0164 if (digi!=0): continue
0165
0166 disk=-3 + int(i-1)/8
0167
0168 if (bld==0 or disk==0): continue
0169
0170 pnl=0
0171
0172 if ((j-1)%4==0 or (j-1)%4==1): pnl=2
0173 if ((j-1)%4==2 or (j-1)%4==3): pnl=1
0174
0175 if (disk < 0 and bld <0):
0176
0177 shell= "BmO"
0178
0179 if ((j-1)%4==0 or (j-1)%4==3): roc= 15-(i-1)%8
0180 if ((j-1)%4==1 or (j-1)%4==2): roc= (i-1)%8
0181
0182 if (disk < 0 and bld >0):
0183
0184 shell= "BmI";
0185
0186 if ((j-1)%4==0 or (j-1)%4==3): roc= 15-(i-1)%8
0187 if ((j-1)%4==1 or (j-1)%4==2): roc= (i-1)%8
0188
0189 if (disk > 0 and bld >0):
0190
0191 shell= "BpI"
0192
0193 if ((j-1)%4==0 or (j-1)%4==3): roc=7-(i-1)%8
0194 if ((j-1)%4==1 or (j-1)%4==2): roc=8+(i-1)%8
0195
0196 if (disk > 0 and bld <0):
0197
0198 shell= "BpO"
0199
0200 if ((j-1)%4==0 or (j-1)%4==3): roc=7-(i-1)%8;
0201 if ((j-1)%4==1 or (j-1)%4==2): roc=8+(i-1)%8;
0202
0203
0204 FPix_Name = "FPix_" + str(shell) + "_D" + str(abs(disk)) + '_BLD'+ str(abs(bld)) + '_PNL' + str(abs(pnl)) + '_RNG'+ str(abs(r)) + "_ROC " + str(roc) +""
0205 FPIXCounter.append(FPix_Name)
0206 return FPIXCounter
0207
0208
0209 deadROCList_online = list(set(BPIX_list(filename_online))) + list(set(FPIX_list(filename_online)))
0210 deadROCList_offline = list(set(BPIX_list(filename_offline))) + list(set(FPIX_list(filename_offline)))
0211
0212 MaskedROC_DuringRun = list(set(deadROCList_online) - set(deadROCList_offline))
0213 print('Number of New Dead ROC: '+ str(len(MaskedROC_DuringRun)))
0214
0215 outFileName = 'DeadROC_Diff.txt'
0216 outFileName_online = 'DeadROC_online.txt'
0217 outFileName_offline = 'DeadROC_offline.txt'
0218 f = open(outFileName,"w")
0219 f1 = open(outFileName_online,"w")
0220 f2 = open(outFileName_offline,"w")
0221
0222 for i in range(len(MaskedROC_DuringRun)):
0223
0224 f.write(MaskedROC_DuringRun[i]+"\n")
0225 f.close()
0226
0227 for i in range(len(deadROCList_online)):
0228
0229 f1.write(deadROCList_online[i]+"\n")
0230 f1.close()
0231
0232 for i in range(len(deadROCList_offline)):
0233
0234 f2.write(deadROCList_offline[i]+"\n")
0235 f2.close()
0236
0237
