#include <TROOT.h>
#include <TSystem.h>
#include <TChain.h>
#include <TFile.h>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <stdlib.h>
#include <sstream>
#include <stdexcept>
#include <TH2.h>
#include <TF1.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <TFile.h>
#include <TObjArray.h>
#include "TLegend.h"
#include "TGraph.h"
#include "stdio.h"
#include <TLatex.h>
#include <TColor.h>
#include <TMarker.h>
#include <TString.h>
#include "TGraphAsymmErrors.h"
#include <vector>
#include <TRandom3.h>
#include <TGaxis.h>
#include <algorithm>
#include <array>


void MyScript(){

TFile* file1 = TFile::Open("1512.03657_7TeV_histoFit_Function1.root"); 

//Canvas: 
TCanvas* c1 = new TCanvas("c1","",800,800);


//Get histograms:
TH1F * h1 = (TH1F*)file1->Get("CombinedJPSIcutPT");
TF1 * f1= (TF1*)h1->GetFunction("fitJPSI_Combined");
TF1 * f1b = (TF1*)f1->Clone("f1b");

f1b->SetParNames("A1","A2","A3","A4","A5");

TH1F * h2 = (TH1F*)file1->Get("CombinedPSI2scutPT");
TF1 * f2 = (TF1*)h2->GetFunction("fitPSI2s_Combined");
TF1 * f2b = (TF1*)f2->Clone("f2b");

f2b->SetParNames("B1","B2","B3","B4","B5");


//Plot f1b multplied by f2b: 
TF1 *Multiply = new TF1("Multiply","f1b*f2b");



c1->cd();
Multiply->Draw();



}