Resolution Map (resolution as a function of eta and E) More...

#include <PFResolutionMap.h>

Inheritance diagram for PFResolutionMap:

Public Member Functions
int	FindBin (double eta, double e, double z=0) override
	extrapolation requires overloading of this function More...

const char *	GetMapFile () const

double	getRes (double eta, double phi, double e, int MapEta=-1)

	PFResolutionMap ()
	default constructor More...

	PFResolutionMap (const char name, const char mapfile)
	create a map from text file mapfile More...

	PFResolutionMap (const char *name, unsigned nbinseta, double mineta, double maxeta, unsigned nbinse, double mine, double maxe, double value=-1)
	create an empty map and initialize it More...

	PFResolutionMap (const TH2D &h)
	create a map from a 2d histogram More...

bool	ReadMapFile (const char *mapfile)
	read text file More...

bool	WriteMapFile (const char *mapfile)

Private Member Functions
double	dCrackPhi (double phi, double eta)

bool	IsInAPhiCrack (double phi, double eta)

double	minimum (double a, double b)

Private Attributes
std::string	mapFile_

Static Private Attributes
static const unsigned	lineSize_ = 10000

Friends
std::ostream &	operator<< (std::ostream &out, const PFResolutionMap &rm)
	print this map More...

Detailed Description

Resolution Map (resolution as a function of eta and E)

Basically just a TH2D with text I/O

Author: Colin Bernet

Date: January 2006

Definition at line 16 of file PFResolutionMap.h.

Constructor & Destructor Documentation

◆ PFResolutionMap() [1/4]

PFResolutionMap::PFResolutionMap ( )

inline

default constructor

Definition at line 19 of file PFResolutionMap.h.

19 : TH2D() {}

◆ PFResolutionMap() [2/4]

PFResolutionMap::PFResolutionMap	(	const char *	name,
		const char *	mapfile
	)

create a map from text file mapfile

Definition at line 38 of file PFResolutionMap.cc.

                                                                       {
   SetTitle(name);
   GetXaxis()->SetTitle("#eta");
   GetYaxis()->SetTitle("E");
   if (!ReadMapFile(mapfile)) {
     string err = "PFResolutionMap::PFResolutionMap : cannot read file ";
     err += mapfile;
     throw invalid_argument(err);
   }
 }

References submitPVResolutionJobs::err, Skims_PA_cff::name, and ReadMapFile().

◆ PFResolutionMap() [3/4]

PFResolutionMap::PFResolutionMap	(	const char *	name,
		unsigned	nbinseta,
		double	mineta,
		double	maxeta,
		unsigned	nbinse,
		double	mine,
		double	maxe,
		double	value = `-1`
	)

create an empty map and initialize it

Definition at line 15 of file PFResolutionMap.cc.

     : TH2D(name, name, nbinseta, mineta, maxeta, nbinse, mine, maxe) {
   GetXaxis()->SetTitle("#eta");
   GetYaxis()->SetTitle("E");
  
   if (value > 0) {
     for (int ie = 1; ie <= GetNbinsY(); ie++) {
       for (int ieta = 1; ieta <= GetNbinsX(); ieta++) {
         SetBinContent(ieta, ie, value);
       }
     }
   }
   // cout<<"creating resolution map "<<endl;
   // Print("all");
 }

References LEDCalibrationChannels::ieta.

◆ PFResolutionMap() [4/4]

PFResolutionMap::PFResolutionMap ( const TH2D & h )

inline

create a map from a 2d histogram

Definition at line 35 of file PFResolutionMap.h.

35 : TH2D(h) {}

Member Function Documentation

◆ dCrackPhi()

double PFResolutionMap::dCrackPhi	(	double	phi,
		double	eta
	)

private

