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 62 of file MuScleFitUtils.cc.

References create_public_lumi_plots::exp.

Referenced by MuScleFitUtils::fitReso().

                                               {
  return par[0]*exp(-0.5*((x[0]-par[1])/par[2])*((x[0]-par[1])/par[2]));
}

void likelihood	(	int &	npar,
		double *	grad,
		double &	fval,
		double *	xval,
		int	flag
	)

Definition at line 1696 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_, create_public_lumi_plots::log, MuScleFitUtils::loopCounter, MuScleFitUtils::massProb(), MuScleFitUtils::massResolution(), MuScleFitUtils::minuitLoop_, MuScleFitUtils::normalizationChanged_, MuScleFitUtils::normalizeLikelihoodByEventNumber_, MuScleFitUtils::oldNormalization_, mix_2012_Summer_inTimeOnly_cff::prob, MuScleFitUtils::ReducedSavedPair, MuScleFitUtils::rminPtr_, MuScleFitUtils::SavedPair, and histoStyle::weight.

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

                                                                                            {

  if (MuScleFitUtils::debug>19) std::cout << "[MuScleFitUtils-likelihood]: In likelihood function" << std::endl;

  const lorentzVector * recMu1;
  const lorentzVector * recMu2;
  lorentzVector corrMu1;
  lorentzVector corrMu2;

  //   if (MuScleFitUtils::debug>19) {
  //     int parnumber = (int)(MuScleFitUtils::parResol.size()+MuScleFitUtils::parScale.size()+
  //                           MuScleFitUtils::parCrossSection.size()+MuScleFitUtils::parBgr.size());
  //     std::cout << "[MuScleFitUtils-likelihood]: Looping on tree with ";
  //     for (int ipar=0; ipar<parnumber; ipar++) {
  //       std::cout << "Parameter #" << ipar << " with value " << xval[ipar] << " ";
  //     }
  //     std::cout << std::endl;
  //   }

  // Loop on the tree
  // ----------------
  double flike = 0;
  int evtsinlik = 0;
  int evtsoutlik = 0;
  // std::cout << "SavedPair.size() = " << MuScleFitUtils::SavedPair.size() << std::endl;
  if( MuScleFitUtils::debug>0 ) {
    std::cout << "SavedPair.size() = " << MuScleFitUtils::SavedPair.size() << std::endl;
    std::cout << "ReducedSavedPair.size() = " << MuScleFitUtils::ReducedSavedPair.size() << std::endl;
  }
  // for( unsigned int nev=0; nev<MuScleFitUtils::SavedPair.size(); ++nev ) {
  for( unsigned int nev=0; nev<MuScleFitUtils::ReducedSavedPair.size(); ++nev ) {

    //     recMu1 = &(MuScleFitUtils::SavedPair[nev].first);
    //     recMu2 = &(MuScleFitUtils::SavedPair[nev].second);
    recMu1 = &(MuScleFitUtils::ReducedSavedPair[nev].first);
    recMu2 = &(MuScleFitUtils::ReducedSavedPair[nev].second);

    // Compute original mass
    // ---------------------
    double mass = MuScleFitUtils::invDimuonMass( *recMu1, *recMu2 );

    // Compute weight and reference mass (from original mass)
    // ------------------------------------------------------
    double weight = MuScleFitUtils::computeWeight(mass, MuScleFitUtils::iev_);
    if( weight!=0. ) {
      // Compute corrected mass (from previous biases) only if we are currently fitting the scale
      // ----------------------------------------------------------------------------------------
      if( MuScleFitUtils::doScaleFit[MuScleFitUtils::loopCounter] ) {
//  std::cout << "Original pt1 = " << corrMu1.Pt() << std::endl;
//  std::cout << "Original pt2 = " << corrMu2.Pt() << std::endl;
        corrMu1 = MuScleFitUtils::applyScale(*recMu1, xval, -1);
        corrMu2 = MuScleFitUtils::applyScale(*recMu2, xval,  1);
        
//         if( (corrMu1.Pt() != corrMu1.Pt()) || (corrMu2.Pt() != corrMu2.Pt()) ) {
//           std::cout << "Rescaled pt1 = " << corrMu1.Pt() << std::endl;
//           std::cout << "Rescaled pt2 = " << corrMu2.Pt() << std::endl;
//         }
//  std::cout << "Rescaled pt1 = " << corrMu1.Pt() << std::endl;
//  std::cout << "Rescaled pt2 = " << corrMu2.Pt() << std::endl;
      }
      else {
        corrMu1 = *recMu1;
        corrMu2 = *recMu2;

//         if( (corrMu1.Pt() != corrMu1.Pt()) || (corrMu2.Pt() != corrMu2.Pt()) ) {
//           std::cout << "Not rescaled pt1 = " << corrMu1.Pt() << std::endl;
//           std::cout << "Not rescaled pt2 = " << corrMu2.Pt() << std::endl;
//         }
      }
      double corrMass = MuScleFitUtils::invDimuonMass(corrMu1, corrMu2);
      double Y = (corrMu1+corrMu2).Rapidity();
      double resEta = (corrMu1+corrMu2).Eta();
      if( MuScleFitUtils::debug>19 ) {
    std::cout << "[MuScleFitUtils-likelihood]: Original/Corrected resonance mass = " << mass
         << " / " << corrMass << std::endl;
      }

      // Compute mass resolution
      // -----------------------
      double massResol = MuScleFitUtils::massResolution(corrMu1, corrMu2, xval);
      if (MuScleFitUtils::debug>19)
    std::cout << "[MuScleFitUtils-likelihood]: Resolution is " << massResol << std::endl;

      // Compute probability of this mass value including background modeling
      // --------------------------------------------------------------------
      if (MuScleFitUtils::debug>1) std::cout << "calling massProb inside likelihood function" << std::endl;

      // double prob = MuScleFitUtils::massProb( corrMass, resEta, Y, massResol, xval );
      double prob = MuScleFitUtils::massProb( corrMass, resEta, Y, massResol, xval, false, corrMu1.eta(), corrMu2.eta() );
      if (MuScleFitUtils::debug>1) std::cout << "likelihood:massProb = " << prob << std::endl;

      // Compute likelihood
      // ------------------
      if( prob>0 ) {
    // flike += log(prob*10000)*weight; // NNBB! x10000 to see if we can recover the problem of boundary
    flike += log(prob)*weight;
    evtsinlik += 1;  // NNBB test: see if likelihood per event is smarter (boundary problem)
      } else {
        if( MuScleFitUtils::debug > 0 ) {
          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;
          std::cout << "Original mass was = " << mass << std::endl;
      std::cout << "WARNING: massResol = " << massResol << " outside window" << std::endl;
        }
    evtsoutlik += 1;
      }
      if (MuScleFitUtils::debug>19)
    std::cout << "[MuScleFitUtils-likelihood]: Mass probability = " << prob << std::endl;
    } // weight!=0

  } // End of loop on tree events

//   // Protection for low statistic. If the likelihood manages to throw out all the signal
//   // events and stays with ~ 10 events in the resonance window it could have a better likelihood
//   // because of ~ uniformly distributed events (a random combination could be good and spoil the fit).
//   // We require that the number of events included in the fit does not change more than 5% in each minuit loop.
//   bool lowStatPenalty = false;
//   if( MuScleFitUtils::minuitLoop_ > 0 ) {
//     double newEventsOutInRatio = double(evtsinlik);
//     // double newEventsOutInRatio = double(evtsoutlik)/double(evtsinlik);
//     double ratio = newEventsOutInRatio/MuScleFitUtils::oldEventsOutInRatio_;
//     MuScleFitUtils::oldEventsOutInRatio_ = newEventsOutInRatio;
//     if( ratio < 0.8 || ratio > 1.2 ) {
//       std::cout << "Warning: too much change from oldEventsInLikelihood to newEventsInLikelihood, ratio is = " << ratio << std::endl;
//       std::cout << "oldEventsInLikelihood = " << MuScleFitUtils::oldEventsOutInRatio_ << ", newEventsInLikelihood = " << newEventsOutInRatio << std::endl;
//       lowStatPenalty = true;
//     }
//   }

  // It is a product of probabilities, we compare the sqrt_N of them. Thus N becomes a denominator of the logarithm.
  if( evtsinlik != 0 ) {

    if( MuScleFitUtils::normalizeLikelihoodByEventNumber_ ) {
      // && !(MuScleFitUtils::duringMinos_) ) {
      if( MuScleFitUtils::rminPtr_ == 0 ) {
        std::cout << "ERROR: rminPtr_ = " << MuScleFitUtils::rminPtr_ << ", code will crash" << std::endl;
      }
      double normalizationArg[] = {1/double(evtsinlik)};
      // Reset the normalizationArg only if it changed
      if( MuScleFitUtils::oldNormalization_ != normalizationArg[0] ) {
        int ierror = 0;
//         if( MuScleFitUtils::likelihoodInLoop_ != 0 ) {
//           // This condition is set only when minimizing. Later calls of hesse and minos will not change the value
//           // This is done to avoid minos being confused by changing the UP parameter during its computation.
//           MuScleFitUtils::rminPtr_->mnexcm("SET ERR", normalizationArg, 1, ierror);
//         }
        MuScleFitUtils::rminPtr_->mnexcm("SET ERR", normalizationArg, 1, ierror);
    std::cout << "oldNormalization = " << MuScleFitUtils::oldNormalization_ << " new = " << normalizationArg[0] << std::endl;
        MuScleFitUtils::oldNormalization_ = normalizationArg[0];
        MuScleFitUtils::normalizationChanged_ += 1;
      }
      fval = -2.*flike/double(evtsinlik);
      // fval = -2.*flike;
      //     if( lowStatPenalty ) {
      //       fval *= 100;
      //     }
    }
    else {
      fval = -2.*flike;
    }
  }
  else {
    std::cout << "Problem: Events in likelihood = " << evtsinlik << std::endl;
    fval = 999999999.;
  }
  // fval = -2.*flike;
  if (MuScleFitUtils::debug>19)
    std::cout << "[MuScleFitUtils-likelihood]: End tree loop with likelihood value = " << fval << std::endl;

//  #ifdef DEBUG

//  if( MuScleFitUtils::minuitLoop_ < 10000 ) {
  if( MuScleFitUtils::likelihoodInLoop_ != 0 ) {
    ++MuScleFitUtils::minuitLoop_;
    MuScleFitUtils::likelihoodInLoop_->SetBinContent(MuScleFitUtils::minuitLoop_, fval);
  }
  //  }
  // else std::cout << "minuitLoop over 10000. Not filling histogram" << std::endl;

  std::cout<<"MINUIT loop number "<<MuScleFitUtils::minuitLoop_<<", likelihood = "<<fval<<std::endl;

  if( MuScleFitUtils::debug > 0 ) {
    //     if( MuScleFitUtils::duringMinos_ ) {
    //       int parnumber = (int)(MuScleFitUtils::parResol.size()+MuScleFitUtils::parScale.size()+
    //                             MuScleFitUtils::parCrossSection.size()+MuScleFitUtils::parBgr.size());
    //       std::cout << "[MuScleFitUtils-likelihood]: Looping on tree with ";
    //       for (int ipar=0; ipar<parnumber; ipar++) {
    //         std::cout << "Parameter #" << ipar << " with value " << xval[ipar] << " ";
    //       }
    //       std::cout << std::endl;
    //       std::cout << "[MuScleFitUtils-likelihood]: likelihood value = " << fval << std::endl;
    //     }
    std::cout << "Events in likelihood = " << evtsinlik << std::endl;
    std::cout << "Events out likelihood = " << evtsoutlik << std::endl;
  }

//  #endif
}

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

