File:Edep.jpg

// // Kyle Snavely June 30, 2008 // // This script reads in rootfiles and makes a plot of // Edep with different cuts // // -------------------------------------------------------- //

// All my standard ROOT includes...

1. include <TString.h>
2. include <string>
3. include <fstream>
4. include <iostream>
5. include <TCanvas.h>
6. include <TH1.h>
7. include <TH2.h>
8. include <TF1.h>
9. include <TH3F.h>
10. include <TFile.h>
11. include <TTree.h>
12. include <TROOT.h>
13. include <TStyle.h>
14. include <TApplication.h>
15. include <TPaveText.h>
16. include <TPaveStats.h>
17. include <TPaveLabel.h>
18. include <vector>
19. include <TProfile.h>
20. include <TF1.h>
21. include <TGraphErrors.h>
22. include <TMultiGraph.h>
23. include <TGraph.h>
24. include <TDatime.h>
25. include <TLegend.h>
26. include <TAxis.h>
27. include <TGaxis.h>
28. include <TPad.h>
29. include <TLatex.h>
30. include <TRandom.h>
31. include <stdio.h>
32. include <stdlib.h>
33. include <math.h>
34. include <TMath.h>
35. include <TMinuit.h>
36. include <TStopwatch.h>
37. include <TSpectrum.h>
38. include "TMath.h"

// Interactive includes

1. include <TRint.h>
2. include <TApplication.h>
3. include <TGX11.h>
4. include <TVirtualX.h>
5. include <TFriendElement.h>

using namespace std ;

// Make tree a global variable so that we can access it outside this function. // You can then also access it from the command line. TTree *tree;

// ************ The main program starts here. ****************//

int main(int argc, char** argv) {

 ///////////////////////////////////////////
 //   Initialize root, get our files and  //
 // make them friends, create our canvas. //
 ///////////////////////////////////////////
 TApplication theApp("App", &argc, argv);
 // gStyle->SetOptFit(0111);
 gROOT->SetStyle("Plain");
 gStyle->SetFillColor(10); 

 // Open Root file
 TFile *f = new TFile("$DATA/proton_500000.root");
 //Get bst tree
 TTree *bst = (TTree*)f->Get("bstT");
 bst->AddFriend("genT", "$DATA/proton_500000.root");

 TFile *g = new TFile("$DATA/pion_500000.root");
 TTree *bst_pi = (TTree*)g->Get("bstT");
 bst_pi->AddFriend("genT", "$DATA/pion_500000.root");
 
 // Make your canvas
 TCanvas* MyCanvas1 = new TCanvas("MyCanvas1", "Some Graphs", 1);
 MyCanvas1->SetFillColor(10);
 TPaveLabel *title1 = NULL;
 title1 = new TPaveLabel(0.1,0.96,0.9,0.99,"A Plot!  Yay!");
 title1->SetFillColor(10);
 //title1->Draw();
 TDatime *now1;
 now1 = new TDatime();
 TPaveLabel *date1 = new TPaveLabel(0.7,0.01,0.9,0.03,now1->AsString());
 date1->SetFillColor(10);
 //date1->Draw();
 //TPad* graphPad1 = new TPad("Graphs","Graphs",0.01,0.05,0.95,0.95);
 TPad* graphPad1 = new TPad("Graphs","Graphs",0.0,0.0,1.0,1.0);
 graphPad1->SetFillColor(10);
 graphPad1->Draw();
 graphPad1->cd();

 ///////////////////////////////////////////
 //   Create and fill our histos          //
 ///////////////////////////////////////////

 TH1F *Edep_h = new TH1F("Edep_h","Edep",100,0,.4);
 bst->Draw("Edep>>Edep_h"); 

 ///////////////////////////////////////////
 //   Draw our graphs                     //
 ///////////////////////////////////////////
 //Edep is in the svtT Tree
 graphPad1->Clear();

 Edep_h->Draw();

 // ------------------------------------------------

 // Main application eventloop. Calls system dependent eventloop
 // -> runs an interactive interface 
 theApp.Run();

 return(0); 

} /********** End of main *************/