calibration/ATRI/makeCalibPlots.cxx
00001 //AraRoot Includes 00002 #include "UsefulAtriStationEvent.h" 00003 #include "RawAtriStationEvent.h" 00004 #include "araSoft.h" 00005 00006 00007 //Root Includes 00008 #include "TTree.h" 00009 #include "TFile.h" 00010 #include "TH1.h" 00011 #include "TTree.h" 00012 #include "TMath.h" 00013 #include "TCanvas.h" 00014 00015 #define USEFUL_RFCHANS 6 00016 00017 int makeCalibPlots(char *calibFileName){ 00018 00019 00020 TFile *inFile = TFile::Open(calibFileName); 00021 if(!inFile){ 00022 printf("Failed to open file %s\n", calibFileName); 00023 return -1; 00024 } 00025 TCanvas *canBinWidths[DDA_PER_ATRI]; 00026 for(int dda=0;dda<DDA_PER_ATRI;dda++){ 00027 00028 char canName[100]; 00029 sprintf(canName, "canBinWidths_DDA%i", dda); 00030 canBinWidths[dda] = new TCanvas(canName); 00031 canBinWidths[dda]->Divide(3,2); 00032 for(int chan=0;chan<USEFUL_RFCHANS;chan++){ 00033 for(int half=0;half<2;half++){ 00034 char histName[100]; 00035 sprintf(histName, "histBinWidth_dda%i_chan%i_%i", dda, chan, half); 00036 TH1* tempHist = (TH1*) inFile->Get(histName); 00037 canBinWidths[dda]->cd(chan+1); 00038 if(half){ 00039 tempHist->SetLineColor(kRed); 00040 tempHist->Draw("SAME"); 00041 } 00042 else{ 00043 tempHist->SetLineColor(kBlue); 00044 tempHist->Draw(); 00045 } 00046 00047 00048 } 00049 } 00050 } 00051 TCanvas *canInterleave[DDA_PER_ATRI]; 00052 for(int dda=0;dda<DDA_PER_ATRI;dda++){ 00053 char canName[100]; 00054 sprintf(canName, "canInterleave_DDA%i", dda); 00055 canInterleave[dda] = new TCanvas(canName); 00056 canInterleave[dda]->Divide(3,2); 00057 for(int chan=0;chan<USEFUL_RFCHANS;chan++){ 00058 char histName[100]; 00059 sprintf(histName, "lag_hist_dda%i_chan%i", dda, chan); 00060 TH1* tempHist = (TH1*) inFile->Get(histName); 00061 canInterleave[dda]->cd(chan+1); 00062 tempHist->Draw(); 00063 } 00064 } 00065 TCanvas *canEpsilon[DDA_PER_ATRI]; 00066 for(int dda=0;dda<DDA_PER_ATRI;dda++){ 00067 char canName[100]; 00068 sprintf(canName, "canEpsilon_DDA%i", dda); 00069 canEpsilon[dda] = new TCanvas(canName); 00070 canEpsilon[dda]->Divide(3,2); 00071 for(int chan=0;chan<USEFUL_RFCHANS;chan++){ 00072 canEpsilon[dda]->cd(chan+1); 00073 char histName[100]; 00074 sprintf(histName, "epsilon_hist_dda%i_chan%i", dda, chan); 00075 TH1* tempHist = (TH1*) inFile->Get(histName); 00076 tempHist->SetLineColor(kBlack); 00077 tempHist->Draw(""); 00078 } 00079 } 00080 00081 00082 return 0; 00083 00084 }
Generated on Tue Jul 16 16:58:02 2013 for ARA ROOT v3.10 Software by
1.4.7