CMS 3D CMS Logo

Functions | Variables
MuScleFitUtils.cc File Reference
#include "MuScleFitUtils.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "SimDataFormats/Track/interface/SimTrack.h"
#include "DataFormats/Candidate/interface/LeafCandidate.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "TString.h"
#include "TFile.h"
#include "TTree.h"
#include "TCanvas.h"
#include "TH2F.h"
#include "TF1.h"
#include "TF2.h"
#include <iostream>
#include <fstream>
#include <memory>
#include "MuonAnalysis/MomentumScaleCalibration/interface/Functions.h"

Go to the source code of this file.

Functions

Double_t Gaussian (Double_t *x, Double_t *par)
 
void likelihood (int &npar, double *grad, double &fval, double *xval, int flag)
 
Double_t lorentzianPeak (Double_t *x, Double_t *par)
 

Variables

double f [11][100]
 
double g [11][100]
 
TF1 * GL
 
TF2 * GL2
 
double isum
 
double mzsum
 

Function Documentation

Double_t Gaussian ( Double_t *  x,
Double_t *  par 
)

Definition at line 63 of file MuScleFitUtils.cc.

References JetChargeProducer_cfi::exp.

Referenced by MuScleFitUtils::fitReso().

63  {
64  return par[0]*exp(-0.5*((x[0]-par[1])/par[2])*((x[0]-par[1])/par[2]));
65 }
T x() const
Cartesian x coordinate.
void likelihood ( int &  npar,
double *  grad,
double &  fval,
double *  xval,
int  flag 
)

Definition at line 1709 of file MuScleFitUtils.cc.

References MuScleFitUtils::applyScale(), MuScleFitUtils::computeWeight(), gather_cfg::cout, MuScleFitUtils::debug, MuScleFitUtils::doScaleFit, reco::tau::disc::Eta(), MuScleFitUtils::iev_, MuScleFitUtils::invDimuonMass(), MuScleFitUtils::likelihoodInLoop_, cmsBatch::log, MuScleFitUtils::loopCounter, ResonanceBuilder::mass, MuScleFitUtils::massProb(), MuScleFitUtils::massResolution(), MuScleFitUtils::minuitLoop_, MuScleFitUtils::normalizationChanged_, MuScleFitUtils::normalizeLikelihoodByEventNumber_, MuScleFitUtils::oldNormalization_, TtFullHadEvtBuilder_cfi::prob, MuScleFitUtils::ReducedSavedPair, MuScleFitUtils::rminPtr_, MuScleFitUtils::SavedPair, and mps_merge::weight.

Referenced by pat::ElectronSelector::filter(), MuScleFitUtils::minimizeLikelihood(), MuScleFitUtils::byPt::operator()(), and PFNuclearProducer::produce().

