Inheritance diagram for SiPixelLorentzAngleCalibration:

Public Member Functions
void	beginOfJob (AlignableTracker tracker, AlignableMuon muon, AlignableExtras *extras) override
	Call at beginning of job: More...

void	beginRun (const edm::Run &, const edm::EventSetup &) override
	Call at beginning of run: More...

unsigned int	derivatives (std::vector< ValuesIndexPair > &outDerivInds, const TransientTrackingRecHit &hit, const TrajectoryStateOnSurface &tsos, const edm::EventSetup &setup, const EventInfo &eventInfo) const override

void	endOfJob () override

double	getParameter (unsigned int index) const override

double	getParameterError (unsigned int index) const override

unsigned int	numParameters () const override
	How many parameters does this calibration define? More...

bool	setParameter (unsigned int index, double value) override

bool	setParameterError (unsigned int index, double error) override

	SiPixelLorentzAngleCalibration (const edm::ParameterSet &cfg)
	Constructor. More...

	~SiPixelLorentzAngleCalibration () override=default
	Destructor. More...

Public Member Functions inherited from IntegratedCalibrationBase
virtual std::vector< Values >	derivatives (const TrackingRecHit &hit, const TrajectoryStateOnSurface &tsos, const edm::EventSetup &setup, const EventInfo &eventInfo) const

virtual unsigned int	derivatives (std::vector< ValuesIndexPair > &outDerivInds, const TrackingRecHit &hit, const TrajectoryStateOnSurface &tsos, const edm::EventSetup &setup, const EventInfo &eventInfo) const =0

virtual void	endOfLoop ()

	IntegratedCalibrationBase (const edm::ParameterSet &cfg)
	Values and their parameter index. More...

const std::string &	name () const
	name of this calibration More...

virtual void	startNewLoop ()

virtual	~IntegratedCalibrationBase ()
	Destructor. More...

Private Member Functions
SiPixelLorentzAngle	createFromTree (const char fileName, const char treeName) const

const SiPixelLorentzAngle *	getLorentzAnglesInput (const align::RunNumber &=0)

double	getParameterForDetId (unsigned int detId, edm::RunNumber_t run) const

void	writeTree (const SiPixelLorentzAngle lorentzAngle, const std::map< unsigned int, TreeStruct > &treeInfo, const char treeName) const

Private Attributes
std::map< align::RunNumber, SiPixelLorentzAngle >	cachedLorentzAngleInputs_

align::RunNumber	currentIOV_ {0}

const std::string	lorentzAngleLabel_

const std::vector< std::string >	mergeFileNames_

const edm::ParameterSet	moduleGroupSelCfg_

std::unique_ptr< TkModuleGroupSelector >	moduleGroupSelector_

const std::string	outFileName_

std::vector< double >	parameters_

std::vector< double >	paramUncertainties_

const std::string	recordNameDBwrite_

const bool	saveToDB_

SiPixelLorentzAngle *	siPixelLorentzAngleInput_ {nullptr}

edm::ESWatcher< SiPixelLorentzAngleRcd >	watchLorentzAngleRcd_

Additional Inherited Members
Public Types inherited from IntegratedCalibrationBase
typedef AlignmentAlgorithmBase::EventInfo	EventInfo

typedef std::pair< double, double >	Values

typedef std::pair< Values, unsigned int >	ValuesIndexPair
	x- and y-values More...

Detailed Description

Calibration of Lorentz angle for the pixel, integrated in the alignment algorithms.

Note that not all algorithms support this...

Author: : Gero Flucke date : September 2012

Revision: 1.4.2.21

Date: 2013/05/31 08:37:12

(last update by

Author: flucke

)

Definition at line 54 of file SiPixelLorentzAngleCalibration.cc.

Constructor & Destructor Documentation

◆ SiPixelLorentzAngleCalibration()

SiPixelLorentzAngleCalibration::SiPixelLorentzAngleCalibration ( const edm::ParameterSet & cfg )

explicit

Constructor.

