Merge pull request cms-sw#234 from hyunyong/GEMCosmicMuon

Gem cosmic muon

Validation/GEMCR/data/GEMRAWID.dat

@@ -0,0 +1 @@
+chamber_1_layer_1_roll_1 : 671170562, chamber_1_layer_1_roll_2 : 671236098, chamber_1_layer_1_roll_3 : 671301634, chamber_1_layer_1_roll_4 : 671367170, chamber_1_layer_1_roll_5 : 671432706, chamber_1_layer_1_roll_6 : 671498242, chamber_1_layer_1_roll_7 : 671563778, chamber_1_layer_1_roll_8 : 671629314, chamber_1_layer_2_roll_1 : 671186946, chamber_1_layer_2_roll_2 : 671252482, chamber_1_layer_2_roll_3 : 671318018, chamber_1_layer_2_roll_4 : 671383554, chamber_1_layer_2_roll_5 : 671449090, chamber_1_layer_2_roll_6 : 671514626, chamber_1_layer_2_roll_7 : 671580162, chamber_1_layer_2_roll_8 : 671645698, chamber_3_layer_1_roll_1 : 671171074, chamber_3_layer_1_roll_2 : 671236610, chamber_3_layer_1_roll_3 : 671302146, chamber_3_layer_1_roll_4 : 671367682, chamber_3_layer_1_roll_5 : 671433218, chamber_3_layer_1_roll_6 : 671498754, chamber_3_layer_1_roll_7 : 671564290, chamber_3_layer_1_roll_8 : 671629826, chamber_3_layer_2_roll_1 : 671187458, chamber_3_layer_2_roll_2 : 671252994, chamber_3_layer_2_roll_3 : 671318530, chamber_3_layer_2_roll_4 : 671384066, chamber_3_layer_2_roll_5 : 671449602, chamber_3_layer_2_roll_6 : 671515138, chamber_3_layer_2_roll_7 : 671580674, chamber_3_layer_2_roll_8 : 671646210, chamber_5_layer_1_roll_1 : 671171586, chamber_5_layer_1_roll_2 : 671237122, chamber_5_layer_1_roll_3 : 671302658, chamber_5_layer_1_roll_4 : 671368194, chamber_5_layer_1_roll_5 : 671433730, chamber_5_layer_1_roll_6 : 671499266, chamber_5_layer_1_roll_7 : 671564802, chamber_5_layer_1_roll_8 : 671630338, chamber_5_layer_2_roll_1 : 671187970, chamber_5_layer_2_roll_2 : 671253506, chamber_5_layer_2_roll_3 : 671319042, chamber_5_layer_2_roll_4 : 671384578, chamber_5_layer_2_roll_5 : 671450114, chamber_5_layer_2_roll_6 : 671515650, chamber_5_layer_2_roll_7 : 671581186, chamber_5_layer_2_roll_8 : 671646722, chamber_7_layer_1_roll_1 : 671172098, chamber_7_layer_1_roll_2 : 671237634, chamber_7_layer_1_roll_3 : 671303170, chamber_7_layer_1_roll_4 : 671368706, chamber_7_layer_1_roll_5 : 671434242, chamber_7_layer_1_roll_6 : 671499778, chamber_7_layer_1_roll_7 : 671565314, chamber_7_layer_1_roll_8 : 671630850, chamber_7_layer_2_roll_1 : 671188482, chamber_7_layer_2_roll_2 : 671254018, chamber_7_layer_2_roll_3 : 671319554, chamber_7_layer_2_roll_4 : 671385090, chamber_7_layer_2_roll_5 : 671450626, chamber_7_layer_2_roll_6 : 671516162, chamber_7_layer_2_roll_7 : 671581698, chamber_7_layer_2_roll_8 : 671647234, chamber_9_layer_1_roll_1 : 671172610, chamber_9_layer_1_roll_2 : 671238146, chamber_9_layer_1_roll_3 : 671303682, chamber_9_layer_1_roll_4 : 671369218, chamber_9_layer_1_roll_5 : 671434754, chamber_9_layer_1_roll_6 : 671500290, chamber_9_layer_1_roll_7 : 671565826, chamber_9_layer_1_roll_8 : 671631362, chamber_9_layer_2_roll_1 : 671188994, chamber_9_layer_2_roll_2 : 671254530, chamber_9_layer_2_roll_3 : 671320066, chamber_9_layer_2_roll_4 : 671385602, chamber_9_layer_2_roll_5 : 671451138, chamber_9_layer_2_roll_6 : 671516674, chamber_9_layer_2_roll_7 : 671582210, chamber_9_layer_2_roll_8 : 671647746, chamber_11_layer_1_roll_1 : 671173122, chamber_11_layer_1_roll_2 : 671238658, chamber_11_layer_1_roll_3 : 671304194, chamber_11_layer_1_roll_4 : 671369730, chamber_11_layer_1_roll_5 : 671435266, chamber_11_layer_1_roll_6 : 671500802, chamber_11_layer_1_roll_7 : 671566338, chamber_11_layer_1_roll_8 : 671631874, chamber_11_layer_2_roll_1 : 671189506, chamber_11_layer_2_roll_2 : 671255042, chamber_11_layer_2_roll_3 : 671320578, chamber_11_layer_2_roll_4 : 671386114, chamber_11_layer_2_roll_5 : 671451650, chamber_11_layer_2_roll_6 : 671517186, chamber_11_layer_2_roll_7 : 671582722, chamber_11_layer_2_roll_8 : 671648258, chamber_13_layer_1_roll_1 : 671173634, chamber_13_layer_1_roll_2 : 671239170, chamber_13_layer_1_roll_3 : 671304706, chamber_13_layer_1_roll_4 : 671370242, chamber_13_layer_1_roll_5 : 671435778, chamber_13_layer_1_roll_6 : 671501314, chamber_13_layer_1_roll_7 : 671566850, chamber_13_layer_1_roll_8 : 671632386, chamber_13_layer_2_roll_1 : 671190018, chamber_13_layer_2_roll_2 : 671255554, chamber_13_layer_2_roll_3 : 671321090, chamber_13_layer_2_roll_4 : 671386626, chamber_13_layer_2_roll_5 : 671452162, chamber_13_layer_2_roll_6 : 671517698, chamber_13_layer_2_roll_7 : 671583234, chamber_13_layer_2_roll_8 : 671648770, chamber_15_layer_1_roll_1 : 671174146, chamber_15_layer_1_roll_2 : 671239682, chamber_15_layer_1_roll_3 : 671305218, chamber_15_layer_1_roll_4 : 671370754, chamber_15_layer_1_roll_5 : 671436290, chamber_15_layer_1_roll_6 : 671501826, chamber_15_layer_1_roll_7 : 671567362, chamber_15_layer_1_roll_8 : 671632898, chamber_15_layer_2_roll_1 : 671190530, chamber_15_layer_2_roll_2 : 671256066, chamber_15_layer_2_roll_3 : 671321602, chamber_15_layer_2_roll_4 : 671387138, chamber_15_layer_2_roll_5 : 671452674, chamber_15_layer_2_roll_6 : 671518210, chamber_15_layer_2_roll_7 : 671583746, chamber_15_layer_2_roll_8 : 671649282, chamber_17_layer_1_roll_1 : 671174658, chamber_17_layer_1_roll_2 : 671240194, chamber_17_layer_1_roll_3 : 671305730, chamber_17_layer_1_roll_4 : 671371266, chamber_17_layer_1_roll_5 : 671436802, chamber_17_layer_1_roll_6 : 671502338, chamber_17_layer_1_roll_7 : 671567874, chamber_17_layer_1_roll_8 : 671633410, chamber_17_layer_2_roll_1 : 671191042, chamber_17_layer_2_roll_2 : 671256578, chamber_17_layer_2_roll_3 : 671322114, chamber_17_layer_2_roll_4 : 671387650, chamber_17_layer_2_roll_5 : 671453186, chamber_17_layer_2_roll_6 : 671518722, chamber_17_layer_2_roll_7 : 671584258, chamber_17_layer_2_roll_8 : 671649794, chamber_19_layer_1_roll_1 : 671175170, chamber_19_layer_1_roll_2 : 671240706, chamber_19_layer_1_roll_3 : 671306242, chamber_19_layer_1_roll_4 : 671371778, chamber_19_layer_1_roll_5 : 671437314, chamber_19_layer_1_roll_6 : 671502850, chamber_19_layer_1_roll_7 : 671568386, chamber_19_layer_1_roll_8 : 671633922, chamber_19_layer_2_roll_1 : 671191554, chamber_19_layer_2_roll_2 : 671257090, chamber_19_layer_2_roll_3 : 671322626, chamber_19_layer_2_roll_4 : 671388162, chamber_19_layer_2_roll_5 : 671453698, chamber_19_layer_2_roll_6 : 671519234, chamber_19_layer_2_roll_7 : 671584770, chamber_19_layer_2_roll_8 : 671650306, chamber_21_layer_1_roll_1 : 671175682, chamber_21_layer_1_roll_2 : 671241218, chamber_21_layer_1_roll_3 : 671306754, chamber_21_layer_1_roll_4 : 671372290, chamber_21_layer_1_roll_5 : 671437826, chamber_21_layer_1_roll_6 : 671503362, chamber_21_layer_1_roll_7 : 671568898, chamber_21_layer_1_roll_8 : 671634434, chamber_21_layer_2_roll_1 : 671192066, chamber_21_layer_2_roll_2 : 671257602, chamber_21_layer_2_roll_3 : 671323138, chamber_21_layer_2_roll_4 : 671388674, chamber_21_layer_2_roll_5 : 671454210, chamber_21_layer_2_roll_6 : 671519746, chamber_21_layer_2_roll_7 : 671585282, chamber_21_layer_2_roll_8 : 671650818, chamber_23_layer_1_roll_1 : 671176194, chamber_23_layer_1_roll_2 : 671241730, chamber_23_layer_1_roll_3 : 671307266, chamber_23_layer_1_roll_4 : 671372802, chamber_23_layer_1_roll_5 : 671438338, chamber_23_layer_1_roll_6 : 671503874, chamber_23_layer_1_roll_7 : 671569410, chamber_23_layer_1_roll_8 : 671634946, chamber_23_layer_2_roll_1 : 671192578, chamber_23_layer_2_roll_2 : 671258114, chamber_23_layer_2_roll_3 : 671323650, chamber_23_layer_2_roll_4 : 671389186, chamber_23_layer_2_roll_5 : 671454722, chamber_23_layer_2_roll_6 : 671520258, chamber_23_layer_2_roll_7 : 671585794, chamber_23_layer_2_roll_8 : 671651330, chamber_25_layer_1_roll_1 : 671176706, chamber_25_layer_1_roll_2 : 671242242, chamber_25_layer_1_roll_3 : 671307778, chamber_25_layer_1_roll_4 : 671373314, chamber_25_layer_1_roll_5 : 671438850, chamber_25_layer_1_roll_6 : 671504386, chamber_25_layer_1_roll_7 : 671569922, chamber_25_layer_1_roll_8 : 671635458, chamber_25_layer_2_roll_1 : 671193090, chamber_25_layer_2_roll_2 : 671258626, chamber_25_layer_2_roll_3 : 671324162, chamber_25_layer_2_roll_4 : 671389698, chamber_25_layer_2_roll_5 : 671455234, chamber_25_layer_2_roll_6 : 671520770, chamber_25_layer_2_roll_7 : 671586306, chamber_25_layer_2_roll_8 : 671651842, chamber_27_layer_1_roll_1 : 671177218, chamber_27_layer_1_roll_2 : 671242754, chamber_27_layer_1_roll_3 : 671308290, chamber_27_layer_1_roll_4 : 671373826, chamber_27_layer_1_roll_5 : 671439362, chamber_27_layer_1_roll_6 : 671504898, chamber_27_layer_1_roll_7 : 671570434, chamber_27_layer_1_roll_8 : 671635970, chamber_27_layer_2_roll_1 : 671193602, chamber_27_layer_2_roll_2 : 671259138, chamber_27_layer_2_roll_3 : 671324674, chamber_27_layer_2_roll_4 : 671390210, chamber_27_layer_2_roll_5 : 671455746, chamber_27_layer_2_roll_6 : 671521282, chamber_27_layer_2_roll_7 : 671586818, chamber_27_layer_2_roll_8 : 671652354, chamber_29_layer_1_roll_1 : 671177730, chamber_29_layer_1_roll_2 : 671243266, chamber_29_layer_1_roll_3 : 671308802, chamber_29_layer_1_roll_4 : 671374338, chamber_29_layer_1_roll_5 : 671439874, chamber_29_layer_1_roll_6 : 671505410, chamber_29_layer_1_roll_7 : 671570946, chamber_29_layer_1_roll_8 : 671636482, chamber_29_layer_2_roll_1 : 671194114, chamber_29_layer_2_roll_2 : 671259650, chamber_29_layer_2_roll_3 : 671325186, chamber_29_layer_2_roll_4 : 671390722, chamber_29_layer_2_roll_5 : 671456258, chamber_29_layer_2_roll_6 : 671521794, chamber_29_layer_2_roll_7 : 671587330, chamber_29_layer_2_roll_8 : 671652866

