ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/Modules/SUSYScan.C

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/SUSYScan.C (file contents):
Revision 1.1 by buchmann, Fri Jul 22 10:15:34 2011 UTC vs.
Revision 1.2 by buchmann, Fri Jul 22 10:25:58 2011 UTC

#	Line 168 | Line 168 | void efficiency_scan_in_susy_space(strin
168		cout << "Doing efficiency scan using " << (scansample.collection)[0].filename << endl;
169		geqleq="geq";
170		automatized=true;
171	+	mcjzbexpression=mcjzb;
172
173		string massgluname="MassGlu";
174		string massLSPname="MassLSP";
#	Line 177 | Line 178 | void efficiency_scan_in_susy_space(strin
178		Double_t g[]    = {0., 0.0, 0.0, 1.0, 1.0};
179		Double_t b[]    = {0., 1.0, 0.0, 0.0, 1.0};
180		Double_t stop[] = {0., .25, .50, .75, 1.0};
181	<	Int_t FI = TColor::CreateGradientColorTable(5, stop, r, g, b, 100);
181	>	Int_t FI = TColor::CreateGradientColorTable(5, stop, r, g, b, 110);
182		for (int i=0;i<100;i++) MyPalette[i] = FI+i;
183
184		gStyle->SetPalette(100, MyPalette);
#	Line 189 | Line 190 | void efficiency_scan_in_susy_space(strin
190
191
192		TH2F *efficiencymap = new TH2F("efficiencymap","",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
193	+	TH2F *Neventsmap    = new TH2F("Neventsmap","",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
194	+
195		float rightmargin=gStyle->GetPadRightMargin();
196	<	gStyle->SetPadRightMargin(0.12);
196	>	gStyle->SetPadRightMargin(0.14);
197	>
198		TCanvas *effcanvas = new TCanvas("effcanvas","Efficiency canvas");
199	+
200		for(int mglu=mglustart;mglu<=mgluend;mglu+=mglustep) {
201		for (int mlsp=mLSPstart;mlsp<=mLSPend&&mlsp<=mglu;mlsp+=mLSPstep)
202		{
#	Line 204 | Line 209 | void efficiency_scan_in_susy_space(strin
209		result=0; // kicking nan and inf errors
210		}
211		efficiencymap->Fill(mglu,mlsp,result);
212	+
213	+	(scansample.collection)[0].events->Draw(mcjzb.c_str(),addcut.str().c_str(),"goff");
214	+	float nevents = (scansample.collection)[0].events->GetSelectedRows();
215	+	Neventsmap->Fill(mglu,mlsp,nevents);
216		}
217		}
218
#	Line 212 | Line 221 | void efficiency_scan_in_susy_space(strin
221		efficiencymap->GetYaxis()->SetTitle("m_{LSP}");
222		efficiencymap->GetYaxis()->CenterTitle();
223
224	+	Neventsmap->GetXaxis()->SetTitle("m_{glu}");
225	+	Neventsmap->GetXaxis()->CenterTitle();
226	+	Neventsmap->GetYaxis()->SetTitle("m_{LSP}");
227	+	Neventsmap->GetYaxis()->CenterTitle();
228
229		efficiencymap->Draw("COLZ");
230		TText *title = write_title("Efficiency in LSP-Glu plane");
231		title->Draw();
232		CompleteSave(effcanvas,"SUSYScan/Efficency");
233	+	Neventsmap->Draw("COLZ");
234	+	TText *title2 = write_title("Number of events in LSP-Glu plane");
235	+	title2->Draw();
236	+	CompleteSave(effcanvas,"SUSYScan/Nevents");
237		gStyle->SetPadRightMargin(rightmargin);
238		}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines