#include <EcalTrivialConditionRetriever.h>

Inheritance diagram for EcalTrivialConditionRetriever:

Public Member Functions
	EcalTrivialConditionRetriever (const edm::ParameterSet &pset)

virtual std::auto_ptr < EcalChannelStatus >	getChannelStatusFromConfiguration (const EcalChannelStatusRcd &)

virtual std::auto_ptr < EcalIntercalibConstants >	getIntercalibConstantsFromConfiguration (const EcalIntercalibConstantsRcd &)

virtual std::auto_ptr < EcalIntercalibConstantsMC >	getIntercalibConstantsMCFromConfiguration (const EcalIntercalibConstantsMCRcd &)

virtual std::auto_ptr < EcalIntercalibErrors >	getIntercalibErrorsFromConfiguration (const EcalIntercalibErrorsRcd &)

virtual std::auto_ptr < EcalMappingElectronics >	getMappingFromConfiguration (const EcalMappingElectronicsRcd &)

virtual std::auto_ptr < EcalTimeCalibConstants >	getTimeCalibConstantsFromConfiguration (const EcalTimeCalibConstantsRcd &)

virtual std::auto_ptr < EcalTimeCalibErrors >	getTimeCalibErrorsFromConfiguration (const EcalTimeCalibErrorsRcd &)

virtual std::auto_ptr < EcalTPGCrystalStatus >	getTrgChannelStatusFromConfiguration (const EcalTPGCrystalStatusRcd &)

virtual std::auto_ptr < EcalADCToGeVConstant >	produceEcalADCToGeVConstant (const EcalADCToGeVConstantRcd &)

virtual std::auto_ptr< Alignments >	produceEcalAlignmentEB (const EBAlignmentRcd &)

virtual std::auto_ptr< Alignments >	produceEcalAlignmentEE (const EEAlignmentRcd &)

virtual std::auto_ptr< Alignments >	produceEcalAlignmentES (const ESAlignmentRcd &)

virtual std::auto_ptr < EcalChannelStatus >	produceEcalChannelStatus (const EcalChannelStatusRcd &)

virtual std::auto_ptr < EcalClusterCrackCorrParameters >	produceEcalClusterCrackCorrParameters (const EcalClusterCrackCorrParametersRcd &)

virtual std::auto_ptr < EcalClusterEnergyCorrectionObjectSpecificParameters >	produceEcalClusterEnergyCorrectionObjectSpecificParameters (const EcalClusterEnergyCorrectionObjectSpecificParametersRcd &)

virtual std::auto_ptr < EcalClusterEnergyCorrectionParameters >	produceEcalClusterEnergyCorrectionParameters (const EcalClusterEnergyCorrectionParametersRcd &)

virtual std::auto_ptr < EcalClusterEnergyUncertaintyParameters >	produceEcalClusterEnergyUncertaintyParameters (const EcalClusterEnergyUncertaintyParametersRcd &)

virtual std::auto_ptr < EcalClusterLocalContCorrParameters >	produceEcalClusterLocalContCorrParameters (const EcalClusterLocalContCorrParametersRcd &)

virtual std::auto_ptr < EcalDAQTowerStatus >	produceEcalDAQTowerStatus (const EcalDAQTowerStatusRcd &)

virtual std::auto_ptr < EcalDCSTowerStatus >	produceEcalDCSTowerStatus (const EcalDCSTowerStatusRcd &)

virtual std::auto_ptr < EcalDQMChannelStatus >	produceEcalDQMChannelStatus (const EcalDQMChannelStatusRcd &)

virtual std::auto_ptr < EcalDQMTowerStatus >	produceEcalDQMTowerStatus (const EcalDQMTowerStatusRcd &)

virtual std::auto_ptr < EcalGainRatios >	produceEcalGainRatios (const EcalGainRatiosRcd &)

virtual std::auto_ptr < EcalIntercalibConstants >	produceEcalIntercalibConstants (const EcalIntercalibConstantsRcd &)

virtual std::auto_ptr < EcalIntercalibConstantsMC >	produceEcalIntercalibConstantsMC (const EcalIntercalibConstantsMCRcd &)

virtual std::auto_ptr < EcalIntercalibErrors >	produceEcalIntercalibErrors (const EcalIntercalibErrorsRcd &)

virtual std::auto_ptr < EcalLaserAlphas >	produceEcalLaserAlphas (const EcalLaserAlphasRcd &)

virtual std::auto_ptr < EcalLaserAPDPNRatios >	produceEcalLaserAPDPNRatios (const EcalLaserAPDPNRatiosRcd &)

virtual std::auto_ptr < EcalLaserAPDPNRatiosRef >	produceEcalLaserAPDPNRatiosRef (const EcalLaserAPDPNRatiosRefRcd &)

virtual std::auto_ptr < EcalLinearCorrections >	produceEcalLinearCorrections (const EcalLinearCorrectionsRcd &)

virtual std::auto_ptr < EcalMappingElectronics >	produceEcalMappingElectronics (const EcalMappingElectronicsRcd &)

virtual std::auto_ptr < EcalPedestals >	produceEcalPedestals (const EcalPedestalsRcd &)

virtual std::auto_ptr < EcalSampleMask >	produceEcalSampleMask (const EcalSampleMaskRcd &)

virtual std::auto_ptr < EcalTBWeights >	produceEcalTBWeights (const EcalTBWeightsRcd &)

virtual std::auto_ptr < EcalTimeBiasCorrections >	produceEcalTimeBiasCorrections (const EcalTimeBiasCorrectionsRcd &)

virtual std::auto_ptr < EcalTimeCalibConstants >	produceEcalTimeCalibConstants (const EcalTimeCalibConstantsRcd &)

virtual std::auto_ptr < EcalTimeCalibErrors >	produceEcalTimeCalibErrors (const EcalTimeCalibErrorsRcd &)

virtual std::auto_ptr < EcalTimeOffsetConstant >	produceEcalTimeOffsetConstant (const EcalTimeOffsetConstantRcd &)

virtual std::auto_ptr < EcalTPGCrystalStatus >	produceEcalTrgChannelStatus (const EcalTPGCrystalStatusRcd &)

virtual std::auto_ptr < EcalWeightXtalGroups >	produceEcalWeightXtalGroups (const EcalWeightXtalGroupsRcd &)

virtual	~EcalTrivialConditionRetriever ()

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

virtual	~ESProducer ()

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

virtual void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval)
	overrides DataProxyProvider method More...

virtual	~ESProxyFactoryProducer ()

Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()

const ComponentDescription &	description () const

bool	isUsingRecord (const EventSetupRecordKey &) const

const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const

void	resetProxies (const EventSetupRecordKey &iRecordType)

void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~DataProxyProvider ()

Public Member Functions inherited from edm::EventSetupRecordIntervalFinder
const eventsetup::ComponentDescription &	descriptionForFinder () const

	EventSetupRecordIntervalFinder ()

std::set < eventsetup::EventSetupRecordKey >	findingForRecords () const

const ValidityInterval &	findIntervalFor (const eventsetup::EventSetupRecordKey &, const IOVSyncValue &)

void	setDescriptionForFinder (const eventsetup::ComponentDescription &iDescription)

virtual	~EventSetupRecordIntervalFinder ()

Protected Member Functions
virtual void	setIntervalFor (const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue &, edm::ValidityInterval &)

Protected Member Functions inherited from edm::ESProducer
template<typename T >
void	setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())

template<typename T >
void	setWhatProduced (T iThis, const char iLabel)

template<typename T >
void	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
void	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord , typename TArg >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())

Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::auto_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::auto_ptr< eventsetup::ProxyFactoryBase > &iFactory, const std::string &iLabel=std::string())

virtual void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList)
	override DataProxyProvider method More...

Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...

void	invalidateProxies (const EventSetupRecordKey &iRecordKey)

template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &)

Protected Member Functions inherited from edm::EventSetupRecordIntervalFinder
template<class T >
void	findingRecord ()

void	findingRecordWithKey (const eventsetup::EventSetupRecordKey &)

Private Member Functions
	EcalTrivialConditionRetriever (const EcalTrivialConditionRetriever &)

void	getWeightsFromConfiguration (const edm::ParameterSet &ps)

const EcalTrivialConditionRetriever &	operator= (const EcalTrivialConditionRetriever &)

Private Attributes
double	adcToGeVEBConstant_

double	adcToGeVEEConstant_

std::vector < ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > >	amplWeights_

std::vector < ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > >	amplWeightsAft_

std::string	amplWeightsAftFile_

std::string	amplWeightsFile_

std::string	channelStatusFile_

std::vector < EcalWeightSet::EcalChi2WeightMatrix >	chi2Matrix_

std::vector < EcalWeightSet::EcalChi2WeightMatrix >	chi2MatrixAft_

std::string	chi2MatrixAftFile_

std::string	chi2MatrixFile_

std::vector< double >	crackCorrParameters_

std::string	EBAlignmentFile_

std::string	EBLaserAlphaFile_

double	EBpedMeanX12_

double	EBpedMeanX1_

double	EBpedMeanX6_

double	EBpedRMSX12_

double	EBpedRMSX1_

double	EBpedRMSX6_

std::vector< double >	EBtimeCorrAmplitudeBins_

std::vector< double >	EBtimeCorrShiftBins_

std::string	EEAlignmentFile_

std::string	EELaserAlphaFile_

double	EEpedMeanX12_

double	EEpedMeanX1_

double	EEpedMeanX6_

double	EEpedRMSX12_

double	EEpedRMSX1_

double	EEpedRMSX6_

std::vector< double >	EEtimeCorrAmplitudeBins_

std::vector< double >	EEtimeCorrShiftBins_

std::vector< double >	energyCorrectionObjectSpecificParameters_

std::vector< double >	energyCorrectionParameters_

std::vector< double >	energyUncertaintyParameters_

std::string	ESAlignmentFile_

double	gainRatio12over6_

double	gainRatio6over1_

bool	getEBAlignmentFromFile_

bool	getEEAlignmentFromFile_

bool	getESAlignmentFromFile_

bool	getLaserAlphaFromFile_

bool	getWeightsFromFile_

double	instLumi_

double	intercalibConstantMean_

double	intercalibConstantMeanMC_

std::string	intercalibConstantsFile_

double	intercalibConstantSigma_

double	intercalibConstantSigmaMC_

std::string	intercalibConstantsMCFile_

double	intercalibErrorMean_

std::string	intercalibErrorsFile_

std::vector < ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > >	jittWeights_

std::vector < ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > >	jittWeightsAft_

std::string	jittWeightsAftFile_

std::string	jittWeightsFile_

double	laserAlphaMean_

double	laserAlphaSigma_

double	laserAPDPNMean_

double	laserAPDPNRefMean_

double	laserAPDPNRefSigma_

double	laserAPDPNSigma_

unsigned long	laserAPDPNTime1_

unsigned long	laserAPDPNTime2_

unsigned long	laserAPDPNTime3_

double	linCorrMean_

double	linCorrSigma_

std::string	linearCorrectionsFile_

unsigned long	linearTime1_

unsigned long	linearTime2_

unsigned long	linearTime3_

std::vector< double >	localContCorrParameters_

std::string	mappingFile_

int	nTDCbins_

std::vector < ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > >	pedWeights_

std::vector < ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > >	pedWeightsAft_

std::string	pedWeightsAftFile_

std::string	pedWeightsFile_

bool	producedEcalADCToGeVConstant_

bool	producedEcalAlignmentEB_

bool	producedEcalAlignmentEE_

bool	producedEcalAlignmentES_

bool	producedEcalChannelStatus_

bool	producedEcalClusterCrackCorrParameters_

bool	producedEcalClusterEnergyCorrectionObjectSpecificParameters_

bool	producedEcalClusterEnergyCorrectionParameters_

bool	producedEcalClusterEnergyUncertaintyParameters_

bool	producedEcalClusterLocalContCorrParameters_

bool	producedEcalDAQTowerStatus_

bool	producedEcalDCSTowerStatus_

bool	producedEcalDQMChannelStatus_

bool	producedEcalDQMTowerStatus_

bool	producedEcalGainRatios_

bool	producedEcalIntercalibConstants_

bool	producedEcalIntercalibConstantsMC_

bool	producedEcalIntercalibErrors_

bool	producedEcalLaserCorrection_

bool	producedEcalLinearCorrections_

bool	producedEcalMappingElectronics_

bool	producedEcalPedestals_

bool	producedEcalSampleMask_

bool	producedEcalTimeBiasCorrections_

bool	producedEcalTimeCalibConstants_

bool	producedEcalTimeCalibErrors_

bool	producedEcalTimeOffsetConstant_

bool	producedEcalTrgChannelStatus_

bool	producedEcalWeights_

unsigned int	sampleMaskEB_

unsigned int	sampleMaskEE_

double	timeCalibConstantMean_

std::string	timeCalibConstantsFile_

double	timeCalibConstantSigma_

double	timeCalibErrorMean_

std::string	timeCalibErrorsFile_

double	timeOffsetEBConstant_

double	timeOffsetEEConstant_

double	totLumi_

std::string	trgChannelStatusFile_

int	verbose_

bool	weightsForAsynchronousRunning_

Additional Inherited Members
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair < DataKey, boost::shared_ptr < DataProxy > > >	KeyedProxies

typedef std::vector < EventSetupRecordKey >	Keys

typedef std::map < EventSetupRecordKey, KeyedProxies >	RecordProxies

Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 115 of file EcalTrivialConditionRetriever.h.

Constructor & Destructor Documentation

EcalTrivialConditionRetriever::EcalTrivialConditionRetriever ( const edm::ParameterSet & pset )

