CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Classes | Functions | Variables
HOCalibAnalyzer.cc File Reference
#include <memory>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/HcalCalibObjects/interface/HOCalibVariables.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "TMath.h"
#include "TFile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TTree.h"
#include "TProfile.h"
#include "TPostScript.h"
#include "TCanvas.h"
#include "TF1.h"
#include "TStyle.h"
#include "TMinuit.h"
#include <string>
#include <iostream>
#include <fstream>
#include <iomanip>

Go to the source code of this file.

Classes

class  HOCalibAnalyzer
 

Functions

 DEFINE_FWK_MODULE (HOCalibAnalyzer)
 
void fcnbg (Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t flag)
 
void fcnsg (Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t flag)
 
Double_t gausX (Double_t *x, Double_t *par)
 
Double_t langaufun (Double_t *x, Double_t *par)
 
void set_mean (double &x, bool mdigi)
 
void set_sigma (double &x, bool mdigi)
 
Double_t totalfunc (Double_t *x, Double_t *par)
 

Variables

vector< float > cro_ssg [netamx][nphimx+1]
 
static const int etamap [4][21]
 
int ietafit
 
int iphifit
 
static const int mapx0m [9][2] ={{17,19}, {14,15}, {13,16}, {9,11}, {0,0}, {8,12}, {5,6}, {4,7}, {1,3}}
 
static const int mapx0p [9][2] ={{3,1}, {7,4}, {6,5}, {12,8}, {0,0}, {11,9}, {16,13}, {15,14}, {19,17}}
 
static const int mapx1 [6][3] ={{1,4,8}, {12,7,3}, {5,9,13}, {11,6,2}, {16,15,14}, {19,18,17}}
 
static const int mapx2 [6][3] ={{1,4,8}, {12,7,3}, {5,9,13}, {11,6,2}, {16,15,14}, {-1,-1,-1}}
 
const int nbgpr = 3
 
static const int netamx =30
 
static const int nphimx =72
 
static const int npixlebt [21] ={0, 0, 0, 0, 0, 1, 2, 3, 0, 4, 0, 6, 7, 8, 9, 0,11,13,14,15,0}
 
static const int npixleft [21] ={0, 0, 1, 2, 0, 4, 5, 6, 0, 8, 0, 0,11, 0,13,14,15, 0,17,18,0}
 
static const int npixleup [21] ={0, 4, 5, 6, 8, 9, 0,11, 0,13, 0,15,16, 0,17,18,19, 0, 0, 0,0}
 
static const int npixribt [21] ={0, 0, 0, 0, 1, 2, 3, 0, 4, 5, 0, 7, 0, 9, 0,11,12,14,15,16,0}
 
static const int npixrigh [21] ={0, 2, 3, 0, 5, 6, 7, 0, 9, 0, 0,12, 0,14,15,16, 0,18,19, 0,0}
 
static const int npixriup [21] ={0, 5, 6, 7, 9, 0,11,12,13,14, 0,16, 0,17,18,19, 0, 0, 0, 0,0}
 
const int nsgpr = 7
 
static const int phimap [4][21]
 
vector< float > sig_reg [netamx][nphimx+1]
 

Function Documentation

DEFINE_FWK_MODULE ( HOCalibAnalyzer  )
void fcnbg ( Int_t &  npar,
Double_t *  gin,
Double_t &  f,
Double_t *  par,
Int_t  flag 
)

Definition at line 174 of file HOCalibAnalyzer.cc.

References cro_ssg, i, ietafit, iphifit, funct::log(), and max().

Referenced by HOCalibAnalyzer::endJob().

174  {
175 
176  double fval = -par[0];
177  for (unsigned i=0; i<cro_ssg[ietafit][iphifit].size(); i++) {
178  double xval = (double)cro_ssg[ietafit][iphifit][i];
179  fval +=log(max(1.e-30,par[0]*TMath::Gaus(xval, par[1], par[2], 1)));
180  // fval +=log(par[0]*TMath::Gaus(xval, par[1], par[2], 1));
181  }
182  f = -fval;
183 }
i
int i
Definition: DBlmapReader.cc:9
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
iphifit
int iphifit
Definition: HOCalibAnalyzer.cc:103
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
funct::log
Log< T >::type log(const T &t)
Definition: Log.h:22
cro_ssg
vector< float > cro_ssg[netamx][nphimx+1]
Definition: HOCalibAnalyzer.cc:105
ietafit
int ietafit
Definition: HOCalibAnalyzer.cc:102
Gflash::par
const double par[8 *NPar][4]
Definition: GflashNameSpace.h:177
void fcnsg ( Int_t &  npar,
Double_t *  gin,
Double_t &  f,
Double_t *  par,
Int_t  flag 
)

Definition at line 185 of file HOCalibAnalyzer.cc.

References i, ietafit, iphifit, funct::log(), sig_reg, and totalfunc().

Referenced by HOCalibAnalyzer::endJob().

