CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Functions
multibsvspvplots.h File Reference

Go to the source code of this file.

Functions

void multibsvspvplots (const char *module, const char *label, const char *postfix, const char *)
 

Function Documentation

void multibsvspvplots ( const char *  module,
const char *  label,
const char *  postfix,
const char *   
)

Definition at line 22 of file multibsvspvplots.cc.

References newFWLiteAna::bin, gather_cfg::cout, mergeVDriftHistosByStation::file, CommonAnalyzer::getObject(), CommonAnalyzer::getRunList(), i, compare_using_db::ifile, gather_cfg::runs, CommonAnalyzer::setPath(), python.multivaluedict::sort(), and AlCaHLTBitMon_QueryRunRegistry::string.

22  {
23 
24  TFile* file[5];
25 
26  file[0] = new TFile("rootfiles/Tracking_PFG_Run2012A_prompt_minbias_v1_190456-194076_muon_relumi_v55_fittedV0.root");
27  file[1] = new TFile("rootfiles/Tracking_PFG_Run2012B_prompt_minbias_v1_190456-196531_muon_relumi_v55_fittedV0.root");
28  file[2] = new TFile("rootfiles/Tracking_PFG_Run2012C_prompt_minbias_v1_190456-199011_muon_relumi_v55_fittedV0.root");
29  file[3] = new TFile("rootfiles/Tracking_PFG_Run2012C_prompt_minbias_v2_190456-203002_muon_relumi_v57_fittedV0.root");
30  file[4] = new TFile("rootfiles/Tracking_PFG_Run2012D_prompt_minbias_v1_190456-208686_muon_relumi_v57_fittedV0.root");
31 
32  char modfull[300];
33  sprintf(modfull,"%s%s",module,postfix);
34 
35  char labfull[300];
36  sprintf(labfull,"%s%s",label,postfix);
37 
38  // Summary histograms
39  TH1D* deltaxsum = new TH1D("deltaxsum","(PV-BS) Fitted X position vs run",10,0.,10.);
40  deltaxsum->SetCanExtend(TH1::kAllAxes);
41  TH1D* deltaysum = new TH1D("deltaysum","(PV-BS) Fitted Y position vs run",10,0.,10.);
42  deltaysum->SetCanExtend(TH1::kAllAxes);
43 
44  TH1D* deltaxmeansum = new TH1D("deltaxmeansum","(PV-BS) Mean X position vs run",10,0.,10.);
45  deltaxmeansum->SetCanExtend(TH1::kAllAxes);
46  TH1D* deltaymeansum = new TH1D("deltaymeansum","(PV-BS) Mean Y position vs run",10,0.,10.);
47  deltaymeansum->SetCanExtend(TH1::kAllAxes);
48 
49  TH1D* deltazmeansum = new TH1D("deltazmeansum","(PV-BS) Mean Z position vs run",10,0.,10.);
50  deltazmeansum->SetCanExtend(TH1::kAllAxes);
51 
52  TF1* fdoubleg = new TF1("doubleg","[1]*exp(-0.5*((x-[0])/[2])**2)+[3]*exp(-0.5*((x-[0])/[4])**2)+[5]*exp(sqrt((x-[0])**2)*[6])",-.2,.2);
53  fdoubleg->SetLineColor(kRed);
54  fdoubleg->SetLineWidth(1);
55 
56  for(unsigned int ifile=0;ifile<5;++ifile) {
57 
58  CommonAnalyzer castat(file[ifile],"",modfull);
59 
60  std::vector<unsigned int> runs = castat.getRunList();
61  std::sort(runs.begin(),runs.end());
62 
63  {
64 
65  std::cout << "Found " << runs.size() << " runs" << std::endl;
66 
67  for(unsigned int i=0;i<runs.size();++i) {
68 
69  char runlabel[100];
70  sprintf(runlabel,"%d",runs[i]);
71  char runpath[100];
72  sprintf(runpath,"run_%d",runs[i]);
73  castat.setPath(runpath);
74  std::cout << runpath << std::endl;
75 
76  TH1F* deltax = (TH1F*)castat.getObject("deltaxrun");
77  if (deltax && deltax->GetEntries()>0) {
78 
79  fdoubleg->SetParameter(0,deltax->GetMean());
80  fdoubleg->SetParameter(2,deltax->GetRMS());
81  fdoubleg->SetParameter(4,deltax->GetRMS());
82  fdoubleg->SetParameter(1,deltax->GetMaximum());
83  fdoubleg->SetParameter(3,0.1*deltax->GetMaximum());
84  fdoubleg->SetParameter(5,0.1*deltax->GetMaximum());
85  /* const int result = */ deltax->Fit(fdoubleg,"q0b","",-.05,.05);
86 
87  int bin = deltaxsum->Fill(runlabel,fdoubleg->GetParameter(0));
88  deltaxsum->SetBinError(bin,fdoubleg->GetParError(0));
89 
90  bin = deltaxmeansum->Fill(runlabel,deltax->GetMean());
91  deltaxmeansum->SetBinError(bin,deltax->GetMeanError());
92 
93  }
94  TH1F* deltay = (TH1F*)castat.getObject("deltayrun");
95  if (deltay && deltay->GetEntries()>0) {
96 
97  fdoubleg->SetParameter(0,deltay->GetMean());
98  fdoubleg->SetParameter(2,deltay->GetRMS());
99  fdoubleg->SetParameter(4,deltay->GetRMS());
100  fdoubleg->SetParameter(1,deltay->GetMaximum());
101  fdoubleg->SetParameter(3,0.1*deltay->GetMaximum());
102  fdoubleg->SetParameter(5,0.1*deltay->GetMaximum());
103  /* const int result = */ deltay->Fit(fdoubleg,"q0b","",-.05,.05);
104  int bin = deltaysum->Fill(runlabel,fdoubleg->GetParameter(0));
105  deltaysum->SetBinError(bin,fdoubleg->GetParError(0));
106 
107  bin = deltaymeansum->Fill(runlabel,deltay->GetMean());
108  deltaymeansum->SetBinError(bin,deltay->GetMeanError());
109 
110  }
111  TH1F* deltaz = (TH1F*)castat.getObject("deltazrun");
112  if (deltaz && deltaz->GetEntries()>0) {
113 
114  int bin = deltazmeansum->Fill(runlabel,deltaz->GetMean());
115  deltazmeansum->SetBinError(bin,deltaz->GetMeanError());
116 
117  }
118 
119 
120 
121  }
122 
123  }
124  }
125 
126  std::string plotfilename;
127 
128  plotfilename = "/afs/cern.ch/cms/tracking/output/";
129  // plotfilename += dirname;
130  // plotfilename += "/";
131  // plotfilename += labfull;
132  plotfilename += "/deltaxsum_";
133  plotfilename += labfull;
134  // plotfilename += "_";
135  // plotfilename += dirname;
136  plotfilename += ".gif";
137 
138  /* TCanvas * cwidedeltax = */ new TCanvas(plotfilename.c_str(),plotfilename.c_str(),1500,500);
139 
140  deltaxsum->SetLineColor(kRed); deltaxsum->SetMarkerColor(kRed);
141  deltaxsum->GetYaxis()->SetRangeUser(-.002,.002);
142  deltaxsum->GetYaxis()->SetTitle("#Delta x (cm)");
143  deltaxsum->Draw();
144  deltaxmeansum->Draw("esame");
145  TLegend deltaxleg(.7,.8,.85,.9,"#Delta(x)");
146  deltaxleg.AddEntry(deltaxsum,"fitted mean","l");
147  deltaxleg.AddEntry(deltaxmeansum,"aritm. mean","l");
148  deltaxleg.Draw();
149  gPad->Print(plotfilename.c_str());
150  // delete cwidedeltax;
151 
152  plotfilename = "/afs/cern.ch/cms/tracking/output/";
153  // plotfilename += dirname;
154  // plotfilename += "/";
155  // plotfilename += labfull;
156  plotfilename += "/deltaysum_";
157  plotfilename += labfull;
158  // plotfilename += "_";
159  // plotfilename += dirname;
160  plotfilename += ".gif";
161 
162  /* TCanvas * cwidedeltay = */ new TCanvas(plotfilename.c_str(),plotfilename.c_str(),1500,500);
163 
164  deltaysum->SetLineColor(kRed); deltaysum->SetMarkerColor(kRed);
165  deltaysum->GetYaxis()->SetRangeUser(-.002,.002);
166  deltaysum->GetYaxis()->SetTitle("#Delta y (cm)");
167  deltaysum->Draw();
168  deltaymeansum->Draw("esame");
169  TLegend deltayleg(.7,.8,.85,.9,"#Delta(y)");
170  deltayleg.AddEntry(deltaysum,"fitted mean","l");
171  deltayleg.AddEntry(deltaymeansum,"aritm. mean","l");
172  deltayleg.Draw();
173  gPad->Print(plotfilename.c_str());
174  // delete cwidedeltay;
175 
176 
177  plotfilename = "/afs/cern.ch/cms/tracking/output/";
178  // plotfilename += dirname;
179  // plotfilename += "/";
180  // plotfilename += labfull;
181  plotfilename += "/deltazsum_";
182  plotfilename += labfull;
183  // plotfilename += "_";
184  // plotfilename += dirname;
185  plotfilename += ".gif";
186 
187  /* TCanvas * cwidedeltaz = */ new TCanvas(plotfilename.c_str(),plotfilename.c_str(),1500,500);
188 
189  deltazmeansum->GetYaxis()->SetRangeUser(-2.,2.);
190  deltazmeansum->GetYaxis()->SetTitle("#Delta z (cm)");
191  deltazmeansum->Draw();
192  gPad->Print(plotfilename.c_str());
193  // delete cwidedeltaz;
194  /*
195  delete deltaxsum;
196  delete deltaysum;
197  delete deltaxmeansum;
198  delete deltaymeansum;
199  delete deltazmeansum;
200  */
201  delete fdoubleg;
202 
203 }
i
int i
Definition: DBlmapReader.cc:9
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
gather_cfg.runs
tuple runs
Definition: gather_cfg.py:87
compare_using_db.ifile
tuple ifile
Definition: compare_using_db.py:250
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
module
Definition: vlib.h:208
Generated for CMSSW Reference Manual by doxygen 1.8.5