Definition at line 30 of file EcalTrivialConditionRetriever.cc.

 {
   // initilize parameters used to produce cond DB objects
   totLumi_=ps.getUntrackedParameter<double>("TotLumi",0.0);
   instLumi_ = ps.getUntrackedParameter<double>("InstLumi",0.0);
 
   // initilize parameters used to produce cond DB objects
   adcToGeVEBConstant_ = ps.getUntrackedParameter<double>("adcToGeVEBConstant",0.035);
   adcToGeVEEConstant_ = ps.getUntrackedParameter<double>("adcToGeVEEConstant",0.060);
 
   intercalibConstantMeanMC_ = ps.getUntrackedParameter<double>("intercalibConstantMeanMC",1.0);
   intercalibConstantSigmaMC_ = ps.getUntrackedParameter<double>("intercalibConstantSigmaMC",0.0);
 
   linCorrMean_ = ps.getUntrackedParameter<double>("linCorrMean",1.0);
   linCorrSigma_ = ps.getUntrackedParameter<double>("linCorrSigma",0.0);
 
   linearTime1_ = (unsigned long)atoi( ps.getUntrackedParameter<std::string>("linearTime1","1").c_str());
   linearTime2_=  (unsigned long)atoi( ps.getUntrackedParameter<std::string>("linearTime2","0").c_str());
   linearTime3_=  (unsigned long)atoi( ps.getUntrackedParameter<std::string>("linearTime3","0").c_str());
 
 
   intercalibConstantMean_ = ps.getUntrackedParameter<double>("intercalibConstantMean",1.0);
   intercalibConstantSigma_ = ps.getUntrackedParameter<double>("intercalibConstantSigma",0.0);
   intercalibErrorMean_ = ps.getUntrackedParameter<double>("IntercalibErrorMean",0.0);
 
   timeCalibConstantMean_ = ps.getUntrackedParameter<double>("timeCalibConstantMean",0.0);
   timeCalibConstantSigma_ = ps.getUntrackedParameter<double>("timeCalibConstantSigma",0.0);
   timeCalibErrorMean_ = ps.getUntrackedParameter<double>("timeCalibErrorMean",0.0);
 
   timeOffsetEBConstant_ = ps.getUntrackedParameter<double>("timeOffsetEBConstant",0.0);
   timeOffsetEEConstant_ = ps.getUntrackedParameter<double>("timeOffsetEEConstant",0.0);
 
   laserAlphaMean_  = ps.getUntrackedParameter<double>("laserAlphaMean",1.55);
   laserAlphaSigma_ = ps.getUntrackedParameter<double>("laserAlphaSigma",0);
 
   laserAPDPNTime1_ = (unsigned long)atoi( ps.getUntrackedParameter<std::string>("laserAPDPNTime1","1").c_str());
   laserAPDPNTime2_= (unsigned long)atoi( ps.getUntrackedParameter<std::string>("laserAPDPNTime2","0").c_str());
   laserAPDPNTime3_= (unsigned long)atoi( ps.getUntrackedParameter<std::string>("laserAPDPNTime3","0").c_str());
 
   laserAPDPNRefMean_ = ps.getUntrackedParameter<double>("laserAPDPNRefMean",1.0);
   laserAPDPNRefSigma_ = ps.getUntrackedParameter<double>("laserAPDPNRefSigma",0.0);
 
   laserAPDPNMean_ = ps.getUntrackedParameter<double>("laserAPDPNMean",1.0);
   laserAPDPNSigma_ = ps.getUntrackedParameter<double>("laserAPDPNSigma",0.0);
 
   localContCorrParameters_ = ps.getUntrackedParameter< std::vector<double> >("localContCorrParameters", std::vector<double>(0) );
   crackCorrParameters_ = ps.getUntrackedParameter< std::vector<double> >("crackCorrParameters", std::vector<double>(0) );
   energyCorrectionParameters_ = ps.getUntrackedParameter< std::vector<double> >("energyCorrectionParameters", std::vector<double>(0) );
   energyUncertaintyParameters_ = ps.getUntrackedParameter< std::vector<double> >("energyUncertaintyParameters", std::vector<double>(0) );
   energyCorrectionObjectSpecificParameters_ = ps.getUntrackedParameter< std::vector<double> >("energyCorrectionObjectSpecificParameters", std::vector<double>(0) );
 
   EBpedMeanX12_ = ps.getUntrackedParameter<double>("EBpedMeanX12", 200.);
   EBpedRMSX12_  = ps.getUntrackedParameter<double>("EBpedRMSX12",  1.10);
   EBpedMeanX6_  = ps.getUntrackedParameter<double>("EBpedMeanX6",  200.);
   EBpedRMSX6_   = ps.getUntrackedParameter<double>("EBpedRMSX6",   0.90);
   EBpedMeanX1_  = ps.getUntrackedParameter<double>("EBpedMeanX1",  200.);
   EBpedRMSX1_   = ps.getUntrackedParameter<double>("EBpedRMSX1",   0.62);
 
   EEpedMeanX12_ = ps.getUntrackedParameter<double>("EEpedMeanX12", 200.);
   EEpedRMSX12_  = ps.getUntrackedParameter<double>("EEpedRMSX12",  2.50);
   EEpedMeanX6_  = ps.getUntrackedParameter<double>("EEpedMeanX6",  200.);
   EEpedRMSX6_   = ps.getUntrackedParameter<double>("EEpedRMSX6",   2.00);
   EEpedMeanX1_  = ps.getUntrackedParameter<double>("EEpedMeanX1",  200.);
   EEpedRMSX1_   = ps.getUntrackedParameter<double>("EEpedRMSX1",   1.40);
 
   gainRatio12over6_ = ps.getUntrackedParameter<double>("gainRatio12over6", 2.0);
   gainRatio6over1_  = ps.getUntrackedParameter<double>("gainRatio6over1",  6.0);
 
   getWeightsFromFile_ = ps.getUntrackedParameter<bool>("getWeightsFromFile",false);
 
   sampleMaskEB_ = ps.getUntrackedParameter<unsigned int>("sampleMaskEB",1023);
   sampleMaskEE_ = ps.getUntrackedParameter<unsigned int>("sampleMaskEB",1023);
 
   EBtimeCorrAmplitudeBins_ = ps.getUntrackedParameter< std::vector<double> >("EBtimeCorrAmplitudeBins", std::vector<double>() );
   EBtimeCorrShiftBins_ = ps.getUntrackedParameter< std::vector<double> >("EBtimeCorrShiftBins", std::vector<double>() );
   EEtimeCorrAmplitudeBins_ = ps.getUntrackedParameter< std::vector<double> >("EEtimeCorrAmplitudeBins", std::vector<double>() );
   EEtimeCorrShiftBins_ = ps.getUntrackedParameter< std::vector<double> >("EEtimeCorrShiftBins", std::vector<double>() );
 
   nTDCbins_ = 1;
 
   weightsForAsynchronousRunning_ = ps.getUntrackedParameter<bool>("weightsForTB",false);
 
   std::cout << " EcalTrivialConditionRetriever " << std::endl;
 
 
 
   if(totLumi_ > 0 ) {
 
     std::cout << " EcalTrivialConditionRetriever going to create conditions based on the damage deu to "<<totLumi_<<
       " fb-1 integrated luminosity" << std::endl;
 
   }
 
   if (weightsForAsynchronousRunning_)
     {
       getWeightsFromFile_ = true; //override user request 
       //nTDCbins_ = 25;
       nTDCbins_ = 50; //modif Alex-21-07-2006
     }
 
   std::string path="CalibCalorimetry/EcalTrivialCondModules/data/";
   std::string weightType;
   std::ostringstream str;
 
   if (!weightsForAsynchronousRunning_)
     str << "_CMS.txt" ;
   else
     str << "_TB.txt" ;
 
   weightType = str.str();
 
   amplWeightsFile_ = ps.getUntrackedParameter<std::string>("amplWeightsFile",path+"ampWeights"+weightType);
   amplWeightsAftFile_ = ps.getUntrackedParameter<std::string>("amplWeightsAftFile",path+"ampWeightsAfterGainSwitch"+weightType);
   pedWeightsFile_ = ps.getUntrackedParameter<std::string>("pedWeightsFile",path+"pedWeights"+weightType);
   pedWeightsAftFile_ = ps.getUntrackedParameter<std::string>("pedWeightsAftFile",path+"pedWeightsAfterGainSwitch"+weightType);
   jittWeightsFile_ = ps.getUntrackedParameter<std::string>("jittWeightsFile",path+"timeWeights"+weightType);
   jittWeightsAftFile_ = ps.getUntrackedParameter<std::string>("jittWeightsAftFile",path+"timeWeightsAfterGainSwitch"+weightType);
   chi2MatrixFile_ = ps.getUntrackedParameter<std::string>("chi2MatrixFile",path+"chi2Matrix"+weightType);
   chi2MatrixAftFile_ = ps.getUntrackedParameter<std::string>("chi2MatrixAftFile",path+"chi2MatrixAfterGainSwitch"+weightType);
 
   amplWeights_.resize(nTDCbins_);  
   amplWeightsAft_.resize(nTDCbins_); 
   pedWeights_.resize(nTDCbins_);  
   pedWeightsAft_.resize(nTDCbins_); 
   jittWeights_.resize(nTDCbins_);  
   jittWeightsAft_.resize(nTDCbins_); 
   chi2Matrix_.resize(nTDCbins_);
   chi2MatrixAft_.resize(nTDCbins_);
 
   // default weights for MGPA shape after pedestal subtraction
   getWeightsFromConfiguration(ps);
 
   producedEcalPedestals_ = ps.getUntrackedParameter<bool>("producedEcalPedestals",true);
   producedEcalWeights_ = ps.getUntrackedParameter<bool>("producedEcalWeights",true);
 
   producedEcalGainRatios_ = ps.getUntrackedParameter<bool>("producedEcalGainRatios",true);
   producedEcalADCToGeVConstant_ = ps.getUntrackedParameter<bool>("producedEcalADCToGeVConstant",true);
 
   producedEcalMappingElectronics_ = ps.getUntrackedParameter<bool>("producedEcalMappingElectronics",true);
   mappingFile_ = ps.getUntrackedParameter<std::string>("mappingFile","");
 
   if ( producedEcalMappingElectronics_ ) {
     if ( mappingFile_ != "" ) { // if file provided read channel map
       setWhatProduced( this, &EcalTrivialConditionRetriever::getMappingFromConfiguration );
     } else { 
       setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalMappingElectronics );
     }
     findingRecord<EcalMappingElectronicsRcd>();
   }
   // Alignment from file
   getEBAlignmentFromFile_ = ps.getUntrackedParameter<bool>("getEBAlignmentFromFile",false);
   if(getEBAlignmentFromFile_) {
     EBAlignmentFile_ = ps.getUntrackedParameter<std::string>("EBAlignmentFile",path+"EBAlignment.txt");
   }
 
   getEEAlignmentFromFile_ = ps.getUntrackedParameter<bool>("getEEAlignmentFromFile",false);
   if(getEEAlignmentFromFile_) {
     EEAlignmentFile_ = ps.getUntrackedParameter<std::string>("EEAlignmentFile",path+"EEAlignment.txt");
   }
 
   getESAlignmentFromFile_ = ps.getUntrackedParameter<bool>("getESAlignmentFromFile",false);
   if(getESAlignmentFromFile_)
     ESAlignmentFile_ = ps.getUntrackedParameter<std::string>("ESAlignmentFile",path+"ESAlignment.txt");
 
   verbose_ = ps.getUntrackedParameter<int>("verbose", 0);
 
   //Tell Producer what we produce
   //setWhatproduce(this);
   if (producedEcalPedestals_)
     setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalPedestals );
 
   if (producedEcalWeights_) {
       setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalWeightXtalGroups );
       setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalTBWeights );
     }
 
   if (producedEcalGainRatios_)
     setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalGainRatios );
 
   if (producedEcalADCToGeVConstant_)
     setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalADCToGeVConstant );
 
   // TimeOffsetConstant
   producedEcalTimeOffsetConstant_ = ps.getUntrackedParameter<bool>("producedEcalTimeOffsetConstant",true);
   //std::cout << " EcalTrivialConditionRetriever : producedEcalTimeOffsetConstant_" << producedEcalTimeOffsetConstant_ << std::endl;
   if (producedEcalTimeOffsetConstant_) {
     setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalTimeOffsetConstant );
     findingRecord<EcalTimeOffsetConstantRcd>();
   }
 
   // linear corrections
   producedEcalLinearCorrections_ = ps.getUntrackedParameter<bool>("producedEcalLinearCorrections",true);
   linearCorrectionsFile_ = ps.getUntrackedParameter<std::string>("linearCorrectionsFile","") ;
 
   if (producedEcalLinearCorrections_) { // user asks to produce constants
     if(linearCorrectionsFile_ != "") {  // if file provided read constants
       setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalLinearCorrections );
     } else { // set all constants to 1. or smear as specified by user
         setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalLinearCorrections ) ;
     }
     findingRecord<EcalLinearCorrectionsRcd> () ;
   }
 
 
 
   // intercalibration constants
   producedEcalIntercalibConstants_ = ps.getUntrackedParameter<bool>("producedEcalIntercalibConstants",true);
   intercalibConstantsFile_ = ps.getUntrackedParameter<std::string>("intercalibConstantsFile","") ;
 
   intercalibConstantsMCFile_ = ps.getUntrackedParameter<std::string>("intercalibConstantsMCFile","") ;
 
   if (producedEcalIntercalibConstants_) { // user asks to produce constants
     if(intercalibConstantsFile_ != "") {  // if file provided read constants
         setWhatProduced (this, &EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration ) ;
     } else { // set all constants to 1. or smear as specified by user
         setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalIntercalibConstants ) ;
     }
     findingRecord<EcalIntercalibConstantsRcd> () ;
   }
   // MC intercalibrations
   producedEcalIntercalibConstantsMC_ = ps.getUntrackedParameter<bool>("producedEcalIntercalibConstantsMC",true);
 
   if (producedEcalIntercalibConstantsMC_) { // user asks to produce constants
     if(intercalibConstantsMCFile_ != "") {  // if file provided read constants
         setWhatProduced (this, &EcalTrivialConditionRetriever::getIntercalibConstantsMCFromConfiguration ) ;
     } else { // set all constants to 1. or smear as specified by user
         setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalIntercalibConstantsMC ) ;
     }
     findingRecord<EcalIntercalibConstantsMCRcd> () ;
   }
 
   // intercalibration constants
   producedEcalIntercalibErrors_ = ps.getUntrackedParameter<bool>("producedEcalIntercalibErrors",true);
   intercalibErrorsFile_ = ps.getUntrackedParameter<std::string>("intercalibErrorsFile","") ;
 
   if (producedEcalIntercalibErrors_) { // user asks to produce constants
     if(intercalibErrorsFile_ != "") {  // if file provided read constants
         setWhatProduced (this, &EcalTrivialConditionRetriever::getIntercalibErrorsFromConfiguration ) ;
     } else { // set all constants to 1. or smear as specified by user
         setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalIntercalibErrors ) ;
     }
     findingRecord<EcalIntercalibErrorsRcd> () ;
   }
 
   // time calibration constants
   producedEcalTimeCalibConstants_ = ps.getUntrackedParameter<bool>("producedEcalTimeCalibConstants",true);
   timeCalibConstantsFile_ = ps.getUntrackedParameter<std::string>("timeCalibConstantsFile","") ;
 
   if (producedEcalTimeCalibConstants_) { // user asks to produce constants
     if(timeCalibConstantsFile_ != "") {  // if file provided read constants
         setWhatProduced (this, &EcalTrivialConditionRetriever::getTimeCalibConstantsFromConfiguration ) ;
     } else { // set all constants to 1. or smear as specified by user
         setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalTimeCalibConstants ) ;
     }
     findingRecord<EcalTimeCalibConstantsRcd> () ;
   }
 
   // time calibration constants
   producedEcalTimeCalibErrors_ = ps.getUntrackedParameter<bool>("producedEcalTimeCalibErrors",true);
   timeCalibErrorsFile_ = ps.getUntrackedParameter<std::string>("timeCalibErrorsFile","") ;
 
   if (producedEcalTimeCalibErrors_) { // user asks to produce constants
     if(timeCalibErrorsFile_ != "") {  // if file provided read constants
         setWhatProduced (this, &EcalTrivialConditionRetriever::getTimeCalibErrorsFromConfiguration ) ;
     } else { // set all constants to 1. or smear as specified by user
         setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalTimeCalibErrors ) ;
     }
     findingRecord<EcalTimeCalibErrorsRcd> () ;
   }
 
   // cluster corrections
   producedEcalClusterLocalContCorrParameters_ = ps.getUntrackedParameter<bool>("producedEcalClusterLocalContCorrParameters", false);
   producedEcalClusterCrackCorrParameters_ = ps.getUntrackedParameter<bool>("producedEcalClusterCrackCorrParameters", false);
   producedEcalClusterEnergyCorrectionParameters_ = ps.getUntrackedParameter<bool>("producedEcalClusterEnergyCorrectionParameters", false);
   producedEcalClusterEnergyUncertaintyParameters_ = ps.getUntrackedParameter<bool>("producedEcalClusterEnergyUncertaintyParameters", false);
   producedEcalClusterEnergyCorrectionObjectSpecificParameters_ = ps.getUntrackedParameter<bool>("producedEcalClusterEnergyCorrectionObjectSpecificParameters", false);
   if ( producedEcalClusterLocalContCorrParameters_ ) {
           setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalClusterLocalContCorrParameters );
           findingRecord<EcalClusterLocalContCorrParametersRcd>();
   }
   if ( producedEcalClusterCrackCorrParameters_ ) {
           setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalClusterCrackCorrParameters );
           findingRecord<EcalClusterCrackCorrParametersRcd>();
   }
   if ( producedEcalClusterEnergyCorrectionParameters_ ) {
           setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalClusterEnergyCorrectionParameters );
           findingRecord<EcalClusterEnergyCorrectionParametersRcd>();
   }
   if ( producedEcalClusterEnergyUncertaintyParameters_ ) {
           setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalClusterEnergyUncertaintyParameters );
           findingRecord<EcalClusterEnergyUncertaintyParametersRcd>();
   }
   if ( producedEcalClusterEnergyCorrectionObjectSpecificParameters_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalClusterEnergyCorrectionObjectSpecificParameters );
     findingRecord<EcalClusterEnergyCorrectionObjectSpecificParametersRcd>();
   }
 
   // laser correction
   producedEcalLaserCorrection_ = ps.getUntrackedParameter<bool>("producedEcalLaserCorrection",true);
   if (producedEcalLaserCorrection_) { // user asks to produce constants
     // set all constants to 1. or smear as specified by user
     setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalLaserAlphas ) ;
     findingRecord<EcalLaserAlphasRcd> () ;
     getLaserAlphaFromFile_ = ps.getUntrackedParameter<bool>("getLaserAlphaFromFile",false);
     std::cout << " getLaserAlphaFromFile_ " <<  getLaserAlphaFromFile_ << std::endl;
     if(getLaserAlphaFromFile_) {
       EBLaserAlphaFile_ = ps.getUntrackedParameter<std::string>("EBLaserAlphaFile",path+"EBLaserAlpha.txt");
       EELaserAlphaFile_ = ps.getUntrackedParameter<std::string>("EELaserAlphaFile",path+"EELaserAlpha.txt");
       std::cout << " EELaserAlphaFile_ " <<  EELaserAlphaFile_.c_str() << std::endl;
     }
     setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatiosRef ) ;
     findingRecord<EcalLaserAPDPNRatiosRefRcd> () ;
     setWhatProduced (this, &EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios) ;
     findingRecord<EcalLaserAPDPNRatiosRcd> () ;
   }
 
   // channel status
   producedEcalChannelStatus_ = ps.getUntrackedParameter<bool>("producedEcalChannelStatus",true);
   channelStatusFile_ = ps.getUntrackedParameter<std::string>("channelStatusFile","");
 
   if ( producedEcalChannelStatus_ ) {
           if ( channelStatusFile_ != "" ) { // if file provided read channel map
                   setWhatProduced( this, &EcalTrivialConditionRetriever::getChannelStatusFromConfiguration );
           } else { // set all channels to working -- FIXME might be changed
                   setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalChannelStatus );
           }
           findingRecord<EcalChannelStatusRcd>();
   }
   // DQM channel status
   producedEcalDQMChannelStatus_ = ps.getUntrackedParameter<bool>("producedEcalDQMChannelStatus",true);
   if ( producedEcalDQMChannelStatus_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalDQMChannelStatus );
     findingRecord<EcalDQMChannelStatusRcd>();
   }
   // DCS Tower status
   producedEcalDCSTowerStatus_ = ps.getUntrackedParameter<bool>("producedEcalDCSTowerStatus",true);
   if ( producedEcalDCSTowerStatus_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalDCSTowerStatus );
     findingRecord<EcalDCSTowerStatusRcd>();
   }
   // DAQ Tower status
   producedEcalDAQTowerStatus_ = ps.getUntrackedParameter<bool>("producedEcalDAQTowerStatus",true);
   if ( producedEcalDAQTowerStatus_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalDAQTowerStatus );
     findingRecord<EcalDAQTowerStatusRcd>();
   }
   // DQM Tower status
   producedEcalDQMTowerStatus_ = ps.getUntrackedParameter<bool>("producedEcalDQMTowerStatus",true);
   if ( producedEcalDQMTowerStatus_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalDQMTowerStatus );
     findingRecord<EcalDQMTowerStatusRcd>();
   }
 
   // trigger channel status
   producedEcalTrgChannelStatus_ = ps.getUntrackedParameter<bool>("producedEcalTrgChannelStatus",true);
   trgChannelStatusFile_ = ps.getUntrackedParameter<std::string>("trgChannelStatusFile","");
 
   if ( producedEcalTrgChannelStatus_ ) {
           if ( trgChannelStatusFile_ != "" ) { // if file provided read channel map
                   setWhatProduced( this, &EcalTrivialConditionRetriever::getTrgChannelStatusFromConfiguration );
           } else { // set all channels to working -- FIXME might be changed
                   setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalTrgChannelStatus );
           }
           findingRecord<EcalTPGCrystalStatusRcd>();
   }
 
   // Alignment
   producedEcalAlignmentEB_ = ps.getUntrackedParameter<bool>("producedEcalAlignmentEB",true);
   if ( producedEcalAlignmentEB_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalAlignmentEB );
     findingRecord<EBAlignmentRcd>();
   }
   producedEcalAlignmentEE_ = ps.getUntrackedParameter<bool>("producedEcalAlignmentEE",true);
   if ( producedEcalAlignmentEE_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalAlignmentEE );
     findingRecord<EEAlignmentRcd>();
   }
   producedEcalAlignmentES_ = ps.getUntrackedParameter<bool>("producedEcalAlignmentES",true);
   if ( producedEcalAlignmentES_ ) {
     setWhatProduced( this, &EcalTrivialConditionRetriever::produceEcalAlignmentES );
     findingRecord<ESAlignmentRcd>();
   }
   //Tell Finder what records we find
   if (producedEcalPedestals_)  findingRecord<EcalPedestalsRcd>();
 
   if (producedEcalWeights_) {
       findingRecord<EcalWeightXtalGroupsRcd>();
       findingRecord<EcalTBWeightsRcd>();
     }
 
   if (producedEcalGainRatios_)  findingRecord<EcalGainRatiosRcd>();
 
   if (producedEcalADCToGeVConstant_)  findingRecord<EcalADCToGeVConstantRcd>();
 
   producedEcalSampleMask_ = ps.getUntrackedParameter<bool>("producedEcalSampleMask",true);
   if (producedEcalSampleMask_) {
     setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalSampleMask );
     findingRecord<EcalSampleMaskRcd>();
   }
   producedEcalTimeBiasCorrections_ = ps.getUntrackedParameter<bool>("producedEcalTimeBiasCorrections", false);
   if (producedEcalTimeBiasCorrections_) {
     setWhatProduced(this, &EcalTrivialConditionRetriever::produceEcalTimeBiasCorrections );
     findingRecord<EcalTimeBiasCorrectionsRcd>();
   }
 }

EcalTrivialConditionRetriever::~EcalTrivialConditionRetriever ( )

virtual

Definition at line 436 of file EcalTrivialConditionRetriever.cc.

437 {

438 }

EcalTrivialConditionRetriever::EcalTrivialConditionRetriever ( const EcalTrivialConditionRetriever & )

private

Member Function Documentation

std::auto_ptr< EcalChannelStatus > EcalTrivialConditionRetriever::getChannelStatusFromConfiguration ( const EcalChannelStatusRcd & )

virtual

Definition at line 1901 of file EcalTrivialConditionRetriever.cc.

