A tool to associate SpaceVoxels with Calibrator objects. More...

#include <SpaceManager.h>

Public Member Functions
void	addEvolution (const DetectorElementPtr &det, Region region, const TF1 &f)

void	assignCalibration (const CalibratorPtr &c, const std::map< DetectorElementPtr, double > &result)

void	clear ()

CalibratorPtr	createCalibrator (const Calibrator &toClone, SpaceVoxelPtr s)

void	createCalibrators (const Calibrator &toClone, const double etaSeg, const double phiSeg, const double energySeg)

void	createCalibrators (const Calibrator &toClone, const unsigned nEta, const double etaMin, const double etaMax, const unsigned nPhi, const double phiMin, const double phiMax, const unsigned nEnergy, const double energyMin, const double energyMax) noexcept(false)

void	createCalibrators (const Calibrator &toClone)

double	evolveCoefficient (DetectorElementPtr det, double energy, double eta, double phi)

TH1 *	extractEvolution (DetectorElementPtr det, Region region, TF1 &f1, bool useTruth=true)

CalibratorPtr	findCalibrator (const double eta, const double phi, const double energy=0) const

std::map< DetectorElementPtr, double >	getCalibration (CalibratorPtr c)

std::map< SpaceVoxelPtr, CalibratorPtr > *	getCalibrators ()

std::string	getName ()

int	getNCalibrations ()

double	interpolateCoefficient (DetectorElementPtr det, double energy, double eta, double phi)

void	makeInverseAddressBook ()

std::ostream &	printCalibrations (std::ostream &stream)

void	setBarrelLimit (double limit)

	SpaceManager (std::string name)

virtual	~SpaceManager ()

Private Attributes
double	barrelLimit_

std::map< DetectorElementPtr, TF1 >	barrelPosEvolutions_

std::vector< SpaceVoxelPtr >	barrelPosRegion_

std::map< CalibratorPtr, std::map< DetectorElementPtr, double > >	calibrationCoeffs_

double	endcapLimit_

std::map< DetectorElementPtr, TF1 >	endcapPosEvolutions_

std::vector< SpaceVoxelPtr >	endcapPosRegion_

std::map< CalibratorPtr, SpaceVoxelPtr >	inverseAddressBook_

std::map< SpaceVoxelPtr, CalibratorPtr >	myAddressBook

std::vector< SpaceVoxelPtr >	myKnownSpaceVoxels

std::string	name_

std::map< Region, std::vector< SpaceVoxelPtr > >	regionsToSVs_

double	transitionLimit_

std::vector< SpaceVoxelPtr >	transitionPosRegion_

Detailed Description

A tool to associate SpaceVoxels with Calibrator objects.

Author: Jamie Ballin

Date: April 2008

Definition at line 23 of file SpaceManager.h.

Constructor & Destructor Documentation

SpaceManager::SpaceManager ( std::string name )

Definition at line 12 of file SpaceManager.cc.

References pftools::BARREL_POS, barrelPosRegion_, pftools::ENDCAP_POS, endcapPosRegion_, and regionsToSVs_.

                                          :
     name_(name), barrelLimit_(1.4), transitionLimit_(1.4), endcapLimit_(5.0) {
     regionsToSVs_[BARREL_POS] = barrelPosRegion_;
     regionsToSVs_[ENDCAP_POS] = endcapPosRegion_;
 }

SpaceManager::~SpaceManager ( )

virtual

Definition at line 18 of file SpaceManager.cc.

                             {
 
 }

Member Function Documentation

void pftools::SpaceManager::addEvolution	(	const DetectorElementPtr &	det,
		Region	region,
		const TF1 &	f
	)

inline

Definition at line 79 of file SpaceManager.h.

References pftools::BARREL_POS, barrelPosEvolutions_, pftools::ENDCAP_POS, endcapPosEvolutions_, evolveCoefficient(), f, and interpolateCoefficient().

                                                                                   {
         if(region == BARREL_POS)
             barrelPosEvolutions_[det] = f;
         if(region == ENDCAP_POS)
             endcapPosEvolutions_[det] = f;
     }

void SpaceManager::assignCalibration	(	const CalibratorPtr &	c,
		const std::map< DetectorElementPtr, double > &	result
	)

Definition at line 246 of file SpaceManager.cc.

References EnergyCorrector::c, calibrationCoeffs_, makeInverseAddressBook(), and mps_fire::result.