Definition at line 55 of file MuScleFitUtils.cc.

References alignCSCRings::e, siStripFEDMonitor_P5_cff::Max, and Pi.

Referenced by MuScleFitUtils::fitMass().

                                                     {
  return (0.5*par[0]*par[1]/TMath::Pi()) /
  TMath::Max(1.e-10,(x[0]-par[2])*(x[0]-par[2]) + .25*par[1]*par[1]);
}

Variable Documentation

double f[11][100]

Definition at line 77 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< PROD_S(F, G), H >::prod0::a(), funct::RatioP1< F, G, H >::prod1::a(), funct::Product< PROD_S(F, G), H >::prod1::a(), funct::RatioP2< F, G, H >::prod0::a(), funct::RatioP2< F, G, H >::prod1::a(), abcd(), abConnect(), edm::helper::Filler< Association< C > >::add(), pftools::SpaceManager::addEvolution(), SiPixelFedCablingTree::addFed(), DDQuery::addFilter(), GroupedCkfTrajectoryBuilder::advanceOneLayer(), SiStripGainFromData::algoEndJob(), PedsFullNoiseAlgorithm::analyse(), ZMuPtScaleAnalyzer::analyze(), EcalSRCondTools::analyze(), EwkMuLumiMonitorDQM::analyze(), GctFibreAnalyzer::analyze(), HcalMonitorModule::analyze(), ZDCMonitorModule::analyze(), DumpSimGeometry::analyze(), HeavyFlavorValidation::analyze(), BTagPerformaceRootProducerFromSQLITE::analyze(), FFTJetPileupAnalyzer::analyze(), CastorMonitorModule::analyze(), MuonTrackValidator::analyze(), ThirdHitPredictionFromCircle::angle(), apply(), npstat::ArrayND< Numeric, StackLen, StackDim >::apply(), FFTJetCorrectionProducer::applyCorrections(), ThreeThresholdAlgorithm::applyGains(), approx_erf(), npstat::ArrayND< Numeric, StackLen, StackDim >::ArrayND(), npstat::ArrayND< Numeric, StackLen, StackDim >::assign(), DTDigitizer::asymGausSmear(), HcalShapeIntegrator::at(), fastmath::atan2r(), SiStripCluster::barycenter(), ProtonTaggerFilter::beginJob(), TSGFromL1Muon::beginRun(), RPCDataCertification::beginRun(), RPCDCSSummary::beginRun(), RPCEventSummary::beginRun(), MaterialForOnline::beginRun(), TrackerOfflineValidation::bookDirHists(), DQMStore::bookTransaction(), MagGeoBuilderFromDDD::bRod::bRod(), MagGeoBuilderFromDDD::bSector::bSector(), FWCaloRecHitDigitSetProxyBuilder::build(), FWTrackHitsDetailView::build(), FWConvTrackHitsDetailView::build(), MagGeoBuilderFromDDD::build(), funct::SimplSumOrd< F, G, H >::prod2::c(), funct::RatioP1< F, G, H >::prod2::c(), funct::Product< PROD_S(F, G), H >::prod2::c(), funct::RatioP2< F, G, H >::prod2::c(), SiPixelDigitizerAlgorithm::calcQ(), HLTEcalResonanceFilter::calcShowerShape(), PFEGammaAlgo::calculate_ele_mva(), Quantile::calculateQ(), edm::ExceptionCollector::call(), reco::HitPattern::call(), GenParticleProducer::chargeTimesThree(), HitPairGeneratorFromLayerPairForPhotonConversion::checkRZCompatibilityWithSeedTrack(), CmsShowEDI::clearPBFrame(), LaserSorter::closeOutStream(), Grid1D::closestNode(), CmsShowMain::CmsShowMain(), TrackerOfflineValidation::collateSummaryHists(), funct::SimplSumOrd< A, B, PROD_S(C, D) >::combine(), funct::combine(), funct::RatioP1< F, G, H >::combine(), funct::Product< PROD_S(F, G), H >::combine(), funct::RatioP2< F, G, H >::combine(), TrackAssociatorByChi2::compareTracksParam(), CompositeTrajectoryFilter::CompositeTrajectoryFilter(), EcalUncalibRecHitRatioMethodAlgo< C >::computeAmplitude(), FWGeoTopNode::ComputeBBox(), MagneticField::computeNominalValue(), EcalUncalibRecHitRatioMethodAlgo< C >::computeTime(), ConnectionManager::connect(), convert64(), edm::ParameterSetConverter::convertParameterSets(), edm::copy_all(), npstat::ArrayND< Numeric, StackLen, StackDim >::copyRangeLoopFunct(), 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(), FFTJetProducer::determinePileup(), FFTJetProducer::determinePileupDensityFromDB(), FWEveDigitSetScalableMarkerGL::DirectDraw(), npstat::HistoND< Numeric, Axis >::dispatch(), DumpGctDigis::doFibres(), DDNsGenerator::doit(), DDFrameGenerator::doit(), CmsShowTaskExecutor::doNextTaskImp(), SiStripRecHitMatcher::doubleMatch(), PixelCPEBase::driftDirection(), DTDigitizer::driftTimeFromParametrization(), edm::JobReport::dumpFiles(), VVIObjFDetails::dzero(), VVIObjDetails::dzero(), sistripvvi::VVIObjDetails::dzero(), EcalSeverityLevelAlgo::EcalSeverityLevelAlgo(), Conv::elec(), EcalTBHodoscopeRawInfoDumper::endJob(), EcalTBTDCRawInfoDumper::endJob(), RPCMonitorRaw::endJob(), miscalibExample::endJob(), EcalSimpleTBAnalyzer::endJob(), EcalSimple2007H4TBAnalyzer::endJob(), PhiSymmetryCalibration::endJob(), TrackerOfflineValidation::endJob(), fireworks::energyScaledBox3DCorners(), BinomialProbability::error(), MagGeoBuilderFromDDD::eSector::eSector(), fireworks::etScaledBox3DCorners(), ESRecHitSimAlgo::evalAmplitude(), fit::RootMinuitFuncEvaluator< Function >::evaluate(), fit::RootMinuitFuncEvaluator< Likelihood< Sample, PDF, Yield > >::evaluate(), edm::JobReport::eventReadFromFile(), edm::JobReport::eventWrittenToFile(), 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(), MTVHistoProducerAlgoForTracker::fill_generic_recoTrack_histos(), MTVHistoProducerAlgoForTracker::fill_recoAssociated_simTrack_histos(), QcdUeDQM::fillChargedJetSpectra(), PhiSymmetryCalibration_step2::fillConstantsHistos(), PhiSymmetryCalibration_step2_SM::fillConstantsHistos(), PhiSymmetryCalibration_step2::fillHistos(), PhiSymmetryCalibration_step2_SM::fillHistos(), StripCPE::fillParams(), PFEGammaAlgo::fillPFCandidates(), LRHelpFunctions::fillToBackgroundHists(), LRHelpFunctions::fillToSignalHists(), edm::AssociationMap< edm::OneToMany< reco::BasicJetCollection, reco::TrackCollection > >::find(), MSLayersAtAngle::findLayer(), muonisolation::IsolatorByNominalEfficiency::findPath(), SymmetryFit::fit(), SymmetryFit::fitfunction(), PhotonFix::fixedEnergy(), extend_argparse::FixedWidthFormatter(), EgammaSCEnergyCorrectionAlgo::fNCrystals(), HiEgammaSCEnergyCorrectionAlgo::fNCrystals(), edm::for_all(), fpeak_convoluted(), 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(), FWFFService::FWFFService(), 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(), TrackAssociatorByChi2::getChi2(), SCRegressionCalculator< VarCalc >::getCorrectionWithErrors(), HitPairGeneratorFromLayerPairForPhotonConversion::getCot(), RegressionHelper::getEcalRegression(), TouchableToHistory::getFilteredView(), HoECalculator::getHoE(), WatcherStreamFileReader::getInputFile(), PhiSymmetryCalibration::getKfactors(), EcalLaserDbService::getLaserCorrection(), EcalTrivialConditionRetriever::getMappingFromConfiguration(), FWMagField::GetMaxFieldMag(), EcalMappingElectronicsHandler::getNewObjects(), popcon::EcalSRPHandler::getNewObjects(), popcon::EcalTPGFineGrainEBIdMapHandler::getNewObjects(), ConversionProducer::getTrackImpactPosition(), EcalClusterLocalContCorrection::getValue(), EcalClusterCrackCorrection::getValue(), CSCAnalogSignal::getValue(), GlobalTrajectoryParameters::GlobalTrajectoryParameters(), CmsAnnotation::Handle(), FWColorPopup::HandleButton(), HcalEEUSMonitor::HcalEEUSMonitor(), HelixBarrelPlaneCrossing2OrderLocal::HelixBarrelPlaneCrossing2OrderLocal(), HitPairGeneratorFromLayerPairForPhotonConversion::hitPairs(), PixelTripletLargeTipGenerator::hitTriplets(), PixelTripletHLTGenerator::hitTriplets(), 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(), Bounds::inside(), PhotonConversionTrajectorySeedProducerFromSingleLegAlgo::inspectTrack(), funct::integral(), LinearGridInterpolator3D::interpolate(), SiStripTemplate::interpolate(), SiPixelTemplate::interpolate(), npstat::LinInterpolatedTableND< Numeric, Axis >::invertRatioResponse(), OAEParametrizedMagneticField::isDefined(), ParabolicParametrizedMagneticField::isDefined(), edm::isTransientEqual(), PerigeeConversions::jacobianParameters2Cartesian(), PerigeeConversions::jacobianPerigee2Curvilinear(), FWGUIValidatingTextEntry::keyPressedInPopup(), EvolutionECAL::LightCollectionEfficiency(), EvolutionECAL::LightCollectionEfficiencyWeighted(), GaussianSumUtilities1D::lnPdf(), HcalDetDiagPedestalMonitor::LoadDataset(), HcalDetDiagLaserMonitor::LoadDataset(), FWFFLooper::loadDefaultGeometryFile(), FFTJetPileupProcessor::loadFlatteningFactors(), MatacqProducer::loadOrbitOffset(), HcalPatternSource::loadPatternFile(), ClusterShapeHitFilter::loadPixelLimits(), HcalDetDiagLEDMonitor::LoadReference(), HcalDetDiagPedestalMonitor::LoadReference(), HcalDetDiagLaserMonitor::LoadReference(), PSimHit::localDirection(), PixelCPETemplateReco::localError(), PixelCPEGeneric::localError(), CSCRadialStripTopology::localError(), StripCPEfromTrackAngle::localParameters(), StripCPEfromTemplate::localParameters(), StripCPE::localParameters(), CSCRadialStripTopology::localPitch(), CSCRadialStripTopology::localStripLength(), GlobalTrajectoryParameters::magneticFieldInInverseGeV(), main(), MuonRPCDetLayerGeometryBuilder::makeBarrelLayers(), MuonRPCDetLayerGeometryBuilder::makeBarrelRods(), fireworks_root_gui::makeHorizontalFrame(), edm::EventSelector::maskTriggerResults(), funct::master(), CSCRadialStripTopology::measurementError(), JetPartonNamespace::UserPartonMixture::mixt(), MuonRoadTrajectoryBuilder::MuonRoadTrajectoryBuilder(), HcalRecHitsMaker::noiseInfCfromDB(), edm::AssociationMap< edm::OneToMany< reco::BasicJetCollection, reco::TrackCollection > >::numberOfAssociations(), LRHelpFunctions::obsFitIncluded(), CaloSamples::offsetTime(), oldComputeElectrons(), GsiFTPStorageMaker::open(), HttpStorageMaker::open(), XrdStorageMaker::open(), LocalStorageMaker::open(), DCacheStorageMaker::open(), StormStorageMaker::open(), StormLcgGtStorageMaker::open(), RFIOStorageMaker::open(), FreeFunctionSelector< T, f >::operator()(), OuterHitPhiPrediction::operator()(), BSpdfsFcn::operator()(), CandMatcherBase< C1, C2 >::operator()(), funct::GaussLegendreIntegrator::operator()(), funct::GaussIntegrator::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::AssociationMap< edm::OneToMany< reco::BasicJetCollection, reco::TrackCollection > >::operator[](), edm::JobReport::outputFileClosed(), PhiSymmetryCalibration_step2::outResidHistos(), PhiSymmetryCalibration_step2_SM::outResidHistos(), npstat::BoxND< Numeric >::overlapFraction(), GenParticlePruner::parse(), parseFFTJetScaleCalculator(), EcalTrigTowerConstituentsMapBuilder::parseTextMap(), PetrukhinFunc(), PFResolutionMapManager::PFResolutionMapManager(), ThirdHitPredictionFromCircle::phi(), CSCSegAlgoSK::phiAtZ(), CSCSegAlgoTC::phiAtZ(), InnerDeltaPhi::phiRange(), SiPixelTemplateReco::PixelTempReco2D(), SiPixelTemplateSplit::PixelTempSplit(), cond::PayLoadInspector< DataT >::plot(), plotMaker(), FWRPZView::populateController(), FWTriggerTableView::populateController(), FWFFService::postBeginJob(), VariableHelperService::postProcess(), InputTagDistributorService::preModule(), VariableHelperService::preModule(), fireworks::prepareTrack(), HcalRawDataMonitor::processEvent(), HcalRecHitMonitor::processEvent_rechit(), reco::modules::MultiTrackSelector::processMVA(), MCTrackMatcher::produce(), L1HLTJetsMatching::produce(), L1HLTTauMatching::produce(), PFJetsMatchedToFilteredCaloJetsProducer::produce(), GenTrackMatcher::produce(), BeamSpotOnlineProducer::produce(), L3TkMuonProducer::produce(), MatchedProbeMaker< T >::produce(), EcalTrivialConditionRetriever::produceEcalAlignmentEB(), EcalTrivialConditionRetriever::produceEcalAlignmentEE(), EcalTrivialConditionRetriever::produceEcalAlignmentES(), perftools::EdmEventSize::produceHistos(), tauImpactParameter::ErrorMatrixPropagator::propagateError(), gen::PtYDistributor::PtYDistributor(), EcalUncalibRecHitRecAnalFitAlgo< EBDataFrame >::pulseShapeFunction(), Histos::put(), SiStripTemplate::qbin(), SiPixelTemplate::qbin(), SiPixelTemplate::qbin_dist(), 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(), MuonResiduals6DOFFitter::readNtuple(), MuonResiduals6DOFrphiFitter::readNtuple(), PasswordReader::readPassword(), TkConvValidator::recalculateMomentumAtFittedVertex(), HcalSimpleRecAlgo::reconstruct(), HcalSimpleRecAlgo::reconstructHFUpgrade(), CmsShowNavigator::removeFilter(), PFEGammaAlgo::removeOrLinkECALClustersToKFTracks(), edm::JobReport::reportDataType(), edm::JobReport::reportFastCopyingStatus(), HcalRawDataMonitor::reset(), ViewerParameterGUI::reset(), CastorSimpleRecAlgo::resetTimeSamples(), EvolutionECAL::ResolutionConstantTermEM50GeV(), LA_Filler_Fitter::result(), fit::RootMinuit< Function >::RootMinuit(), LocalError::rotate(), LocalError::rotateCosSin(), PixelFitterByHelixProjections::run(), DQMStore::save(), TkHistoMap::saveAsCanvas(), npstat::ArrayND< Numeric, StackLen, StackDim >::scanInPlace(), FWViewContextMenuHandlerGL::select(), reco::modules::MultiTrackSelector::select(), edm::EventSelector::selectionDecision(), BaselinePFSCRegression::set(), VariableHelperService::set(), FWGeometryTableManagerBase::NodeInfo::setBit(), TrackInformation::setCaloIDChecked(), CmsShowMainBase::setConfigFilename(), UnbinnedLikelihoodFit::setFunction(), PFECALSuperClusterAlgo::setMajorityFraction(), hcalCalib::SetMaxProbeJetEmFrac(), hcalCalib::SetMaxTagJetEmFrac(), BscG4Hit::setPhiAtEntry(), FP420G4Hit::setPhiAtEntry(), TotemG4Hit::setPhiAtEntry(), BaseParticlePropagator::setPropagationConditions(), FWMET3DProxyBuilder::Arrow::setScale(), HcalPulseShape::setShapeBin(), CastorPulseShapes::Shape::setShapeBin(), HcalEEUSMonitor::setup(), HcalRawDataMonitor::setup(), FWEveView::setupEnergyScale(), sicif(), PhotonFix::sigmaEnergy(), SiPixelTemplate::simpletemplate2D(), XHistogram::splitSegment(), FWFFLooper::startingNewLoop(), stringToNumber(), cscdqm::stringToNumber(), SiStripTemplateReco::StripTempReco1D(), SiStripTemplateSplit::StripTempSplit(), FWGUIManager::subviewIsBeingDestroyed(), MSLayersAtAngle::sum2RmRn(), SummaryPlotXmlParser::SummaryPlotXmlParser(), MSLayersAtAngle::sumX0D(), MSLayer::sumX0D(), SiStripTemplate::sxtemp(), SiPixelTemplate::temperrors(), FWGeometryTableManagerBase::NodeInfo::testBit(), edm::EventSelector::testSelectionOverlap(), root::tf1(), TFileAdaptor::TFileAdaptor(), GlobalTrajectoryParameters::transverseCurvature(), funct::trapezoid_integral(), PrimaryVertexAnalyzer4PU::truthMatchedTrack(), reco::typeCode(), PFEGammaAlgo::unlinkRefinableObjectKFandECALMatchedToHCAL(), PFEGammaAlgo::unlinkRefinableObjectKFandECALWithBadEoverP(), HcalUnpacker::unpack(), unsafe_expf_impl(), spu::Untar(), HcalRawDataMonitor::UpdateMEs(), FWPFRhoPhiRecHit::updateScale(), root::helper::RootVarsAdapter< F, 1 >::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(), sistrip::MeasureLA::write_report_text_ms(), edm::JobReport::JobReportImpl::writeInputFile(), VirtualJetProducer::writeJets(), edm::JobReport::JobReportImpl::writeOutputFile(), MSLayer::x0(), SiStripTemplate::xflcorr(), SiPixelTemplate::xflcorr(), geometryXMLparser.MuonGeometry::xml(), SiStripTemplate::xsigma2(), SiPixelTemplate::xsigma2(), SiPixelTemplate2D::xysigma2(), SiPixelTemplate2D::xytemp(), SiPixelTemplate::yflcorr(), SiPixelTemplate::ysigma2(), and edm::service::InitRootHandlers::~InitRootHandlers().

double g[11][100]

Definition at line 78 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 82 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 86 of file MuScleFitUtils.cc.

double isum

Definition at line 76 of file MuScleFitUtils.cc.

Referenced by SiPixelActionExecutor::fillFEDErrorSummary(), SiStripSummaryCreator::fillGrandSummaryHistos(), SiPixelActionExecutor::fillSummary(), SiStripSummaryCreator::fillSummaryHistos(), and SiStripSummaryCreator::setSummaryMENames().

double mzsum

Definition at line 75 of file MuScleFitUtils.cc.