References edm::hlt::Exception, contentValuesFiles::fullPath, edm::FileInPath::fullPath(), EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, ntuplemaker::status, AlCaHLTBitMon_QueryRunRegistry::string, EBDetId::unhashIndex(), EEDetId::unhashIndex(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalChannelStatus> ecalStatus = std::auto_ptr<EcalChannelStatus>( new EcalChannelStatus() );
   
 
   // start by setting all statuses to 0
   
   // barrel
   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) {
       if (EBDetId::validDetId(ieta,iphi)) {
     EBDetId ebid(ieta,iphi);
     ecalStatus->setValue( ebid, 0 );
       }
     }
   }
   // endcap
   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)) {
     EEDetId eedetidpos(iX,iY,1);
     ecalStatus->setValue( eedetidpos, 0 );
       }
       if (EEDetId::validDetId(iX,iY,-1)) {
     EEDetId eedetidneg(iX,iY,-1);
     ecalStatus->setValue( eedetidneg, 0 );
       }
     }
   }
   
   
   
   // overwrite the statuses which are in the file
   
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading channel statuses from file " << edm::FileInPath(channelStatusFile_).fullPath().c_str() ;
   std::ifstream statusFile(edm::FileInPath(channelStatusFile_).fullPath().c_str());
   if ( !statusFile.good() ) {
     edm::LogError ("EcalTrivialConditionRetriever") 
       << "*** Problems opening file: " << channelStatusFile_ ;
     throw cms::Exception ("Cannot open ECAL channel status file") ;
   }
   
   std::string EcalSubDet;
   std::string str;
   int hashIndex(0);
   int status(0);
   
   while (!statusFile.eof()) 
     {
       statusFile >> EcalSubDet;
       if (EcalSubDet!=std::string("EB") && EcalSubDet!=std::string("EE"))
     {
       std::getline(statusFile,str);
       continue;
     }
       else
     {
       statusFile>> hashIndex >> status;
     }
       // std::cout << EcalSubDet << " " << hashIndex << " " << status;
       
       if(EcalSubDet==std::string("EB"))
     {
       EBDetId ebid = EBDetId::unhashIndex(hashIndex);
       ecalStatus->setValue( ebid, status );
     }
       else if(EcalSubDet==std::string("EE"))
     {
       EEDetId eedetid = EEDetId::unhashIndex(hashIndex);
       ecalStatus->setValue( eedetid, status );
     }
       else
     {
       edm::LogError ("EcalTrivialConditionRetriever") 
         << " *** " << EcalSubDet << " is neither EB nor EE ";
     }
       
     }
   // the file is supposed to be in the form  -- FIXME
   
   
   statusFile.close();
   return ecalStatus;
 }

std::auto_ptr< EcalIntercalibConstants > EcalTrivialConditionRetriever::getIntercalibConstantsFromConfiguration ( const EcalIntercalibConstantsRcd & )

virtual

Definition at line 2313 of file EcalTrivialConditionRetriever.cc.

References EBDetId::approxEta(), EnergyResolutionVsLumi::calcresolutitonConstantTerm(), calib, gather_cfg::cout, EcalCondObjectContainer< T >::end(), eta(), edm::hlt::Exception, EcalCondObjectContainer< T >::find(), EcalCondObjectContainer< T >::getMap(), h, i, cuy::ii, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, geometryCSVtoXML::line, fff_deleter::log, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MAX_SM, EBDetId::MIN_IPHI, EBDetId::MIN_SM, alignCSCRings::r, DetId::rawId(), EcalFloatCondObjectContainerXMLTranslator::readXML(), EnergyResolutionVsLumi::setInstLumi(), EnergyResolutionVsLumi::setLumi(), EcalCondObjectContainer< T >::setValue(), EBDetId::SMCRYSTALMODE, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalIntercalibConstants>  ical;
   //        std::auto_ptr<EcalIntercalibConstants>( new EcalIntercalibConstants() );
 
   // Read the values from a txt file
   // -------------------------------
 
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading intercalibration constants from file "
                                                 << intercalibConstantsFile_.c_str() ;
 
   if(intercalibConstantsFile_.find(".xml")!= std::string::npos) {
 
     std::cout<<"generating Intercalib from xml file"<<std::endl; 
   
     EcalCondHeader h;
     EcalIntercalibConstants * rcd = new EcalIntercalibConstants;
     EcalIntercalibConstantsXMLTranslator::readXML(intercalibConstantsFile_,h,*rcd);
     
 
     if(totLumi_ !=0 || instLumi_!=0) {
       std::cout<<"implementing ageing for intercalib"<<std::endl; 
 
       EcalIntercalibConstantsMC * rcdMC = new EcalIntercalibConstantsMC;
 
       if(intercalibConstantsMCFile_.find(".xml")!= std::string::npos) {
 
     std::cout<<"generating IntercalibMC from xml file"<<std::endl; 
   
     EcalCondHeader h;
     EcalIntercalibConstantsMCXMLTranslator::readXML(intercalibConstantsMCFile_,h,*rcdMC);
 
       } else {
     std::cout<<"please provide the xml file of the EcalIntercalibConstantsMC"<<std::endl;
       }
 
 
       TRandom3 * gRandom = new TRandom3();
 
       EnergyResolutionVsLumi ageing;
       ageing.setLumi(totLumi_);
       ageing.setInstLumi(instLumi_);
       
 
       const EcalIntercalibConstantMap& mymap= rcd->getMap();
       const EcalIntercalibConstantMCMap& mymapMC= rcdMC->getMap();
 
       for(int ieta=-EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA ;++ieta) {
     if(ieta==0) continue;
     
     double eta=EBDetId::approxEta(EBDetId(ieta,1));
     eta = fabs(eta);
     double constantTerm= ageing.calcresolutitonConstantTerm(eta);
     std::cout<<"EB at eta="<<eta<<" constant term is "<<constantTerm<<std::endl;
     
     for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId::MAX_IPHI; ++iphi) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
         {
           EBDetId ebid(ieta,iphi);
           EcalIntercalibConstants::const_iterator idref=mymap.find(ebid);
           EcalIntercalibConstant icalconstant=1;
           if(idref!=mymap.end())icalconstant=(*idref);
 
           EcalIntercalibConstantsMC::const_iterator idrefMC=mymapMC.find(ebid);
           EcalIntercalibConstantMC icalconstantMC=1;
           if(idrefMC!=mymapMC.end())icalconstantMC=(*idrefMC);
           
           double r = gRandom->Gaus(0,constantTerm); 
 
           if(iphi==10) std::cout<<"EB at eta="<<eta<<" IC="<<icalconstant<<" ICMC="<<icalconstantMC<<" smear="<<r<<" ";
 
           icalconstant = icalconstant + r*1.29*icalconstantMC;
           rcd->setValue( ebid.rawId(), icalconstant );
 
           if(iphi==10) std::cout<<"newIC="<<icalconstant<<std::endl;
 
           EcalIntercalibConstant icalconstant2=(*idref);
           if(icalconstant !=icalconstant2) std::cout<<">>>> error in smearing intercalib"<<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))
         {
           
           EEDetId eedetidpos(iX,iY,1);
           double eta= -log(tan(0.5*atan(sqrt((iX-50.0)*(iX-50.0)+(iY-50.0)*(iY-50.0))*2.98/328.)));
           eta = fabs(eta);
           double constantTerm=ageing.calcresolutitonConstantTerm(eta);
           if(iX==50) std::cout<<"EE at eta="<<eta<<" constant term is "<<constantTerm<<std::endl;
 
 
 
           
           EcalIntercalibConstants::const_iterator idref=mymap.find(eedetidpos);
           EcalIntercalibConstant icalconstant=1;
           if(idref!=mymap.end())icalconstant=(*idref);
 
           EcalIntercalibConstantsMC::const_iterator idrefMC=mymapMC.find(eedetidpos);
           EcalIntercalibConstantMC icalconstantMC=1;
           if(idrefMC!=mymapMC.end())icalconstantMC=(*idrefMC);
           
           double r = gRandom->Gaus(0,constantTerm); 
 
           if(iX==10) std::cout<<"EE at eta="<<eta<<" IC="<<icalconstant<<" ICMC="<<icalconstantMC<<" smear="<<r<<" ";
           icalconstant = icalconstant + r*1.29*icalconstantMC;
           rcd->setValue( eedetidpos.rawId(), icalconstant );
           if(iX==10) std::cout<<"newIC="<<icalconstant<<std::endl;
 
 
           
         }
       if(EEDetId::validDetId(iX,iY,-1))
         {
           
           EEDetId eedetidneg(iX,iY,-1);
           double eta= -log(tan(0.5*atan(sqrt((iX-50.0)*(iX-50.0)+(iY-50.0)*(iY-50.0))*2.98/328.)));
           eta = fabs(eta);
           double constantTerm=ageing.calcresolutitonConstantTerm(eta);
           EcalIntercalibConstant icalconstant=1;
 
 
           EcalIntercalibConstantsMC::const_iterator idrefMC=mymapMC.find(eedetidneg);
           EcalIntercalibConstantMC icalconstantMC=1;
           if(idrefMC!=mymapMC.end())icalconstantMC=(*idrefMC);
           
           double r = gRandom->Gaus(0,constantTerm); 
           icalconstant = icalconstant + r*1.29*icalconstantMC;
           rcd->setValue( eedetidneg.rawId(), icalconstant );
 
           
         }
     }
       }
       
       ical = std::auto_ptr<EcalIntercalibConstants> (rcd);
 
       delete gRandom;
 
     }
 
 
   } else {
 
     ical =std::auto_ptr<EcalIntercalibConstants>( new EcalIntercalibConstants() );
 
     FILE *inpFile ;
     inpFile = fopen (intercalibConstantsFile_.c_str (),"r") ;
     if (!inpFile) 
       {
     edm::LogError ("EcalTrivialConditionRetriever") 
       << "*** Can not open file: " << intercalibConstantsFile_ ;
     throw cms::Exception ("Cannot open inter-calibration coefficients txt file") ;
       }
     
     char line[256] ;
     std::ostringstream str ;
     fgets (line,255,inpFile) ;
     int sm_number=atoi (line) ;
     str << "sm: " << sm_number ;  
     
     fgets (line,255,inpFile) ;
     //int nevents=atoi (line) ; // not necessary here just for online conddb
     
     fgets (line,255,inpFile) ;
     std::string gen_tag = line ;
     str << "gen tag: " << gen_tag ;  // should I use this? 
     
     fgets (line,255,inpFile) ;
     std::string cali_method = line ;
     str << "cali method: " << cali_method << std::endl ; // not important 
     
     fgets (line,255,inpFile) ;
     std::string cali_version = line ;
     str << "cali version: " << cali_version << std::endl ; // not important 
     
     fgets (line,255,inpFile) ;
     std::string cali_type = line ;
     str << "cali type: " << cali_type ; // not important
     
     edm::LogInfo("EcalTrivialConditionRetriever")
       << "[PIETRO] Intercalibration file - " 
       << str.str () << std::endl ;
     
     float calib[1700]={1} ;
     int calib_status[1700]={0} ;
     
     int ii = 0 ;
     
     while (fgets (line,255,inpFile)) 
       {
     ii++;
     int dmy_num = 0 ;
     float dmy_calib = 0. ;
     float dmy_RMS = 0. ;
     int dmy_events = 0 ;
     int dmy_status = 0 ;
     sscanf (line, "%d %f %f %d %d", &dmy_num, &dmy_calib,
         &dmy_RMS, &dmy_events,
         &dmy_status) ;
     assert (dmy_num >= 1) ;
     assert (dmy_num <= 1700) ;
     calib[dmy_num-1] = dmy_calib ; 
     calib_status[dmy_num-1] = dmy_status ;
     
     //       edm::LogInfo ("EcalTrivialConditionRetriever")
     //                 << "[PIETRO] cry = " << dmy_num 
     //                 << " calib = " << calib[dmy_num-1] 
     //                 << " RMS = " << dmy_RMS
     //                 << " events = " << dmy_events
     //                 << " status = " << calib_status[dmy_num-1] 
     //                 << std::endl ;
       }
 
     fclose (inpFile) ;           // close inp. file
     edm::LogInfo ("EcalTrivialConditionRetriever") << "Read intercalibrations for " << ii << " xtals " ; 
     if (ii!=1700) edm::LogWarning ("StoreEcalCondition") 
       << "Some crystals missing, set to 1" << std::endl ;
     
     // Transfer the data to the inter-calibration coefficients container
     // -----------------------------------------------------------------
     
     // DB supermodule always set to 1 for the TestBeam FIXME
     int sm_db=1 ;
     // loop over channels 
     for (int i=0 ; i<1700 ; i++)
       {
     //if (EBDetId::validDetId(iEta,iPhi)) {
     // CANNOT be used -- validDetId only checks ETA PHI method
     // doing a check by hand, here is the only place in CMSSW 
     // outside TB code and EcalRawToDigi where needed
     // => no need for changing the DetId interface
     //
     // checking only sm_db -- guess can change with the above FIXME
     if (sm_db >= EBDetId::MIN_SM && sm_db <= EBDetId::MAX_SM) {
       EBDetId ebid (sm_db,i+1,EBDetId::SMCRYSTALMODE) ;
       if (calib_status[i]) ical->setValue (ebid.rawId (), calib[i]) ;
       else ical->setValue (ebid.rawId (), 1.) ;
     }
     //}
       } // loop over channels 
 
 
 
 
   }
 
     return ical;
 
 
 }

std::auto_ptr< EcalIntercalibConstantsMC > EcalTrivialConditionRetriever::getIntercalibConstantsMCFromConfiguration ( const EcalIntercalibConstantsMCRcd & )

virtual

Definition at line 2571 of file EcalTrivialConditionRetriever.cc.

References gather_cfg::cout, h, and EcalFloatCondObjectContainerXMLTranslator::readXML().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalIntercalibConstantsMC>  ical;
   //        std::auto_ptr<EcalIntercalibConstants>( new EcalIntercalibConstants() );
 
   // Read the values from a xml file
   // -------------------------------
 
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading intercalibration constants MC from file "
                                                 << intercalibConstantsMCFile_.c_str() ;
 
   if(intercalibConstantsMCFile_.find(".xml")!= std::string::npos) {
 
     std::cout<<"generating Intercalib MC from xml file"<<std::endl; 
   
     EcalCondHeader h;
     EcalIntercalibConstantsMC * rcd = new EcalIntercalibConstantsMC;
     EcalIntercalibConstantsMCXMLTranslator::readXML(intercalibConstantsMCFile_,h,*rcd);
 
       ical = std::auto_ptr<EcalIntercalibConstants> (rcd);
 
   } else {
 
     std::cout <<"ERROR>>> please provide a xml file"<<std::endl;
   }
 
 
   return ical;
 
 }

std::auto_ptr< EcalIntercalibErrors > EcalTrivialConditionRetriever::getIntercalibErrorsFromConfiguration ( const EcalIntercalibErrorsRcd & )

virtual

Definition at line 2605 of file EcalTrivialConditionRetriever.cc.

References calib, edm::hlt::Exception, i, cuy::ii, geometryCSVtoXML::line, EBDetId::MAX_SM, EBDetId::MIN_SM, DetId::rawId(), EBDetId::SMCRYSTALMODE, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalIntercalibErrors>  ical = 
       std::auto_ptr<EcalIntercalibErrors>( new EcalIntercalibErrors() );
 
   // Read the values from a txt file
   // -------------------------------
 
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading intercalibration constants from file "
                                                 << intercalibErrorsFile_.c_str() ;
 
   FILE *inpFile ;
   inpFile = fopen (intercalibErrorsFile_.c_str (),"r") ;
   if (!inpFile) 
     {
       edm::LogError ("EcalTrivialConditionRetriever") 
          << "*** Can not open file: " << intercalibErrorsFile_ ;
       throw cms::Exception ("Cannot open inter-calibration coefficients txt file") ;
     }
 
   char line[256] ;
   std::ostringstream str ;
   fgets (line,255,inpFile) ;
   int sm_number=atoi (line) ;
   str << "sm: " << sm_number ;  
 
   fgets (line,255,inpFile) ;
   //int nevents=atoi (line) ; // not necessary here just for online conddb
 
   fgets (line,255,inpFile) ;
   std::string gen_tag = line ;
   str << "gen tag: " << gen_tag ;  // should I use this? 
 
   fgets (line,255,inpFile) ;
   std::string cali_method = line ;
   str << "cali method: " << cali_method << std::endl ; // not important 
 
   fgets (line,255,inpFile) ;
   std::string cali_version = line ;
   str << "cali version: " << cali_version << std::endl ; // not important 
 
   fgets (line,255,inpFile) ;
   std::string cali_type = line ;
   str << "cali type: " << cali_type ; // not important
 
   edm::LogInfo("EcalTrivialConditionRetriever")
             << "[PIETRO] Intercalibration file - " 
             << str.str () << std::endl ;
 
   float calib[1700]={1} ;
   int calib_status[1700]={0} ;
 
   int ii = 0 ;
 
   while (fgets (line,255,inpFile)) 
     {
       ii++;
       int dmy_num = 0 ;
       float dmy_calib = 0. ;
       float dmy_RMS = 0. ;
       int dmy_events = 0 ;
       int dmy_status = 0 ;
       sscanf (line, "%d %f %f %d %d", &dmy_num, &dmy_calib,
                                       &dmy_RMS, &dmy_events,
                                       &dmy_status) ;
       assert (dmy_num >= 1) ;
       assert (dmy_num <= 1700) ;
       calib[dmy_num-1] = dmy_calib ; 
       calib_status[dmy_num-1] = dmy_status ;
 
 //       edm::LogInfo ("EcalTrivialConditionRetriever")
 //                 << "[PIETRO] cry = " << dmy_num 
 //                 << " calib = " << calib[dmy_num-1] 
 //                 << " RMS = " << dmy_RMS
 //                 << " events = " << dmy_events
 //                 << " status = " << calib_status[dmy_num-1] 
 //                 << std::endl ;
     }
 
   fclose (inpFile) ;           // close inp. file
   edm::LogInfo ("EcalTrivialConditionRetriever") << "Read intercalibrations for " << ii << " xtals " ; 
   if (ii!=1700) edm::LogWarning ("StoreEcalCondition") 
                 << "Some crystals missing, set to 1" << std::endl ;
 
   // Transfer the data to the inter-calibration coefficients container
   // -----------------------------------------------------------------
 
   // DB supermodule always set to 1 for the TestBeam FIXME
   int sm_db=1 ;
   // loop over channels 
   for (int i=0 ; i<1700 ; i++)
     {
       //if (EBDetId::validDetId(iEta,iPhi)) {
       // CANNOT be used -- validDetId only checks ETA PHI method
       // doing a check by hand, here is the only place in CMSSW 
       // outside TB code and EcalRawToDigi where needed
       // => no need for changing the DetId interface
       //
       // checking only sm_db -- guess can change with the above FIXME
       if (sm_db >= EBDetId::MIN_SM && sm_db <= EBDetId::MAX_SM) {
         EBDetId ebid (sm_db,i+1,EBDetId::SMCRYSTALMODE) ;
         if (calib_status[i]) ical->setValue (ebid.rawId (), calib[i]) ;
         else ical->setValue (ebid.rawId (), 1.) ;
       }
       //}
     } // loop over channels 
     
 //  edm::LogInfo ("EcalTrivialConditionRetriever") << "INTERCALIBRATION DONE" ; 
   return ical;
 
 
 
 
 
 }

std::auto_ptr< EcalMappingElectronics > EcalTrivialConditionRetriever::getMappingFromConfiguration ( const EcalMappingElectronicsRcd & )

virtual

Definition at line 2954 of file EcalTrivialConditionRetriever.cc.