Referenced by getCalibrators().

                                                           {
     calibrationCoeffs_[c] = result;
     makeInverseAddressBook();
 }

void SpaceManager::clear ( void )

Definition at line 22 of file SpaceManager.cc.

References EnergyCorrector::c, myAddressBook, and alignCSCRings::s.

Referenced by createCalibrators(), and getNCalibrations().

                          {
     for (std::map<SpaceVoxelPtr, CalibratorPtr>::iterator it =
             myAddressBook.begin(); it!= myAddressBook.end(); ++it) {
         SpaceVoxelPtr s = (*it).first;
         CalibratorPtr c = (*it).second;
     }
 }

CalibratorPtr SpaceManager::createCalibrator	(	const Calibrator &	toClone,
		SpaceVoxelPtr	s
	)

Definition at line 213 of file SpaceManager.cc.

References EnergyCorrector::c, pftools::Calibrator::clone(), KineDebug3::count(), myAddressBook, myKnownSpaceVoxels, and alignCSCRings::s.

Referenced by getCalibrators().

                          {
     CalibratorPtr c;
     int known = count(myKnownSpaceVoxels.begin(), myKnownSpaceVoxels.end(), s);
     if (known == 0) {
         myKnownSpaceVoxels.push_back(s);
         c.reset(toClone.clone());
         myAddressBook[s] = c;
     } else {
         c = myAddressBook[s];
     }
     assert(c != 0);
     return c;
 
 }

void SpaceManager::createCalibrators	(	const Calibrator &	toClone,
		const double	etaSeg,
		const double	phiSeg,
		const double	energySeg
	)

Definition at line 30 of file SpaceManager.cc.

References gather_cfg::cout.

Referenced by getName().

                                                                           {
     std::cout << __PRETTY_FUNCTION__
             << ": this method has not yet been implemented!\n";
     PFToolsException me("Unimplemented method! Sorry!");
     throw me;
 
 }

void SpaceManager::createCalibrators	(	const Calibrator &	toClone,
		const unsigned	nEta,
		const double	etaMin,
		const double	etaMax,
		const unsigned	nPhi,
		const double	phiMin,
		const double	phiMax,
		const unsigned	nEnergy,
		const double	energyMin,
		const double	energyMax
	)

noexcept

Definition at line 168 of file SpaceManager.cc.

References EnergyCorrector::c, clear(), gather_cfg::cout, eleEcalExtractorBlocks_cff::energyMin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, gen::k, checklumidiff::l, funct::m, makeInverseAddressBook(), myAddressBook, myKnownSpaceVoxels, DetIdAssociatorESProducer_cff::nEta, DetIdAssociatorESProducer_cff::nPhi, AlignmentTrackSelector_cfi::phiMax, and AlignmentTrackSelector_cfi::phiMin.

                         {
     clear();
 
     if (nEta == 0|| nPhi ==0|| nEnergy == 0) {
         PFToolsException
                 me("Can't create calibrators with zero values for nEta, nPhi or nEnergy!");
         throw me;
     }
 
     double etaSeg = (etaMax - etaMin) / nEta;
     double phiSeg = (phiMax - phiMin) / nPhi;
     double energySeg = (energyMax - energyMin) / nEnergy;
 
     double eta1, eta2, phi1, phi2, energy1, energy2;
     for (unsigned k(0); k < nEta; ++k) {
         for (unsigned l(0); l < nPhi; ++l) {
             for (unsigned m(0); m < nEnergy; ++m) {
                 eta1 = etaMin + k * etaSeg;
                 eta2 = eta1 + etaSeg;
 
                 phi1 = phiMin + l * phiSeg;
                 phi2 = phi1 + phiSeg;
 
                 energy1 = energyMin + m * energySeg;
                 energy2 = energy1 + energySeg;
                 SpaceVoxelPtr sv(new SpaceVoxel(eta1, eta2, phi1, phi2, energy1, energy2));
                 myKnownSpaceVoxels.push_back(sv);
                 CalibratorPtr c(toClone.clone());
                 myAddressBook[sv] = c;
             }
         }
     }
     unsigned nCalibrators = nEta * nPhi * nEnergy;
     std::cout << "Created "<< nCalibrators << " calibrators.\n";
     std::cout << "Address book size: \t\t"<< myAddressBook.size() << "\n";
     std::cout << "Known space voxels size: \t"<< myKnownSpaceVoxels.size()
             << "\n";
     assert(myAddressBook.size() == myKnownSpaceVoxels.size());
     makeInverseAddressBook();
 
 }