Validation/GEMCR/test/MakePlot.py

@@ -30,7 +30,7 @@
 oDir = "Run%06d_Plots/"%run
 
 
-maskPlot = True
+maskPlot = runConfig.runWithMasking
 
 
 SLOTLIST=[]
@@ -58,7 +58,7 @@
     chi+=1
 
 import csv
-tmpID = [x for x in csv.reader(open("GEMRAWID.dat","r"))][0]
+tmpID = [x for x in csv.reader(open("../data/GEMRAWID.dat","r"))][0]
 GEMRAWID = {}
 GEMNAME = {}
 for x in tmpID:
@@ -113,8 +113,8 @@ def writeMask(hist):
       if v>cut*2.0: hotStrip.write("%d %d\n"%(findRawId(name), x))
 
 if maskPlot:
-  tmpM = open("GEMMaskVec.dat", "r")
-  tmpH = open("GEMHotVec.dat", "r")
+  tmpM = open("../data/GEMMaskVecRun%6d.dat"%run, "r")
+  tmpH = open("../data/GEMHotVecRun%6d.dat"%run, "r")
   maskL = {}
   for x in tmpM:
     tmpL = x.split()
@@ -705,8 +705,8 @@ def draw_plot( file, tDir,oDir ) :
 
   os.system("mkdir -p "+oDir )
   outRoot = TFile(oDir+oDir[:-1]+".root", "recreate")
