|
|
|
#include <SiPixelDigiSource.h>
Description: <one line="" class="" summary>="">
Usage: <usage>
Definition at line 49 of file SiPixelDigiSource.h.
| typedef edm::DetSet<PixelDigi>::const_iterator SiPixelDigiSource::DigiIterator |
Definition at line 54 of file SiPixelDigiSource.h.
| SiPixelDigiSource::SiPixelDigiSource | ( | const edm::ParameterSet & | conf | ) | [explicit] |
Definition at line 46 of file SiPixelDigiSource.cc.
References cmsCodeRules::cppFunctionSkipper::operator, and theDMBE.
: conf_(iConfig), src_( conf_.getParameter<edm::InputTag>( "src" ) ), saveFile( conf_.getUntrackedParameter<bool>("saveFile",false) ), isPIB( conf_.getUntrackedParameter<bool>("isPIB",false) ), slowDown( conf_.getUntrackedParameter<bool>("slowDown",false) ), modOn( conf_.getUntrackedParameter<bool>("modOn",true) ), twoDimOn( conf_.getUntrackedParameter<bool>("twoDimOn",true) ), twoDimModOn( conf_.getUntrackedParameter<bool>("twoDimModOn",true) ), twoDimOnlyLayDisk( conf_.getUntrackedParameter<bool>("twoDimOnlyLayDisk",false) ), hiRes( conf_.getUntrackedParameter<bool>("hiRes",false) ), reducedSet( conf_.getUntrackedParameter<bool>("reducedSet",false) ), ladOn( conf_.getUntrackedParameter<bool>("ladOn",false) ), layOn( conf_.getUntrackedParameter<bool>("layOn",false) ), phiOn( conf_.getUntrackedParameter<bool>("phiOn",false) ), ringOn( conf_.getUntrackedParameter<bool>("ringOn",false) ), bladeOn( conf_.getUntrackedParameter<bool>("bladeOn",false) ), diskOn( conf_.getUntrackedParameter<bool>("diskOn",false) ), bigEventSize( conf_.getUntrackedParameter<int>("bigEventSize",1000) ) { theDMBE = edm::Service<DQMStore>().operator->(); LogInfo ("PixelDQM") << "SiPixelDigiSource::SiPixelDigiSource: Got DQM BackEnd interface"<<endl; }
| SiPixelDigiSource::~SiPixelDigiSource | ( | ) |
Definition at line 71 of file SiPixelDigiSource.cc.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
LogInfo ("PixelDQM") << "SiPixelDigiSource::~SiPixelDigiSource: Destructor"<<endl;
}
| void SiPixelDigiSource::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 134 of file SiPixelDigiSource.cc.
References Reference_intrackfit_cff::barrel, bigEventSize, PixelEndcapName::bladeName(), bladeOn, edm::EventBase::bunchCrossing(), PixelEndcapName::diskName(), diskOn, Reference_intrackfit_cff::endcap, eventNo, MonitorElement::Fill(), DQMStore::get(), edm::Event::getByLabel(), PixelEndcapName::halfCylinder(), i, I_detId, I_fedId, I_linkId1, I_linkId2, collect_tpl::input, edm::HandleBase::isValid(), j, ladOn, layOn, edm::EventBase::luminosityBlock(), PixelEndcapName::mI, PixelEndcapName::mO, modOn, nBPIXDigis, nDigisA, nDigisB, nDigisPerChan, nDigisPerDisk, nDigisPerFed, nFPIXDigis, numberOfDigis, cmsCodeRules::cppFunctionSkipper::operator, PixelEndcapName::pannelName(), phiOn, PixelEndcapName::pI, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, PixelEndcapName::pO, reducedSet, MonitorElement::Reset(), ringOn, MonitorElement::setBinContent(), slowDown, src_, DetId::subdetId(), theDMBE, thePixelStructure, twoDimModOn, twoDimOn, and twoDimOnlyLayDisk.
{
eventNo++;
//cout<<"BIGFATEVENTNUMBER: "<<eventNo<<endl;
// get input data
edm::Handle< edm::DetSetVector<PixelDigi> > input;
iEvent.getByLabel( src_, input );
if (!input.isValid()) return;
// Get DQM interface
DQMStore* theDMBE = edm::Service<DQMStore>().operator->();
//float iOrbitSec = iEvent.orbitNumber()/11223.;
int bx = iEvent.bunchCrossing();
//long long tbx = (long long)iEvent.orbitNumber() * 3564 + bx;
int lumiSection = (int)iEvent.luminosityBlock();
int nEventDigis = 0; int nActiveModules = 0;
//int nEventBPIXDigis = 0; int nEventFPIXDigis = 0;
if(modOn){
MonitorElement* meReset = theDMBE->get("Pixel/averageDigiOccupancy");
//if(meReset && eventNo%1000==0){
if(meReset && lumiSection%8==0){
meReset->Reset();
nBPIXDigis = 0;
nFPIXDigis = 0;
for(int i=0; i!=40; i++) nDigisPerFed[i]=0;
}
}
std::map<uint32_t,SiPixelDigiModule*>::iterator struct_iter;
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;
for (struct_iter = thePixelStructure.begin() ; struct_iter != thePixelStructure.end() ; struct_iter++) {
int numberOfDigisMod = (*struct_iter).second->fill(*input, modOn,
ladOn, layOn, phiOn,
bladeOn, diskOn, ringOn,
twoDimOn, reducedSet, twoDimModOn, twoDimOnlyLayDisk,
nDigisA, nDigisB);
if(numberOfDigisMod>0){
nEventDigis = nEventDigis + numberOfDigisMod;
nActiveModules++;
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((*struct_iter).first >= 302055684 && (*struct_iter).first <= 302197792 ){ // Barrel
if(barrel){ // Barrel
//cout<<"AAbpix: "<<numberOfDigisMod<<" + "<<nBPIXDigis<<" = ";
nBPIXDigis = nBPIXDigis + numberOfDigisMod;
//cout<<nBPIXDigis<<endl;
for(int i=0; i!=768; i++){
//cout<<"\t\t\t bpix: "<<i<<" , "<<(*struct_iter).first<<" , "<<I_detId[i]<<endl;
if((*struct_iter).first == I_detId[i]){
//if(I_fedId[i]>=32&&I_fedId[i]<=39) std::cout<<"Attention: a BPIX module matched to an FPIX FED!"<<std::endl;
nDigisPerFed[I_fedId[i]]=nDigisPerFed[I_fedId[i]]+numberOfDigisMod;
//cout<<"BPIX: "<<i<<" , "<<I_fedId[i]<<" , "<<numberOfDigisMod<<" , "<<nDigisPerFed[I_fedId[i]]<<endl;
int index1 = I_fedId[i]*36+(I_linkId1[i]-1); int index2 = I_fedId[i]*36+(I_linkId2[i]-1);
if(nDigisA>0) nDigisPerChan[index1]=nDigisPerChan[index1]+nDigisA;
if(nDigisB>0) nDigisPerChan[index2]=nDigisPerChan[index2]+nDigisB;
//cout<<"BPIX: "<<i<<" , "<<I_linkId1[i]<<" , "<<I_linkId2[i]<<" , "<<nDigisA<<" , "<<nDigisB<<endl;
i=767;
}
}
//}else if((*struct_iter).first >= 343999748 && (*struct_iter).first <= 352477708 ){ // Endcap
}else if(endcap){ // Endcap
//cout<<"AAfpix: "<<nFPIXDigis<<" = ";
nFPIXDigis = nFPIXDigis + numberOfDigisMod;
//cout<<nFPIXDigis<<endl;
PixelEndcapName::HalfCylinder side = PixelEndcapName(DetId((*struct_iter).first)).halfCylinder();
int disk = PixelEndcapName(DetId((*struct_iter).first)).diskName();
int blade = PixelEndcapName(DetId((*struct_iter).first)).bladeName();
int panel = PixelEndcapName(DetId((*struct_iter).first)).pannelName();
//std::cout<<"Endcap: "<<side<<" , "<<disk<<" , "<<blade<<" , "<<panel<<" , "<<std::endl;
int iter=-1; int i=0;
if(side==PixelEndcapName::mI){
if(disk==1){
i=0;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}else if(disk==2){
i=24;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}
}else if(side==PixelEndcapName::mO){
if(disk==1){
i=48;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}else if(disk==2){
i=72;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}
}else if(side==PixelEndcapName::pI){
if(disk==1){
i=96;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}else if(disk==2){
i=120;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}
}else if(side==PixelEndcapName::pO){
if(disk==1){
i=144;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}else if(disk==2){
i=168;
if(blade==1){ if(panel==1) iter=i; else if(panel==2) iter=i+1; }
if(blade==2){ if(panel==1) iter=i+2; else if(panel==2) iter=i+3; }
if(blade==3){ if(panel==1) iter=i+4; else if(panel==2) iter=i+5; }
if(blade==4){ if(panel==1) iter=i+6; else if(panel==2) iter=i+7; }
if(blade==5){ if(panel==1) iter=i+8; else if(panel==2) iter=i+9; }
if(blade==6){ if(panel==1) iter=i+10; else if(panel==2) iter=i+11; }
if(blade==7){ if(panel==1) iter=i+12; else if(panel==2) iter=i+13; }
if(blade==8){ if(panel==1) iter=i+14; else if(panel==2) iter=i+15; }
if(blade==9){ if(panel==1) iter=i+16; else if(panel==2) iter=i+17; }
if(blade==10){ if(panel==1) iter=i+18; else if(panel==2) iter=i+19; }
if(blade==11){ if(panel==1) iter=i+20; else if(panel==2) iter=i+21; }
if(blade==12){ if(panel==1) iter=i+22; else if(panel==2) iter=i+23; }
}
}
numberOfDigis[iter]=numberOfDigis[iter]+numberOfDigisMod;
//if(side==PixelEndcapName::pO||side==PixelEndcapName::pI){
// if(disk==2){
// std::cout<<"status: "<<iter<<","<<side<<","<<disk<<","<<blade<<","<<panel<<","<<numberOfDigisMod<<","<<numberOfDigis[iter]<<std::endl;
//}}
for(int i=768; i!=1440; i++){
//cout<<"\t\t\t fpix: "<<i<<" , "<<(*struct_iter).first<<" , "<<I_detId[i]<<endl;
if((*struct_iter).first == I_detId[i]){
//if(I_fedId[i]<32||I_fedId[i]>39) std::cout<<"Attention: an FPIX module matched to a BPIX FED!"<<std::endl;
nDigisPerFed[I_fedId[i]]=nDigisPerFed[I_fedId[i]]+numberOfDigisMod;
//cout<<"FPIX: "<<i<<" , "<<I_fedId[i]<<" , "<<nDigisPerFed[I_fedId[i]]<<endl;
i=1439;
}
}
} //endif Barrel/Endcap
//cout<<"numberOfDigis: "<<numberOfDigisMod<<" , nBPIXDigis: "<<nBPIXDigis<<" , nFPIXDigis: "<<nFPIXDigis<<endl;
// digi occupancy per individual FED channel:
} // endif any digis in this module
} // endfor loop over all modules
// if(lumiSection>lumSec){ lumSec = lumiSection; nLumiSecs++; }
// if(nEventDigis>bigEventSize) nBigEvents++;
// if(nLumiSecs%5==0){
MonitorElement* meE; MonitorElement* meE1; MonitorElement* meE2; MonitorElement* meE3; MonitorElement* meE4;
meE=theDMBE->get("Pixel/Endcap/ALLMODS_ndigisCHAN_Endcap");
if(meE){ for(int j=0; j!=192; j++) if(numberOfDigis[j]>0) meE->Fill((float)numberOfDigis[j]);}
meE1=theDMBE->get("Pixel/Endcap/ALLMODS_ndigisCHAN_EndcapDm1");
if(meE1){ for(int j=0; j!=72; j++) if((j<24||j>47)&&numberOfDigis[j]>0) meE1->Fill((float)numberOfDigis[j]);}
meE2=theDMBE->get("Pixel/Endcap/ALLMODS_ndigisCHAN_EndcapDm2");
if(meE2){ for(int j=24; j!=96; j++) if((j<48||j>71)&&numberOfDigis[j]>0) meE2->Fill((float)numberOfDigis[j]);}
meE3=theDMBE->get("Pixel/Endcap/ALLMODS_ndigisCHAN_EndcapDp1");
if(meE3){ for(int j=96; j!=168; j++) if((j<120||j>143)&&numberOfDigis[j]>0) meE3->Fill((float)numberOfDigis[j]);}
meE4=theDMBE->get("Pixel/Endcap/ALLMODS_ndigisCHAN_EndcapDp2");
if(meE4){ for(int j=120; j!=192; j++) if((j<144||j>167)&&numberOfDigis[j]>0) meE4->Fill((float)numberOfDigis[j]);}
MonitorElement* me1; MonitorElement* me2; MonitorElement* me3; MonitorElement* me4; MonitorElement* me5;
MonitorElement* me6; MonitorElement* me7; MonitorElement* me8; MonitorElement* me9; MonitorElement* me10; MonitorElement* me11;
MonitorElement* me12; MonitorElement* me13; MonitorElement* me14; MonitorElement* me15; MonitorElement* me16; MonitorElement* me17;
MonitorElement* me18; MonitorElement* me19; MonitorElement* me20; MonitorElement* me21; MonitorElement* me22; MonitorElement* me23;
MonitorElement* me24; MonitorElement* me25; MonitorElement* me26; MonitorElement* me27; MonitorElement* me28; MonitorElement* me29;
MonitorElement* me30; MonitorElement* me31; MonitorElement* me32; MonitorElement* me33; MonitorElement* me34; MonitorElement* me35;
MonitorElement* me36;
me1=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh1");
if(me1){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+0]>0) me1->Fill((float)nDigisPerChan[i*36+0]);}
me2=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh2");
if(me2){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+1]>0) me2->Fill((float)nDigisPerChan[i*36+1]);}
me3=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh3");
if(me3){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+2]>0) me3->Fill((float)nDigisPerChan[i*36+2]);}
me4=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh4");
if(me4){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+3]>0) me4->Fill((float)nDigisPerChan[i*36+3]);}
me5=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh5");
if(me5){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+4]>0) me5->Fill((float)nDigisPerChan[i*36+4]);}
me6=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh6");
if(me6){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+5]>0) me6->Fill((float)nDigisPerChan[i*36+5]);}
me7=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh7");
if(me7){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+6]>0) me7->Fill((float)nDigisPerChan[i*36+6]);}
me8=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh8");
if(me8){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+7]>0) me8->Fill((float)nDigisPerChan[i*36+7]);}
me9=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh9");
if(me9){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+8]>0) me9->Fill((float)nDigisPerChan[i*36+8]);}
me10=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh10");
if(me10){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+9]>0) me10->Fill((float)nDigisPerChan[i*36+9]);}
me11=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh11");
if(me11){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+10]>0) me11->Fill((float)nDigisPerChan[i*36+10]);}
me12=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh12");
if(me12){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+11]>0) me12->Fill((float)nDigisPerChan[i*36+11]);}
me13=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh13");
if(me13){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+12]>0) me13->Fill((float)nDigisPerChan[i*36+12]);}
me14=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh14");
if(me14){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+13]>0) me14->Fill((float)nDigisPerChan[i*36+13]);}
me15=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh15");
if(me15){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+14]>0) me15->Fill((float)nDigisPerChan[i*36+14]);}
me16=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh16");
if(me16){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+15]>0) me16->Fill((float)nDigisPerChan[i*36+15]);}
me17=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh17");
if(me17){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+16]>0) me17->Fill((float)nDigisPerChan[i*36+16]);}
me18=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh18");
if(me18){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+17]>0) me18->Fill((float)nDigisPerChan[i*36+17]);}
me19=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh19");
if(me19){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+18]>0) me19->Fill((float)nDigisPerChan[i*36+18]);}
me20=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh20");
if(me20){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+19]>0) me20->Fill((float)nDigisPerChan[i*36+19]);}
me21=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh21");
if(me21){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+20]>0) me21->Fill((float)nDigisPerChan[i*36+20]);}
me22=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh22");
if(me22){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+21]>0) me22->Fill((float)nDigisPerChan[i*36+21]);}
me23=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh23");
if(me23){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+22]>0) me23->Fill((float)nDigisPerChan[i*36+22]);}
me24=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh24");
if(me24){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+23]>0) me24->Fill((float)nDigisPerChan[i*36+23]);}
me25=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh25");
if(me25){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+24]>0) me25->Fill((float)nDigisPerChan[i*36+24]);}
me26=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh26");
if(me26){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+25]>0) me26->Fill((float)nDigisPerChan[i*36+25]);}
me27=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh27");
if(me27){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+26]>0) me27->Fill((float)nDigisPerChan[i*36+26]);}
me28=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh28");
if(me28){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+27]>0) me28->Fill((float)nDigisPerChan[i*36+27]);}
me29=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh29");
if(me29){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+28]>0) me29->Fill((float)nDigisPerChan[i*36+28]);}
me30=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh30");
if(me30){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+29]>0) me30->Fill((float)nDigisPerChan[i*36+29]);}
me31=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh31");
if(me31){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+30]>0) me31->Fill((float)nDigisPerChan[i*36+30]);}
me32=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh32");
if(me32){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+31]>0) me32->Fill((float)nDigisPerChan[i*36+31]);}
me33=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh33");
if(me33){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+32]>0) me33->Fill((float)nDigisPerChan[i*36+32]);}
me34=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh34");
if(me34){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+33]>0) me34->Fill((float)nDigisPerChan[i*36+33]);}
me35=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh35");
if(me35){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+34]>0) me35->Fill((float)nDigisPerChan[i*36+34]);}
me36=theDMBE->get("Pixel/Barrel/ALLMODS_ndigisCHAN_BarrelCh36");
if(me36){ for(int i=0; i!=32; i++) if(nDigisPerChan[i*36+35]>0) me36->Fill((float)nDigisPerChan[i*36+35]);}
// Rate of events with >N digis:
MonitorElement* meX1;
if(nEventDigis>bigEventSize){
meX1 = theDMBE->get("Pixel/bigEventRate");
if(meX1) meX1->Fill(lumiSection,1./23.);
}
//std::cout<<"nEventDigis: "<<nEventDigis<<" , nLumiSecs: "<<nLumiSecs<<" , nBigEvents: "<<nBigEvents<<std::endl;
// Rate of pixel events and total number of pixel events per BX:
MonitorElement* meX2; MonitorElement* meX3;
if(nActiveModules>=4){
meX2 = theDMBE->get("Pixel/pixEvtsPerBX");
if(meX2) meX2->Fill(float(bx));
meX3 = theDMBE->get("Pixel/pixEventRate");
if(meX3) meX3->Fill(lumiSection, 1./23.);
}
// Actual digi occupancy in a FED compared to average digi occupancy per FED
MonitorElement* meX4; MonitorElement* meX5;
meX4 = theDMBE->get("Pixel/averageDigiOccupancy");
meX5 = theDMBE->get("Pixel/avgfedDigiOccvsLumi");
if(meX4){
int maxfed=0;
for(int i=0; i!=32; i++){
if(nDigisPerFed[i]>maxfed) maxfed=nDigisPerFed[i];
}
for(int i=0; i!=40; i++){
float averageOcc = 0.;
if(i<32){
float averageBPIXFed = float(nBPIXDigis-maxfed)/31.;
if(averageBPIXFed>0.) averageOcc = nDigisPerFed[i]/averageBPIXFed;
//cout<<"\t BPIX i: "<<i<<" , "<<nBPIXDigis<<" , "<<averageBPIXFed<<" , "<<nDigisPerFed[i]<<" , "<<averageOcc<<endl;
}else{
float averageFPIXFed = float(nFPIXDigis)/8.;
if(averageFPIXFed>0.) averageOcc = nDigisPerFed[i]/averageFPIXFed;
//cout<<"\t FPIX i: "<<i<<" , "<<nFPIXDigis<<" , "<<averageFPIXFed<<" , "<<nDigisPerFed[i]<<" , "<<averageOcc<<endl;
}
meX4->setBinContent(i+1,averageOcc);
int lumiSections8 = int(lumiSection/8);
if (modOn){
if (meX5){
meX5->setBinContent(1+lumiSections8, i+1, averageOcc);
}//endif meX5
}//endif modOn
}
}
// slow down...
if(slowDown) usleep(10000);
}
| void SiPixelDigiSource::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 79 of file SiPixelDigiSource.cc.
References firstRun, I_detId, I_fedId, I_linkId1, I_linkId2, I_name, recoMuon::in, and EdgesToViz::infile.
{
firstRun = true;
// find a FED# for the current detId:
ifstream infile(edm::FileInPath("DQM/SiPixelMonitorClient/test/detId.dat").fullPath().c_str(),ios::in);
int nModsInFile=0;
assert(!infile.fail());
while(!infile.eof()&&nModsInFile<1440) {
infile >> I_name[nModsInFile] >> I_detId[nModsInFile] >> I_fedId[nModsInFile] >> I_linkId1[nModsInFile] >> I_linkId2[nModsInFile];
//cout<<nModsInFile<<" , "<<I_name[nModsInFile]<<" , "<<I_detId[nModsInFile]<<" , "<<I_fedId[nModsInFile]<<" , "<<I_linkId[nModsInFile]<<endl; ;
nModsInFile++;
}
infile.close();
}
| void SiPixelDigiSource::beginRun | ( | const edm::Run & | r, |
| edm::EventSetup const & | iSetup | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 93 of file SiPixelDigiSource.cc.
References bladeOn, bookMEs(), buildStructure(), diskOn, eventNo, firstRun, hiRes, i, ladOn, layOn, lumSec, modOn, nBigEvents, nBPIXDigis, nDigisPerDisk, nDigisPerFed, nFPIXDigis, nLumiSecs, phiOn, ringOn, and twoDimOn.
{
LogInfo ("PixelDQM") << " SiPixelDigiSource::beginJob - Initialisation ... " << std::endl;
LogInfo ("PixelDQM") << "Mod/Lad/Lay/Phi " << modOn << "/" << ladOn << "/"
<< 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){
eventNo = 0;
lumSec = 0;
nLumiSecs = 0;
nBigEvents = 0;
nBPIXDigis = 0;
nFPIXDigis = 0;
for(int i=0; i!=40; i++) nDigisPerFed[i]=0;
for(int i=0; i!=4; i++) nDigisPerDisk[i]=0;
// Build map
buildStructure(iSetup);
// Book Monitoring Elements
bookMEs();
firstRun = false;
}
}
| void SiPixelDigiSource::bookMEs | ( | ) | [virtual] |
Create folder tree and book histograms
Definition at line 563 of file SiPixelDigiSource.cc.
References averageDigiOccupancy, avgfedDigiOccvsLumi, bigEventRate, bigEventSize, bladeOn, DQMStore::book1D(), DQMStore::book2D(), DQMStore::cd(), conf_, diskOn, Exception, hiRes, isPIB, ladOn, layOn, LogDebug, meNDigisCHANBarrel_, meNDigisCHANBarrelCh10_, meNDigisCHANBarrelCh11_, meNDigisCHANBarrelCh12_, meNDigisCHANBarrelCh13_, meNDigisCHANBarrelCh14_, meNDigisCHANBarrelCh15_, meNDigisCHANBarrelCh16_, meNDigisCHANBarrelCh17_, meNDigisCHANBarrelCh18_, meNDigisCHANBarrelCh19_, meNDigisCHANBarrelCh1_, meNDigisCHANBarrelCh20_, meNDigisCHANBarrelCh21_, meNDigisCHANBarrelCh22_, meNDigisCHANBarrelCh23_, meNDigisCHANBarrelCh24_, meNDigisCHANBarrelCh25_, meNDigisCHANBarrelCh26_, meNDigisCHANBarrelCh27_, meNDigisCHANBarrelCh28_, meNDigisCHANBarrelCh29_, meNDigisCHANBarrelCh2_, meNDigisCHANBarrelCh30_, meNDigisCHANBarrelCh31_, meNDigisCHANBarrelCh32_, meNDigisCHANBarrelCh33_, meNDigisCHANBarrelCh34_, meNDigisCHANBarrelCh35_, meNDigisCHANBarrelCh36_, meNDigisCHANBarrelCh3_, meNDigisCHANBarrelCh4_, meNDigisCHANBarrelCh5_, meNDigisCHANBarrelCh6_, meNDigisCHANBarrelCh7_, meNDigisCHANBarrelCh8_, meNDigisCHANBarrelCh9_, meNDigisCHANBarrelL1_, meNDigisCHANBarrelL2_, meNDigisCHANBarrelL3_, meNDigisCHANEndcap_, meNDigisCHANEndcapDm1_, meNDigisCHANEndcapDm2_, meNDigisCHANEndcapDp1_, meNDigisCHANEndcapDp2_, meNDigisCOMBBarrel_, meNDigisCOMBEndcap_, modOn, cmsCodeRules::cppFunctionSkipper::operator, phiOn, pixEventRate, pixEvtsPerBX, DQMStore::pwd(), reducedSet, ringOn, MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), SiPixelFolderOrganizer::setModuleFolder(), theDMBE, thePixelStructure, indexGen::title, twoDimModOn, twoDimOn, and twoDimOnlyLayDisk.
Referenced by beginRun().
{
// Get DQM interface
DQMStore* theDMBE = edm::Service<DQMStore>().operator->();
theDMBE->setCurrentFolder("Pixel");
char title[80]; sprintf(title, "Rate of events with >%i digis;LumiSection;Rate [Hz]",bigEventSize);
bigEventRate = theDMBE->book1D("bigEventRate",title,5000,0.,5000.);
char title1[80]; sprintf(title1, "Pixel events vs. BX;BX;# events");
pixEvtsPerBX = theDMBE->book1D("pixEvtsPerBX",title1,3565,0.,3565.);
char title2[80]; sprintf(title2, "Rate of Pixel events;LumiSection;Rate [Hz]");
pixEventRate = theDMBE->book1D("pixEventRate",title2,5000,0.,5000.);
char title3[80]; sprintf(title3, "Average digi occupancy per FED;FED;NDigis/<NDigis>");
averageDigiOccupancy = theDMBE->book1D("averageDigiOccupancy",title3,40,-0.5,39.5);
if(modOn){
char title4[80]; sprintf(title4, "FED Digi Occupancy (NDigis/<NDigis>) vs LumiSections;Lumi Section;FED");
avgfedDigiOccvsLumi = theDMBE->book2D ("avgfedDigiOccvsLumi", title4, 400,0., 3200., 40, -0.5, 39.5);
}
std::map<uint32_t,SiPixelDigiModule*>::iterator struct_iter;
SiPixelFolderOrganizer theSiPixelFolder;
for(struct_iter = thePixelStructure.begin(); struct_iter != thePixelStructure.end(); struct_iter++){
if(modOn){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first)){
(*struct_iter).second->book( conf_,0,twoDimOn,hiRes, reducedSet, twoDimModOn);
} else {
if(!isPIB) throw cms::Exception("LogicError")
<< "[SiPixelDigiSource::bookMEs] Creation of DQM folder failed";
}
}
if(ladOn){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first,1)){
(*struct_iter).second->book( conf_,1,twoDimOn,hiRes, reducedSet);
} else {
LogDebug ("PixelDQM") << "PROBLEM WITH LADDER-FOLDER\n";
}
}
if(layOn || twoDimOnlyLayDisk){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first,2)){
(*struct_iter).second->book( conf_,2,twoDimOn,hiRes, reducedSet, twoDimOnlyLayDisk);
} else {
LogDebug ("PixelDQM") << "PROBLEM WITH LAYER-FOLDER\n";
}
}
if(phiOn){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first,3)){
(*struct_iter).second->book( conf_,3,twoDimOn,hiRes, reducedSet);
} else {
LogDebug ("PixelDQM") << "PROBLEM WITH PHI-FOLDER\n";
}
}
if(bladeOn){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first,4)){
(*struct_iter).second->book( conf_,4,twoDimOn,hiRes, reducedSet);
} else {
LogDebug ("PixelDQM") << "PROBLEM WITH BLADE-FOLDER\n";
}
}
if(diskOn || twoDimOnlyLayDisk){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first,5)){
(*struct_iter).second->book( conf_,5,twoDimOn,hiRes, reducedSet, twoDimOnlyLayDisk);
} else {
LogDebug ("PixelDQM") << "PROBLEM WITH DISK-FOLDER\n";
}
}
if(ringOn){
if(theSiPixelFolder.setModuleFolder((*struct_iter).first,6)){
(*struct_iter).second->book( conf_,6,twoDimOn,hiRes, reducedSet);
} else {
LogDebug ("PixelDQM") << "PROBLEM WITH RING-FOLDER\n";
}
}
}
std::string currDir = theDMBE->pwd();
theDMBE->cd("Pixel/Barrel");
meNDigisCOMBBarrel_ = theDMBE->book1D("ALLMODS_ndigisCOMB_Barrel","Number of Digis",200,0.,400.);
meNDigisCOMBBarrel_->setAxisTitle("Number of digis per module per event",1);
meNDigisCHANBarrel_ = theDMBE->book1D("ALLMODS_ndigisCHAN_Barrel","Number of Digis",100,0.,1000.);
meNDigisCHANBarrel_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelL1_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelL1","Number of Digis L1",100,0.,1000.);
meNDigisCHANBarrelL1_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelL2_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelL2","Number of Digis L2",100,0.,1000.);
meNDigisCHANBarrelL2_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelL3_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelL3","Number of Digis L3",100,0.,1000.);
meNDigisCHANBarrelL3_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh1_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh1","Number of Digis Ch1",100,0.,1000.);
meNDigisCHANBarrelCh1_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh2_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh2","Number of Digis Ch2",100,0.,1000.);
meNDigisCHANBarrelCh2_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh3_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh3","Number of Digis Ch3",100,0.,1000.);
meNDigisCHANBarrelCh3_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh4_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh4","Number of Digis Ch4",100,0.,1000.);
meNDigisCHANBarrelCh4_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh5_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh5","Number of Digis Ch5",100,0.,1000.);
meNDigisCHANBarrelCh5_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh6_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh6","Number of Digis Ch6",100,0.,1000.);
meNDigisCHANBarrelCh6_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh7_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh7","Number of Digis Ch7",100,0.,1000.);
meNDigisCHANBarrelCh7_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh8_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh8","Number of Digis Ch8",100,0.,1000.);
meNDigisCHANBarrelCh8_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh9_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh9","Number of Digis Ch9",100,0.,1000.);
meNDigisCHANBarrelCh9_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh10_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh10","Number of Digis Ch10",100,0.,1000.);
meNDigisCHANBarrelCh10_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh11_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh11","Number of Digis Ch11",100,0.,1000.);
meNDigisCHANBarrelCh11_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh12_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh12","Number of Digis Ch12",100,0.,1000.);
meNDigisCHANBarrelCh12_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh13_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh13","Number of Digis Ch13",100,0.,1000.);
meNDigisCHANBarrelCh13_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh14_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh14","Number of Digis Ch14",100,0.,1000.);
meNDigisCHANBarrelCh14_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh15_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh15","Number of Digis Ch15",100,0.,1000.);
meNDigisCHANBarrelCh15_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh16_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh16","Number of Digis Ch16",100,0.,1000.);
meNDigisCHANBarrelCh16_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh17_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh17","Number of Digis Ch17",100,0.,1000.);
meNDigisCHANBarrelCh17_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh18_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh18","Number of Digis Ch18",100,0.,1000.);
meNDigisCHANBarrelCh18_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh19_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh19","Number of Digis Ch19",100,0.,1000.);
meNDigisCHANBarrelCh19_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh20_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh20","Number of Digis Ch20",100,0.,1000.);
meNDigisCHANBarrelCh20_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh21_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh21","Number of Digis Ch21",100,0.,1000.);
meNDigisCHANBarrelCh21_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh22_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh22","Number of Digis Ch22",100,0.,1000.);
meNDigisCHANBarrelCh22_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh23_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh23","Number of Digis Ch23",100,0.,1000.);
meNDigisCHANBarrelCh23_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh24_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh24","Number of Digis Ch24",100,0.,1000.);
meNDigisCHANBarrelCh24_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh25_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh25","Number of Digis Ch25",100,0.,1000.);
meNDigisCHANBarrelCh25_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh26_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh26","Number of Digis Ch26",100,0.,1000.);
meNDigisCHANBarrelCh26_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh27_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh27","Number of Digis Ch27",100,0.,1000.);
meNDigisCHANBarrelCh27_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh28_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh28","Number of Digis Ch28",100,0.,1000.);
meNDigisCHANBarrelCh28_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh29_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh29","Number of Digis Ch29",100,0.,1000.);
meNDigisCHANBarrelCh29_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh30_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh30","Number of Digis Ch30",100,0.,1000.);
meNDigisCHANBarrelCh30_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh31_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh31","Number of Digis Ch31",100,0.,1000.);
meNDigisCHANBarrelCh31_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh32_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh32","Number of Digis Ch32",100,0.,1000.);
meNDigisCHANBarrelCh32_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh33_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh33","Number of Digis Ch33",100,0.,1000.);
meNDigisCHANBarrelCh33_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh34_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh34","Number of Digis Ch34",100,0.,1000.);
meNDigisCHANBarrelCh34_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh35_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh35","Number of Digis Ch35",100,0.,1000.);
meNDigisCHANBarrelCh35_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANBarrelCh36_ = theDMBE->book1D("ALLMODS_ndigisCHAN_BarrelCh36","Number of Digis Ch36",100,0.,1000.);
meNDigisCHANBarrelCh36_->setAxisTitle("Number of digis per FED channel per event",1);
theDMBE->cd("Pixel/Endcap");
meNDigisCOMBEndcap_ = theDMBE->book1D("ALLMODS_ndigisCOMB_Endcap","Number of Digis",200,0.,400.);
meNDigisCOMBEndcap_->setAxisTitle("Number of digis per module per event",1);
meNDigisCHANEndcap_ = theDMBE->book1D("ALLMODS_ndigisCHAN_Endcap","Number of Digis",100,0.,1000.);
meNDigisCHANEndcap_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANEndcapDp1_ = theDMBE->book1D("ALLMODS_ndigisCHAN_EndcapDp1","Number of Digis Disk p1",100,0.,1000.);
meNDigisCHANEndcapDp1_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANEndcapDp2_ = theDMBE->book1D("ALLMODS_ndigisCHAN_EndcapDp2","Number of Digis Disk p2",100,0.,1000.);
meNDigisCHANEndcapDp2_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANEndcapDm1_ = theDMBE->book1D("ALLMODS_ndigisCHAN_EndcapDm1","Number of Digis Disk m1",100,0.,1000.);
meNDigisCHANEndcapDm1_->setAxisTitle("Number of digis per FED channel per event",1);
meNDigisCHANEndcapDm2_ = theDMBE->book1D("ALLMODS_ndigisCHAN_EndcapDm2","Number of Digis Disk m2",100,0.,1000.);
meNDigisCHANEndcapDm2_->setAxisTitle("Number of digis per FED channel per event",1);
theDMBE->cd(currDir);
}
| void SiPixelDigiSource::buildStructure | ( | edm::EventSetup const & | iSetup | ) | [virtual] |
Definition at line 500 of file SiPixelDigiSource.cc.
References PixelEndcapName::bladeName(), PixelEndcapName::diskName(), edm::EventSetup::get(), PixelEndcapName::halfCylinder(), isPIB, LogDebug, PixelEndcapName::pannelName(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, PixelEndcapName::plaquetteName(), DetId::rawId(), PixelGeomDetUnit::specificTopology(), DetId::subdetId(), and thePixelStructure.
Referenced by beginRun().
{
LogInfo ("PixelDQM") <<" SiPixelDigiSource::buildStructure" ;
edm::ESHandle<TrackerGeometry> pDD;
iSetup.get<TrackerDigiGeometryRecord>().get( pDD );
LogVerbatim ("PixelDQM") << " *** Geometry node for TrackerGeom is "<<&(*pDD)<<std::endl;
LogVerbatim ("PixelDQM") << " *** I have " << pDD->dets().size() <<" detectors"<<std::endl;
LogVerbatim ("PixelDQM") << " *** I have " << pDD->detTypes().size() <<" types"<<std::endl;
for(TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++){
if(dynamic_cast<PixelGeomDetUnit*>((*it))!=0){
DetId detId = (*it)->geographicalId();
const GeomDetUnit * geoUnit = pDD->idToDetUnit( detId );
const PixelGeomDetUnit * pixDet = dynamic_cast<const PixelGeomDetUnit*>(geoUnit);
int nrows = (pixDet->specificTopology()).nrows();
int ncols = (pixDet->specificTopology()).ncolumns();
if(detId.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) {
if(isPIB) continue;
LogDebug ("PixelDQM") << " ---> Adding Barrel Module " << detId.rawId() << endl;
uint32_t id = detId();
SiPixelDigiModule* theModule = new SiPixelDigiModule(id, ncols, nrows);
thePixelStructure.insert(pair<uint32_t,SiPixelDigiModule*> (id,theModule));
} else if(detId.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
LogDebug ("PixelDQM") << " ---> Adding Endcap Module " << detId.rawId() << endl;
uint32_t id = detId();
SiPixelDigiModule* theModule = new SiPixelDigiModule(id, ncols, nrows);
PixelEndcapName::HalfCylinder side = PixelEndcapName(DetId(id)).halfCylinder();
int disk = PixelEndcapName(DetId(id)).diskName();
int blade = PixelEndcapName(DetId(id)).bladeName();
int panel = PixelEndcapName(DetId(id)).pannelName();
int module = PixelEndcapName(DetId(id)).plaquetteName();
char sside[80]; sprintf(sside, "HalfCylinder_%i",side);
char sdisk[80]; sprintf(sdisk, "Disk_%i",disk);
char sblade[80]; sprintf(sblade, "Blade_%02i",blade);
char spanel[80]; sprintf(spanel, "Panel_%i",panel);
char smodule[80];sprintf(smodule,"Module_%i",module);
std::string side_str = sside;
std::string disk_str = sdisk;
bool mask = side_str.find("HalfCylinder_1")!=string::npos||
side_str.find("HalfCylinder_2")!=string::npos||
side_str.find("HalfCylinder_4")!=string::npos||
disk_str.find("Disk_2")!=string::npos;
// clutch to take all of FPIX, but no BPIX:
mask = false;
if(isPIB && mask) continue;
thePixelStructure.insert(pair<uint32_t,SiPixelDigiModule*> (id,theModule));
}
}
}
LogInfo ("PixelDQM") << " *** Pixel Structure Size " << thePixelStructure.size() << endl;
}
| void SiPixelDigiSource::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 121 of file SiPixelDigiSource.cc.
References conf_, edm::ParameterSet::getParameter(), download_sqlite_cfg::outputFile, DQMStore::save(), saveFile, and theDMBE.
{
if(saveFile) {
LogInfo ("PixelDQM") << " SiPixelDigiSource::endJob - Saving Root File " << std::endl;
std::string outputFile = conf_.getParameter<std::string>("outputFile");
theDMBE->save( outputFile.c_str() );
}
}
Definition at line 92 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 93 of file SiPixelDigiSource.h.
Referenced by bookMEs().
MonitorElement* SiPixelDigiSource::bigEventRate [private] |
Definition at line 89 of file SiPixelDigiSource.h.
Referenced by bookMEs().
int SiPixelDigiSource::bigEventSize [private] |
Definition at line 142 of file SiPixelDigiSource.h.
bool SiPixelDigiSource::bladeOn [private] |
Definition at line 79 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
edm::ParameterSet SiPixelDigiSource::conf_ [private] |
Definition at line 65 of file SiPixelDigiSource.h.
bool SiPixelDigiSource::diskOn [private] |
Definition at line 79 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
int SiPixelDigiSource::eventNo [private] |
Definition at line 80 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginRun().
bool SiPixelDigiSource::firstRun [private] |
Definition at line 144 of file SiPixelDigiSource.h.
Referenced by beginJob(), and beginRun().
bool SiPixelDigiSource::hiRes [private] |
Definition at line 74 of file SiPixelDigiSource.h.
Referenced by beginRun(), and bookMEs().
unsigned int SiPixelDigiSource::I_detId[1440] [private] |
Definition at line 147 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginJob().
int SiPixelDigiSource::I_fedId[1440] [private] |
Definition at line 148 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginJob().
int SiPixelDigiSource::I_linkId1[1440] [private] |
Definition at line 149 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginJob().
int SiPixelDigiSource::I_linkId2[1440] [private] |
Definition at line 150 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginJob().
std::string SiPixelDigiSource::I_name[1440] [private] |
Definition at line 146 of file SiPixelDigiSource.h.
Referenced by beginJob().
bool SiPixelDigiSource::isPIB [private] |
Definition at line 68 of file SiPixelDigiSource.h.
Referenced by bookMEs(), and buildStructure().
bool SiPixelDigiSource::ladOn [private] |
Definition at line 77 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
bool SiPixelDigiSource::layOn [private] |
Definition at line 77 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
int SiPixelDigiSource::lumSec [private] |
Definition at line 81 of file SiPixelDigiSource.h.
Referenced by beginRun().
Definition at line 96 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 109 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 110 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 111 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 112 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 113 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 114 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 115 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 116 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 117 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 118 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 100 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 119 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 120 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 121 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 122 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 123 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 124 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 125 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 126 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 127 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 128 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 101 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 129 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 130 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 131 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 132 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 133 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 134 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 135 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 102 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 103 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 104 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 105 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 106 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 107 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 108 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 97 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 98 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 99 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 136 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 139 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 140 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 137 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 138 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 94 of file SiPixelDigiSource.h.
Referenced by bookMEs().
Definition at line 95 of file SiPixelDigiSource.h.
Referenced by bookMEs().
bool SiPixelDigiSource::modOn [private] |
Definition at line 70 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
int SiPixelDigiSource::nBigEvents [private] |
Definition at line 86 of file SiPixelDigiSource.h.
Referenced by beginRun().
int SiPixelDigiSource::nBPIXDigis [private] |
Definition at line 87 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginRun().
int SiPixelDigiSource::nDigisA [private] |
Definition at line 155 of file SiPixelDigiSource.h.
Referenced by analyze().
int SiPixelDigiSource::nDigisB [private] |
Definition at line 156 of file SiPixelDigiSource.h.
Referenced by analyze().
int SiPixelDigiSource::nDigisPerChan[1152] [private] |
Definition at line 152 of file SiPixelDigiSource.h.
Referenced by analyze().
int SiPixelDigiSource::nDigisPerDisk[4] [private] |
Definition at line 153 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginRun().
int SiPixelDigiSource::nDigisPerFed[40] [private] |
Definition at line 151 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginRun().
int SiPixelDigiSource::nFPIXDigis [private] |
Definition at line 88 of file SiPixelDigiSource.h.
Referenced by analyze(), and beginRun().
int SiPixelDigiSource::nLumiSecs [private] |
Definition at line 82 of file SiPixelDigiSource.h.
Referenced by beginRun().
int SiPixelDigiSource::numberOfDigis[192] [private] |
Definition at line 154 of file SiPixelDigiSource.h.
Referenced by analyze().
bool SiPixelDigiSource::phiOn [private] |
Definition at line 77 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
MonitorElement* SiPixelDigiSource::pixEventRate [private] |
Definition at line 91 of file SiPixelDigiSource.h.
Referenced by bookMEs().
MonitorElement* SiPixelDigiSource::pixEvtsPerBX [private] |
Definition at line 90 of file SiPixelDigiSource.h.
Referenced by bookMEs().
bool SiPixelDigiSource::reducedSet [private] |
Definition at line 75 of file SiPixelDigiSource.h.
bool SiPixelDigiSource::ringOn [private] |
Definition at line 79 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
bool SiPixelDigiSource::saveFile [private] |
Definition at line 67 of file SiPixelDigiSource.h.
Referenced by endJob().
bool SiPixelDigiSource::slowDown [private] |
Definition at line 69 of file SiPixelDigiSource.h.
Referenced by analyze().
edm::InputTag SiPixelDigiSource::src_ [private] |
Definition at line 66 of file SiPixelDigiSource.h.
Referenced by analyze().
DQMStore* SiPixelDigiSource::theDMBE [private] |
Definition at line 83 of file SiPixelDigiSource.h.
Referenced by analyze(), bookMEs(), endJob(), and SiPixelDigiSource().
std::map<uint32_t,SiPixelDigiModule*> SiPixelDigiSource::thePixelStructure [private] |
Definition at line 84 of file SiPixelDigiSource.h.
Referenced by analyze(), bookMEs(), and buildStructure().
bool SiPixelDigiSource::twoDimModOn [private] |
Definition at line 72 of file SiPixelDigiSource.h.
bool SiPixelDigiSource::twoDimOn [private] |
Definition at line 71 of file SiPixelDigiSource.h.
Referenced by analyze(), beginRun(), and bookMEs().
bool SiPixelDigiSource::twoDimOnlyLayDisk [private] |
Definition at line 73 of file SiPixelDigiSource.h.
1.7.3