SiPixelLorentzAngle/test/plotLorentzAngle.C

0001 {
0002     setTDRStyle();
0003     TH1F * h_lorentzAngle_layer1 = new TH1F("lorentzAnglelayer1","lorentzAngle", 8 , 0.5, 8.5);
0004     TH1F * h_lorentzAngle_layer2 = new TH1F("lorentzAnglelayer2","lorentzAngle", 8 , 0.5, 8.5);
0005     TH1F * h_lorentzAngle_layer3 = new TH1F("lorentzAnglelayer3","lorentzAngle", 8 , 0.5, 8.5);
0006     ifstream fin;
0007     fin.open("lorentzFit.txt");
0008     int module,  layer;
0009     double offset,  error,   slope,   error,   rel, pull,    chi2,    prob;
0010     fin.ignore(100,'\n');
0011 //  fin >> module >>  layer >> offset >> error >> slope >>   error >>   rel >> pull >>    chi2 >>    prob;
0012 //  for(int i_layer = 1; i_layer <= 3; i_layer++){
0013     for(int i_module = 1; i_module <= 8; i_module++){
0014         fin >> module >>  layer >> offset >> error >> slope >>   error >>   rel >> pull >>    chi2 >>    prob;
0015         cout << slope << endl;
0016         h_lorentzAngle_layer1->SetBinContent(i_module, slope);
0017         h_lorentzAngle_layer1->SetBinError(i_module, error);
0018     }
0019     for(int i_module = 1; i_module <= 8; i_module++){
0020         fin >> module >>  layer >> offset >> error >> slope >>   error >>   rel >> pull >>    chi2 >>    prob;
0021         cout << slope << endl;
0022         h_lorentzAngle_layer2->SetBinContent(i_module, slope);
0023         h_lorentzAngle_layer2->SetBinError(i_module, error);
0024     }
0025     for(int i_module = 1; i_module <= 8; i_module++){
0026         fin >> module >>  layer >> offset >> error >> slope >>   error >>   rel >> pull >>    chi2 >>    prob;
0027         cout << slope << endl;
0028         h_lorentzAngle_layer3->SetBinContent(i_module, slope);
0029         h_lorentzAngle_layer3->SetBinError(i_module, error);
0030     }
0031
0032     h_lorentzAngle_layer3->SetMarkerColor(7);
0033     h_lorentzAngle_layer3->SetLineColor(7);
0034     h_lorentzAngle_layer3->SetMarkerStyle(20);
0035     h_lorentzAngle_layer3->SetMarkerSize(0.9);
0036     h_lorentzAngle_layer3->SetMaximum(0.55);
0037     h_lorentzAngle_layer3->SetMinimum(0.35);
0038     h_lorentzAngle_layer3->GetXaxis()->SetTitle("module");
0039     h_lorentzAngle_layer3->GetYaxis()->SetTitle("tan(#Theta_{L})   ");
0040
0041     h_lorentzAngle_layer3->Draw("");
0042
0043     TLine* l = new TLine(0.5, 0.424, 8.5, 0.424);
0044     l->SetLineColor(2);
0045     l->SetLineWidth(2);
0046     l->Draw();
0047
0048     h_lorentzAngle_layer2->SetMarkerColor(4);
0049     h_lorentzAngle_layer2->SetLineColor(4);
0050     h_lorentzAngle_layer2->SetMarkerStyle(21);
0051     h_lorentzAngle_layer2->SetMarkerSize(0.7);
0052     h_lorentzAngle_layer2->Draw("same");
0053
0054     h_lorentzAngle_layer1->SetMarkerColor(6);
0055     h_lorentzAngle_layer1->SetLineColor(6);
0056     h_lorentzAngle_layer1->SetMarkerStyle(22);
0057     h_lorentzAngle_layer1->SetMarkerSize(1.1);
0058     h_lorentzAngle_layer1->Draw("same");
0059
0060     TLegend *legend1=new TLegend(0.7,0.75,0.94,0.94);
0061 //  legend1->SetTextSize(0.045);
0062     legend1->AddEntry(h_lorentzAngle_layer1,"layer 1","lp");
0063     legend1->AddEntry(h_lorentzAngle_layer2,"layer 2","lp");
0064     legend1->AddEntry(h_lorentzAngle_layer3,"layer 3","lp");
0065     legend1->Draw("");
0066
0067 }