Definition at line 216 of file PFResolutionMap.cc.

                                                         {
   constexpr double pi = M_PI;  // 3.14159265358979323846;
   constexpr double twopi = 2. * pi;
   constexpr double twopiO18 = pi / 9;
  
   //Location of the 18 phi-cracks
   constexpr double c0 = 2.97025;
   constexpr std::array<double, 18> cPhi{{c0,
                                          c0 - twopiO18,
                                          c0 - 2 * twopiO18,
                                          c0 - 3 * twopiO18,
                                          c0 - 4 * twopiO18,
                                          c0 - 5 * twopiO18,
                                          c0 - 6 * twopiO18,
                                          c0 - 7 * twopiO18,
                                          c0 - 8 * twopiO18,
                                          c0 - 9 * twopiO18,
                                          c0 - 10 * twopiO18,
                                          c0 - 11 * twopiO18,
                                          c0 - 12 * twopiO18,
                                          c0 - 13 * twopiO18,
                                          c0 - 14 * twopiO18,
                                          c0 - 15 * twopiO18,
                                          c0 - 16 * twopiO18,
                                          c0 - 17 * twopiO18}};
  
   //Shift of this location if eta<0
   constexpr double delta_cPhi = 0.00638;
  
   double m;  //the result
  
   //the location is shifted
   if (eta < 0) {
     phi += delta_cPhi;
     if (phi > pi)
       phi -= twopi;
   }
   if (phi >= -pi && phi <= pi) {
     //the problem of the extrema
     if (phi < cPhi[17] || phi >= cPhi[0]) {
       if (phi < 0)
         phi += twopi;
       m = minimum(phi - cPhi[0], phi - cPhi[17] - twopi);
     }
  
     //between these extrema...
     else {
       bool OK = false;
       unsigned i = 16;
       while (!OK) {
         if (phi < cPhi[i]) {
           m = minimum(phi - cPhi[i + 1], phi - cPhi[i]);
           OK = true;
         } else
           i -= 1;
       }
     }
   } else {
     m = 0.;  //if there is a problem, we assum that we are in a crack
     edm::LogWarning("PFResolutionMap:Problem") << "Problem in dminphi";
   }
   if (eta < 0)
     m = -m;  //because of the disymetry
   return m;
 }

References fftjetpileupestimator_calo_uncalib_cfi::c0, PVValHelper::eta, mps_fire::i, visualization-live-secondInstance_cfg::m, M_PI, minimum(), phi, and pi.

Referenced by IsInAPhiCrack().

◆ FindBin()

int PFResolutionMap::FindBin	(	double	eta,
		double	e,
		double	z = `0`
	)

override

extrapolation requires overloading of this function

Definition at line 168 of file PFResolutionMap.cc.

                                                            {
   if (e >= GetYaxis()->GetXmax())
     e = GetYaxis()->GetXmax() - 0.001;
  
   return TH2D::FindBin(eta, e);
 }

References MillePedeFileConverter_cfg::e, and PVValHelper::eta.

Referenced by getRes().

◆ GetMapFile()

const char* PFResolutionMap::GetMapFile ( ) const

inline

Definition at line 49 of file PFResolutionMap.h.

49 { return mapFile_.c_str(); }

References mapFile_.

◆ getRes()

double PFResolutionMap::getRes	(	double	eta,
		double	phi,
		double	e,
		int	MapEta = `-1`
	)

Definition at line 138 of file PFResolutionMap.cc.

                                                                            {
   constexpr double fMinEta = -2.95;
   constexpr double fMaxEta = 2.95;
   constexpr double fMinE = 0;
   constexpr double fMaxE = 100;
  
   if (eta < fMinEta)
     eta = fMinEta + 0.001;
   if (eta > fMaxEta)
     eta = fMaxEta - 0.001;
  
   if (e < fMinE)
     e = fMinE + 0.001;
   if (e > fMaxE)
     e = fMaxE - 0.001;
  
   unsigned bin = FindBin(TMath::Abs(eta), e);
  
   double res = GetBinContent(bin);
   if (MapEta > -1) {
     if ((eta < 1.48) && IsInAPhiCrack(phi, eta)) {
       if (MapEta == 1)
         res *= 1.88;
       else
         res *= 1.65;
     }
   }
   return res;
 }

References Abs(), newFWLiteAna::bin, MillePedeFileConverter_cfg::e, PVValHelper::eta, FindBin(), IsInAPhiCrack(), and phi.

◆ IsInAPhiCrack()

bool PFResolutionMap::IsInAPhiCrack	(	double	phi,
		double	eta
	)

private

Definition at line 198 of file PFResolutionMap.cc.

                                                           {
   double dminPhi = dCrackPhi(phi, eta);
   bool Is = (TMath::Abs(dminPhi) < 0.005);
   return Is;
 }