185  {
186 
187  double xval[2];
188  double fval = -(par[0]+par[5]);
189  for (unsigned i=0; i<sig_reg[ietafit][iphifit].size(); i++) {
190  xval[0] = (double)sig_reg[ietafit][iphifit][i];
191  fval +=log(totalfunc(xval, par));
192  }
193  f = -fval;
194 }
i
int i
Definition: DBlmapReader.cc:9
totalfunc
Double_t totalfunc(Double_t *x, Double_t *par)
Definition: HOCalibAnalyzer.cc:170
sig_reg
vector< float > sig_reg[netamx][nphimx+1]
Definition: HOCalibAnalyzer.cc:104
iphifit
int iphifit
Definition: HOCalibAnalyzer.cc:103
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
funct::log
Log< T >::type log(const T &t)
Definition: Log.h:22
ietafit
int ietafit
Definition: HOCalibAnalyzer.cc:102
Gflash::par
const double par[8 *NPar][4]
Definition: GflashNameSpace.h:177
Double_t gausX ( Double_t *  x,
Double_t *  par 
)

Definition at line 113 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::endJob(), and totalfunc().

113  {
114  return par[0]*(TMath::Gaus(x[0], par[1], par[2], kTRUE));
115 }
x
Definition: DDAxes.h:10
Gflash::par
const double par[8 *NPar][4]
Definition: GflashNameSpace.h:177
Double_t langaufun ( Double_t *  x,
Double_t *  par 
)

Definition at line 117 of file HOCalibAnalyzer.cc.

References i, runTheMatrix::np, and launcher::step.

Referenced by fitUtilities::doLanGaussFit(), HOCalibAnalyzer::endJob(), langaupro(), and totalfunc().

117  {
118 
119  //Fit parameters:
120  //par[0]*par[1]=Width (scale) parameter of Landau density
121  //par[1]=Most Probable (MP, location) parameter of Landau density
122  //par[2]=Total area (integral -inf to inf, normalization constant)
123  //par[3]=Width (sigma) of convoluted Gaussian function
124  //
125  //In the Landau distribution (represented by the CERNLIB approximation),
126  //the maximum is located at x=-0.22278298 with the location parameter=0.
127  //This shift is corrected within this function, so that the actual
128  //maximum is identical to the MP parameter.
129  // /*
130  // Numeric constants
131  Double_t invsq2pi = 0.3989422804014; // (2 pi)^(-1/2)
132  Double_t mpshift = -0.22278298; // Landau maximum location
133 
134  // Control constants
135  Double_t np = 100.0; // number of convolution steps
136  Double_t sc = 5.0; // convolution extends to +-sc Gaussian sigmas
137 
138  // Variables
139  Double_t xx;
140  Double_t mpc;
141  Double_t fland;
142  Double_t sum = 0.0;
143  Double_t xlow,xupp;
144  Double_t step;
145  Double_t i;
146 
147 
148  // MP shift correction
149  mpc = par[1] - mpshift * par[0]*par[1];
150 
151  // Range of convolution integral
152  xlow = x[0] - sc * par[3];
153  xupp = x[0] + sc * par[3];
154 
155  step = (xupp-xlow) / np;
156 
157  // Convolution integral of Landau and Gaussian by sum
158  for(i=1.0; i<=np/2; i++) {
159  xx = xlow + (i-.5) * step;
160  fland = TMath::Landau(xx,mpc,par[0]*par[1], kTRUE); // / par[0];
161  sum += fland * TMath::Gaus(x[0],xx,par[3]);
162  xx = xupp - (i-.5) * step;
163  fland = TMath::Landau(xx,mpc,par[0]*par[1], kTRUE); // / par[0];
164  sum += fland * TMath::Gaus(x[0],xx,par[3]);
165  }
166 
167  return (par[2] * step * sum * invsq2pi / par[3]);
168 }
i
int i
Definition: DBlmapReader.cc:9
launcher.step
list step
Definition: launcher.py:15
x
Definition: DDAxes.h:10
Gflash::par
const double par[8 *NPar][4]
Definition: GflashNameSpace.h:177
void set_mean ( double &  x,
bool  mdigi 
)

Definition at line 196 of file HOCalibAnalyzer.cc.

References max(), and min.

Referenced by HOCalibAnalyzer::endJob().

196  {
197  if(mdigi) {
198  x = min(x, 0.5);
199  x = max(x, -0.5);
200  } else {
201  x = min(x, 0.1);
202  x = max(x, -0.1);
203  }
204 }
min
#define min(a, b)
Definition: mlp_lapack.h:161
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
x
Definition: DDAxes.h:10
void set_sigma ( double &  x,
bool  mdigi 
)

Definition at line 206 of file HOCalibAnalyzer.cc.

References max(), and min.

Referenced by HOCalibAnalyzer::endJob().

206  {
207  if(mdigi) {
208  x = min(x, 1.2);
209  x = max(x, -1.2);
210  } else {
211  x = min(x, 0.24);
212  x = max(x, -0.24);
213  }
214 }
min
#define min(a, b)
Definition: mlp_lapack.h:161
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
x
Definition: DDAxes.h:10
Double_t totalfunc ( Double_t *  x,
Double_t *  par 
)

