Skip to content

Commit

Permalink
FixZeroOccupancyOnline
Browse files Browse the repository at this point in the history
  • Loading branch information
@fioriNTU
fioriNTU committed Oct 9, 2016
1 parent d1a8ec9 commit 22d450c
Show file tree
Hide file tree
Showing 3 changed files with 185 additions and 29 deletions.
38 changes: 37 additions & 1 deletion DQM/SiPixelMonitorDigi/interface/SiPixelDigiSource.h
View file
Expand Up @@ -62,6 +62,8 @@
virtual void buildStructure(edm::EventSetup const&);
virtual void bookMEs(DQMStore::IBooker &, const edm::EventSetup& iSetup);

virtual void CountZeroROCsInSubstructure (bool , bool&, SiPixelDigiModule *);

std::string topFolderName_;

private:
Expand Down Expand Up @@ -186,10 +188,44 @@
std::vector<MonitorElement*> meNDigisCHANEndcapDps_;
std::vector<MonitorElement*> meNDigisCHANEndcapDms_;

int NzeroROCs[2];
int NloEffROCs[2];

bool ROCMapToReset;
//the following long list of bools is to patch the ZeroOccupancy ROC filling in a way that a substructure (like BPix/BmO/Layer1) is counted only once as it should be (in the past for each module in the substructure the same number of ZeroOccupancy rocs was added)

bool DoZeroRocsBMO1;
bool DoZeroRocsBMO2;
bool DoZeroRocsBMO3;

bool DoZeroRocsBMI1;
bool DoZeroRocsBMI2;
bool DoZeroRocsBMI3;

bool DoZeroRocsBPO1;
bool DoZeroRocsBPO2;
bool DoZeroRocsBPO3;

bool DoZeroRocsBPI1;
bool DoZeroRocsBPI2;
bool DoZeroRocsBPI3;

bool DoZeroRocsFPO1;
bool DoZeroRocsFPO2;

bool DoZeroRocsFMO1;
bool DoZeroRocsFMO2;

bool DoZeroRocsFPI1;
bool DoZeroRocsFPI2;

bool DoZeroRocsFMI1;
bool DoZeroRocsFMI2;

int bigEventSize;
bool isUpgrade;
bool firstRun;

