(r14) SiPixelQuality < CMS

CMS Web>NewPhysics>CmsSusy>SiPixelQuality (revision 14) (raw view)~~EditAttach~~

<style type="text/css" media="all"> pre { text-align: left; padding: 10px;margin-left: 20px; color: black; } pre.command {background-color: lightgrey;} pre.cfg {background-color: lightblue;} pre.code {background-color: lightpink;} pre.output {background-color: lightgreen;} pre.tofix {background-color: thistle;} </style>

%TOC{}%

---++ Tracking effect

What happens with modules marked bad that come alive from time to time?

Studying the valid/missing/inactive modules in three scenarios

---+++ Scenario 1)
---++++ We have a new bad component but it is not included in the DB yet <br/>
--> in SiPixelQuality_v34_data we had a new bad half module (BmO_Sec8_Lyr1_Ldr9F_Mod4 (TBM-B) -- 302059012) <br/>
   * Run range 275891 (the moduel is bad from this run) - 275914 (hlt tag was accepted and uploaded, in the meantime there was a half module that was bad but the DB did not contain it)
<pre class="command">
TChain chain("trajTree");
chain.Add("nTuplePilotBlade_All*.root")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.type==0&&traj.onEdge==0","COLZ") // validhit
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.type==1&&traj.onEdge==0","COLZ") // missing
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.type==2&&traj.onEdge==0","COLZ") // inactive

chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302059012&&traj.type==0&&traj.onEdge==0","COLZ") // valid pixels

TimingStudy
TChain chain("trajTree");
chain.Add("Ntuple*.root")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.validhit==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.missing==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.inactive==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302059012&&traj.badhit==1&&traj.onedge==0","COLZ")

chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302059012&&traj.validhit==1&&traj.onedge==0","COLZ")

TChain chain("clustTree");
chain.Add("Ntuple*.root")
chain.Draw("lx:ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302059012","COLZ") // all clusters
</pre>

---++++ Result
   * half of the module does not contain valid hits
   * there are no missing hits
   * there are no inactive hits

---+++ Scenario 2)
---++++ We have a bad component and it is included in the DB already <br/>
--> in SiPixelQuality_v34_data an old half module (BPix_BmO_SEC7_LYR2_LDR13F_MOD3 (ROC0-7)  -- 302126344 [this was added in v33])
   * Run range 275891-275914 (same range is OK)
<pre class="command">
TChain chain("trajTree");
chain.Add("nTuplePilotBlade_All*.root")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.type==0&&traj.onEdge==0","COLZ") // validhit
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.type==1&&traj.onEdge==0","COLZ") // missing
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.type==2&&traj.onEdge==0","COLZ") // inactive

chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302126344&&traj.type==0&&traj.onEdge==0","COLZ") // valid pixels
chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302126344&&traj.type==1&&traj.onEdge==0","COLZ") // missing pixels

TimingStudy
TChain chain("trajTree");
chain.Add("Ntuple*.root")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.validhit==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.missing==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.inactive==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302126344&&traj.badhit==1&&traj.onedge==0","COLZ")

chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302126344&&traj.validhit==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302126344&&traj.missing==1&&traj.onedge==0","COLZ")

TChain chain("clustTree");
chain.Add("Ntuple*.root")
chain.Draw("lx:ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302126344","COLZ") // all clusters

</pre>

---++++ Result
   * half of the module is does not contain the valid hits
   * the other half  contain the missing hits
   * there are no inactive hits

---+++ Scenario 3)
---++++ A module is marked bad although it is working <br/>
--> v35_data (302123800)<br/>
   * Run range 276460-276870 (after this the mistake was corrected)<br/>

<pre class="command">
TChain chain("trajTree");
chain.Add("nTuplePilotBlade_All*.root")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.type==0&&traj.onEdge==0","COLZ") // validhit
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.type==1&&traj.onEdge==0","COLZ") // missing
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.type==2&&traj.onEdge==0","COLZ") // inactive

TimingStudy
TChain chain("trajTree");
chain.Add("Ntuple*.root")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.validhit==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.missing==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.inactive==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(8,-3.0,3.0,2,-.55,.55)","module.rawid==302123800&&traj.badhit==1&&traj.onedge==0","COLZ")

chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302123800&&traj.validhit==1&&traj.onedge==0","COLZ")
chain.Draw("traj.lx:traj.ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302123800&&traj.missing==1&&traj.onedge==0","COLZ")

TChain chain("clustTree");
chain.Add("Ntuple*.root")
chain.Draw("lx:ly>>h1(416,-3.0,3.0,160,-.55,.55)","module.rawid==302123800","COLZ") // all clusters
</pre>


---++++ Result
   * the whole module contains valid hits, the lower part has twice of the number of hits as the upper half
   * there are missing hits in the upper half 
   * there are no inactive hits

---++ Dateset boarders

