|
|
|
#include <EELaserTask.h>
Definition at line 20 of file EELaserTask.h.
| EELaserTask::EELaserTask | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 33 of file EELaserTask.cc.
References dqmStore_, EcalPnDiodeDigiCollection_, EcalRawDataCollection_, EcalUncalibratedRecHitCollection_, EEDigiCollection_, enableCleanup_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, init_, laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, mergeRuns_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, cppFunctionSkipper::operator, and prefixME_.
{
init_ = false;
dqmStore_ = edm::Service<DQMStore>().operator->();
prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
EEDigiCollection_ = ps.getParameter<edm::InputTag>("EEDigiCollection");
EcalPnDiodeDigiCollection_ = ps.getParameter<edm::InputTag>("EcalPnDiodeDigiCollection");
EcalUncalibratedRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection");
// vector of enabled wavelengths (Default to all 4)
laserWavelengths_.reserve(4);
for ( unsigned int i = 1; i <= 4; i++ ) laserWavelengths_.push_back(i);
laserWavelengths_ = ps.getUntrackedParameter<std::vector<int> >("laserWavelengths", laserWavelengths_);
for (int i = 0; i < 18; i++) {
meShapeMapL1_[i] = 0;
meAmplMapL1_[i] = 0;
meTimeMapL1_[i] = 0;
meAmplPNMapL1_[i] = 0;
mePnAmplMapG01L1_[i] = 0;
mePnPedMapG01L1_[i] = 0;
mePnAmplMapG16L1_[i] = 0;
mePnPedMapG16L1_[i] = 0;
meShapeMapL2_[i] = 0;
meAmplMapL2_[i] = 0;
meTimeMapL2_[i] = 0;
meAmplPNMapL2_[i] = 0;
mePnAmplMapG01L2_[i] = 0;
mePnPedMapG01L2_[i] = 0;
mePnAmplMapG16L2_[i] = 0;
mePnPedMapG16L2_[i] = 0;
meShapeMapL3_[i] = 0;
meAmplMapL3_[i] = 0;
meTimeMapL3_[i] = 0;
meAmplPNMapL3_[i] = 0;
mePnAmplMapG01L3_[i] = 0;
mePnPedMapG01L3_[i] = 0;
mePnAmplMapG16L3_[i] = 0;
mePnPedMapG16L3_[i] = 0;
meShapeMapL4_[i] = 0;
meAmplMapL4_[i] = 0;
meTimeMapL4_[i] = 0;
meAmplPNMapL4_[i] = 0;
mePnAmplMapG01L4_[i] = 0;
mePnPedMapG01L4_[i] = 0;
mePnAmplMapG16L4_[i] = 0;
mePnPedMapG16L4_[i] = 0;
}
for(int i=0; i<2; i++){
meAmplSummaryMapL1_[i] = 0;
meAmplSummaryMapL2_[i] = 0;
meAmplSummaryMapL3_[i] = 0;
meAmplSummaryMapL4_[i] = 0;
}
}
| EELaserTask::~EELaserTask | ( | ) | [virtual] |
| void EELaserTask::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Analyze.
Implements edm::EDAnalyzer.
Definition at line 726 of file EELaserTask.cc.
References ecalMGPA::adc(), EcalMGPASample::adc(), EcalEndcap, EcalPnDiodeDigiCollection_, EcalRawDataCollection_, EcalUncalibratedRecHitCollection_, EEDigiCollection_, MonitorElement::Fill(), edm::Event::getByLabel(), NumbersPn::getPNs(), i, Numbers::icEE(), ievt_, init_, NumbersPn::ipnEE(), Numbers::iSM(), ecalpyutils::ism(), EEDetId::ix(), EcalDCCHeaderBlock::LASER_GAP, EcalDCCHeaderBlock::LASER_STD, LogDebug, max(), meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, min, Numbers::RtHalf(), environment_file_cff::runType, EcalDataFrame::sample(), setup(), Numbers::subDet(), and dtDQMClient_cfg::threshold.
{
bool enable = false;
int runType[18];
for (int i=0; i<18; i++) runType[i] = -1;
unsigned rtHalf[18];
for (int i=0; i<18; i++) rtHalf[i] = -1;
int waveLength[18];
for (int i=0; i<18; i++) waveLength[i] = -1;
edm::Handle<EcalRawDataCollection> dcchs;
if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
int ism = Numbers::iSM( *dcchItr, EcalEndcap );
runType[ism-1] = dcchItr->getRunType();
rtHalf[ism-1] = dcchItr->getRtHalf();
waveLength[ism-1] = dcchItr->getEventSettings().wavelength;
if ( dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_STD ||
dcchItr->getRunType() == EcalDCCHeaderBlock::LASER_GAP ) enable = true;
}
} else {
edm::LogWarning("EELaserTask") << EcalRawDataCollection_ << " not available";
}
if ( ! enable ) return;
if ( ! init_ ) this->setup();
ievt_++;
bool numPN[80];
float adcPN[80];
for ( int i = 0; i < 80; i++ ) {
numPN[i] = false;
adcPN[i] = 0.;
}
std::vector<int> PNs;
PNs.reserve(12);
edm::Handle<EEDigiCollection> digis;
if ( e.getByLabel(EEDigiCollection_, digis) ) {
int maxpos[10];
for(int i(0); i < 10; i++)
maxpos[i] = 0;
int nReadouts(0);
for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
EEDetId id = digiItr->id();
int ism = Numbers::iSM( id );
if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
nReadouts++;
EEDataFrame dataframe = (*digiItr);
int iMax(-1);
float max(0.);
float min(4096.);
for (int i = 0; i < 10; i++) {
int adc = dataframe.sample(i).adc();
if(adc > max){
max = adc;
iMax = i;
}
if(adc < min)
min = adc;
}
if(iMax >= 0 && max - min > 20.)
maxpos[iMax] += 1;
}
int threshold(nReadouts / 2);
enable = false;
for(int i(0); i < 10; i++){
if(maxpos[i] > threshold){
enable = true;
break;
}
}
if(!enable) return;
int need = digis->size();
LogDebug("EELaserTask") << "event " << ievt_ << " digi collection size " << need;
for ( EEDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
EEDetId id = digiItr->id();
int ix = id.ix();
int iy = id.iy();
int ism = Numbers::iSM( id );
if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
int ic = Numbers::icEE(ism, ix, iy);
EEDataFrame dataframe = (*digiItr);
for (int i = 0; i < 10; i++) {
int adc = dataframe.sample(i).adc();
MonitorElement* meShapeMap = 0;
if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
if ( waveLength[ism-1] == 0 ) meShapeMap = meShapeMapL1_[ism-1];
if ( waveLength[ism-1] == 1 ) meShapeMap = meShapeMapL2_[ism-1];
if ( waveLength[ism-1] == 2 ) meShapeMap = meShapeMapL3_[ism-1];
if ( waveLength[ism-1] == 3 ) meShapeMap = meShapeMapL4_[ism-1];
} else {
edm::LogWarning("EELaserTask") << " RtHalf = " << rtHalf[ism-1];
}
float xval = float(adc);
if ( meShapeMap ) meShapeMap->Fill(ic - 0.5, i + 0.5, xval);
}
NumbersPn::getPNs( ism, ix, iy, PNs );
for (unsigned int i=0; i<PNs.size(); i++) {
int ipn = PNs[i];
if ( ipn >= 0 && ipn < 80 ) numPN[ipn] = true;
}
}
} else {
edm::LogWarning("EELaserTask") << EEDigiCollection_ << " not available";
}
edm::Handle<EcalPnDiodeDigiCollection> pns;
if ( e.getByLabel(EcalPnDiodeDigiCollection_, pns) ) {
int nep = pns->size();
LogDebug("EELaserTask") << "event " << ievt_ << " pns collection size " << nep;
for ( EcalPnDiodeDigiCollection::const_iterator pnItr = pns->begin(); pnItr != pns->end(); ++pnItr ) {
if ( Numbers::subDet( pnItr->id() ) != EcalEndcap ) continue;
int ism = Numbers::iSM( pnItr->id() );
int num = pnItr->id().iPnId();
if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
int ipn = NumbersPn::ipnEE( ism, num );
if ( ipn >= 0 && ipn < 80 && numPN[ipn] == false ) continue;
float xvalped = 0.;
for (int i = 0; i < 4; i++) {
int adc = pnItr->sample(i).adc();
MonitorElement* mePNPed = 0;
if ( pnItr->sample(i).gainId() == 0 ) {
if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG01L1_[ism-1];
if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG01L2_[ism-1];
if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG01L3_[ism-1];
if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG01L4_[ism-1];
}
if ( pnItr->sample(i).gainId() == 1 ) {
if ( waveLength[ism-1] == 0 ) mePNPed = mePnPedMapG16L1_[ism-1];
if ( waveLength[ism-1] == 1 ) mePNPed = mePnPedMapG16L2_[ism-1];
if ( waveLength[ism-1] == 2 ) mePNPed = mePnPedMapG16L3_[ism-1];
if ( waveLength[ism-1] == 3 ) mePNPed = mePnPedMapG16L4_[ism-1];
}
float xval = float(adc);
if ( mePNPed ) mePNPed->Fill(num - 0.5, xval);
xvalped = xvalped + xval;
}
xvalped = xvalped / 4;
float xvalmax = 0.;
MonitorElement* mePN = 0;
for (int i = 0; i < 50; i++) {
int adc = pnItr->sample(i).adc();
float xval = float(adc);
if ( xval >= xvalmax ) xvalmax = xval;
}
xvalmax = xvalmax - xvalped;
if ( pnItr->sample(0).gainId() == 0 ) {
if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG01L1_[ism-1];
if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG01L2_[ism-1];
if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG01L3_[ism-1];
if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG01L4_[ism-1];
}
if ( pnItr->sample(0).gainId() == 1 ) {
if ( waveLength[ism-1] == 0 ) mePN = mePnAmplMapG16L1_[ism-1];
if ( waveLength[ism-1] == 1 ) mePN = mePnAmplMapG16L2_[ism-1];
if ( waveLength[ism-1] == 2 ) mePN = mePnAmplMapG16L3_[ism-1];
if ( waveLength[ism-1] == 3 ) mePN = mePnAmplMapG16L4_[ism-1];
}
if ( mePN ) mePN->Fill(num - 0.5, xvalmax);
if ( ipn >= 0 && ipn < 80 ) adcPN[ipn] = xvalmax;
}
} else {
edm::LogWarning("EELaserTask") << EcalPnDiodeDigiCollection_ << " not available";
}
edm::Handle<EcalUncalibratedRecHitCollection> hits;
if ( e.getByLabel(EcalUncalibratedRecHitCollection_, hits) ) {
int neh = hits->size();
LogDebug("EELaserTask") << "event " << ievt_ << " hits collection size " << neh;
for ( EcalUncalibratedRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
EEDetId id = hitItr->id();
int ix = id.ix();
int iy = id.iy();
int iz;
int ism = Numbers::iSM( id );
if ( ism >= 1 && ism <= 9 ){
ix = 101 - ix;
iz = 0;
}else{
iz = 1;
}
float xix = ix - 0.5;
float xiy = iy - 0.5;
if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::LASER_STD ||
runType[ism-1] == EcalDCCHeaderBlock::LASER_GAP ) ) continue;
if ( rtHalf[ism-1] != Numbers::RtHalf(id) ) continue;
MonitorElement* meAmplMap = 0;
MonitorElement* meTimeMap = 0;
MonitorElement* meAmplPNMap = 0;
MonitorElement* meAmplSummaryMap = 0;
if ( rtHalf[ism-1] == 0 || rtHalf[ism-1] == 1 ) {
if ( waveLength[ism-1] == 0 ) {
meAmplMap = meAmplMapL1_[ism-1];
meTimeMap = meTimeMapL1_[ism-1];
meAmplPNMap = meAmplPNMapL1_[ism-1];
meAmplSummaryMap = meAmplSummaryMapL1_[iz];
}
if ( waveLength[ism-1] == 1 ) {
meAmplMap = meAmplMapL2_[ism-1];
meTimeMap = meTimeMapL2_[ism-1];
meAmplPNMap = meAmplPNMapL2_[ism-1];
meAmplSummaryMap = meAmplSummaryMapL2_[iz];
}
if ( waveLength[ism-1] == 2 ) {
meAmplMap = meAmplMapL3_[ism-1];
meTimeMap = meTimeMapL3_[ism-1];
meAmplPNMap = meAmplPNMapL3_[ism-1];
meAmplSummaryMap = meAmplSummaryMapL3_[iz];
}
if ( waveLength[ism-1] == 3 ) {
meAmplMap = meAmplMapL4_[ism-1];
meTimeMap = meTimeMapL4_[ism-1];
meAmplPNMap = meAmplPNMapL4_[ism-1];
meAmplSummaryMap = meAmplSummaryMapL4_[iz];
}
} else {
edm::LogWarning("EELaserTask") << " RtHalf = " << rtHalf[ism-1];
}
float xval = hitItr->amplitude();
if ( xval <= 0. ) xval = 0.0;
float yval = hitItr->jitter() + 5.0;
if ( yval <= 0. ) yval = 0.0;
float zval = hitItr->pedestal();
if ( zval <= 0. ) zval = 0.0;
if ( meAmplMap ) meAmplMap->Fill(xix, xiy, xval);
if ( xval > 16. ) {
if ( meTimeMap ) meTimeMap->Fill(xix, xiy, yval);
}
float wval = 0.;
NumbersPn::getPNs( ism, ix, iy, PNs );
if ( PNs.size() > 0 ) {
int ipn = PNs[0];
if ( ipn >= 0 && ipn < 80 ) {
if ( adcPN[ipn] != 0. ) wval = xval / adcPN[ipn];
}
}
if ( meAmplPNMap ) meAmplPNMap->Fill(xix, xiy, wval);
if( meAmplSummaryMap ) meAmplSummaryMap->Fill(xix, xiy, xval);
}
} else {
edm::LogWarning("EELaserTask") << EcalUncalibratedRecHitCollection_ << " not available";
}
}
| void EELaserTask::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 106 of file EELaserTask.cc.
References dqmStore_, ievt_, prefixME_, DQMStore::rmdir(), and DQMStore::setCurrentFolder().
| void EELaserTask::beginRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
BeginRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 117 of file EELaserTask.cc.
References Numbers::initGeometry(), mergeRuns_, and reset().
{
Numbers::initGeometry(c, false);
if ( ! mergeRuns_ ) this->reset();
}
| void EELaserTask::cleanup | ( | void | ) | [protected] |
Cleanup.
Definition at line 542 of file EELaserTask.cc.
References dqmStore_, spr::find(), edm::getName(), i, init_, laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().
Referenced by endJob().
{
if ( ! init_ ) return;
if ( dqmStore_ ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1");
for (int i = 0; i < 18; i++) {
if ( meShapeMapL1_[i] ) dqmStore_->removeElement( meShapeMapL1_[i]->getName() );
meShapeMapL1_[i] = 0;
if ( meAmplMapL1_[i] ) dqmStore_->removeElement( meAmplMapL1_[i]->getName() );
meAmplMapL1_[i] = 0;
if ( meTimeMapL1_[i] ) dqmStore_->removeElement( meTimeMapL1_[i]->getName() );
meTimeMapL1_[i] = 0;
if ( meAmplPNMapL1_[i] ) dqmStore_->removeElement( meAmplPNMapL1_[i]->getName() );
meAmplPNMapL1_[i] = 0;
}
for(int i=0; i<2; i++)
if( meAmplSummaryMapL1_[i] ) dqmStore_->removeElement( meAmplSummaryMapL1_[i]->getName() );
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2");
for (int i = 0; i < 18; i++) {
if ( meShapeMapL2_[i] ) dqmStore_->removeElement( meShapeMapL2_[i]->getName() );
meShapeMapL2_[i] = 0;
if ( meAmplMapL2_[i] ) dqmStore_->removeElement( meAmplMapL2_[i]->getName() );
meAmplMapL2_[i] = 0;
if ( meTimeMapL2_[i] ) dqmStore_->removeElement( meTimeMapL2_[i]->getName() );
meTimeMapL2_[i] = 0;
if ( meAmplPNMapL2_[i] ) dqmStore_->removeElement( meAmplPNMapL2_[i]->getName() );
meAmplPNMapL2_[i] = 0;
}
for(int i=0; i<2; i++)
if( meAmplSummaryMapL2_[i] ) dqmStore_->removeElement( meAmplSummaryMapL2_[i]->getName() );
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3");
for (int i = 0; i < 18; i++) {
if ( meShapeMapL3_[i] ) dqmStore_->removeElement( meShapeMapL3_[i]->getName() );
meShapeMapL3_[i] = 0;
if ( meAmplMapL3_[i] ) dqmStore_->removeElement( meAmplMapL3_[i]->getName() );
meAmplMapL3_[i] = 0;
if ( meTimeMapL3_[i] ) dqmStore_->removeElement( meTimeMapL3_[i]->getName() );
meTimeMapL3_[i] = 0;
if ( meAmplPNMapL3_[i] ) dqmStore_->removeElement( meAmplPNMapL3_[i]->getName() );
meAmplPNMapL3_[i] = 0;
}
for(int i=0; i<2; i++)
if( meAmplSummaryMapL3_[i] ) dqmStore_->removeElement( meAmplSummaryMapL3_[i]->getName() );
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4");
for (int i = 0; i < 18; i++) {
if ( meShapeMapL4_[i] ) dqmStore_->removeElement( meShapeMapL4_[i]->getName() );
meShapeMapL4_[i] = 0;
if ( meAmplMapL4_[i] ) dqmStore_->removeElement( meAmplMapL4_[i]->getName() );
meAmplMapL4_[i] = 0;
if ( meTimeMapL4_[i] ) dqmStore_->removeElement( meTimeMapL4_[i]->getName() );
meTimeMapL4_[i] = 0;
if ( meAmplPNMapL4_[i] ) dqmStore_->removeElement( meAmplPNMapL4_[i]->getName() );
meAmplPNMapL4_[i] = 0;
}
for(int i=0; i<2; i++)
if( meAmplSummaryMapL4_[i] ) dqmStore_->removeElement( meAmplSummaryMapL4_[i]->getName() );
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN");
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN/Gain01");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG01L1_[i] ) dqmStore_->removeElement( mePnAmplMapG01L1_[i]->getName() );
mePnAmplMapG01L1_[i] = 0;
if ( mePnPedMapG01L1_[i] ) dqmStore_->removeElement( mePnPedMapG01L1_[i]->getName() );
mePnPedMapG01L1_[i] = 0;
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser1/PN/Gain16");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG16L1_[i] ) dqmStore_->removeElement( mePnAmplMapG16L1_[i]->getName() );
mePnAmplMapG16L1_[i] = 0;
if ( mePnPedMapG16L1_[i] ) dqmStore_->removeElement( mePnPedMapG16L1_[i]->getName() );
mePnPedMapG16L1_[i] = 0;
}
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN");
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN/Gain01");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG01L2_[i] ) dqmStore_->removeElement( mePnAmplMapG01L2_[i]->getName() );
mePnAmplMapG01L2_[i] = 0;
if ( mePnPedMapG01L2_[i] ) dqmStore_->removeElement( mePnPedMapG01L2_[i]->getName() );
mePnPedMapG01L2_[i] = 0;
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser2/PN/Gain16");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG16L2_[i] ) dqmStore_->removeElement( mePnAmplMapG16L2_[i]->getName() );
mePnAmplMapG16L2_[i] = 0;
if ( mePnPedMapG16L2_[i] ) dqmStore_->removeElement( mePnPedMapG16L2_[i]->getName() );
mePnPedMapG16L2_[i] = 0;
}
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN");
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN/Gain01");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG01L3_[i] ) dqmStore_->removeElement( mePnAmplMapG01L3_[i]->getName() );
mePnAmplMapG01L3_[i] = 0;
if ( mePnPedMapG01L3_[i] ) dqmStore_->removeElement( mePnPedMapG01L3_[i]->getName() );
mePnPedMapG01L3_[i] = 0;
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser3/PN/Gain16");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG16L3_[i] ) dqmStore_->removeElement( mePnAmplMapG16L3_[i]->getName() );
mePnAmplMapG16L3_[i] = 0;
if ( mePnPedMapG16L3_[i] ) dqmStore_->removeElement( mePnPedMapG16L3_[i]->getName() );
mePnPedMapG16L3_[i] = 0;
}
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN");
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN/Gain01");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG01L4_[i] ) dqmStore_->removeElement( mePnAmplMapG01L4_[i]->getName() );
mePnAmplMapG01L4_[i] = 0;
if ( mePnPedMapG01L4_[i] ) dqmStore_->removeElement( mePnPedMapG01L4_[i]->getName() );
mePnPedMapG01L4_[i] = 0;
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/Laser4/PN/Gain16");
for (int i = 0; i < 18; i++) {
if ( mePnAmplMapG16L4_[i] ) dqmStore_->removeElement( mePnAmplMapG16L4_[i]->getName() );
mePnAmplMapG16L4_[i] = 0;
if ( mePnPedMapG16L4_[i] ) dqmStore_->removeElement( mePnPedMapG16L4_[i]->getName() );
mePnPedMapG16L4_[i] = 0;
}
}
}
init_ = false;
}
| void EELaserTask::endJob | ( | void | ) | [protected, virtual] |
EndJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 718 of file EELaserTask.cc.
References cleanup(), enableCleanup_, and ievt_.
{
edm::LogInfo("EELaserTask") << "analyzed " << ievt_ << " events";
if ( enableCleanup_ ) this->cleanup();
}
| void EELaserTask::endRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
| void EELaserTask::reset | ( | void | ) | [protected] |
Reset.
Definition at line 129 of file EELaserTask.cc.
References spr::find(), i, laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, and MonitorElement::Reset().
Referenced by beginRun().
{
for (int i = 0; i < 18; i++) {
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
if ( meShapeMapL1_[i] ) meShapeMapL1_[i]->Reset();
if ( meAmplMapL1_[i] ) meAmplMapL1_[i]->Reset();
if ( meTimeMapL1_[i] ) meTimeMapL1_[i]->Reset();
if ( meAmplPNMapL1_[i] ) meAmplPNMapL1_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
if ( meShapeMapL2_[i] ) meShapeMapL2_[i]->Reset();
if ( meAmplMapL2_[i] ) meAmplMapL2_[i]->Reset();
if ( meTimeMapL2_[i] ) meTimeMapL2_[i]->Reset();
if ( meAmplPNMapL2_[i] ) meAmplPNMapL2_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
if ( meShapeMapL3_[i] ) meShapeMapL3_[i]->Reset();
if ( meAmplMapL3_[i] ) meAmplMapL3_[i]->Reset();
if ( meTimeMapL3_[i] ) meTimeMapL3_[i]->Reset();
if ( meAmplPNMapL3_[i] ) meAmplPNMapL3_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
if ( meShapeMapL4_[i] ) meShapeMapL4_[i]->Reset();
if ( meAmplMapL4_[i] ) meAmplMapL4_[i]->Reset();
if ( meTimeMapL4_[i] ) meTimeMapL4_[i]->Reset();
if ( meAmplPNMapL4_[i] ) meAmplPNMapL4_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
if ( mePnAmplMapG01L1_[i] ) mePnAmplMapG01L1_[i]->Reset();
if ( mePnPedMapG01L1_[i] ) mePnPedMapG01L1_[i]->Reset();
if ( mePnAmplMapG16L1_[i] ) mePnAmplMapG16L1_[i]->Reset();
if ( mePnPedMapG16L1_[i] ) mePnPedMapG16L1_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
if ( mePnAmplMapG01L2_[i] ) mePnAmplMapG01L2_[i]->Reset();
if ( mePnPedMapG01L2_[i] ) mePnPedMapG01L2_[i]->Reset();
if ( mePnAmplMapG16L2_[i] ) mePnAmplMapG16L2_[i]->Reset();
if ( mePnPedMapG16L2_[i] ) mePnPedMapG16L2_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
if ( mePnAmplMapG01L3_[i] ) mePnAmplMapG01L3_[i]->Reset();
if ( mePnPedMapG01L3_[i] ) mePnPedMapG01L3_[i]->Reset();
if ( mePnAmplMapG16L3_[i] ) mePnAmplMapG16L3_[i]->Reset();
if ( mePnPedMapG16L3_[i] ) mePnPedMapG16L3_[i]->Reset();
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
if ( mePnAmplMapG01L4_[i] ) mePnAmplMapG01L4_[i]->Reset();
if ( mePnPedMapG01L4_[i] ) mePnPedMapG01L4_[i]->Reset();
if ( mePnAmplMapG16L4_[i] ) mePnAmplMapG16L4_[i]->Reset();
if ( mePnPedMapG16L4_[i] ) mePnPedMapG16L4_[i]->Reset();
}
}
for(int i=0; i<2; i++){
if( meAmplSummaryMapL1_[i] ) meAmplSummaryMapL1_[i]->Reset();
if( meAmplSummaryMapL2_[i] ) meAmplSummaryMapL2_[i]->Reset();
if( meAmplSummaryMapL3_[i] ) meAmplSummaryMapL3_[i]->Reset();
if( meAmplSummaryMapL4_[i] ) meAmplSummaryMapL4_[i]->Reset();
}
}
| void EELaserTask::setup | ( | void | ) | [protected] |
Setup.
Definition at line 204 of file EELaserTask.cc.
References DQMStore::bookProfile(), DQMStore::bookProfile2D(), dqmStore_, spr::find(), i, init_, Numbers::ix0EE(), Numbers::iy0EE(), laserWavelengths_, meAmplMapL1_, meAmplMapL2_, meAmplMapL3_, meAmplMapL4_, meAmplPNMapL1_, meAmplPNMapL2_, meAmplPNMapL3_, meAmplPNMapL4_, meAmplSummaryMapL1_, meAmplSummaryMapL2_, meAmplSummaryMapL3_, meAmplSummaryMapL4_, mePnAmplMapG01L1_, mePnAmplMapG01L2_, mePnAmplMapG01L3_, mePnAmplMapG01L4_, mePnAmplMapG16L1_, mePnAmplMapG16L2_, mePnAmplMapG16L3_, mePnAmplMapG16L4_, mePnPedMapG01L1_, mePnPedMapG01L2_, mePnPedMapG01L3_, mePnPedMapG01L4_, mePnPedMapG16L1_, mePnPedMapG16L2_, mePnPedMapG16L3_, mePnPedMapG16L4_, meShapeMapL1_, meShapeMapL2_, meShapeMapL3_, meShapeMapL4_, meTimeMapL1_, meTimeMapL2_, meTimeMapL3_, meTimeMapL4_, mergeVDriftHistosByStation::name, prefixME_, Numbers::sEE(), MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), and DQMStore::tag().
Referenced by analyze().
{
init_ = true;
std::string name;
std::stringstream LaserN, LN;
if ( dqmStore_ ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask");
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 1) != laserWavelengths_.end() ) {
LaserN.str("");
LaserN << "Laser" << 1;
LN.str("");
LN << "L" << 1;
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
for (int i = 0; i < 18; i++) {
name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
meShapeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
meShapeMapL1_[i]->setAxisTitle("channel", 1);
meShapeMapL1_[i]->setAxisTitle("sample", 2);
meShapeMapL1_[i]->setAxisTitle("amplitude", 3);
dqmStore_->tag(meShapeMapL1_[i], i+1);
name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
meAmplMapL1_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplMapL1_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL1_[i]->setAxisTitle("101-ix", 1);
meAmplMapL1_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplMapL1_[i], i+1);
name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
meTimeMapL1_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
meTimeMapL1_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL1_[i]->setAxisTitle("101-ix", 1);
meTimeMapL1_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meTimeMapL1_[i], i+1);
name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
meAmplPNMapL1_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplPNMapL1_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL1_[i]->setAxisTitle("101-ix", 1);
meAmplPNMapL1_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplPNMapL1_[i], i+1);
}
name = "EELT amplitude map " + LN.str() + " EE -";
meAmplSummaryMapL1_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL1_[0]->setAxisTitle("ix", 1);
meAmplSummaryMapL1_[0]->setAxisTitle("iy", 2);
name = "EELT amplitude map " + LN.str() + " EE +";
meAmplSummaryMapL1_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL1_[1]->setAxisTitle("ix", 1);
meAmplSummaryMapL1_[1]->setAxisTitle("iy", 2);
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN");
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain01");
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnAmplMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG01L1_[i]->setAxisTitle("channel", 1);
mePnAmplMapG01L1_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG01L1_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnPedMapG01L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG01L1_[i]->setAxisTitle("channel", 1);
mePnPedMapG01L1_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG01L1_[i], i+1);
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnAmplMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG16L1_[i]->setAxisTitle("channel", 1);
mePnAmplMapG16L1_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG16L1_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnPedMapG16L1_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG16L1_[i]->setAxisTitle("channel", 1);
mePnPedMapG16L1_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG16L1_[i], i+1);
}
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 2) != laserWavelengths_.end() ) {
LaserN.str("");
LaserN << "Laser" << 2;
LN.str("");
LN << "L" << 2;
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
for (int i = 0; i < 18; i++) {
name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
meShapeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
meShapeMapL2_[i]->setAxisTitle("channel", 1);
meShapeMapL2_[i]->setAxisTitle("sample", 2);
meShapeMapL2_[i]->setAxisTitle("amplitude", 3);
dqmStore_->tag(meShapeMapL2_[i], i+1);
name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
meAmplMapL2_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplMapL2_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL2_[i]->setAxisTitle("101-ix", 1);
meAmplMapL2_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplMapL2_[i], i+1);
name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
meTimeMapL2_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
meTimeMapL2_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL2_[i]->setAxisTitle("101-ix", 1);
meTimeMapL2_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meTimeMapL2_[i], i+1);
name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
meAmplPNMapL2_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplPNMapL2_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL2_[i]->setAxisTitle("101-ix", 1);
meAmplPNMapL2_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplPNMapL2_[i], i+1);
}
name = "EELT amplitude map " + LN.str() + " EE -";
meAmplSummaryMapL2_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL2_[0]->setAxisTitle("ix", 1);
meAmplSummaryMapL2_[0]->setAxisTitle("iy", 2);
name = "EELT amplitude map " + LN.str() + " EE +";
meAmplSummaryMapL2_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL2_[1]->setAxisTitle("ix", 1);
meAmplSummaryMapL2_[1]->setAxisTitle("iy", 2);
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnAmplMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG01L2_[i]->setAxisTitle("channel", 1);
mePnAmplMapG01L2_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG01L2_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnPedMapG01L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG01L2_[i]->setAxisTitle("channel", 1);
mePnPedMapG01L2_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG01L2_[i], i+1);
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnAmplMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG16L2_[i]->setAxisTitle("channel", 1);
mePnAmplMapG16L2_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG16L2_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnPedMapG16L2_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG16L2_[i]->setAxisTitle("channel", 1);
mePnPedMapG16L2_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG16L2_[i], i+1);
}
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 3) != laserWavelengths_.end() ) {
LaserN.str("");
LaserN << "Laser" << 3;
LN.str("");
LN << "L" << 3;
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
for (int i = 0; i < 18; i++) {
name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
meShapeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
meShapeMapL3_[i]->setAxisTitle("channel", 1);
meShapeMapL3_[i]->setAxisTitle("sample", 2);
meShapeMapL3_[i]->setAxisTitle("amplitude", 3);
dqmStore_->tag(meShapeMapL3_[i], i+1);
name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
meAmplMapL3_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplMapL3_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL3_[i]->setAxisTitle("101-ix", 1);
meAmplMapL3_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplMapL3_[i], i+1);
name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
meTimeMapL3_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
meTimeMapL3_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL3_[i]->setAxisTitle("101-ix", 1);
meTimeMapL3_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meTimeMapL3_[i], i+1);
name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
meAmplPNMapL3_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplPNMapL3_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL3_[i]->setAxisTitle("101-ix", 1);
meAmplPNMapL3_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplPNMapL3_[i], i+1);
}
name = "EELT amplitude map " + LN.str() + " EE -";
meAmplSummaryMapL3_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL3_[0]->setAxisTitle("ix", 1);
meAmplSummaryMapL3_[0]->setAxisTitle("iy", 2);
name = "EELT amplitude map " + LN.str() + " EE +";
meAmplSummaryMapL3_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL3_[1]->setAxisTitle("ix", 1);
meAmplSummaryMapL3_[1]->setAxisTitle("iy", 2);
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnAmplMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG01L3_[i]->setAxisTitle("channel", 1);
mePnAmplMapG01L3_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG01L3_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnPedMapG01L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG01L3_[i]->setAxisTitle("channel", 1);
mePnPedMapG01L3_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG01L3_[i], i+1);
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnAmplMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG16L3_[i]->setAxisTitle("channel", 1);
mePnAmplMapG16L3_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG16L3_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnPedMapG16L3_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG16L3_[i]->setAxisTitle("channel", 1);
mePnPedMapG16L3_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG16L3_[i], i+1);
}
}
if ( find(laserWavelengths_.begin(), laserWavelengths_.end(), 4) != laserWavelengths_.end() ) {
LaserN.str("");
LaserN << "Laser" << 4;
LN.str("");
LN << "L" << 4;
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str());
for (int i = 0; i < 18; i++) {
name = "EELT shape " + Numbers::sEE(i+1) + " " + LN.str();
meShapeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 850, 0., 850., 10, 0., 10., 4096, 0., 4096., "s");
meShapeMapL4_[i]->setAxisTitle("channel", 1);
meShapeMapL4_[i]->setAxisTitle("sample", 2);
meShapeMapL4_[i]->setAxisTitle("amplitude", 3);
dqmStore_->tag(meShapeMapL4_[i], i+1);
name = "EELT amplitude " + Numbers::sEE(i+1) + " " + LN.str();
meAmplMapL4_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplMapL4_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplMapL4_[i]->setAxisTitle("101-ix", 1);
meAmplMapL4_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplMapL4_[i], i+1);
name = "EELT timing " + Numbers::sEE(i+1) + " " + LN.str();
meTimeMapL4_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 250, 0., 10., "s");
meTimeMapL4_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meTimeMapL4_[i]->setAxisTitle("101-ix", 1);
meTimeMapL4_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meTimeMapL4_[i], i+1);
name = "EELT amplitude over PN " + Numbers::sEE(i+1) + " " + LN.str();
meAmplPNMapL4_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., 4096, 0., 4096.*12., "s");
meAmplPNMapL4_[i]->setAxisTitle("ix", 1);
if ( i+1 >= 1 && i+1 <= 9 ) meAmplPNMapL4_[i]->setAxisTitle("101-ix", 1);
meAmplPNMapL4_[i]->setAxisTitle("iy", 2);
dqmStore_->tag(meAmplPNMapL4_[i], i+1);
}
name = "EELT amplitude map " + LN.str() + " EE -";
meAmplSummaryMapL4_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL4_[0]->setAxisTitle("ix", 1);
meAmplSummaryMapL4_[0]->setAxisTitle("iy", 2);
name = "EELT amplitude map " + LN.str() + " EE +";
meAmplSummaryMapL4_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., 0., 4096.);
meAmplSummaryMapL4_[1]->setAxisTitle("ix", 1);
meAmplSummaryMapL4_[1]->setAxisTitle("iy", 2);
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnAmplMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG01L4_[i]->setAxisTitle("channel", 1);
mePnAmplMapG01L4_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG01L4_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G01 " + LN.str();
mePnPedMapG01L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG01L4_[i]->setAxisTitle("channel", 1);
mePnPedMapG01L4_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG01L4_[i], i+1);
}
dqmStore_->setCurrentFolder(prefixME_ + "/EELaserTask/" + LaserN.str() + "/PN/Gain16");
for (int i = 0; i < 18; i++) {
name = "EELT PNs amplitude " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnAmplMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnAmplMapG16L4_[i]->setAxisTitle("channel", 1);
mePnAmplMapG16L4_[i]->setAxisTitle("amplitude", 2);
dqmStore_->tag(mePnAmplMapG16L4_[i], i+1);
name = "EELT PNs pedestal " + Numbers::sEE(i+1) + " G16 " + LN.str();
mePnPedMapG16L4_[i] = dqmStore_->bookProfile(name, name, 10, 0., 10., 4096, 0., 4096., "s");
mePnPedMapG16L4_[i]->setAxisTitle("channel", 1);
mePnPedMapG16L4_[i]->setAxisTitle("pedestal", 2);
dqmStore_->tag(mePnPedMapG16L4_[i], i+1);
}
}
}
}
DQMStore* EELaserTask::dqmStore_ [private] |
Definition at line 60 of file EELaserTask.h.
Referenced by beginJob(), cleanup(), EELaserTask(), and setup().
Definition at line 70 of file EELaserTask.h.
Referenced by analyze(), and EELaserTask().
Definition at line 68 of file EELaserTask.h.
Referenced by analyze(), and EELaserTask().
Definition at line 71 of file EELaserTask.h.
Referenced by analyze(), and EELaserTask().
edm::InputTag EELaserTask::EEDigiCollection_ [private] |
Definition at line 69 of file EELaserTask.h.
Referenced by analyze(), and EELaserTask().
bool EELaserTask::enableCleanup_ [private] |
Definition at line 64 of file EELaserTask.h.
Referenced by EELaserTask(), and endJob().
int EELaserTask::ievt_ [private] |
Definition at line 58 of file EELaserTask.h.
Referenced by analyze(), beginJob(), and endJob().
bool EELaserTask::init_ [private] |
Definition at line 115 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), and setup().
std::vector<int> EELaserTask::laserWavelengths_ [private] |
Definition at line 72 of file EELaserTask.h.
Referenced by cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplMapL1_[18] [private] |
Definition at line 75 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplMapL2_[18] [private] |
Definition at line 84 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplMapL3_[18] [private] |
Definition at line 93 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplMapL4_[18] [private] |
Definition at line 102 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplPNMapL1_[18] [private] |
Definition at line 77 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplPNMapL2_[18] [private] |
Definition at line 86 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplPNMapL3_[18] [private] |
Definition at line 95 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplPNMapL4_[18] [private] |
Definition at line 104 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplSummaryMapL1_[2] [private] |
Definition at line 110 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplSummaryMapL2_[2] [private] |
Definition at line 111 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplSummaryMapL3_[2] [private] |
Definition at line 112 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meAmplSummaryMapL4_[2] [private] |
Definition at line 113 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG01L1_[18] [private] |
Definition at line 78 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG01L2_[18] [private] |
Definition at line 87 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG01L3_[18] [private] |
Definition at line 96 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG01L4_[18] [private] |
Definition at line 105 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG16L1_[18] [private] |
Definition at line 80 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG16L2_[18] [private] |
Definition at line 89 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG16L3_[18] [private] |
Definition at line 98 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnAmplMapG16L4_[18] [private] |
Definition at line 107 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG01L1_[18] [private] |
Definition at line 79 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG01L2_[18] [private] |
Definition at line 88 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG01L3_[18] [private] |
Definition at line 97 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG01L4_[18] [private] |
Definition at line 106 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG16L1_[18] [private] |
Definition at line 81 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG16L2_[18] [private] |
Definition at line 90 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG16L3_[18] [private] |
Definition at line 99 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::mePnPedMapG16L4_[18] [private] |
Definition at line 108 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
bool EELaserTask::mergeRuns_ [private] |
Definition at line 66 of file EELaserTask.h.
Referenced by beginRun(), and EELaserTask().
MonitorElement* EELaserTask::meShapeMapL1_[18] [private] |
Definition at line 74 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meShapeMapL2_[18] [private] |
Definition at line 83 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meShapeMapL3_[18] [private] |
Definition at line 92 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meShapeMapL4_[18] [private] |
Definition at line 101 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meTimeMapL1_[18] [private] |
Definition at line 76 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meTimeMapL2_[18] [private] |
Definition at line 85 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meTimeMapL3_[18] [private] |
Definition at line 94 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
MonitorElement* EELaserTask::meTimeMapL4_[18] [private] |
Definition at line 103 of file EELaserTask.h.
Referenced by analyze(), cleanup(), EELaserTask(), reset(), and setup().
std::string EELaserTask::prefixME_ [private] |
Definition at line 62 of file EELaserTask.h.
Referenced by beginJob(), cleanup(), EELaserTask(), and setup().
1.7.3