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

add jets, report #8910

Closed
wants to merge 15 commits into from
Closed

add jets, report #8910

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
6bbb5d7
add var.bins histograms
popovvp Feb 12, 2015
9e45771
add var.bins histograms
popovvp Feb 12, 2015
5dec9ef
add var.bins histograms
popovvp Feb 12, 2015
152dc0f
add var.bins histograms
popovvp Feb 12, 2015
bd2a166
add tower collection
popovvp Feb 16, 2015
9f9d3e5
add tower collection
popovvp Feb 16, 2015
bbf7f6b
add tower collection
popovvp Feb 16, 2015
58e5c9f
add tower collection
popovvp Feb 16, 2015
cdd1c90
Merged refs/pull/7693/head from repository cms-sw
popovvp Apr 29, 2015
76072e3
add jets
popovvp Apr 29, 2015
f48bd43
add jets
popovvp Apr 29, 2015
03fb58e
add report
popovvp Apr 29, 2015
c7c614c
add report
popovvp Apr 29, 2015
11223ce
remove bunch hist
popovvp Apr 29, 2015
6a60584
remove bunch hist
popovvp Apr 29, 2015
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 12 additions & 4 deletions DQM/CastorMonitor/interface/CastorDigiMonitor.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -19,23 +19,31 @@ class CastorDigiMonitor: public CastorBaseMonitor {
CastorDigiMonitor(const edm::ParameterSet& ps);
~CastorDigiMonitor();

void setup(const edm::ParameterSet& ps);
void setup(const edm::ParameterSet& ps);
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &);
void processEvent(const CastorDigiCollection& cast,const CastorDbService& cond, int bunch);

void processEvent(const CastorDigiCollection& cast,const CastorDbService& cond);
int ModSecToIndex(int module, int sector);
private:
std::string subsystemname_;
int fVerbosity;
int ievt_;
float Qrms_DEAD;

MonitorElement* h2QrmsTSvsCh;
MonitorElement* hQIErms[10];
MonitorElement* hTSratio;
MonitorElement* h2TSratio;
MonitorElement* h2status;
MonitorElement* h2digierr;
MonitorElement* h2reportMap;
MonitorElement* hReport;
MonitorElement* h2QtsvsCh;
MonitorElement *h2QmeantsvsCh;
MonitorElement *h2QmeanMap;
MonitorElement *hModule;
MonitorElement *hSector;
MonitorElement* hdigisize;
MonitorElement* hBunchOcc;
// MonitorElement* hBunchOcc;
};

#endif
8 changes: 8 additions & 0 deletions DQM/CastorMonitor/interface/CastorMonitorModule.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,6 +32,8 @@
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"

#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Common/interface/TriggerNames.h"

