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.37 by buchmann, Mon Sep 26 15:43:37 2011 UTC vs.
Revision 1.38 by buchmann, Tue Sep 27 08:09:32 2011 UTC

#	Line 59 | Line 59 | void do_limit_wrapper(float mceff,float
59		int counterinterval=5;
60		sleep(1); //waiting a second for the process to become active
61		while(counter<limitpatience*60 && isThreadActive) {
62	<	std::cout << "Limits are being calculated; Checking round " << counter << " ( corresponds to " << seconds_to_time(counter) << " ) , patience will end in " << seconds_to_time(60*limitpatience-counter) << std::endl;
62	>	std::cout << "Limits are being calculated; Checking round " << counter/counterinterval << " ( corresponds to " << seconds_to_time(counter) << " ) , patience will end in " << seconds_to_time(60*limitpatience-counter) << std::endl;
63		counter+=counterinterval;
64		sleep(counterinterval);
65		}
#	Line 232 | Line 232 | void establish_SUSY_limits(string mcjzb,
232		exp1mlimitmap->Fill(mglu,mlsp,sigmas[3]);
233		exp2plimitmap->Fill(mglu,mlsp,sigmas[4]);
234		exp2mlimitmap->Fill(mglu,mlsp,sigmas[5]);
235	–	} else {
236	–	write_warning(__FUNCTION__,"Watch out, limit results did not contain expected limits!");
235		}
236
237		dout << "An upper limit has been added for this point ( " << massgluname << " = " << mglu << " and " << massLSPname << " = " << mlsp << " ) at " << sigmas[0] << endl;
238	<	if(1) {
238	>	if(mSUGRA) { // for SMS this is a bit easier at the moment - we have a reference XS file which we use when plotting
239		dout << "Computing exclusion status" << endl;
240		float rel_limit=0;
241		float xs=get_xs(mglu,mlsp,massgluname,massLSPname,xsec,mcjzb,requireZ);
#	Line 348 | Line 346 | cout << "FILE USED (name): " << (scansam
346		TH2F *efficiencymap =  new TH2F((prefix+"efficiencymap"+any2string(jzbSel)).c_str(),"",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
347		TH2F *noscefficiencymap = new TH2F((prefix+"nosc_efficiencymap"+any2string(jzbSel)).c_str(),"",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
348		TH2F *Neventsmap    =  new TH2F((prefix+"Neventsmap"+any2string(jzbSel)).c_str(),"",   (mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
349	<	TH2F *ipointmap    =   new TH2F((prefix+"Neventsmap"+any2string(jzbSel)).c_str(),"",   (mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
349	>	TH2F *ipointmap    =   new TH2F((prefix+"ipointmap"+any2string(jzbSel)).c_str(),"",   (mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
350		TH2F *syspdfmap =      new TH2F((prefix+"syspdfmap"+any2string(jzbSel)).c_str(),"",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
351		TH2F *systotmap =      new TH2F((prefix+"systotmap"+any2string(jzbSel)).c_str(),"",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
352		TH2F *sysstatmap =     new TH2F((prefix+"sysstatmap"+any2string(jzbSel)).c_str(),"",(mgluend-mglustart)/mglustep+1,mglustart-0.5*mglustep,mgluend+0.5*mglustep,(mLSPend-mLSPstart)/mLSPstep+1,mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
#	Line 437 | Line 435 | write_warning(__FUNCTION__,"CURRENTLY SW
435		exp1mlimitmap->Fill(mglu,mlsp,sigmas[3]);
436		exp2plimitmap->Fill(mglu,mlsp,sigmas[4]);
437		exp2mlimitmap->Fill(mglu,mlsp,sigmas[5]);
438	<	} else {
441	<	write_warning(__FUNCTION__,"Watch out, limit results did not contain expected limits!");
442	<	}
438	>	}
439
440		sysjesmap->Fill(mglu,mlsp,sys_jes);
441		sysjsumap->Fill(mglu,mlsp,sys_jsu);

Diff Legend

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