|
|
|
#include <EBBeamHodoTask.h>
Definition at line 21 of file EBBeamHodoTask.h.
| EBBeamHodoTask::EBBeamHodoTask | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 33 of file EBBeamHodoTask.cc.
References cryInBeam_, dqmStore_, EcalRawDataCollection_, EcalTBEventHeader_, EcalTBHodoscopeRawInfo_, EcalTBHodoscopeRecInfo_, EcalTBTDCRawInfo_, EcalTBTDCRecInfo_, EcalUncalibratedRecHitCollection_, enableCleanup_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, init_, meCaloVsHodoXPos_, meCaloVsHodoYPos_, meCaloVsTDCTime_, meEvsXRecHis_, meEvsXRecProf_, meEvsYRecHis_, meEvsYRecProf_, meHodoOcc_, meHodoPosRecX_, meHodoPosRecXY_, meHodoPosRecY_, meHodoPosXMinusCaloPosXVsCry_, meHodoPosYMinusCaloPosYVsCry_, meHodoQuaXRec_, meHodoQuaYRec_, meHodoRaw_, meHodoSloXRec_, meHodoSloYRec_, meMissingCollections_, mergeRuns_, meTDCRec_, meTDCTimeMinusCaloTimeVsCry_, cmsCodeRules::cppFunctionSkipper::operator, prefixME_, previousCryInBeam_, and tableIsMoving_.
{
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);
EcalTBEventHeader_ = ps.getParameter<edm::InputTag>("EcalTBEventHeader");
EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
EcalUncalibratedRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection");
EcalTBTDCRawInfo_ = ps.getParameter<edm::InputTag>("EcalTBTDCRawInfo");
EcalTBHodoscopeRawInfo_ = ps.getParameter<edm::InputTag>("EcalTBHodoscopeRawInfo");
EcalTBTDCRecInfo_ = ps.getParameter<edm::InputTag>("EcalTBTDCRecInfo");
EcalTBHodoscopeRecInfo_ = ps.getParameter<edm::InputTag>("EcalTBHodoscopeRecInfo");
tableIsMoving_ = false;
cryInBeam_ =0;
previousCryInBeam_ = -99999;
// to be filled all the time
for (int i=0; i<4; i++) {
meHodoOcc_[i] =0;
meHodoRaw_[i] =0;
}
meTDCRec_ =0;
// filled only when: the table does not move
meHodoPosRecX_ =0;
meHodoPosRecY_ =0;
meHodoPosRecXY_ =0;
meHodoSloXRec_ =0;
meHodoSloYRec_ =0;
meHodoQuaXRec_ =0;
meHodoQuaYRec_ =0;
meHodoPosXMinusCaloPosXVsCry_ =0;
meHodoPosYMinusCaloPosYVsCry_ =0;
meTDCTimeMinusCaloTimeVsCry_ =0;
meMissingCollections_ =0;
meEvsXRecProf_ =0;
meEvsYRecProf_ =0;
meEvsXRecHis_ =0;
meEvsYRecHis_ =0;
// and matrix 5x5 available
meCaloVsHodoXPos_ =0;
meCaloVsHodoYPos_ =0;
meCaloVsTDCTime_ =0;
}
| EBBeamHodoTask::~EBBeamHodoTask | ( | ) | [virtual] |
| void EBBeamHodoTask::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Analyze.
Implements edm::EDAnalyzer.
Definition at line 334 of file EBBeamHodoTask.cc.
References EcalDCCHeaderBlock::BEAMH2, EcalDCCHeaderBlock::BEAMH4, EcalTBHodoscopePlaneRawHits::channels(), EcalTBHodoscopeRawInfo::channels(), cryInBeam_, cryInBeamCounter_, EcalTBEventHeader::crystalInBeam(), EcalBarrel, EcalRawDataCollection_, EcalTBEventHeader_, EcalTBHodoscopeRawInfo_, EcalTBHodoscopeRecInfo_, EcalTBTDCRawInfo_, EcalTBTDCRecInfo_, EcalUncalibratedRecHitCollection_, EBDetId::ETAPHIMODE, EcalTBEventHeader::eventNumber(), MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::Event::getByLabel(), EcalTBHodoscopeRawInfo::getPlaneRawHits(), i, EBDetId::ieta(), ievt_, init_, EcalTBHodoscopePlaneRawHits::isChannelFired(), funct::log(), LogDebug, LV1_, meCaloVsHodoXPos_, meCaloVsHodoYPos_, meCaloVsTDCTime_, meEvsXRecHis_, meEvsXRecProf_, meEvsYRecHis_, meEvsYRecProf_, meHodoOcc_, meHodoPosRecX_, meHodoPosRecXY_, meHodoPosRecY_, meHodoPosXMinusCaloPosXVsCry_, meHodoPosYMinusCaloPosYVsCry_, meHodoQuaXRec_, meHodoQuaYRec_, meHodoRaw_, meHodoSloXRec_, meHodoSloYRec_, meMissingCollections_, meTDCRec_, meTDCTimeMinusCaloTimeVsCry_, EcalTBHodoscopePlaneRawHits::numberOfFiredHits(), EcalTBTDCRecInfo::offset(), EcalTBHodoscopeRawInfo::planes(), EcalTBHodoscopeRecInfo::posX(), EcalTBHodoscopeRecInfo::posY(), previousCryInBeam_, edm::Handle< T >::product(), EcalTBHodoscopeRecInfo::qualX(), EcalTBHodoscopeRecInfo::qualY(), MonitorElement::Reset(), resetNow_, setup(), EcalTBHodoscopeRecInfo::slopeX(), EcalTBHodoscopeRecInfo::slopeY(), Numbers::subDet(), EcalTBEventHeader::tableIsMoving(), tableIsMoving_, EBDetId::validDetId(), and CommonMethods::weight().
{
bool enable = false;
edm::Handle<EcalTBEventHeader> pHeader;
const EcalTBEventHeader* Header =0;
if ( e.getByLabel(EcalTBEventHeader_, pHeader) ) {
Header = pHeader.product(); // get a ptr to the product
if (!Header) {
edm::LogWarning("EBBeamHodoTask") << "Event header not found. Returning. ";
meMissingCollections_-> Fill(0); // bin1: missing CMSSW edm::Event header
return;
}
tableIsMoving_ = Header->tableIsMoving();
cryInBeam_ = Header->crystalInBeam(); // cryInBeam_ = Header->nominalCrystalInBeam();
if (previousCryInBeam_ == -99999 )
{ previousCryInBeam_ = cryInBeam_ ; }
LogDebug("EBBeamHodoTask") << "event: " << ievt_ << " event header found ";
if (tableIsMoving_){
LogDebug("EBBeamHodoTask") << "Table is moving. "; }
else
{ LogDebug("EBBeamHodoTask") << "Table is not moving. "; }
} else {
edm::LogWarning("EBBeamHodoTask") << "Event header not found (exception caught). Returning. ";
return;
}
edm::Handle<EcalRawDataCollection> dcchs;
if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
if ( Numbers::subDet( *dcchItr ) != EcalBarrel ) continue;
if ( dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH4 ||
dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH2 ) enable = true;
}
} else {
edm::LogWarning("EcalBeamTask") << EcalRawDataCollection_ << " not available";
meMissingCollections_-> Fill(1); // bin2: missing DCC headers
return;
// see bottom of cc file for compatibility to 2004 data [***]
}
if ( ! enable ) return;
if ( ! init_ ) this->setup();
ievt_++;
LV1_ = Header->eventNumber();
edm::Handle<EcalUncalibratedRecHitCollection> pUncalRH;
const EcalUncalibratedRecHitCollection* uncalRecH =0;
if ( e.getByLabel(EcalUncalibratedRecHitCollection_, pUncalRH) ) {
uncalRecH = pUncalRH.product(); // get a ptr to the product
int neh = pUncalRH->size();
LogDebug("EBBeamHodoTask") << EcalUncalibratedRecHitCollection_ << " found in event " << ievt_ << "; hits collection size " << neh;
} else {
edm::LogWarning("EBBeamHodoTask") << EcalUncalibratedRecHitCollection_ << " not available";
meMissingCollections_-> Fill(2); // bin3: missing uncalibRecHits
return;
}
edm::Handle<EcalTBTDCRawInfo> pTDCRaw;
const EcalTBTDCRawInfo* rawTDC=0;
if ( e.getByLabel(EcalTBTDCRawInfo_, pTDCRaw) ) {
rawTDC = pTDCRaw.product();
} else {
edm::LogError("EcalBeamTask") << "Error! Can't get the product EcalTBTDCRawInfo. Returning.";
meMissingCollections_-> Fill(4); // bin5: missing raw TDC
return;
}
edm::Handle<EcalTBHodoscopeRawInfo> pHodoRaw;
const EcalTBHodoscopeRawInfo* rawHodo=0;
if ( e.getByLabel(EcalTBHodoscopeRawInfo_, pHodoRaw) ) {
rawHodo = pHodoRaw.product();
if(rawHodo->planes() ){
LogDebug("EcalBeamTask") << "hodoscopeRaw: num planes: " << rawHodo->planes()
<< " channels in plane 1: " << rawHodo->channels(0);
}
} else {
edm::LogError("EcalBeamTask") << "Error! Can't get the product EcalTBHodoscopeRawInfo. Returning";
meMissingCollections_-> Fill(3); // bin4: missing raw hodo hits collection
return;
}
if ( !rawTDC ||!rawHodo || !uncalRecH || !( rawHodo->planes() )) {
edm::LogWarning("EcalBeamTask") << "analyze: missing a needed collection or hodo collection empty. Returning.";
return;
}
LogDebug("EBBeamHodoTask") << " TDC raw, Hodo raw, uncalRecH and DCCheader found.";
// table has come to a stop is identified by new value of cry_in_beam
// - increase counter of crystals that have been on beam
// - set flag for resetting
if (cryInBeam_ != previousCryInBeam_ ) {
previousCryInBeam_ = cryInBeam_ ;
cryInBeamCounter_++;
resetNow_ = true;
// since flag "tableIsMoving==false" is reliable (as we can tell, so far),
// operations due when "table has started moving"
// can be done after the change in crystal in beam
LogDebug("EcalBeamTask") << "At event number : " << LV1_ << " switching table status: from moving to still. "
<< " cry in beam is: " << cryInBeam_ << ", step being: " << cryInBeamCounter_ ;
// fill here plots which keep history of beamed crystals
float HodoPosXMinusCaloPosXVsCry_mean =0;
float HodoPosXMinusCaloPosXVsCry_rms =0;
float HodoPosYMinusCaloPosYVsCry_mean =0;
float HodoPosYMinusCaloPosYVsCry_rms =0;
float TDCTimeMinusCaloTimeVsCry_mean =0;
float TDCTimeMinusCaloTimeVsCry_rms =0;
// min number of entries chosen assuming:
// prescaling = 100 X 2FU
// that we want at leas 2k events per crystal
if (meCaloVsHodoXPos_-> getEntries() > 10){
HodoPosXMinusCaloPosXVsCry_mean = meCaloVsHodoXPos_-> getMean(1);
HodoPosXMinusCaloPosXVsCry_rms = meCaloVsHodoXPos_-> getRMS(1);
meHodoPosXMinusCaloPosXVsCry_-> setBinContent( cryInBeamCounter_, HodoPosXMinusCaloPosXVsCry_mean);
meHodoPosXMinusCaloPosXVsCry_-> setBinError( cryInBeamCounter_, HodoPosXMinusCaloPosXVsCry_rms);
LogDebug("EcalBeamTask") << "At event number: " << LV1_ << " step: " << cryInBeamCounter_
<< " DeltaPosX is: " << (meCaloVsHodoXPos_-> getMean(1))
<< " +-" << ( meCaloVsHodoXPos_-> getRMS(1));
}
if (meCaloVsHodoYPos_-> getEntries() > 10){
HodoPosYMinusCaloPosYVsCry_mean = meCaloVsHodoYPos_-> getMean(1);
HodoPosYMinusCaloPosYVsCry_rms = meCaloVsHodoYPos_-> getRMS(1);
meHodoPosYMinusCaloPosYVsCry_-> setBinContent( cryInBeamCounter_, HodoPosYMinusCaloPosYVsCry_mean);
meHodoPosYMinusCaloPosYVsCry_-> setBinError( cryInBeamCounter_, HodoPosYMinusCaloPosYVsCry_rms);
LogDebug("EcalBeamTask") << "At event number: " << LV1_ << " step: " << cryInBeamCounter_
<< " DeltaPosY is: " << (meCaloVsHodoYPos_-> getMean(1))
<< " +-" << ( meCaloVsHodoYPos_-> getRMS(1));
}
if (meCaloVsTDCTime_-> getEntries() > 10){
TDCTimeMinusCaloTimeVsCry_mean = meCaloVsTDCTime_-> getMean(1);
TDCTimeMinusCaloTimeVsCry_rms = meCaloVsTDCTime_-> getRMS(1);
meTDCTimeMinusCaloTimeVsCry_-> setBinContent(cryInBeamCounter_, TDCTimeMinusCaloTimeVsCry_mean);
meTDCTimeMinusCaloTimeVsCry_-> setBinError(cryInBeamCounter_, TDCTimeMinusCaloTimeVsCry_rms);
LogDebug("EcalBeamTask") << "At event number: " << LV1_ << " step: " << cryInBeamCounter_
<< " DeltaT is: " << (meCaloVsTDCTime_-> getMean(1))
<< " +-" << ( meCaloVsTDCTime_-> getRMS(1));
}
LogDebug("EcalBeamTask") << "At event number: " << LV1_ << " trace histos filled ( cryInBeamCounter_="
<< cryInBeamCounter_ << ")";
}
// if table has come to rest (from movement), reset concerned ME's
if (resetNow_)
{
meEvsXRecProf_->Reset();
meEvsYRecProf_->Reset();
meEvsXRecHis_->Reset();
meEvsYRecHis_->Reset();
meCaloVsHodoXPos_->Reset();
meCaloVsHodoYPos_->Reset();
meCaloVsTDCTime_->Reset();
resetNow_ = false;
}
/**************************************
// handling histos type I:
**************************************/
for (unsigned int planeId=0; planeId <4; planeId++){
const EcalTBHodoscopePlaneRawHits& planeRaw = rawHodo->getPlaneRawHits(planeId);
LogDebug("EcalBeamTask") << "\t plane: " << (planeId+1)
<< "\t number of fibers: " << planeRaw.channels()
<< "\t number of hits: " << planeRaw.numberOfFiredHits();
meHodoOcc_[planeId]-> Fill( planeRaw.numberOfFiredHits() );
for (unsigned int i=0;i<planeRaw.channels();i++)
{
if (planeRaw.isChannelFired(i))
{
edm::LogInfo("EcalBeamTask")<< " channel " << (i+1) << " has fired";
meHodoRaw_[planeId]-> Fill(i+0.5);
}
}
}
edm::Handle<EcalTBTDCRecInfo> pTDC;
const EcalTBTDCRecInfo* recTDC=0;
if ( e.getByLabel(EcalTBTDCRecInfo_, pTDC) ) {
recTDC = pTDC.product();
LogDebug("EBBeamHodoTask") << " TDC offset is: " << recTDC->offset();
} else {
edm::LogError("EcalBeamTask") << "Error! Can't get the product EcalTBTDCRecInfo. Returning";
meMissingCollections_-> Fill(5); // bin6: missing reconstructed TDC
return;
}
edm::Handle<EcalTBHodoscopeRecInfo> pHodo;
const EcalTBHodoscopeRecInfo* recHodo=0;
if ( e.getByLabel(EcalTBHodoscopeRecInfo_, pHodo) ) {
recHodo = pHodo.product();
LogDebug("EcalBeamTask") << "hodoscopeReco: x: " << recHodo->posX()
<< "\ty: " << recHodo->posY()
<< "\t sx: " << recHodo->slopeX() << "\t qualx: " << recHodo->qualX()
<< "\t sy: " << recHodo->slopeY() << "\t qualy: " << recHodo->qualY();
} else {
edm::LogError("EcalBeamTask") << "Error! Can't get the product EcalTBHodoscopeRecInfo";
meMissingCollections_-> Fill(6); // bin7: missing reconstructed hodoscopes
return;
}
if ( (!recHodo) || (!recTDC) ) {
edm::LogWarning("EcalBeamTask") << "analyze: missing a needed collection, recHodo or recTDC. Returning.";
return;
}
LogDebug("EBBeamHodoTask") << " Hodo reco and TDC reco found.";
meTDCRec_->Fill( recTDC->offset());
meHodoPosRecXY_->Fill( recHodo->posX(), recHodo->posY() );
meHodoPosRecX_->Fill( recHodo->posX());
meHodoPosRecY_->Fill( recHodo->posY() );
meHodoSloXRec_->Fill( recHodo->slopeX());
meHodoSloYRec_->Fill( recHodo->slopeY());
meHodoQuaXRec_->Fill( recHodo->qualX());
meHodoQuaYRec_->Fill( recHodo->qualY());
/**************************************
// handling histos type II:
**************************************/
if (tableIsMoving_) {
LogDebug("EcalBeamTask")<< "At event number:" << LV1_ << " table is moving. Not filling concerned monitoring elements. ";
return;
} else {
LogDebug("EcalBeamTask")<< "At event number:" << LV1_ << " table is not moving - thus filling alos monitoring elements requiring so.";
}
float maxE =0;
EBDetId maxHitId(0);
for ( EBUncalibratedRecHitCollection::const_iterator uncalHitItr = pUncalRH->begin(); uncalHitItr!= pUncalRH->end(); uncalHitItr++ ) {
double e = (*uncalHitItr).amplitude();
if ( e <= 0. ) e = 0.0;
if ( e > maxE ) {
maxE = e;
maxHitId = (*uncalHitItr).id();
}
}
if ( maxHitId == EBDetId(0) ) {
edm::LogError("EBBeamHodoTask") << "No positive UncalRecHit found in ECAL in event " << ievt_ << " - returning.";
return;
}
meEvsXRecProf_-> Fill(recHodo->posX(), maxE);
meEvsYRecProf_-> Fill(recHodo->posY(), maxE);
meEvsXRecHis_-> Fill(recHodo->posX(), maxE);
meEvsYRecHis_-> Fill(recHodo->posY(), maxE);
edm::LogInfo("EcalBeamTask")<< " channel with max is " << maxHitId;
bool mat5x5 =true;
int ietaMax = maxHitId.ieta();
int iphiMax = (maxHitId.iphi() % 20);
if (ietaMax ==1 || ietaMax ==2 || ietaMax ==84 || ietaMax == 85 ||
iphiMax ==1 || iphiMax ==2 || iphiMax ==19 || iphiMax == 20 )
{mat5x5 =false;}
if (!mat5x5) return;
EBDetId Xtals5x5[25];
double ene5x5[25];
double e25 =0;
for (unsigned int icry=0;icry<25;icry++)
{
unsigned int row = icry / 5;
unsigned int column= icry %5;
if ( EBDetId::validDetId(maxHitId.ieta()+column-2,maxHitId.iphi()+row-2) ) {
Xtals5x5[icry]=EBDetId(maxHitId.ieta()+column-2,maxHitId.iphi()+row-2,EBDetId::ETAPHIMODE);
double e = ( *pUncalRH->find( Xtals5x5[icry] ) ).amplitude();
if ( e <= 0. ) e = 0.0;
ene5x5[icry] =e;
e25 +=e;
} else {
LogDebug("EcalBeamTask")<< "Cannot construct 5x5 matrix around EBDetId " << maxHitId;
mat5x5 =false;
}
}
if (!mat5x5) return;
LogDebug("EcalBeamTask")<< "Could construct 5x5 matrix around EBDetId " << maxHitId;
// im mm
float sideX =24.06;
float sideY =22.02;
float caloX =0;
float caloY =0;
float weight=0;
float sumWeight=0;
// X and Y calculated from log-weighted energy center of mass
for (unsigned int icry=0;icry<25;icry++)
{
weight = log( ene5x5[icry] / e25) + 3.8;
if (weight>0)
{
unsigned int row = icry / 5;
unsigned int column = icry %5;
caloX += (column-2) * sideX * weight;
caloY -= (row-2) * sideY * weight;
sumWeight += weight;
}
}
caloX /=sumWeight;
caloY /=sumWeight;
meCaloVsHodoXPos_->Fill( recHodo->posX()-caloX );
meCaloVsHodoYPos_->Fill( recHodo->posY()-caloY );
meCaloVsTDCTime_->Fill( (*pUncalRH->find( maxHitId ) ).jitter() - recTDC->offset() - 3);
LogDebug("EcalBeamTask")<< "jiitter from uncalRecHit: " << (*pUncalRH->find( maxHitId ) ).jitter();
}
| void EBBeamHodoTask::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 93 of file EBBeamHodoTask.cc.
References cryInBeamCounter_, dqmStore_, ievt_, LV1_, prefixME_, resetNow_, DQMStore::rmdir(), and DQMStore::setCurrentFolder().
| void EBBeamHodoTask::beginRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
BeginRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 108 of file EBBeamHodoTask.cc.
References Numbers::initGeometry(), mergeRuns_, and reset().
{
Numbers::initGeometry(c, false);
if ( ! mergeRuns_ ) this->reset();
}
| void EBBeamHodoTask::cleanup | ( | void | ) | [protected] |
Cleanup.
Definition at line 267 of file EBBeamHodoTask.cc.
References dqmStore_, edm::getName(), MonitorElement::getName(), i, init_, meCaloVsHodoXPos_, meCaloVsHodoYPos_, meCaloVsTDCTime_, meEvsXRecHis_, meEvsXRecProf_, meEvsYRecHis_, meEvsYRecProf_, meHodoOcc_, meHodoPosRecX_, meHodoPosRecXY_, meHodoPosRecY_, meHodoPosXMinusCaloPosXVsCry_, meHodoPosYMinusCaloPosYVsCry_, meHodoQuaXRec_, meHodoQuaYRec_, meHodoRaw_, meHodoSloXRec_, meHodoSloYRec_, meMissingCollections_, meTDCRec_, meTDCTimeMinusCaloTimeVsCry_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().
Referenced by endJob().
{
if ( ! init_ ) return;
if ( dqmStore_ ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EBBeamHodoTask");
for (int i=0; i<4; i++) {
if ( meHodoOcc_[i] ) dqmStore_->removeElement( meHodoOcc_[i]->getName() );
meHodoOcc_[i] = 0;
if ( meHodoRaw_[i] ) dqmStore_->removeElement( meHodoRaw_[i]->getName() );
meHodoRaw_[i] = 0;
}
if ( meHodoPosRecX_ ) dqmStore_->removeElement( meHodoPosRecX_->getName() );
meHodoPosRecX_ = 0;
if ( meHodoPosRecY_ ) dqmStore_->removeElement( meHodoPosRecY_->getName() );
meHodoPosRecY_ = 0;
if ( meHodoPosRecXY_ ) dqmStore_->removeElement( meHodoPosRecXY_->getName() );
meHodoPosRecXY_ = 0;
if ( meHodoSloXRec_ ) dqmStore_->removeElement( meHodoSloXRec_->getName() );
meHodoSloXRec_ = 0;
if ( meHodoSloYRec_ ) dqmStore_->removeElement( meHodoSloYRec_->getName() );
meHodoSloYRec_ = 0;
if ( meHodoQuaXRec_ ) dqmStore_->removeElement( meHodoQuaXRec_->getName() );
meHodoQuaXRec_ = 0;
if ( meHodoQuaYRec_ ) dqmStore_->removeElement( meHodoQuaYRec_->getName() );
meHodoQuaYRec_ = 0;
if ( meTDCRec_ ) dqmStore_->removeElement( meTDCRec_->getName() );
meTDCRec_ = 0;
if ( meEvsXRecProf_ ) dqmStore_->removeElement( meEvsXRecProf_->getName() );
meEvsXRecProf_ = 0;
if ( meEvsYRecProf_ ) dqmStore_->removeElement( meEvsYRecProf_->getName() );
meEvsYRecProf_ = 0;
if ( meEvsXRecHis_ ) dqmStore_->removeElement( meEvsXRecHis_->getName() );
meEvsXRecHis_ = 0;
if ( meEvsYRecHis_ ) dqmStore_->removeElement( meEvsYRecHis_->getName() );
meEvsYRecHis_ = 0;
if ( meCaloVsHodoXPos_ ) dqmStore_->removeElement( meCaloVsHodoXPos_->getName() );
meCaloVsHodoXPos_ = 0;
if ( meCaloVsHodoYPos_ ) dqmStore_->removeElement( meCaloVsHodoYPos_->getName() );
meCaloVsHodoYPos_ = 0;
if ( meCaloVsTDCTime_ ) dqmStore_->removeElement( meCaloVsTDCTime_->getName() );
meCaloVsTDCTime_ = 0;
if ( meHodoPosXMinusCaloPosXVsCry_ ) dqmStore_->removeElement( meHodoPosXMinusCaloPosXVsCry_->getName() );
meHodoPosXMinusCaloPosXVsCry_ = 0;
if ( meHodoPosYMinusCaloPosYVsCry_ ) dqmStore_->removeElement( meHodoPosYMinusCaloPosYVsCry_->getName() );
meHodoPosYMinusCaloPosYVsCry_ = 0;
if ( meTDCTimeMinusCaloTimeVsCry_ ) dqmStore_->removeElement( meTDCTimeMinusCaloTimeVsCry_->getName() );
meTDCTimeMinusCaloTimeVsCry_ = 0;
if ( meMissingCollections_ ) dqmStore_->removeElement( meMissingCollections_->getName() );
meMissingCollections_ = 0;
}
init_ = false;
}
| void EBBeamHodoTask::endJob | ( | void | ) | [protected, virtual] |
EndJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 326 of file EBBeamHodoTask.cc.
References cleanup(), enableCleanup_, and ievt_.
{
edm::LogInfo("EBBeamHodoTask") << "analyzed " << ievt_ << " events";
if ( enableCleanup_ ) this->cleanup();
}
| void EBBeamHodoTask::endRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
| void EBBeamHodoTask::reset | ( | void | ) | [protected] |
Reset.
Definition at line 120 of file EBBeamHodoTask.cc.
References i, meCaloVsHodoXPos_, meCaloVsHodoYPos_, meCaloVsTDCTime_, meEvsXRecHis_, meEvsXRecProf_, meEvsYRecHis_, meEvsYRecProf_, meHodoOcc_, meHodoPosRecX_, meHodoPosRecXY_, meHodoPosRecY_, meHodoPosXMinusCaloPosXVsCry_, meHodoPosYMinusCaloPosYVsCry_, meHodoQuaXRec_, meHodoQuaYRec_, meHodoRaw_, meHodoSloXRec_, meHodoSloYRec_, meMissingCollections_, meTDCRec_, meTDCTimeMinusCaloTimeVsCry_, and MonitorElement::Reset().
Referenced by beginRun().
{
for (int i=0; i<4; i++) {
if ( meHodoOcc_[i] ) meHodoOcc_[i]->Reset();
if ( meHodoRaw_[i] ) meHodoRaw_[i]->Reset();
}
if ( meHodoPosRecX_ ) meHodoPosRecX_->Reset();
if ( meHodoPosRecY_ ) meHodoPosRecY_->Reset();
if ( meHodoPosRecXY_ ) meHodoPosRecXY_->Reset();
if ( meHodoSloXRec_ ) meHodoSloXRec_->Reset();
if ( meHodoSloYRec_ ) meHodoSloYRec_->Reset();
if ( meHodoQuaXRec_ ) meHodoQuaXRec_->Reset();
if ( meHodoQuaYRec_ ) meHodoQuaYRec_->Reset();
if ( meTDCRec_ ) meTDCRec_->Reset();
if ( meEvsXRecProf_ ) meEvsXRecProf_->Reset();
if ( meEvsYRecProf_ ) meEvsYRecProf_->Reset();
if ( meEvsXRecHis_ ) meEvsXRecHis_->Reset();
if ( meEvsYRecHis_ ) meEvsYRecHis_->Reset();
if ( meCaloVsHodoXPos_ ) meCaloVsHodoXPos_->Reset();
if ( meCaloVsHodoYPos_ ) meCaloVsHodoYPos_->Reset();
if ( meCaloVsTDCTime_ ) meCaloVsTDCTime_->Reset();
if ( meHodoPosXMinusCaloPosXVsCry_ ) meHodoPosXMinusCaloPosXVsCry_->Reset();
if ( meHodoPosYMinusCaloPosYVsCry_ ) meHodoPosYMinusCaloPosYVsCry_->Reset();
if ( meTDCTimeMinusCaloTimeVsCry_ ) meTDCTimeMinusCaloTimeVsCry_->Reset();
if ( meMissingCollections_ ) meMissingCollections_->Reset();
}
| void EBBeamHodoTask::setup | ( | void | ) | [protected] |
Setup.
Definition at line 149 of file EBBeamHodoTask.cc.
References DQMStore::book1D(), DQMStore::book2D(), dqmStore_, i, init_, meCaloVsHodoXPos_, meCaloVsHodoYPos_, meCaloVsTDCTime_, meEvsXRecHis_, meEvsXRecProf_, meEvsYRecHis_, meEvsYRecProf_, meHodoOcc_, meHodoPosRecX_, meHodoPosRecXY_, meHodoPosRecY_, meHodoPosXMinusCaloPosXVsCry_, meHodoPosYMinusCaloPosYVsCry_, meHodoQuaXRec_, meHodoQuaYRec_, meHodoRaw_, meHodoSloXRec_, meHodoSloYRec_, meMissingCollections_, meTDCRec_, meTDCTimeMinusCaloTimeVsCry_, mergeVDriftHistosByStation::name, prefixME_, Numbers::sEB(), MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), and smId.
Referenced by analyze().
{
init_ = true;
smId =1;
if ( dqmStore_ ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EBBeamHodoTask");
std::stringstream ss;
std::string name;
// following ME (type I):
// *** do not need to be ever reset
// *** can be filled regardless of the moving/notMoving status of the table
for (int i=0; i<4; i++) {
ss << std::setw(2) << std::setfill('0') << i+1;
name = "EBBHT occup " + Numbers::sEB(smId) + " " + ss.str();
meHodoOcc_[i] = dqmStore_->book1D(name, name, 30, 0., 30.);
meHodoOcc_[i]->setAxisTitle("hits per event", 1);
name = "EBBHT raw " + Numbers::sEB(smId) + " " + ss.str();
meHodoRaw_[i] = dqmStore_->book1D(name, name, 64, 0., 64.);
meHodoRaw_[i]->setAxisTitle("hodo fiber number", 1);
}
name = "EBBHT PosX rec " + Numbers::sEB(smId);
meHodoPosRecX_ = dqmStore_->book1D(name, name, 100, -20, 20);
meHodoPosRecX_->setAxisTitle("reconstructed position (mm)", 1);
name = "EBBHT PosY rec " + Numbers::sEB(smId);
meHodoPosRecY_ = dqmStore_->book1D(name, name, 100, -20, 20);
meHodoPosRecY_->setAxisTitle("reconstructed position (mm)", 1);
name = "EBBHT PosYX rec " + Numbers::sEB(smId);
meHodoPosRecXY_ = dqmStore_->book2D(name, name, 100, -20, 20,100, -20, 20);
meHodoPosRecXY_->setAxisTitle("reconstructed X position (mm)", 1);
meHodoPosRecXY_->setAxisTitle("reconstructed Y position (mm)", 2);
name = "EBBHT SloX " + Numbers::sEB(smId);
meHodoSloXRec_ = dqmStore_->book1D(name, name, 50, -0.005, 0.005);
meHodoSloXRec_->setAxisTitle("reconstructed track slope", 1);
name = "EBBHT SloY " + Numbers::sEB(smId);
meHodoSloYRec_ = dqmStore_->book1D(name, name, 50, -0.005, 0.005);
meHodoSloYRec_->setAxisTitle("reconstructed track slope", 1);
name = "EBBHT QualX " + Numbers::sEB(smId);
meHodoQuaXRec_ = dqmStore_->book1D(name, name, 50, 0, 5);
meHodoQuaXRec_->setAxisTitle("track fit quality", 1);
name = "EBBHT QualY " + Numbers::sEB(smId);
meHodoQuaYRec_ = dqmStore_->book1D(name, name, 50, 0, 5);
meHodoQuaYRec_->setAxisTitle("track fit quality", 1);
name = "EBBHT TDC rec " + Numbers::sEB(smId);
meTDCRec_ = dqmStore_->book1D(name, name, 25, 0, 1);
meTDCRec_->setAxisTitle("offset", 1);
name = "EBBHT Hodo-Calo X vs Cry " + Numbers::sEB(smId);
meHodoPosXMinusCaloPosXVsCry_ = dqmStore_->book1D(name, name, 50, 0, 50);
meHodoPosXMinusCaloPosXVsCry_->setAxisTitle("scan step number", 1);
meHodoPosXMinusCaloPosXVsCry_->setAxisTitle("PosX_{hodo} - PosX_{calo} (mm)", 2);
name = "EBBHT Hodo-Calo Y vs Cry " + Numbers::sEB(smId);
meHodoPosYMinusCaloPosYVsCry_ = dqmStore_->book1D(name, name, 50, 0, 50);
meHodoPosYMinusCaloPosYVsCry_->setAxisTitle("scan step number", 1);
meHodoPosYMinusCaloPosYVsCry_->setAxisTitle("PosY_{hodo} - PosY_{calo} (mm)", 2);
name = "EBBHT TDC-Calo vs Cry " + Numbers::sEB(smId);
meTDCTimeMinusCaloTimeVsCry_ = dqmStore_->book1D(name, name, 50, 0, 50);
meTDCTimeMinusCaloTimeVsCry_->setAxisTitle("scan step number", 1);
meTDCTimeMinusCaloTimeVsCry_->setAxisTitle("Time_{TDC} - Time_{calo} (sample)", 2);
name = "EBBHT Missing Collections " + Numbers::sEB(smId);
meMissingCollections_ = dqmStore_->book1D(name, name, 7, 0, 7);
meMissingCollections_->setAxisTitle("missing collection", 1);
// following ME (type II):
// *** can be filled only when table is **not** moving
// *** need to be reset once table goes from 'moving'->notMoving
name = "EBBHT prof E1 vs X " + Numbers::sEB(smId);
meEvsXRecProf_ = dqmStore_-> bookProfile(name, name, 100, -20, 20, 500, 0, 5000, "s");
meEvsXRecProf_->setAxisTitle("PosX (mm)", 1);
meEvsXRecProf_->setAxisTitle("E1 (ADC)", 2);
name = "EBBHT prof E1 vs Y " + Numbers::sEB(smId);
meEvsYRecProf_ = dqmStore_-> bookProfile(name, name, 100, -20, 20, 500, 0, 5000, "s");
meEvsYRecProf_->setAxisTitle("PosY (mm)", 1);
meEvsYRecProf_->setAxisTitle("E1 (ADC)", 2);
name = "EBBHT his E1 vs X " + Numbers::sEB(smId);
meEvsXRecHis_ = dqmStore_-> book2D(name, name, 100, -20, 20, 500, 0, 5000);
meEvsXRecHis_->setAxisTitle("PosX (mm)", 1);
meEvsXRecHis_->setAxisTitle("E1 (ADC)", 2);
name = "EBBHT his E1 vs Y " + Numbers::sEB(smId);
meEvsYRecHis_ = dqmStore_-> book2D(name, name, 100, -20, 20, 500, 0, 5000);
meEvsYRecHis_->setAxisTitle("PosY (mm)", 1);
meEvsYRecHis_->setAxisTitle("E1 (ADC)", 2);
name = "EBBHT PosX Hodo-Calo " + Numbers::sEB(smId);
meCaloVsHodoXPos_ = dqmStore_->book1D(name, name, 40, -20, 20);
meCaloVsHodoXPos_->setAxisTitle("PosX_{hodo} - PosX_{calo} (mm)", 1);
name = "EBBHT PosY Hodo-Calo " + Numbers::sEB(smId);
meCaloVsHodoYPos_ = dqmStore_->book1D(name, name, 40, -20, 20);
meCaloVsHodoYPos_->setAxisTitle("PosY_{hodo} - PosY_{calo} (mm)", 1);
name = "EBBHT TimeMax TDC-Calo " + Numbers::sEB(smId);
meCaloVsTDCTime_ = dqmStore_->book1D(name, name, 100, -1, 1);//tentative
meCaloVsTDCTime_->setAxisTitle("Time_{TDC} - Time_{calo} (samples)", 1);
}
}
int EBBeamHodoTask::cryInBeam_ [private] |
Definition at line 80 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
int EBBeamHodoTask::cryInBeamCounter_ [private] |
Definition at line 82 of file EBBeamHodoTask.h.
Referenced by analyze(), and beginJob().
DQMStore* EBBeamHodoTask::dqmStore_ [private] |
Definition at line 61 of file EBBeamHodoTask.h.
Referenced by beginJob(), cleanup(), EBBeamHodoTask(), and setup().
Definition at line 70 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
Definition at line 69 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
Definition at line 73 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
Definition at line 75 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
Definition at line 72 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
Definition at line 74 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
Definition at line 71 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
bool EBBeamHodoTask::enableCleanup_ [private] |
Definition at line 65 of file EBBeamHodoTask.h.
Referenced by EBBeamHodoTask(), and endJob().
int EBBeamHodoTask::ievt_ [private] |
Definition at line 59 of file EBBeamHodoTask.h.
Referenced by analyze(), beginJob(), and endJob().
bool EBBeamHodoTask::init_ [private] |
Definition at line 109 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), and setup().
int EBBeamHodoTask::LV1_ [private] |
Definition at line 77 of file EBBeamHodoTask.h.
Referenced by analyze(), and beginJob().
MonitorElement* EBBeamHodoTask::meCaloVsHodoXPos_ [private] |
Definition at line 105 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meCaloVsHodoYPos_ [private] |
Definition at line 106 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meCaloVsTDCTime_ [private] |
Definition at line 107 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meEvsXRecHis_ [private] |
Definition at line 103 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meEvsXRecProf_ [private] |
Definition at line 101 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meEvsYRecHis_ [private] |
Definition at line 104 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meEvsYRecProf_ [private] |
Definition at line 102 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoOcc_[4] [private] |
Definition at line 85 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoPosRecX_ [private] |
Definition at line 88 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoPosRecXY_ [private] |
Definition at line 87 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoPosRecY_ [private] |
Definition at line 89 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
Definition at line 95 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
Definition at line 96 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoQuaXRec_ [private] |
Definition at line 92 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoQuaYRec_ [private] |
Definition at line 93 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoRaw_[4] [private] |
Definition at line 86 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoSloXRec_ [private] |
Definition at line 90 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
MonitorElement* EBBeamHodoTask::meHodoSloYRec_ [private] |
Definition at line 91 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
Definition at line 98 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
bool EBBeamHodoTask::mergeRuns_ [private] |
Definition at line 67 of file EBBeamHodoTask.h.
Referenced by beginRun(), and EBBeamHodoTask().
MonitorElement* EBBeamHodoTask::meTDCRec_ [private] |
Definition at line 94 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
Definition at line 97 of file EBBeamHodoTask.h.
Referenced by analyze(), cleanup(), EBBeamHodoTask(), reset(), and setup().
std::string EBBeamHodoTask::prefixME_ [private] |
Definition at line 63 of file EBBeamHodoTask.h.
Referenced by beginJob(), cleanup(), EBBeamHodoTask(), and setup().
int EBBeamHodoTask::previousCryInBeam_ [private] |
Definition at line 81 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
bool EBBeamHodoTask::resetNow_ [private] |
Definition at line 79 of file EBBeamHodoTask.h.
Referenced by analyze(), and beginJob().
int EBBeamHodoTask::smId [private] |
Definition at line 111 of file EBBeamHodoTask.h.
Referenced by setup().
bool EBBeamHodoTask::tableIsMoving_ [private] |
Definition at line 78 of file EBBeamHodoTask.h.
Referenced by analyze(), and EBBeamHodoTask().
1.7.3