1709  {
1710 
1711  if (MuScleFitUtils::debug>19) std::cout << "[MuScleFitUtils-likelihood]: In likelihood function" << std::endl;
1712 
1713  const lorentzVector * recMu1;
1714  const lorentzVector * recMu2;
1715  lorentzVector corrMu1;
1716  lorentzVector corrMu2;
1717 
1718  // if (MuScleFitUtils::debug>19) {
1719  // int parnumber = (int)(MuScleFitUtils::parResol.size()+MuScleFitUtils::parScale.size()+
1720  // MuScleFitUtils::parCrossSection.size()+MuScleFitUtils::parBgr.size());
1721  // std::cout << "[MuScleFitUtils-likelihood]: Looping on tree with ";
1722  // for (int ipar=0; ipar<parnumber; ipar++) {
1723  // std::cout << "Parameter #" << ipar << " with value " << xval[ipar] << " ";
1724  // }
1725  // std::cout << std::endl;
1726  // }
1727 
1728  // Loop on the tree
1729  // ----------------
1730  double flike = 0;
1731  int evtsinlik = 0;
1732  int evtsoutlik = 0;
1733  // std::cout << "SavedPair.size() = " << MuScleFitUtils::SavedPair.size() << std::endl;
1734  if( MuScleFitUtils::debug>0 ) {
1735  std::cout << "SavedPair.size() = " << MuScleFitUtils::SavedPair.size() << std::endl;
1736  std::cout << "ReducedSavedPair.size() = " << MuScleFitUtils::ReducedSavedPair.size() << std::endl;
1737  }
1738  // for( unsigned int nev=0; nev<MuScleFitUtils::SavedPair.size(); ++nev ) {
1739  for( unsigned int nev=0; nev<MuScleFitUtils::ReducedSavedPair.size(); ++nev ) {
1740 
1741  // recMu1 = &(MuScleFitUtils::SavedPair[nev].first);
1742  // recMu2 = &(MuScleFitUtils::SavedPair[nev].second);
1743  recMu1 = &(MuScleFitUtils::ReducedSavedPair[nev].first);
1744  recMu2 = &(MuScleFitUtils::ReducedSavedPair[nev].second);
1745 
1746  // Compute original mass
1747  // ---------------------
1748  double mass = MuScleFitUtils::invDimuonMass( *recMu1, *recMu2 );
1749 
1750  // Compute weight and reference mass (from original mass)
1751  // ------------------------------------------------------
1753  if( weight!=0. ) {
1754  // Compute corrected mass (from previous biases) only if we are currently fitting the scale
1755  // ----------------------------------------------------------------------------------------
1757 // std::cout << "Original pt1 = " << corrMu1.Pt() << std::endl;
1758 // std::cout << "Original pt2 = " << corrMu2.Pt() << std::endl;
1759  corrMu1 = MuScleFitUtils::applyScale(*recMu1, xval, -1);
1760  corrMu2 = MuScleFitUtils::applyScale(*recMu2, xval, 1);
1761 
1762 // if( (corrMu1.Pt() != corrMu1.Pt()) || (corrMu2.Pt() != corrMu2.Pt()) ) {
1763 // std::cout << "Rescaled pt1 = " << corrMu1.Pt() << std::endl;
1764 // std::cout << "Rescaled pt2 = " << corrMu2.Pt() << std::endl;
1765 // }
1766 // std::cout << "Rescaled pt1 = " << corrMu1.Pt() << std::endl;
1767 // std::cout << "Rescaled pt2 = " << corrMu2.Pt() << std::endl;
1768  }
1769  else {
1770  corrMu1 = *recMu1;
1771  corrMu2 = *recMu2;
1772 
1773 // if( (corrMu1.Pt() != corrMu1.Pt()) || (corrMu2.Pt() != corrMu2.Pt()) ) {
1774 // std::cout << "Not rescaled pt1 = " << corrMu1.Pt() << std::endl;
1775 // std::cout << "Not rescaled pt2 = " << corrMu2.Pt() << std::endl;
1776 // }
1777  }
1778  double corrMass = MuScleFitUtils::invDimuonMass(corrMu1, corrMu2);
1779  double Y = (corrMu1+corrMu2).Rapidity();
1780  double resEta = (corrMu1+corrMu2).Eta();
1781  if( MuScleFitUtils::debug>19 ) {
1782  std::cout << "[MuScleFitUtils-likelihood]: Original/Corrected resonance mass = " << mass
1783  << " / " << corrMass << std::endl;
1784  }
1785 
1786  // Compute mass resolution
1787  // -----------------------
1788  double massResol = MuScleFitUtils::massResolution(corrMu1, corrMu2, xval);
1789  if (MuScleFitUtils::debug>19)
1790  std::cout << "[MuScleFitUtils-likelihood]: Resolution is " << massResol << std::endl;
1791 
1792  // Compute probability of this mass value including background modeling
1793  // --------------------------------------------------------------------
1794  if (MuScleFitUtils::debug>1) std::cout << "calling massProb inside likelihood function" << std::endl;
1795 
1796  // double prob = MuScleFitUtils::massProb( corrMass, resEta, Y, massResol, xval );
1797  double prob = MuScleFitUtils::massProb( corrMass, resEta, Y, massResol, xval, false, corrMu1.eta(), corrMu2.eta() );
1798  if (MuScleFitUtils::debug>1) std::cout << "likelihood:massProb = " << prob << std::endl;
1799 
1800  // Compute likelihood
1801  // ------------------
1802  if( prob>0 ) {
1803  // flike += log(prob*10000)*weight; // NNBB! x10000 to see if we can recover the problem of boundary
1804  flike += log(prob)*weight;
1805  evtsinlik += 1; // NNBB test: see if likelihood per event is smarter (boundary problem)
1806  } else {
1807  if( MuScleFitUtils::debug > 0 ) {
1808  std::cout << "WARNING: corrMass = " << corrMass << " outside window, this will cause a discontinuity in the likelihood. Consider increasing the safety bands which are now set to 90% of the normalization window to avoid this problem" << std::endl;
1809  std::cout << "Original mass was = " << mass << std::endl;
1810  std::cout << "WARNING: massResol = " << massResol << " outside window" << std::endl;
1811  }
1812  evtsoutlik += 1;
1813  }
1814  if (MuScleFitUtils::debug>19)
1815  std::cout << "[MuScleFitUtils-likelihood]: Mass probability = " << prob << std::endl;
1816  } // weight!=0
1817 
1818  } // End of loop on tree events
1819 
1820 // // Protection for low statistic. If the likelihood manages to throw out all the signal
1821 // // events and stays with ~ 10 events in the resonance window it could have a better likelihood
1822 // // because of ~ uniformly distributed events (a random combination could be good and spoil the fit).
1823 // // We require that the number of events included in the fit does not change more than 5% in each minuit loop.
1824 // bool lowStatPenalty = false;
1825 // if( MuScleFitUtils::minuitLoop_ > 0 ) {
1826 // double newEventsOutInRatio = double(evtsinlik);
1827 // // double newEventsOutInRatio = double(evtsoutlik)/double(evtsinlik);
1828 // double ratio = newEventsOutInRatio/MuScleFitUtils::oldEventsOutInRatio_;
1829 // MuScleFitUtils::oldEventsOutInRatio_ = newEventsOutInRatio;
1830 // if( ratio < 0.8 || ratio > 1.2 ) {
1831 // std::cout << "Warning: too much change from oldEventsInLikelihood to newEventsInLikelihood, ratio is = " << ratio << std::endl;
1832 // std::cout << "oldEventsInLikelihood = " << MuScleFitUtils::oldEventsOutInRatio_ << ", newEventsInLikelihood = " << newEventsOutInRatio << std::endl;
1833 // lowStatPenalty = true;
1834 // }
1835 // }
1836 
1837  // It is a product of probabilities, we compare the sqrt_N of them. Thus N becomes a denominator of the logarithm.
1838  if( evtsinlik != 0 ) {
1839 
1841  // && !(MuScleFitUtils::duringMinos_) ) {
1842  if( MuScleFitUtils::rminPtr_ == 0 ) {
1843  std::cout << "ERROR: rminPtr_ = " << MuScleFitUtils::rminPtr_ << ", code will crash" << std::endl;
1844  }
1845  double normalizationArg[] = {1/double(evtsinlik)};
1846  // Reset the normalizationArg only if it changed
1847  if( MuScleFitUtils::oldNormalization_ != normalizationArg[0] ) {
1848  int ierror = 0;
1849 // if( MuScleFitUtils::likelihoodInLoop_ != 0 ) {
1850 // // This condition is set only when minimizing. Later calls of hesse and minos will not change the value
1851 // // This is done to avoid minos being confused by changing the UP parameter during its computation.
1852 // MuScleFitUtils::rminPtr_->mnexcm("SET ERR", normalizationArg, 1, ierror);
1853 // }
1854  MuScleFitUtils::rminPtr_->mnexcm("SET ERR", normalizationArg, 1, ierror);
1855  std::cout << "oldNormalization = " << MuScleFitUtils::oldNormalization_ << " new = " << normalizationArg[0] << std::endl;
1856  MuScleFitUtils::oldNormalization_ = normalizationArg[0];
1858  }
1859  fval = -2.*flike/double(evtsinlik);
1860  // fval = -2.*flike;
1861  // if( lowStatPenalty ) {
1862  // fval *= 100;
1863  // }
1864  }
1865  else {
1866  fval = -2.*flike;
1867  }
1868  }
1869  else {
1870  std::cout << "Problem: Events in likelihood = " << evtsinlik << std::endl;
1871  fval = 999999999.;
1872  }
1873  // fval = -2.*flike;
1874  if (MuScleFitUtils::debug>19)
1875  std::cout << "[MuScleFitUtils-likelihood]: End tree loop with likelihood value = " << fval << std::endl;
1876 
1877 // #ifdef DEBUG
1878 
1879 // if( MuScleFitUtils::minuitLoop_ < 10000 ) {
1883  }
1884  // }
1885  // else std::cout << "minuitLoop over 10000. Not filling histogram" << std::endl;
1886 
1887  std::cout<<"MINUIT loop number "<<MuScleFitUtils::minuitLoop_<<", likelihood = "<<fval<<std::endl;
1888 
1889  if( MuScleFitUtils::debug > 0 ) {
1890  // if( MuScleFitUtils::duringMinos_ ) {
1891  // int parnumber = (int)(MuScleFitUtils::parResol.size()+MuScleFitUtils::parScale.size()+
1892  // MuScleFitUtils::parCrossSection.size()+MuScleFitUtils::parBgr.size());
1893  // std::cout << "[MuScleFitUtils-likelihood]: Looping on tree with ";
1894  // for (int ipar=0; ipar<parnumber; ipar++) {
1895  // std::cout << "Parameter #" << ipar << " with value " << xval[ipar] << " ";
1896  // }
1897  // std::cout << std::endl;
1898  // std::cout << "[MuScleFitUtils-likelihood]: likelihood value = " << fval << std::endl;
1899  // }
1900  std::cout << "Events in likelihood = " << evtsinlik << std::endl;
1901  std::cout << "Events out likelihood = " << evtsoutlik << std::endl;
1902  }
1903 
1904 // #endif
1905 }
static std::vector< int > doScaleFit
static unsigned int loopCounter
static int debug
static unsigned int normalizationChanged_
Definition: weight.py:1
static double massProb(const double &mass, const double &rapidity, const int ires, const double &massResol)
reco::Particle::LorentzVector lorentzVector
Definition: GenMuonPair.h:9
static std::vector< std::pair< lorentzVector, lorentzVector > > ReducedSavedPair
static double massResolution(const lorentzVector &mu1, const lorentzVector &mu2)
static int minuitLoop_
static double computeWeight(const double &mass, const int iev, const bool doUseBkgrWindow=false)
static std::vector< std::pair< lorentzVector, lorentzVector > > SavedPair
static double invDimuonMass(const lorentzVector &mu1, const lorentzVector &mu2)
static lorentzVector applyScale(const lorentzVector &muon, const std::vector< double > &parval, const int charge)
static TMinuit * rminPtr_
static TH1D * likelihoodInLoop_
static double oldNormalization_
static int iev_
static bool normalizeLikelihoodByEventNumber_
Double_t lorentzianPeak ( Double_t *  x,
Double_t *  par 
)

