#include <CovarianceParameterization.h>

Classes
struct	CompressionSchema

Public Member Functions
	CovarianceParameterization ()

bool	isValid () const

void	load (int version)

int	loadedVersion () const

float	meanValue (int i, int j, int sign, float pt, float eta, int nHits, int pixelHits, float cii=1., float cjj=1.) const

float	pack (float value, int schema, int i, int j, float pt, float eta, int nHits, int pixelHits, float cii=1., float cjj=1.) const

float	unpack (uint16_t packed, int schema, int i, int j, float pt, float eta, int nHits, int pixelHits, float cii=1., float cjj=1.) const

Static Public Member Functions
static int	index (int i, int j)

Private Member Functions
void	addTheHistogram (std::vector< TH3D * > *HistoVector, std::string StringToAddInTheName, int i, int j, TFile &fileToRead)

void	readFile (TFile &)

Private Attributes
std::vector< TH3D * >	cov_elements_noPixelHit

std::vector< TH3D * >	cov_elements_pixelHit

TFile *	fileToRead_

int	loadedVersion_

std::unordered_map< uint16_t, CompressionSchema >	schemas

Detailed Description

Definition at line 24 of file CovarianceParameterization.h.

Constructor & Destructor Documentation

CovarianceParameterization::CovarianceParameterization ( )

inline

Definition at line 38 of file CovarianceParameterization.h.

38 : loadedVersion_(-1) {}

CovarianceParameterization::loadedVersion_

int loadedVersion_

Definition: CovarianceParameterization.h:69

Member Function Documentation

void CovarianceParameterization::addTheHistogram	(	std::vector< TH3D * > *	HistoVector,
		std::string	StringToAddInTheName,
		int	i,
		int	j,
		TFile &	fileToRead
	)

private

Definition at line 137 of file CovarianceParameterization.cc.

References mps_fire::i, dqmiolumiharvest::j, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                      {
   std::string List_covName[5] = {"qoverp", "lambda", "phi", "dxy", "dsz"};
 
   std::string histoNameString = "covariance_" + List_covName[i] + "_" + List_covName[j] + StringToAddInTheName +
                                 "parametrization";  // + "_entries";
   TH3D *matrixElememtHistogramm = (TH3D *)fileToRead.Get(histoNameString.c_str());
   HistoVector->push_back(matrixElememtHistogramm);
 }

static int CovarianceParameterization::index	(	int	i,
		int	j
	)

inlinestatic

Definition at line 26 of file CovarianceParameterization.h.

Referenced by BeautifulSoup.PageElement::insert().

                                  {
     if (i >= j)
       return j + i * (i + 1) / 2;
     else
       return i + j * (j + 1) / 2;
   }

bool CovarianceParameterization::isValid ( void ) const

inline

Definition at line 39 of file CovarianceParameterization.h.

Referenced by ntupleDataFormat._Object::_checkIsValid(), core.AutoHandle.AutoHandle::ReallyLoad(), and pat::PackedCandidate::unpackCovariance().

39 { return loadedVersion_ != -1; }

CovarianceParameterization::loadedVersion_

int loadedVersion_

Definition: CovarianceParameterization.h:69

void CovarianceParameterization::load ( int version )

Definition at line 78 of file CovarianceParameterization.cc.

References GetRecoTauVFromDQM_MC_cff::cl, CompressionElement::CompressionElement(), printsummarytable::folder, dqm-mbProfile::format, edm::FileInPath::fullPath(), mps_fire::i, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, crabWrapper::key, GetRecoTauVFromDQM_MC_cff::next, AlCaHLTBitMon_ParallelJobs::p, dtResolutionTest_cfi::readFile, str, AlCaHLTBitMon_QueryRunRegistry::string, and BeamSplash_cfg::version.

                                                  {
   edm::FileInPath fip(
       (boost::format("DataFormats/PatCandidates/data/CovarianceParameterization_version%d.root") % version).str());
   fileToRead_ = TFile::Open(fip.fullPath().c_str());
   TFile &fileToRead = *fileToRead_;
   //Read files from here fip.fullPath().c_str();
   if (fileToRead.IsOpen()) {
     readFile(fileToRead);
 
     TIter next(((TDirectoryFile *)fileToRead.Get("schemas"))->GetListOfKeys());
     TKey *key;
     while ((key = (TKey *)next())) {
       TClass *cl = gROOT->GetClass(key->GetClassName());
       if (!cl->InheritsFrom("TDirectoryFile"))
         continue;
       std::string schemaNumber = key->ReadObj()->GetName();
       uint16_t schemaN = std::stoi(schemaNumber);
       //for (int folderNumber = 0; folderNumber < 6 ; folderNumber++) {
       CompressionSchema schema;
       for (int i = 0; i < 5; i++) {
         for (int j = i; j < 5; j++) {  //FILLING ONLY THE SCHEMA OF SOME ELEMENTS
           std::string folder = "schemas/" + schemaNumber + "/" + char(48 + i) + char(48 + j);
           std::string methodString = folder + "/method";
           std::string targetString = folder + "/target";
           std::string bitString = folder + "/bit";
           std::vector<float> vParams;
           TVector *p = (TVector *)fileToRead.Get((folder + "/param").c_str());
           for (int k = 0; k < p->GetNoElements(); k++) {
             vParams.push_back((*p)[k]);
           }
 
           schema(i, j) = CompressionElement(
               (CompressionElement::Method)((TParameter<int> *)fileToRead.Get(methodString.c_str()))->GetVal(),
               (CompressionElement::Target)((TParameter<int> *)fileToRead.Get(targetString.c_str()))->GetVal(),
               (int)((TParameter<int> *)fileToRead.Get(bitString.c_str()))->GetVal(),
               vParams);
         }
       }
       schemas[schemaN] = schema;
     }
 
     loadedVersion_ = version;
   } else {
     loadedVersion_ = -1;
   }
 }