Definition at line 137 of file SiPixelLorentzAngleCalibration.cc.

     : IntegratedCalibrationBase(cfg),
       saveToDB_(cfg.getParameter<bool>("saveToDB")),
       recordNameDBwrite_(cfg.getParameter<std::string>("recordNameDBwrite")),
       outFileName_(cfg.getParameter<std::string>("treeFile")),
       mergeFileNames_(cfg.getParameter<std::vector<std::string> >("mergeTreeFiles")),
       lorentzAngleLabel_(cfg.getParameter<std::string>("lorentzAngleLabel")),
       moduleGroupSelCfg_(cfg.getParameter<edm::ParameterSet>("LorentzAngleModuleGroups")) {}

◆ ~SiPixelLorentzAngleCalibration()

SiPixelLorentzAngleCalibration::~SiPixelLorentzAngleCalibration ( )

overridedefault

Destructor.

Member Function Documentation

◆ beginOfJob()

void SiPixelLorentzAngleCalibration::beginOfJob	(	AlignableTracker *	tracker,
		AlignableMuon *	muon,
		AlignableExtras *	extras
	)

overridevirtual

Call at beginning of job:

Reimplemented from IntegratedCalibrationBase.

Definition at line 259 of file SiPixelLorentzAngleCalibration.cc.

                                                                      {
   //specify the sub-detectors for which the LA is determined
   const std::vector<int> sdets = {PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap};
  
   moduleGroupSelector_ = std::make_unique<TkModuleGroupSelector>(aliTracker, moduleGroupSelCfg_, sdets);
  
   parameters_.resize(moduleGroupSelector_->getNumberOfParameters(), 0.);
   paramUncertainties_.resize(moduleGroupSelector_->getNumberOfParameters(), 0.);
  
   edm::LogInfo("Alignment") << "@SUB=SiPixelLorentzAngleCalibration"
                             << "Created with name " << this->name() << "',\n"
                             << this->numParameters() << " parameters to be determined,"
                             << "\n saveToDB = " << saveToDB_ << "\n outFileName = " << outFileName_
                             << "\n N(merge files) = " << mergeFileNames_.size()
                             << "\n number of IOVs = " << moduleGroupSelector_->numIovs();
  
   if (!mergeFileNames_.empty()) {
     edm::LogInfo("Alignment") << "@SUB=SiPixelLorentzAngleCalibration"
                               << "First file to merge: " << mergeFileNames_[0];
   }
 }

References mergeFileNames_, moduleGroupSelCfg_, moduleGroupSelector_, IntegratedCalibrationBase::name(), numParameters(), outFileName_, parameters_, paramUncertainties_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, and saveToDB_.

◆ beginRun()

void SiPixelLorentzAngleCalibration::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

overridevirtual

Call at beginning of run:

Reimplemented from IntegratedCalibrationBase.

