|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
|
#! /usr/bin/env python3
import ROOT
import optparse
import re
import os
def main():
parser = optparse.OptionParser()
(options, args) = parser.parse_args()
ROOT.gROOT.SetBatch(True)
for filename in args:
if "vDrift_segment_" in filename:
#name="DQM"
m=re.search("vDrift_segment_(\d*)",filename)
run=m.group(1)
#path="DQMData/Run "+run+"/DT/Run summary/DTCalibValidation"
f = open('dtVDriftAnalyzer_cfg.py','w')
print("from CalibMuon.DTCalibration.Workflow.addPoolDBESSource import addPoolDBESSource", file=f)
print("from CalibMuon.DTCalibration.dtVDriftAnalyzer_cfg import process", file=f)
print("addPoolDBESSource(process = process, moduleName = 'vDriftDB',record = 'DTMtimeRcd',tag = 'vDrift', connect = 'sqlite_file:"+filename+"')", file=f)
print("process.dtVDriftAnalyzer.rootFileName = 'dtVDriftAnalyzer_dtVDriftCalibration"+run+".root'", file=f)
f.close()
os.system("cmsRun dtVDriftAnalyzer_cfg.py")
name="vdrift"
runvdrift(name, run, "dtVDriftAnalyzer_dtVDriftCalibration"+run+".root")
if "vDrift_meantimer_" in filename:
#name="DQM"
m=re.search("vDrift_meantimer_(\d*)",filename)
run=m.group(1)
#path="DQMData/Run "+run+"/DT/Run summary/DTCalibValidation"
f = open('dtVDriftAnalyzer_cfg.py','w')
print("from CalibMuon.DTCalibration.Workflow.addPoolDBESSource import addPoolDBESSource", file=f)
print("from CalibMuon.DTCalibration.dtVDriftAnalyzer_cfg import process", file=f)
print("addPoolDBESSource(process = process, moduleName = 'vDriftDB',record = 'DTMtimeRcd',tag = 'vDrift', connect = 'sqlite_file:"+filename+"')", file=f)
print("process.dtVDriftAnalyzer.rootFileName = 'dtVDriftAnalyzer_dtVDriftCalibration"+run+".root'", file=f)
f.close()
os.system("cmsRun dtVDriftAnalyzer_cfg.py")
name="vdrift"
runvdrift(name, run, "dtVDriftAnalyzer_dtVDriftCalibration"+run+".root")
elif "DQM" in filename:
name="DQM"
m=re.search("R000(\d*)__",filename)
run=m.group(1)
path="DQMData/Run "+run+"/DT/Run summary/DTCalibValidation"
runttrig(name, filename, path, run)
elif "TestPulses" in filename:
name="T0"
runt0(name,filename,123456)
else:
name="DTRV"
m=re.search("DTResidualValidation_(\d*)\.root",filename)
run=m.group(1)
path="DTResiduals"
runttrig(name, filename, path, run)
def runttrig(name, filename, path, run):
from CalibMuon.DTCalibration.PlottingTools.plotResiduals import *
for SL in [1,2,3]:
mean,sigma = plot(filename, SL,dir=path)
mean[0].Print(name+run+"-SL"+str(SL)+"-mean.pdf")
mean[0].SaveAs(name+run+"-SL"+str(SL)+"-mean.root")
sigma[0].Print(name+run+"-SL"+str(SL)+"-sigma.pdf")
sigma[0].SaveAs(name+run+"-SL"+str(SL)+"-sigma.root")
def runvdrift(name, run, filename):
from CalibMuon.DTCalibration.PlottingTools.plotVDriftFromHistos import *
for SL in [1,2,3]:
mean = plot(filename, SL)
mean[0].Print(name+run+"-SL"+str(SL)+"-mean.pdf")
mean[0].SaveAs(name+run+"-SL"+str(SL)+"-mean.root")
#sigma[0].Print(name+run+"-SL"+str(SL)+"-sigma.pdf")
#plot("DTResidualValidation_210614.root", 1,dir="DTResiduals")
#plot("DQM_V0001_R000210634__StreamExpress__HIRun2013-DtCalib-Express-v1-dtTtrigCalibrationFromResiduals-NEW-usedb-rev1__ALCARECO.root", 1,dir="DQMData/Run 210634/DT/Run summary/DTCalibValidation")
def runt0(name,filename,run):
from CalibMuon.DTCalibration.PlottingTools.plotT0FromHistos import *
for SL in [1,2,3]:
mean = plot(filename, SL ,run)
mean[0].Print(name+run+"-SL"+str(SL)+"-mean.pdf")
mean[0].SaveAs(name+run+"-SL"+str(SL)+"-mean.root")
if __name__=="__main__":
main()
|