-  mask = open("GEMMaskVec.dat", "w")
-  hotStrip = open("GEMHotVec.dat", "w")
+  mask = open("../data/GEMMaskVecRun%06d.dat"%run, "w")
+  hotStrip = open("../data/GEMHotVecRun%06d.dat"%run, "w")
   for c in chamber:
     outRoot.mkdir(c.replace("/",""))
     os.system("mkdir -p "+oDir+"/"+c.replace("/",""))  

Validation/GEMCR/test/configureRun_cfi.py

@@ -6,9 +6,22 @@
 #RAWFileName="/afs/cern.ch/work/d/dorney/CMS_GEM/Data/QC8/run000036_Test_TIF_2016-11-28.dat"
 #RAWFileName="/afs/cern.ch/work/d/dorney/CMS_GEM/Data/QC8/run000037_Test_TIF_2016-11-28.dat"
 RAWFileName="run000044_Cosmics_TIF_2016-12-03.dat"
-
+RAWFileName="/afs/cern.ch/user/h/hyunyong/public/run000080_Cosmics_TIF_2016-12-05.dat"
 RunNumber=int(RAWFileName.split("/")[-1].split("_")[0][3:])
 
+GEMMask = "Validation/GEMCR/data/GEMMaskVecRun%06d.dat"%RunNumber
+GEMHot = "Validation/GEMCR/data/GEMHotVecRun%06d.dat"%RunNumber
+
+if os.path.isfile(os.environ['SRT_CMSSW_BASE_SCRAMRTDEL']+"/src/"+GEMHot): runWithMasking = True
+else : runWithMasking = False
+
+if not runWithMasking :
+  GEMMask = "Validation/GEMCR/data/GEMMaskEmp.dat" 
+  GEMHot = "Validation/GEMCR/data/GEMMaskEmp.dat"
+
+if runWithMasking : print "This run is running with masking"
+else : print  "This run dose not have mask lists" 
+
 makeTrack = False
 if makeTrack : OutputFileName='Reco_Run%06d.root'%RunNumber
 else : OutputFileName='Reco_Run%06d_noTr.root'%RunNumber
@@ -22,7 +35,7 @@
 
 #MaxEvents=-1
 #MaxEvents=92000
-MaxEvents=9600
+MaxEvents=1000
 
 #makeTrack = False
 

Validation/GEMCR/test/runGEMCosmicStand_data.py

@@ -109,7 +109,8 @@
 process.MuonServiceProxy.ServiceParameters.Propagators.append('StraightLinePropagator')
 
 process.gemRecHits.gemDigiLabel = cms.InputTag("muonGEMDigis","","RECO")
-
+process.gemRecHits.maskvecfile = cms.FileInPath(runConfig.GEMMask)
+process.gemRecHits.deadvecfile = cms.FileInPath(runConfig.GEMHot)
 process.GEMCosmicMuon = cms.EDProducer("GEMCosmicMuon",
                                        process.MuonServiceProxy,
                                        gemRecHitLabel = cms.InputTag("gemRecHits"),