|
|
|
#include <EcalLaserDbService.h>
Public Member Functions | |
| EcalLaserDbService () | |
| EcalLaserDbService (const edm::ParameterSet &) | |
| const EcalLaserAlphas * | getAlphas () const |
| const EcalLaserAPDPNRatios * | getAPDPNRatios () const |
| const EcalLaserAPDPNRatiosRef * | getAPDPNRatiosRef () const |
| float | getLaserCorrection (DetId const &xid, edm::Timestamp const &iTime) const |
| const EcalLinearCorrections * | getLinearCorrections () const |
| void | setAlphaData (const EcalLaserAlphas *fItem) |
| void | setAPDPNData (const EcalLaserAPDPNRatios *fItem) |
| void | setAPDPNRefData (const EcalLaserAPDPNRatiosRef *fItem) |
| void | setLinearCorrectionsData (const EcalLinearCorrections *fItem) |
Private Attributes | |
| const EcalLaserAlphas * | mAlphas_ |
| const EcalLaserAPDPNRatios * | mAPDPNRatios_ |
| const EcalLaserAPDPNRatiosRef * | mAPDPNRatiosRef_ |
| const EcalLinearCorrections * | mLinearCorrections_ |
Definition at line 25 of file EcalLaserDbService.h.
| EcalLaserDbService::EcalLaserDbService | ( | ) |
Definition at line 16 of file EcalLaserDbService.cc.
: mAlphas_ (0), mAPDPNRatiosRef_ (0), mAPDPNRatios_ (0), mLinearCorrections_ (0) {}
| EcalLaserDbService::EcalLaserDbService | ( | const edm::ParameterSet & | ) |
| const EcalLaserAlphas * EcalLaserDbService::getAlphas | ( | ) | const |
| const EcalLaserAPDPNRatios * EcalLaserDbService::getAPDPNRatios | ( | ) | const |
Definition at line 34 of file EcalLaserDbService.cc.
References mAPDPNRatios_.
{
return mAPDPNRatios_;
}
| const EcalLaserAPDPNRatiosRef * EcalLaserDbService::getAPDPNRatiosRef | ( | ) | const |
Definition at line 30 of file EcalLaserDbService.cc.
References mAPDPNRatiosRef_.
{
return mAPDPNRatiosRef_;
}
| float EcalLaserDbService::getLaserCorrection | ( | DetId const & | xid, |
| edm::Timestamp const & | iTime | ||
| ) | const |
Definition at line 43 of file EcalLaserDbService.cc.
References alpha, DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, EcalCondObjectContainer< T >::end(), EcalCondObjectContainer< T >::find(), EcalLaserAPDPNRatios::getLaserMap(), EcalCondObjectContainer< T >::getMap(), EcalTimeDependentCorrections::getTimeMap(), EcalLaserAPDPNRatios::getTimeMap(), EcalTimeDependentCorrections::getValueMap(), MEEEGeom::lmr(), MEEBGeom::lmr(), mAlphas_, mAPDPNRatios_, mAPDPNRatiosRef_, mLinearCorrections_, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p1, EcalTimeDependentCorrections::Values::p1, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p2, EcalTimeDependentCorrections::Values::p2, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p3, EcalTimeDependentCorrections::Values::p3, funct::pow(), DetId::rawId(), DetId::subdetId(), lumiQTWidget::t, EcalTimeDependentCorrections::Times::t1, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t1, EcalTimeDependentCorrections::Times::t2, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t2, EcalTimeDependentCorrections::Times::t3, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t3, cond::timestamp, and edm::Timestamp::value().
Referenced by EcalHitResponse::findLaserConstant().
{
float correctionFactor = 1.0;
const EcalLaserAPDPNRatios::EcalLaserAPDPNRatiosMap& laserRatiosMap = mAPDPNRatios_->getLaserMap();
const EcalLaserAPDPNRatios::EcalLaserTimeStampMap& laserTimeMap = mAPDPNRatios_->getTimeMap();
const EcalLaserAPDPNRatiosRefMap& laserRefMap = mAPDPNRatiosRef_->getMap();
const EcalLaserAlphaMap& laserAlphaMap = mAlphas_->getMap();
const EcalLinearCorrections::EcalValueMap& linearValueMap = mLinearCorrections_->getValueMap();
const EcalLinearCorrections::EcalTimeMap& linearTimeMap = mLinearCorrections_->getTimeMap();
EcalLaserAPDPNRatios::EcalLaserAPDPNpair apdpnpair;
EcalLaserAPDPNRatios::EcalLaserTimeStamp timestamp;
EcalLaserAPDPNref apdpnref;
EcalLaserAlpha alpha;
EcalLinearCorrections::Values linValues;
EcalLinearCorrections::Times linTimes;
if (xid.det()==DetId::Ecal) {
// std::cout << " XID is in Ecal : ";
} else {
// std::cout << " XID is NOT in Ecal : ";
edm::LogError("EcalLaserDbService") << " DetId is NOT in ECAL" << endl;
return correctionFactor;
}
// int hi = -1;
// if (xid.subdetId()==EcalBarrel) {
// // std::cout << "EcalBarrel" << std::endl;
// // std::cout << "--> rawId() = " << xid.rawId() << " id() = " << EBDetId( xid ).hashedIndex() << std::endl;
// hi = EBDetId( xid ).hashedIndex();
// } else if (xid.subdetId()==EcalEndcap) {
// // std::cout << "EcalEndcap" << std::endl;
// hi = EEDetId( xid ).hashedIndex() + EBDetId::MAX_HASH + 1;
//
// } else {
// // std::cout << "NOT EcalBarrel or EcalEndCap" << std::endl;
// edm::LogError("EcalLaserDbService") << " DetId is NOT in ECAL Barrel or Endcap" << endl;
// return correctionFactor;
// }
int iLM;
if (xid.subdetId()==EcalBarrel) {
EBDetId ebid( xid.rawId() );
iLM = MEEBGeom::lmr(ebid.ieta(), ebid.iphi());
} else if (xid.subdetId()==EcalEndcap) {
EEDetId eeid( xid.rawId() );
// SuperCrystal coordinates
MEEEGeom::SuperCrysCoord iX = (eeid.ix()-1)/5 + 1;
MEEEGeom::SuperCrysCoord iY = (eeid.iy()-1)/5 + 1;
iLM = MEEEGeom::lmr(iX, iY, eeid.zside());
} else {
edm::LogError("EcalLaserDbService") << " DetId is NOT in ECAL Barrel or Endcap" << endl;
return correctionFactor;
}
// std::cout << " LM num ====> " << iLM << endl;
// get alpha, apd/pn ref, apd/pn pairs and timestamps for interpolation
EcalLaserAPDPNRatios::EcalLaserAPDPNRatiosMap::const_iterator itratio = laserRatiosMap.find(xid);
if (itratio != laserRatiosMap.end()) {
apdpnpair = (*itratio);
} else {
edm::LogError("EcalLaserDbService") << "error with laserRatiosMap!" << endl;
return correctionFactor;
}
if (iLM-1< (int)laserTimeMap.size()) {
timestamp = laserTimeMap[iLM-1];
} else {
edm::LogError("EcalLaserDbService") << "error with laserTimeMap!" << endl;
return correctionFactor;
}
EcalLinearCorrections::EcalValueMap::const_iterator itlin = linearValueMap.find(xid);
if (itlin != linearValueMap.end()) {
linValues = (*itlin);
} else {
edm::LogError("EcalLaserDbService") << "error with linearValueMap!" << endl;
return correctionFactor;
}
if (iLM-1< (int)linearTimeMap.size()) {
linTimes = linearTimeMap[iLM-1];
} else {
edm::LogError("EcalLaserDbService") << "error with laserTimeMap!" << endl;
return correctionFactor;
}
EcalLaserAPDPNRatiosRefMap::const_iterator itref = laserRefMap.find(xid);
if ( itref != laserRefMap.end() ) {
apdpnref = (*itref);
} else {
edm::LogError("EcalLaserDbService") << "error with laserRefMap!" << endl;
return correctionFactor;
}
EcalLaserAlphaMap::const_iterator italpha = laserAlphaMap.find(xid);
if ( italpha != laserAlphaMap.end() ) {
alpha = (*italpha);
} else {
edm::LogError("EcalLaserDbService") << "error with laserAlphaMap!" << endl;
return correctionFactor;
}
// std::cout << " APDPN pair " << apdpnpair.p1 << " , " << apdpnpair.p2 << std::endl;
// std::cout << " TIME pair " << timestamp.t1.value() << " , " << timestamp.t2.value() << " iLM " << iLM << std::endl;
// std::cout << " LM module " << iLM << std::endl;
// std::cout << " APDPN ref " << apdpnref << std::endl;
// std::cout << " ALPHA " << alpha << std::endl;
// should implement some default in case of error...
// should do some quality checks first
// ...
// we will need to treat time differently...
// is time in DB same format as in MC? probably not...
// interpolation
edm::TimeValue_t t = iTime.value();
edm::TimeValue_t t_i = 0, t_f = 0;
float p_i = 0, p_f = 0;
edm::TimeValue_t lt_i = 0, lt_f = 0;
float lp_i = 0, lp_f = 0;
if ( t >= timestamp.t1.value() && t < timestamp.t2.value() ) {
t_i = timestamp.t1.value();
t_f = timestamp.t2.value();
p_i = apdpnpair.p1;
p_f = apdpnpair.p2;
} else if ( t >= timestamp.t2.value() && t <= timestamp.t3.value() ) {
t_i = timestamp.t2.value();
t_f = timestamp.t3.value();
p_i = apdpnpair.p2;
p_f = apdpnpair.p3;
} else if ( t < timestamp.t1.value() ) {
t_i = timestamp.t1.value();
t_f = timestamp.t2.value();
p_i = apdpnpair.p1;
p_f = apdpnpair.p2;
//edm::LogWarning("EcalLaserDbService") << "The event timestamp t=" << t
// << " is lower than t1=" << t_i << ". Extrapolating...";
} else if ( t > timestamp.t3.value() ) {
t_i = timestamp.t2.value();
t_f = timestamp.t3.value();
p_i = apdpnpair.p2;
p_f = apdpnpair.p3;
//edm::LogWarning("EcalLaserDbService") << "The event timestamp t=" << t
// << " is greater than t3=" << t_f << ". Extrapolating...";
}
if ( t >= linTimes.t1.value() && t < linTimes.t2.value() ) {
lt_i = linTimes.t1.value();
lt_f = linTimes.t2.value();
lp_i = linValues.p1;
lp_f = linValues.p2;
} else if ( t >= linTimes.t2.value() && t <= linTimes.t3.value() ) {
lt_i = linTimes.t2.value();
lt_f = linTimes.t3.value();
lp_i = linValues.p2;
lp_f = linValues.p3;
} else if ( t < linTimes.t1.value() ) {
lt_i = linTimes.t1.value();
lt_f = linTimes.t2.value();
lp_i = linValues.p1;
lp_f = linValues.p2;
//edm::LogWarning("EcalLaserDbService") << "The event timestamp t=" << t
// << " is lower than t1=" << t_i << ". Extrapolating...";
} else if ( t > linTimes.t3.value() ) {
lt_i = linTimes.t2.value();
lt_f = linTimes.t3.value();
lp_i = linValues.p2;
lp_f = linValues.p3;
//edm::LogWarning("EcalLaserDbService") << "The event timestamp t=" << t
// << " is greater than t3=" << t_f << ". Extrapolating...";
}
if ( apdpnref != 0 && (t_i - t_f) != 0 && (lt_i - lt_f) != 0) {
float interpolatedLaserResponse = p_i/apdpnref + (t-t_i)*(p_f-p_i)/apdpnref/(t_f-t_i);
float interpolatedLinearResponse = lp_i/apdpnref + (t-lt_i)*(lp_f-lp_i)/apdpnref/(lt_f-lt_i); // FIXED BY FC
if(interpolatedLinearResponse >2 || interpolatedLinearResponse <0.1) interpolatedLinearResponse=1;
if ( interpolatedLaserResponse <= 0 ) {
edm::LogWarning("EcalLaserDbService") << "The interpolated laser correction is <= zero! ("
<< interpolatedLaserResponse << "). Using 1. as correction factor.";
return correctionFactor;
} else {
float interpolatedTransparencyResponse = interpolatedLaserResponse / interpolatedLinearResponse;
correctionFactor = 1/( pow(interpolatedTransparencyResponse,alpha) *interpolatedLinearResponse );
}
} else {
edm::LogError("EcalLaserDbService")
<< "apdpnref (" << apdpnref << ") "
<< "or t_i-t_f (" << (t_i - t_f) << " is zero!";
return correctionFactor;
}
return correctionFactor;
}
| const EcalLinearCorrections * EcalLaserDbService::getLinearCorrections | ( | ) | const |
Definition at line 38 of file EcalLaserDbService.cc.
References mLinearCorrections_.
{
return mLinearCorrections_;
}
| void EcalLaserDbService::setAlphaData | ( | const EcalLaserAlphas * | fItem | ) | [inline] |
| void EcalLaserDbService::setAPDPNData | ( | const EcalLaserAPDPNRatios * | fItem | ) | [inline] |
Definition at line 38 of file EcalLaserDbService.h.
References mAPDPNRatios_.
{mAPDPNRatios_ = fItem;}
| void EcalLaserDbService::setAPDPNRefData | ( | const EcalLaserAPDPNRatiosRef * | fItem | ) | [inline] |
Definition at line 37 of file EcalLaserDbService.h.
References mAPDPNRatiosRef_.
{mAPDPNRatiosRef_ = fItem;}
| void EcalLaserDbService::setLinearCorrectionsData | ( | const EcalLinearCorrections * | fItem | ) | [inline] |
Definition at line 39 of file EcalLaserDbService.h.
References mLinearCorrections_.
{mLinearCorrections_ = fItem;}
const EcalLaserAlphas* EcalLaserDbService::mAlphas_ [private] |
Definition at line 43 of file EcalLaserDbService.h.
Referenced by getAlphas(), getLaserCorrection(), and setAlphaData().
const EcalLaserAPDPNRatios* EcalLaserDbService::mAPDPNRatios_ [private] |
Definition at line 45 of file EcalLaserDbService.h.
Referenced by getAPDPNRatios(), getLaserCorrection(), and setAPDPNData().
const EcalLaserAPDPNRatiosRef* EcalLaserDbService::mAPDPNRatiosRef_ [private] |
Definition at line 44 of file EcalLaserDbService.h.
Referenced by getAPDPNRatiosRef(), getLaserCorrection(), and setAPDPNRefData().
const EcalLinearCorrections* EcalLaserDbService::mLinearCorrections_ [private] |
Definition at line 46 of file EcalLaserDbService.h.
Referenced by getLaserCorrection(), getLinearCorrections(), and setLinearCorrectionsData().
1.7.3