pftools::SpaceManager Class Reference

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

#include <SpaceManager.h>

Public Member Functions
void	addEvolution (DetectorElementPtr det, Region region, TF1 f)
void	assignCalibration (CalibratorPtr c, 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) throw (PFToolsException&)
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 ( DetectorElementPtr  det, Region  region, TF1  f  )

inline

Definition at line 79 of file SpaceManager.h.

References pftools::BARREL_POS, barrelPosEvolutions_, pftools::ENDCAP_POS, endcapPosEvolutions_, and f.

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

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

Definition at line 246 of file SpaceManager.cc.

References trackerHits::c, calibrationCoeffs_, makeInverseAddressBook(), and query::result.

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

void SpaceManager::clear

(

void

)

Definition at line 22 of file SpaceManager.cc.

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

Referenced by createCalibrators().

                         {
    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 trackerHits::c, pftools::Calibrator::clone(), prof2calltree::count, myAddressBook, myKnownSpaceVoxels, and alignCSCRings::s.

                         {
    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.

                                                                          {
    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
	)
throw	(	PFToolsException &
	)

Definition at line 168 of file SpaceManager.cc.

References trackerHits::c, hitfit::clear(), gather_cfg::cout, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, gen::k, prof2calltree::l, m, jptDQMConfig_cff::phiMax, and jptDQMConfig_cff::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_, trackerHits::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_, and endcapPosEvolutions_.

                                {
    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_, trackerHits::c, calibrationCoeffs_, pftools::DetElNames, pftools::ENDCAP_POS, endcapPosRegion_, i, myAddressBook, mergeVDriftHistosByStation::name, profilereader::profile, pftools::RegionNames, alignCSCRings::s, and python.multivaluedict::sort().

                       {

    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, prof2calltree::count, myAddressBook, myKnownSpaceVoxels, alignCSCRings::s, and edm::second().

Referenced by 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 trackerHits::c, and calibrationCoeffs_.

                                                                               {
    return calibrationCoeffs_[c];
}

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

inline

Definition at line 52 of file SpaceManager.h.

References myAddressBook.

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

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

inline

Definition at line 29 of file SpaceManager.h.

References name_.

                        {
        return name_;
    }

int pftools::SpaceManager::getNCalibrations ( )

inline

Definition at line 90 of file SpaceManager.h.

References calibrationCoeffs_.

                           {
        return calibrationCoeffs_.size();
    }

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

Definition at line 330 of file SpaceManager.cc.

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

                                               {
    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 trackerHits::c, inverseAddressBook_, myAddressBook, and alignCSCRings::s.

Referenced by assignCalibration().

                                          {
    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, trackerHits::c, calibrationCoeffs_, myAddressBook, myKnownSpaceVoxels, query::result, and alignCSCRings::s.

                                                              {
    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_cff::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.