pftools::SpaceManager Class Reference

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 != nullptr);
    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, PostProcessor_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 != nullptr) {
            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 != nullptr) {
            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.