Definition at line 170 of file HOCalibAnalyzer.cc.

References gausX(), and langaufun().

Referenced by HOCalibAnalyzer::endJob(), and fcnsg().

170  {
171  return gausX(x, par) + langaufun(x, &par[3]);
172 }
gausX
Double_t gausX(Double_t *x, Double_t *par)
Definition: HOCalibAnalyzer.cc:113
x
Definition: DDAxes.h:10
langaufun
Double_t langaufun(Double_t *x, Double_t *par)
Definition: HOCalibAnalyzer.cc:117
Gflash::par
const double par[8 *NPar][4]
Definition: GflashNameSpace.h:177

Variable Documentation

vector<float> cro_ssg[netamx][nphimx+1]

Definition at line 105 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze(), and fcnbg().

const int etamap[4][21]
static
Initial value:
={{-1, 0,3,1, 0,2,3, 1,0,2, -1, 3,1,2, 4,4,4, -1,-1,-1, -1},
{-1, 0,3,1, 0,2,3, 1,0,2, -1, 3,1,2, 4,4,4, 5,5,5, -1},
{-1, 0,-1,0, 1,2,2, 1,3,5, -1, 5,3,6, 7,7,6, 8,-1,8, -1},
{-1, 8,-1,8, 7,6,6, 7,5,3, -1, 3,5,2, 1,1,2, 0,-1,0, -1}}

Definition at line 79 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

int ietafit

Definition at line 102 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::endJob(), fcnbg(), and fcnsg().

int iphifit

Definition at line 103 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::endJob(), fcnbg(), and fcnsg().

const int mapx0m[9][2] ={{17,19}, {14,15}, {13,16}, {9,11}, {0,0}, {8,12}, {5,6}, {4,7}, {1,3}}
static

Definition at line 77 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int mapx0p[9][2] ={{3,1}, {7,4}, {6,5}, {12,8}, {0,0}, {11,9}, {16,13}, {15,14}, {19,17}}
static

Definition at line 76 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int mapx1[6][3] ={{1,4,8}, {12,7,3}, {5,9,13}, {11,6,2}, {16,15,14}, {19,18,17}}
static

Definition at line 72 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int mapx2[6][3] ={{1,4,8}, {12,7,3}, {5,9,13}, {11,6,2}, {16,15,14}, {-1,-1,-1}}
static

Definition at line 74 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int nbgpr = 3

Definition at line 99 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::endJob().

const int netamx =30
static

Definition at line 97 of file HOCalibAnalyzer.cc.

const int nphimx =72
static

Definition at line 98 of file HOCalibAnalyzer.cc.

const int npixlebt[21] ={0, 0, 0, 0, 0, 1, 2, 3, 0, 4, 0, 6, 7, 8, 9, 0,11,13,14,15,0}
static

Definition at line 92 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int npixleft[21] ={0, 0, 1, 2, 0, 4, 5, 6, 0, 8, 0, 0,11, 0,13,14,15, 0,17,18,0}
static

Definition at line 90 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int npixleup[21] ={0, 4, 5, 6, 8, 9, 0,11, 0,13, 0,15,16, 0,17,18,19, 0, 0, 0,0}
static

Definition at line 94 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int npixribt[21] ={0, 0, 0, 0, 1, 2, 3, 0, 4, 5, 0, 7, 0, 9, 0,11,12,14,15,16,0}
static

Definition at line 93 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int npixrigh[21] ={0, 2, 3, 0, 5, 6, 7, 0, 9, 0, 0,12, 0,14,15,16, 0,18,19, 0,0}
static

Definition at line 91 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int npixriup[21] ={0, 5, 6, 7, 9, 0,11,12,13,14, 0,16, 0,17,18,19, 0, 0, 0, 0,0}
static

Definition at line 95 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

const int nsgpr = 7

Definition at line 100 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::endJob().

const int phimap[4][21]
static
Initial value:
={{-1, 0,2,2, 1,0,1, 1,2,1, -1, 0,0,2, 2,1,0, 2,1,0, -1},
{-1, 0,2,2, 1,0,1, 1,2,1, -1, 0,0,2, 2,1,0, 2,1,0, -1},
{-1, 1,-1,0, 1,1,0, 0,1,1, -1, 0,0,1, 1,0,0, 1,-1,0, -1},
{-1, 0,-1,1, 0,0,1, 1,0,0, -1, 1,1,0, 0,1,1, 0,-1,1, -1}}

Definition at line 84 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze().

vector<float> sig_reg[netamx][nphimx+1]

Definition at line 104 of file HOCalibAnalyzer.cc.

Referenced by HOCalibAnalyzer::analyze(), HOCalibAnalyzer::endJob(), and fcnsg().

Generated for CMSSW Reference Manual by doxygen 1.8.5