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.26 by buchmann, Thu Sep 1 09:27:46 2011 UTC vs.
Revision 1.29 by buchmann, Wed Sep 7 06:42:21 2011 UTC

#	Line 328 | Line 328 | void MCefficiency(TTree *events,float &r
328		delete effh;
329		}
330
331	+
332	+	//____________________________________________________________________________________
333	+	// Total selection efficiency (MC)
334	+	vector<float> processMCefficiency(TTree *events,string mcjzb,bool requireZ,int Neventsinfile, string addcut) {
335	+	vector<float> process_efficiencies;
336	+	for(int iprocess=0;iprocess<=10;iprocess++) {
337	+	float this_process_efficiency,efferr;
338	+	stringstream addcutplus;
339	+	addcutplus<<addcut<<"&&(process=="<<iprocess<<")";
340	+	MCefficiency(events,this_process_efficiency, efferr,mcjzb,requireZ,Neventsinfile, addcutplus.str(),-1);
341	+	process_efficiencies.push_back(this_process_efficiency);
342	+	}
343	+	return process_efficiencies;
344	+	}
345	+
346	+
347		void JZBefficiency(TTree *events, string informalname, float &jzbeff, float &jzbefferr, bool requireZ, string addcut="") {
348		TCut kbase("abs(genMll-91.2)<20&&genNjets>2&&genZPt>0&&abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
349		if(addcut!="") kbase=kbase&&addcut.c_str();//this is mostly for SUSY scans (adding requirements on masses)
#	Line 459 | Line 475 | void do_systematics_for_one_file(TTree *
475		float JetEnergyScaleUncert=0.1;
476		float JZBScaleUncert=0.1;
477		mcjzbexpression=mcjzb;
478	<	float triggereff=4.0/100;// in range [0,1]
478	>	float triggereff=5.0/100;// in range [0,1]
479		dout << "Trigger efficiency not implemented in this script  yet, still using external one" << endl;
480		float leptonseleff=2.0/100;// in range [0,1]
481		dout << "Lepton selection efficiency not implemented in this script  yet, still using external one" << endl;
482
483	<	int NPdfs=get_npdfs(events);
484	<
483	>	int NPdfs=0;
484	>	if(ismSUGRA) NPdfs = get_npdfs(events);
485	>
486		float mceff,mcefferr,jzbeff,jzbefferr;
487		if(!automatized) dout << "MC efficiencies:" << endl;
488		MCefficiency(events,mceff,mcefferr,mcjzb,requireZ,Neventsinfile,addcut,-1);

Diff Legend

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