void merge()
{
    
    TFile* f1 = new TFile("RBBhistogram.root");
    
    TH1::AddDirectory(kFALSE) ;
    
    
    TH1D* hEdepTheta0RBB = (TH1D*)f1->Clone("EdepTheta0");
    hEdepTheta0RBB -> Scale(0.0005);

    TH1D* hEdepTheta68RBB = (TH1D*)f1->Clone("EdepTheta68");
    hEdepTheta68RBB -> Scale(0.0005);

    TH2D* hEdepPos0RBB = (TH2D*)f1->Clone("EdepPosTheta0");
    hEdepPos0RBB -> Scale(0.0005);
    
    TH2D* hEdepPos68RBB = (TH2D*)f1->Clone("EdepPosTheta68");
    hEdepPos68RBB -> Scale(0.0005);

    TH1D* hEdepCell0RBB = (TH1D*)f1->Clone("EdepPerRelevantCell1");
    hEdepCell0RBB -> Scale(0.0005);
    
    TH1D* hEdepCell68RBB = (TH1D*)f1->Clone("EdepPerRelevantCell2");
    hEdepCell68RBB -> Scale(0.0005);

    TH1D* hNumHits0RBB = (TH1D*)f1->Clone("NumOfHitsRelevantCell1");
    hNumHits0RBB -> Scale(0.0005);

    TH1D* hNumHits68RBB = (TH1D*)f1->Clone("NumOfHitsRelevantCell2");
    hNumHits68RBB -> Scale(0.0005);

    TH2D* hEdepPosCell0RBB = (TH2D*)f1->Clone("EdepPostitionRelevantCell1");
    hEdepPosCell0RBB -> Scale(0.0005);
    

    TH2D* hEdepPosCell68RBB = (TH2D*)f1->Clone("EdepPostitionRelevantCell2");
    hEdepPosCell68RBB -> Scale(0.0005);
 
    f1->Close();

    TFile* f2 = new TFile("AllButRBBhistogram.root");

    TH1D* hEdepTheta0All = (TH1D*)f2->Clone("EdepTheta0");
    hEdepTheta0All -> Scale(0.00005);
    TH1D* hEdepTheta68All = (TH1D*)f2->Clone("EdepTheta68");
    hEdepTheta68All -> Scale(0.00005);
    TH2D* hEdepPos0All = (TH2D*)f2->Clone("EdepPosTheta0");
    hEdepPos0All -> Scale(0.00005);
    TH2D* hEdepPos68All = (TH2D*)f2->Clone("EdepPosTheta68");
    hEdepPos68All -> Scale(0.00005);
    TH1D* hEdepCell0All = (TH1D*)f2->Clone("EdepPerRelevantCell1");
    hEdepCell0All -> Scale(0.00005);
    TH1D* hEdepCell68All = (TH1D*)f2->Clone("EdepPerRelevantCell2");
    hEdepCell68All -> Scale(0.00005);
    TH1D* hNumHits0All = (TH1D*)f2->Clone("NumOfHitsRelevantCell1");
    hNumHits0All -> Scale(0.00005);
    TH1D* hNumHits68All = (TH1D*)f2->Clone("NumOfHitsRelevantCell2");
    hNumHits68All -> Scale(0.00005);
    

    TH2D* hEdepPosCell0All = (TH2D*)f2->Clone("EdepPostitionRelevantCell1");
    hEdepPosCell0All -> Scale(0.00005);
    TH2D* hEdepPosCell68All = (TH2D*)f2->Clone("EdepPostitionRelevantCell2");
    hEdepPosCell68All -> Scale(0.00005);
   
    f2->Close();
    
    TH1D* hEdepTheta0 = new TH1D();
    TH1D* hEdepTheta68= new TH1D();
    TH2D* hEdepPos0= new TH2D();
    TH2D* hEdepPos68= new TH2D();
    TH1D* hEdepCell0= new TH1D();
    TH1D* hEdepCell68= new TH1D();
    TH1D* hNumHits0= new TH1D();
    TH1D* hNumHits68= new TH1D();
    TH2D* hEdepPosCell0= new TH2D();
    TH2D* hEdepPosCell68= new TH2D();

    hEdepTheta0 ->Add(hEdepTheta0RBB,hEdepTheta0All);

    hEdepTheta68 -> Add(hEdepTheta68RBB, hEdepTheta68All);
    hEdepPos0 -> Add(hEdepPos0RBB, hEdepPos0All);
    hEdepPos68 -> Add(hEdepPos68RBB, hEdepPos68All);
    hEdepCell0 -> Add(hEdepCell0RBB, hEdepCell0All);
    hEdepCell68 -> Add(hEdepCell68RBB, hEdepCell68All, 1, 1);
    hNumHits0 -> Add(hNumHits0RBB, hNumHits0All, 1 ,1);
    hNumHits68 -> Add(hNumHits68RBB, hNumHits68All, 1 ,1);
    hEdepPosCell0 -> Add(hEdepPosCell0RBB, hEdepPosCell0All, 1 ,1);
    hEdepPosCell68 -> Add(hEdepPosCell68RBB, hEdepPosCell68All, 1, 1);
    

    

}