#include <SiStripGainRandomCalculator.h>

Inheritance diagram for SiStripGainRandomCalculator:

Public Member Functions
	SiStripGainRandomCalculator (const edm::ParameterSet &)

	~SiStripGainRandomCalculator ()

Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

virtual	~ConditionDBWriter ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Member Functions
void	algoAnalyze (const edm::Event &, const edm::EventSetup &)

SiStripApvGain *	getNewObject ()

Private Attributes
std::vector< std::pair < uint32_t, unsigned short > >	detid_apvs_

unsigned long long	m_cacheID_

double	meanGain_

double	minimumPosValue_

bool	printdebug_

double	sigmaGain_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from ConditionDBWriter< SiStripApvGain >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 27 of file SiStripGainRandomCalculator.h.

Constructor & Destructor Documentation

SiStripGainRandomCalculator::SiStripGainRandomCalculator ( const edm::ParameterSet & iConfig )

explicit

Definition at line 26 of file SiStripGainRandomCalculator.cc.

References detid_apvs_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), meanGain_, minimumPosValue_, printdebug_, and sigmaGain_.

                                                                                        : ConditionDBWriter<SiStripApvGain>(iConfig), m_cacheID_(0){
 
   
   edm::LogInfo("SiStripGainRandomCalculator::SiStripGainRandomCalculator");
 
 //   std::string Mode=iConfig.getParameter<std::string>("Mode");
 //   if (Mode==std::string("Gaussian")) GaussianMode_=true;
 //   else if (IOVMode==std::string("Constant")) ConstantMode_=true;
 //   else  edm::LogError("SiStripGainRandomCalculator::SiStripGainRandomCalculator(): ERROR - unknown generation mode...will not store anything on the DB") << std::endl;
 
   detid_apvs_.clear();
 
   meanGain_=iConfig.getParameter<double>("MeanGain");
   sigmaGain_=iConfig.getParameter<double>("SigmaGain");
   minimumPosValue_=iConfig.getParameter<double>("MinPositiveGain");
   printdebug_ = iConfig.getUntrackedParameter<bool>("printDebug", false);
 
 
 }

SiStripGainRandomCalculator::~SiStripGainRandomCalculator ( )

Definition at line 47 of file SiStripGainRandomCalculator.cc.

                                                          {
 
    edm::LogInfo("SiStripGainRandomCalculator::~SiStripGainRandomCalculator");
 }

Member Function Documentation

void SiStripGainRandomCalculator::algoAnalyze	(	const edm::Event &	event,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 54 of file SiStripGainRandomCalculator.cc.

References cond::rpcobgas::detid, detid_apvs_, cmsRelvalreport::exit, edm::EventSetup::get(), m_cacheID_, StripTopology::nstrips(), AlCaHLTBitMon_ParallelJobs::p, and printdebug_.

                                                                                                 {
 
 
   unsigned long long cacheID = iSetup.get<TrackerDigiGeometryRecord>().cacheIdentifier();
   
   if (m_cacheID_ != cacheID) {
     
     m_cacheID_ = cacheID; 
 
     edm::ESHandle<TrackerGeometry> pDD;
 
     iSetup.get<TrackerDigiGeometryRecord>().get( pDD );
     edm::LogInfo("SiStripGainRandomCalculator::algoAnalyze - got new geometry  ")<<std::endl;
 
     detid_apvs_.clear();
     
     edm::LogInfo("SiStripGainCalculator") <<" There are "<<pDD->detUnits().size() <<" detectors"<<std::endl;
     
     for(TrackerGeometry::DetUnitContainer::const_iterator it = pDD->detUnits().begin(); it != pDD->detUnits().end(); it++){
   
       if( dynamic_cast<StripGeomDetUnit*>((*it))!=0){
     uint32_t detid=((*it)->geographicalId()).rawId();            
     const StripTopology& p = dynamic_cast<StripGeomDetUnit*>((*it))->specificTopology();
     unsigned short NAPVs = p.nstrips()/128;
     if(NAPVs<1 || NAPVs>6 ) {
       edm::LogError("SiStripGainCalculator")<<" Problem with Number of strips in detector.. "<< p.nstrips() <<" Exiting program"<<endl;
       exit(1);
     }
     detid_apvs_.push_back( pair<uint32_t,unsigned short>(detid,NAPVs) );
     if (printdebug_)
       edm::LogInfo("SiStripGainCalculator")<< "detid " << detid << " apvs " << NAPVs;
       }
     }
   }
 
 
 }

SiStripApvGain * SiStripGainRandomCalculator::getNewObject ( )

privatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 93 of file SiStripGainRandomCalculator.cc.

References gather_cfg::cout, detid_apvs_, j, meanGain_, minimumPosValue_, getGTfromDQMFile::obj, printdebug_, SiStripApvGain::put(), and sigmaGain_.

                                                            {
 
   std::cout<<"SiStripGainRandomCalculator::getNewObject called"<<std::endl;
 
   SiStripApvGain * obj = new SiStripApvGain();
 
   for(std::vector< pair<uint32_t,unsigned short> >::const_iterator it = detid_apvs_.begin(); it != detid_apvs_.end(); it++){
     //Generate Gain for det detid
     std::vector<float> theSiStripVector;
     for(unsigned short j=0; j<it->second; j++){
       float gain;
 
       //      if(sigmaGain_/meanGain_ < 0.00001) gain = meanGain_;
       //      else{
       gain = CLHEP::RandGauss::shoot(meanGain_, sigmaGain_);
       if(gain<=minimumPosValue_) gain=minimumPosValue_;
       //      }
 
       if (printdebug_)
     edm::LogInfo("SiStripGainCalculator") << "detid " << it->first << " \t"
                           << " apv " << j << " \t"
                           << gain    << " \t" 
                           << std::endl;         
       theSiStripVector.push_back(gain);
     }
     
     
     SiStripApvGain::Range range(theSiStripVector.begin(),theSiStripVector.end());
     if ( ! obj->put(it->first,range) )
       edm::LogError("SiStripGainCalculator")<<"[SiStripGainCalculator::beginJob] detid already exists"<<std::endl;
 
   }
   
   return obj;
 
 }