Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

add pixel roc occupancy dqm plots #15820

Merged
merged 6 commits into from Nov 15, 2016
Merged

add pixel roc occupancy dqm plots #15820

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
98cad9d
add pixel roc occupancy dqm plots
gouskos Sep 12, 2016
7d6a50d
clean up some comments
gouskos Sep 12, 2016
943cd4d
fix compiler error
gouskos Sep 24, 2016
b99ecbd
fix override
gouskos Oct 6, 2016
933fc73
wteaks and add some comments
gouskos Oct 27, 2016
a9699dc
remove if statements
gouskos Nov 7, 2016
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
13 changes: 11 additions & 2 deletions DQM/SiPixelMonitorCluster/interface/SiPixelClusterSource.h
View file
Open in desktop
Expand Up @@ -64,8 +64,8 @@

typedef edmNew::DetSet<SiPixelCluster>::const_iterator ClusterIterator;

virtual void analyze(const edm::Event&, const edm::EventSetup&);
virtual void dqmBeginRun(const edm::Run&, edm::EventSetup const&) ;
virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
virtual void dqmBeginRun(const edm::Run&, edm::EventSetup const&) override;
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;

virtual void buildStructure(edm::EventSetup const&);
Expand All @@ -76,6 +76,7 @@
private:
edm::ParameterSet conf_;
edm::InputTag src_;
edm::InputTag digisrc_;
bool saveFile;
bool isPIB;
bool slowDown;
Expand Down Expand Up @@ -105,11 +106,19 @@
MonitorElement* meClusFpixPProf;
MonitorElement* meClusFpixMProf;

std::vector<MonitorElement*> meZeroRocBPIX;
MonitorElement* meZeroRocFPIX;

int noOfLayers;
int noOfDisks;

void getrococcupancy(DetId detId,const edm::DetSetVector<PixelDigi> diginp,const TrackerTopology* const tTopo,std::vector<MonitorElement*> meinput);
void getrococcupancye(DetId detId,const edmNew::DetSetVector<SiPixelCluster> & clustColl,const TrackerTopology* const pTT,
edm::ESHandle<TrackerGeometry> pDD,MonitorElement* meinput);

//define Token(-s)
edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > srcToken_;
edm::EDGetTokenT<edm::DetSetVector<PixelDigi> > digisrcToken_;
};

#endif
1 change: 1 addition & 0 deletions DQM/SiPixelMonitorCluster/python/SiPixelMonitorCluster_cfi.py
View file
Open in desktop
Expand Up @@ -8,6 +8,7 @@
SiPixelClusterSource = cms.EDAnalyzer("SiPixelClusterSource",
TopFolderName = cms.string('Pixel'),
src = cms.InputTag("siPixelClusters"),
digisrc = cms.InputTag("siPixelDigis"),
outputFile = cms.string('Pixel_DQM_Cluster.root'),
saveFile = cms.untracked.bool(False),
slowDown = cms.untracked.bool(False),
Expand Down
156 changes: 155 additions & 1 deletion DQM/SiPixelMonitorCluster/src/SiPixelClusterSource.cc
View file
Open in desktop
Expand Up @@ -48,6 +48,7 @@ using namespace edm;
SiPixelClusterSource::SiPixelClusterSource(const edm::ParameterSet& iConfig) :
conf_(iConfig),
src_( conf_.getParameter<edm::InputTag>( "src" ) ),
digisrc_( conf_.getParameter<edm::InputTag>( "digisrc" ) ),
saveFile( conf_.getUntrackedParameter<bool>("saveFile",false) ),
isPIB( conf_.getUntrackedParameter<bool>("isPIB",false) ),
slowDown( conf_.getUntrackedParameter<bool>("slowDown",false) ),
Expand All @@ -70,8 +71,10 @@ SiPixelClusterSource::SiPixelClusterSource(const edm::ParameterSet& iConfig) :

//set Token(-s)
srcToken_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(conf_.getParameter<edm::InputTag>("src"));
digisrcToken_ = consumes<edm::DetSetVector<PixelDigi> >(conf_.getParameter<edm::InputTag>("digisrc"));
firstRun = true;
topFolderName_ = conf_.getParameter<std::string>("TopFolderName");

}


Expand Down Expand Up @@ -152,6 +155,32 @@ void SiPixelClusterSource::bookHistograms(DQMStore::IBooker & iBooker, edm::Run
meClusFpixMProf = iBooker.bookProfile(ss1.str(),ss2.str(),2400,0.,150,0,0,"");
meClusFpixMProf->getTH1()->SetCanExtend(TH1::kAllAxes);

iBooker.setCurrentFolder(topFolderName_+"/Barrel");
for (int i = 1; i <= noOfLayers; i++) {
int ybins = -1; float ymin = 0.; float ymax = 0.;
if (i==1) { ybins = 42; ymin = -10.5; ymax = 10.5; }
if (i==2) { ybins = 66; ymin = -16.5; ymax = 16.5; }
if (i==3) { ybins = 90; ymin = -22.5; ymax = 22.5; }
ss1.str(std::string()); ss1 << "pix_bar Occ_roc_online_" + digisrc_. label() + "_layer_" << i;
ss2.str(std::string()); ss2 << "Pixel Barrel Occupancy, ROC level (Online): Layer " << i;
meZeroRocBPIX.push_back(iBooker.book2D(ss1.str(),ss2.str(),72,-4.5,4.5,ybins,ymin,ymax));
meZeroRocBPIX.at(i-1)->setAxisTitle("ROC / Module",1);
meZeroRocBPIX.at(i-1)->setAxisTitle("ROC / Ladder",2);
}

iBooker.setCurrentFolder(topFolderName_+"/Endcap");
meZeroRocFPIX = iBooker.book2D("ROC_endcap_occupancy","Pixel Endcap Occupancy, ROC level (Online)",72, -4.5, 4.5,288,-12.5,12.5);
meZeroRocFPIX->setBinLabel(1, "Disk-2 Pnl2",1);
meZeroRocFPIX->setBinLabel(9, "Disk-2 Pnl1",1);
meZeroRocFPIX->setBinLabel(19, "Disk-1 Pnl2",1);
meZeroRocFPIX->setBinLabel(27, "Disk-1 Pnl1",1);
meZeroRocFPIX->setBinLabel(41, "Disk+1 Pnl1",1);
meZeroRocFPIX->setBinLabel(49, "Disk+1 Pnl2",1);
meZeroRocFPIX->setBinLabel(59, "Disk+2 Pnl1",1);
meZeroRocFPIX->setBinLabel(67, "Disk+2 Pnl2",1);
meZeroRocFPIX->setAxisTitle("Blades in Inner (>0) / Outer(<) Halves",2);
meZeroRocFPIX->setAxisTitle("ROC occupancy",3);

}

//------------------------------------------------------------------
Expand All @@ -176,11 +205,21 @@ void SiPixelClusterSource::analyze(const edm::Event& iEvent, const edm::EventSet
// get input data
edm::Handle< edmNew::DetSetVector<SiPixelCluster> > input;
iEvent.getByToken(srcToken_, input);
auto const & clustColl = *(input.product());

edm::ESHandle<TrackerGeometry> pDD;
iSetup.get<TrackerDigiGeometryRecord> ().get (pDD);
const TrackerGeometry* tracker = &(* pDD);

edm::Handle< edm::DetSetVector<PixelDigi> > digiinput;
iEvent.getByToken( digisrcToken_, digiinput );
const edm::DetSetVector<PixelDigi> diginp = *(digiinput.product());

edm::ESHandle<TrackerTopology> tTopoHandle;
iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
const TrackerTopology *pTT = tTopoHandle.product();


int lumiSection = (int)iEvent.luminosityBlock();
int nEventFpixClusters = 0;

Expand Down Expand Up @@ -218,7 +257,23 @@ void SiPixelClusterSource::analyze(const edm::Event& iEvent, const edm::EventSet

//std::cout<<"nEventFpixClusters: "<<nEventFpixClusters<<" , nLumiSecs: "<<nLumiSecs<<" , nBigEvents: "<<nBigEvents<<std::endl;


for(TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++){

DetId detId = (*it)->geographicalId();

// fill barrel
if (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) {
getrococcupancy(detId,diginp,pTT,meZeroRocBPIX);
}

// fill endcap
if (detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
getrococcupancye(detId,clustColl,pTT,pDD,meZeroRocFPIX);
}

}



// slow down...
if(slowDown) usleep(10000);
Expand Down Expand Up @@ -373,5 +428,104 @@ void SiPixelClusterSource::bookMEs(DQMStore::IBooker & iBooker, const edm::Event

}


void SiPixelClusterSource::getrococcupancy(DetId detId,const edm::DetSetVector<PixelDigi> diginp,const TrackerTopology* const tTopo,std::vector<MonitorElement*> meinput) {

edm::DetSetVector<PixelDigi>::const_iterator ipxsearch = diginp.find(detId);
if( ipxsearch != diginp.end() ) {

// Look at digis now
edm::DetSet<PixelDigi>::const_iterator pxdi;
for (pxdi = ipxsearch->begin(); pxdi != ipxsearch->end(); pxdi++) {

bool isHalfModule = PixelBarrelName(DetId(detId),tTopo,isUpgrade).isHalfModule();
int DBlayer = PixelBarrelName(DetId(detId),tTopo,isUpgrade).layerName();
int DBmodule = PixelBarrelName(DetId(detId),tTopo,isUpgrade).moduleName();
int DBladder = PixelBarrelName(DetId(detId),tTopo,isUpgrade).ladderName();
int DBshell = PixelBarrelName(DetId(detId),tTopo,isUpgrade).shell();

// add sign to the modules
if (DBshell==1 || DBshell==2) { DBmodule = -DBmodule; }
if (DBshell==1 || DBshell==3) { DBladder = -DBladder; }

int col = pxdi->column();
int row = pxdi->row();

float modsign = (float)DBmodule/(abs((float)DBmodule));
float ladsign = (float)DBladder/(abs((float)DBladder));
float rocx = ((float)col/(52.*8.))*modsign + ((float)DBmodule-(modsign)*0.5);
float rocy = ((float)row/(80.*2.))*ladsign + ((float)DBladder-(ladsign)*0.5);

// do the flip where need
bool flip = false;
if ( (DBladder%2==0) && (!isHalfModule) ) { flip = true; }
if ((flip) && (DBladder>0)) {
if ( ( ((float)DBladder-(ladsign)*0.5)<=rocy) && (rocy<(float)DBladder)) { rocy = rocy + ladsign*0.5; }
else if ( ( ((float)DBladder)<=rocy) && (rocy<((float)DBladder+(ladsign)*0.5)) ) { rocy = rocy - ladsign*0.5; }
}

// tweak border effect for negative modules/ladders
if (modsign<0) { rocx = rocx -0.0001; }
if (ladsign<0) { rocy = rocy -0.0001; } else { rocy = rocy +0.0001; }
if (abs(DBladder)==1) { rocy = rocy + ladsign*0.5; } //take care of the half module
meinput[DBlayer-1]->Fill(rocx,rocy);
} // end of looping over pxdi
}
}



void SiPixelClusterSource::getrococcupancye(DetId detId, const edmNew::DetSetVector<SiPixelCluster> & clustColl,
const TrackerTopology* const pTT,edm::ESHandle<TrackerGeometry> pDD,MonitorElement* meinput) {

edmNew::DetSetVector<SiPixelCluster>::const_iterator ipxsearch = clustColl.find(detId);
if( ipxsearch != clustColl.end() ) {

// Look at clusters now
edmNew::DetSet<SiPixelCluster>::const_iterator pxclust;
for (pxclust = ipxsearch->begin(); pxclust != ipxsearch->end(); pxclust++) {

const GeomDetUnit * geoUnit = pDD->idToDetUnit( detId );
const PixelGeomDetUnit * pixDet = dynamic_cast<const PixelGeomDetUnit*>(geoUnit);
const PixelTopology * topol = &(pixDet->specificTopology());
LocalPoint clustlp = topol->localPosition( MeasurementPoint(pxclust->x(),pxclust->y()) );
GlobalPoint clustgp = geoUnit->surface().toGlobal( clustlp );

float xclust = pxclust->x();
float yclust = pxclust->y();
float z = clustgp.z();

int pxfside = PixelEndcapName(detId,pTT,isUpgrade).halfCylinder();
int pxfpanel = PixelEndcapName(detId,pTT,isUpgrade).pannelName();
int pxfmodule = PixelEndcapName(detId,pTT,isUpgrade).plaquetteName();
int pxfdisk = PixelEndcapName(detId,pTT,isUpgrade).diskName();
int pxfblade = PixelEndcapName(detId,pTT,isUpgrade).bladeName();

if ( (pxfside==1) || (pxfside==3) ) { pxfblade = -1.*pxfblade; }

if (z<0.) { pxfdisk = -1.*pxfdisk; }

int clu_sdpx = ((pxfdisk>0) ? 1 : -1) * (2 * (abs(pxfdisk) - 1) + pxfpanel);
int binselx = (pxfpanel==1&&(pxfmodule==1||pxfmodule==4)) ? (pxfmodule==1) : ((pxfpanel==1&& xclust<80.0)||(pxfpanel==2&&xclust>=80.0));
int nperpan = 2 * pxfmodule + pxfpanel - 1 + binselx;
int clu_roc_binx = ((pxfdisk>0) ? nperpan : 9 - nperpan) + (clu_sdpx + 4) * 8 - 2 * ((abs(pxfdisk)==1) ? pxfdisk : 0);

int clu_roc_biny = -99.;
int nrocly = pxfmodule + pxfpanel;
for (int i=0; i<nrocly; i++) {
int j = (pxfdisk<0) ? i : nrocly - 1 - i;
if (yclust>=(j*52.0)&& yclust<((j+1)*52.0))
clu_roc_biny = 6 - nrocly + 2 * i + ((pxfblade>0) ? pxfblade-1 : pxfblade + 12)*12 + 1;
}
if (pxfblade>0) { clu_roc_biny = clu_roc_biny+144; }

meinput->setBinContent(clu_roc_binx,clu_roc_biny, meinput->getBinContent(clu_roc_binx,clu_roc_biny)+1);
meinput->setBinContent(clu_roc_binx,clu_roc_biny+1, meinput->getBinContent(clu_roc_binx,clu_roc_biny+1)+1);

}
}

}

//define this as a plug-in
DEFINE_FWK_MODULE(SiPixelClusterSource);
4 changes: 2 additions & 2 deletions DQM/SiPixelMonitorTrack/interface/SiPixelHitEfficiencySource.h
View file
Open in desktop
Expand Up @@ -46,9 +46,9 @@ class SiPixelHitEfficiencySource : public DQMEDAnalyzer {
explicit SiPixelHitEfficiencySource(const edm::ParameterSet&);
~SiPixelHitEfficiencySource();

virtual void dqmBeginRun(const edm::Run& r, edm::EventSetup const& iSetup);
virtual void dqmBeginRun(const edm::Run& r, edm::EventSetup const& iSetup) override;
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
virtual void analyze(const edm::Event&, const edm::EventSetup&);
virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
virtual void fillClusterProbability(int , int, bool, double );

private:
Expand Down
12 changes: 9 additions & 3 deletions DQM/SiPixelMonitorTrack/interface/SiPixelTrackResidualSource.h
View file
Open in desktop
Expand Up @@ -51,9 +51,9 @@ class SiPixelTrackResidualSource : public DQMEDAnalyzer {
explicit SiPixelTrackResidualSource(const edm::ParameterSet&);
~SiPixelTrackResidualSource();

virtual void dqmBeginRun(const edm::Run& r, edm::EventSetup const& iSetup);
virtual void dqmBeginRun(const edm::Run& r, edm::EventSetup const& iSetup) override;
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
virtual void analyze(const edm::Event&, const edm::EventSetup&);
virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
void getrococcupancy(DetId detId,const edm::DetSetVector<PixelDigi> diginp,const TrackerTopology* const tTopo,std::vector<MonitorElement*> meinput);
void triplets(double x1,double y1,double z1,double x2,double y2,double z2,double x3,double y3,double z3,
double ptsig, double & dc,double & dz, double kap);
Expand Down Expand Up @@ -181,7 +181,13 @@ class SiPixelTrackResidualSource : public DQMEDAnalyzer {
std::vector<MonitorElement*> meZeroRocLadvsModOffTrackBarrel;

MonitorElement* meHitProbability;
MonitorElement* meRocBladevsDiskEndcap;
MonitorElement* meRocBladevsDiskEndcapOnTrk;
MonitorElement* meRocBladevsDiskEndcapOffTrk;

void getepixrococcupancyontrk(const TrackerTopology* const tTopo, TransientTrackingRecHit::ConstRecHitPointer hit,
float xclust, float yclust, float z, MonitorElement* meinput);
void getepixrococcupancyofftrk(DetId detId, const TrackerTopology* const tTopo,
float xclust, float yclust, float z, MonitorElement* meinput);

int noOfLayers;
int noOfDisks;
Expand Down