int CovarianceParameterization::loadedVersion ( ) const

inline

Definition at line 40 of file CovarianceParameterization.h.

References PVValHelper::eta, mps_fire::i, dqmiolumiharvest::j, svgfig::load(), CompressionElement::pack(), DiDispStaMuonMonitor_cfi::pt, dtResolutionTest_cfi::readFile, Validation_hcalonly_cfi::sign, AlCaHLTBitMon_QueryRunRegistry::string, CompressionElement::unpack(), and BeamSplash_cfg::version.

40 { return loadedVersion_; }

CovarianceParameterization::loadedVersion_

int loadedVersion_

Definition: CovarianceParameterization.h:69

float CovarianceParameterization::meanValue	(	int	i,
		int	j,
		int	sign,
		float	pt,
		float	eta,
		int	nHits,
		int	pixelHits,
		float	cii = `1.`,
		float	cjj = `1.`
	)		const

Definition at line 147 of file CovarianceParameterization.cc.

References funct::abs(), muonRecoAnalyzer_cfi::etaBin, mps_fire::i, dqmiolumiharvest::j, and muonRecoAnalyzer_cfi::ptBin.

                                                                                                        {
   int hitNumberToUse = nHits;
   if (hitNumberToUse < 2)
     hitNumberToUse = 2;
   if (hitNumberToUse > 32)
     hitNumberToUse = 32;
   int ptBin = cov_elements_pixelHit[0]->GetXaxis()->FindBin(pt);
   int etaBin = cov_elements_pixelHit[0]->GetYaxis()->FindBin(std::abs(eta));
   int hitBin = cov_elements_pixelHit[0]->GetZaxis()->FindBin(hitNumberToUse);
   int min_idx = i;
   int max_idx = j;
 
   if (i > j) {
     min_idx = j;
     max_idx = i;
   }
 
   int indexOfTheHitogramInTheList = ((9 - min_idx) * min_idx) / 2 + max_idx;
 
   double meanValue = 0.;
   if (pixelHits > 0) {
     meanValue = sign * cov_elements_pixelHit[indexOfTheHitogramInTheList]->GetBinContent(ptBin, etaBin, hitBin);
   } else {
     meanValue = sign * cov_elements_noPixelHit[indexOfTheHitogramInTheList]->GetBinContent(ptBin, etaBin, hitBin);
   }
   return meanValue;
 }

float CovarianceParameterization::pack	(	float	value,
		int	schema,
		int	i,
		int	j,
		float	pt,
		float	eta,
		int	nHits,
		int	pixelHits,
		float	cii = `1.`,
		float	cjj = `1.`
	)		const

Definition at line 176 of file CovarianceParameterization.cc.

References DiDispStaMuonMonitor_cfi::pt, edm::second(), and std::swap().

Referenced by pat::PackedCandidate::packCovarianceElement().

                                                                                                                       {
   if (i > j)
     std::swap(i, j);
   float ref = meanValue(i, j, 1., pt, eta, nHits, pixelHits, cii, cjj);
   if (ref == 0) {
     schema = 0;
   }
   if (schema == 0 && i == j && (i == 2 || i == 0))
     ref = 1. / (pt * pt);
   /*  //Used for debugging, to be later removed  
     uint16_t p=(*schemas.find(schema)).second(i,j).pack(value,ref);
     float up=(*schemas.find(schema)).second(i,j).unpack(p,ref);
     std::cout << "check " << i << " " << j << " " << value << " " << up << " " << p << " " << ref << " " << schema<< std::endl;*/
   return (*schemas.find(schema)).second(i, j).pack(value, ref);
 }

void CovarianceParameterization::readFile ( TFile & f )

private

Definition at line 125 of file CovarianceParameterization.cc.

References mps_fire::i, dqmiolumiharvest::j, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by Vispa.Plugins.EventBrowser.EventBrowserTabController.EventBrowserTabController::navigate(), Vispa.Main.TabController.TabController::open(), and Vispa.Main.TabController.TabController::refresh().

                                                   {
   for (int i = 0; i < 5; i++) {
     for (int j = i; j < 5; j++) {
       std::string String_first_positive = "_pixel_";
       std::string String_second_positive = "_noPixel_";
 
       addTheHistogram(&cov_elements_pixelHit, String_first_positive, i, j, f);
       addTheHistogram(&cov_elements_noPixelHit, String_second_positive, i, j, f);
     }
   }
 }

float CovarianceParameterization::unpack	(	uint16_t	packed,
		int	schema,
		int	i,
		int	j,
		float	pt,
		float	eta,
		int	nHits,
		int	pixelHits,
		float	cii = `1.`,
		float	cjj = `1.`
	)		const

Definition at line 192 of file CovarianceParameterization.cc.

References MillePedeFileConverter_cfg::e, DiDispStaMuonMonitor_cfi::pt, edm::second(), and std::swap().

Referenced by pat::PackedCandidate::unpackCovarianceElement().

           {
   if (i > j)
     std::swap(i, j);
   float ref = meanValue(i, j, 1., pt, eta, nHits, pixelHits, cii, cjj);
   if (ref == 0) {
     schema = 0;
   }
   if (schema == 0 && i == j && (i == 2 || i == 0))
     ref = 1. / (pt * pt);
   if (i == j && (*schemas.find(schema)).second(i, j).unpack(packed, ref) == 0)
     return 1e-9;
   else
     return (*schemas.find(schema)).second(i, j).unpack(packed, ref);
 }