Definition at line 150 of file SiPixelLorentzAngleCalibration.cc.

                                                                                            {
   // no action needed if the LA record didn't change
   if (!(watchLorentzAngleRcd_.check(setup)))
     return;
  
   const auto runNumber = run.run();
   auto firstRun = cond::timeTypeSpecs[cond::runnumber].beginValue;
  
   // avoid infinite loop due to wrap-around of unsigned variable 'i' including
   // arrow from i to zero and a nice smiley ;)
   for (unsigned int i = moduleGroupSelector_->numIovs(); i-- > 0;) {
     const auto firstRunOfIOV = moduleGroupSelector_->firstRunOfIOV(i);
     if (runNumber >= firstRunOfIOV) {
       firstRun = firstRunOfIOV;
       break;
     }
   }
  
   edm::ESHandle<SiPixelLorentzAngle> lorentzAngleHandle;
   const auto &lorentzAngleRcd = setup.get<SiPixelLorentzAngleRcd>();
   lorentzAngleRcd.get(lorentzAngleLabel_, lorentzAngleHandle);
   if (cachedLorentzAngleInputs_.find(firstRun) == cachedLorentzAngleInputs_.end()) {
     cachedLorentzAngleInputs_.emplace(firstRun, SiPixelLorentzAngle(*lorentzAngleHandle));
   } else {
     if (lorentzAngleRcd.validityInterval().first().eventID().run() > firstRun &&
         lorentzAngleHandle->getLorentzAngles()                            // only bad if non-identical values
             != cachedLorentzAngleInputs_[firstRun].getLorentzAngles()) {  // (comparing maps)
       // Maps are containers sorted by key, but comparison problems may arise from
       // 'floating point comparison' problems (FIXME?)
       throw cms::Exception("BadInput") << "Trying to cache SiPixelLorentzAngle payload for a run (" << runNumber
                                        << ") in an IOV (" << firstRun << ") that was already cached.\n"
                                        << "The following record in your input database tag has an IOV "
                                        << "boundary that does not match your IOV definition:\n"
                                        << " - SiPixelLorentzAngleRcd '" << lorentzAngleRcd.key().name() << "' (since "
                                        << lorentzAngleRcd.validityInterval().first().eventID().run() << ")\n";
     }
   }
  
   siPixelLorentzAngleInput_ = &(cachedLorentzAngleInputs_[firstRun]);
   currentIOV_ = firstRun;
 }

References cond::TimeTypeSpecs::beginValue, cachedLorentzAngleInputs_, edm::ESWatcher< T >::check(), currentIOV_, Exception, dataset::firstRun, edm::eventsetup::EventSetupRecordImplementation< T >::get(), SiPixelLorentzAngle::getLorentzAngles(), mps_fire::i, lorentzAngleLabel_, moduleGroupSelector_, writedatasetfile::run, cond::runnumber, convertSQLiteXML::runNumber, singleTopDQM_cfi::setup, siPixelLorentzAngleInput_, cond::timeTypeSpecs, and watchLorentzAngleRcd_.

◆ createFromTree()

SiPixelLorentzAngle SiPixelLorentzAngleCalibration::createFromTree	(	const char *	fileName,
		const char *	treeName
	)		const

private

Definition at line 440 of file SiPixelLorentzAngleCalibration.cc.

                                                                                                                    {
   // Check for file existence on your own to work around
   // https://hypernews.cern.ch/HyperNews/CMS/get/swDevelopment/2715.html:
   TFile *file = nullptr;
   FILE *testFile = fopen(fileName, "r");
   if (testFile) {
     fclose(testFile);
     file = TFile::Open(fileName, "READ");
   }  // else not existing, see error below
  
   TTree *tree = nullptr;
   if (file)
     file->GetObject(treeName, tree);
  
   SiPixelLorentzAngle result{};
   if (tree) {
     unsigned int id = 0;
     float value = 0.;
     tree->SetBranchAddress("detId", &id);
     tree->SetBranchAddress("value", &value);
  
     const Long64_t nEntries = tree->GetEntries();
     for (Long64_t iEntry = 0; iEntry < nEntries; ++iEntry) {
       tree->GetEntry(iEntry);
       result.putLorentzAngle(id, value);
     }
   } else {  // Warning only since could be parallel job on no events.
     edm::LogWarning("Alignment") << "@SUB=SiPixelLorentzAngleCalibration::createFromTree"
                                  << "Could not get TTree '" << treeName << "' from file '" << fileName
                                  << (file ? "'." : "' (file does not exist).");
   }
  
   delete file;  // tree will vanish with file
   return result;
 }

References FrontierConditions_GlobalTag_cff::file, MillePedeFileConverter_cfg::fileName, and mps_fire::result.

Referenced by getLorentzAnglesInput().

◆ derivatives()

unsigned int SiPixelLorentzAngleCalibration::derivatives	(	std::vector< ValuesIndexPair > &	outDerivInds,
		const TransientTrackingRecHit &	hit,
		const TrajectoryStateOnSurface &	tsos,
		const edm::EventSetup &	setup,
		const EventInfo &	eventInfo
	)		const

override

Return non-zero derivatives for x- and y-measurements with their indices by reference. Return value is their number.

Definition at line 193 of file SiPixelLorentzAngleCalibration.cc.

                                                                                            {
   outDerivInds.clear();
  
   if (hit.det()) {  // otherwise 'constraint hit' or whatever
  
     const int index =
         moduleGroupSelector_->getParameterIndexFromDetId(hit.det()->geographicalId(), eventInfo.eventId().run());
     if (index >= 0) {  // otherwise not treated
       edm::ESHandle<MagneticField> magneticField;
       setup.get<IdealMagneticFieldRecord>().get(magneticField);
       const GlobalVector bField(magneticField->inTesla(hit.det()->surface().position()));
       const LocalVector bFieldLocal(hit.det()->surface().toLocal(bField));
       const double dZ = hit.det()->surface().bounds().thickness();  // it is a float only...
       // shift due to LA: dx = tan(LA) * dz/2 = mobility * B_y * dz/2,
       // shift due to LA: dy = - mobility * B_x * dz/2,
       // '-' since we have derivative of the residual r = trk -hit
       const double xDerivative = bFieldLocal.y() * dZ * -0.5;  // parameter is mobility!
       const double yDerivative = bFieldLocal.x() * dZ * 0.5;   // parameter is mobility!
       if (xDerivative || yDerivative) {                        // If field is zero, this is zero: do not return it
         const Values derivs{xDerivative, yDerivative};
         outDerivInds.push_back(ValuesIndexPair(derivs, index));
       }
     }
   } else {
     edm::LogWarning("Alignment") << "@SUB=SiPixelLorentzAngleCalibration::derivatives2"
                                  << "Hit without GeomDet, skip!";
   }
  
   return outDerivInds.size();
 }

References Calorimetry_cff::bField, ZMuMuAnalysisNtupler_cff::eventInfo, get, HLT_2018_cff::magneticField, moduleGroupSelector_, singleTopDQM_cfi::setup, and hit::y.

◆ endOfJob()

void SiPixelLorentzAngleCalibration::endOfJob ( )

overridevirtual

Called at end of a the job of the AlignmentProducer. Write out determined parameters.

Implements IntegratedCalibrationBase.

Definition at line 284 of file SiPixelLorentzAngleCalibration.cc.

                                               {
   // loginfo output
   std::ostringstream out;
   out << "Parameter results\n";
   for (unsigned int iPar = 0; iPar < parameters_.size(); ++iPar) {
     out << iPar << ": " << parameters_[iPar] << " +- " << paramUncertainties_[iPar] << "\n";
   }
   edm::LogInfo("Alignment") << "@SUB=SiPixelLorentzAngleCalibration::endOfJob" << out.str();
  
   std::map<unsigned int, TreeStruct> treeInfo;  // map of TreeStruct for each detId
  
   // now write 'input' tree(s)
   const std::string treeName{this->name() + '_'};
   std::vector<const SiPixelLorentzAngle *> inputs{};
   inputs.reserve(moduleGroupSelector_->numIovs());
   for (unsigned int iIOV = 0; iIOV < moduleGroupSelector_->numIovs(); ++iIOV) {
     const auto firstRunOfIOV = moduleGroupSelector_->firstRunOfIOV(iIOV);
     inputs.push_back(this->getLorentzAnglesInput(firstRunOfIOV));  // never NULL
     this->writeTree(inputs.back(),
                     treeInfo,
                     (treeName + "input_" + std::to_string(firstRunOfIOV)).c_str());  // empty treeInfo for input...
  
     if (inputs.back()->getLorentzAngles().empty()) {
       edm::LogError("Alignment") << "@SUB=SiPixelLorentzAngleCalibration::endOfJob"
                                  << "Input Lorentz angle map is empty, skip writing output!";
       return;
     }
   }
  
   const unsigned int nonZeroParamsOrErrors =  // Any determined value?
       count_if(parameters_.begin(), parameters_.end(), [](auto c) { return c != 0.; }) +
       count_if(paramUncertainties_.begin(), paramUncertainties_.end(), [](auto c) { return c != 0.; });
  
   for (unsigned int iIOV = 0; iIOV < moduleGroupSelector_->numIovs(); ++iIOV) {
     auto firstRunOfIOV = static_cast<cond::Time_t>(moduleGroupSelector_->firstRunOfIOV(iIOV));
     SiPixelLorentzAngle output{};
     // Loop on map of values from input and add (possible) parameter results
     for (const auto &iterIdValue : inputs[iIOV]->getLorentzAngles()) {
       // type of 'iterIdValue' is pair<unsigned int, float>
       const auto detId = iterIdValue.first;  // key of map is DetId
       const auto param = this->getParameterForDetId(detId, firstRunOfIOV);
       // put result in output, i.e. sum of input and determined parameter, but the nasty
       // putLorentzAngle(..) takes float by reference - not even const reference:
       auto value = iterIdValue.second + static_cast<float>(param);
       output.putLorentzAngle(detId, value);
       const int paramIndex = moduleGroupSelector_->getParameterIndexFromDetId(detId, firstRunOfIOV);
       treeInfo[detId] = TreeStruct(param, this->getParameterError(paramIndex), paramIndex);
     }
  
     if (saveToDB_ || nonZeroParamsOrErrors != 0) {  // Skip writing mille jobs...
       this->writeTree(&output, treeInfo, (treeName + Form("result_%lld", firstRunOfIOV)).c_str());
     }
  
     if (saveToDB_) {  // If requested, write out to DB
       edm::Service<cond::service::PoolDBOutputService> dbService;
       if (dbService.isAvailable()) {
         dbService->writeOne(&output, firstRunOfIOV, recordNameDBwrite_);
       } else {
         edm::LogError("BadConfig") << "@SUB=SiPixelLorentzAngleCalibration::endOfJob"
                                    << "No PoolDBOutputService available, but saveToDB true!";
       }
     }
   }  // end loop on IOVs
 }

References HltBtagPostValidation_cff::c, getLorentzAnglesInput(), getParameterError(), getParameterForDetId(), PixelMapPlotter::inputs, edm::Service< T >::isAvailable(), moduleGroupSelector_, IntegratedCalibrationBase::name(), MillePedeFileConverter_cfg::out, convertSQLitetoXML_cfg::output, parameters_, paramUncertainties_, recordNameDBwrite_, saveToDB_, AlCaHLTBitMon_QueryRunRegistry::string, cond::service::PoolDBOutputService::writeOne(), and writeTree().

◆ getLorentzAnglesInput()

const SiPixelLorentzAngle * SiPixelLorentzAngleCalibration::getLorentzAnglesInput ( const align::RunNumber & run = 0 )

private

Input LorentzAngle values:

either from EventSetup of first call to derivatives(..)
or created from files of passed by configuration (i.e. from parallel processing)

Definition at line 350 of file SiPixelLorentzAngleCalibration.cc.

                                                                                                           {
   const auto &resolvedRun = run > 0 ? run : currentIOV_;
   // For parallel processing in Millepede II, create SiPixelLorentzAngle
   // from info stored in files of parallel jobs and check that they are identical.
   // If this job has run on data, still check that LA is identical to the ones
   // from mergeFileNames_.
   const std::string treeName{this->name() + "_input_" + std::to_string(resolvedRun)};
   for (const auto &iFile : mergeFileNames_) {
     auto la = this->createFromTree(iFile.c_str(), treeName.c_str());
     // siPixelLorentzAngleInput_ could be non-null from previous file of this loop
     // or from checkLorentzAngleInput(..) when running on data in this job as well
     if (!siPixelLorentzAngleInput_ || siPixelLorentzAngleInput_->getLorentzAngles().empty()) {
       cachedLorentzAngleInputs_[resolvedRun] = la;
       siPixelLorentzAngleInput_ = &(cachedLorentzAngleInputs_[resolvedRun]);
       currentIOV_ = resolvedRun;
     } else {
       // FIXME: about comparison of maps see comments in checkLorentzAngleInput
       if (!la.getLorentzAngles().empty() &&  // single job might not have got events
           la.getLorentzAngles() != siPixelLorentzAngleInput_->getLorentzAngles()) {
         // Throw exception instead of error?
         edm::LogError("NoInput") << "@SUB=SiPixelLorentzAngleCalibration::getLorentzAnglesInput"
                                  << "Different input values from tree " << treeName << " in file " << iFile << ".";
       }
     }
   }
  
   if (!siPixelLorentzAngleInput_) {  // no files nor ran on events
     // [] operator default-constructs an empty SiPixelLorentzAngle object in place:
     siPixelLorentzAngleInput_ = &(cachedLorentzAngleInputs_[resolvedRun]);
     currentIOV_ = resolvedRun;
     edm::LogError("NoInput") << "@SUB=SiPixelLorentzAngleCalibration::getLorentzAnglesInput"
                              << "No input, create an empty one!";
   } else if (siPixelLorentzAngleInput_->getLorentzAngles().empty()) {
     edm::LogError("NoInput") << "@SUB=SiPixelLorentzAngleCalibration::getLorentzAnglesInput"
                              << "Empty result!";
   }
  
   return siPixelLorentzAngleInput_;
 }

References cachedLorentzAngleInputs_, createFromTree(), currentIOV_, SiPixelLorentzAngle::getLorentzAngles(), mergeFileNames_, IntegratedCalibrationBase::name(), writedatasetfile::run, siPixelLorentzAngleInput_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by endOfJob().

◆ getParameter()

double SiPixelLorentzAngleCalibration::getParameter ( unsigned int index ) const

overridevirtual

Return current value of parameter identified by index. Return 0. if index out-of-bounds.

Implements IntegratedCalibrationBase.

Definition at line 249 of file SiPixelLorentzAngleCalibration.cc.

                                                                             {
   return (index >= parameters_.size() ? 0. : parameters_[index]);
 }

References parameters_.

Referenced by Mixins._Parameterizable::hasParameter().

◆ getParameterError()

double SiPixelLorentzAngleCalibration::getParameterError ( unsigned int index ) const

overridevirtual

Return current value of parameter identified by index. Returns 0. if index out-of-bounds or if errors undetermined.

Implements IntegratedCalibrationBase.

Definition at line 254 of file SiPixelLorentzAngleCalibration.cc.

                                                                                  {
   return (index >= paramUncertainties_.size() ? 0. : paramUncertainties_[index]);
 }

References paramUncertainties_.

Referenced by endOfJob().

◆ getParameterForDetId()

double SiPixelLorentzAngleCalibration::getParameterForDetId	(	unsigned int	detId,
		edm::RunNumber_t	run
	)		const

private

Determined parameter value for this detId (detId not treated => 0.) and the given run.

Definition at line 391 of file SiPixelLorentzAngleCalibration.cc.

                                                                                                         {
   const int index = moduleGroupSelector_->getParameterIndexFromDetId(detId, run);
   return (index < 0 ? 0. : parameters_[index]);
 }

References moduleGroupSelector_, parameters_, and writedatasetfile::run.

Referenced by endOfJob().

◆ numParameters()

unsigned int SiPixelLorentzAngleCalibration::numParameters ( ) const

overridevirtual

How many parameters does this calibration define?

Implements IntegratedCalibrationBase.

Definition at line 147 of file SiPixelLorentzAngleCalibration.cc.

147 { return parameters_.size(); }

References parameters_.

Referenced by beginOfJob().

◆ setParameter()

bool SiPixelLorentzAngleCalibration::setParameter	(	unsigned int	index,
		double	value
	)

overridevirtual

Setting the determined parameter identified by index, returns false if out-of-bounds, true otherwise.

Implements IntegratedCalibrationBase.

Definition at line 229 of file SiPixelLorentzAngleCalibration.cc.

                                                                                   {
   if (index >= parameters_.size()) {
     return false;
   } else {
     parameters_[index] = value;
     return true;
   }
 }

