scripts/Beam/ReadBeamSpot.cpp

#include "ReadBeamSpot.h"
#include <TFile.h>
#include <TH1F.h>
#include <TCanvas.h>
#include <TROOT.h>
#include <iostream>
#include <vector>
#include <map>
#include <string>
#include <TStyle.h>
#include <TPaveText.h>
#include <TLegend.h>
#include <TGaxis.h>

bool myci_equal(char ch1, char ch2){
  return toupper((unsigned char)ch1) == toupper((unsigned char)ch2);
}

//------------------------------------------------------------
ReadBeamSpot::ReadBeamSpot(){
  string mthn = "[ReadBeamSpot::ReadBeamSpot()]\t";
  text_ = "";
  dataDir_ = "./";
}

//------------------------------------------------------------
ReadBeamSpot::~ReadBeamSpot(){
}

//------------------------------------------------------------
void ReadBeamSpot::openHtmlFile(string fileName){
  text_ = "";
  theFile_.open(fileName.c_str());
  if(!theFile_.is_open()){
    cout << "[ReadCsvFiles::openCsvFile()]\t" << "Can't open file: " << fileName << endl;
  }
}

//------------------------------------------------------------
void ReadBeamSpot::closeHtmlFile(void){
  if(theFile_.is_open()){
    theFile_.close();
  }  
  text_ = "";
  theFile_.clear();
}

//------------------------------------------------------------
void ReadBeamSpot::makePlots(){
}

void ReadBeamSpot::setStyle(){
  gStyle->SetFrameBorderMode(0);
  gStyle->SetCanvasBorderMode(0);
  gStyle->SetPadBorderMode(0);
  gStyle->SetFrameFillStyle(0);
  gStyle->SetLegendBorderSize(0);
  
  gStyle->SetPadColor(0);
  gStyle->SetCanvasColor(0);
  gStyle->SetTitleColor(1);
  gStyle->SetStatColor(0);
  
  // set the paper & margin sizes
  gStyle->SetPaperSize(20,26);
  gStyle->SetPadTopMargin(0.12);
  //gStyle->SetPadRightMargin(0.05);
  gStyle->SetPadBottomMargin(0.2);
  gStyle->SetPadLeftMargin(0.2);
  gStyle->SetPadTickX(1);
  gStyle->SetPadTickY(1);
  
  //
  gStyle->SetTextFont(132);
  gStyle->SetTextSize(0.09);
  gStyle->SetLabelFont(42,"xyz");
  gStyle->SetTitleFont(42,"xyz");
  //gStyle->SetLabelSize(0.032,"xyz");
  gStyle->SetLabelSize(0.035,"xyz");
  gStyle->SetTitleSize(0.045,"xyz");
  gStyle->SetTitleOffset(1.2,"xy");
  
  // use bold lines and markers
  gStyle->SetMarkerStyle(8);
  gStyle->SetHistLineWidth(2.3);
  gStyle->SetLineWidth(1.5);
  //gStyle->SetLineStyleString(2,"[12 12]"); // postscript dashes

  // do not display any of the standard histogram decorations
  gStyle->SetOptTitle(0); //0
  gStyle->SetOptStat(0); //0
  gStyle->SetOptFit(0); //0

  gStyle->SetPalette(1,(Int_t*)0);
  gStyle->SetErrorX(0);
  //  gStyle->SetEndErrorSize(0.5);
}

//------------------------------------------------------------
void ReadBeamSpot::findString(string what){
  setStyle();
  if(text_ == ""){
    char line[65536];
    while (theFile_.good()){
      theFile_.getline(line, 65536);
      text_ += string(line) + string(" ");
    }
  }
  int stringPos = 0;
  string runNumber;
  float x0,x0Error;
  float y0,y0Error;
  float z0,z0Error;
  float sigmaZ0,sigmaZ0Error;
  int valBegin = 0;
  int valSplit = 0;
  int valEnd = 0;

  int numberOfRuns = 0;
  while(stringPos >= 0){
    stringPos = ci_find(text_,what,stringPos);
    if(stringPos < 0) break;
    ++numberOfRuns;
    ++stringPos;
  }


  map<string,string> runToFillMap;
  runToFillMap["123815"] = "904";
  runToFillMap["123908"] = "906";
  runToFillMap["123970"] = "907";
  runToFillMap["123976"] = "907";
  runToFillMap["124023"] = "911";
  runToFillMap["124024"] = "911";
  runToFillMap["124025"] = "911";
  runToFillMap["124027"] = "912";
  runToFillMap["124030"] = "912";
  runToFillMap["124120"] = "916 (2.36 TeV)";
  runToFillMap["124230"] = "919";
  runToFillMap["124275"] = "923 (2.36 TeV)";

  TFile rootFile("Output.root","RECREATE");
  TH1F* hX0VsRun      = new TH1F("hX0VsRun",     "X0 vs Run number",     numberOfRuns,0,numberOfRuns);
  TH1F* hY0VsRun      = new TH1F("hY0VsRun",     "Y0 vs Run number",     numberOfRuns,0,numberOfRuns);
  TH1F* hZ0VsRun      = new TH1F("hZ0VsRun",     "Z0 vs Run number",     numberOfRuns,0,numberOfRuns);
  TH1F* hSigmaZ0VsRun = new TH1F("hSigmaZ0VsRun","SigmaZ0 vs Run number",numberOfRuns,0,numberOfRuns);
  
  hX0VsRun->GetXaxis()->SetTitle("Run");
  hX0VsRun->GetYaxis()->SetTitle("X position of beam spots (cm)");

  hY0VsRun->GetXaxis()->SetTitle("Run");
  hY0VsRun->GetYaxis()->SetTitle("Y position of beam spots (cm)");

  hZ0VsRun->GetXaxis()->SetTitle("Run");
  hZ0VsRun->GetYaxis()->SetTitle("Z position of beam spots (cm)");

  hSigmaZ0VsRun->GetXaxis()->SetTitle("Run");
  hSigmaZ0VsRun->GetYaxis()->SetTitle("#sigma_{z_{0}} of beam spots (cm)");

  TH1F* hX0VsFill      = new TH1F("hX0VsFill",     "X0 vs Fill number",     numberOfRuns,0,numberOfRuns);
  TH1F* hY0VsFill      = new TH1F("hY0VsFill",     "Y0 vs Fill number",     numberOfRuns,0,numberOfRuns);
  TH1F* hZ0VsFill      = new TH1F("hZ0VsFill",     "Z0 vs Fill number",     numberOfRuns,0,numberOfRuns);
  TH1F* hSigmaZ0VsFill = new TH1F("hSigmaZ0VsFill","SigmaZ0 vs Fill number",numberOfRuns,0,numberOfRuns);
  
  hX0VsFill->GetXaxis()->SetTitle("LHC Fill");
  hX0VsFill->GetYaxis()->SetTitle("X position of beam spots (cm)");

  hY0VsFill->GetXaxis()->SetTitle("LHC Fill");
  hY0VsFill->GetYaxis()->SetTitle("Y position of beam spots (cm)");

  hZ0VsFill->GetXaxis()->SetTitle("LHC Fill");
  hZ0VsFill->GetYaxis()->SetTitle("Z position of beam spots (cm)");

  hSigmaZ0VsFill->GetXaxis()->SetTitle("LHC Fill");
  hSigmaZ0VsFill->GetYaxis()->SetTitle("#sigma_{z_{0}} of beam spots (cm)");


  stringPos=0;
  while(stringPos >= 0){
    //          stringPos = text_.find(what,stringPos);
    //          stringPos = ci_find(text_.substr(stringPos,text_.length()-stringPos),what);
    stringPos = ci_find(text_,what,stringPos);
    //    stringPos = ci_find(text_,what,stringPos);
    if(stringPos < 0) break;
    runNumber = text_.substr((text_.substr(0,stringPos)).rfind("<a name=\"Run_")+13,6);

    valBegin = ci_find(text_,"X0",stringPos)+13;
    valSplit = ci_find(text_,"+/-",stringPos);
    valEnd   = ci_find(text_,"&#91;cm]",stringPos);
    x0      = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
    x0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
    stringPos = valEnd+8;
    
    valBegin = ci_find(text_,"Y0",stringPos)+13;
    valSplit = ci_find(text_,"+/-",stringPos);
    valEnd   = ci_find(text_,"&#91;cm]",stringPos);
    y0      = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
    y0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
    stringPos = valEnd+8;

    valBegin = ci_find(text_,"Z0",stringPos)+13;
    valSplit = ci_find(text_,"+/-",stringPos);
    valEnd   = ci_find(text_,"&#91;cm]",stringPos);
    z0      = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
    z0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
    stringPos = valEnd+8;

    valBegin = ci_find(text_,"Sigma Z0",stringPos)+19;
    valSplit = ci_find(text_,"+/-",stringPos);
    valEnd   = ci_find(text_,"&#91;cm]",stringPos);
    sigmaZ0      = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
    sigmaZ0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
    stringPos = valEnd+8;

    cout << __PRETTY_FUNCTION__ 
         << "Pos: " << stringPos 
         << " run: " << runNumber << endl
      //         << "-X0: "   << x0
      //         << "-X0error: " << x0Error << "-" << endl
      //         << "-Y0: "   << y0
      //         << "-Y0error: " << y0Error  << "-" << endl
      //         << "-Z0: "   << z0
      //         << "-Z0error: " << z0Error  << "-" << endl
      //         << "-sigmaZ0: "   << sigmaZ0
      //         << "-sigmaZ0error: " << sigmaZ0Error
         << "-" << endl;
    TString runLabel;
    if (runNumber == "124120" || runNumber == "124275")
      runLabel = TString(runNumber.c_str())+" (2.36 TeV)";
    else
      runLabel = runNumber.c_str();

    hX0VsRun->SetBinContent(numberOfRuns,x0);
    hX0VsRun->SetBinError(numberOfRuns,x0Error);
    hX0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);

    hY0VsRun->SetBinContent(numberOfRuns,y0);
    hY0VsRun->SetBinError(numberOfRuns,y0Error);
    hY0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);

    hZ0VsRun->SetBinContent(numberOfRuns,z0);
    hZ0VsRun->SetBinError(numberOfRuns,z0Error);
    hZ0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);

    hSigmaZ0VsRun->SetBinContent(numberOfRuns,sigmaZ0);
    hSigmaZ0VsRun->SetBinError(numberOfRuns,sigmaZ0Error);
    hSigmaZ0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);
    
    hX0VsFill->SetBinContent(numberOfRuns,x0);
    hX0VsFill->SetBinError(numberOfRuns,x0Error);
    hX0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());

    hY0VsFill->SetBinContent(numberOfRuns,y0);
    hY0VsFill->SetBinError(numberOfRuns,y0Error);
    hY0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());

    hZ0VsFill->SetBinContent(numberOfRuns,z0);
    hZ0VsFill->SetBinError(numberOfRuns,z0Error);
    hZ0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());

    hSigmaZ0VsFill->SetBinContent(numberOfRuns,sigmaZ0);
    hSigmaZ0VsFill->SetBinError(numberOfRuns,sigmaZ0Error);
    hSigmaZ0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());
    ++stringPos;
    --numberOfRuns;
  }
  rootFile.Write();

  TPaveText *ppt= new TPaveText(.36,.72,.76,.78,"NDC");
  ppt->AddText("CMS Preliminary 2009");
  ppt->SetTextAlign(12);
  ppt->SetTextFont(62);
  ppt->SetTextSize(0.045);
  ppt->SetBorderSize(0);
  ppt->SetFillStyle(0);
  
  TCanvas cTmp("cTmp","Drawing canvas",800,600);
  hX0VsRun->GetYaxis()->SetRangeUser(0.085,0.135);
  hX0VsRun->GetXaxis()->LabelsOption("u");
  hX0VsRun->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hX0VsRun->GetTitle())+".png");
  hZ0VsRun->GetXaxis()->LabelsOption("u");
  hZ0VsRun->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hZ0VsRun->GetTitle())+".png");
  hSigmaZ0VsRun->GetXaxis()->LabelsOption("u");
  hSigmaZ0VsRun->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hSigmaZ0VsRun->GetTitle())+".png");
  hX0VsFill->GetYaxis()->SetRangeUser(0.085,0.135);
  hX0VsFill->GetXaxis()->LabelsOption("u");
  hX0VsFill->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hX0VsFill->GetTitle())+".png");
  hZ0VsFill->GetXaxis()->LabelsOption("u");
  hZ0VsFill->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hZ0VsFill->GetTitle())+".png");
  hSigmaZ0VsFill->GetXaxis()->LabelsOption("u");
  hSigmaZ0VsFill->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hSigmaZ0VsFill->GetTitle())+".png");

  ppt->SetX1NDC(0.45);
  ppt->SetY1NDC(0.72);
  ppt->SetX2NDC(0.85);
  ppt->SetY2NDC(0.78);
  hY0VsRun->GetYaxis()->SetRangeUser(0.01,0.06);
  hY0VsRun->GetXaxis()->LabelsOption("u");
  hY0VsRun->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hY0VsRun->GetTitle())+".png");
  hY0VsFill->GetYaxis()->SetRangeUser(0.01,0.06);
  hY0VsFill->GetXaxis()->LabelsOption("u");
  hY0VsFill->Draw("e1");
  ppt->Draw();
  cTmp.Print(TString(hY0VsFill->GetTitle())+".png");

  // Superimpose Transverse
  Double_t ymin,ymax;
  ymin=0.085;ymax=0.15;

  hX0VsRun->GetYaxis()->SetTitle("Transverse position of beam spots (cm)");
  hX0VsRun->GetXaxis()->LabelsOption("u");
  hX0VsRun->GetXaxis()->CenterLabels();
  hX0VsRun->Draw("e1");
  hX0VsRun->GetYaxis()->SetRangeUser(ymin,ymax);
  gPad->Update();
//   hY0VsRun->SetMarkerStyle(21);
//   hY0VsRun->SetMarkerColor(2);
//   hY0VsRun->SetLineColor(2);
//   hY0VsRun->Draw("samee1][");
  int nbins = hY0VsRun->GetEntries();
  TH1F*hY0tmp=new TH1F("y0_tmp","y0_tmp",nbins,0,nbins);
  double level=0.08;
  for (int i=0;i<nbins;i++){
    hY0tmp->SetBinContent(i+1,(level+hY0VsRun->GetBinContent(i+1)));
    hY0tmp->SetBinError(i+1,hY0VsRun->GetBinError(i+1));
  }
  hY0tmp->SetMarkerStyle(21);
  hY0tmp->SetMarkerColor(2);
  hY0tmp->SetLineColor(2);
  hY0tmp->Draw("samee1][");
//   cout << gPad->GetUymin() << endl;
  TGaxis*axis=new TGaxis(gPad->GetUxmax(),gPad->GetUymin(),gPad->GetUxmax(),gPad->GetUymax(),ymin-level,ymax-level,510,"+L");
  axis->SetLineColor(2);
  axis->SetLabelColor(2);
  axis->SetLabelFont(42);
  axis->SetLabelSize(0.035);
  axis->Draw();
  
  ppt->SetX1NDC(0.36);
  ppt->SetY1NDC(0.78);
  ppt->SetX2NDC(0.76);
  ppt->SetY2NDC(0.86);
  ppt->Draw();

  TLegend*l0=new TLegend(0.55,0.6,0.63,0.7);
  l0->SetFillStyle(0);
  l0->AddEntry(hX0VsRun,"  X_{0}","p");
//   l0->AddEntry(hY0VsRun,"  Y_{0}","p");
  l0->AddEntry(hY0tmp,"  Y_{0}","p");
  l0->Draw();
//   cTmp.Print(TString(hY0VsRun->GetTitle())+".png");
  cTmp.Print("Transverse Beam Potition vs Run number.png");

  hX0VsFill->GetYaxis()->SetTitle("Transverse position of beam spots (cm)");
  hX0VsFill->GetXaxis()->LabelsOption("u");
  hX0VsFill->Draw("e1");
  hX0VsFill->GetYaxis()->SetRangeUser(ymin,ymax);
//   hY0VsFill->SetMarkerStyle(21);
//   hY0VsFill->SetMarkerColor(2);
//   hY0VsFill->SetLineColor(2);
//   hY0VsFill->Draw("samee1][");
  hY0tmp->Draw("samee1][");
  axis->Draw();
  ppt->Draw();
//   TLegend*l1=new TLegend(0.6,0.5,0.68,0.6);
//   l1->SetFillStyle(0);
//   l1->AddEntry(hX0VsFill,"  X_{0}","p");
//   l1->AddEntry(hY0VsFill,"  Y_{0}","p");
  l0->Draw();
  cTmp.Print("Transverse Beam Potition vs Fill Number.png");

  rootFile.Close();
}

ostream& operator<<(ostream& out, const ReadBeamSpot& reader){
  return (out << reader.textFound_.str().c_str());
}

size_t ReadBeamSpot::ci_find(string& str1, string str2, int stringPos){
  string::iterator pos = search(str1.begin()+stringPos, str1.end(), str2.begin(), str2.end(), myci_equal);
  if (pos == str1.end())
    return string::npos;
  else
    return pos - str1.begin();
}
Revision:	1.1
Committed:	Fri Feb 12 18:03:59 2010 UTC (15 years, 2 months ago) by jengbou
Branch:	MAIN
CVS Tags:	HEAD
Log Message:	Add scripts for plotting beam positions vs run number (modified from Lorenzo's code)
#	Content
1	#include "ReadBeamSpot.h"
2	#include <TFile.h>
3	#include <TH1F.h>
4	#include <TCanvas.h>
5	#include <TROOT.h>
6	#include <iostream>
7	#include <vector>
8	#include <map>
9	#include <string>
10	#include <TStyle.h>
11	#include <TPaveText.h>
12	#include <TLegend.h>
13	#include <TGaxis.h>
14
15	bool myci_equal(char ch1, char ch2){
16	return toupper((unsigned char)ch1) == toupper((unsigned char)ch2);
17	}
18
19	//------------------------------------------------------------
20	ReadBeamSpot::ReadBeamSpot(){
21	string mthn = "[ReadBeamSpot::ReadBeamSpot()]\t";
22	text_ = "";
23	dataDir_ = "./";
24	}
25
26	//------------------------------------------------------------
27	ReadBeamSpot::~ReadBeamSpot(){
28	}
29
30	//------------------------------------------------------------
31	void ReadBeamSpot::openHtmlFile(string fileName){
32	text_ = "";
33	theFile_.open(fileName.c_str());
34	if(!theFile_.is_open()){
35	cout << "[ReadCsvFiles::openCsvFile()]\t" << "Can't open file: " << fileName << endl;
36	}
37	}
38
39	//------------------------------------------------------------
40	void ReadBeamSpot::closeHtmlFile(void){
41	if(theFile_.is_open()){
42	theFile_.close();
43	}
44	text_ = "";
45	theFile_.clear();
46	}
47
48	//------------------------------------------------------------
49	void ReadBeamSpot::makePlots(){
50	}
51
52	void ReadBeamSpot::setStyle(){
53	gStyle->SetFrameBorderMode(0);
54	gStyle->SetCanvasBorderMode(0);
55	gStyle->SetPadBorderMode(0);
56	gStyle->SetFrameFillStyle(0);
57	gStyle->SetLegendBorderSize(0);
58
59	gStyle->SetPadColor(0);
60	gStyle->SetCanvasColor(0);
61	gStyle->SetTitleColor(1);
62	gStyle->SetStatColor(0);
63
64	// set the paper & margin sizes
65	gStyle->SetPaperSize(20,26);
66	gStyle->SetPadTopMargin(0.12);
67	//gStyle->SetPadRightMargin(0.05);
68	gStyle->SetPadBottomMargin(0.2);
69	gStyle->SetPadLeftMargin(0.2);
70	gStyle->SetPadTickX(1);
71	gStyle->SetPadTickY(1);
72
73	//
74	gStyle->SetTextFont(132);
75	gStyle->SetTextSize(0.09);
76	gStyle->SetLabelFont(42,"xyz");
77	gStyle->SetTitleFont(42,"xyz");
78	//gStyle->SetLabelSize(0.032,"xyz");
79	gStyle->SetLabelSize(0.035,"xyz");
80	gStyle->SetTitleSize(0.045,"xyz");
81	gStyle->SetTitleOffset(1.2,"xy");
82
83	// use bold lines and markers
84	gStyle->SetMarkerStyle(8);
85	gStyle->SetHistLineWidth(2.3);
86	gStyle->SetLineWidth(1.5);
87	//gStyle->SetLineStyleString(2,"[12 12]"); // postscript dashes
88
89	// do not display any of the standard histogram decorations
90	gStyle->SetOptTitle(0); //0
91	gStyle->SetOptStat(0); //0
92	gStyle->SetOptFit(0); //0
93
94	gStyle->SetPalette(1,(Int_t*)0);
95	gStyle->SetErrorX(0);
96	// gStyle->SetEndErrorSize(0.5);
97	}
98
99	//------------------------------------------------------------
100	void ReadBeamSpot::findString(string what){
101	setStyle();
102	if(text_ == ""){
103	char line[65536];
104	while (theFile_.good()){
105	theFile_.getline(line, 65536);
106	text_ += string(line) + string(" ");
107	}
108	}
109	int stringPos = 0;
110	string runNumber;
111	float x0,x0Error;
112	float y0,y0Error;
113	float z0,z0Error;
114	float sigmaZ0,sigmaZ0Error;
115	int valBegin = 0;
116	int valSplit = 0;
117	int valEnd = 0;
118
119	int numberOfRuns = 0;
120	while(stringPos >= 0){
121	stringPos = ci_find(text_,what,stringPos);
122	if(stringPos < 0) break;
123	++numberOfRuns;
124	++stringPos;
125	}
126
127
128	map<string,string> runToFillMap;
129	runToFillMap["123815"] = "904";
130	runToFillMap["123908"] = "906";
131	runToFillMap["123970"] = "907";
132	runToFillMap["123976"] = "907";
133	runToFillMap["124023"] = "911";
134	runToFillMap["124024"] = "911";
135	runToFillMap["124025"] = "911";
136	runToFillMap["124027"] = "912";
137	runToFillMap["124030"] = "912";
138	runToFillMap["124120"] = "916 (2.36 TeV)";
139	runToFillMap["124230"] = "919";
140	runToFillMap["124275"] = "923 (2.36 TeV)";
141
142	TFile rootFile("Output.root","RECREATE");
143	TH1F* hX0VsRun = new TH1F("hX0VsRun", "X0 vs Run number", numberOfRuns,0,numberOfRuns);
144	TH1F* hY0VsRun = new TH1F("hY0VsRun", "Y0 vs Run number", numberOfRuns,0,numberOfRuns);
145	TH1F* hZ0VsRun = new TH1F("hZ0VsRun", "Z0 vs Run number", numberOfRuns,0,numberOfRuns);
146	TH1F* hSigmaZ0VsRun = new TH1F("hSigmaZ0VsRun","SigmaZ0 vs Run number",numberOfRuns,0,numberOfRuns);
147
148	hX0VsRun->GetXaxis()->SetTitle("Run");
149	hX0VsRun->GetYaxis()->SetTitle("X position of beam spots (cm)");
150
151	hY0VsRun->GetXaxis()->SetTitle("Run");
152	hY0VsRun->GetYaxis()->SetTitle("Y position of beam spots (cm)");
153
154	hZ0VsRun->GetXaxis()->SetTitle("Run");
155	hZ0VsRun->GetYaxis()->SetTitle("Z position of beam spots (cm)");
156
157	hSigmaZ0VsRun->GetXaxis()->SetTitle("Run");
158	hSigmaZ0VsRun->GetYaxis()->SetTitle("#sigma_{z_{0}} of beam spots (cm)");
159
160	TH1F* hX0VsFill = new TH1F("hX0VsFill", "X0 vs Fill number", numberOfRuns,0,numberOfRuns);
161	TH1F* hY0VsFill = new TH1F("hY0VsFill", "Y0 vs Fill number", numberOfRuns,0,numberOfRuns);
162	TH1F* hZ0VsFill = new TH1F("hZ0VsFill", "Z0 vs Fill number", numberOfRuns,0,numberOfRuns);
163	TH1F* hSigmaZ0VsFill = new TH1F("hSigmaZ0VsFill","SigmaZ0 vs Fill number",numberOfRuns,0,numberOfRuns);
164
165	hX0VsFill->GetXaxis()->SetTitle("LHC Fill");
166	hX0VsFill->GetYaxis()->SetTitle("X position of beam spots (cm)");
167
168	hY0VsFill->GetXaxis()->SetTitle("LHC Fill");
169	hY0VsFill->GetYaxis()->SetTitle("Y position of beam spots (cm)");
170
171	hZ0VsFill->GetXaxis()->SetTitle("LHC Fill");
172	hZ0VsFill->GetYaxis()->SetTitle("Z position of beam spots (cm)");
173
174	hSigmaZ0VsFill->GetXaxis()->SetTitle("LHC Fill");
175	hSigmaZ0VsFill->GetYaxis()->SetTitle("#sigma_{z_{0}} of beam spots (cm)");
176
177
178	stringPos=0;
179	while(stringPos >= 0){
180	// stringPos = text_.find(what,stringPos);
181	// stringPos = ci_find(text_.substr(stringPos,text_.length()-stringPos),what);
182	stringPos = ci_find(text_,what,stringPos);
183	// stringPos = ci_find(text_,what,stringPos);
184	if(stringPos < 0) break;
185	runNumber = text_.substr((text_.substr(0,stringPos)).rfind("<a name=\"Run_")+13,6);
186
187	valBegin = ci_find(text_,"X0",stringPos)+13;
188	valSplit = ci_find(text_,"+/-",stringPos);
189	valEnd = ci_find(text_,"[cm]",stringPos);
190	x0 = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
191	x0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
192	stringPos = valEnd+8;
193
194	valBegin = ci_find(text_,"Y0",stringPos)+13;
195	valSplit = ci_find(text_,"+/-",stringPos);
196	valEnd = ci_find(text_,"[cm]",stringPos);
197	y0 = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
198	y0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
199	stringPos = valEnd+8;
200
201	valBegin = ci_find(text_,"Z0",stringPos)+13;
202	valSplit = ci_find(text_,"+/-",stringPos);
203	valEnd = ci_find(text_,"[cm]",stringPos);
204	z0 = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
205	z0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
206	stringPos = valEnd+8;
207
208	valBegin = ci_find(text_,"Sigma Z0",stringPos)+19;
209	valSplit = ci_find(text_,"+/-",stringPos);
210	valEnd = ci_find(text_,"[cm]",stringPos);
211	sigmaZ0 = atof(text_.substr(valBegin, valSplit-valBegin-1).c_str());
212	sigmaZ0Error = atof(text_.substr(valSplit+4, valEnd-valSplit-5).c_str());
213	stringPos = valEnd+8;
214
215	cout << __PRETTY_FUNCTION__
216	<< "Pos: " << stringPos
217	<< " run: " << runNumber << endl
218	// << "-X0: " << x0
219	// << "-X0error: " << x0Error << "-" << endl
220	// << "-Y0: " << y0
221	// << "-Y0error: " << y0Error << "-" << endl
222	// << "-Z0: " << z0
223	// << "-Z0error: " << z0Error << "-" << endl
224	// << "-sigmaZ0: " << sigmaZ0
225	// << "-sigmaZ0error: " << sigmaZ0Error
226	<< "-" << endl;
227	TString runLabel;
228	if (runNumber == "124120" \|\| runNumber == "124275")
229	runLabel = TString(runNumber.c_str())+" (2.36 TeV)";
230	else
231	runLabel = runNumber.c_str();
232
233	hX0VsRun->SetBinContent(numberOfRuns,x0);
234	hX0VsRun->SetBinError(numberOfRuns,x0Error);
235	hX0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);
236
237	hY0VsRun->SetBinContent(numberOfRuns,y0);
238	hY0VsRun->SetBinError(numberOfRuns,y0Error);
239	hY0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);
240
241	hZ0VsRun->SetBinContent(numberOfRuns,z0);
242	hZ0VsRun->SetBinError(numberOfRuns,z0Error);
243	hZ0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);
244
245	hSigmaZ0VsRun->SetBinContent(numberOfRuns,sigmaZ0);
246	hSigmaZ0VsRun->SetBinError(numberOfRuns,sigmaZ0Error);
247	hSigmaZ0VsRun->GetXaxis()->SetBinLabel(numberOfRuns,runLabel);
248
249	hX0VsFill->SetBinContent(numberOfRuns,x0);
250	hX0VsFill->SetBinError(numberOfRuns,x0Error);
251	hX0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());
252
253	hY0VsFill->SetBinContent(numberOfRuns,y0);
254	hY0VsFill->SetBinError(numberOfRuns,y0Error);
255	hY0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());
256
257	hZ0VsFill->SetBinContent(numberOfRuns,z0);
258	hZ0VsFill->SetBinError(numberOfRuns,z0Error);
259	hZ0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());
260
261	hSigmaZ0VsFill->SetBinContent(numberOfRuns,sigmaZ0);
262	hSigmaZ0VsFill->SetBinError(numberOfRuns,sigmaZ0Error);
263	hSigmaZ0VsFill->GetXaxis()->SetBinLabel(numberOfRuns,runToFillMap[runNumber].c_str());
264	++stringPos;
265	--numberOfRuns;
266	}
267	rootFile.Write();
268
269	TPaveText *ppt= new TPaveText(.36,.72,.76,.78,"NDC");
270	ppt->AddText("CMS Preliminary 2009");
271	ppt->SetTextAlign(12);
272	ppt->SetTextFont(62);
273	ppt->SetTextSize(0.045);
274	ppt->SetBorderSize(0);
275	ppt->SetFillStyle(0);
276
277	TCanvas cTmp("cTmp","Drawing canvas",800,600);
278	hX0VsRun->GetYaxis()->SetRangeUser(0.085,0.135);
279	hX0VsRun->GetXaxis()->LabelsOption("u");
280	hX0VsRun->Draw("e1");
281	ppt->Draw();
282	cTmp.Print(TString(hX0VsRun->GetTitle())+".png");
283	hZ0VsRun->GetXaxis()->LabelsOption("u");
284	hZ0VsRun->Draw("e1");
285	ppt->Draw();
286	cTmp.Print(TString(hZ0VsRun->GetTitle())+".png");
287	hSigmaZ0VsRun->GetXaxis()->LabelsOption("u");
288	hSigmaZ0VsRun->Draw("e1");
289	ppt->Draw();
290	cTmp.Print(TString(hSigmaZ0VsRun->GetTitle())+".png");
291	hX0VsFill->GetYaxis()->SetRangeUser(0.085,0.135);
292	hX0VsFill->GetXaxis()->LabelsOption("u");
293	hX0VsFill->Draw("e1");
294	ppt->Draw();
295	cTmp.Print(TString(hX0VsFill->GetTitle())+".png");
296	hZ0VsFill->GetXaxis()->LabelsOption("u");
297	hZ0VsFill->Draw("e1");
298	ppt->Draw();
299	cTmp.Print(TString(hZ0VsFill->GetTitle())+".png");
300	hSigmaZ0VsFill->GetXaxis()->LabelsOption("u");
301	hSigmaZ0VsFill->Draw("e1");
302	ppt->Draw();
303	cTmp.Print(TString(hSigmaZ0VsFill->GetTitle())+".png");
304
305	ppt->SetX1NDC(0.45);
306	ppt->SetY1NDC(0.72);
307	ppt->SetX2NDC(0.85);
308	ppt->SetY2NDC(0.78);
309	hY0VsRun->GetYaxis()->SetRangeUser(0.01,0.06);
310	hY0VsRun->GetXaxis()->LabelsOption("u");
311	hY0VsRun->Draw("e1");
312	ppt->Draw();
313	cTmp.Print(TString(hY0VsRun->GetTitle())+".png");
314	hY0VsFill->GetYaxis()->SetRangeUser(0.01,0.06);
315	hY0VsFill->GetXaxis()->LabelsOption("u");
316	hY0VsFill->Draw("e1");
317	ppt->Draw();
318	cTmp.Print(TString(hY0VsFill->GetTitle())+".png");
319
320	// Superimpose Transverse
321	Double_t ymin,ymax;
322	ymin=0.085;ymax=0.15;
323
324	hX0VsRun->GetYaxis()->SetTitle("Transverse position of beam spots (cm)");
325	hX0VsRun->GetXaxis()->LabelsOption("u");
326	hX0VsRun->GetXaxis()->CenterLabels();
327	hX0VsRun->Draw("e1");
328	hX0VsRun->GetYaxis()->SetRangeUser(ymin,ymax);
329	gPad->Update();
330	// hY0VsRun->SetMarkerStyle(21);
331	// hY0VsRun->SetMarkerColor(2);
332	// hY0VsRun->SetLineColor(2);
333	// hY0VsRun->Draw("samee1][");
334	int nbins = hY0VsRun->GetEntries();
335	TH1F*hY0tmp=new TH1F("y0_tmp","y0_tmp",nbins,0,nbins);
336	double level=0.08;
337	for (int i=0;i<nbins;i++){
338	hY0tmp->SetBinContent(i+1,(level+hY0VsRun->GetBinContent(i+1)));
339	hY0tmp->SetBinError(i+1,hY0VsRun->GetBinError(i+1));
340	}
341	hY0tmp->SetMarkerStyle(21);
342	hY0tmp->SetMarkerColor(2);
343	hY0tmp->SetLineColor(2);
344	hY0tmp->Draw("samee1][");
345	// cout << gPad->GetUymin() << endl;
346	TGaxis*axis=new TGaxis(gPad->GetUxmax(),gPad->GetUymin(),gPad->GetUxmax(),gPad->GetUymax(),ymin-level,ymax-level,510,"+L");
347	axis->SetLineColor(2);
348	axis->SetLabelColor(2);
349	axis->SetLabelFont(42);
350	axis->SetLabelSize(0.035);
351	axis->Draw();
352
353	ppt->SetX1NDC(0.36);
354	ppt->SetY1NDC(0.78);
355	ppt->SetX2NDC(0.76);
356	ppt->SetY2NDC(0.86);
357	ppt->Draw();
358
359	TLegend*l0=new TLegend(0.55,0.6,0.63,0.7);
360	l0->SetFillStyle(0);
361	l0->AddEntry(hX0VsRun," X_{0}","p");
362	// l0->AddEntry(hY0VsRun," Y_{0}","p");
363	l0->AddEntry(hY0tmp," Y_{0}","p");
364	l0->Draw();
365	// cTmp.Print(TString(hY0VsRun->GetTitle())+".png");
366	cTmp.Print("Transverse Beam Potition vs Run number.png");
367
368	hX0VsFill->GetYaxis()->SetTitle("Transverse position of beam spots (cm)");
369	hX0VsFill->GetXaxis()->LabelsOption("u");
370	hX0VsFill->Draw("e1");
371	hX0VsFill->GetYaxis()->SetRangeUser(ymin,ymax);
372	// hY0VsFill->SetMarkerStyle(21);
373	// hY0VsFill->SetMarkerColor(2);
374	// hY0VsFill->SetLineColor(2);
375	// hY0VsFill->Draw("samee1][");
376	hY0tmp->Draw("samee1][");
377	axis->Draw();
378	ppt->Draw();
379	// TLegend*l1=new TLegend(0.6,0.5,0.68,0.6);
380	// l1->SetFillStyle(0);
381	// l1->AddEntry(hX0VsFill," X_{0}","p");
382	// l1->AddEntry(hY0VsFill," Y_{0}","p");
383	l0->Draw();
384	cTmp.Print("Transverse Beam Potition vs Fill Number.png");
385
386	rootFile.Close();
387	}
388
389	ostream& operator<<(ostream& out, const ReadBeamSpot& reader){
390	return (out << reader.textFound_.str().c_str());
391	}
392
393	size_t ReadBeamSpot::ci_find(string& str1, string str2, int stringPos){
394	string::iterator pos = search(str1.begin()+stringPos, str1.end(), str2.begin(), str2.end(), myci_equal);
395	if (pos == str1.end())
396	return string::npos;
397	else
398	return pos - str1.begin();
399	}