See header file for a class description

Author
S. Bolognesi - INFN Torino / T. Dorigo, M. De Mattia - INFN Padova Revised S. Casasso, E. Migliore - UniTo & INFN Torino

Definition at line 56 of file MuScleFitUtils.cc.

References MillePedeFileConverter_cfg::e, Max(), and Pi.

Referenced by MuScleFitUtils::fitMass().

56  {
57  return (0.5*par[0]*par[1]/TMath::Pi()) /
58  TMath::Max(1.e-10,(x[0]-par[2])*(x[0]-par[2]) + .25*par[1]*par[1]);
59 }
const double Pi
T x() const
Cartesian x coordinate.
T Max(T a, T b)
Definition: MathUtil.h:44

Variable Documentation

double f[11][100]

Definition at line 78 of file MuScleFitUtils.cc.

Referenced by funct::SimplSumOrd< F, G, H >::prod0::a(), funct::SimplSumOrd< F, G, H >::prod1::a(), funct::RatioP1< F, G, H >::prod0::a(), funct::Product< ProductStruct< F, G >, H >::prod0::a(), funct::RatioP1< F, G, H >::prod1::a(), funct::Product< ProductStruct< F, G >, H >::prod1::a(), funct::RatioP2< F, G, H >::prod0::a(), funct::RatioP2< F, G, H >::prod1::a(), abcd(), abConnect(), abNum(), LocalTrajectoryParameters::absdz(), HGCDigitizer::accumulate(), ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::accumulate(), edm::helper::Filler< Association< C > >::add(), FWRecoGeometryESProducer::ADD_PIXEL_TOPOLOGY(), pftools::SpaceManager::addEvolution(), SiPixelFedCablingTree::addFed(), DDQuery::addFilter(), GEMDigiModel::addLinks(), GroupedCkfTrajectoryBuilder::advanceOneLayer(), SiStripGainFromData::algoEndJob(), PedsFullNoiseAlgorithm::analyse(), ZMuPtScaleAnalyzer::analyze(), EcalSRCondTools::analyze(), JetResolutionDBReader::analyze(), EwkMuLumiMonitorDQM::analyze(), GctFibreAnalyzer::analyze(), SiPixelHitEfficiencySource::analyze(), JetResolutionScaleFactorDBReader::analyze(), DumpSimGeometry::analyze(), HeavyFlavorValidation::analyze(), BTagPerformaceRootProducerFromSQLITE::analyze(), FFTJetPileupAnalyzer::analyze(), MuonTrackValidator::analyze(), ThirdHitPredictionFromCircle::angle(), apply(), npstat::ArrayND< Numeric, StackLen, StackDim >::apply(), FFTJetCorrectionProducer::applyCorrections(), approx_asin_P< 11 >(), approx_asin_P< 7 >(), approx_asin_P< 9 >(), approx_erf(), npstat::ArrayND< Numeric, StackLen, StackDim >::ArrayND(), npstat::ArrayND< Numeric, StackLen, StackDim >::assign(), OuterHitPhiPrediction::asym(), DTDigitizer::asymGausSmear(), HcalShapeIntegrator::at(), HcalItemArrayCollById< Item, N >::at(), fastmath::atan2r(), QGLikelihoodDBWriter::beginJob(), ProtonTaggerFilter::beginJob(), EcalPedestalHistory::beginRun(), TrackerOfflineValidation::bookDirHists(), DQMStore::bookTransaction(), MagGeoBuilderFromDDD::bRod::bRod(), MagGeoBuilderFromDDD::bSector::bSector(), FWCaloRecHitDigitSetProxyBuilder::build(), FWPFCandidateDetailView::build(), FWTrackHitsDetailView::build(), FWConvTrackHitsDetailView::build(), MagGeoBuilderFromDDD::build(), GenericSimClusterMapper::buildClusters(), FWBoxRecHit::buildTower(), funct::SimplSumOrd< F, G, H >::prod2::c(), funct::RatioP1< F, G, H >::prod2::c(), funct::Product< ProductStruct< F, G >, H >::prod2::c(), funct::RatioP2< F, G, H >::prod2::c(), SiTrackerMultiRecHitUpdator::calcParameters(), Phase2TrackerDigitizerAlgorithm::calcQ(), SiPixelDigitizerAlgorithm::calcQ(), HLTEcalResonanceFilter::calcShowerShape(), HLTRegionalEcalResonanceFilter::calcShowerShape(), PFEGammaAlgo::calculate_ele_mva(), Basic2DGenericPFlowPositionCalc::calculateAndSetPositionActual(), tkDetUtil::calculatePhiWindow(), l1t::EmtfPtAssignment::calculatePt(), CalculatePt_FullPrecision(), Quantile::calculateQ(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::calibFit(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::calibFitErr(), do_iterate< TupleType, N >::call(), edm::ExceptionCollector::call(), reco::HitPattern::call(), HFPreRecHit::chargeAsymmetry(), GenParticleProducer::chargeTimesThree(), check_file(), RPCDCSSummary::checkDCSbit(), RPCDataCertification::checkFED(), checkPhiInRange(), HitPairGeneratorFromLayerPairForPhotonConversion::checkRZCompatibilityWithSeedTrack(), RectangularEtaPhiTrackingRegion::checkRZOld(), FastTrajectoryCleaner::clean(), CmsShowEDI::clearPBFrame(), LaserSorter::closeOutStream(), Grid1D::closestNode(), DTTrigGeom::CMSPosition(), CmsShowCommonPopup::CmsShowCommonPopup(), CmsShowMain::CmsShowMain(), funct::SimplSumOrd< A, B, ProductStruct< C, D > >::combine(), funct::combine(), funct::RatioP1< A, B, Numerical< n > >::combine(), funct::Product< ProductStruct< F, G >, H >::combine(), funct::RatioP2< A, B, ProductStruct< C, D > >::combine(), TrackAssociatorByChi2Impl::compareTracksParam(), ForwardDetLayer::compatible(), GsfBetheHeitlerUpdator::compute(), EcalUncalibRecHitRatioMethodAlgo< C >::computeAmplitudeImpl(), FWGeoTopNode::ComputeBBox(), VolumeEnergyLossEstimator::computeElectrons(), TBPLayer::computeIndexes(), MagneticField::computeNominalValue(), EcalUncalibRecHitRatioMethodAlgo< C >::computeTime(), HBHENegativeFlagSetter::configFilter(), ConnectionManager::connect(), OpticalObject::constructMaterial(), convert64(), convertPackedEtaToPackedY(), edm::ParameterSetConverter::convertParameterSets(), edm::copy_all(), npstat::ArrayND< Numeric, StackLen, StackDim >::copyRangeLoopFunct(), EcalClusterToolsT< noZS >::covariances(), StripCPE::Param::coveredStrips(), spu::create_file(), EveService::createEventNavigationGUI(), muonisolation::IsolatorByNominalEfficiency::cuts(), GaussianSumUtilities1D::d1LnPdf(), GaussianSumUtilities1D::d2LnPdf(), dd_html_frameset(), dd_html_menu(), dd_html_menu_frameset(), dd_to_html(), DDDefaultNumberingScheme::DDDefaultNumberingScheme(), BTagLikeDeDxDiscriminator::dedx(), HcalTimeSlewSim::delay(), DDMaterial::density(), FFTJetProducer::determinePileup(), FFTJetProducer::determinePileupDensityFromDB(), FWEveDigitSetScalableMarkerGL::DirectDraw(), LocalTrajectoryParameters::direction(), MuonTagger::discriminator(), ElectronTagger::discriminator(), npstat::HistoND< Numeric, Axis >::dispatch(), DumpGctDigis::doFibres(), DDNsGenerator::doit(), DDFrameGenerator::doit(), CmsShowTaskExecutor::doNextTaskImp(), RPCEventSummary::dqmEndLuminosityBlock(), FP420G4Hit::Draw(), BscG4Hit::Draw(), PixelCPEBase::driftDirection(), DTDigitizer::driftTimeFromParametrization(), edm::JobReport::dumpFiles(), VVIObjFDetails::dzero(), VVIObjDetails::dzero(), sistripvvi::VVIObjDetails::dzero(), EcalSeverityLevelAlgo::EcalSeverityLevelAlgo(), Conv::elec(), EcalTBHodoscopeRawInfoDumper::endJob(), EcalTBTDCRawInfoDumper::endJob(), miscalibExample::endJob(), EcalSimpleTBAnalyzer::endJob(), EcalSimple2007H4TBAnalyzer::endJob(), PhiSymmetryCalibration::endJob(), HFPreRecHit::energyAsymmetry(), fireworks::energyScaledBox3DCorners(), AbsPlan1RechitCombiner::energyWeightedAverage(), BinomialProbability::error(), MagGeoBuilderFromDDD::eSector::eSector(), ForwardMeasurementEstimator::estimate(), BarrelMeasurementEstimator::estimate(), RectangularEtaPhiTrackingRegion::estimator(), fireworks::etScaledBox3DCorners(), ESRecHitSimAlgo::evalAmplitude(), fit::RootMinuitFuncEvaluator< Function >::evaluate(), fit::RootMinuitFuncEvaluator< Likelihood< Sample, PDF, Yield > >::evaluate(), PixelQuadrupletGenerator::QuantityDependsPt::evaluator(), CAHitTripletGenerator::QuantityDependsPt::evaluator(), CAHitQuadrupletGenerator::QuantityDependsPt::evaluator(), edm::JobReport::eventReadFromFile(), edm::JobReport::eventWrittenToFile(), FWGUIManager::exportAllViews(), TFParams::f3deg(), HitQuadrupletGeneratorFromLayerPairForPhotonConversion::failCheckRZCompatibility(), EEDetId::fastValidDetId(), HiEgammaSCEnergyCorrectionAlgo::fBrem(), VVIObj::fcn(), VVIObjF::fcn(), sistripvvi::VVIObj::fcn(), HiEgammaSCEnergyCorrectionAlgo::fEta(), HiEgammaSCEnergyCorrectionAlgo::fEtEta(), PixelRecoUtilities::FieldAt0::FieldAt0(), MultipleScatteringX0Data::fileName(), XHistogram::fill(), big::bigStore::fill(), QcdUeDQM::fillChargedJetSpectra(), PhiSymmetryCalibration_step2::fillConstantsHistos(), PhiSymmetryCalibration_step2_SM::fillConstantsHistos(), PhiSymmetryCalibration_step2::fillHistos(), PhiSymmetryCalibration_step2_SM::fillHistos(), StripCPE::fillParams(), GEDPhotonProducer::fillPhotonCollection(), LRHelpFunctions::fillToBackgroundHists(), LRHelpFunctions::fillToSignalHists(), edm::AssociationMap< edm::OneToMany< reco::BasicJetCollection, reco::TrackCollection > >::find(), MSLayersAtAngle::findLayer(), reco::findMethod(), muonisolation::IsolatorByNominalEfficiency::findPath(), CSCXonStrip_MatchGatti::findXOnStrip(), SymmetryFit::fit(), CosmicMuonSmoother::fit(), SymmetryFit::fitfunction(), CSCFindPeakTime::fivePoleFitTime(), PhotonFix::fixedEnergy(), extend_argparse::FixedWidthFormatter(), EgammaSCEnergyCorrectionAlgo::fNCrystals(), HiEgammaSCEnergyCorrectionAlgo::fNCrystals(), edm::for_all(), fpeak_convoluted(), cond::SmallWORMDict::IterHelp::frame(), lumi::from_string(), HcalDbOmds::from_string(), root::helper::RootFunctionHelper< F, args, Tag >::fun(), root::function(), npstat::ArrayND< Numeric, StackLen, StackDim >::functorFill(), npstat::ArrayND< Numeric, StackLen, StackDim >::functorFillLoop(), FWEveViewManager::FWEveViewManager(), FWGeometryTableViewManager::FWGeometryTableViewManager(), FWGUIEventFilter::FWGUIEventFilter(), FWLegoCandidate::FWLegoCandidate(), FWLegoEvePFCandidate::FWLegoEvePFCandidate(), FWTableViewManager::FWTableViewManager(), FWTriggerTableViewManager::FWTriggerTableViewManager(), EcalSimRawData::genFeData(), EcalSimRawData::genSrData(), funct::NthDerivative< n, X, F >::get(), funct::NthDerivative< 0, X, F >::get(), funct::ConstPrimitive< X, F, true >::get(), FTSFromVertexToPointFactory::get(), Phase2ITPixelClusterShapeCache::get(), SiPixelClusterShapeCache::get(), edm::AssociationMap< edm::OneToMany< reco::BasicJetCollection, reco::TrackCollection > >::get(), StripCPE::getAlgoParam(), HcalItemArrayCollById< Item, N >::getByIndex(), EcalBarrelGeometry::getCells(), track_associator::getChi2(), hgcal::ClusterTools::getClusterHadronFraction(), fireworks::GetColorValuesForPaletteExtra(), SCRegressionCalculator< VarCalc >::getCorrectionWithErrors(), HitPairGeneratorFromLayerPairForPhotonConversion::getCot(), TouchableToHistory::getFilteredView(), WatcherStreamFileReader::getInputFile(), PhiSymmetryCalibration::getKfactors(), EcalLaserDbService::getLaserCorrection(), EcalTrivialConditionRetriever::getMappingFromConfiguration(), FWMagField::GetMaxFieldMag(), hgcal::ClusterTools::getMultiClusterPosition(), EcalMappingElectronicsHandler::getNewObjects(), popcon::EcalSRPHandler::getNewObjects(), popcon::EcalTPGFineGrainEBIdMapHandler::getNewObjects(), EcalRecHit::getPower10(), ConversionProducer::getTrackImpactPosition(), EcalClusterLocalContCorrection::getValue(), EcalClusterCrackCorrection::getValue(), CSCAnalogSignal::getValue(), hcaldqm::quantity::FlagQuantity::getValue(), gen::Hydjet2Hadronizer::GetWeakDecayLimit(), GlobalTrajectoryParameters::GlobalTrajectoryParameters(), TIBRing::groupedCompatibleDetsV(), FWConfigurationManager::guessAndReadFromFile(), CmsAnnotation::Handle(), FWColorPopup::HandleButton(), CACell::haveSimilarCurvature(), HelixBarrelCylinderCrossing::HelixBarrelCylinderCrossing(), HelixBarrelPlaneCrossing2OrderLocal::HelixBarrelPlaneCrossing2OrderLocal(), HGCalRecHitWorkerSimple::HGCalRecHitWorkerSimple(), CSCMake2DRecHit::hitFromStripAndWire(), HitPairGeneratorFromLayerPairForPhotonConversion::hitPairs(), PixelQuadrupletGenerator::hitQuadruplets(), RectangularEtaPhiTrackingRegion::hits(), MultiHitGeneratorFromChi2::hitSets(), PixelTripletLargeTipGenerator::hitTriplets(), PixelTripletHLTGenerator::hitTriplets(), HLTCommonConfig::HLTCommonConfig(), IEEE32toDSP(), pat::strbitset::index(), SiTrivialInduceChargeOnStrips::induceVector(), ThirdHitCorrection::init(), PixelRecoLineRZ::initCot(), fireworks::Context::initEveElements(), PFEGammaAlgo::initializeProtoCands(), CandMatcherBase< C1, C2 >::initMaps(), reco::MustacheKernel::inMustache(), edm::JobReport::inputFileClosed(), edm::helper::Filler< Association< C > >::insert(), SimpleDiskBounds::inside(), DiskSectorBounds::inside(), GeneralNSurfaceDelimitedBounds::inside(), RectangularPlaneBounds::inside(), Bounds::inside(), PhotonConversionTrajectorySeedProducerFromSingleLegAlgo::inspectTrack(), funct::integral(), LinearGridInterpolator3D::interpolate(), SiStripTemplate::interpolate(), SiPixelTemplate::interpolate(), npstat::LinInterpolatedTableND< Numeric, Axis >::invertRatioResponse(), cond::persistency::iovSorter(), ParabolicParametrizedMagneticField::isDefined(), OAEParametrizedMagneticField::isDefined(), heppy::IsolationComputer::isoSumNeutralsWeighted(), edm::isTransientEqual(), EcalSimRawData::iTtPhi2cIndex(), PerigeeConversions::jacobianParameters2Cartesian(), PerigeeConversions::jacobianPerigee2Curvilinear(), JME::JetResolutionObject::JetResolutionObject(), FWGUIValidatingTextEntry::keyPressedInPopup(), hitfit::Gentop_Args::kt_res_str(), StripCPEfromTrackAngle::legacyStripErrorSquared(), EvolutionECAL::LightCollectionEfficiency(), EvolutionECAL::LightCollectionEfficiencyWeighted(), PFMultiDepthClusterizer::link(), GaussianSumUtilities1D::lnPdf(), FWFFLooper::loadDefaultGeometryFile(), FFTJetPileupProcessor::loadFlatteningFactors(), MatacqProducer::loadOrbitOffset(), HcalPatternSource::loadPatternFile(), ClusterShapeHitFilter::loadPixelLimits(), SteppingHelixPropagator::loadState(), EcalClusterToolsT< noZS >::localCovariances(), PSimHit::localDirection(), RectangularStripTopology::localError(), PixelCPETemplateReco::localError(), CSCRadialStripTopology::localError(), TkRadialStripTopology::localError(), PixelCPEGeneric::localError(), StripCPEfromTemplate::localParameters(), StripCPE::localParameters(), TrapezoidalStripTopology::localPitch(), CSCRadialStripTopology::localPitch(), RectangularStripTopology::localPosition(), TrapezoidalStripTopology::localStripLength(), CSCRadialStripTopology::localStripLength(), TkRadialStripTopology::localStripLength(), GroupedTrajCandLess::looperPenalty(), Lut::Lut(), SimpleHBHEPhase1Algo::m0Time(), GlobalTrajectoryParameters::magneticFieldInInverseGeV(), main(), SiPixelDigitizerAlgorithm::make_digis(), MuonRPCDetLayerGeometryBuilder::makeBarrelLayers(), MuonRPCDetLayerGeometryBuilder::makeBarrelRods(), HcalDbHardcode::makeGain(), fireworks_root_gui::makeHorizontalFrame(), HcalDbHardcode::makeMCParam(), HcalDbHardcode::makePedestal(), reco::MaskCollection< Object >::mask(), edm::EventSelector::maskTriggerResults(), funct::master(), SiStripRecHitMatcher::match(), CSCRadialStripTopology::measurementError(), TkRadialStripTopology::measurementError(), DQMStore::meBookerGetter(), DQMStore::meGetter(), JetPartonNamespace::UserPartonMixture::mixt(), LocalTrajectoryParameters::momentum(), HelixArbitraryPlaneCrossing::notAtSurface(), edm::AssociationMap< edm::OneToMany< reco::BasicJetCollection, reco::TrackCollection > >::numberOfAssociations(), LRHelpFunctions::obsFitIncluded(), CaloSamples::offsetTime(), oldComputeElectrons(), GsiFTPStorageMaker::open(), HttpStorageMaker::open(), DavixStorageMaker::open(), LocalStorageMaker::open(), DCacheStorageMaker::open(), XrdStorageMaker::open(), StormStorageMaker::open(), StormLcgGtStorageMaker::open(), RFIOStorageMaker::open(), GsfEleEffAreaPFIsoCut::operator()(), FreeFunctionSelector< T, f >::operator()(), GsfEleDeltaBetaIsoCut::operator()(), GsfEleDeltaBetaIsoCutStandalone::operator()(), GsfEleCalPFClusterIsoCut::operator()(), TkClonerImpl::operator()(), BSpdfsFcn::operator()(), CandMatcherBase< C1, C2 >::operator()(), funct::GaussLegendreIntegrator::operator()(), funct::GaussIntegrator::operator()(), reco::operator<<(), funct::operator<<(), reco::parser::operator<<(), operator<<(), cscdqm::LogCout::operator<<(), edm::LogWarning::operator<<(), edm::LogError::operator<<(), edm::LogSystem::operator<<(), edm::LogInfo::operator<<(), cms::operator<<(), edm::LogVerbatim::operator<<(), edm::LogPrint::operator<<(), edm::LogProblem::operator<<(), edm::LogImportant::operator<<(), edm::LogAbsolute::operator<<(), edm::LogDebug_::operator<<(), edm::LogTrace_::operator<<(), edm::edmmltest::LogWarningThatSuppressesLikeLogInfo::operator<<(), MSLayer::operator==(), edm::JobReport::outputFileClosed(), PhiSymmetryCalibration_step2_SM::outResidHistos(), PhiSymmetryCalibration_step2::outResidHistos(), npstat::BoxND< Numeric >::overlapFraction(), Phase2EndcapLayer::overlapInR(), TIDLayer::overlapInR(), pat::PackedCandidate::packVtx(), CSCFindPeakTime::parabolaFitTime(), GenParticlePruner::parse(), parseFFTJetScaleCalculator(), EcalTrigTowerConstituentsMapBuilder::parseTextMap(), StraightLinePlaneCrossing::pathLength(), PetrukhinFunc(), PFResolutionMapManager::PFResolutionMapManager(), ThirdHitPredictionFromCircle::phi(), CSCSegAlgoSK::phiAtZ(), CSCSegAlgoRU::phiAtZ(), CSCSegAlgoTC::phiAtZ(), InnerDeltaPhi::phiRange(), SiPixelTemplateReco::PixelTempReco2D(), SiPixelTemplateSplit::PixelTempSplit(), CmsShowMainFrame::playEventsBackwardsAction(), HIPplots::plotAlignableChiSquare(), HIPplots::plotAlignParams(), HIPplots::plotAlignParamsAtIter(), plotMaker(), TwoTrackMinimumDistance::pointsHelixHelix(), reco::exprEvalDetails::popenCPP(), FWRPZView::populateController(), FWTriggerTableView::populateController(), HelixBarrelPlaneCrossing2OrderLocal::positionOnly(), FWFFService::postBeginJob(), InputTagDistributorService::preModule(), VariableHelperService::preModule(), TrackerOfflineValidation::prepareSummaryHists(), fireworks::prepareTrack(), DeDxHitInfoProducer::processHit(), MultiTrackSelector::processMVA(), HIMultiTrackSelector::processMVA(), RawDataUnpacker::ProcessVFATDataParallel(), edm::ExpoRandomPGunProducer::produce(), PFJetsMatchedToFilteredCaloJetsProducer::produce(), L1HLTJetsMatching::produce(), MCTrackMatcher::produce(), GenTrackMatcher::produce(), GenParticleProducer::produce(), BeamSpotOnlineProducer::produce(), L3TkMuonProducer::produce(), MatchedProbeMaker< T >::produce(), L1TMicroGMTInputProducer::produce(), GoodSeedProducer::produce(), EcalTrivialConditionRetriever::produceEcalAlignmentEB(), EcalTrivialConditionRetriever::produceEcalAlignmentEE(), EcalTrivialConditionRetriever::produceEcalAlignmentES(), EcalTrivialConditionRetriever::produceEcalSamplesCorrelation(), perftools::EdmEventSize::produceHistos(), TkClonerImpl::project(), tauImpactParameter::ErrorMatrixPropagator::propagateError(), Geant4ePropagator::propagateGeneric(), AnalyticalPropagator::propagateParametersOnCylinder(), AnalyticalPropagator::propagateParametersOnPlane(), gen::PtYDistributor::PtYDistributor(), EcalUncalibRecHitRecAnalFitAlgo< EBDataFrame >::pulseShapeFunction(), Histos::put(), SiPixelGenError::qbin(), SiStripTemplate::qbin(), SiPixelTemplate::qbin(), SiPixelTemplate::qbin_dist(), cms::Exception::raise(), PixelRecoLineRZ::rAtZ(), edm::ValueMap< edm::Ref< std::vector< reco::PFCandidate >, reco::PFCandidate, edm::refhelper::FindUsingAdvance< std::vector< reco::PFCandidate >, reco::PFCandidate > > >::rawIndexOf(), hitfit::Defaults_Textrep::read_file(), TStorageFactoryFile::ReadBufferAsync(), DQMStore::readFile(), FWConfigurationManager::readFromFile(), reco::details::readGzipFile(), EcalFEtoDigi::readInput(), MuonResiduals5DOFFitter::readNtuple(), MuonResiduals6DOFrphiFitter::readNtuple(), MuonResiduals6DOFFitter::readNtuple(), PasswordReader::readPassword(), TkConvValidator::recalculateMomentumAtFittedVertex(), HcalSimpleRecAlgoImpl::recoHFTime(), HFSimpleTimeCheck::reconstruct(), MiniFloatConverter::reduceMantissaToNbits(), CmsShowNavigator::removeFilter(), PFEGammaAlgo::removeOrLinkECALClustersToKFTracks(), edm::JobReport::reportDataType(), edm::JobReport::reportFastCopyingStatus(), ViewerParameterGUI::reset(), CastorSimpleRecAlgo::resetTimeSamples(), EvolutionECAL::ResolutionConstantTermEM50GeV(), LA_Filler_Fitter::result(), edm::RootEmbeddedFileSequence::RootEmbeddedFileSequence(), fit::RootMinuit< Function >::RootMinuit(), LocalError::rotate(), LocalError::rotateCosSin(), RPAlignmentCorrectionsMethods::RPAlignmentCorrectionsMethods(), PixelFitterByHelixProjections::run(), SimpleElectronicsSimInMIPs::run(), HGCHEbackDigitizer::runCaliceLikeDigitizer(), HGCFEElectronics< DFr >::runShaperWithToT(), HGCDigitizerBase< DFr >::runSimple(), HGCFEElectronics< DFr >::runSimpleShaper(), safe_atan2f(), DQMStore::save(), TkHistoMap::saveAsCanvas(), npstat::ArrayND< Numeric, StackLen, StackDim >::scanInPlace(), EcalClusterToolsT< noZS >::scLocalCovariances(), GroupedTrajCandLess::score(), FWViewContextMenuHandlerGL::select(), reco::SelectInCollection< Object >::select(), reco::SelectIndecesInCollection< Object >::select(), MultiTrackSelector::select(), HIMultiTrackSelector::select(), siStripClusterTools::sensorThicknessInverse(), VariableHelperService::set(), external::HEPTopTaggerV2::set_optimalR_calc_fun(), npstat::HistoND< Numeric, Axis >::setAxisLabel(), FWGeometryTableManagerBase::NodeInfo::setBit(), TrackInformation::setCaloIDChecked(), CmsShowMainBase::setConfigFilename(), EcalRecHit::setEnergyError(), TotemG4Hit::setEntry(), CompositeTrajectoryFilter::setEvent(), UnbinnedLikelihoodFit::setFunction(), PFECALSuperClusterAlgo::setMajorityFraction(), hcalCalib::SetMaxProbeJetEmFrac(), hcalCalib::SetMaxTagJetEmFrac(), BscG4Hit::setPhiAtEntry(), FP420G4Hit::setPhiAtEntry(), TotemG4Hit::setPhiAtEntry(), BaseParticlePropagator::setPropagationConditions(), FWMET3DProxyBuilder::Arrow::setScale(), HcalPulseShape::setShapeBin(), CastorPulseShapes::Shape::setShapeBin(), FWEveView::setupEnergyScale(), DTConfigLUTs::setWHEEL(), sicif(), PhotonFix::sigmaEnergy(), SimHitPrinter::SimHitPrinter(), SimHitMatcher::simHitsMeanStrip(), SiPixelTemplate::simpletemplate2D(), sortNtupleByEvent(), XHistogram::splitSegment(), FWFFLooper::startingNewLoop(), npstat::ArrayND< Num1, Len1, Dim1 >::strides(), stringToNumber(), cscdqm::stringToNumber(), SiStripTemplateReco::StripTempReco1D(), SiStripTemplateSplit::StripTempSplit(), FWGUIManager::subviewIsBeingDestroyed(), MSLayersAtAngle::sum2RmRn(), SummaryPlotXmlParser::SummaryPlotXmlParser(), MSLayersAtAngle::sumX0D(), MSLayer::sumX0D(), SiStripTemplate::sxtemp(), SiPixelTemplate::temperrors(), MinPtTrajectoryFilter::test(), FWGeometryTableManagerBase::NodeInfo::testBit(), StripSubClusterShapeFilterBase::testLastHit(), edm::EventSelector::testSelectionOverlap(), TkStripMeasurementDet::testStrips(), root::tf1(), TFileAdaptor::TFileAdaptor(), CaloParticle::threeCharge(), SimCluster::threeCharge(), TrackingParticle::threeCharge(), timeout::timed_out(), timeout::TimedOutFn(), Traj2TrackHits::toChi2x5(), GlobalTrajectoryParameters::transverseCurvature(), ThirdHitPredictionFromCircle::transverseIP(), funct::trapezoid_integral(), reco::typeCode(), PFEGammaAlgo::unlinkRefinableObjectKFandECALMatchedToHCAL(), PFEGammaAlgo::unlinkRefinableObjectKFandECALWithBadEoverP(), HcalUnpacker::unpack(), unsafe_acos(), unsafe_asin(), spu::Untar(), edm::service::CondorStatusService::updateImpl(), LocalTrajectoryParameters::updateP(), FWPFRhoPhiRecHit::updateScale(), root::helper::RootVarsAdapter< F, 1 >::value(), GsfEleDeltaBetaIsoCutStandalone::value(), GsfEleEffAreaPFIsoCut::value(), GsfEleDeltaBetaIsoCut::value(), GsfEleCalPFClusterIsoCut::value(), root::helper::RootVarsAdapter< F, 2 >::value(), JetPartonNamespace::ParametrizationJetParton::value(), SiStripTemplate::vavilov2_pars(), SiPixelTemplate::vavilov2_pars(), SiStripTemplate::vavilov_pars(), SiPixelTemplate::vavilov_pars(), FWCaloRecHitDigitSetProxyBuilder::viewContextBoxScale(), FWPFCandidateWithHitsProxyBuilder::viewContextBoxScale(), VVIObjF::VVIObjF(), WatcherStreamFileReader::WatcherStreamFileReader(), reco::parser::Exception::what(), sistrip::MeasureLA::write_report_text_ms(), edm::JobReport::JobReportImpl::writeInputFile(), edm::JobReport::JobReportImpl::writeOutputFile(), MSLayer::x0(), SiStripTemplate::xflcorr(), SiPixelTemplate::xflcorr(), geometryXMLparser.MuonGeometry::xml(), TkRadialStripTopology::xOfStrip(), SiStripTemplate::xsigma2(), SiPixelTemplate::xsigma2(), SiPixelTemplate2D::xysigma2(), SiPixelTemplate2D::xytemp(), SiPixelTemplate::yflcorr(), CSCLayerGeometry::yLimitsOfStripPlane(), SiPixelTemplate::ysigma2(), hitfit::Base_Constrainer::~Base_Constrainer(), cscdqm::Configuration::~Configuration(), tauImpactParameter::ErrorMatrixPropagator::~ErrorMatrixPropagator(), edm::service::InitRootHandlers::~InitRootHandlers(), SiPixelFedCablingTree::~SiPixelFedCablingTree(), and TestPythiaDecays::~TestPythiaDecays().

double g[11][100]

Definition at line 79 of file MuScleFitUtils.cc.

TF1* GL
Initial value:
= new TF1 ("GL",
"0.5/3.1415926*[0]/(pow(x-[1],2)+pow(0.5*[0],2))*exp(-0.5*pow((x-[2])/[3],2))/([3]*sqrt(6.283185))",
0, 1000)

Definition at line 83 of file MuScleFitUtils.cc.

Referenced by MuScleFitUtils::massProb(), and MuScleFitBase::readProbabilityDistributionsFromFile().

TF2* GL2
Initial value:
= new TF2 ("GL2",
"0.5/3.1415926*[0]/(pow(x-[1],2)+pow(0.5*[0],2))*exp(-0.5*pow((x-y)/[2],2))/([2]*sqrt(6.283185))",
0, 200, 0, 200)

Definition at line 87 of file MuScleFitUtils.cc.

double isum
double mzsum

Definition at line 76 of file MuScleFitUtils.cc.