References cond::rpcobgas::detid, EcalMappingElement::electronicsid, edm::hlt::Exception, f, contentValuesFiles::fullPath, edm::FileInPath::fullPath(), EcalElectronicsId::rawId(), EcalTriggerElectronicsId::rawId(), EcalMappingElement::triggerid, and EEDetId::XYMODE.

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalMappingElectronics> mapping = std::auto_ptr<EcalMappingElectronics>( new EcalMappingElectronics() );
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading mapping from file " << edm::FileInPath(mappingFile_).fullPath().c_str() ;
   
   std::ifstream f(edm::FileInPath(mappingFile_).fullPath().c_str());
   if (!f.good())
     {
       edm::LogError("EcalTrivialConditionRetriever") << "File not found";
       throw cms::Exception("FileNotFound");
     }
   
   // uint32_t detid, elecid, triggerid;
   
   int ix, iy, iz, CL;
   // int dccid, towerid, stripid, xtalid;
   // int tccid, tower, ipseudostrip, xtalinps;
   int dccid, towerid, pseudostrip_in_SC, xtal_in_pseudostrip;
   int tccid, tower, pseudostrip_in_TCC, pseudostrip_in_TT;
   
   while ( ! f.eof()) 
     {
       // f >> detid >> elecid >> triggerid; 
       f >> ix >> iy >> iz >> CL >> dccid >> towerid >> pseudostrip_in_SC >> xtal_in_pseudostrip >> tccid >> tower >> 
     pseudostrip_in_TCC >> pseudostrip_in_TT ;
       
 //       if (!EEDetId::validDetId(ix,iy,iz))
 //    continue;
       
       EEDetId detid(ix,iy,iz,EEDetId::XYMODE);
       // std::cout << " dcc tower ps_in_SC xtal_in_ps " << dccid << " " << towerid << " " << pseudostrip_in_SC << " " << xtal_in_pseudostrip << std::endl;
       EcalElectronicsId elecid(dccid,towerid, pseudostrip_in_SC, xtal_in_pseudostrip);
       // std::cout << " tcc tt ps_in_TT xtal_in_ps " << tccid << " " << tower << " " << pseudostrip_in_TT << " " << xtal_in_pseudostrip << std::endl;
       EcalTriggerElectronicsId triggerid(tccid, tower, pseudostrip_in_TT, xtal_in_pseudostrip);
       EcalMappingElement aElement;
       aElement.electronicsid = elecid.rawId();
       aElement.triggerid = triggerid.rawId();
       (*mapping).setValue(detid, aElement);
     }
   
   f.close();
   return mapping;
 }

std::auto_ptr< EcalTimeCalibConstants > EcalTrivialConditionRetriever::getTimeCalibConstantsFromConfiguration ( const EcalTimeCalibConstantsRcd & )

virtual

Definition at line 2726 of file EcalTrivialConditionRetriever.cc.

References calib, edm::hlt::Exception, i, cuy::ii, geometryCSVtoXML::line, EBDetId::MAX_SM, EBDetId::MIN_SM, DetId::rawId(), EBDetId::SMCRYSTALMODE, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalTimeCalibConstants>  ical = 
       std::auto_ptr<EcalTimeCalibConstants>( new EcalTimeCalibConstants() );
 
   // Read the values from a txt file
   // -------------------------------
 
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading time calibration constants from file "
                                                 << timeCalibConstantsFile_.c_str() ;
 
   FILE *inpFile ;
   inpFile = fopen (timeCalibConstantsFile_.c_str (),"r") ;
   if (!inpFile) 
     {
       edm::LogError ("EcalTrivialConditionRetriever") 
          << "*** Can not open file: " << timeCalibConstantsFile_ ;
       throw cms::Exception ("Cannot open inter-calibration coefficients txt file") ;
     }
 
   char line[256] ;
   std::ostringstream str ;
   fgets (line,255,inpFile) ;
   int sm_number=atoi (line) ;
   str << "sm: " << sm_number ;  
 
   fgets (line,255,inpFile) ;
   //int nevents=atoi (line) ; // not necessary here just for online conddb
 
   fgets (line,255,inpFile) ;
   std::string gen_tag = line ;
   str << "gen tag: " << gen_tag ;  // should I use this? 
 
   fgets (line,255,inpFile) ;
   std::string cali_method = line ;
   str << "cali method: " << cali_method << std::endl ; // not important 
 
   fgets (line,255,inpFile) ;
   std::string cali_version = line ;
   str << "cali version: " << cali_version << std::endl ; // not important 
 
   fgets (line,255,inpFile) ;
   std::string cali_type = line ;
   str << "cali type: " << cali_type ; // not important
 
   edm::LogInfo("EcalTrivialConditionRetriever")
             << "TimeCalibration file - " 
             << str.str () << std::endl ;
 
   float calib[1700]={1} ;
   int calib_status[1700]={0} ;
 
   int ii = 0 ;
 
   while (fgets (line,255,inpFile)) 
     {
       ii++;
       int dmy_num = 0 ;
       float dmy_calib = 0. ;
       float dmy_RMS = 0. ;
       int dmy_events = 0 ;
       int dmy_status = 0 ;
       sscanf (line, "%d %f %f %d %d", &dmy_num, &dmy_calib,
                                       &dmy_RMS, &dmy_events,
                                       &dmy_status) ;
       assert (dmy_num >= 1) ;
       assert (dmy_num <= 1700) ;
       calib[dmy_num-1] = dmy_calib ; 
       calib_status[dmy_num-1] = dmy_status ;
 
 //       edm::LogInfo ("EcalTrivialConditionRetriever")
 //                 << "[PIETRO] cry = " << dmy_num 
 //                 << " calib = " << calib[dmy_num-1] 
 //                 << " RMS = " << dmy_RMS
 //                 << " events = " << dmy_events
 //                 << " status = " << calib_status[dmy_num-1] 
 //                 << std::endl ;
     }
 
   fclose (inpFile) ;           // close inp. file
   edm::LogInfo ("EcalTrivialConditionRetriever") << "Read timeCalibrations for " << ii << " xtals " ; 
   if (ii!=1700) edm::LogWarning ("StoreEcalCondition") 
                 << "Some crystals missing, set to 1" << std::endl ;
 
   // Transfer the data to the inter-calibration coefficients container
   // -----------------------------------------------------------------
 
   // DB supermodule always set to 1 for the TestBeam FIXME
   int sm_db=1 ;
   // loop over channels 
   for (int i=0 ; i<1700 ; i++)
     {
       //if (EBDetId::validDetId(iEta,iPhi)) {
       // CANNOT be used -- validDetId only checks ETA PHI method
       // doing a check by hand, here is the only place in CMSSW 
       // outside TB code and EcalRawToDigi where needed
       // => no need for changing the DetId interface
       //
       // checking only sm_db -- guess can change with the above FIXME
       if (sm_db >= EBDetId::MIN_SM && sm_db <= EBDetId::MAX_SM) {
         EBDetId ebid (sm_db,i+1,EBDetId::SMCRYSTALMODE) ;
         if (calib_status[i]) ical->setValue (ebid.rawId (), calib[i]) ;
         else ical->setValue (ebid.rawId (), 1.) ;
       }
       //}
     } // loop over channels 
     
 //  edm::LogInfo ("EcalTrivialConditionRetriever") << "INTERCALIBRATION DONE" ; 
   return ical;
 }

std::auto_ptr< EcalTimeCalibErrors > EcalTrivialConditionRetriever::getTimeCalibErrorsFromConfiguration ( const EcalTimeCalibErrorsRcd & )

virtual

Definition at line 2840 of file EcalTrivialConditionRetriever.cc.

References calib, edm::hlt::Exception, i, cuy::ii, geometryCSVtoXML::line, EBDetId::MAX_SM, EBDetId::MIN_SM, DetId::rawId(), EBDetId::SMCRYSTALMODE, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalTimeCalibErrors>  ical = 
       std::auto_ptr<EcalTimeCalibErrors>( new EcalTimeCalibErrors() );
 
   // Read the values from a txt file
   // -------------------------------
 
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading timeCalibration constants from file "
                                                 << timeCalibErrorsFile_.c_str() ;
 
   FILE *inpFile ;
   inpFile = fopen (timeCalibErrorsFile_.c_str (),"r") ;
   if (!inpFile) 
     {
       edm::LogError ("EcalTrivialConditionRetriever") 
          << "*** Can not open file: " << timeCalibErrorsFile_ ;
       throw cms::Exception ("Cannot open inter-calibration coefficients txt file") ;
     }
 
   char line[256] ;
   std::ostringstream str ;
   fgets (line,255,inpFile) ;
   int sm_number=atoi (line) ;
   str << "sm: " << sm_number ;  
 
   fgets (line,255,inpFile) ;
   //int nevents=atoi (line) ; // not necessary here just for online conddb
 
   fgets (line,255,inpFile) ;
   std::string gen_tag = line ;
   str << "gen tag: " << gen_tag ;  // should I use this? 
 
   fgets (line,255,inpFile) ;
   std::string cali_method = line ;
   str << "cali method: " << cali_method << std::endl ; // not important 
 
   fgets (line,255,inpFile) ;
   std::string cali_version = line ;
   str << "cali version: " << cali_version << std::endl ; // not important 
 
   fgets (line,255,inpFile) ;
   std::string cali_type = line ;
   str << "cali type: " << cali_type ; // not important
 
   edm::LogInfo("EcalTrivialConditionRetriever")
             << "TimeCalibration file - " 
             << str.str () << std::endl ;
 
   float calib[1700]={1} ;
   int calib_status[1700]={0} ;
 
   int ii = 0 ;
 
   while (fgets (line,255,inpFile)) 
     {
       ii++;
       int dmy_num = 0 ;
       float dmy_calib = 0. ;
       float dmy_RMS = 0. ;
       int dmy_events = 0 ;
       int dmy_status = 0 ;
       sscanf (line, "%d %f %f %d %d", &dmy_num, &dmy_calib,
                                       &dmy_RMS, &dmy_events,
                                       &dmy_status) ;
       assert (dmy_num >= 1) ;
       assert (dmy_num <= 1700) ;
       calib[dmy_num-1] = dmy_calib ; 
       calib_status[dmy_num-1] = dmy_status ;
 
 //       edm::LogInfo ("EcalTrivialConditionRetriever")
 //                 << "[PIETRO] cry = " << dmy_num 
 //                 << " calib = " << calib[dmy_num-1] 
 //                 << " RMS = " << dmy_RMS
 //                 << " events = " << dmy_events
 //                 << " status = " << calib_status[dmy_num-1] 
 //                 << std::endl ;
     }
 
   fclose (inpFile) ;           // close inp. file
   edm::LogInfo ("EcalTrivialConditionRetriever") << "Read time calibrations for " << ii << " xtals " ; 
   if (ii!=1700) edm::LogWarning ("StoreEcalCondition") 
                 << "Some crystals missing, set to 1" << std::endl ;
 
   // Transfer the data to the inter-calibration coefficients container
   // -----------------------------------------------------------------
 
   // DB supermodule always set to 1 for the TestBeam FIXME
   int sm_db=1 ;
   // loop over channels 
   for (int i=0 ; i<1700 ; i++)
     {
       //if (EBDetId::validDetId(iEta,iPhi)) {
       // CANNOT be used -- validDetId only checks ETA PHI method
       // doing a check by hand, here is the only place in CMSSW 
       // outside TB code and EcalRawToDigi where needed
       // => no need for changing the DetId interface
       //
       // checking only sm_db -- guess can change with the above FIXME
       if (sm_db >= EBDetId::MIN_SM && sm_db <= EBDetId::MAX_SM) {
         EBDetId ebid (sm_db,i+1,EBDetId::SMCRYSTALMODE) ;
         if (calib_status[i]) ical->setValue (ebid.rawId (), calib[i]) ;
         else ical->setValue (ebid.rawId (), 1.) ;
       }
       //}
     } // loop over channels 
     
 //  edm::LogInfo ("EcalTrivialConditionRetriever") << "INTERCALIBRATION DONE" ; 
   return ical;
 }

std::auto_ptr< EcalTPGCrystalStatus > EcalTrivialConditionRetriever::getTrgChannelStatusFromConfiguration ( const EcalTPGCrystalStatusRcd & )

virtual

Definition at line 2187 of file EcalTrivialConditionRetriever.cc.

References edm::hlt::Exception, contentValuesFiles::fullPath, edm::FileInPath::fullPath(), EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, ntuplemaker::status, AlCaHLTBitMon_QueryRunRegistry::string, EBDetId::unhashIndex(), EEDetId::unhashIndex(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalTPGCrystalStatus> ecalStatus = std::auto_ptr<EcalTPGCrystalStatus>( new EcalTPGCrystalStatus() );
   
 
   // start by setting all statuses to 0
   
   // barrel
   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) {
       if (EBDetId::validDetId(ieta,iphi)) {
     EBDetId ebid(ieta,iphi);
     ecalStatus->setValue( ebid, 0 );
       }
     }
   }
   // endcap
   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)) {
     EEDetId eedetidpos(iX,iY,1);
     ecalStatus->setValue( eedetidpos, 0 );
       }
       if (EEDetId::validDetId(iX,iY,-1)) {
     EEDetId eedetidneg(iX,iY,-1);
     ecalStatus->setValue( eedetidneg, 0 );
       }
     }
   }
   
   
   
   // overwrite the statuses which are in the file
   
   edm::LogInfo("EcalTrivialConditionRetriever") << "Reading channel statuses from file " << edm::FileInPath(channelStatusFile_).fullPath().c_str() ;
   std::ifstream statusFile(edm::FileInPath(channelStatusFile_).fullPath().c_str());
   if ( !statusFile.good() ) {
     edm::LogError ("EcalTrivialConditionRetriever") 
       << "*** Problems opening file: " << channelStatusFile_ ;
     throw cms::Exception ("Cannot open ECAL channel status file") ;
   }
   
   std::string EcalSubDet;
   std::string str;
   int hashIndex(0);
   int status(0);
   
   while (!statusFile.eof()) 
     {
       statusFile >> EcalSubDet;
       if (EcalSubDet!=std::string("EB") && EcalSubDet!=std::string("EE"))
     {
       std::getline(statusFile,str);
       continue;
     }
       else
     {
       statusFile>> hashIndex >> status;
     }
       // std::cout << EcalSubDet << " " << hashIndex << " " << status;
       
       if(EcalSubDet==std::string("EB"))
     {
       EBDetId ebid = EBDetId::unhashIndex(hashIndex);
       ecalStatus->setValue( ebid, status );
     }
       else if(EcalSubDet==std::string("EE"))
     {
       EEDetId eedetid = EEDetId::unhashIndex(hashIndex);
       ecalStatus->setValue( eedetid, status );
     }
       else
     {
       edm::LogError ("EcalTrivialConditionRetriever") 
         << " *** " << EcalSubDet << " is neither EB nor EE ";
     }
       
     }
   // the file is supposed to be in the form  -- FIXME
   
   
   statusFile.close();
   return ecalStatus;
 }

void EcalTrivialConditionRetriever::getWeightsFromConfiguration ( const edm::ParameterSet & ps )

private

Definition at line 1293 of file EcalTrivialConditionRetriever.cc.