References parameters_, and relativeConstraints::value.

Referenced by coreTools.RunOnData::__call__(), trackTools.MakeAODTrackCandidates::__call__(), runJetUncertainties.RunJetUncertainties::__call__(), metTools.AddMETCollection::__call__(), cmsswVersionTools.PickRelValInputFiles::__call__(), coreTools.RemoveMCMatching::__call__(), trackTools.MakePATTrackCandidates::__call__(), trigTools.SwitchOnTrigger::__call__(), runMETCorrectionsAndUncertainties.RunMETCorrectionsAndUncertainties::__call__(), trigTools.SwitchOnTriggerStandAlone::__call__(), tauTools.AddTauCollection::__call__(), trackTools.MakeTrackCandidates::__call__(), trigTools.SwitchOnTriggerMatching::__call__(), trigTools.SwitchOnTriggerMatchingStandAlone::__call__(), trigTools.SwitchOnTriggerMatchEmbedding::__call__(), jetTools.AddJetCollection::__call__(), jetTools.SwitchJetCollection::__call__(), jetTools.UpdateJetCollection::__call__(), jetTools.AddJetID::__call__(), jetTools.SetTagInfos::__call__(), and runMETCorrectionsAndUncertainties.RunMETCorrectionsAndUncertainties::jetConfiguration().

◆ setParameterError()

bool SiPixelLorentzAngleCalibration::setParameterError	(	unsigned int	index,
		double	error
	)

overridevirtual

Setting the determined parameter uncertainty identified by index, returns false if out-of-bounds, true otherwise.

Implements IntegratedCalibrationBase.

Definition at line 239 of file SiPixelLorentzAngleCalibration.cc.

                                                                                        {
   if (index >= paramUncertainties_.size()) {
     return false;
   } else {
     paramUncertainties_[index] = error;
     return true;
   }
 }

References relativeConstraints::error, and paramUncertainties_.

◆ writeTree()

void SiPixelLorentzAngleCalibration::writeTree	(	const SiPixelLorentzAngle *	lorentzAngle,
		const std::map< unsigned int, TreeStruct > &	treeInfo,
		const char *	treeName
	)		const

private

Definition at line 397 of file SiPixelLorentzAngleCalibration.cc.

                                                                            {
   if (!lorentzAngle)
     return;
  
   TFile *file = TFile::Open(outFileName_.c_str(), "UPDATE");
   if (!file) {
     edm::LogError("BadConfig") << "@SUB=SiPixelLorentzAngleCalibration::writeTree"
                                << "Could not open file '" << outFileName_ << "'.";
     return;
   }
  
   TTree *tree = new TTree(treeName, treeName);
   unsigned int id = 0;
   float value = 0.;
   TreeStruct treeStruct;
   tree->Branch("detId", &id, "detId/i");
   tree->Branch("value", &value, "value/F");
   tree->Branch("treeStruct", &treeStruct, TreeStruct::LeafList());
  
   for (auto iterIdValue = lorentzAngle->getLorentzAngles().begin();
        iterIdValue != lorentzAngle->getLorentzAngles().end();
        ++iterIdValue) {
     // type of (*iterIdValue) is pair<unsigned int, float>
     id = iterIdValue->first;  // key of map is DetId
     value = iterIdValue->second;
     // type of (*treeStructIter) is pair<unsigned int, TreeStruct>
     auto treeStructIter = treeInfo.find(id);  // find info for this id
     if (treeStructIter != treeInfo.end()) {
       treeStruct = treeStructIter->second;  // info from input map
     } else {                                // if none found, fill at least parameter index (using 1st IOV...)
       const cond::Time_t run1of1stIov = moduleGroupSelector_->firstRunOfIOV(0);
       const int ind = moduleGroupSelector_->getParameterIndexFromDetId(id, run1of1stIov);
       treeStruct = TreeStruct(ind);
     }
     tree->Fill();
   }
   tree->Write();
   delete file;  // tree vanishes with the file... (?)
 }