pimpCOLZ.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////
//                                                            //
//    SECOND ANALYSIS                                         //
//                                                            //
//    Detector response                                       //
//    TH2 COLZ plots PECorr/cm vs W vs Cell                   //
//                                                            //
//    Created by: Diana Patricia Mendez                       //
//    10 March 2016                                           //
//                                                            //
//    Input files generated with COLZ.cxx                     //
////////////////////////////////////////////////////////////////

#include <iostream>
#include <TFile.h>
#include <TPad.h>
#include <TProfile.h>
#include <TProfile2D.h>
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TH3D.h>
#include <TCanvas.h>
#include <TVector.h>
#include <TColor.h>
#include <TLegend.h>
#include <TLatex.h>
#include <TMath.h>
#include <TStyle.h>
#include <TROOT.h>
#include <TPad.h>
#include <TFile.h>
#include <TTree.h>
#include <TSystem.h>
#include <TColor.h>

#include "Utilities/rootlogon.C"

// canvas size
#define c1_x 700
#define c1_y 500 

#define c2_x 1000
#define c2_y 500


double GetHistMax(TH1* h0, TH1* h1, TH1* h2, TH1* h3){
  double max[4];
  max[0] = h0->GetMaximum();
  max[1] = h1->GetMaximum();
  max[2] = h2->GetMaximum();
  max[3] = h3->GetMaximum();

  double hiy = 0;
  for(int i=1;i<4;i++){
    if(max[i]>hiy)
      hiy = max[i];
  }
  return (hiy * 1.1);
}


int main(){

  gROOT->ForceStyle();


  //Normalization constants = MEUs for each set
  double const fdmcg100  = 1;
  // double const fdmcg100  = 40.2581;
  
  TFile *fmc_gain100 = new TFile("FD_WCellResponse_mcg100.root");
  TH2D *hfdmc_gain100_PECorrWX1 = (TH2D*)fmc_gain100->Get("hPECorrCellPlaneX");
  TH2D *hfdmc_gain100_PECorrWY1 = (TH2D*)fmc_gain100->Get("hPECorrCellPlaneY");


  
  ////////////////////////////////////////////
  //   W vs Cell vs Response                //
  //   PROF COLZ Plots                      //
  ////////////////////////////////////////////
  
 
  TCanvas *cFDMC_PlaneCellX_mcg100 = new TCanvas("cFDMC_PlaneCellX","cFDMC_PlaneCellX", 700, 500);
  cFDMC_PlaneCellX_mcg100->SetRightMargin(0.15);
  hfdmc_gain100_PECorrWX1->GetXaxis()->SetTitle("W, x-view (cm)");
  hfdmc_gain100_PECorrWX1->GetYaxis()->SetTitle("Cell");
  hfdmc_gain100_PECorrWX1->GetZaxis()->SetRangeUser(0,80);
  hfdmc_gain100_PECorrWX1->GetZaxis()->SetTitle("PECorr/cm");
  CenterTitles(hfdmc_gain100_PECorrWX1);
  hfdmc_gain100_PECorrWX1->Draw("colz");
  //  NovaFD();
  Simulation();

  cFDMC_PlaneCellX_mcg100 -> Print("cFDmcgain100_WCellPECorrX.pdf");
  cFDMC_PlaneCellX_mcg100 -> Print("cFDmcgain100_WCellPECorrX.png");
  
}

int respCOLZ_WCell(){
  return main();
}