References begin, edm::hlt::Exception, contentValuesFiles::fullPath, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), i, j, ConfigFiles::l, and findQualityFiles::size.

 {
 
   std::vector < std::vector<double> > amplwgtv(nTDCbins_);
 
   if (!getWeightsFromFile_ && nTDCbins_ == 1)
     {
       std::vector<double> vampl;
       //As default using simple 3+1 weights
       vampl.push_back( -0.33333 );
       vampl.push_back( -0.33333 );
       vampl.push_back( -0.33333 );
       vampl.push_back(  0. );
       vampl.push_back(  0. );
       vampl.push_back(  1. );
       vampl.push_back(  0. );
       vampl.push_back(  0. );
       vampl.push_back(  0. );
       vampl.push_back(  0. );
       amplwgtv[0]= ps.getUntrackedParameter< std::vector<double> >("amplWeights", vampl);
     }
   else if (getWeightsFromFile_)
     {
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading amplitude weights from file " << edm::FileInPath(amplWeightsFile_).fullPath().c_str() ;
       std::ifstream amplFile(edm::FileInPath(amplWeightsFile_).fullPath().c_str());
       int tdcBin=0;
       while (!amplFile.eof() && tdcBin < nTDCbins_) 
     {
       for(int j = 0; j < 10; ++j) {
         float ww;
         amplFile >> ww;
         amplwgtv[tdcBin].push_back(ww);
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
       //Read from file
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
 
   
   for (int i=0;i<nTDCbins_;i++)
     {
       assert(amplwgtv[i].size() == 10);
       int j=0;
       for(std::vector<double>::const_iterator it = amplwgtv[i].begin(); it != amplwgtv[i].end(); ++it) 
     {
       (amplWeights_[i])[j]=*it;
       j++;
     }
     }
   
 
   std::vector < std::vector<double> > amplwgtvAftGain(nTDCbins_);
 
   if (!getWeightsFromFile_ && nTDCbins_ == 1 )
     {
       std::vector<double> vamplAftGain;
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  1. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       vamplAftGain.push_back(  0. );
       amplwgtvAftGain[0] = ps.getUntrackedParameter< std::vector<double> >("amplWeightsAftGain", vamplAftGain);
     }
   else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading amplitude weights aftre gain switch from file " << edm::FileInPath(amplWeightsAftFile_).fullPath().c_str() ;
       std::ifstream amplFile(edm::FileInPath(amplWeightsAftFile_).fullPath().c_str());
       int tdcBin=0;
       while (!amplFile.eof() && tdcBin < nTDCbins_) 
     {
       for(int j = 0; j < 10; ++j) {
         float ww;
         amplFile >> ww;
         amplwgtvAftGain[tdcBin].push_back(ww);
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
 
   for (int i=0;i<nTDCbins_;i++)
     {
       assert(amplwgtvAftGain[i].size() == 10);
       int j=0;
       for(std::vector<double>::const_iterator it = amplwgtvAftGain[i].begin(); it != amplwgtvAftGain[i].end(); ++it) {
     (amplWeightsAft_[i])[j]=*it;
     j++;
       }
     }
       
   // default weights to reco amplitude w/o pedestal subtraction
 
   std::vector< std::vector<double> > pedwgtv(nTDCbins_);
 
   if (!getWeightsFromFile_ && nTDCbins_ == 1)
     {
       std::vector<double> vped;
       vped.push_back( 0.33333 );
       vped.push_back( 0.33333 );
       vped.push_back( 0.33333 );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       pedwgtv[0] = ps.getUntrackedParameter< std::vector<double> >("pedWeights", vped);
     }
   else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading pedestal weights from file " << edm::FileInPath(pedWeightsFile_).fullPath().c_str() ;
       std::ifstream pedFile(edm::FileInPath(pedWeightsFile_).fullPath().c_str());
       int tdcBin=0;
       while (!pedFile.eof() && tdcBin < nTDCbins_) 
     {
       for(int j = 0; j < 10; ++j) {
         float ww;
         pedFile >> ww;
         pedwgtv[tdcBin].push_back(ww);
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
 
   for (int i=0;i<nTDCbins_;i++)
     {
       assert(pedwgtv[i].size() == 10);
       int j=0;
       for(std::vector<double>::const_iterator it = pedwgtv[i].begin(); it != pedwgtv[i].end(); ++it) {
     (pedWeights_[i])[j] = *it;
     j++;
       }
     }
   
   std::vector< std::vector<double> > pedwgtvaft(nTDCbins_);
 
   if (!getWeightsFromFile_ && nTDCbins_ == 1)
     {
       std::vector<double> vped;
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       vped.push_back( 0. );
       pedwgtvaft[0] = ps.getUntrackedParameter< std::vector<double> >("pedWeightsAft", vped);
     }
   else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading pedestal after gain switch weights from file " << edm::FileInPath(pedWeightsAftFile_).fullPath().c_str() ;
       std::ifstream pedFile(edm::FileInPath(pedWeightsAftFile_).fullPath().c_str());
       int tdcBin=0;
       while (!pedFile.eof() && tdcBin < nTDCbins_) 
     {
       for(int j = 0; j < 10; ++j) {
         float ww;
         pedFile >> ww;
         pedwgtvaft[tdcBin].push_back(ww);
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
   
   for (int i=0;i<nTDCbins_;i++)
     {
       assert(pedwgtvaft[i].size() == 10);
       int j=0;
       for(std::vector<double>::const_iterator it = pedwgtvaft[i].begin(); it != pedwgtvaft[i].end(); ++it) {
     (pedWeightsAft_[i])[j]=*it;
     j++;
       }
     }
   
   
   
   // default weights to reco jitter
 
   std::vector< std::vector<double> > jittwgtv(nTDCbins_);
 
   if (!getWeightsFromFile_  && nTDCbins_ == 1 )
     {
       std::vector<double> vjitt;
       vjitt.push_back( 0.04066309 );
       vjitt.push_back( 0.04066309 );
       vjitt.push_back( 0.04066309 );
       vjitt.push_back( 0.000 );
       vjitt.push_back( 1.325176 );
       vjitt.push_back( -0.04997078 );
       vjitt.push_back( -0.504338 );
       vjitt.push_back( -0.5024844 );
       vjitt.push_back( -0.3903718 );
       vjitt.push_back( 0.000 );
       jittwgtv[0] = ps.getUntrackedParameter< std::vector<double> >("jittWeights", vjitt);
     }
   else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading jitter weights from file " << edm::FileInPath(jittWeightsFile_).fullPath().c_str() ;
       std::ifstream jittFile(edm::FileInPath(jittWeightsFile_).fullPath().c_str());
       int tdcBin=0;
       while (!jittFile.eof() && tdcBin < nTDCbins_) 
     {
       for(int j = 0; j < 10; ++j) {
         float ww;
         jittFile >> ww;
         jittwgtv[tdcBin].push_back(ww);
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
   
   for (int i=0;i<nTDCbins_;i++)
     {
       assert(jittwgtv[i].size() == 10);
       int j=0;
       for(std::vector<double>::const_iterator it = jittwgtv[i].begin(); it != jittwgtv[i].end(); ++it) {
     (jittWeights_[i])[j]= *it;
     j++;
       }
     }
   
   std::vector< std::vector<double> > jittwgtvaft(nTDCbins_);
 
   if (!getWeightsFromFile_ && nTDCbins_ == 1)
     {
       std::vector<double> vjitt;
       vjitt.push_back( 0. );
       vjitt.push_back( 0. );
       vjitt.push_back( 0. );
       vjitt.push_back( 0. );
       vjitt.push_back( 1.097871 );
       vjitt.push_back( -0.04551035 );
       vjitt.push_back( -0.4159156 );
       vjitt.push_back( -0.4185352 );
       vjitt.push_back( -0.3367127 );
       vjitt.push_back( 0. );
       jittwgtvaft[0] = ps.getUntrackedParameter< std::vector<double> >("jittWeightsAft", vjitt);
     }
   else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading jitter after gain switch weights from file " << edm::FileInPath(jittWeightsAftFile_).fullPath().c_str() ;
       std::ifstream jittFile(edm::FileInPath(jittWeightsAftFile_).fullPath().c_str());
       int tdcBin=0;
       while (!jittFile.eof() && tdcBin < nTDCbins_) 
     {
       for(int j = 0; j < 10; ++j) {
         float ww;
         jittFile >> ww;
         jittwgtvaft[tdcBin].push_back(ww);
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
   
   for (int i=0;i<nTDCbins_;i++)
     {
       assert(jittwgtvaft[i].size() == 10);
       int j=0;
       for(std::vector<double>::const_iterator it = jittwgtvaft[i].begin(); it != jittwgtvaft[i].end(); ++it) {
     (jittWeightsAft_[i])[j]= *it;
     j++;
       }
     }
 
    
    std::vector<  EcalWeightSet::EcalChi2WeightMatrix > chi2Matrix(nTDCbins_);
    if (!getWeightsFromFile_  && nTDCbins_ == 1 )
      {
        //        chi2Matrix[0].resize(10);
        //        for (int i=0;i<10;i++)
      //  chi2Matrix[0][i].resize(10);
 
        chi2Matrix[0](0,0) = 0.694371;
        chi2Matrix[0](0,1) = -0.305629;  
        chi2Matrix[0](0,2) = -0.305629;
        chi2Matrix[0](0,3) = 0.;
        chi2Matrix[0](0,4) = 0.;
        chi2Matrix[0](0,5) = 0.;
        chi2Matrix[0](0,6) = 0.;
        chi2Matrix[0](0,7) = 0.;
        chi2Matrix[0](0,8) = 0.;
        chi2Matrix[0](0,9) = 0.;
        chi2Matrix[0](1,0) = -0.305629;
        chi2Matrix[0](1,1) = 0.694371;
        chi2Matrix[0](1,2) = -0.305629;
        chi2Matrix[0](1,3) = 0.;
        chi2Matrix[0](1,4) = 0.;
        chi2Matrix[0](1,5) = 0.;
        chi2Matrix[0](1,6) = 0.;
        chi2Matrix[0](1,7) = 0.;
        chi2Matrix[0](1,8) = 0.;
        chi2Matrix[0](1,9) = 0.;
        chi2Matrix[0](2,0) = -0.305629;
        chi2Matrix[0](2,1) = -0.305629;
        chi2Matrix[0](2,2) = 0.694371;
        chi2Matrix[0](2,3) = 0.;
        chi2Matrix[0](2,4) = 0.;
        chi2Matrix[0](2,5) = 0.;
        chi2Matrix[0](2,6) = 0.;
        chi2Matrix[0](2,7) = 0.;
        chi2Matrix[0](2,8) = 0.;
        chi2Matrix[0](2,9) = 0.;
        chi2Matrix[0](3,0) = 0.;
        chi2Matrix[0](3,1) = 0.;
        chi2Matrix[0](3,2) = 0.;
        chi2Matrix[0](3,3) = 0.;
        chi2Matrix[0](3,4) = 0.;
        chi2Matrix[0](3,5) = 0.;
        chi2Matrix[0](3,6) = 0.;
        chi2Matrix[0](3,7) = 0.;
        chi2Matrix[0](3,8) = 0.;
        chi2Matrix[0](3,9) = 0.; 
        chi2Matrix[0](4,0) = 0.;
        chi2Matrix[0](4,1) = 0.;
        chi2Matrix[0](4,2) = 0.;
        chi2Matrix[0](4,3) = 0.;
        chi2Matrix[0](4,4) = 0.8027116;
        chi2Matrix[0](4,5) = -0.2517103;
        chi2Matrix[0](4,6) = -0.2232882;
        chi2Matrix[0](4,7) = -0.1716192;
        chi2Matrix[0](4,8) = -0.1239006;
        chi2Matrix[0](4,9) = 0.; 
        chi2Matrix[0](5,0) = 0.;
        chi2Matrix[0](5,1) = 0.;
        chi2Matrix[0](5,2) = 0.;
        chi2Matrix[0](5,3) = 0.;
        chi2Matrix[0](5,4) = -0.2517103;
        chi2Matrix[0](5,5) = 0.6528964;
        chi2Matrix[0](5,6) = -0.2972839;
        chi2Matrix[0](5,7) = -0.2067162;
        chi2Matrix[0](5,8) = -0.1230729;
        chi2Matrix[0](5,9) = 0.;
        chi2Matrix[0](6,0) = 0.;
        chi2Matrix[0](6,1) = 0.;
        chi2Matrix[0](6,2) = 0.;
        chi2Matrix[0](6,3) = 0.;
        chi2Matrix[0](6,4) = -0.2232882;
        chi2Matrix[0](6,5) = -0.2972839;
        chi2Matrix[0](6,6) = 0.7413607;
        chi2Matrix[0](6,7) = -0.1883866;
        chi2Matrix[0](6,8) = -0.1235052;
        chi2Matrix[0](6,9) = 0.; 
        chi2Matrix[0](7,0) = 0.;
        chi2Matrix[0](7,1) = 0.;
        chi2Matrix[0](7,2) = 0.;
        chi2Matrix[0](7,3) = 0.;
        chi2Matrix[0](7,4) = -0.1716192;
        chi2Matrix[0](7,5) = -0.2067162;
        chi2Matrix[0](7,6) = -0.1883866;
        chi2Matrix[0](7,7) = 0.844935;
        chi2Matrix[0](7,8) = -0.124291;
        chi2Matrix[0](7,9) = 0.; 
        chi2Matrix[0](8,0) = 0.;
        chi2Matrix[0](8,1) = 0.;
        chi2Matrix[0](8,2) = 0.;
        chi2Matrix[0](8,3) = 0.;
        chi2Matrix[0](8,4) = -0.1239006;
        chi2Matrix[0](8,5) = -0.1230729;
        chi2Matrix[0](8,6) = -0.1235052;
        chi2Matrix[0](8,7) = -0.124291;
        chi2Matrix[0](8,8) = 0.8749833;
        chi2Matrix[0](8,9) = 0.;
        chi2Matrix[0](9,0) = 0.;
        chi2Matrix[0](9,1) = 0.;
        chi2Matrix[0](9,2) = 0.;
        chi2Matrix[0](9,3) = 0.;
        chi2Matrix[0](9,4) = 0.;
        chi2Matrix[0](9,5) = 0.;
        chi2Matrix[0](9,6) = 0.;
        chi2Matrix[0](9,7) = 0.;
        chi2Matrix[0](9,8) = 0.;
        chi2Matrix[0](9,9) = 0.; 
      }
    else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading chi2Matrix from file " << edm::FileInPath(chi2MatrixFile_).fullPath().c_str() ;
       std::ifstream chi2MatrixFile(edm::FileInPath(chi2MatrixFile_).fullPath().c_str());
       int tdcBin=0;
       while (!chi2MatrixFile.eof() && tdcBin < nTDCbins_) 
     {
       //      chi2Matrix[tdcBin].resize(10);
       for(int j = 0; j < 10; ++j) {
         for(int l = 0; l < 10; ++l) {
           float ww;
           chi2MatrixFile >> ww;
           chi2Matrix[tdcBin](j,l)=ww;
         }
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
    
 //    for (int i=0;i<nTDCbins_;i++)
 //      {
        //       assert(chi2Matrix[i].size() == 10);
    chi2Matrix_ =  chi2Matrix;
 //      }
 
    std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2MatrixAft(nTDCbins_);
    if (!getWeightsFromFile_  && nTDCbins_ == 1 )
      {
        //       chi2MatrixAft[0].resize(10);
        //       for (int i=0;i<10;i++)
      //  chi2MatrixAft[0][i].resize(10);
 
        chi2MatrixAft[0](0,0) = 0.;
        chi2MatrixAft[0](0,1) = 0.;  
        chi2MatrixAft[0](0,2) = 0.;
        chi2MatrixAft[0](0,3) = 0.;
        chi2MatrixAft[0](0,4) = 0.;
        chi2MatrixAft[0](0,5) = 0.;
        chi2MatrixAft[0](0,6) = 0.;
        chi2MatrixAft[0](0,7) = 0.;
        chi2MatrixAft[0](0,8) = 0.;
        chi2MatrixAft[0](0,9) = 0.;
        chi2MatrixAft[0](1,0) = 0.;
        chi2MatrixAft[0](1,1) = 0.;
        chi2MatrixAft[0](1,2) = 0.;
        chi2MatrixAft[0](1,3) = 0.;
        chi2MatrixAft[0](1,4) = 0.;
        chi2MatrixAft[0](1,5) = 0.;
        chi2MatrixAft[0](1,6) = 0.;
        chi2MatrixAft[0](1,7) = 0.;
        chi2MatrixAft[0](1,8) = 0.;
        chi2MatrixAft[0](1,9) = 0.;
        chi2MatrixAft[0](2,0) = 0.;
        chi2MatrixAft[0](2,1) = 0.;
        chi2MatrixAft[0](2,2) = 0.;
        chi2MatrixAft[0](2,3) = 0.;
        chi2MatrixAft[0](2,4) = 0.;
        chi2MatrixAft[0](2,5) = 0.;
        chi2MatrixAft[0](2,6) = 0.;
        chi2MatrixAft[0](2,7) = 0.;
        chi2MatrixAft[0](2,8) = 0.;
        chi2MatrixAft[0](2,9) = 0.;
        chi2MatrixAft[0](3,0) = 0.;
        chi2MatrixAft[0](3,1) = 0.;
        chi2MatrixAft[0](3,2) = 0.;
        chi2MatrixAft[0](3,3) = 0.;
        chi2MatrixAft[0](3,4) = 0.;
        chi2MatrixAft[0](3,5) = 0.;
        chi2MatrixAft[0](3,6) = 0.;
        chi2MatrixAft[0](3,7) = 0.;
        chi2MatrixAft[0](3,8) = 0.;
        chi2MatrixAft[0](3,9) = 0.; 
        chi2MatrixAft[0](4,0) = 0.;
        chi2MatrixAft[0](4,1) = 0.;
        chi2MatrixAft[0](4,2) = 0.;
        chi2MatrixAft[0](4,3) = 0.;
        chi2MatrixAft[0](4,4) = 0.8030884;
        chi2MatrixAft[0](4,5) = -0.2543541;
        chi2MatrixAft[0](4,6) = -0.2243544;
        chi2MatrixAft[0](4,7) = -0.1698177;
        chi2MatrixAft[0](4,8) = -0.1194506;
        chi2MatrixAft[0](4,9) = 0.; 
        chi2MatrixAft[0](5,0) = 0.;
        chi2MatrixAft[0](5,1) = 0.;
        chi2MatrixAft[0](5,2) = 0.;
        chi2MatrixAft[0](5,3) = 0.;
        chi2MatrixAft[0](5,4) = -0.2543541;
        chi2MatrixAft[0](5,5) = 0.6714465;
        chi2MatrixAft[0](5,6) = -0.2898025;
        chi2MatrixAft[0](5,7) = -0.2193564;
        chi2MatrixAft[0](5,8) = -0.1542964;
        chi2MatrixAft[0](5,9) = 0.;
        chi2MatrixAft[0](6,0) = 0.;
        chi2MatrixAft[0](6,1) = 0.;
        chi2MatrixAft[0](6,2) = 0.;
        chi2MatrixAft[0](6,3) = 0.;
        chi2MatrixAft[0](6,4) = -0.2243544;
        chi2MatrixAft[0](6,5) = -0.2898025;
        chi2MatrixAft[0](6,6) = 0.7443781;
        chi2MatrixAft[0](6,7) = -0.1934846;
        chi2MatrixAft[0](6,8) = -0.136098;
        chi2MatrixAft[0](6,9) = 0.; 
        chi2MatrixAft[0](7,0) = 0.;
        chi2MatrixAft[0](7,1) = 0.;
        chi2MatrixAft[0](7,2) = 0.;
        chi2MatrixAft[0](7,3) = 0.;
        chi2MatrixAft[0](7,4) = -0.1698177;
        chi2MatrixAft[0](7,5) = -0.2193564;
        chi2MatrixAft[0](7,6) = -0.1934846;
        chi2MatrixAft[0](7,7) = 0.8535482;
        chi2MatrixAft[0](7,8) = -0.1030149;
        chi2MatrixAft[0](7,9) = 0.; 
        chi2MatrixAft[0](8,0) = 0.;
        chi2MatrixAft[0](8,1) = 0.;
        chi2MatrixAft[0](8,2) = 0.;
        chi2MatrixAft[0](8,3) = 0.;
        chi2MatrixAft[0](8,4) = -0.1194506;
        chi2MatrixAft[0](8,5) = -0.1542964;
        chi2MatrixAft[0](8,6) = -0.136098;
        chi2MatrixAft[0](8,7) = -0.1030149;
        chi2MatrixAft[0](8,8) = 0.9275388;
        chi2MatrixAft[0](8,9) = 0.;
        chi2MatrixAft[0](9,0) = 0.;
        chi2MatrixAft[0](9,1) = 0.;
        chi2MatrixAft[0](9,2) = 0.;
        chi2MatrixAft[0](9,3) = 0.;
        chi2MatrixAft[0](9,4) = 0.;
        chi2MatrixAft[0](9,5) = 0.;
        chi2MatrixAft[0](9,6) = 0.;
        chi2MatrixAft[0](9,7) = 0.;
        chi2MatrixAft[0](9,8) = 0.;
        chi2MatrixAft[0](9,9) = 0.; 
      }
    else if (getWeightsFromFile_)
     {
       //Read from file
       edm::LogInfo("EcalTrivialConditionRetriever") << "Reading chi2MatrixAft from file " << edm::FileInPath(chi2MatrixAftFile_).fullPath().c_str() ;
       std::ifstream chi2MatrixAftFile(edm::FileInPath(chi2MatrixAftFile_).fullPath().c_str());
       int tdcBin=0;
       while (!chi2MatrixAftFile.eof() && tdcBin < nTDCbins_) 
     {
       //      chi2MatrixAft[tdcBin].resize(10);
       for(int j = 0; j < 10; ++j) {
         for(int l = 0; l < 10; ++l) {
           float ww;
           chi2MatrixAftFile >> ww;
           chi2MatrixAft[tdcBin](j,l)=ww;
         }
       }
       ++tdcBin;
     }
       assert (tdcBin == nTDCbins_);
     }
   else
     {
       //Not supported
       edm::LogError("EcalTrivialConditionRetriever") << "Configuration not supported. Exception is raised ";
       throw cms::Exception("WrongConfig");
     }
    
 //    for (int i=0;i<nTDCbins_;i++)
 //      {
        //       assert(chi2MatrixAft[i].size() == 10);
    chi2MatrixAft_ =  chi2MatrixAft;
    //      }
 
 }

const EcalTrivialConditionRetriever& EcalTrivialConditionRetriever::operator= ( const EcalTrivialConditionRetriever & )

private

std::auto_ptr< EcalADCToGeVConstant > EcalTrivialConditionRetriever::produceEcalADCToGeVConstant ( const EcalADCToGeVConstantRcd & )

virtual

Definition at line 878 of file EcalTrivialConditionRetriever.cc.

Referenced by EcalTrivialConditionRetriever().

 {
   return std::auto_ptr<EcalADCToGeVConstant>( new EcalADCToGeVConstant(adcToGeVEBConstant_,adcToGeVEEConstant_) );
 }

std::auto_ptr< Alignments > EcalTrivialConditionRetriever::produceEcalAlignmentEB ( const EBAlignmentRcd & )

virtual

Definition at line 3011 of file EcalTrivialConditionRetriever.cc.

References a, gather_cfg::cout, f, contentValuesFiles::fullPath, Alignments::m_align, and create_public_lumi_plots::transform.

Referenced by EcalTrivialConditionRetriever().

                                                                              {
   double mytrans[3] = {0., 0., 0.};
   double myeuler[3] = {0., 0., 0.};
   std::ifstream f;
   if(getEBAlignmentFromFile_)
     f.open(edm::FileInPath(EBAlignmentFile_).fullPath().c_str());
   std::vector<AlignTransform> my_align;
   int ieta = 1;
   int iphi = 0;
   for(int SM = 1 ; SM < 37; SM++ ) {
     // make an EBDetId since we need EBDetId::rawId()
     if(SM > 18) { 
       iphi = 1 + (SM - 19) * 20;
       ieta = -1;
     }
     else
       iphi = SM * 20;
     EBDetId ebdetId(ieta,iphi);
     if(getEBAlignmentFromFile_) {
       f >> myeuler[0] >> myeuler[1] >> myeuler[2] >> mytrans[0] >> mytrans[1] >> mytrans[2];
       std::cout << " translation " << mytrans[0] << " " << mytrans[1] << " " << mytrans[2] << "\n" 
        << " euler " << myeuler[0] << " " << myeuler[1] << " " << myeuler[2] << std::endl;
     }
     CLHEP::Hep3Vector translation( mytrans[0], mytrans[1], mytrans[2]);
     CLHEP::HepEulerAngles euler( myeuler[0], myeuler[1], myeuler[2]);
     AlignTransform transform( translation, euler, ebdetId );
     my_align.push_back(transform);
   }
   /*
   // check
   uint32_t m_rawid;
   double m_x, m_y, m_z;
   double m_phi, m_theta, m_psi;
   int SM = 0;
   for (std::vector<AlignTransform>::const_iterator i = my_align.begin(); i != my_align.end(); ++i){ 
     m_rawid = i->rawId();
     CLHEP::Hep3Vector translation = i->translation();
     m_x = translation.x();
     m_y = translation.y();
     m_z = translation.z();
     CLHEP::HepRotation rotation = i->rotation();
     m_phi = rotation.getPhi();
     m_theta = rotation.getTheta();
     m_psi = rotation.getPsi();
     SM++;
     std::cout << " SM " << SM << " Id " << m_rawid << " x " << m_x << " y " << m_y << " z " << m_z
           << " phi " << m_phi << " theta " << m_theta << " psi " << m_psi << std::endl;
   }
   */
   Alignments a; 
   a.m_align = my_align; 
 
   std::auto_ptr<Alignments> ical = std::auto_ptr<Alignments>( new Alignments(a) );
   return ical;
 }

std::auto_ptr< Alignments > EcalTrivialConditionRetriever::produceEcalAlignmentEE ( const EEAlignmentRcd & )

virtual

Definition at line 3068 of file EcalTrivialConditionRetriever.cc.

References a, gather_cfg::cout, f, contentValuesFiles::fullPath, Alignments::m_align, and create_public_lumi_plots::transform.

Referenced by EcalTrivialConditionRetriever().

                                                                              {
   double mytrans[3] = {0., 0., 0.};
   double myeuler[3] = {0., 0., 0.};
   std::ifstream f;
   if(getEEAlignmentFromFile_)
     f.open(edm::FileInPath(EEAlignmentFile_).fullPath().c_str());
   std::vector<AlignTransform> my_align;
   int ix = 20;
   int iy = 50;
   int side = -1;
   for(int Dee = 0 ; Dee < 4; Dee++ ) {
     // make an EEDetId since we need EEDetId::rawId()
     if(Dee == 1 || Dee == 3) 
       ix = 70;
     else ix = 20;
     if(Dee == 2)
       side = 1;
     EEDetId eedetId(ix, iy, side);
     if(getEEAlignmentFromFile_) {
       f >> myeuler[0] >> myeuler[1] >> myeuler[2] >> mytrans[0] >> mytrans[1] >> mytrans[2];
       std::cout << " translation " << mytrans[0] << " " << mytrans[1] << " " << mytrans[2] << "\n" 
        << " euler " << myeuler[0] << " " << myeuler[1] << " " << myeuler[2] << std::endl;
     }
     CLHEP::Hep3Vector translation( mytrans[0], mytrans[1], mytrans[2]);
     CLHEP::HepEulerAngles euler( myeuler[0], myeuler[1], myeuler[2]);
     AlignTransform transform( translation, euler, eedetId );
     my_align.push_back(transform);
   }
   Alignments a; 
   a.m_align = my_align; 
   std::auto_ptr<Alignments> ical = std::auto_ptr<Alignments>( new Alignments(a) );
   return ical;
 }

std::auto_ptr< Alignments > EcalTrivialConditionRetriever::produceEcalAlignmentES ( const ESAlignmentRcd & )

virtual

Definition at line 3103 of file EcalTrivialConditionRetriever.cc.

References a, gather_cfg::cout, f, contentValuesFiles::fullPath, Alignments::m_align, and create_public_lumi_plots::transform.

Referenced by EcalTrivialConditionRetriever().

                                                                              {
   double mytrans[3] = {0., 0., 0.};
   double myeuler[3] = {0., 0., 0.};
   std::ifstream f;
   if(getESAlignmentFromFile_)
     f.open(edm::FileInPath(ESAlignmentFile_).fullPath().c_str());
   std::vector<AlignTransform> my_align;
   //  int ix_vect[10] = {10, 30, 30, 50, 10, 30, 10, 30};
   int pl_vect[10] = {2, 2, 1, 1, 1, 1, 2, 2};
   int iy = 10;
   int side = -1;
   int strip = 1;
   for(int layer = 0 ; layer < 8; layer++ ) {
     // make an ESDetId since we need ESDetId::rawId()
     int ix = 10 + (layer%2) * 20;
     int plane = pl_vect[layer];
     if(layer > 3) side = 1;
     ESDetId esdetId(strip, ix, iy, plane, side);
     if(getESAlignmentFromFile_) {
       f >> myeuler[0] >> myeuler[1] >> myeuler[2] >> mytrans[0] >> mytrans[1] >> mytrans[2];
       std::cout << " translation " << mytrans[0] << " " << mytrans[1] << " " << mytrans[2] << "\n" 
         << " euler " << myeuler[0] << " " << myeuler[1] << " " << myeuler[2] << std::endl;
     }
     CLHEP::Hep3Vector translation( mytrans[0], mytrans[1], mytrans[2]);
     CLHEP::HepEulerAngles euler( myeuler[0], myeuler[1], myeuler[2]);
     AlignTransform transform( translation, euler, esdetId );
     my_align.push_back(transform);
   }
   Alignments a; 
   a.m_align = my_align; 
   std::auto_ptr<Alignments> ical = std::auto_ptr<Alignments>( new Alignments(a) );
   return ical;
 }

std::auto_ptr< EcalChannelStatus > EcalTrivialConditionRetriever::produceEcalChannelStatus ( const EcalChannelStatusRcd & )

virtual

Definition at line 1991 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
 
         std::auto_ptr<EcalChannelStatus>  ical = std::auto_ptr<EcalChannelStatus>( new EcalChannelStatus() );
         // barrel
         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) {
                         if (EBDetId::validDetId(ieta,iphi)) {
                                 EBDetId ebid(ieta,iphi);
                                 ical->setValue( ebid, 0 );
                         }
                 }
         }
         // endcap
         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)) {
                                 EEDetId eedetidpos(iX,iY,1);
                                 ical->setValue( eedetidpos, 0 );
                         }
                         if (EEDetId::validDetId(iX,iY,-1)) {
                                 EEDetId eedetidneg(iX,iY,-1);
                                 ical->setValue( eedetidneg, 0 );
                         }
                 }
         }
         return ical;
 }

std::auto_ptr< EcalClusterCrackCorrParameters > EcalTrivialConditionRetriever::produceEcalClusterCrackCorrParameters ( const EcalClusterCrackCorrParametersRcd & )

virtual

Definition at line 982 of file EcalTrivialConditionRetriever.cc.

References i.

Referenced by EcalTrivialConditionRetriever().

 {
         std::auto_ptr<EcalClusterCrackCorrParameters> ipar = std::auto_ptr<EcalClusterCrackCorrParameters>( new EcalClusterCrackCorrParameters() );
         for (size_t i = 0; i < crackCorrParameters_.size(); ++i ) {
                 ipar->params().push_back( crackCorrParameters_[i] );
         }
         return ipar;
 }

std::auto_ptr< EcalClusterEnergyCorrectionObjectSpecificParameters > EcalTrivialConditionRetriever::produceEcalClusterEnergyCorrectionObjectSpecificParameters ( const EcalClusterEnergyCorrectionObjectSpecificParametersRcd & )

virtual

Definition at line 1009 of file EcalTrivialConditionRetriever.cc.

References i.

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalClusterEnergyCorrectionObjectSpecificParameters> ipar = std::auto_ptr<EcalClusterEnergyCorrectionObjectSpecificParameters>( new EcalClusterEnergyCorrectionObjectSpecificParameters() );
   for (size_t i = 0; i < energyCorrectionObjectSpecificParameters_.size(); ++i ) {
     ipar->params().push_back( energyCorrectionObjectSpecificParameters_[i] );
   }
   return ipar;
 }

std::auto_ptr< EcalClusterEnergyCorrectionParameters > EcalTrivialConditionRetriever::produceEcalClusterEnergyCorrectionParameters ( const EcalClusterEnergyCorrectionParametersRcd & )

virtual

Definition at line 991 of file EcalTrivialConditionRetriever.cc.

References i.

Referenced by EcalTrivialConditionRetriever().

 {
         std::auto_ptr<EcalClusterEnergyCorrectionParameters> ipar = std::auto_ptr<EcalClusterEnergyCorrectionParameters>( new EcalClusterEnergyCorrectionParameters() );
         for (size_t i = 0; i < energyCorrectionParameters_.size(); ++i ) {
                 ipar->params().push_back( energyCorrectionParameters_[i] );
         }
         return ipar;
 }

std::auto_ptr< EcalClusterEnergyUncertaintyParameters > EcalTrivialConditionRetriever::produceEcalClusterEnergyUncertaintyParameters ( const EcalClusterEnergyUncertaintyParametersRcd & )

virtual

Definition at line 1000 of file EcalTrivialConditionRetriever.cc.

References i.

Referenced by EcalTrivialConditionRetriever().

 {
         std::auto_ptr<EcalClusterEnergyUncertaintyParameters> ipar = std::auto_ptr<EcalClusterEnergyUncertaintyParameters>( new EcalClusterEnergyUncertaintyParameters() );
         for (size_t i = 0; i < energyUncertaintyParameters_.size(); ++i ) {
                 ipar->params().push_back( energyUncertaintyParameters_[i] );
         }
         return ipar;
 }

std::auto_ptr< EcalClusterLocalContCorrParameters > EcalTrivialConditionRetriever::produceEcalClusterLocalContCorrParameters ( const EcalClusterLocalContCorrParametersRcd & )

virtual

Definition at line 973 of file EcalTrivialConditionRetriever.cc.

References i.

Referenced by EcalTrivialConditionRetriever().

 {
         std::auto_ptr<EcalClusterLocalContCorrParameters> ipar = std::auto_ptr<EcalClusterLocalContCorrParameters>( new EcalClusterLocalContCorrParameters() );
         for (size_t i = 0; i < localContCorrParameters_.size(); ++i ) {
                 ipar->params().push_back( localContCorrParameters_[i] );
         }
         return ipar;
 }

std::auto_ptr< EcalDAQTowerStatus > EcalTrivialConditionRetriever::produceEcalDAQTowerStatus ( const EcalDAQTowerStatusRcd & )

virtual

Definition at line 2144 of file EcalTrivialConditionRetriever.cc.

References EcalBarrel, i, j, gen::k, ntuplemaker::status, EcalTrigTowerDetId::validDetId(), and EcalScDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
 
         std::auto_ptr<EcalDAQTowerStatus>  ical = std::auto_ptr<EcalDAQTowerStatus>( new EcalDAQTowerStatus() );
 
     int status(0);
     
         // barrel
     int iz=0;
         for(int k=0 ; k<2; k++ ) {
       if(k==0) iz=-1;
       if(k==1) iz=+1;
       for(int i=1 ; i<73; i++) {
         for(int j=1 ; j<18; j++) {
           if (EcalTrigTowerDetId::validDetId(iz,EcalBarrel,j,i )){
         EcalTrigTowerDetId ebid(iz,EcalBarrel,j,i);
 
         ical->setValue( ebid, status );
           }
         }
       }
     }
 
 
         // endcap
         for(int k=0 ; k<2; k++ ) {
       if(k==0) iz=-1;
       if(k==1) iz=+1;
       for(int i=1 ; i<21; i++) {
         for(int j=1 ; j<21; j++) {
           if (EcalScDetId::validDetId(i,j,iz )){
         EcalScDetId eeid(i,j,iz);
         ical->setValue( eeid, status );
           }
         }
       }
     }
 
         return ical;
 }

std::auto_ptr< EcalDCSTowerStatus > EcalTrivialConditionRetriever::produceEcalDCSTowerStatus ( const EcalDCSTowerStatusRcd & )

virtual

Definition at line 2101 of file EcalTrivialConditionRetriever.cc.

References EcalBarrel, i, j, gen::k, ntuplemaker::status, EcalTrigTowerDetId::validDetId(), and EcalScDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
 
         std::auto_ptr<EcalDCSTowerStatus>  ical = std::auto_ptr<EcalDCSTowerStatus>( new EcalDCSTowerStatus() );
 
     int status(0);
     
         // barrel
     int iz=0;
         for(int k=0 ; k<2; k++ ) {
       if(k==0) iz=-1;
       if(k==1) iz=+1;
       for(int i=1 ; i<73; i++) {
         for(int j=1 ; j<18; j++) {
           if (EcalTrigTowerDetId::validDetId(iz,EcalBarrel,j,i )){
         EcalTrigTowerDetId ebid(iz,EcalBarrel,j,i);
 
         ical->setValue( ebid, status );
           }
         }
       }
     }
 
 
         // endcap
         for(int k=0 ; k<2; k++ ) {
       if(k==0) iz=-1;
       if(k==1) iz=+1;
       for(int i=1 ; i<21; i++) {
         for(int j=1 ; j<21; j++) {
           if (EcalScDetId::validDetId(i,j,iz )){
         EcalScDetId eeid(i,j,iz);
         ical->setValue( eeid, status );
           }
         }
       }
     }
 
         return ical;
 }

std::auto_ptr< EcalDQMChannelStatus > EcalTrivialConditionRetriever::produceEcalDQMChannelStatus ( const EcalDQMChannelStatusRcd & )

virtual

Definition at line 2024 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   uint32_t sta(0);
 
         std::auto_ptr<EcalDQMChannelStatus>  ical = std::auto_ptr<EcalDQMChannelStatus>( new EcalDQMChannelStatus() );
         // barrel
         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) {
                         if (EBDetId::validDetId(ieta,iphi)) {
                                 EBDetId ebid(ieta,iphi);
                                 ical->setValue( ebid, sta );
                         }
                 }
         }
         // endcap
         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)) {
                                 EEDetId eedetidpos(iX,iY,1);
                                 ical->setValue( eedetidpos, sta );
                         }
                         if (EEDetId::validDetId(iX,iY,-1)) {
                                 EEDetId eedetidneg(iX,iY,-1);
                                 ical->setValue( eedetidneg, sta );
                         }
                 }
         }
         return ical;
 }

