append.C File Reference

#include "TROOT.h"
#include "TObject.h"
#include "TH2.h"
#include "TH1.h"
#include "TMatrix.h"
#include "TMatrixD.h"
#include "TFile.h"
#include <iostream>
#include <algorithm>
#include <string>

Go to the source code of this file.

Functions
void	append ()

Function Documentation

void append

(

)

Definition at line 24 of file append.C.

References MECModelEnuComparisons::i, inFile, calib::j, string, and febshutoff_auto::val.

Referenced by pclistAnalyser::add_to_children(), raw2rootAnalyser::analyse(), pclistAnalyser::analyse(), lemmon_tools::appendStatus(), boost::python::vector_indexing_suite< Container, NoProxy, final_vector_derived_policies< Container, NoProxy > >::base_append(), runNovaSAM::checkAndMoveFiles(), ProductionTestTools::crontabToGantt(), log_trawl.Project::decomposeModes(), ProductionTestTools::directoryToGantt(), ControllerView.ControllerView::draw_component_plots(), bp::FindChanges(), bp::FindFiles(), ProductionTestTools::findMissingFiles(), ViewHTML::foldersToButtonGroups(), generic.tools.HistogramTools::getPlotStats(), common_tools::getPlotStats(), lem_server::handle_get_work(), lem_server::handle_pid(), log_trawl.Project::inspectProcs(), PandAna.utils.sadd_pandana::main(), FirstAnaPageMakers::makeDatasetPage(), PageMakers::makeResultsPage(), runNovaSAM::MakeTransposeJson(), PageMakers::makeValidationPage(), SupernovaUtilities::ParseSNDFile(), PandAna.core.core.loader::readData(), SupernovaUtilities::ReadGarchingLegacyFluxes(), SamTimer.SamTimer::record_time(), ViewHTML::samProjectCollectionToTable(), and validateGRLvsSAM::validate().

{

  // file where xsec is scaled up 
  std::string inName = 
  "<MY_INPUTS>/C_PiMinus_fit.root";

  TFile *inFile = new TFile(inName.c_str(), "UPDATE");

 
  TMatrixT<double> *FitCovariance = (TMatrixT<double>*)inFile->Get("FitCovariance");

  FitCovariance->ResizeTo(6,6);

  // Set diagonal entry to (sigma)**2, all the others stay as zero
  double cov_value = 0.2*0.2;

  (*FitCovariance)(4,4) = cov_value;
  (*FitCovariance)(5,5) = cov_value;   

  FitCovariance->Write("FitCovariance", TObject::kOverwrite);
    

  std::vector existing_ints{"fAbs", "fCex", "fInel", "fReac"};
  TH1D* parsHist = (TH1D*)inFile->Get("parsHist");
  TH1D new_parsHist("new_parsHist", "", 6, 0, 6);

  TH2D* FitCovHist = (TH2D*)inFile->Get("FitCovHist");
  TH2D new_FitCovHist("new_FitCovHist", "", 6, 0, 6, 6, 0, 6);


  for(uint i=0; i<existing_ints.size(); ++i){
    double val = parsHist->GetBinContent(i+1);
    double err = parsHist->GetBinError(i+1);
    //std::cout << i+1 << " = " << val << " +/- " << err << std::endl;
    new_parsHist.SetBinContent(i+1, val);
    new_parsHist.GetXaxis()->SetBinLabel( i+1, existing_ints.at(i));
    new_parsHist.SetBinError(i+1, err);

    new_FitCovHist.GetXaxis()->SetBinLabel( i+1, existing_ints.at(i));
    new_FitCovHist.GetYaxis()->SetBinLabel( i+1, existing_ints.at(i));
    
    for(uint j=0; j<existing_ints.size(); ++j){
      double cov_val = FitCovHist->GetBinContent(i+1, j+1);
      new_FitCovHist.SetBinContent(i+1, j+1, cov_val);
    }

  }

  new_parsHist.SetBinContent(5, 1.);
  new_parsHist.GetXaxis()->SetBinLabel( 5, "fDcex");
  new_parsHist.SetBinError(5, 0.2);

  new_parsHist.SetBinContent(6, 1.);
  new_parsHist.GetXaxis()->SetBinLabel( 6, "fProd");
  new_parsHist.SetBinError(6, 0.2);

  new_parsHist.Write("parsHist", TObject::kOverwrite);


  new_FitCovHist.SetBinContent(5, 5, cov_value);
  new_FitCovHist.GetXaxis()->SetBinLabel(5, "fDcex");
  new_FitCovHist.GetYaxis()->SetBinLabel(5, "fDcex");

  new_FitCovHist.SetBinContent(6, 6, cov_value);
  new_FitCovHist.GetXaxis()->SetBinLabel(6, "fProd");
  new_FitCovHist.GetYaxis()->SetBinLabel(6, "fProd");

  new_FitCovHist.Write("FitCovHist", TObject::kOverwrite);


  inFile->Close();

}