Xephyr
0.0
|
implements 2D pdf profiled likelihood fitting. More...
#include <XeLikelihoods.h>
Public Member Functions | |
pdfLikelihood (TString name, double wimpMass) | |
void | initialize () |
void | setData (int dataType) |
double | computeTheLogLikelihood () |
void | generateAsimov (double mu_prime) |
void | generateAndUseAsimov (double mu_prime) |
void | generateToyDataset (double seed, double mu_prime) |
double | getSignalDefaultNorm () |
double | getSignalMultiplier () |
void | setSignalMultiplier (double val) |
void | setTreeIndex (int index) |
void | setSignalDefaultNorm (double val) |
double | getCurrentNs () |
void | checkInputs () |
void | addBkgPdfComponent (pdfComponent *addMe, bool Safeguarded=false) |
void | setSignalPdf (pdfComponent *signal) |
void | setDataHandler (dataHandler *d) |
void | setAsimovData (dataHandler *d) |
int | useAsimovData () |
int | useDMData () |
bool | isNegativeAnywhere (TH2F histo) |
void | setCalibrationData (dataHandler *d) |
void | setWithSafeGuard (bool b) |
void | setSafeGuardPosDef (bool b) |
Set whether to use safeguard or not. | |
void | drawAllOnProjection (bool isS1Projection) |
Set whether safeguard should be forced to be possitive. | |
void | printEventSummary (bool isForWiki=false) |
prints a summary of all bkg and signal events with current parameter choice More... | |
vector< string > | getTrueParamsNames () |
vector< double > | getTrueParams () |
double | getSafeguardValue () |
pdfComponent * | getBkgComponent (TString name) |
returns current value of safeguard More... | |
vector< pdfComponent * > | getBkgComponents () |
Returns vector of all background components. | |
histoCompare | getModelCompare () |
Get a handler for plotting comparisons between models. | |
TH2F | getOverallBkgHisto () |
TH2F | getSafeguardedBkgPdf () |
TH2F | getSafeguardedBkgPdfOnly () |
void | setAdditionalSafeGuardComponent (TH2F *h) |
void | setFixedValueForSafeguard (double fv) |
double | LLsafeGuard () |
int | numberOfSafeguarded () |
returns the number of bkg components that ask for safeguard | |
double | getWimpMass () |
![]() | |
ProfileLikelihood (TString name) | |
void | printFlagsAndParameters () |
void | estimateCrossSection () |
bool | initialize () |
TGraph * | getGraphOfLogLikelihood (int n_points) |
TGraph * | getGraphAtQval (double qval, bool forceMuhatComp=true) |
double | getSigmaAtQval (double qval, bool forceMuhatComp=true) |
double | returnLimitHagar (double cl, bool forceMuhatComp=true) |
vector< TGraph * > | getGraphOfParameters (int n_points) |
TGraph * | getLikelihoodScanOfParameter (int n_points, LKParameter *par, double mu=0) |
![]() | |
double | computeTheConstraint () |
Likelihood (TString name) | |
double | getSigmaHat () |
double | getLogD () |
bool | parameterExists (int p) |
void | addParameter (int id, int type, TString nam, double initialVal, double step, double mi, double ma) |
void | addParameter (LKParameter *param, int id=SAME) |
void | replaceParameter (LKParameter *param) |
void | addParameterTolerant (LKParameter *param) |
void | activateParameter (LKParameter *param, bool active=true) |
void | removeParameter (LKParameter *param, bool tolerant=false) |
void | removeParameter (int id, bool tolerant=false) |
void | listParameters () |
void | printInitialParameters () |
void | printResultParameters () |
void | printCurrentParameters () |
void | printCurrentParameters (bool with_err) |
void | ignoreParameter (int id) |
void | setParameterType (int id, int type) |
void | setParameterValue (int id, double v) |
void | setInitialValue (int id, double v) |
double | getParameterValue (int id) |
double | maximize (bool freezeParametersOfInterest) |
double | maximizeNumerically (int numberOfToys, bool freezeParametersOfInterest) |
void | setSeed (double Inputseed) |
int | getNTotalParameters () |
int | getNActiveParameters () |
int | getNParametersOfInterest () |
int | getNMinuitParameters () |
int | getNNuisanceParameters () |
int | getNParameters (int type) |
void | resetParameters () |
LKParameter * | getParameter (int id) |
LKParameter * | getPOI () |
map< int, LKParameter * > * | getParameters () |
get Parameter Of Interest | |
TH1F | getPullsHisto () |
void | setCurrentValues (const double *v, const double *ers=NULL) |
void | setCurrentValuesInMinuitUnits (const double *v, const double *e=NULL) |
void | setCombinedMode (bool mode=true) |
void | printInitialHeader () |
void | printCurrentHeader () |
bool | checkConsistency () |
bool | inCombinedMode () |
int | mapMinuitParameters (bool freezeParametersOfInterest) |
int | getNParametersForChi2 () |
void | forceNParametersOfInterest (int nF) |
void | clearTheParameters () |
![]() | |
XeStat (TString name) | |
void | setName (TString nam="No_Name") |
helper method to set the name of this instance. | |
TString | getName () |
helper method to return the name of this instance. | |
void | setExperiment (int exp) |
set the experiment number, usefull for combination | |
int | getExperiment () |
set the experiment number, usefull for combination | |
![]() | |
errorHandler (TString name) | |
void | Error (TString functionName, TString message) |
Send an Error message to the user. More... | |
void | Warning (TString functionName, TString message) |
Send a Warning message to the user. More... | |
void | Info (TString functionName, TString message) |
Send an Info message to the user. More... | |
void | Debug (TString functionName, TString message) |
Send a Debug message to the user. More... | |
void | setPrintLevel (int level) |
Set the local printing level (for this instance only) More... | |
int | getPrintLevel () |
Return the current print level for this instance. remind that if you set with setPrintLevel this will override the global print level. | |
![]() | |
printTools () | |
---------— Print Toools -----------— /// | |
Public Attributes | |
vector< pdfComponent * > | bkg_components |
vector< bool > | safeguarded_bkg_components |
says if a bkg component needs to be safeguarded. | |
TH2F * | safeguardAdditionalComponent |
is an additional component to be added in the fit to calibration for | |
pdfComponent * | signal_component |
dataHandler * | data |
dataHandler * | calibrationData |
dataHandler * | asimovData |
dataHandler * | dmData |
double | siganlDefaultNorm |
double | XsecMultiplier |
bool | withSafeGuard |
bool | safeGuardPosDef |
double | wimp_mass |
Force safeguard parameter to be positive. | |
double | safeguard_fixValue |
LKParameter * | safeGuardParam |
double | safeguard_scaling |
![]() | |
double | seed |
![]() | |
int | localPrintLevel |
TString | className |
Additional Inherited Members | |
![]() | |
static double | shapeLikelihood (vector< int > *bins, int nDist, double **dists, double *norm) |
![]() | |
static void | setAnalysisMode (int mode) |
static int | getAnalysisMode () |
static bool | checkAnalysisMode (TString name, int requested) |
static bool | isAnalysisDefined (bool verbose=false) |
static bool | isCutsBased () |
static bool | isPL () |
static int | getSigmaLinLog (int unit) |
static TString | getAnalysisModeName () |
static TString | getAnalysisModeName (int mode) |
static TString | getSystematicModeName (int syst) |
static TString | getSigmaUnitName (int unit) |
static TString | getSigmaLabel (int unit) |
static TString | getSigmaModeName (int mode) |
static TString | getUpperSigmaLabel (int unit) |
![]() | |
static TString | upperCase (TString s) |
static TString | lowerCase (TString s) |
static TString | trim (TString s) |
static TString | justify (TString s, int w, bool left, bool trim) |
static TString | rightJustify (TString s, int w, bool trim=false) |
static TString | leftJustify (TString s, int w, bool trim=false) |
static TString | format0I (int v, int w=1) |
static TString | formatLI (int v, int w=1, int p=1) |
static TString | formatRI (int v, int w=1, int p=1) |
static TString | formatF (double v, int w=1, int p=1) |
static TString | formatR (double v, int w=1, int p=1) |
static TString | formatLF (double v, int w=1, int p=1) |
static TString | formatRF (double v, int w=1, int p=1) |
static TString | formatG (double v, int w=1, int p=1) |
static TString | formatLG (double v, int w=1, int p=1) |
static TString | formatRG (double v, int w=1, int p=1) |
static TString | formatI (int v, int w=1, int trailer=0) |
static TString | doOrDont (bool b) |
![]() | |
static int | globalPrintLevel = 1 |
![]() | |
void | setup () |
![]() | |
void | setup () |
void | clear () |
bool | checkParameter (int p, bool shouldExist) |
![]() | |
LKParameter * | sigPar |
![]() | |
int | currentId |
int | nParametersOfInterest |
int | nActiveParameters |
int | nNuisanceParameters |
int | forcedNParametersOfInterest |
bool | combinedMode |
vector< LKParameter * > | MinuitParameters |
map< int, LKParameter * > | parameters |
double | sigmaHat |
double | LogD |
![]() | |
TString | name |
int | experiment |
![]() | |
static int | analysisMode = NO_ANALYSIS |
implements 2D pdf profiled likelihood fitting.
It expects as inputs pdfObjects. It also implements safeguard.
|
virtual |
none | All parameters values set thru LKParameter |
Implements Likelihood.
pdfComponent * pdfLikelihood::getBkgComponent | ( | TString | name | ) |
returns current value of safeguard
Returns pointer to background component by name, first match is returned.
TH2F pdfLikelihood::getOverallBkgHisto | ( | ) |
returns the current overall background hostogram (with all background added up). Current means that is for the current set of parameters, so each bkg model is interpolated, scaled and then added to the others.
TH2F pdfLikelihood::getSafeguardedBkgPdf | ( | ) |
returns (1-epsilon)Fb(x,y) + epsilon*Fs(x,y) with all uncertainties interpolated and added. this term is a full bkg pdf including also the non safeguarded bkgs TO BE USED IN THE "physics" likelihood
TH2F pdfLikelihood::getSafeguardedBkgPdfOnly | ( | ) |
returns (1-epsilon)Fb(x,y) + epsilon*Fs(x,y) for all bkg components that are considered for dafeguard, with all uncertainties interpolated and added. TO BE USED in the fit to calibration
|
virtual |
prints a summary of all bkg and signal events with current parameter choice
Implements ProfileLikelihood.