XeLikelihoods.h

#ifndef XELIKELIHOODS_H
#define XELIKELIHOODS_H

#include "XeStat.h"
#include <TString.h>
#include "TGraph2D.h"
#include "TNtuple.h"
#include <cmath>
#include <csignal>
#include <iostream>
#include <cstdio>
#include <vector>
#include <utility>      // std::pair, std::make_pair
#include "dataHandler.h"

#include "XePdfObjects.h"
#include "plotHelpers.h"

using namespace std;

class pdfLikelihood : public ProfileLikelihood {

  public :

    //FIXME  TODO put some comments

    ~pdfLikelihood();
    pdfLikelihood(TString name, double wimpMass);


    void initialize();

    void setData(int dataType);

    double computeTheLogLikelihood();

    void generateAsimov(double mu_prime) ;

    void generateAndUseAsimov(double mu_prime) {generateAsimov(mu_prime); useAsimovData() ;}

    void generateToyDataset(double seed, double mu_prime);

    double getSignalDefaultNorm() { return siganlDefaultNorm; } ;

    double getSignalMultiplier()  { return XsecMultiplier;    } ;

    void   setSignalMultiplier(double val) { XsecMultiplier = val; } ;

    void   setTreeIndex(int index);

    void   setSignalDefaultNorm(double val) { siganlDefaultNorm = val; } ;

        double getCurrentNs();

    void   checkInputs();

    // add a bkg pdf, if its shape should be considered in safeguard method
    // Safeguarded should be true, default is false.
    void addBkgPdfComponent(pdfComponent *addMe , bool Safeguarded = false);

    void setSignalPdf(pdfComponent *signal)  {signal_component = signal; };

    void setDataHandler(dataHandler *d )     {dmData = d; useDMData(); } ;

    void setAsimovData(dataHandler *d) {asimovData = d;  useAsimovData(); } ;

    int useAsimovData();

    int useDMData();

    bool isNegativeAnywhere(TH2F histo);

    void setCalibrationData(dataHandler *d )     {calibrationData = d; } ;

    void setWithSafeGuard(bool b)  {withSafeGuard = b;} ; 

  void setSafeGuardPosDef(bool b)  {safeGuardPosDef = b;} ; 

  void drawAllOnProjection(bool isS1Projection);

    void printEventSummary(bool isForWiki=false);

    vector<string> getTrueParamsNames() { return data->getTrueParamsNames(); };

    vector<double> getTrueParams()      { return data->getTrueParams(); };

  double getSafeguardValue(); 


    pdfComponent* getBkgComponent(TString name);

  vector<pdfComponent *> getBkgComponents();

    // used by drawAllOnProjection()
    histoCompare getModelCompare();

    TH2F getOverallBkgHisto();

    TH2F getSafeguardedBkgPdf();

    TH2F getSafeguardedBkgPdfOnly();

    vector <pdfComponent*> bkg_components;

    vector <bool>    safeguarded_bkg_components;

    // safeguard. Used in this case for the AC that is different between data and Rn
    // this component has to be scaled to Bkg_data ---> Nbkg/Ncal
    TH2F *safeguardAdditionalComponent;

    void setAdditionalSafeGuardComponent(TH2F *h){ safeguardAdditionalComponent = h;};

    void setFixedValueForSafeguard(double fv) { safeguard_fixValue = fv; } ;

    double LLsafeGuard();

    int   numberOfSafeguarded();

    pdfComponent           *signal_component;

    dataHandler             *data;

    dataHandler            *calibrationData;

    dataHandler            *asimovData;

    dataHandler            *dmData;

    double                 siganlDefaultNorm;   //corresponding cross section of histos

    double                 XsecMultiplier;   // handy multiplier to keep sigma [0,50]

    bool                   withSafeGuard;

  bool                   safeGuardPosDef; 

    double                 wimp_mass;

    double                 safeguard_fixValue;

    LKParameter            *safeGuardParam;

    double                  safeguard_scaling;

    //This is needed for compatibility, ancestral xephyr roots.
    //FIXME: move getWimpMass to Asymptotics
    double getWimpMass() {return wimp_mass;};
     //------------------------------//



};





#endif