#include "DataFormats/CastorReco/interface/CastorTower.h"
#include "DataFormats/CastorReco/interface/CastorCluster.h"
Expand Down Expand Up @@ -115,6 +117,12 @@ class CastorMonitorModule : public DQMEDAnalyzer{
edm::InputTag inputLabelDigi_;
edm::InputTag inputLabelRecHitCASTOR_;

edm::InputTag inputLabelCastorTowers_;
typedef std::vector<reco::CastorTower> CastorTowerCollection;

edm::InputTag JetAlgorithm;
// edm::InputTag trigResultsSource;

CastorRecHitMonitor* RecHitMon_;
CastorDigiMonitor* DigiMon_;
CastorLEDMonitor* LedMon_;
Expand Down
28 changes: 25 additions & 3 deletions DQM/CastorMonitor/interface/CastorRecHitMonitor.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,6 +6,18 @@
#include "DQMServices/Core/interface/MonitorElement.h"
//#include "DataFormats/HcalRecHit/interface/CastorRecHit.h"

#include "DataFormats/CastorReco/interface/CastorTower.h"
#include "DataFormats/CastorReco/interface/CastorCluster.h"
#include "DataFormats/CastorReco/interface/CastorJet.h"
#include "DataFormats/JetReco/interface/CastorJetID.h"
#include "RecoJets/JetProducers/interface/CastorJetIDHelper.h"
#include "RecoJets/JetProducers/plugins/CastorJetIDProducer.h"
#include "DataFormats/JetReco/interface/BasicJet.h"
#include "DataFormats/JetReco/interface/BasicJetCollection.h"
#include "DataFormats/JetReco/interface/Jet.h"

#include "DataFormats/Common/interface/TriggerResults.h"

class CastorRecHitMonitor: public CastorBaseMonitor {
public:
CastorRecHitMonitor(const edm::ParameterSet& ps);
Expand All @@ -14,15 +26,25 @@ class CastorRecHitMonitor: public CastorBaseMonitor {
void setup(const edm::ParameterSet& ps);
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &);
void processEvent(const CastorRecHitCollection& castorHits);

void processEventTowers(const reco::CastorTowerCollection& castorTowers);
void processEventJets(const reco::BasicJetCollection& Jets);
private:
int ievt_;
float energyInEachChannel[14][16];
std::string subsystemname;

MonitorElement *hTowerE;
MonitorElement *hTowerDepth;
MonitorElement *h2TowerEMhad;
MonitorElement *hTowerMultipl;
MonitorElement *hJetsMultipl;
MonitorElement *hJetEnergy;
MonitorElement *hJetEta;
MonitorElement *hJetPhi;

TH2F *h2RecHitMap;
//MonitorElement* h2RHchan;
//MonitorElement* h2RHvsSec;
MonitorElement* h2RHchan;
MonitorElement* h2RHvsSec;
MonitorElement* h2RHmap;
MonitorElement* h2RHoccmap;
MonitorElement* h2RHentriesMap;
Expand Down
185 changes: 144 additions & 41 deletions DQM/CastorMonitor/src/CastorDigiMonitor.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,12 +11,24 @@
//// last revision: 31.05.2011 (Panos Katsas) to remove selecting N events for filling the histograms
//****************************************************//
//---- critical revision 26.06.2014 (Vladimir Popov)
// add rms check 15.04.2015 (Vladimir Popov)
//==================================================================//

static int TS_MAX = 10;
static float RatioThresh1 = 0.;
static float QIEerrThreshold = 0.0001;
static double QrmsTS[224][10], QmeanTS[224][10];
const int TSped = 0;

//======================= Constructor ==============================//
CastorDigiMonitor::CastorDigiMonitor(const edm::ParameterSet& ps)
{
subsystemname_=ps.getUntrackedParameter<std::string>("subSystemFolder","Castor");
fVerbosity = ps.getUntrackedParameter<int>("debug",0);
subsystemname_=ps.getUntrackedParameter<std::string>("subSystemFolder","Castor");
RatioThresh1 = ps.getUntrackedParameter<double>("ratioThreshold",0.9);
Qrms_DEAD = ps.getUntrackedParameter<double>("QrmsDead",0.01); //fC
Qrms_DEAD = Qrms_DEAD*Qrms_DEAD;
TS_MAX = ps.getUntrackedParameter<double>("qieTSmax",6);
}

//======================= Destructor ===============================//
Expand All @@ -38,21 +50,35 @@ void CastorDigiMonitor::bookHistograms(DQMStore::IBooker& ibooker,
{
char s[60];
if(fVerbosity>0) std::cout << "CastorDigiMonitor::beginRun (start)" << std::endl;
char sTileIndex[50];
sprintf(sTileIndex,"Tile(=module*16+sector)");

ievt_=0;
// pset_.getUntrackedParameter<std::string>("subSystemFolder","Castor");
// baseFolder_ = rootFolder_+"CastorDigiMonitor";

ibooker.setCurrentFolder(subsystemname_ + "/CastorDigiMonitor");
std::string s2 = "QIE_capID+er+dv";
h2digierr = ibooker.book2D(s2,s2,14,0.,14., 16,0.,16.);
h2digierr->getTH2F()->GetXaxis()->SetTitle("ModuleZ");
h2digierr->getTH2F()->GetYaxis()->SetTitle("SectorPhi");
h2digierr->getTH2F()->SetOption("colz");

sprintf(s,"BunchOccupancy(fC)_all_TS");
hBunchOcc = ibooker.bookProfile(s,s,4000,0,4000, 100,0,1.e10,"");
hBunchOcc->getTProfile()->GetXaxis()->SetTitle("BX");
hBunchOcc->getTProfile()->GetYaxis()->SetTitle("QIE(fC)");
std::string s2 = "CASTOR QIE_capID+er+dv";
h2digierr=ibooker.bookProfile2D(s2,s2,14,0.,14., 16,0.,16.,100,0,1.e10,"");
h2digierr->getTProfile2D()->GetXaxis()->SetTitle("ModuleZ");
h2digierr->getTProfile2D()->GetYaxis()->SetTitle("SectorPhi");
h2digierr->getTProfile2D()->SetMaximum(1.);
h2digierr->getTProfile2D()->SetMinimum(QIEerrThreshold);
h2digierr->getTProfile2D()->SetOption("colz");

sprintf(s,"CASTOR DeadChannelsMap");
h2status = ibooker.book2D(s,s,14,0.,14., 16,0.,16.);
h2status->getTH2F()->GetXaxis()->SetTitle("ModuleZ");
h2status->getTH2F()->GetYaxis()->SetTitle("SectorPhi");
h2status->getTH2F()->SetOption("colz");

sprintf(s,"CASTOR AverageToMaxRatioMap");
h2TSratio = ibooker.book2D(s,s,14,0.,14., 16,0.,16.);
h2TSratio->getTH2F()->GetXaxis()->SetTitle("ModuleZ");
h2TSratio->getTH2F()->GetYaxis()->SetTitle("SectorPhi");
h2TSratio->getTH2F()->SetOption("colz");

sprintf(s,"CASTOR AverageToMaxRatio");
hTSratio = ibooker.book1D(s,s,100,0.,1.);

sprintf(s,"DigiSize");
hdigisize = ibooker.book1D(s,s,20,0.,20.);
Expand All @@ -64,37 +90,61 @@ void CastorDigiMonitor::bookHistograms(DQMStore::IBooker& ibooker,
hSector = ibooker.book1D(s,s,16,0.,16.);
hSector->getTH1F()->GetXaxis()->SetTitle("SectorPhi");
hSector->getTH1F()->GetYaxis()->SetTitle("QIE(fC)");
sprintf(s,"QfC=f(Tile TS) (cumulative)");
h2QtsvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
h2QtsvsCh->getTH2F()->GetXaxis()->SetTitle("Tile(=sector*14+module)");
h2QtsvsCh->getTH2F()->GetYaxis()->SetTitle("TS");
h2QtsvsCh->getTH2F()->SetOption("colz");

sprintf(s,"QfC=f(x=Tile y=TS) (cumulative)");
h2QtsvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
h2QtsvsCh->getTH2F()->GetXaxis()->SetTitle(sTileIndex);
h2QtsvsCh->getTH2F()->GetYaxis()->SetTitle("TS");
h2QtsvsCh->getTH2F()->SetOption("colz");

sprintf(s,"QmeanfC=f(Tile TS)");
h2QmeantsvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
h2QmeantsvsCh->getTH2F()->GetXaxis()->SetTitle("Tile(=sector*14+module)");
h2QmeantsvsCh->getTH2F()->GetYaxis()->SetTitle("TS");
h2QmeantsvsCh->getTH2F()->GetXaxis()->SetTitle(sTileIndex);
h2QmeantsvsCh->getTH2F()->GetYaxis()->SetTitle("Time Slice");
h2QmeantsvsCh->getTH2F()->SetOption("colz");
sprintf(s,"QmeanfC_map(allTS)");

sprintf(s,"QrmsfC=f(Tile TS)");
h2QrmsTSvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
h2QrmsTSvsCh->getTH2F()->GetXaxis()->SetTitle(sTileIndex);
h2QrmsTSvsCh->getTH2F()->GetYaxis()->SetTitle("TS");
h2QrmsTSvsCh->getTH2F()->SetOption("colz");

sprintf(s,"CASTORreportSummaryMap");
h2reportMap = ibooker.book2D(s,s,14, 0,14, 16, 0,16);
h2reportMap->getTH2F()->GetXaxis()->SetTitle("moduleZ");
h2reportMap->getTH2F()->GetYaxis()->SetTitle("sectorPhi");
h2reportMap->getTH2F()->SetOption("colz");

hReport = ibooker.bookFloat("CASTOR reportSummary");

sprintf(s,"QmeanfC_map(allTS)");
h2QmeanMap = ibooker.book2D(s,s,14,0.,14., 16,0.,16.);
h2QmeanMap->getTH2F()->GetXaxis()->SetTitle("ModuleZ");
h2QmeanMap->getTH2F()->GetYaxis()->SetTitle("SectorPhi");
h2QmeanMap->getTH2F()->SetOption("textcolz");

for(int ts=0; ts<=1; ts++) {
sprintf(s,"QIErms_TS=%d",ts);
hQIErms[ts] = ibooker.book1D(s,s,1000,0.,100.);
hQIErms[ts]->getTH1F()->GetXaxis()->SetTitle("QIErms(fC)");
}

for(int ind=0; ind<224; ind++) for(int ts=0; ts<10; ts++)
QrmsTS[ind][ts] = QmeanTS[ind][ts]= 0.;

if(fVerbosity>0) std::cout<<"CastorDigiMonitor::beginRun(end)"<<std::endl;
return;
}


//=============== processEvent =========
void CastorDigiMonitor::processEvent(const CastorDigiCollection& castorDigis,
const CastorDbService& cond, int iBunch) {
const CastorDbService& cond) {
if(fVerbosity>0) std::cout << "CastorDigiMonitor::processEvent (begin)"<< std::endl;

// CaloSamples tool;

if(castorDigis.size() <= 0) {
if(fVerbosity>0) std::cout<<"CastorPSMonitor::processEvent NO Castor Digis"<<std::endl;
return;
if(fVerbosity>0) std::cout<<"CastorPSMonitor::processEvent NO Castor Digis"<<std::endl;
return;
}

for(CastorDigiCollection::const_iterator j=castorDigis.begin();
Expand All @@ -105,42 +155,57 @@ void CastorDigiMonitor::processEvent(const CastorDigiCollection& castorDigis,
int capid1 = digi.sample(0).capid();
hdigisize->Fill(digi.size());
double sum = 0.;
// for (int i=1; i<digi.size(); i++) {
for (int i=0; i<digi.size(); i++) {
int module = digi.id().module()-1;
int sector = digi.id().sector()-1;
if(capid1 < 3) capid1++;
else capid1 = 0;
int capid = digi.sample(i).capid();
int dv = digi.sample(i).dv();
int er = digi.sample(i).er();
int rawd = digi.sample(i).adc();
rawd = rawd&0x7F;
int err = (capid != capid1) | er<<1 | (!dv)<<2; // =0
if(err !=0) h2digierr->Fill(module,sector);
int ind = sector*14 + module;
h2digierr->Fill(module,sector,err);
// if(err !=0) continue;
int ind = ModSecToIndex(module,sector);
h2QtsvsCh->Fill(ind,i,LedMonAdc2fc[rawd]);
sum += LedMonAdc2fc[rawd];
// if(err != 0 && fVerbosity>0)
// std::cout<<"event/idigi=" <<ievt_<<"/" <<i<< " cap_cap1_dv_er: " <<
// capid <<"="<< capid1 <<" "<< dv <<" "<< er<<" "<< err << std::endl;
capid1 = capid;
float q = LedMonAdc2fc[rawd];
sum += q; // sum += LedMonAdc2fc[rawd];
QrmsTS[ind][i] += (q*q);
QmeanTS[ind][i] += q;
if(err != 0 && fVerbosity>0)
std::cout<<"event/idigi=" <<ievt_<<"/"<<i<<" cap=cap1_dv_er_err: "<<
capid <<"="<< capid1 <<" "<< dv <<" "<< er<<" "<< err << std::endl;
if(capid1 < 3) capid1 = capid+1;
else capid1 = 0;
}
hBunchOcc->Fill(iBunch,sum);
// hBunchOcc->Fill(iBunch,sum);
} //end for(CastorDigiCollection::const_iterator ...

ievt_++;
if(ievt_ %100 == 0) {

const float repChanBAD = 0.9;
const float repChanWarning = 0.95;
if(ievt_ %100 != 0) return;
float ModuleSum[14], SectorSum[16];
for(int m=0; m<14; m++) ModuleSum[m]=0.;
for(int s=0; s<16; s++) SectorSum[s]=0.;
for(int mod=0; mod<14; mod++) for(int sec=0; sec<16; sec++) {
for(int ts=0; ts<=1; ts++) {
int ind = ModSecToIndex(mod,sec);
double Qmean = QmeanTS[ind][ts]/ievt_;
double Qrms = sqrt(QrmsTS[ind][ts]/ievt_ - Qmean*Qmean);
hQIErms[ts]->Fill(Qrms);
}

double sum=0.;
for(int ts=1; ts<=10; ts++) {
int ind = sec*14 + mod +1;
for(int ts=1; ts<=TS_MAX; ts++) {
int ind = ModSecToIndex(mod,sec) + 1;
double a=h2QtsvsCh->getTH2F()->GetBinContent(ind,ts);
h2QmeantsvsCh->getTH2F()->SetBinContent(ind,ts,a/double(ievt_));
sum += a;
double Qmean = QmeanTS[ind-1][ts-1]/ievt_;
double Qrms = QrmsTS[ind-1][ts-1]/ievt_ - Qmean*Qmean;
h2QrmsTSvsCh->getTH2F()->SetBinContent(ind,ts,sqrt(Qrms));
}
sum /= double(ievt_);
ModuleSum[mod] += sum;
Expand All @@ -153,8 +218,46 @@ void CastorDigiMonitor::processEvent(const CastorDigiCollection& castorDigis,
hModule->getTH1F()->SetBinContent(mod+1,ModuleSum[mod]);
for(int sec=0; sec<16; sec++)
hSector->getTH1F()->SetBinContent(sec+1,SectorSum[sec]);
} //end if(ievt_ %100 == 0) {

if(fVerbosity>0) std::cout << "CastorDigiMonitor::processEvent (end)"<< std::endl;
int nGoodCh = 0;
for(int mod=0; mod<14; mod++) for(int sec=0; sec<16;sec++) {
int ind = ModSecToIndex(mod,sec);
double Qmean = QmeanTS[ind][TSped]/ievt_;
double Qrms = QrmsTS[ind][TSped]/ievt_ - Qmean*Qmean;
float ChanStatus = 0.;
if(Qrms < Qrms_DEAD) ChanStatus = 1.;
h2status->getTH2F()->SetBinContent(mod+1,sec+1,ChanStatus);

int tsm = 0;
float am=0.;
for(int ts=0; ts<TS_MAX; ts++) {
float a = h2QmeantsvsCh->getTH2F()->GetBinContent(ind+1,ts+1);
if(am < a) {am = a; tsm = ts;}
}

double sum = 0.;
for(int ts=0; ts<TS_MAX; ts++) if(ts != tsm)
sum += h2QmeantsvsCh->getTH2F()->GetBinContent(ind+1,ts+1);
float r = 0.;
if(am > 0.) r = sum/(TS_MAX-1)/am;
h2TSratio->getTH2F()->SetBinContent(mod+1,sec+1,r);
hTSratio->Fill(r);
float statusTS = 1.0;
if(r > RatioThresh1) statusTS = repChanWarning;
else if(r > 0.99) statusTS = repChanBAD;
float gChanStatus = statusTS;
if(ChanStatus > 0.) gChanStatus = repChanBAD; // RMS
if(h2digierr->getTProfile2D()->GetBinContent(mod+1,sec+1)>QIEerrThreshold)
gChanStatus = repChanBAD;
h2reportMap->getTH2F()->SetBinContent(mod+1,sec+1,gChanStatus);
if(gChanStatus > repChanBAD) ++nGoodCh;
}
hReport->Fill(float(nGoodCh)/224.);
return;
}

int CastorDigiMonitor::ModSecToIndex(int module, int sector) {
int ind = sector + module*16;
if(ind>223) ind=223;
return(ind);
}
Loading