References Abs(), dCrackPhi(), PVValHelper::eta, and phi.

Referenced by getRes().

◆ minimum()

double PFResolutionMap::minimum	(	double	a,
		double	b
	)

private

Definition at line 205 of file PFResolutionMap.cc.

                                                   {
   if (TMath::Abs(b) < TMath::Abs(a))
     a = b;
   return a;
 }

References a, Abs(), and b.

Referenced by dCrackPhi().

◆ ReadMapFile()

bool PFResolutionMap::ReadMapFile ( const char * mapfile )

read text file

Definition at line 70 of file PFResolutionMap.cc.

                                                      {
   // open the file
  
   std::ifstream inf(mapfile);
   if (!inf.good()) {
     // cout<<"PFResolutionMap::Read : cannot open file "<<mapfile<<endl;
     return false;
   }
  
   // first data describes the map
  
   int nbinseta = 0;
   double mineta = 0;
   double maxeta = 0;
  
   int nbinse = 0;
   double mine = 0;
   double maxe = 0;
  
   inf >> nbinseta;
   inf >> mineta;
   inf >> maxeta;
  
   inf >> nbinse;
   inf >> mine;
   inf >> maxe;
  
   SetBins(nbinseta, mineta, maxeta, nbinse, mine, maxe);
  
   char s[lineSize_];
   int pos = inf.tellg();
  
   // parse map data
   int i = 1;
   do {
     inf.seekg(pos);
     inf.getline(s, lineSize_);
  
     pos = inf.tellg();
  
     if (string(s).empty()) {
       continue;  // remove empty lines
     }
  
     istringstream lin(s);
  
     double dataw;
     int j = 1;
     do {
       lin >> dataw;
       SetBinContent(j, i, dataw);
       j++;
     } while (lin.good());
     i++;
   } while (inf.good());
  
   if (inf.eof()) {
     mapFile_ = mapfile;
     return true;
   } else {
     // cout<<"error"<<endl;
     return false;
   }
  
   mapFile_ = mapfile;
   return true;
 }

References relativeConstraints::empty, mps_fire::i, dqmiodatasetharvest::inf, dqmiolumiharvest::j, lineSize_, mapFile_, and alignCSCRings::s.

Referenced by PFResolutionMap().

◆ WriteMapFile()

bool PFResolutionMap::WriteMapFile ( const char * mapfile )

write text file is not const because mapFile_ will be updated

Definition at line 49 of file PFResolutionMap.cc.

                                                       {
   //  assert(fData.size() == fNbinsEta*fNbinsE);
  
   // open the file
  
   std::ofstream outf(mapfile);
   if (!outf.good()) {
     edm::LogWarning("PFResolutionMap::Write") << " : cannot open file " << mapfile;
     return false;
   }
  
   outf << (*this) << endl;
   if (!outf.good()) {
     edm::LogError("PFResolutionMap::Write") << " : corrupted file " << mapfile;
     return false;
   } else {
     mapFile_ = mapfile;
     return true;
   }
 }

References mapFile_.

Friends And Related Function Documentation

◆ operator<<

std::ostream& operator<<	(	std::ostream &	out,
		const PFResolutionMap &	rm
	)

friend

print this map

Member Data Documentation

◆ lineSize_

const unsigned PFResolutionMap::lineSize_ = 10000

staticprivate

Definition at line 58 of file PFResolutionMap.h.

Referenced by ReadMapFile().

◆ mapFile_

std::string PFResolutionMap::mapFile_

private

Definition at line 59 of file PFResolutionMap.h.

Referenced by GetMapFile(), ReadMapFile(), and WriteMapFile().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ PFResolutionMap() [1/4]

◆ PFResolutionMap() [2/4]

◆ PFResolutionMap() [3/4]

◆ PFResolutionMap() [4/4]

Member Function Documentation

◆ dCrackPhi()

◆ FindBin()

◆ GetMapFile()

◆ getRes()

◆ IsInAPhiCrack()

◆ minimum()

◆ ReadMapFile()

◆ WriteMapFile()

Friends And Related Function Documentation

◆ operator<<

Member Data Documentation

◆ lineSize_

◆ mapFile_