void SpaceManager::createCalibrators ( const Calibrator & toClone )

Definition at line 39 of file SpaceManager.cc.

References barrelLimit_, barrelPosRegion_, EnergyCorrector::c, clear(), pftools::Calibrator::clone(), gather_cfg::cout, endcapLimit_, endcapPosRegion_, myAddressBook, and myKnownSpaceVoxels.

                                                               {
     clear();
     std::cout << __PRETTY_FUNCTION__
             << ": creating default calibration schema.\n";
 
     SpaceVoxelPtr sv(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 0, 3.0));
     barrelPosRegion_.push_back(sv);
     SpaceVoxelPtr sv1(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 3, 9.0));
     barrelPosRegion_.push_back(sv1);
     SpaceVoxelPtr sv2(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 9.0, 16.0));
     barrelPosRegion_.push_back(sv2);
     SpaceVoxelPtr sv3(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 16.0, 25.0));
     barrelPosRegion_.push_back(sv3);
     SpaceVoxelPtr sv4(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 25.0, 100.0));
     barrelPosRegion_.push_back(sv4);
     SpaceVoxelPtr sv5(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 100.0, 200.0));
     barrelPosRegion_.push_back(sv5);
     SpaceVoxelPtr sv6(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 200.0, 400.0));
     barrelPosRegion_.push_back(sv6);
 
     //  SpaceVoxelPtr sv(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 0, 2.0));
     //  barrelPosRegion_.push_back(sv);
     //  SpaceVoxelPtr sv1(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 1, 1.5));
     //  barrelPosRegion_.push_back(sv1);
     //  SpaceVoxelPtr sv2(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 1.5, 2.0));
     //  barrelPosRegion_.push_back(sv2);
     //  SpaceVoxelPtr sv3(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 2.0, 2.5));
     //  barrelPosRegion_.push_back(sv3);
     //  SpaceVoxelPtr sv4(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 2.5, 4.0));
     //  barrelPosRegion_.push_back(sv4);
     //  SpaceVoxelPtr sv41(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 4.0, 5.0));
     //  barrelPosRegion_.push_back(sv41);
     //  SpaceVoxelPtr sv5(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 5.0, 6.5));
     //  barrelPosRegion_.push_back(sv5);
     //  SpaceVoxelPtr sv51(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 6.5, 8.0));
     //  barrelPosRegion_.push_back(sv51);
     //  SpaceVoxelPtr sv6(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 8.0, 10.0));
     //  barrelPosRegion_.push_back(sv6);
     //  SpaceVoxelPtr sv61(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 10.0, 12.0));
     //  barrelPosRegion_.push_back(sv61);
     //  SpaceVoxelPtr sv7(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 12.0, 16.0));
     //  barrelPosRegion_.push_back(sv7);
     //  SpaceVoxelPtr sv8(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 16.0, 25.0));
     //  barrelPosRegion_.push_back(sv8);
     //  SpaceVoxelPtr sv9(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 25.0, 40.0));
     //  barrelPosRegion_.push_back(sv9);
     //  SpaceVoxelPtr sv10(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 40.0, 60.0));
     //  barrelPosRegion_.push_back(sv10);
     //  SpaceVoxelPtr sv11(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 60.0, 100.0));
     //  barrelPosRegion_.push_back(sv11);
     //  SpaceVoxelPtr sv12(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 100.0, 200.0));
     //  barrelPosRegion_.push_back(sv12);
     //  SpaceVoxelPtr sv13(new SpaceVoxel(-1.0*barrelLimit_, barrelLimit_, -3.2, 3.2, 200.0, 400.0));
     //  barrelPosRegion_.push_back(sv13);
 
     SpaceVoxelPtr sve(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 0, 3.0));
     endcapPosRegion_.push_back(sve);
     SpaceVoxelPtr sve0(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 3.0, 9.0));
     endcapPosRegion_.push_back(sve0);
     SpaceVoxelPtr sve1(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 9.0, 16.0));
     endcapPosRegion_.push_back(sve1);
     SpaceVoxelPtr sve2(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 16.0, 25.0));
     endcapPosRegion_.push_back(sve2);
     SpaceVoxelPtr sve3(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 25.0, 100.0));
     endcapPosRegion_.push_back(sve3);
     SpaceVoxelPtr sve4(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 100.0, 200.0));
     endcapPosRegion_.push_back(sve4);
     SpaceVoxelPtr sve5(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 200.0, 400.0));
     endcapPosRegion_.push_back(sve5);
 
     //  SpaceVoxelPtr sve(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 0, 0.5));
     //  endcapPosRegion_.push_back(sve);
     //  SpaceVoxelPtr sve0(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 0.5, 1.0));
     //  endcapPosRegion_.push_back(sve0);
     //  SpaceVoxelPtr sve1(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 1, 1.5));
     //  endcapPosRegion_.push_back(sve1);
     //  SpaceVoxelPtr sve2(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 1.5, 2.0));
     //  endcapPosRegion_.push_back(sve2);
     //  SpaceVoxelPtr sve3(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 2.0, 2.5));
     //  endcapPosRegion_.push_back(sve3);
     //  SpaceVoxelPtr sve4(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 2.5, 4.0));
     //  endcapPosRegion_.push_back(sve4);
     //  SpaceVoxelPtr sve5(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 4.0, 5.0));
     //  endcapPosRegion_.push_back(sve5);
     //  SpaceVoxelPtr sve51(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 5.0, 6.5));
     //  endcapPosRegion_.push_back(sve51);
     //  SpaceVoxelPtr sve6(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 6.5, 8.0));
     //  endcapPosRegion_.push_back(sve6);
     //  SpaceVoxelPtr sve61(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 8.0, 10.0));
     //  endcapPosRegion_.push_back(sve61);
     //  SpaceVoxelPtr sve62(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 10.0, 12.0));
     //  endcapPosRegion_.push_back(sve62);
     //  SpaceVoxelPtr sve7(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 12.0, 16.0));
     //  endcapPosRegion_.push_back(sve7);
     //  SpaceVoxelPtr sve8(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 16.0, 25.0));
     //  endcapPosRegion_.push_back(sve8);
     //  SpaceVoxelPtr sve9(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 25.0, 40.0));
     //  endcapPosRegion_.push_back(sve9);
     //  SpaceVoxelPtr sve10(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 40.0, 60.0));
     //  endcapPosRegion_.push_back(sve10);
     //  SpaceVoxelPtr sve11(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 60.0, 100.0));
     //  endcapPosRegion_.push_back(sve11);
     //  SpaceVoxelPtr sve12(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 100.0, 200.0));
     //  endcapPosRegion_.push_back(sve12);
     //  SpaceVoxelPtr sve13(new SpaceVoxel(barrelLimit_, endcapLimit_, -3.2, 3.2, 200.0, 400.0));
     //  endcapPosRegion_.push_back(sve13);
 
     for (std::vector<SpaceVoxelPtr>::iterator it = barrelPosRegion_.begin(); it
             != barrelPosRegion_.end(); ++it) {
         myKnownSpaceVoxels.push_back(*it);
     }
 
     for (std::vector<SpaceVoxelPtr>::iterator it = endcapPosRegion_.begin(); it
             != endcapPosRegion_.end(); ++it) {
         myKnownSpaceVoxels.push_back(*it);
     }
 
     for (std::vector<SpaceVoxelPtr>::iterator it = myKnownSpaceVoxels.begin(); it
             != myKnownSpaceVoxels.end(); ++it) {
         CalibratorPtr c(toClone.clone());
         myAddressBook[*it] = c;
     }
     std::cout << "Address book size: \t\t"<< myAddressBook.size() << "\n";
     std::cout << "Known space voxels size: \t"<< myKnownSpaceVoxels.size()
             << "\n";
     assert(myAddressBook.size() == myKnownSpaceVoxels.size());
 
 }

double SpaceManager::evolveCoefficient	(	DetectorElementPtr	det,
		double	energy,
		double	eta,
		double	phi
	)

Definition at line 320 of file SpaceManager.cc.

References barrelLimit_, barrelPosEvolutions_, endcapPosEvolutions_, and RecoEcal_EventContent_cff::func.

Referenced by addEvolution().

                                 {
     if (eta < barrelLimit_) {
         TF1& func = barrelPosEvolutions_[det];
         return func.Eval(energy);
     }
     TF1& func = endcapPosEvolutions_[det];
     return func.Eval(energy);
 }

TH1 * SpaceManager::extractEvolution	(	DetectorElementPtr	det,
		Region	region,
		TF1 &	f1,
		bool	useTruth = `true`
	)

Definition at line 256 of file SpaceManager.cc.

References pftools::BARREL_POS, barrelPosRegion_, EnergyCorrector::c, calibrationCoeffs_, pftools::DetElNames, pftools::ENDCAP_POS, endcapPosRegion_, mps_fire::i, particleFlow_cfi::maxE, myAddressBook, dataset::name, HadronAndPartonSelector_cfi::particles, NewPostProcessor_cff::profile, pftools::RegionNames, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by getCalibrators().

                        {
 
     std::vector<SpaceVoxelPtr> region;
     if (r == BARREL_POS)
         region = barrelPosRegion_;
     if (r == ENDCAP_POS)
         region = endcapPosRegion_;
     //region = regionsToSVs_[r];
 
     std::sort(region.begin(), region.end(), SpaceVoxel());
 
     std::string detElName = DetElNames[det->getType()];
     std::string name("hDist_");
     name.append(RegionNames[r]);
     name.append("_");
     name.append(DetElNames[det->getType()]);
 
     double minE(1000);
     double maxE(0);
 
     TH2F hDist(name.c_str(), name.c_str(), 100, 0, 300, 50, 0.0, 2.5);
     //  TH3F hSurf(nameSurf.c_str(), nameSurf.c_str(), 30, 0, 50, 10, 0.0, 3.0, 30,
     //          0.0, 2.5);
     for (std::vector<SpaceVoxelPtr>::iterator i = region.begin(); i
             != region.end(); ++i) {
         SpaceVoxelPtr s = *i;
         //double midE = (s->maxEnergy() + s->minEnergy()) / 2.0;
         if (s->maxEnergy() > maxE)
             maxE = s->maxEnergy();
         if (s->minEnergy() < minE)
             minE = s->minEnergy();
         CalibratorPtr c = myAddressBook[s];
         double coeff = calibrationCoeffs_[c][det];
         if (coeff != 0.0) {
             std::vector<ParticleDepositPtr> particles = c->getParticles();
             for (std::vector<ParticleDepositPtr>::iterator it =
                     particles.begin(); it != particles.end(); ++it) {
                 if (useTruth)
                     hDist.Fill((*it)->getTruthEnergy(), coeff);
                 else
                     hDist.Fill((*it)->getRecEnergy(), coeff);
             }
         }
     }
 
     hDist.FitSlicesY();
     hDist.ProfileX();
     hDist.Write();
     std::string nameProfile(name);
     nameProfile.append("_pfx");
     name.append("_1");
 
     TH1D* slices = (TH1D*) gDirectory->Get(name.c_str());
     //TH2D* slicesSurf = (TH2D*) gDirectory->Get(nameSurf.c_str());
     TProfile* profile = (TProfile*) gDirectory->Get(nameProfile.c_str());
     profile->Fit(&f1);
     slices->Fit(&f1);
     profile->Write();
 
     return slices;
 
 }

CalibratorPtr SpaceManager::findCalibrator	(	const double	eta,
		const double	phi,
		const double	energy = `0`
	)		const

Definition at line 229 of file SpaceManager.cc.

References submit::answer, KineDebug3::count(), myAddressBook, myKnownSpaceVoxels, alignCSCRings::s, and edm::second().

Referenced by getCalibrators(), and interpolateCoefficient().

                                    {
     CalibratorPtr answer;
     for (std::vector<SpaceVoxelPtr>::const_iterator cit =
             myKnownSpaceVoxels.begin(); cit != myKnownSpaceVoxels.end(); ++cit) {
         SpaceVoxelPtr s = *cit;
         if (s->contains(eta, phi, energy)) {
             assert(count(myKnownSpaceVoxels.begin(), myKnownSpaceVoxels.end(), s) != 0);
             answer = (*myAddressBook.find(s)).second;
             break;
         } else {
             //assert(count(myKnownSpaceVoxels.begin(), myKnownSpaceVoxels.end(), s) == 0);
         }
     }
     return answer;
 }

std::map< DetectorElementPtr, double > SpaceManager::getCalibration ( CalibratorPtr c )

Definition at line 252 of file SpaceManager.cc.

References EnergyCorrector::c, and calibrationCoeffs_.

Referenced by getCalibrators().

                                                                                {
     return calibrationCoeffs_[c];
 }

