Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 11:56:23

0001 #include <string>
0002 #include <sstream>
0003     
0004 #include "TFile.h"
0005 #include "TList.h"
0006 #include "TNtuple.h"
0007 #include "TTree.h"
0008 
0009 void compareGeomSerial(){
0010     
0011     
0012     //gStyle->SetOptStat(000000000);
0013     gStyle->SetOptStat("emr");
0014     gStyle->SetPadLeftMargin(0.15);
0015     gStyle->SetPadTopMargin(0.20);
0016     gStyle->SetTitleFontSize(0.08);
0017     
0018     
0019 
0020     TFile fin("<PATH>/comparisonV3_<N>.root");
0021     
0022     fin.cd();
0023     TTree* data = alignTree;
0024     
0025     // subdetectors ("sublevel"): PXB (1), PXF, TIB (3), TID, TOB, TEC (6)
0026     TCut levelCut = "((level == 1) && (sublevel >= 3))";
0027     
0028     TCanvas* c = new TCanvas("c", "c", 200, 10, 800, 800);
0029     c->SetFillColor(0);
0030     data->SetMarkerStyle(6);
0031     
0032     TH2D* hist2d = new TH2D("hist2D", "#Delta r vs. r; #Delta r (cm); r (cm)", 500, 0, 130, 1000, -1.0, 1.0);
0033     data->Project("hist2D", "r*dphi:r", levelCut);
0034     
0035     
0036     hist2d->Draw();
0037     TAxis *xaxis = hist2d->GetXaxis();
0038     TAxis *yaxis = hist2d->GetYaxis();
0039     xaxis->SetTitle("");
0040     yaxis->SetTitle("");
0041     xaxis->CenterTitle(true);
0042     yaxis->CenterTitle(true);
0043     xaxis->SetTitleSize(0.06);
0044     yaxis->SetTitleSize(0.06);
0045     yaxis->SetTitleOffset(-0.25);
0046     yaxis->SetLabelSize(.06);
0047     xaxis->SetLabelSize(.06);
0048     
0049     
0050     c->Print("<PATH>/comparisonV3_<N>.jpg");
0051 
0052 
0053     
0054 }