|  Dataset name  |  From #Run  |  Until #Run  |  Comment  |
|  Run2016B-v2   |  273150  |  275376  | /ZeroBias/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias1/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias2/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias3/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias4/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias5/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias6/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias7/Run2016B-v2/RAW |
|  Run2016B-v2   |  273511  |  274100  | /ZeroBias8/Run2016B-v2/RAW |
|  Run2016B-v1   |  272023  |  273146  |  /ZeroBias/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias1/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias2/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias3/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias4/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias5/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias6/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias7/Run2016B-v1/RAW  |
|  Run2016B-v1   |  272007  |  273138  |  /ZeroBias8/Run2016B-v1/RAW  |
|  Run2016B-v1   |  271884  |  273138  |  /ZeroBias9/Run2016B-v1/RAW  |
---++ Recipe
---+++ CMSSW 80X
---++++ Enviroment
<pre class="command">
cmsrel CMSSW_8_0_2
cd CMSSW_8_0_2/src
cmsenv

git cms-addpkg CondTools/SiPixel
git clone https://github.com/jkarancs/PixelTimingStudy DPGAnalysis/PixelTimingStudy
vim DPGAnalysis/PixelTimingStudy/plugins/TimingStudy.cc 
#define COMPLETE 0
#define SPLIT 2
scram b -j
</pre>

---++++ New Quality
<pre class="command">
cd CondTools/SiPixel/test
wget http://www.grid.kfki.hu/twiki/pub/CMS/SiPixelQuality/0SiPixelBadModuleByHandBuilder_cfg_v29.py.txt 0SiPixelBadModuleByHandBuilder_cfg_v29.py
vim 0SiPixelBadModuleByHandBuilder_cfg_*.py
# edit accordingly
</pre>

---++++ Validation using MC
<pre class="command">
cd DPGAnalysis/PixelTimingStudy/test/
wget http://www.grid.kfki.hu/twiki/pub/CMS/SiPixelQuality/2TimingStudy_RawData_CMSSW802_cfg.py.txt 2TimingStudy_RawData_CMSSW802_cfg.py
cmsrun 2TimingStudy_RawData_CMSSW802_cfg.py
# make sure this is True: process.siPixelDigis.UseQualityInfo = cms.bool(True)
# creates Ntuple_IOV30_5k.root

cp Ntuple_IOV30_5k.root  ../../3Ploting@/MCPlot/.
cd  ../../3Ploting@/MCPlot/
# edit ROC_Occupancy.C accordingly
root -l 
.L ROC_Occupancy.C+g
.X ROC_Occupancy.C
</pre>

---++++ Occupancy map from Data
<pre class="command">
cd DPGAnalysis/PixelTimingStudy/test/
wget http://www.grid.kfki.hu/twiki/pub/CMS/SiPixelQuality/3TimingStudy_RawData_CMSSW802_cfg.py.txt 3TimingStudy_RawData_CMSSW802_cfg.py
cmsrun 3TimingStudy_RawData_CMSSW802_cfg.py
#creates Ntuple_IOV30_5k.root

cp Ntuple_IOV30_5k.root  ../../3Ploting@/DataPlot/.
cd  ../../3Ploting@/DataPlot/
# edit ROC_Occupancy.C accordingly
root -l 
.L ROC_Occupancy.C+g
.X ROC_Occupancy.C
</pre>

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
txt	0SiPixelBadModuleByHandBuilder_cfg_v29.py.txt	r1	manage	5.2 K	2016-03-22 - 10:28	TamasVami	SiPixelQuality_byHand_cfg_v29
txt	0SiPixelBadModuleByHandBuilder_cfg_v33.py.txt	r1	manage	3.6 K	2016-05-18 - 16:53	TamasVami	0SiPixelBadModuleByHandBuilder_cfg_v33
txt	0SiPixelBadModuleByHandBuilder_cfg_v36.py.txt	r1	manage	3.9 K	2016-07-13 - 04:07	TamasVami	0SiPixelBadModuleByHandBuilder_cfg_v36
txt	2TimingStudy_RawData_CMSSW802_cfg.py.txt	r1	manage	9.5 K	2016-05-09 - 13:33	TamasVami	2TimingStudy_RawData_CMSSW802_cfg
png	OneCoolingLoop_Layer1.png	r1	manage	11.7 K	2016-07-13 - 03:16	TamasVami	OneCoolingLoop _Layer1
png	OneCoolingLoop_Layer2.png	r1	manage	11.5 K	2016-07-13 - 03:17	TamasVami	OneCoolingLoop _Layer2
png	OneCoolingLoop_Layer3.png	r1	manage	46.6 K	2016-07-13 - 03:17	TamasVami	OneCoolingLoop _Layer3
png	OneCoolingLoop_Layer4.png	r1	manage	27.0 K	2016-07-13 - 03:17	TamasVami	OneCoolingLoop _Layer4
c	ROC_Occupancy.C	r1	manage	17.0 K	2016-07-22 - 19:24	TamasVami	ROC_Occupancy -- 2016-07-22

Topic revision: r14 - 2016-07-30 - TamasVami

CMS (CERN)

Copyright &© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback