Go to the documentation of this file.00001 #include "CalibTracker/SiStripESProducers/interface/SiStripApvGainGenerator.h"
00002 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00003 #include <boost/cstdint.hpp>
00004 #include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
00005 #include "FWCore/ParameterSet/interface/FileInPath.h"
00006
00007 #include "CLHEP/Random/RandFlat.h"
00008 #include "CLHEP/Random/RandGauss.h"
00009
00010 SiStripApvGainGenerator::SiStripApvGainGenerator(const edm::ParameterSet& iConfig,const edm::ActivityRegistry& aReg):
00011 SiStripCondObjBuilderBase<SiStripApvGain>::SiStripCondObjBuilderBase(iConfig)
00012 {
00013 edm::LogInfo("SiStripApvGainGenerator") << "[SiStripApvGainGenerator::SiStripApvGainGenerator]";
00014 }
00015
00016 SiStripApvGainGenerator::~SiStripApvGainGenerator() {
00017 edm::LogInfo("SiStripApvGainGenerator") << "[SiStripApvGainGenerator::~SiStripApvGainGenerator]";
00018 }
00019
00020 void SiStripApvGainGenerator::createObject(){
00021
00022 obj_ = new SiStripApvGain();
00023
00024 std::string genMode = _pset.getParameter<std::string>("genMode");
00025
00026 double meanGain_=_pset.getParameter<double>("MeanGain");
00027 double sigmaGain_=_pset.getParameter<double>("SigmaGain");
00028 double minimumPosValue_=_pset.getParameter<double>("MinPositiveGain");
00029
00030 edm::FileInPath fp_ = _pset.getParameter<edm::FileInPath>("file");
00031 uint32_t printdebug_ = _pset.getUntrackedParameter<uint32_t>("printDebug", 5);
00032 uint32_t count=0;
00033
00034 SiStripDetInfoFileReader reader(fp_.fullPath());
00035
00036 const std::map<uint32_t, SiStripDetInfoFileReader::DetInfo > DetInfos = reader.getAllData();
00037 float gainValue;
00038 for(std::map<uint32_t, SiStripDetInfoFileReader::DetInfo >::const_iterator it = DetInfos.begin(); it != DetInfos.end(); it++){
00039
00040 std::vector<float> theSiStripVector;
00041 for(unsigned short j=0; j<it->second.nApvs; j++){
00042
00043 if(genMode=="default")
00044 gainValue=meanGain_;
00045 else if (genMode=="gaussian") {
00046 gainValue = CLHEP::RandGauss::shoot(meanGain_, sigmaGain_);
00047 if(gainValue<=minimumPosValue_) gainValue=minimumPosValue_;
00048 }
00049 else {
00050 LogDebug("SiStripApvGain") << "ERROR: wrong genMode specifier : " << genMode << ", please select one of \"default\" or \"gaussian\"" << std::endl;
00051 exit(1);
00052 }
00053
00054 if (count<printdebug_) {
00055 edm::LogInfo("SiStripApvGainGenerator") << "detid: " << it->first << " Apv: " << j << " gain: " << gainValue << std::endl;
00056 }
00057 theSiStripVector.push_back(gainValue);
00058 }
00059 count++;
00060 SiStripApvGain::Range range(theSiStripVector.begin(),theSiStripVector.end());
00061 if ( ! obj_->put(it->first,range) )
00062 edm::LogError("SiStripApvGainGenerator")<<" detid already exists"<<std::endl;
00063 }
00064 }
1.7.3