|
|
|
#include <EcalTrivialConditionRetriever.h>
Definition at line 114 of file EcalTrivialConditionRetriever.h.
| EcalTrivialConditionRetriever::EcalTrivialConditionRetriever | ( | const edm::ParameterSet & | pset | ) |
Definition at line 31 of file EcalTrivialConditionRetriever.cc.
References gather_cfg::cout, getChannelStatusFromConfiguration(), getIntercalibConstantsFromConfiguration(), getIntercalibConstantsMCFromConfiguration(), getIntercalibErrorsFromConfiguration(), getMappingFromConfiguration(), getTimeCalibConstantsFromConfiguration(), getTimeCalibErrorsFromConfiguration(), getTrgChannelStatusFromConfiguration(), edm::ParameterSet::getUntrackedParameter(), getHLTPrescaleColumns::path, produceEcalADCToGeVConstant(), produceEcalAlignmentEB(), produceEcalAlignmentEE(), produceEcalAlignmentES(), produceEcalChannelStatus(), produceEcalClusterCrackCorrParameters(), produceEcalClusterEnergyCorrectionObjectSpecificParameters(), produceEcalClusterEnergyCorrectionParameters(), produceEcalClusterEnergyUncertaintyParameters(), produceEcalClusterLocalContCorrParameters(), produceEcalDAQTowerStatus(), produceEcalDCSTowerStatus(), produceEcalDQMChannelStatus(), produceEcalDQMTowerStatus(), produceEcalGainRatios(), produceEcalIntercalibConstants(), produceEcalIntercalibConstantsMC(), produceEcalIntercalibErrors(), produceEcalLaserAlphas(), produceEcalLaserAPDPNRatios(), produceEcalLaserAPDPNRatiosRef(), produceEcalLinearCorrections(), produceEcalMappingElectronics(), produceEcalPedestals(), produceEcalSampleMask(), produceEcalTBWeights(), produceEcalTimeCalibConstants(), produceEcalTimeCalibErrors(), produceEcalTimeOffsetConstant(), produceEcalTrgChannelStatus(), produceEcalWeightXtalGroups(), and AlCaHLTBitMon_QueryRunRegistry::string.
{
// 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);
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>();
}
}
| EcalTrivialConditionRetriever::~EcalTrivialConditionRetriever | ( | ) | [virtual] |
Definition at line 427 of file EcalTrivialConditionRetriever.cc.
{
}
| EcalTrivialConditionRetriever::EcalTrivialConditionRetriever | ( | const EcalTrivialConditionRetriever & | ) | [private] |
| std::auto_ptr< EcalChannelStatus > EcalTrivialConditionRetriever::getChannelStatusFromConfiguration | ( | const EcalChannelStatusRcd & | ) | [virtual] |
Definition at line 1892 of file EcalTrivialConditionRetriever.cc.
References Exception, 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(), EEDetId::validDetId(), and EBDetId::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 2304 of file EcalTrivialConditionRetriever.cc.
References EBDetId::approxEta(), EnergyResolutionVsLumi::calcresolutitonConstantTerm(), calib, gather_cfg::cout, EcalCondObjectContainer< T >::end(), eta(), 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, create_public_lumi_plots::log, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MAX_SM, EBDetId::MIN_IPHI, EBDetId::MIN_SM, alignCSCRings::r, DetId::rawId(), insertMaterial::readXML(), EnergyResolutionVsLumi::setInstLumi(), EnergyResolutionVsLumi::setLumi(), EcalCondObjectContainer< T >::setValue(), EBDetId::SMCRYSTALMODE, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), EEDetId::validDetId(), and EBDetId::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);
EcalIntercalibConstants::const_iterator idref=mymap.find(eedetidneg);
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 2563 of file EcalTrivialConditionRetriever.cc.
References gather_cfg::cout, h, and insertMaterial::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 2597 of file EcalTrivialConditionRetriever.cc.
References calib, 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 2946 of file EcalTrivialConditionRetriever.cc.
References cond::rpcobgas::detid, EcalMappingElement::electronicsid, Exception, f, 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 2718 of file EcalTrivialConditionRetriever.cc.
References calib, 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 2832 of file EcalTrivialConditionRetriever.cc.
References calib, 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 2178 of file EcalTrivialConditionRetriever.cc.
References Exception, 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(), EEDetId::validDetId(), and EBDetId::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 1284 of file EcalTrivialConditionRetriever.cc.
References begin, Exception, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), i, j, prof2calltree::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 869 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 3003 of file EcalTrivialConditionRetriever.cc.
References a, gather_cfg::cout, f, Alignments::m_align, and create_public_pileup_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 3060 of file EcalTrivialConditionRetriever.cc.
References a, gather_cfg::cout, f, Alignments::m_align, and create_public_pileup_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 3095 of file EcalTrivialConditionRetriever.cc.
References a, gather_cfg::cout, f, Alignments::m_align, and create_public_pileup_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 1982 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, EEDetId::validDetId(), and EBDetId::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 973 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 1000 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 982 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 991 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 964 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 2135 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 2092 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 2015 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, EEDetId::validDetId(), and EBDetId::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 2049 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 831 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(), EEDetId::validDetId(), and EBDetId::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 634 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, DetId::rawId(), EEDetId::validDetId(), and EBDetId::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 673 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, DetId::rawId(), EEDetId::validDetId(), and EBDetId::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 712 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(), EEDetId::validDetId(), and EBDetId::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 1012 of file EcalTrivialConditionRetriever.cc.
References abs, alpha, python::rootplot::core::batch, edm::check(), gather_cfg::cout, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, pos, alignCSCRings::r, diffTwoXMLs::r1, EBDetId::SMCRYSTALMODE, AlCaHLTBitMon_QueryRunRegistry::string, EEDetId::validDetId(), EBDetId::validDetId(), x, detailsBasic3DVector::y, and 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 1180 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, create_public_lumi_plots::log, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p1, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p2, EcalLaserAPDPNRatios::EcalLaserAPDPNpair::p3, alignCSCRings::r, diffTwoXMLs::r1, EnergyResolutionVsLumi::setInstLumi(), EnergyResolutionVsLumi::setLumi(), mathSSE::sqrt(), EcalLaserAPDPNRatios::EcalLaserTimeStamp::t1, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t2, EcalLaserAPDPNRatios::EcalLaserTimeStamp::t3, funct::tan(), EEDetId::validDetId(), EBDetId::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 1144 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, EEDetId::validDetId(), and EBDetId::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 560 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, EcalTimeDependentCorrections::Times::t1, EcalTimeDependentCorrections::Times::t2, EcalTimeDependentCorrections::Times::t3, EEDetId::validDetId(), EBDetId::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 2993 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 446 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, EBDetId::MIN_IPHI, DetId::rawId(), EnergyResolutionVsLumi::setInstLumi(), EnergyResolutionVsLumi::setLumi(), EEDetId::validDetId(), and EBDetId::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 3130 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 875 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< EcalTimeCalibConstants > EcalTrivialConditionRetriever::produceEcalTimeCalibConstants | ( | const EcalTimeCalibConstantsRcd & | ) | [virtual] |
Definition at line 748 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, DetId::rawId(), EEDetId::validDetId(), and EBDetId::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 787 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(), EEDetId::validDetId(), and EBDetId::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 823 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 2268 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, EEDetId::validDetId(), and EBDetId::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 523 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(), EEDetId::validDetId(), and EBDetId::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 | ||
| ) | [protected, virtual] |
Implements edm::EventSetupRecordIntervalFinder.
Definition at line 435 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() );
}
double EcalTrivialConditionRetriever::adcToGeVEBConstant_ [private] |
Definition at line 187 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::adcToGeVEEConstant_ [private] |
Definition at line 188 of file EcalTrivialConditionRetriever.h.
std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::amplWeights_ [private] |
Definition at line 253 of file EcalTrivialConditionRetriever.h.
std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::amplWeightsAft_ [private] |
Definition at line 254 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::amplWeightsAftFile_ [private] |
Definition at line 266 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::amplWeightsFile_ [private] |
Definition at line 265 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::channelStatusFile_ [private] |
Definition at line 279 of file EcalTrivialConditionRetriever.h.
std::vector< EcalWeightSet::EcalChi2WeightMatrix > EcalTrivialConditionRetriever::chi2Matrix_ [private] |
Definition at line 262 of file EcalTrivialConditionRetriever.h.
std::vector< EcalWeightSet::EcalChi2WeightMatrix > EcalTrivialConditionRetriever::chi2MatrixAft_ [private] |
Definition at line 263 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::chi2MatrixAftFile_ [private] |
Definition at line 272 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::chi2MatrixFile_ [private] |
Definition at line 271 of file EcalTrivialConditionRetriever.h.
std::vector<double> EcalTrivialConditionRetriever::crackCorrParameters_ [private] |
Definition at line 209 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::EBAlignmentFile_ [private] |
Definition at line 282 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::EBLaserAlphaFile_ [private] |
Definition at line 285 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EBpedMeanX12_ [private] |
Definition at line 236 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EBpedMeanX1_ [private] |
Definition at line 240 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EBpedMeanX6_ [private] |
Definition at line 238 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EBpedRMSX12_ [private] |
Definition at line 237 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EBpedRMSX1_ [private] |
Definition at line 241 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EBpedRMSX6_ [private] |
Definition at line 239 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::EEAlignmentFile_ [private] |
Definition at line 283 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::EELaserAlphaFile_ [private] |
Definition at line 286 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EEpedMeanX12_ [private] |
Definition at line 243 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EEpedMeanX1_ [private] |
Definition at line 247 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EEpedMeanX6_ [private] |
Definition at line 245 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EEpedRMSX12_ [private] |
Definition at line 244 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EEpedRMSX1_ [private] |
Definition at line 248 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::EEpedRMSX6_ [private] |
Definition at line 246 of file EcalTrivialConditionRetriever.h.
std::vector<double> EcalTrivialConditionRetriever::energyCorrectionObjectSpecificParameters_ [private] |
Definition at line 212 of file EcalTrivialConditionRetriever.h.
std::vector<double> EcalTrivialConditionRetriever::energyCorrectionParameters_ [private] |
Definition at line 210 of file EcalTrivialConditionRetriever.h.
std::vector<double> EcalTrivialConditionRetriever::energyUncertaintyParameters_ [private] |
Definition at line 211 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::ESAlignmentFile_ [private] |
Definition at line 284 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::gainRatio12over6_ [private] |
Definition at line 250 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::gainRatio6over1_ [private] |
Definition at line 251 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::getEBAlignmentFromFile_ [private] |
Definition at line 321 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::getEEAlignmentFromFile_ [private] |
Definition at line 322 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::getESAlignmentFromFile_ [private] |
Definition at line 323 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::getLaserAlphaFromFile_ [private] |
Definition at line 324 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::getWeightsFromFile_ [private] |
Definition at line 292 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::instLumi_ [private] |
Definition at line 216 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::intercalibConstantMean_ [private] |
Definition at line 191 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::intercalibConstantMeanMC_ [private] |
Definition at line 194 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::intercalibConstantsFile_ [private] |
Definition at line 274 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::intercalibConstantSigma_ [private] |
Definition at line 192 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::intercalibConstantSigmaMC_ [private] |
Definition at line 195 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::intercalibConstantsMCFile_ [private] |
Definition at line 275 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::intercalibErrorMean_ [private] |
Definition at line 198 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::intercalibErrorsFile_ [private] |
Definition at line 276 of file EcalTrivialConditionRetriever.h.
std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::jittWeights_ [private] |
Definition at line 259 of file EcalTrivialConditionRetriever.h.
std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::jittWeightsAft_ [private] |
Definition at line 260 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::jittWeightsAftFile_ [private] |
Definition at line 270 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::jittWeightsFile_ [private] |
Definition at line 269 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::laserAlphaMean_ [private] |
Definition at line 219 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::laserAlphaSigma_ [private] |
Definition at line 220 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::laserAPDPNMean_ [private] |
Definition at line 223 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::laserAPDPNRefMean_ [private] |
Definition at line 221 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::laserAPDPNRefSigma_ [private] |
Definition at line 222 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::laserAPDPNSigma_ [private] |
Definition at line 224 of file EcalTrivialConditionRetriever.h.
unsigned long EcalTrivialConditionRetriever::laserAPDPNTime1_ [private] |
Definition at line 225 of file EcalTrivialConditionRetriever.h.
unsigned long EcalTrivialConditionRetriever::laserAPDPNTime2_ [private] |
Definition at line 226 of file EcalTrivialConditionRetriever.h.
unsigned long EcalTrivialConditionRetriever::laserAPDPNTime3_ [private] |
Definition at line 227 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::linCorrMean_ [private] |
Definition at line 229 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::linCorrSigma_ [private] |
Definition at line 230 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::linearCorrectionsFile_ [private] |
Definition at line 273 of file EcalTrivialConditionRetriever.h.
unsigned long EcalTrivialConditionRetriever::linearTime1_ [private] |
Definition at line 232 of file EcalTrivialConditionRetriever.h.
unsigned long EcalTrivialConditionRetriever::linearTime2_ [private] |
Definition at line 233 of file EcalTrivialConditionRetriever.h.
unsigned long EcalTrivialConditionRetriever::linearTime3_ [private] |
Definition at line 234 of file EcalTrivialConditionRetriever.h.
std::vector<double> EcalTrivialConditionRetriever::localContCorrParameters_ [private] |
Definition at line 208 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::mappingFile_ [private] |
Definition at line 281 of file EcalTrivialConditionRetriever.h.
int EcalTrivialConditionRetriever::nTDCbins_ [private] |
Definition at line 290 of file EcalTrivialConditionRetriever.h.
std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::pedWeights_ [private] |
Definition at line 256 of file EcalTrivialConditionRetriever.h.
std::vector< ROOT::Math::SVector<double,EcalDataFrame::MAXSAMPLES> > EcalTrivialConditionRetriever::pedWeightsAft_ [private] |
Definition at line 257 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::pedWeightsAftFile_ [private] |
Definition at line 268 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::pedWeightsFile_ [private] |
Definition at line 267 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalADCToGeVConstant_ [private] |
Definition at line 304 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalAlignmentEB_ [private] |
Definition at line 318 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalAlignmentEE_ [private] |
Definition at line 319 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalAlignmentES_ [private] |
Definition at line 320 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalChannelStatus_ [private] |
Definition at line 306 of file EcalTrivialConditionRetriever.h.
Definition at line 313 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalClusterEnergyCorrectionObjectSpecificParameters_ [private] |
Definition at line 316 of file EcalTrivialConditionRetriever.h.
Definition at line 314 of file EcalTrivialConditionRetriever.h.
Definition at line 315 of file EcalTrivialConditionRetriever.h.
Definition at line 312 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalDAQTowerStatus_ [private] |
Definition at line 310 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalDCSTowerStatus_ [private] |
Definition at line 309 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalDQMChannelStatus_ [private] |
Definition at line 308 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalDQMTowerStatus_ [private] |
Definition at line 307 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalGainRatios_ [private] |
Definition at line 303 of file EcalTrivialConditionRetriever.h.
Definition at line 297 of file EcalTrivialConditionRetriever.h.
Definition at line 298 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalIntercalibErrors_ [private] |
Definition at line 299 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalLaserCorrection_ [private] |
Definition at line 305 of file EcalTrivialConditionRetriever.h.
Definition at line 296 of file EcalTrivialConditionRetriever.h.
Definition at line 317 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalPedestals_ [private] |
Definition at line 294 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalSampleMask_ [private] |
Definition at line 325 of file EcalTrivialConditionRetriever.h.
Definition at line 300 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalTimeCalibErrors_ [private] |
Definition at line 301 of file EcalTrivialConditionRetriever.h.
Definition at line 302 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalTrgChannelStatus_ [private] |
Definition at line 311 of file EcalTrivialConditionRetriever.h.
bool EcalTrivialConditionRetriever::producedEcalWeights_ [private] |
Definition at line 295 of file EcalTrivialConditionRetriever.h.
unsigned int EcalTrivialConditionRetriever::sampleMaskEB_ [private] |
Definition at line 287 of file EcalTrivialConditionRetriever.h.
unsigned int EcalTrivialConditionRetriever::sampleMaskEE_ [private] |
Definition at line 288 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::timeCalibConstantMean_ [private] |
Definition at line 200 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::timeCalibConstantsFile_ [private] |
Definition at line 277 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::timeCalibConstantSigma_ [private] |
Definition at line 201 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::timeCalibErrorMean_ [private] |
Definition at line 202 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::timeCalibErrorsFile_ [private] |
Definition at line 278 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::timeOffsetEBConstant_ [private] |
Definition at line 204 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::timeOffsetEEConstant_ [private] |
Definition at line 205 of file EcalTrivialConditionRetriever.h.
double EcalTrivialConditionRetriever::totLumi_ [private] |
Definition at line 215 of file EcalTrivialConditionRetriever.h.
std::string EcalTrivialConditionRetriever::trgChannelStatusFile_ [private] |
Definition at line 280 of file EcalTrivialConditionRetriever.h.
int EcalTrivialConditionRetriever::verbose_ [private] |
Definition at line 327 of file EcalTrivialConditionRetriever.h.
Definition at line 293 of file EcalTrivialConditionRetriever.h.
1.7.3