ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/Systematics.C (file contents):
Revision 1.29 by buchmann, Wed Sep 7 06:42:21 2011 UTC vs.
Revision 1.38 by buchmann, Wed Sep 28 16:00:49 2011 UTC

#	Line 1 | Line 1
1		#include <iostream>
2		#include <vector>
3		#include <sys/stat.h>
4	+	#include <algorithm>
5	+	#include <cmath>
6
7		#include <TMath.h>
8		#include <TColor.h>
#	Line 276 | Line 278 | void PeakError(TTree *events,float &resu
278
279		//____________________________________________________________________________________
280		// Total selection efficiency (MC)
281	<	void MCefficiency(TTree *events,float &result, float &resulterr,string mcjzb,bool requireZ,int Neventsinfile, string addcut="", int k = 0) {
281	>	//returns the efficiency WITHOUT signal contamination, and the result and resulterr contain the result and the corresponding error
282	>	Value MCefficiency(TTree *events,float &result, float &resulterr,string mcjzb,bool requireZ,int Neventsinfile, string addcut="", int k = 0) {
283	>
284	>	if(!events) {
285	>	write_error(__FUNCTION__,"Tree passed for efficiency calculation is invalid!");
286	>	return Value(0,0);
287	>	}
288
289		char jzbSelStr[256]; sprintf(jzbSelStr,"%f",jzbSel);
290		// All acceptance cuts at gen. level
#	Line 317 | Line 325 | void MCefficiency(TTree *events,float &r
325		//      Float_t tot = events->GetSelectedRows();
326		Float_t tot = Neventsinfile;
327
328	+	Value result_wo_signalcont;
329	+
330		if(ConsiderSignalContaminationForLimits) {
331		result=(sel-nsel)/tot;
332		resulterr=(1.0/tot)*TMath::Sqrt(sel+nsel+(sel-nsel)*(sel-nsel)/tot);
333	+	result_wo_signalcont=Value(sel/tot,TMath::Sqrt(sel/tot*(1+sel/tot)/tot));
334		} else {//no signal contamination considered:
335		result=(sel)/tot;
336		resulterr=TMath::Sqrt(sel/tot*(1+sel/tot)/tot);
337	+	result_wo_signalcont=Value(result,resulterr);
338		}
339	+	if(!automatized && k>0 ) dout << "PDF assessment: ";
340		if(!automatized) dout << "  MC efficiency: " << result << "+-" << resulterr << "  ( JZB>" << jzbSel << " : " << sel << " , JZB<-" << jzbSel << " : " << nsel << " and nevents=" << tot << ") with normFact=" << normFactor << std::endl;
341		delete effh;
342	+	return result_wo_signalcont;
343		}
344
345
346		//____________________________________________________________________________________
347	<	// Total selection efficiency (MC)
347	>	// Selection efficiency for one process (MC)
348		vector<float> processMCefficiency(TTree *events,string mcjzb,bool requireZ,int Neventsinfile, string addcut) {
349		vector<float> process_efficiencies;
350		for(int iprocess=0;iprocess<=10;iprocess++) {
#	Line 412 | Line 426 | void doJZBscale(TTree *events, float &do
426		Float_t eff  = Interpolate(jzbSel,hist);
427		Float_t effp = Interpolate(jzbSel*(1.+systematic),hist);
428		Float_t effm = Interpolate(jzbSel*(1.-systematic),hist);
429	<	if(!automatized) dout << "  efficiency at JZB==" << jzbSel*(1.+systematic)  << "(-"<<syst*100<<"%) : " << effp << " (" << ((effp-eff)/eff)*100. << "%)"  << std::endl;
429	>	if(!automatized) dout << "  efficiency at JZB==" << jzbSel*(1.+systematic)  << "(-"<<systematic*100<<"%) : " << effp << " (" << ((effp-eff)/eff)*100. << "%)"  << std::endl;
430		if(!automatized) dout << "  efficiency at JZB==" << jzbSel  << ": " << eff << std::endl;
431	<	if(!automatized) dout << "  efficiency at JZB==" << jzbSel*(1.-systematic)  << "(-"<<syst*100<<"%) : " << effm << " (" << ((effm-eff)/eff)*100. << "%)"  << std::endl;
431	>	if(!automatized) dout << "  efficiency at JZB==" << jzbSel*(1.-systematic)  << "(-"<<systematic*100<<"%) : " << effm << " (" << ((effm-eff)/eff)*100. << "%)"  << std::endl;
432		up=((effp-eff)/eff);
433		down=((effm-eff)/eff);
434		}
#	Line 438 | Line 452 | void JZBresponse(TTree *events, bool req
452		hJzbResp->SetMinimum(0.2);
453		hJzbResp->Fit("pol0","Q");
454		TF1 *fittedfunction = hJzbResp->GetFunction("pol0");
455	<	resp=fittedfunction->GetParameter(0);
456	<	resperr=fittedfunction->GetParError(0);
457	<	if(!automatized) dout << "  Response: " << resp << " +/- " << resperr << endl;
455	>	if(!fittedfunction) {
456	>	// in case there are not enough points passing our selection
457	>	cout << "OOPS response function invalid, assuming 100% error !!!!" << endl;
458	>	resp=1;
459	>	resperr=1;
460	>	} else {
461	>	resp=fittedfunction->GetParameter(0);
462	>	resperr=fittedfunction->GetParError(0);
463	>	if(!automatized) dout << "  Response: " << resp << " +/- " << resperr << endl;
464	>	}
465		delete hJzbResp;
466		}
467
#	Line 478 | Line 499 | void do_systematics_for_one_file(TTree *
499		float triggereff=5.0/100;// in range [0,1]
500		dout << "Trigger efficiency not implemented in this script  yet, still using external one" << endl;
501		float leptonseleff=2.0/100;// in range [0,1]
502	+	leptonseleff=TMath::Sqrt(leptonseleff*leptonseleff+leptonseleff*leptonseleff); // because the 2% is per lepton
503		dout << "Lepton selection efficiency not implemented in this script  yet, still using external one" << endl;
504
505		int NPdfs=0;
#	Line 485 | Line 507 | void do_systematics_for_one_file(TTree *
507
508		float mceff,mcefferr,jzbeff,jzbefferr;
509		if(!automatized) dout << "MC efficiencies:" << endl;
510	<	MCefficiency(events,mceff,mcefferr,mcjzb,requireZ,Neventsinfile,addcut,-1);
511	<	JZBefficiency(events,informalname,jzbeff,jzbefferr,requireZ,addcut);
510	>	Value mceff_nosigcont = MCefficiency(events,mceff,mcefferr,mcjzb,requireZ,Neventsinfile,addcut,-1);
511	>	if(!automatized) cout << "   Without signal contamination, we find an efficiency of " << mceff_nosigcont << endl;
512	>
513	>	if(PlottingSetup::computeJZBefficiency) JZBefficiency(events,informalname,jzbeff,jzbefferr,requireZ,addcut);
514		if(!automatized) dout << "JZB efficiency: " << jzbeff << "+/-" << jzbefferr << endl;
515
516		if(!automatized) dout << "Error from Peak position:" << endl;
#	Line 503 | Line 527 | void do_systematics_for_one_file(TTree *
527
528		if(!automatized) dout << "JZB response: " << std::endl;
529		float resp,resperr;
530	<	JZBresponse(events,requireZ,resp,resperr,addcut);
530	>	if(PlottingSetup::computeJZBresponse) {
531	>	if(!automatized) dout << "JZB response: " << std::endl;
532	>	JZBresponse(events,requireZ,resp,resperr,addcut);
533	>	}
534
535		if(!automatized) dout << "Pileup: " << std::endl;
536	<	float resolution=pileup(events,requireZ,informalname,addcut);
536	>	float resolution;
537	>	resolution=pileup(events,requireZ,informalname,addcut);
538
539		float PDFuncert=0;
540	+	if(!automatized) dout << "Assessing PDF uncertainty: " << std::endl;
541		if(ismSUGRA) PDFuncert = get_pdf_uncertainty(events, mcjzb, requireZ, Neventsinfile, NPdfs, addcut);
542
543		dout << "_______________________________________________" << endl;
#	Line 523 | Line 552 | void do_systematics_for_one_file(TTree *
552		dout << "Resolution : " << resolution << endl; // in range [0,1]
553		dout << "From peak : " << sysfrompeak << endl; // in range [0,1]
554		if(ismSUGRA) dout << "PDF uncertainty  : " << PDFuncert << endl; // in range [0,1]
555	<	dout << "JZB efficiency: " << jzbeff << "+/-" << jzbefferr << " (not yet included below) " << endl; // in range [0,1]
556	<	dout << "JZB response  : " << resp << " +/-" << resperr << " (not yet included below) " << endl; // in range [0,1]
555	>	if(PlottingSetup::computeJZBefficiency) dout << "JZB efficiency: " << jzbeff << "+/-" << jzbefferr << " (not yet included below) " << endl; // in range [0,1]
556	>	if(PlottingSetup::computeJZBresponse)dout << "JZB response  : " << resp << " +/-" << resperr << " (not yet included below) " << endl; // in range [0,1]
557
558		float toterr=0;
559		toterr+=(triggereff)*(triggereff);
#	Line 540 | Line 569 | void do_systematics_for_one_file(TTree *
569
570		dout << "FINAL RESULT : " << 100*mceff << " +/- "<< 100*mcefferr << " (stat) +/- " << 100*systerr << " (syst)   %" << endl;
571		dout << "     we thus use the sqrt of the sum of the squares of the stat & syst err, which is : " << 100*toterr << endl;
572	+	dout << "_______________________________________________" << endl;
573
574		//Do not modify the lines below or mess with the order; this order is expected by all limit calculating functions!
575		vector<float> res;
#	Line 551 | Line 581 | void do_systematics_for_one_file(TTree *
581		if(fabs(jesup)>fabs(jesdown)) res.push_back(fabs(jesup)); else res.push_back(fabs(jesdown));
582		if(fabs(scaleup)>fabs(scaledown)) res.push_back(fabs(scaleup)); else res.push_back(fabs(scaledown));
583		res.push_back(fabs(resolution));
584	+	res.push_back(mceff_nosigcont.getValue());
585	+	res.push_back(mceff_nosigcont.getError());
586		if(ismSUGRA) res.push_back(PDFuncert);
587		results.push_back(res);
588		}

Diff Legend

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