std::string I_name[1856];
unsigned int I_detId[1856];
int I_fedId[1856];
Expand Down
10 changes: 5 additions & 5 deletions DQM/SiPixelMonitorDigi/src/SiPixelDigiModule.cc
View file
Expand Up @@ -543,8 +543,8 @@ std::pair<int,int> SiPixelDigiModule::getZeroLoEffROCs(){
float SF = 1.0;
if (mePixRocsDisk_ && meZeroOccRocsDisk_){
if (mePixRocsDisk_->getEntries() > 0) SF = float(mePixRocsDisk_->getNbinsX()*mePixRocsDisk_->getNbinsY()/mePixRocsDisk_->getEntries());
for (int i = 1; i < mePixRocsDisk_->getNbinsX(); ++i){
for (int j = 1; j < mePixRocsDisk_->getNbinsY(); ++j){
for (int i = 1; i < mePixRocsDisk_->getNbinsX()+1; ++i){
for (int j = 1; j < mePixRocsDisk_->getNbinsY()+1; ++j){
float localX = float(i) - 0.5;
float localY = float(j)/2.0 + 0.75;
if (mePixRocsDisk_->getBinContent(i,j) < 1 ) {nZeroROC++; meZeroOccRocsDisk_->Fill(localX,localY);}
Expand All @@ -555,11 +555,11 @@ std::pair<int,int> SiPixelDigiModule::getZeroLoEffROCs(){
}
if (mePixRocsLay_ && meZeroOccRocsLay_){
if (mePixRocsLay_->getEntries() > 0) SF = float(mePixRocsLay_->getNbinsX()*mePixRocsLay_->getNbinsY()/mePixRocsLay_->getEntries());
for (int i = 1; i < mePixRocsLay_->getNbinsX(); ++i){
for (int j = 1; j < mePixRocsLay_->getNbinsY(); ++j){
for (int i = 1; i < mePixRocsLay_->getNbinsX()+1; ++i){
for (int j = 1; j < mePixRocsLay_->getNbinsY()+1; ++j){
float localX = float(i) - 0.5;
float localY = float(j)/2.0 + 1.25;
if (mePixRocsLay_->getBinContent(i,j) < 1 ) {nZeroROC++; meZeroOccRocsLay_->Fill(localX,localY);}
if (mePixRocsLay_->getBinContent(i,j) < 8 ) {nZeroROC++; meZeroOccRocsLay_->Fill(localX,localY);} //in some regions of pixel there are modules with no HV but enabled ROCs that sometime give a fake hit, so the dead rocs have to be counted to have less than 8 hits in 10 LS
if (mePixRocsLay_->getBinContent(i,j)*SF < 0.25){nLoEffROC++;}
}
}
Expand Down
166 changes: 143 additions & 23 deletions DQM/SiPixelMonitorDigi/src/SiPixelDigiSource.cc
View file
Expand Up @@ -101,8 +101,10 @@ SiPixelDigiSource::~SiPixelDigiSource()
void
SiPixelDigiSource::beginLuminosityBlock(const edm::LuminosityBlock& lb, edm::EventSetup const&)
{
int thisls = lb.id().luminosityBlock();


int thisls = lb.id().luminosityBlock();

if(modOn && thisls % 10 == 0 && averageDigiOccupancy){
nBPIXDigis = 0;
nFPIXDigis = 0;
Expand All @@ -113,7 +115,47 @@ SiPixelDigiSource::beginLuminosityBlock(const edm::LuminosityBlock& lb, edm::Eve
nFPIXDigis = 0;
for(int i=0; i!=40; i++) nDigisPerFed[i]=0;
}


if (modOn && thisls % 10 == 0) {

ROCMapToReset=true; //the ROC map is reset each 10 lumisections

for (int i=0; i<2;i++) NzeroROCs[i]=0;
for (int i=0; i<2;i++) NloEffROCs[i]=0; //resetting also Zero and low eff. ROC counters

NzeroROCs[1]=-672;
NloEffROCs[1]=-672;//this magic number derives by the way the endcap occupancy is filled, there are always 672 empty bins by construction

//these bools are needed to count zero occupancy plots in the substructure only once each 10 LS
DoZeroRocsBMO1=true;
DoZeroRocsBMO2=true;
DoZeroRocsBMO3=true;

DoZeroRocsBMI1=true;
DoZeroRocsBMI2=true;
DoZeroRocsBMI3=true;

DoZeroRocsBPO1=true;
DoZeroRocsBPO2=true;
DoZeroRocsBPO3=true;

DoZeroRocsBPI1=true;
DoZeroRocsBPI2=true;
DoZeroRocsBPI3=true;

DoZeroRocsFPO1=true;
DoZeroRocsFPO2=true;

DoZeroRocsFMO1=true;
DoZeroRocsFMO2=true;

DoZeroRocsFPI1=true;
DoZeroRocsFPI2=true;

DoZeroRocsFMI1=true;
DoZeroRocsFMI2=true;
}

}

void
Expand Down Expand Up @@ -155,7 +197,6 @@ void SiPixelDigiSource::dqmBeginRun(const edm::Run& r, const edm::EventSetup& iS
<< layOn << "/" << phiOn << std::endl;
LogInfo ("PixelDQM") << "Blade/Disk/Ring" << bladeOn << "/" << diskOn << "/"
<< ringOn << std::endl;

LogInfo ("PixelDQM") << "2DIM IS " << twoDimOn << " and set to high resolution? " << hiRes << "\n";

if(firstRun){
Expand Down Expand Up @@ -213,6 +254,39 @@ void SiPixelDigiSource::dqmBeginRun(const edm::Run& r, const edm::EventSetup& iS
nL4M3 = 0;
nL4M4 = 0;

ROCMapToReset=false;

DoZeroRocsBMO1=false;
DoZeroRocsBMO2=false;
DoZeroRocsBMO3=false;

DoZeroRocsBMI1=false;
DoZeroRocsBMI2=false;
DoZeroRocsBMI3=false;

DoZeroRocsBPO1=false;
DoZeroRocsBPO2=false;
DoZeroRocsBPO3=false;

DoZeroRocsBPI1=false;
DoZeroRocsBPI2=false;
DoZeroRocsBPI3=false;

DoZeroRocsFPO1=false;
DoZeroRocsFPO2=false;

DoZeroRocsFMO1=false;
DoZeroRocsFMO2=false;

DoZeroRocsFPI1=false;
DoZeroRocsFPI2=false;

DoZeroRocsFMI1=false;
DoZeroRocsFMI2=false;

for (int i=0; i<2; i++) NzeroROCs[i]=0;
for (int i=0; i<2;i++) NloEffROCs[i]=0;

// Build map
buildStructure(iSetup);
// Book Monitoring Elements
Expand All @@ -229,6 +303,7 @@ void SiPixelDigiSource::bookHistograms(DQMStore::IBooker & iBooker, edm::Run con
//------------------------------------------------------------------
void SiPixelDigiSource::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{

edm::ESHandle<TrackerTopology> tTopoHandle;
iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
const TrackerTopology *pTT = tTopoHandle.product();
Expand All @@ -237,7 +312,6 @@ void SiPixelDigiSource::analyze(const edm::Event& iEvent, const edm::EventSetup&
edm::Handle< edm::DetSetVector<PixelDigi> > input;
iEvent.getByToken(srcToken_, input);
if (!input.isValid()) return;


int bx = iEvent.bunchCrossing();

Expand All @@ -248,20 +322,20 @@ void SiPixelDigiSource::analyze(const edm::Event& iEvent, const edm::EventSetup&
for(int i=0; i!=192; i++) numberOfDigis[i]=0;
for(int i=0; i!=1152; i++) nDigisPerChan[i]=0;
for(int i=0; i!=4; i++) nDigisPerDisk[i]=0;
int NzeroROCs[2] = {0,-672};
int NloEffROCs[2] = {0,-672};

for (struct_iter = thePixelStructure.begin() ; struct_iter != thePixelStructure.end() ; struct_iter++) {

int numberOfDigisMod = (*struct_iter).second->fill(*input, iSetup,
meNDigisCOMBBarrel_, meNDigisCHANBarrel_,meNDigisCHANBarrelLs_,meNDigisCOMBEndcap_,
modOn, ladOn, layOn, phiOn,
bladeOn, diskOn, ringOn,
twoDimOn, reducedSet, twoDimModOn, twoDimOnlyLayDisk,
nDigisA, nDigisB, isUpgrade);

if (modOn && twoDimOnlyLayDisk && lumiSection%10 == 0) (*struct_iter).second->resetRocMap();



bool barrel = DetId((*struct_iter).first).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
bool endcap = DetId((*struct_iter).first).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);

if(numberOfDigisMod>0){
nEventDigis = nEventDigis + numberOfDigisMod;
nActiveModules++;
Expand All @@ -275,6 +349,28 @@ void SiPixelDigiSource::analyze(const edm::Event& iEvent, const edm::EventSetup&
nTOTmodules=1440;
}
if(barrel){ // Barrel
int layer=PixelBarrelName(DetId((*struct_iter).first),pTT,isUpgrade).layerName();
PixelBarrelName::Shell shell=PixelBarrelName(DetId((*struct_iter).first),pTT,isUpgrade).shell();

//Count Zero Occ Rocs in Barrel in the first event after each 10 Ls
if (ROCMapToReset && shell==PixelBarrelName::mO && twoDimOnlyLayDisk){
if (DoZeroRocsBMO1 && layer==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsBMO1,(*struct_iter).second);
if (DoZeroRocsBMO2 && layer==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsBMO2,(*struct_iter).second);
if (DoZeroRocsBMO3 && layer==3) CountZeroROCsInSubstructure(barrel, DoZeroRocsBMO3,(*struct_iter).second);
} else if (ROCMapToReset && shell==PixelBarrelName::mI && twoDimOnlyLayDisk){
if (DoZeroRocsBMI1 && layer==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsBMI1,(*struct_iter).second);
if (DoZeroRocsBMI2 && layer==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsBMI2,(*struct_iter).second);
if (DoZeroRocsBMI3 && layer==3) CountZeroROCsInSubstructure(barrel, DoZeroRocsBMI3,(*struct_iter).second);
} else if (ROCMapToReset && shell==PixelBarrelName::pO && twoDimOnlyLayDisk){
if (DoZeroRocsBPO1 && layer==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsBPO1,(*struct_iter).second);
if (DoZeroRocsBPO2 && layer==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsBPO2,(*struct_iter).second);
if (DoZeroRocsBPO3 && layer==3) CountZeroROCsInSubstructure(barrel, DoZeroRocsBPO3,(*struct_iter).second);
} else if (ROCMapToReset && shell==PixelBarrelName::pI && twoDimOnlyLayDisk){
if (DoZeroRocsBPI1 && layer==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsBPI1,(*struct_iter).second);
if (DoZeroRocsBPI2 && layer==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsBPI2,(*struct_iter).second);
if (DoZeroRocsBPI3 && layer==3) CountZeroROCsInSubstructure(barrel, DoZeroRocsBPI3,(*struct_iter).second);
}

nBPIXDigis = nBPIXDigis + numberOfDigisMod;
for(int i=0; i!=nBPiXmodules; ++i){
if((*struct_iter).first == I_detId[i]){
Expand All @@ -295,6 +391,22 @@ void SiPixelDigiSource::analyze(const edm::Event& iEvent, const edm::EventSetup&
int panel = PixelEndcapName(DetId((*struct_iter).first),pTT,isUpgrade).pannelName();
int module = PixelEndcapName(DetId((*struct_iter).first),pTT,isUpgrade).plaquetteName();
int iter=0; int i=0;

//count Zero Occupancy ROCs in Endcap in the first event after each 10 Ls
if (ROCMapToReset && side==PixelEndcapName::mO && twoDimOnlyLayDisk){
if (DoZeroRocsFMO1 && disk==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsFMO1,(*struct_iter).second);
if (DoZeroRocsFMO2 && disk==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsFMO2,(*struct_iter).second);
} else if (ROCMapToReset && side==PixelEndcapName::mI && twoDimOnlyLayDisk){
if (DoZeroRocsFMI1 && disk==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsFMI1,(*struct_iter).second);
if (DoZeroRocsFMI2 && disk==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsFMI2,(*struct_iter).second);
} else if (ROCMapToReset && side==PixelEndcapName::pO && twoDimOnlyLayDisk){
if (DoZeroRocsFPO1 && disk==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsFPO1,(*struct_iter).second);
if (DoZeroRocsFPO2 && disk==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsFPO2,(*struct_iter).second);
} else if (ROCMapToReset && side==PixelEndcapName::pI && twoDimOnlyLayDisk){
if (DoZeroRocsFPI1 && disk==1) CountZeroROCsInSubstructure(barrel, DoZeroRocsFPI1,(*struct_iter).second);
if (DoZeroRocsFPI2 && disk==2) CountZeroROCsInSubstructure(barrel, DoZeroRocsFPI2,(*struct_iter).second);
}

if(side==PixelEndcapName::mI){
if(disk==1){
i=0;
Expand Down Expand Up @@ -480,28 +592,19 @@ void SiPixelDigiSource::analyze(const edm::Event& iEvent, const edm::EventSetup&
}
}//endif(Endcap && isUpgrade)
} // endif any digis in this module

if (twoDimOnlyLayDisk && lumiSection%10 > 2){
std::pair<int,int> tempPair = (*struct_iter).second->getZeroLoEffROCs();
if (barrel){
NzeroROCs[0] += tempPair.first;
NloEffROCs[0] += tempPair.second;
}
else if (endcap){
NzeroROCs[1] += tempPair.first;
NloEffROCs[1] += tempPair.second;
}
}
} // endfor loop over all modules

if (lumiSection%10> 2){


if (lumiSection%10==0 && ROCMapToReset) {
for (int i =0; i < 2; ++i) NloEffROCs[i] = NloEffROCs[i] - NzeroROCs[i];
if(noOccROCsBarrel) noOccROCsBarrel->setBinContent(1+lumiSection/10, NzeroROCs[0]);
if(loOccROCsBarrel) loOccROCsBarrel->setBinContent(1+lumiSection/10, NloEffROCs[0]);
if(noOccROCsEndcap) noOccROCsEndcap->setBinContent(1+lumiSection/10, NzeroROCs[1]);
if(loOccROCsEndcap) loOccROCsEndcap->setBinContent(1+lumiSection/10, NloEffROCs[1]);
ROCMapToReset=false; // in this way the ROC maps are reset for one event only (the first event in LS multiple of 10
}

if (noOfDisks == 2) { // if (!isUpgrade)
if(meNDigisCHANEndcap_){ for(int j=0; j!=192; j++) if(numberOfDigis[j]>0) meNDigisCHANEndcap_->Fill((float)numberOfDigis[j]);}
if(meNDigisCHANEndcapDms_.at(0)){ for(int j=0; j!=72; j++) if((j<24||j>47)&&numberOfDigis[j]>0) meNDigisCHANEndcapDms_.at(0)->Fill((float)numberOfDigis[j]);}
Expand Down Expand Up @@ -864,5 +967,22 @@ void SiPixelDigiSource::bookMEs(DQMStore::IBooker & iBooker, const edm::EventSet
iBooker.cd(topFolderName_);
}

void SiPixelDigiSource::CountZeroROCsInSubstructure(bool barrel, bool& DoZeroRocs, SiPixelDigiModule* mod){

std::pair<int,int> tempPair = mod->getZeroLoEffROCs();

if (barrel){
NzeroROCs[0] += tempPair.first;
NloEffROCs[0] += tempPair.second;
}
else {
NzeroROCs[1] += tempPair.first;
NloEffROCs[1] += tempPair.second;
}

DoZeroRocs=false;
mod->resetRocMap(); //once got the number of ZeroOccupancy Rocs, reset the ROC map of the corresponding Pixel substructure
}

//define this as a plug-in
DEFINE_FWK_MODULE(SiPixelDigiSource);

0 comments on commit 22d450c

Please sign in to comment.