std::map<SpaceVoxelPtr, CalibratorPtr>* pftools::SpaceManager::getCalibrators ( )

inline

Definition at line 52 of file SpaceManager.h.

References assignCalibration(), EnergyCorrector::c, createCalibrator(), PVValHelper::eta, extractEvolution(), connectstrParser::f1, findCalibrator(), getCalibration(), myAddressBook, printCalibrations(), mps_fire::result, and alignCSCRings::s.

                                                            {
         std::map<SpaceVoxelPtr, CalibratorPtr>* ptr = &myAddressBook;
         return ptr;
     }

std::string pftools::SpaceManager::getName ( void )

inline

Definition at line 29 of file SpaceManager.h.

References createCalibrators(), eleEcalExtractorBlocks_cff::energyMin, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, name_, DetIdAssociatorESProducer_cff::nEta, noexcept, DetIdAssociatorESProducer_cff::nPhi, AlignmentTrackSelector_cfi::phiMax, and AlignmentTrackSelector_cfi::phiMin.

Referenced by plotting.Plot::draw().

                         {
         return name_;
     }

int pftools::SpaceManager::getNCalibrations ( )

inline

Definition at line 90 of file SpaceManager.h.

References calibrationCoeffs_, clear(), and makeInverseAddressBook().

                            {
         return calibrationCoeffs_.size();
     }

double SpaceManager::interpolateCoefficient	(	DetectorElementPtr	det,
		double	energy,
		double	eta,
		double	phi
	)

Definition at line 330 of file SpaceManager.cc.

References EnergyCorrector::c, calibrationCoeffs_, findCalibrator(), inverseAddressBook_, and alignCSCRings::s.

Referenced by addEvolution().

                                                {
     CalibratorPtr c = findCalibrator(eta, phi, energy);
 
     SpaceVoxelPtr s = inverseAddressBook_[c];
 
     double midEnergy = (s->minEnergy() + s->maxEnergy())/2.0;
     //interpolate left or right?
     double diffEnergy = energy - midEnergy;
     double thisCoeff = calibrationCoeffs_[c][det];
 
     double interpolatedCoeff = thisCoeff;
     double adjacentCoeff = 0.0;
     double adjacentEnergy = 0.0;
     if (diffEnergy > 0) {
         //look to higher energy calibrators
         CalibratorPtr adjC = findCalibrator(eta, phi, s->maxEnergy() + 0.1);
         if (adjC != 0) {
             SpaceVoxelPtr adjS = inverseAddressBook_[adjC];
             adjacentCoeff = calibrationCoeffs_[adjC][det];
             adjacentEnergy = (adjS->minEnergy() + adjS->maxEnergy()) / 2.0;
         }
     } else {
         //look to lower energy calibrations
         CalibratorPtr adjC = findCalibrator(eta, phi, s->minEnergy() - 0.1);
         if (adjC != 0) {
             SpaceVoxelPtr adjS = inverseAddressBook_[adjC];
             adjacentCoeff = calibrationCoeffs_[adjC][det];
             adjacentEnergy = (adjS->minEnergy() + adjS->maxEnergy()) / 2.0;
         }
     }
     if (adjacentCoeff != 0) {
         interpolatedCoeff = thisCoeff + diffEnergy* (adjacentCoeff - thisCoeff)
                 / (adjacentEnergy - midEnergy);
     }
     return interpolatedCoeff;
 }

void SpaceManager::makeInverseAddressBook ( )

Definition at line 393 of file SpaceManager.cc.

References EnergyCorrector::c, inverseAddressBook_, myAddressBook, and alignCSCRings::s.

Referenced by assignCalibration(), createCalibrators(), and getNCalibrations().

                                           {
     inverseAddressBook_.clear();
     for (std::map<SpaceVoxelPtr, CalibratorPtr>::iterator it =
             myAddressBook.begin(); it != myAddressBook.end(); ++it) {
         SpaceVoxelPtr s = (*it).first;
         CalibratorPtr c = (*it).second;
         inverseAddressBook_[c] = s;
     }
 }

std::ostream & SpaceManager::printCalibrations ( std::ostream & stream )

Definition at line 368 of file SpaceManager.cc.

References b, EnergyCorrector::c, calibrationCoeffs_, edmIntegrityCheck::d, myAddressBook, myKnownSpaceVoxels, mps_fire::result, and alignCSCRings::s.

Referenced by getCalibrators().

                                                               {
     stream << "Calibration results: \n";
     //  std::sort(myKnownSpaceVoxels.begin(), myKnownSpaceVoxels.end(),
     //          SpaceVoxel());
     stream << "WARNING! Haven't sorted space voxels properly!\n";
     for (std::vector<SpaceVoxelPtr>::iterator it = myKnownSpaceVoxels.begin(); it
             != myKnownSpaceVoxels.end(); ++it) {
         SpaceVoxelPtr s = *it;
         CalibratorPtr c = myAddressBook[s];
         stream << *s << "\n";
         stream << "\t[";
         std::map<DetectorElementPtr, double> result = calibrationCoeffs_[c];
         for (std::map<DetectorElementPtr, double>::iterator b = result.begin(); b
                 != result.end(); ++b) {
             DetectorElementPtr d = (*b).first;
             stream << *d << ": ";
             double ans = (*b).second;
             stream << ans << ", ";
         }
         stream << "]\n";
     }
 
     return stream;
 }

void pftools::SpaceManager::setBarrelLimit ( double limit )

inline

Definition at line 97 of file SpaceManager.h.

References barrelLimit_, and MessageLogger_cfi::limit.

                                       {
         barrelLimit_ = limit;
     }

Member Data Documentation

double pftools::SpaceManager::barrelLimit_

private

Definition at line 105 of file SpaceManager.h.

Referenced by createCalibrators(), evolveCoefficient(), and setBarrelLimit().

std::map<DetectorElementPtr, TF1> pftools::SpaceManager::barrelPosEvolutions_

private

Definition at line 118 of file SpaceManager.h.

Referenced by addEvolution(), and evolveCoefficient().

std::vector<SpaceVoxelPtr> pftools::SpaceManager::barrelPosRegion_

private

Definition at line 114 of file SpaceManager.h.

Referenced by createCalibrators(), extractEvolution(), and SpaceManager().

std::map<CalibratorPtr, std::map<DetectorElementPtr, double> > pftools::SpaceManager::calibrationCoeffs_

private

Definition at line 111 of file SpaceManager.h.

Referenced by assignCalibration(), extractEvolution(), getCalibration(), getNCalibrations(), interpolateCoefficient(), and printCalibrations().

double pftools::SpaceManager::endcapLimit_

private

Definition at line 107 of file SpaceManager.h.

Referenced by createCalibrators().

std::map<DetectorElementPtr, TF1> pftools::SpaceManager::endcapPosEvolutions_

private

Definition at line 119 of file SpaceManager.h.

Referenced by addEvolution(), and evolveCoefficient().

std::vector<SpaceVoxelPtr> pftools::SpaceManager::endcapPosRegion_

private

Definition at line 116 of file SpaceManager.h.

Referenced by createCalibrators(), extractEvolution(), and SpaceManager().

std::map<CalibratorPtr, SpaceVoxelPtr> pftools::SpaceManager::inverseAddressBook_

private

Definition at line 110 of file SpaceManager.h.

Referenced by interpolateCoefficient(), and makeInverseAddressBook().

std::map<SpaceVoxelPtr, CalibratorPtr> pftools::SpaceManager::myAddressBook

private

Definition at line 109 of file SpaceManager.h.

Referenced by clear(), createCalibrator(), createCalibrators(), extractEvolution(), findCalibrator(), getCalibrators(), makeInverseAddressBook(), and printCalibrations().

std::vector<SpaceVoxelPtr> pftools::SpaceManager::myKnownSpaceVoxels

private

Definition at line 112 of file SpaceManager.h.

Referenced by createCalibrator(), createCalibrators(), findCalibrator(), and printCalibrations().

std::string pftools::SpaceManager::name_

private

Definition at line 103 of file SpaceManager.h.

Referenced by getName().

std::map<Region, std::vector<SpaceVoxelPtr> > pftools::SpaceManager::regionsToSVs_

private

Definition at line 121 of file SpaceManager.h.

Referenced by SpaceManager().

double pftools::SpaceManager::transitionLimit_

private

Definition at line 106 of file SpaceManager.h.

std::vector<SpaceVoxelPtr> pftools::SpaceManager::transitionPosRegion_

private

Definition at line 115 of file SpaceManager.h.

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation