|
|
|
#include <EcalTestDevDB.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &evt, const edm::EventSetup &evtSetup) |
| EcalTestDevDB (const edm::ParameterSet &iConfig) | |
| EcalADCToGeVConstant * | generateEcalADCToGeVConstant () |
| EcalGainRatios * | generateEcalGainRatios () |
| EcalIntercalibConstants * | generateEcalIntercalibConstants () |
| EcalLaserAlphas * | generateEcalLaserAlphas () |
| EcalLaserAPDPNRatios * | generateEcalLaserAPDPNRatios (uint32_t i_run) |
| EcalLaserAPDPNRatiosRef * | generateEcalLaserAPDPNRatiosRef () |
| EcalPedestals * | generateEcalPedestals () |
| EcalTBWeights * | generateEcalTBWeights () |
| EcalWeightXtalGroups * | generateEcalWeightXtalGroups () |
| ~EcalTestDevDB () | |
Private Attributes | |
| std::map< std::string, unsigned long long > | m_cacheIDs |
| unsigned long | m_firstRun |
| unsigned int | m_interval |
| unsigned long | m_lastRun |
| std::map< std::string, std::string > | m_records |
| std::string | m_timetype |
Definition at line 38 of file EcalTestDevDB.h.
| EcalTestDevDB::EcalTestDevDB | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 25 of file EcalTestDevDB.cc.
References edm::ParameterSet::getParameter(), i, m_cacheIDs, m_firstRun, m_interval, m_lastRun, m_records, record, GlobalPosition_Frontier_DevDB_cff::tag, and copyFileAlignEB_cfg::toCopy.
: m_timetype(iConfig.getParameter<std::string>("timetype")), m_cacheIDs(), m_records() { std::string container; std::string tag; std::string record; m_firstRun=static_cast<unsigned long>(atoi( iConfig.getParameter<std::string>("firstRun").c_str())); m_lastRun=static_cast<unsigned long>(atoi( iConfig.getParameter<std::string>("lastRun").c_str())); m_interval=static_cast<unsigned long>(atoi( iConfig.getParameter<std::string>("interval").c_str())); typedef std::vector< edm::ParameterSet > Parameters; Parameters toCopy = iConfig.getParameter<Parameters>("toCopy"); for(Parameters::iterator i = toCopy.begin(); i != toCopy.end(); ++i) { container = i->getParameter<std::string>("container"); record = i->getParameter<std::string>("record"); m_cacheIDs.insert( std::make_pair(container, 0) ); m_records.insert( std::make_pair(container, record) ); } }
| EcalTestDevDB::~EcalTestDevDB | ( | ) |
Definition at line 51 of file EcalTestDevDB.cc.
{
}
| void EcalTestDevDB::analyze | ( | const edm::Event & | evt, |
| const edm::EventSetup & | evtSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 56 of file EcalTestDevDB.cc.
References gather_cfg::cout, edm::IOVSyncValue::endOfTime(), Exception, generateEcalADCToGeVConstant(), generateEcalGainRatios(), generateEcalIntercalibConstants(), generateEcalLaserAlphas(), generateEcalLaserAPDPNRatios(), generateEcalLaserAPDPNRatiosRef(), generateEcalPedestals(), generateEcalTBWeights(), generateEcalWeightXtalGroups(), i, edm::Service< T >::isAvailable(), m_firstRun, m_interval, m_lastRun, m_records, record, and ecalTPGAnalyzer_cfg::recordName.
{
edm::Service<cond::service::PoolDBOutputService> dbOutput;
if ( !dbOutput.isAvailable() ) {
throw cms::Exception("PoolDBOutputService is not available");
}
std::string container;
std::string record;
typedef std::map<std::string, std::string>::const_iterator recordIter;
for (recordIter i = m_records.begin(); i != m_records.end(); ++i) {
container = (*i).first;
record = (*i).second;
std::string recordName = m_records[container];
// Loop through each of the runs
unsigned long nrec=(m_lastRun-m_firstRun)/m_interval+1;
unsigned long nstart=0;
if (m_firstRun == 0 && m_lastRun == 0) {
// it should do at least once the loop
nstart=0;
nrec=1;
}
for(unsigned long i=nstart; i<nrec; i++) {
unsigned long irun=m_firstRun+i*m_interval;
// Arguments 0 0 mean infinite IOV
if (m_firstRun == 0 && m_lastRun == 0) {
std::cout << "Infinite IOV mode" << std::endl;
irun = edm::IOVSyncValue::endOfTime().eventID().run();
}
std::cout << "Starting Transaction for run " << irun << "..." << std::flush;
if (container == "EcalPedestals") {
EcalPedestals* condObject= generateEcalPedestals();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalPedestals>( condObject, dbOutput->beginOfTime(),dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalPedestals>( condObject, irun , recordName);
}
} else if (container == "EcalADCToGeVConstant") {
EcalADCToGeVConstant* condObject= generateEcalADCToGeVConstant();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalADCToGeVConstant>( condObject, dbOutput->beginOfTime(),dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalADCToGeVConstant>( condObject, irun , recordName);
}
} else if (container == "EcalIntercalibConstants") {
EcalIntercalibConstants* condObject= generateEcalIntercalibConstants();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalIntercalibConstants>( condObject, dbOutput->beginOfTime(),dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalIntercalibConstants>( condObject, irun , recordName);
}
} else if (container == "EcalGainRatios") {
EcalGainRatios* condObject= generateEcalGainRatios();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalGainRatios>( condObject,dbOutput->beginOfTime(), dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalGainRatios>( condObject, irun , recordName);
}
} else if (container == "EcalWeightXtalGroups") {
EcalWeightXtalGroups* condObject= generateEcalWeightXtalGroups();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalWeightXtalGroups>( condObject, dbOutput->beginOfTime(), dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalWeightXtalGroups>( condObject, irun , recordName);
}
} else if (container == "EcalTBWeights") {
EcalTBWeights* condObject= generateEcalTBWeights();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalTBWeights>( condObject, dbOutput->beginOfTime(), dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalTBWeights>( condObject, irun , recordName);
}
} else if (container == "EcalLaserAPDPNRatios") {
EcalLaserAPDPNRatios* condObject= generateEcalLaserAPDPNRatios(irun);
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalLaserAPDPNRatios>( condObject,dbOutput->beginOfTime(), dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalLaserAPDPNRatios>( condObject, irun , recordName);
}
} else if (container == "EcalLaserAPDPNRatiosRef") {
EcalLaserAPDPNRatiosRef* condObject= generateEcalLaserAPDPNRatiosRef();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalLaserAPDPNRatiosRef>( condObject, dbOutput->beginOfTime(), dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalLaserAPDPNRatiosRef>( condObject, irun , recordName);
}
} else if (container == "EcalLaserAlphas") {
EcalLaserAlphas* condObject= generateEcalLaserAlphas();
if(irun==m_firstRun && dbOutput->isNewTagRequest(recordName)) {
// create new
std::cout<<"First One "<<std::endl;
dbOutput->createNewIOV<const EcalLaserAlphas>( condObject,dbOutput->beginOfTime(), dbOutput->endOfTime() ,recordName);
} else {
// append
std::cout<<"Old One "<<std::endl;
dbOutput->appendSinceTime<const EcalLaserAlphas>( condObject, irun , recordName);
}
} else {
std::cout << "it does not work yet for " << container << "..." << std::flush;
}
}
}
}
| EcalADCToGeVConstant * EcalTestDevDB::generateEcalADCToGeVConstant | ( | ) |
Definition at line 247 of file EcalTestDevDB.cc.
References csvReporter::r.
Referenced by analyze().
{
//-------------------------------------------------------------
double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
EcalADCToGeVConstant* agc = new EcalADCToGeVConstant(36.+r*4., 60.+r*4);
return agc;
}
| EcalGainRatios * EcalTestDevDB::generateEcalGainRatios | ( | ) |
Definition at line 277 of file EcalTestDevDB.cc.
References EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, csvReporter::r, DetId::rawId(), EcalMGPAGainRatio::setGain12Over6(), EcalMGPAGainRatio::setGain6Over1(), and EcalCondObjectContainer< T >::setValue().
Referenced by analyze().
{
//-------------------------------------------------------------
// create gain ratios
EcalGainRatios* gratio = new EcalGainRatios;
for(int ieta=-EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA; ++ieta) {
if(ieta==0) continue;
for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId:: MAX_IPHI; ++iphi) {
EBDetId ebid(ieta,iphi);
double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
EcalMGPAGainRatio gr;
gr.setGain12Over6( 1.9 + r*0.2 );
gr.setGain6Over1( 5.9 + r*0.2 );
gratio->setValue( ebid.rawId(), gr );
} // loop over phi
} // loop over eta
return gratio;
}
| EcalIntercalibConstants * EcalTestDevDB::generateEcalIntercalibConstants | ( | ) |
Definition at line 257 of file EcalTestDevDB.cc.
References EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, csvReporter::r, DetId::rawId(), and EcalCondObjectContainer< T >::setValue().
Referenced by analyze().
{
//-------------------------------------------------------------
EcalIntercalibConstants* ical = new EcalIntercalibConstants();
for(int ieta=-EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA; ++ieta) {
if(ieta==0) continue;
for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId:: MAX_IPHI; ++iphi) {
EBDetId ebid(ieta,iphi);
double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
ical->setValue( ebid.rawId(), 0.85 + r*0.3 );
} // loop over phi
} // loop over eta
return ical;
}
| EcalLaserAlphas * EcalTestDevDB::generateEcalLaserAlphas | ( | ) |
Definition at line 585 of file EcalTestDevDB.cc.
References EcalCondObjectContainer< T >::getMap(), EEDetId::hashedIndex(), EBDetId::hashedIndex(), EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_HASH, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalCondObjectContainer< T >::setValue(), and EEDetId::validDetId().
Referenced by analyze().
{
//--------------------------------------------------------------
EcalLaserAlphas* laser = new EcalLaserAlphas();
EcalLaserAlpha Alpha;
for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
if(iEta==0) continue;
for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
Alpha = double(1.55);
EBDetId ebid(iEta,iPhi);
int hi = ebid.hashedIndex();
if (hi< static_cast<int>(laser->getMap().size())) {
laser->setValue(hi, Alpha);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (alpha)!" << std::endl;
}
}
}
for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) {
for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) {
if (!EEDetId::validDetId(iX,iY,1))
continue;
EEDetId eedetidpos(iX,iY,1);
Alpha = double(1.55);
int hi = eedetidpos.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getMap().size())) {
laser->setValue(hi, Alpha);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (alpha)!" << std::endl;
}
if (!EEDetId::validDetId(iX,iY,-1))
continue;
EEDetId eedetidneg(iX,iY,-1);
Alpha = double(1.55);
hi = eedetidneg.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getMap().size())) {
laser->setValue(hi, Alpha);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (alpha)!" << std::endl;
}
}
}
return laser;
}
| EcalLaserAPDPNRatios * EcalTestDevDB::generateEcalLaserAPDPNRatios | ( | uint32_t | i_run | ) |
Definition at line 370 of file EcalTestDevDB.cc.
References gather_cfg::cout, funct::exp(), EcalLaserAPDPNRatios::getLaserMap(), EcalLaserAPDPNRatios::getTimeMap(), EEDetId::hashedIndex(), EBDetId::hashedIndex(), i, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, funct::log(), m_firstRun, m_lastRun, EBDetId::MAX_HASH, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p1, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p2, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p3, EcalLaserAPDPNRatios::setTime(), EcalLaserAPDPNRatios::setValue(), EcalLaserAPDPNRatios::EcalLaserTimeStamp::t1, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t2, and EEDetId::validDetId().
Referenced by analyze().
{
//--------------------------------------------------------------
EcalLaserAPDPNRatios* laser = new EcalLaserAPDPNRatios();
EcalLaserAPDPNRatios::EcalLaserAPDPNpair APDPNpair;
EcalLaserAPDPNRatios::EcalLaserTimeStamp TimeStamp;
// if((m_firstRun == 0 && i_run == 0) || (m_firstRun == 1 && i_run == 1)){
std::cout << "First & last run: " << i_run << " " << m_firstRun << " " << m_lastRun << " " << std::endl;
if (m_firstRun == i_run && (i_run == 0 || i_run == 1) ) {
APDPNpair.p1= (double(1)+1/double(log(exp(1)+double((i_run-m_firstRun)*10))))/double(2);
APDPNpair.p2= (double(1)+1/double(log(exp(1)+double((i_run-m_firstRun)*10)+double(10))))/double(2);
APDPNpair.p3= double(0);
std::cout << i_run << " " << m_firstRun << " " << APDPNpair.p1 << " " << APDPNpair.p2 << std::endl;
for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
if(iEta==0) continue;
for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
//APDPNpair.p1= double(1);
//APDPNpair.p2= double(1);
EBDetId ebid(iEta,iPhi);
int hi = ebid.hashedIndex();
if (hi< static_cast<int>(laser->getLaserMap().size())) {
laser->setValue(hi, APDPNpair);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ratio)!" << std::endl;
continue;
}
}
}
for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) {
for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) {
// make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
if (!EEDetId::validDetId(iX,iY,1))
continue;
EEDetId eedetidpos(iX,iY,1);
//APDPNpair.p1 = double(1);
//APDPNpair.p2 = double(1);
int hi = eedetidpos.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getLaserMap().size())) {
laser->setValue(hi, APDPNpair);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ratio)!" << std::endl;
continue;
}
if (!EEDetId::validDetId(iX,iY,-1))
continue;
EEDetId eedetidneg(iX,iY,-1);
//APDPNpair.p1 = double(1);
//APDPNpair.p2 = double(1);
hi = eedetidneg.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getLaserMap().size())) {
laser->setValue(hi, APDPNpair);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ratio)!" << std::endl;
continue;
}
}
}
std::cout << std::endl;
for(int i=0; i<92; i++){
if (i< static_cast<int>(laser->getTimeMap().size())) {
TimeStamp.t1 = edm::Timestamp(1380*(i_run-m_firstRun) + 15*i);
TimeStamp.t2 = edm::Timestamp(1380*(i_run-m_firstRun + 1) + 15*i);
laser->setTime(i, TimeStamp);
//std::cout << " Timestamp for " << i << " : " << TimeStamp.t1.value() << " , " << TimeStamp.t2.value() << std::endl;
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (time)!" << std::endl;
continue;
}
}
}else{
APDPNpair.p1= (double(1)+1/double(log(exp(1)+double((i_run-m_firstRun)*10))))/double(2);
APDPNpair.p2= (double(1)+1/double(log(exp(1)+double((i_run-m_firstRun)*10)+double(10))))/double(2);
APDPNpair.p3= double(0);
std::cout << i_run << " " << m_firstRun << " " << APDPNpair.p1 << " " << APDPNpair.p2 << std::endl;
for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
if(iEta==0) continue;
for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
EBDetId ebid(iEta,iPhi);
int hi = ebid.hashedIndex();
if (hi< static_cast<int>(laser->getLaserMap().size())) {
laser->setValue(hi, APDPNpair);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ratio)!" << std::endl;
}
}
}
for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) {
for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) {
// make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
if (!EEDetId::validDetId(iX,iY,1))
continue;
EEDetId eedetidpos(iX,iY,1);
int hi = eedetidpos.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getLaserMap().size())) {
laser->setValue(hi, APDPNpair);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ratio)!" << std::endl;
continue;
}
if (!EEDetId::validDetId(iX,iY,-1))
continue;
EEDetId eedetidneg(iX,iY,-1);
hi = eedetidneg.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getLaserMap().size())) {
laser->setValue(hi, APDPNpair);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ratio)!" << std::endl;
continue;
}
}
}
std::cout << std::endl;
for(int i=0; i<92; i++){
if (i< static_cast<int>(laser->getTimeMap().size())) {
TimeStamp.t1 = edm::Timestamp(1380*(i_run-m_firstRun) + 15*i);
TimeStamp.t2 = edm::Timestamp(1380*(i_run-m_firstRun + 1) + 15*i);
laser->setTime(i, TimeStamp);
//std::cout << " Timestamp for " << i << " : " << TimeStamp.t1.value() << " , " << TimeStamp.t2.value() << std::endl;
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (time)!" << std::endl;
continue;
}
}
}
return laser;
}
| EcalLaserAPDPNRatiosRef * EcalTestDevDB::generateEcalLaserAPDPNRatiosRef | ( | ) |
Definition at line 525 of file EcalTestDevDB.cc.
References EcalCondObjectContainer< T >::getMap(), EEDetId::hashedIndex(), EBDetId::hashedIndex(), EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_HASH, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalCondObjectContainer< T >::setValue(), and EEDetId::validDetId().
Referenced by analyze().
{
//--------------------------------------------------------------
EcalLaserAPDPNRatiosRef* laser = new EcalLaserAPDPNRatiosRef();
EcalLaserAPDPNref APDPNref;
for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
if(iEta==0) continue;
for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
APDPNref=double(1.5);
EBDetId ebid(iEta,iPhi);
int hi = ebid.hashedIndex();
if (hi< static_cast<int>(laser->getMap().size())) {
laser->setValue(hi, APDPNref);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ref)!" << std::endl;
}
}
}
for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) {
for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) {
if (!EEDetId::validDetId(iX,iY,1))
continue;
EEDetId eedetidpos(iX,iY,1);
APDPNref=double(1.5);
int hi = eedetidpos.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getMap().size())) {
laser->setValue(hi, APDPNref);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ref)!" << std::endl;
}
if (!EEDetId::validDetId(iX,iY,-1))
continue;
EEDetId eedetidneg(iX,iY,-1);
APDPNref=double(1.5);
hi = eedetidneg.hashedIndex() + EBDetId::MAX_HASH + 1;
if (hi< static_cast<int>(laser->getMap().size())) {
laser->setValue(hi, APDPNref);
} else {
edm::LogError("EcalTestDevDB") << "error with laser Map (ref)!" << std::endl;
}
}
}
return laser;
}
| EcalPedestals * EcalTestDevDB::generateEcalPedestals | ( | ) |
Definition at line 223 of file EcalTestDevDB.cc.
References EcalCondObjectContainer< T >::insert(), EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, and DetId::rawId().
Referenced by analyze().
{
//-------------------------------------------------------------
EcalPedestals* peds = new EcalPedestals();
EcalPedestals::Item item;
for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
if(iEta==0) continue;
for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
item.mean_x1 = 200.*( (double)std::rand()/(double(RAND_MAX)+double(1)) );
item.rms_x1 = (double)std::rand()/(double(RAND_MAX)+double(1));
item.mean_x6 = 1200.*( (double)std::rand()/(double(RAND_MAX)+double(1)) );
item.rms_x6 = 6.*( (double)std::rand()/(double(RAND_MAX)+double(1)) );
item.mean_x12 = 2400.*( (double)std::rand()/(double(RAND_MAX)+double(1)) );
item.rms_x12 = 12.*( (double)std::rand()/(double(RAND_MAX)+double(1)) );
EBDetId ebdetid(iEta,iPhi);
peds->insert(std::make_pair(ebdetid.rawId(),item));
}
}
return peds;
}
| EcalTBWeights * EcalTestDevDB::generateEcalTBWeights | ( | ) |
Definition at line 320 of file EcalTestDevDB.cc.
References EcalWeightSet::getChi2WeightsAfterGainSwitch(), EcalWeightSet::getChi2WeightsBeforeGainSwitch(), EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), i, j, EBDetId::MAX_IETA, csvReporter::r, and EcalTBWeights::setValue().
Referenced by analyze().
{
//-------------------------------------------------------------
EcalTBWeights* tbwgt = new EcalTBWeights();
// create weights for each distinct group ID
int nMaxTDC = 10;
for(int igrp=-EBDetId::MAX_IETA; igrp<=EBDetId::MAX_IETA; ++igrp) {
if(igrp==0) continue;
for(int itdc=1; itdc<=nMaxTDC; ++itdc) {
// generate random number
double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
// make a new set of weights
EcalWeightSet wgt;
EcalWeightSet::EcalWeightMatrix& mat1 = wgt.getWeightsBeforeGainSwitch();
EcalWeightSet::EcalWeightMatrix& mat2 = wgt.getWeightsAfterGainSwitch();
for(size_t i=0; i<3; ++i) {
for(size_t j=0; j<10; ++j) {
double ww = igrp*itdc*r + i*10. + j;
//std::cout << "row: " << i << " col: " << j << " - val: " << ww << std::endl;
mat1(i,j)=ww;
mat2(i,j)=100+ww;
}
}
// fill the chi2 matrcies
r = (double)std::rand()/( double(RAND_MAX)+double(1) );
EcalWeightSet::EcalChi2WeightMatrix& mat3 = wgt.getChi2WeightsBeforeGainSwitch();
EcalWeightSet::EcalChi2WeightMatrix& mat4 = wgt.getChi2WeightsAfterGainSwitch();
for(size_t i=0; i<10; ++i) {
for(size_t j=0; j<10; ++j) {
double ww = igrp*itdc*r + i*10. + j;
mat3(i,j)=1000+ww;
mat4(i,j)=1000+100+ww;
}
}
// put the weight in the container
tbwgt->setValue(std::make_pair(igrp,itdc), wgt);
}
}
return tbwgt;
}
| EcalWeightXtalGroups * EcalTestDevDB::generateEcalWeightXtalGroups | ( | ) |
Definition at line 304 of file EcalTestDevDB.cc.
References EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, DetId::rawId(), and EcalCondObjectContainer< T >::setValue().
Referenced by analyze().
{
//-------------------------------------------------------------
EcalWeightXtalGroups* xtalGroups = new EcalWeightXtalGroups();
for(int ieta=-EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA; ++ieta) {
if(ieta==0) continue;
for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId:: MAX_IPHI; ++iphi) {
EBDetId ebid(ieta,iphi);
xtalGroups->setValue(ebid.rawId(), EcalXtalGroupId(ieta) ); // define rings in eta
}
}
return xtalGroups;
}
std::map<std::string, unsigned long long> EcalTestDevDB::m_cacheIDs [private] |
Definition at line 60 of file EcalTestDevDB.h.
Referenced by EcalTestDevDB().
unsigned long EcalTestDevDB::m_firstRun [private] |
Definition at line 62 of file EcalTestDevDB.h.
Referenced by analyze(), EcalTestDevDB(), and generateEcalLaserAPDPNRatios().
unsigned int EcalTestDevDB::m_interval [private] |
Definition at line 64 of file EcalTestDevDB.h.
Referenced by analyze(), and EcalTestDevDB().
unsigned long EcalTestDevDB::m_lastRun [private] |
Definition at line 63 of file EcalTestDevDB.h.
Referenced by analyze(), EcalTestDevDB(), and generateEcalLaserAPDPNRatios().
std::map<std::string, std::string> EcalTestDevDB::m_records [private] |
Definition at line 61 of file EcalTestDevDB.h.
Referenced by analyze(), and EcalTestDevDB().
std::string EcalTestDevDB::m_timetype [private] |
Definition at line 59 of file EcalTestDevDB.h.
1.7.3