std::auto_ptr< EcalDQMTowerStatus > EcalTrivialConditionRetriever::produceEcalDQMTowerStatus ( const EcalDQMTowerStatusRcd & )

virtual

Definition at line 2058 of file EcalTrivialConditionRetriever.cc.

References EcalBarrel, i, j, gen::k, EcalTrigTowerDetId::validDetId(), and EcalScDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
 
         std::auto_ptr<EcalDQMTowerStatus>  ical = std::auto_ptr<EcalDQMTowerStatus>( new EcalDQMTowerStatus() );
 
     uint32_t sta(0);
     
         // barrel
     int iz=0;
         for(int k=0 ; k<2; k++ ) {
       if(k==0) iz=-1;
       if(k==1) iz=+1;
       for(int i=1 ; i<73; i++) {
         for(int j=1 ; j<18; j++) {
           if (EcalTrigTowerDetId::validDetId(iz,EcalBarrel,j,i )){
         EcalTrigTowerDetId ebid(iz,EcalBarrel,j,i);
 
         ical->setValue( ebid, sta );
           }
         }
       }
     }
 
 
         // endcap
         for(int k=0 ; k<2; k++ ) {
       if(k==0) iz=-1;
       if(k==1) iz=+1;
       for(int i=1 ; i<21; i++) {
         for(int j=1 ; j<21; j++) {
           if (EcalScDetId::validDetId(i,j,iz )){
         EcalScDetId eeid(i,j,iz);
         ical->setValue( eeid, sta );
           }
         }
       }
     }
 
         return ical;
 }

std::auto_ptr< EcalGainRatios > EcalTrivialConditionRetriever::produceEcalGainRatios ( const EcalGainRatiosRcd & )

virtual

Definition at line 840 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, DetId::rawId(), EcalMGPAGainRatio::setGain12Over6(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalGainRatios> gratio = std::auto_ptr<EcalGainRatios>( new EcalGainRatios() );
   EcalMGPAGainRatio gr;
   gr.setGain12Over6( gainRatio12over6_ );
   gr.setGain6Over1( gainRatio6over1_ );
 
   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) {
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       gratio->setValue( ebid.rawId(), gr );
     }
     }
   }
   
   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))
     {
       EEDetId eedetidpos(iX,iY,1);
       gratio->setValue( eedetidpos.rawId(), gr );
     }
       if (EEDetId::validDetId(iX,iY,-1))
     {
       EEDetId eedetidneg(iX,iY,-1);
       gratio->setValue( eedetidneg.rawId(), gr );
     }
     }
   }
   
   return gratio;
 }

std::auto_ptr< EcalIntercalibConstants > EcalTrivialConditionRetriever::produceEcalIntercalibConstants ( const EcalIntercalibConstantsRcd & )

virtual

Definition at line 643 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, alignCSCRings::r, diffTwoXMLs::r1, rand(), DetId::rawId(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalIntercalibConstants>  ical = std::auto_ptr<EcalIntercalibConstants>( 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) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       ical->setValue( ebid.rawId(), intercalibConstantMean_ + r*intercalibConstantSigma_ );
     }
     }
   }
 
   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))
     {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos.rawId(), intercalibConstantMean_ + r*intercalibConstantSigma_ );
     }
       if(EEDetId::validDetId(iX,iY,-1))
         {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg.rawId(), intercalibConstantMean_ + r1*intercalibConstantSigma_ );
     }
     }
   }
   
   return ical;
 }

std::auto_ptr< EcalIntercalibConstantsMC > EcalTrivialConditionRetriever::produceEcalIntercalibConstantsMC ( const EcalIntercalibConstantsMCRcd & )

virtual

Definition at line 682 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, alignCSCRings::r, diffTwoXMLs::r1, rand(), DetId::rawId(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalIntercalibConstantsMC>  ical = std::auto_ptr<EcalIntercalibConstantsMC>( new EcalIntercalibConstantsMC() );
 
   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) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       ical->setValue( ebid.rawId(), intercalibConstantMeanMC_ + r*intercalibConstantSigmaMC_ );
     }
     }
   }
 
   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))
     {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos.rawId(), intercalibConstantMeanMC_ + r*intercalibConstantSigmaMC_ );
     }
       if(EEDetId::validDetId(iX,iY,-1))
         {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg.rawId(), intercalibConstantMeanMC_ + r1*intercalibConstantSigmaMC_ );
     }
     }
   }
   
   return ical;
 }

std::auto_ptr< EcalIntercalibErrors > EcalTrivialConditionRetriever::produceEcalIntercalibErrors ( const EcalIntercalibErrorsRcd & )

virtual

Definition at line 721 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, DetId::rawId(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalIntercalibErrors>  ical = std::auto_ptr<EcalIntercalibErrors>( new EcalIntercalibErrors() );
 
   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) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       ical->setValue( ebid.rawId(), intercalibErrorMean_);
     }
     }
   }
 
   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))
     {
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos.rawId(), intercalibErrorMean_ );
     }
       if(EEDetId::validDetId(iX,iY,-1))
         {
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg.rawId(), intercalibErrorMean_ );
     }
     }
   }
   
   return ical;
 }

std::auto_ptr< EcalLaserAlphas > EcalTrivialConditionRetriever::produceEcalLaserAlphas ( const EcalLaserAlphasRcd & )

virtual

Definition at line 1021 of file EcalTrivialConditionRetriever.cc.

References funct::abs(), alpha, python.rootplot.core::batch, edm::check(), gather_cfg::cout, contentValuesFiles::fullPath, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, alignCSCRings::r, diffTwoXMLs::r1, rand(), EBDetId::SMCRYSTALMODE, AlCaHLTBitMon_QueryRunRegistry::string, EBDetId::validDetId(), EEDetId::validDetId(), x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by EcalTrivialConditionRetriever().

 {
 
   std::cout << " produceEcalLaserAlphas " << std::endl;
   std::auto_ptr<EcalLaserAlphas> ical = std::auto_ptr<EcalLaserAlphas>( new EcalLaserAlphas() );
   if(getLaserAlphaFromFile_) {
     std::ifstream fEB(edm::FileInPath(EBLaserAlphaFile_).fullPath().c_str());
     int SMpos[36] = {-10, 4, -7, -16, 6, -9, 11, -17, 5, 18, 3, -8, 1, -3, -13, 14, -6, 2,
              15, -18, 8, 17, -2, 9, -1, 10, -5, 7, -12, -11, 16, -4, -15, -14, 12, 13};
     // check!
     int SMCal[36] = {12,17,10, 1, 8, 4,27,20,23,25, 6,34,35,15,18,30,21, 9,
              24,22,13,31,26,16, 2,11, 5, 0,29,28,14,33,32, 3, 7,19};
     /*
   int slot_to_constr[37]={-1,12,17,10,1,8,4,27,20,23,25,6,34,35,15,18,30,21,9
               ,24,22,13,31,26,16,2,11,5,0,29,28,14,33,32,3,7,19};
   int constr_to_slot[36]={28,4,25,34,6,27,11,35,5,18,3,26,1,21,31,14,24,2,15,
               36,8,17,20,9,19,10,23,7,30,29,16,22,33,32,12,13  };
     */
     for(int SMcons = 0; SMcons < 36; SMcons++) {
       int SM = SMpos[SMcons];
       if(SM < 0) SM = 17 + abs(SM);
       else SM--;
       if(SMCal[SM] != SMcons)
      std::cout << " SM pb : read SM " <<  SMcons<< " SMpos " << SM
            << " SMCal " << SMCal[SM] << std::endl;
     }
     // check
     std::string type, batch;
     int readSM, pos, bar, bar2;
     float alpha = 0;
     for(int SMcons = 0; SMcons < 36; SMcons++) {
       int SM = SMpos[SMcons];
       for(int ic = 0; ic < 1700; ic++) {
     fEB >> readSM >> pos >> bar >>  bar2 >> type >> batch;
     //  if(ic == 0) std::cout << readSM << " " << pos << " " << bar << " " << bar2 << " " 
     //                << type << " " << batch << std::endl;
     if(readSM != SMcons || pos != ic + 1) 
       std::cout << " barrel read pb read SM " << readSM << " const SM " << SMcons
             << " read pos " << pos << " ic " << ic << std::endl;
     if(SM < 0) SM = 18 + abs(SM);
     EBDetId ebdetid(SM, pos, EBDetId::SMCRYSTALMODE);
     if(bar == 33101 || bar == 30301 )
       alpha = 1.52;
     else if(bar == 33106) {
       if(bar2 <= 2000)
         alpha = 1.0;
       else {
         std::cout << " problem with barcode first " << bar << " last " << bar2 
               << " read SM " << readSM << " read pos " << pos << std::endl;
         alpha = 0.0;
       }
     }
     ical->setValue( ebdetid, alpha );
       }
     }  // loop over SMcons
   }   // laserAlpha from a file
   else {
     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) {
     if (EBDetId::validDetId(ieta,iphi)) {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       ical->setValue( ebid, laserAlphaMean_ + r*laserAlphaSigma_ );
     }
       }  // loop over iphi
     }  // loop over ieta
   }   // do not read a file
 
   std::cout << " produceEcalLaserAlphas EE" << std::endl;
   if(getLaserAlphaFromFile_) {
     std::ifstream fEE(edm::FileInPath(EELaserAlphaFile_).fullPath().c_str());
     int check[101][101];
     for(int x = 1; x < 101; x++)
       for(int y = 1; y < 101; y++)
     check[x][y] = -1;
     for(int crystal = 0; crystal < 14648; crystal++) {
       int x, y ,z, bid, bar, bar2;
       float LY, alpha = 0;
       fEE >> z >> x >> y >> LY >> bid >> bar >> bar2;
       if(x < 1 || x > 100 || y < 1 || y > 100)
     std::cout << " wrong coordinates for barcode " << bar 
           << " x " << x << " y " << y << " z " << z << std::endl;
       else {
     if(z == 1) check[x][y] = 1;
     else check[x][y] = 0;
     if(bar == 33201 || (bar == 30399 && bar2 < 568))
       alpha = 1.52;
     else if((bar == 33106 && bar2 > 2000 && bar2 < 4669) 
         || (bar == 30399 && bar2 > 567))
       alpha = 1.0;
     else {
       std::cout << " problem with barcode " << bar << " " << bar2 
             << " x " << x << " y " << y << " z " << z << std::endl;
       alpha = 0.0;
     }
       } 
       if (EEDetId::validDetId(x, y, z)) {
     EEDetId eedetidpos(x, y, z);
     ical->setValue( eedetidpos, alpha );
       }
       else // should not occur
     std::cout << " problem with EEDetId " << " x " << x << " y " << y << " z " << z << std::endl;
     }  // loop over crystal in file
     for(int x = 1; x < 101; x++)
       for(int y = 1; y < 101; y++)
     if(check[x][y] == 1) std::cout << " missing x " << x << " y " << y << std::endl;
   }  // laserAlpha from a file
   else {
     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)) {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos, laserAlphaMean_ + r*laserAlphaSigma_ );
     }
 
     if (EEDetId::validDetId(iX,iY,-1)) {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg, laserAlphaMean_ + r1*laserAlphaSigma_ );
     }
       } // loop over iY
     } // loop over iX
   }  // do not read a file
   
   return ical;
 }

std::auto_ptr< EcalLaserAPDPNRatios > EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatios ( const EcalLaserAPDPNRatiosRcd & )

virtual

Definition at line 1189 of file EcalTrivialConditionRetriever.cc.

References EBDetId::approxEta(), EnergyResolutionVsLumi::calcampDropTotal(), gather_cfg::cout, drop, edm::Timestamp::endOfTime(), eta(), i, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, fff_deleter::log, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p1, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p2, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p3, alignCSCRings::r, diffTwoXMLs::r1, rand(), EnergyResolutionVsLumi::setInstLumi(), EnergyResolutionVsLumi::setLumi(), mathSSE::sqrt(), EcalLaserAPDPNRatios::EcalLaserTimeStamp::t1, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t2, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t3, funct::tan(), EBDetId::validDetId(), EEDetId::validDetId(), and relativeConstraints::value.

Referenced by EcalTrivialConditionRetriever().

 {
   EnergyResolutionVsLumi ageing;
   ageing.setLumi(totLumi_);
   ageing.setInstLumi(instLumi_);
 
 
   std::auto_ptr<EcalLaserAPDPNRatios>  ical = std::auto_ptr<EcalLaserAPDPNRatios>( new EcalLaserAPDPNRatios() );
   for(int ieta=-EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA; ++ieta) {
     if(ieta==0) continue;
 
 
     
     double eta=EBDetId::approxEta(EBDetId(ieta,1));
     
     eta = fabs(eta);
     double drop=ageing.calcampDropTotal(eta);
     std::cout<<"EB at eta="<<eta<<" dropping by "<<drop<<std::endl;
     
     for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId::MAX_IPHI; ++iphi) {
       if (EBDetId::validDetId(ieta,iphi)) {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
 
       EcalLaserAPDPNRatios::EcalLaserAPDPNpair pairAPDPN;
       pairAPDPN.p1 = laserAPDPNMean_*drop + r*laserAPDPNSigma_;
       pairAPDPN.p2 = laserAPDPNMean_*drop + r*laserAPDPNSigma_;
       pairAPDPN.p3 = laserAPDPNMean_*drop + r*laserAPDPNSigma_;
       ical->setValue( ebid, pairAPDPN );
       }
      }
   }
 
 
   std::cout<<"----- EE -----"<<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)) {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
 
       double eta= -log(tan(0.5*atan(sqrt((iX-50.0)
                          *(iX-50.0)+
                          (iY-50.0)*
                          (iY-50.0))
                     *2.98/328.)));
       eta = fabs(eta);
       double drop=ageing.calcampDropTotal(eta);
       if(iX==50) std::cout<<"EE at eta="<<eta<<" dropping by "<<drop<<std::endl;
       
 
       EcalLaserAPDPNRatios::EcalLaserAPDPNpair pairAPDPN;
       pairAPDPN.p1 = laserAPDPNMean_*drop + r*laserAPDPNSigma_;
       pairAPDPN.p2 = laserAPDPNMean_*drop + r*laserAPDPNSigma_;
       pairAPDPN.p3 = laserAPDPNMean_*drop + r*laserAPDPNSigma_;
       ical->setValue( eedetidpos, pairAPDPN );
        }
 
        if (EEDetId::validDetId(iX,iY,-1)) {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
 
       double eta= -log(tan(0.5*atan(sqrt((iX-50.0)*(iX-50.0)+(iY-50.0)*(iY-50.0))*2.98/328.)));
       eta = fabs(eta);
       double drop=ageing.calcampDropTotal(eta);
       if(iX==50) std::cout<<"EE at eta="<<eta<<" dropping by "<<drop<<std::endl;
 
 
       EcalLaserAPDPNRatios::EcalLaserAPDPNpair pairAPDPN;
       pairAPDPN.p1 = laserAPDPNMean_*drop + r1*laserAPDPNSigma_;
       pairAPDPN.p2 = laserAPDPNMean_*drop + r1*laserAPDPNSigma_;
       pairAPDPN.p3 = laserAPDPNMean_*drop + r1*laserAPDPNSigma_;
       ical->setValue( eedetidneg, pairAPDPN );
        }
      }
    }
 
   EcalLaserAPDPNRatios::EcalLaserTimeStamp TimeStamp;
   //  for(int i=1; i<=92; i++){
   for(int i=0; i<92; i++){
     TimeStamp.t1 = Timestamp(laserAPDPNTime1_);
     if(laserAPDPNTime2_ == 0 ){ 
       TimeStamp.t2 = Timestamp(edm::Timestamp::endOfTime().value());
     } else {
       TimeStamp.t2 = Timestamp(laserAPDPNTime2_);
     }
     if(laserAPDPNTime3_ == 0 ){ 
       TimeStamp.t3 = Timestamp(edm::Timestamp::endOfTime().value());
     } else {
       TimeStamp.t3 = Timestamp(laserAPDPNTime3_);
     }
 
     ical->setTime( i, TimeStamp );
   }
   
   return ical;
 
 }

std::auto_ptr< EcalLaserAPDPNRatiosRef > EcalTrivialConditionRetriever::produceEcalLaserAPDPNRatiosRef ( const EcalLaserAPDPNRatiosRefRcd & )

virtual

Definition at line 1153 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, alignCSCRings::r, diffTwoXMLs::r1, rand(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalLaserAPDPNRatiosRef>  ical = std::auto_ptr<EcalLaserAPDPNRatiosRef>( new EcalLaserAPDPNRatiosRef() );
   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) {
       if (EBDetId::validDetId(ieta,iphi)) {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       ical->setValue( ebid, laserAPDPNRefMean_ + r*laserAPDPNRefSigma_ );
       }
     }
   }
 
    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)) {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos, laserAPDPNRefMean_ + r*laserAPDPNRefSigma_ );
        }
 
        if (EEDetId::validDetId(iX,iY,-1)) {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg, laserAPDPNRefMean_ + r1*laserAPDPNRefSigma_ );
        }
      }
    }
   
   return ical;
 }

std::auto_ptr< EcalLinearCorrections > EcalTrivialConditionRetriever::produceEcalLinearCorrections ( const EcalLinearCorrectionsRcd & )

virtual

Definition at line 569 of file EcalTrivialConditionRetriever.cc.

References edm::Timestamp::endOfTime(), i, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalTimeDependentCorrections::Values::p1, EcalTimeDependentCorrections::Values::p2, EcalTimeDependentCorrections::Values::p3, alignCSCRings::r, diffTwoXMLs::r1, rand(), EcalTimeDependentCorrections::Times::t1, EcalTimeDependentCorrections::Times::t2, EcalTimeDependentCorrections::Times::t3, EBDetId::validDetId(), EEDetId::validDetId(), and relativeConstraints::value.

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalLinearCorrections>  ical = std::auto_ptr<EcalLinearCorrections>( new EcalLinearCorrections() );
 
   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) {
       if (EBDetId::validDetId(ieta,iphi)) {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
 
       EcalLinearCorrections::Values pairAPDPN;
       pairAPDPN.p1 = linCorrMean_ + r*linCorrSigma_;
       pairAPDPN.p2 = linCorrMean_ + r*linCorrSigma_;
       pairAPDPN.p3 = linCorrMean_ + r*linCorrSigma_;
       ical->setValue( ebid, pairAPDPN );
       }
      }
   }
 
    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)) {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
 
       EcalLinearCorrections::Values pairAPDPN;
       pairAPDPN.p1 = linCorrMean_ + r*linCorrSigma_;
       pairAPDPN.p2 = linCorrMean_ + r*linCorrSigma_;
       pairAPDPN.p3 = linCorrMean_ + r*linCorrSigma_;
 
       ical->setValue( eedetidpos, pairAPDPN );
        }
 
        if (EEDetId::validDetId(iX,iY,-1)) {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
 
       EcalLinearCorrections::Values pairAPDPN;
       pairAPDPN.p1 = linCorrMean_ + r1*linCorrSigma_;
       pairAPDPN.p2 = linCorrMean_ + r1*linCorrSigma_;
       pairAPDPN.p3 = linCorrMean_ + r1*linCorrSigma_;
 
       ical->setValue( eedetidneg, pairAPDPN );
        }
      }
    }
 
   EcalLinearCorrections::Times TimeStamp;
   //  for(int i=1; i<=92; i++){
   for(int i=0; i<92; i++){
     TimeStamp.t1 = Timestamp(linearTime1_);
     if(linearTime2_ == 0 ){ 
       TimeStamp.t2 = Timestamp(edm::Timestamp::endOfTime().value());
     } else {
       TimeStamp.t2 = Timestamp(linearTime2_);
     }
     if(linearTime3_ == 0 ){ 
       TimeStamp.t3 = Timestamp(edm::Timestamp::endOfTime().value());
     } else {
       TimeStamp.t3 = Timestamp(linearTime3_);
     }
 
     ical->setTime( i, TimeStamp );
   }
   
   return ical;
 
 }

std::auto_ptr< EcalMappingElectronics > EcalTrivialConditionRetriever::produceEcalMappingElectronics ( const EcalMappingElectronicsRcd & )

virtual

Definition at line 3001 of file EcalTrivialConditionRetriever.cc.

Referenced by EcalTrivialConditionRetriever().

 {
 
         std::auto_ptr<EcalMappingElectronics>  ical = std::auto_ptr<EcalMappingElectronics>( new EcalMappingElectronics() );
         return ical;
 }

std::auto_ptr< EcalPedestals > EcalTrivialConditionRetriever::produceEcalPedestals ( const EcalPedestalsRcd & )

virtual

Definition at line 455 of file EcalTrivialConditionRetriever.cc.

References EBDetId::approxEta(), EnergyResolutionVsLumi::calcnoiseIncreaseADC(), gather_cfg::cout, eta(), EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, EBDetId::MIN_IPHI, DetId::rawId(), EcalPedestal::rms_x1, EcalPedestal::rms_x12, EcalPedestal::rms_x6, EnergyResolutionVsLumi::setInstLumi(), EnergyResolutionVsLumi::setLumi(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

                                                                              {
   std::auto_ptr<EcalPedestals>  peds = std::auto_ptr<EcalPedestals>( new EcalPedestals() );
   EcalPedestals::Item EBitem;
   EcalPedestals::Item EEitem;
   
   EBitem.mean_x1  = EBpedMeanX1_;
   EBitem.rms_x1   = EBpedRMSX1_;
   EBitem.mean_x6  = EBpedMeanX6_;
   EBitem.rms_x6   = EBpedRMSX6_;
   EBitem.mean_x12 = EBpedMeanX12_;
   EBitem.rms_x12  = EBpedRMSX12_;
 
   EEitem.mean_x1  = EEpedMeanX1_;
   EEitem.rms_x1   = EEpedRMSX1_;
   EEitem.mean_x6  = EEpedMeanX6_;
   EEitem.rms_x6   = EEpedRMSX6_;
   EEitem.mean_x12 = EEpedMeanX12_;
   EEitem.rms_x12  = EEpedRMSX12_;
   
 
 
 
 
   for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
     if(iEta==0) continue;
 
     if(totLumi_>0) {
       double eta=EBDetId::approxEta(EBDetId(iEta,1));
 
       EnergyResolutionVsLumi ageing; 
       ageing.setLumi(totLumi_);
       ageing.setInstLumi(instLumi_);
       eta = fabs(eta);
       double noisefactor= ageing.calcnoiseIncreaseADC(eta);
 
 
       EBitem.rms_x1   = EBpedRMSX1_*noisefactor;
       EBitem.rms_x6   = EBpedRMSX6_*noisefactor;
       EBitem.rms_x12  = EBpedRMSX12_*noisefactor;
       std::cout << "rms ped at eta:"<< eta<<" ="<< EBitem.rms_x12 << std::endl;
     }
 
 
 
     for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(iEta,iPhi))
     {
       EBDetId ebdetid(iEta,iPhi);
       peds->insert(std::make_pair(ebdetid.rawId(),EBitem));
     }
     }
   }
   
   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))
     {
       EEDetId eedetidpos(iX,iY,1);
       peds->insert(std::make_pair(eedetidpos.rawId(),EEitem));
     }
       if(EEDetId::validDetId(iX,iY,-1))
     {
       EEDetId eedetidneg(iX,iY,-1);
       peds->insert(std::make_pair(eedetidneg.rawId(),EEitem));
     }
     }
   }
 
   //return std::auto_ptr<EcalPedestals>( peds );
   return peds;
 }

std::auto_ptr< EcalSampleMask > EcalTrivialConditionRetriever::produceEcalSampleMask ( const EcalSampleMaskRcd & )

virtual

Definition at line 3138 of file EcalTrivialConditionRetriever.cc.

Referenced by EcalTrivialConditionRetriever().

 {
   return std::auto_ptr<EcalSampleMask>( new EcalSampleMask(sampleMaskEB_, sampleMaskEE_) );
 }

std::auto_ptr< EcalTBWeights > EcalTrivialConditionRetriever::produceEcalTBWeights ( const EcalTBWeightsRcd & )

virtual

   for(size_t i=0; i<3; ++i) {
   std::vector<EcalWeight> tv1, tv2;
   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; tv1.push_back( EcalWeight(ww) ); tv2.push_back( EcalWeight(100+ww) ); } mat1.push_back(tv1); mat2.push_back(tv2); }

Definition at line 884 of file EcalTrivialConditionRetriever.cc.

References EcalWeightSet::getChi2WeightsAfterGainSwitch(), EcalWeightSet::getChi2WeightsBeforeGainSwitch(), EcalWeightSet::getWeightsAfterGainSwitch(), and EcalWeightSet::getWeightsBeforeGainSwitch().

Referenced by EcalTrivialConditionRetriever().

 {
   // create weights for the test-beam
   std::auto_ptr<EcalTBWeights> tbwgt = std::auto_ptr<EcalTBWeights>( 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; 
   int igrp=1;
   for(int itdc=1; itdc<=nTDCbins_; ++itdc) {
     // generate random number
     //    double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
     
     // make a new set of weights
     EcalWeightSet wgt;
     //typedef std::vector< std::vector<EcalWeight> > EcalWeightSet::EcalWeightMatrix;
     EcalWeightSet::EcalWeightMatrix& mat1 = wgt.getWeightsBeforeGainSwitch();
     EcalWeightSet::EcalWeightMatrix& mat2 = wgt.getWeightsAfterGainSwitch();
     
 //     if(verbose_>=1) {
 //       std::cout << "initial size of mat1: " << mat1.size() << std::endl;
 //       std::cout << "initial size of mat2: " << mat2.size() << std::endl;
 //     }
     
     // generate random numbers to use as weights
     // use values provided by user
     mat1.Place_in_row(amplWeights_[itdc-1],0,0);
     mat1.Place_in_row(pedWeights_[itdc-1],1,0);
     mat1.Place_in_row(jittWeights_[itdc-1],2,0);
     
     // wdights after gain switch 
     mat2.Place_in_row(amplWeightsAft_[itdc-1],0,0);
     mat2.Place_in_row(pedWeightsAft_[itdc-1],1,0);
     mat2.Place_in_row(jittWeightsAft_[itdc-1],2,0);
     
     // fill the chi2 matrcies with random numbers
     //    r = (double)std::rand()/( double(RAND_MAX)+double(1) );
     EcalWeightSet::EcalChi2WeightMatrix& mat3 = wgt.getChi2WeightsBeforeGainSwitch();
     EcalWeightSet::EcalChi2WeightMatrix& mat4 = wgt.getChi2WeightsAfterGainSwitch();
     mat3=chi2Matrix_[itdc-1];
     mat4=chi2MatrixAft_[itdc-1];
     
     //     for(size_t i=0; i<10; ++i) 
     //       {
     //  mat3.push_back(chi2Matrix_[itdc-1][i]);
     //  mat4.push_back(chi2MatrixAft_[itdc-1][i]);
     //       }
     //       std::vector<EcalWeight> tv1, tv2;
     //       for(size_t j=0; j<10; ++j) {
     //  double ww = igrp*itdc*r + i*10. + j;
     //  tv1.push_back( EcalWeight(1000+ww) );
     //  tv2.push_back( EcalWeight(1000+100+ww) );
     //       }
     
   
     
     
 //     if(verbose_>=1) {
 //       std::cout << "group: " << igrp << " TDC: " << itdc 
 //      << " mat1: " << mat1.size() << " mat2: " << mat2.size()
 //      << " mat3: " << mat3.size() << " mat4: " << mat4.size()
 //      << std::endl;
 //     }
     
     // put the weight in the container
     tbwgt->setValue(std::make_pair(igrp,itdc), wgt);
   } 
   //   }
   return tbwgt;
 }

std::auto_ptr< EcalTimeBiasCorrections > EcalTrivialConditionRetriever::produceEcalTimeBiasCorrections ( const EcalTimeBiasCorrectionsRcd & )

virtual

Definition at line 3144 of file EcalTrivialConditionRetriever.cc.

References filterCSVwithJSON::copy.

Referenced by EcalTrivialConditionRetriever().

                                                                                                  {
   std::auto_ptr<EcalTimeBiasCorrections> ipar = std::auto_ptr<EcalTimeBiasCorrections>( new EcalTimeBiasCorrections() );
   copy(EBtimeCorrAmplitudeBins_.begin(), EBtimeCorrAmplitudeBins_.end(),
        back_inserter(ipar->EBTimeCorrAmplitudeBins));
   copy(EBtimeCorrShiftBins_.begin(), EBtimeCorrShiftBins_.end(),
        back_inserter(ipar->EBTimeCorrShiftBins));
   copy(EEtimeCorrAmplitudeBins_.begin(), EEtimeCorrAmplitudeBins_.end(),
        back_inserter(ipar->EETimeCorrAmplitudeBins));
   copy(EEtimeCorrShiftBins_.begin(), EEtimeCorrShiftBins_.end(),
        back_inserter(ipar->EETimeCorrShiftBins));
   return ipar;
 }

std::auto_ptr< EcalTimeCalibConstants > EcalTrivialConditionRetriever::produceEcalTimeCalibConstants ( const EcalTimeCalibConstantsRcd & )

virtual

Definition at line 757 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, alignCSCRings::r, diffTwoXMLs::r1, rand(), DetId::rawId(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalTimeCalibConstants>  ical = std::auto_ptr<EcalTimeCalibConstants>( new EcalTimeCalibConstants() );
 
   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) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       ical->setValue( ebid.rawId(), timeCalibConstantMean_ + r*timeCalibConstantSigma_ );
     }
     }
   }
 
   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))
     {
       double r = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos.rawId(), timeCalibConstantMean_ + r*timeCalibConstantSigma_ );
     }
       if(EEDetId::validDetId(iX,iY,-1))
         {
       double r1 = (double)std::rand()/( double(RAND_MAX)+double(1) );
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg.rawId(), timeCalibConstantMean_ + r1*timeCalibConstantSigma_ );
     }
     }
   }
   
   return ical;
 }

std::auto_ptr< EcalTimeCalibErrors > EcalTrivialConditionRetriever::produceEcalTimeCalibErrors ( const EcalTimeCalibErrorsRcd & )

virtual

Definition at line 796 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, DetId::rawId(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalTimeCalibErrors>  ical = std::auto_ptr<EcalTimeCalibErrors>( new EcalTimeCalibErrors() );
 
   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) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       ical->setValue( ebid.rawId(), timeCalibErrorMean_);
     }
     }
   }
 
   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))
     {
       EEDetId eedetidpos(iX,iY,1);
       ical->setValue( eedetidpos.rawId(), timeCalibErrorMean_ );
     }
       if(EEDetId::validDetId(iX,iY,-1))
         {
       EEDetId eedetidneg(iX,iY,-1);
       ical->setValue( eedetidneg.rawId(), timeCalibErrorMean_ );
     }
     }
   }
   
   return ical;
 }

std::auto_ptr< EcalTimeOffsetConstant > EcalTrivialConditionRetriever::produceEcalTimeOffsetConstant ( const EcalTimeOffsetConstantRcd & )

virtual

Definition at line 832 of file EcalTrivialConditionRetriever.cc.

References gather_cfg::cout.

Referenced by EcalTrivialConditionRetriever().

 {
   std::cout << " produceEcalTimeOffsetConstant: " << std::endl;
   std::cout << "  EB " << timeOffsetEBConstant_ << " EE " <<  timeOffsetEEConstant_<< std::endl;
   return std::auto_ptr<EcalTimeOffsetConstant>( new EcalTimeOffsetConstant(timeOffsetEBConstant_,timeOffsetEEConstant_) );
 }

std::auto_ptr< EcalTPGCrystalStatus > EcalTrivialConditionRetriever::produceEcalTrgChannelStatus ( const EcalTPGCrystalStatusRcd & )

virtual

Definition at line 2277 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
 
         std::auto_ptr<EcalTPGCrystalStatus>  ical = std::auto_ptr<EcalTPGCrystalStatus>( new EcalTPGCrystalStatus() );
         // barrel
         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) {
                         if (EBDetId::validDetId(ieta,iphi)) {
                                 EBDetId ebid(ieta,iphi);
                                 ical->setValue( ebid, 0 );
                         }
                 }
         }
         // endcap
         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)) {
                                 EEDetId eedetidpos(iX,iY,1);
                                 ical->setValue( eedetidpos, 0 );
                         }
                         if (EEDetId::validDetId(iX,iY,-1)) {
                                 EEDetId eedetidneg(iX,iY,-1);
                                 ical->setValue( eedetidneg, 0 );
                         }
                 }
         }
         return ical;
 }

std::auto_ptr< EcalWeightXtalGroups > EcalTrivialConditionRetriever::produceEcalWeightXtalGroups ( const EcalWeightXtalGroupsRcd & )

virtual

Definition at line 532 of file EcalTrivialConditionRetriever.cc.

References EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, DetId::rawId(), EBDetId::validDetId(), and EEDetId::validDetId().

Referenced by EcalTrivialConditionRetriever().

 {
   std::auto_ptr<EcalWeightXtalGroups> xtalGroups = std::auto_ptr<EcalWeightXtalGroups>( new EcalWeightXtalGroups() );
   EcalXtalGroupId defaultGroupId(1);
   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) {
       // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
       if (EBDetId::validDetId(ieta,iphi))
     {
       EBDetId ebid(ieta,iphi);
       //      xtalGroups->setValue(ebid.rawId(), EcalXtalGroupId(ieta) ); // define rings in eta
       xtalGroups->setValue(ebid.rawId(), defaultGroupId ); // define rings in eta
     }
     }
   }
 
   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))
     {
       EEDetId eedetidpos(iX,iY,1);
       xtalGroups->setValue(eedetidpos.rawId(), defaultGroupId ); 
     }
       if(EEDetId::validDetId(iX,iY,-1))
     {
       EEDetId eedetidneg(iX,iY,-1);
       xtalGroups->setValue(eedetidneg.rawId(), defaultGroupId ); 
     }
     }
   }
   return xtalGroups;
 }

void EcalTrivialConditionRetriever::setIntervalFor	(	const edm::eventsetup::EventSetupRecordKey &	rk,
		const edm::IOVSyncValue &	iTime,
		edm::ValidityInterval &	oValidity
	)

protectedvirtual

Implements edm::EventSetupRecordIntervalFinder.

Definition at line 444 of file EcalTrivialConditionRetriever.cc.

References edm::IOVSyncValue::beginOfTime(), gather_cfg::cout, edm::IOVSyncValue::endOfTime(), edm::eventsetup::EventSetupRecordKey::name(), edm::IOVSyncValue::time(), and edm::Timestamp::value().

 {
   if(verbose_>=1) std::cout << "EcalTrivialConditionRetriever::setIntervalFor(): record key = " << rk.name() << "\ttime: " << iTime.time().value() << std::endl;
   //For right now, we will just use an infinite interval of validity
   oValidity = edm::ValidityInterval( edm::IOVSyncValue::beginOfTime(),edm::IOVSyncValue::endOfTime() );
 }

Member Data Documentation

double EcalTrivialConditionRetriever::adcToGeVEBConstant_

private

Definition at line 190 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::adcToGeVEEConstant_

private

Definition at line 191 of file EcalTrivialConditionRetriever.h.

std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::amplWeights_

private

Definition at line 256 of file EcalTrivialConditionRetriever.h.

std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::amplWeightsAft_

private

Definition at line 257 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::amplWeightsAftFile_

private

Definition at line 269 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::amplWeightsFile_

private

Definition at line 268 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::channelStatusFile_

private

Definition at line 282 of file EcalTrivialConditionRetriever.h.

std::vector< EcalWeightSet::EcalChi2WeightMatrix > EcalTrivialConditionRetriever::chi2Matrix_

private

Definition at line 265 of file EcalTrivialConditionRetriever.h.

std::vector< EcalWeightSet::EcalChi2WeightMatrix > EcalTrivialConditionRetriever::chi2MatrixAft_

private

Definition at line 266 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::chi2MatrixAftFile_

private

Definition at line 275 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::chi2MatrixFile_

private

Definition at line 274 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::crackCorrParameters_

private

Definition at line 212 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::EBAlignmentFile_

private

Definition at line 285 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::EBLaserAlphaFile_

private

Definition at line 288 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EBpedMeanX12_

private

Definition at line 239 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EBpedMeanX1_

private

Definition at line 243 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EBpedMeanX6_

private

Definition at line 241 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EBpedRMSX12_

private

Definition at line 240 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EBpedRMSX1_

private

Definition at line 244 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EBpedRMSX6_

private

Definition at line 242 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::EBtimeCorrAmplitudeBins_

private

Definition at line 292 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::EBtimeCorrShiftBins_

private

Definition at line 293 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::EEAlignmentFile_

private

Definition at line 286 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::EELaserAlphaFile_

private

Definition at line 289 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EEpedMeanX12_

private

Definition at line 246 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EEpedMeanX1_

private

Definition at line 250 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EEpedMeanX6_

private

Definition at line 248 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EEpedRMSX12_

private

Definition at line 247 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EEpedRMSX1_

private

Definition at line 251 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::EEpedRMSX6_

private

Definition at line 249 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::EEtimeCorrAmplitudeBins_

private

Definition at line 294 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::EEtimeCorrShiftBins_

private

Definition at line 295 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::energyCorrectionObjectSpecificParameters_

private

Definition at line 215 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::energyCorrectionParameters_

private

Definition at line 213 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::energyUncertaintyParameters_

private

Definition at line 214 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::ESAlignmentFile_

private

Definition at line 287 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::gainRatio12over6_

private

Definition at line 253 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::gainRatio6over1_

private

Definition at line 254 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::getEBAlignmentFromFile_

private

Definition at line 328 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::getEEAlignmentFromFile_

private

Definition at line 329 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::getESAlignmentFromFile_

private

Definition at line 330 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::getLaserAlphaFromFile_

private

Definition at line 331 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::getWeightsFromFile_

private

Definition at line 299 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::instLumi_

private

Definition at line 219 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::intercalibConstantMean_

private

Definition at line 194 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::intercalibConstantMeanMC_

private

Definition at line 197 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::intercalibConstantsFile_

private

Definition at line 277 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::intercalibConstantSigma_

private

Definition at line 195 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::intercalibConstantSigmaMC_

private

Definition at line 198 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::intercalibConstantsMCFile_

private

Definition at line 278 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::intercalibErrorMean_

private

Definition at line 201 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::intercalibErrorsFile_

private

Definition at line 279 of file EcalTrivialConditionRetriever.h.

std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::jittWeights_

private

Definition at line 262 of file EcalTrivialConditionRetriever.h.

std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::jittWeightsAft_

private

Definition at line 263 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::jittWeightsAftFile_

private

Definition at line 273 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::jittWeightsFile_

private

Definition at line 272 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::laserAlphaMean_

private

Definition at line 222 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::laserAlphaSigma_

private

Definition at line 223 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::laserAPDPNMean_

private

Definition at line 226 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::laserAPDPNRefMean_

private

Definition at line 224 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::laserAPDPNRefSigma_

private

Definition at line 225 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::laserAPDPNSigma_

private

Definition at line 227 of file EcalTrivialConditionRetriever.h.

unsigned long EcalTrivialConditionRetriever::laserAPDPNTime1_

private

Definition at line 228 of file EcalTrivialConditionRetriever.h.

unsigned long EcalTrivialConditionRetriever::laserAPDPNTime2_

private

Definition at line 229 of file EcalTrivialConditionRetriever.h.

unsigned long EcalTrivialConditionRetriever::laserAPDPNTime3_

private

Definition at line 230 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::linCorrMean_

private

Definition at line 232 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::linCorrSigma_

private

Definition at line 233 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::linearCorrectionsFile_

private

Definition at line 276 of file EcalTrivialConditionRetriever.h.

unsigned long EcalTrivialConditionRetriever::linearTime1_

private

Definition at line 235 of file EcalTrivialConditionRetriever.h.

unsigned long EcalTrivialConditionRetriever::linearTime2_

private

Definition at line 236 of file EcalTrivialConditionRetriever.h.

unsigned long EcalTrivialConditionRetriever::linearTime3_

private

Definition at line 237 of file EcalTrivialConditionRetriever.h.

std::vector<double> EcalTrivialConditionRetriever::localContCorrParameters_

private

Definition at line 211 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::mappingFile_

private

Definition at line 284 of file EcalTrivialConditionRetriever.h.

int EcalTrivialConditionRetriever::nTDCbins_

private

Definition at line 297 of file EcalTrivialConditionRetriever.h.

std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::pedWeights_

private

Definition at line 259 of file EcalTrivialConditionRetriever.h.

std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::pedWeightsAft_

private

Definition at line 260 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::pedWeightsAftFile_

private

Definition at line 271 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::pedWeightsFile_

private

Definition at line 270 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalADCToGeVConstant_

private

Definition at line 311 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalAlignmentEB_

private

Definition at line 325 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalAlignmentEE_

private

Definition at line 326 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalAlignmentES_

private

Definition at line 327 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalChannelStatus_

private

Definition at line 313 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalClusterCrackCorrParameters_

private

Definition at line 320 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalClusterEnergyCorrectionObjectSpecificParameters_

private

Definition at line 323 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalClusterEnergyCorrectionParameters_

private

Definition at line 321 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalClusterEnergyUncertaintyParameters_

private

Definition at line 322 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalClusterLocalContCorrParameters_

private

Definition at line 319 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalDAQTowerStatus_

private

Definition at line 317 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalDCSTowerStatus_

private

Definition at line 316 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalDQMChannelStatus_

private

Definition at line 315 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalDQMTowerStatus_

private

Definition at line 314 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalGainRatios_

private

Definition at line 310 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalIntercalibConstants_

private

Definition at line 304 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalIntercalibConstantsMC_

private

Definition at line 305 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalIntercalibErrors_

private

Definition at line 306 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalLaserCorrection_

private

Definition at line 312 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalLinearCorrections_

private

Definition at line 303 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalMappingElectronics_

private

Definition at line 324 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalPedestals_

private

Definition at line 301 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalSampleMask_

private

Definition at line 332 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalTimeBiasCorrections_

private

Definition at line 333 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalTimeCalibConstants_

private

Definition at line 307 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalTimeCalibErrors_

private

Definition at line 308 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalTimeOffsetConstant_

private

Definition at line 309 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalTrgChannelStatus_

private

Definition at line 318 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::producedEcalWeights_

private

Definition at line 302 of file EcalTrivialConditionRetriever.h.

unsigned int EcalTrivialConditionRetriever::sampleMaskEB_

private

Definition at line 290 of file EcalTrivialConditionRetriever.h.

unsigned int EcalTrivialConditionRetriever::sampleMaskEE_

private

Definition at line 291 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::timeCalibConstantMean_

private

Definition at line 203 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::timeCalibConstantsFile_

private

Definition at line 280 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::timeCalibConstantSigma_

private

Definition at line 204 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::timeCalibErrorMean_

private

Definition at line 205 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::timeCalibErrorsFile_

private

Definition at line 281 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::timeOffsetEBConstant_

private

Definition at line 207 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::timeOffsetEEConstant_

private

Definition at line 208 of file EcalTrivialConditionRetriever.h.

double EcalTrivialConditionRetriever::totLumi_

private

Definition at line 218 of file EcalTrivialConditionRetriever.h.

std::string EcalTrivialConditionRetriever::trgChannelStatusFile_

private

Definition at line 283 of file EcalTrivialConditionRetriever.h.

int EcalTrivialConditionRetriever::verbose_

private

Definition at line 335 of file EcalTrivialConditionRetriever.h.

bool EcalTrivialConditionRetriever::weightsForAsynchronousRunning_

private

Definition at line 300 of file EcalTrivialConditionRetriever.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation