| 16 |
|
#include <TProfile.h> |
| 17 |
|
#include <TPaveStats.h> |
| 18 |
|
|
| 19 |
+ |
#include "GeneratorLevelStudyModule.C" |
| 20 |
+ |
|
| 21 |
|
//#include "TTbar_stuff.C" |
| 22 |
|
using namespace std; |
| 23 |
|
|
| 761 |
|
} |
| 762 |
|
|
| 763 |
|
|
| 764 |
+ |
void generator_study_plots() { |
| 765 |
+ |
|
| 766 |
+ |
//uncomment whichever one you want to see :-) |
| 767 |
+ |
|
| 768 |
+ |
|
| 769 |
+ |
GenLevelStudy::X_vs_generation_lm4(); |
| 770 |
+ |
/* |
| 771 |
+ |
GenLevelStudy::compare_sms_lm4(); |
| 772 |
+ |
GenLevelStudy::MomentumFraction(); |
| 773 |
+ |
GenLevelStudy::AngleMETsumLSP(); |
| 774 |
+ |
GenLevelStudy::RatioMETsumLSP(); |
| 775 |
+ |
GenLevelStudy::AngleLSPLSP(); |
| 776 |
+ |
GenLevelStudy::AngleLSPZ(); |
| 777 |
+ |
GenLevelStudy::WidthIllustration(); |
| 778 |
+ |
GenLevelStudy::ZDecayIllustration(); |
| 779 |
+ |
GenLevelStudy::DeltaLSPmomentum(); |
| 780 |
+ |
GenLevelStudy::pStarIllustration(0.5); |
| 781 |
+ |
GenLevelStudy::pStarIllustration(0.25); |
| 782 |
+ |
GenLevelStudy::pStarIllustration(0.75); |
| 783 |
+ |
GenLevelStudy::pStarDistributions(); |
| 784 |
+ |
GenLevelStudy::DrawJetBand(); |
| 785 |
+ |
GenLevelStudy::DrawOnly100to150inJetBand(); |
| 786 |
+ |
GenLevelStudy::ImpactOfGluinoChi2MassDifference(); |
| 787 |
+ |
*/ |
| 788 |
+ |
} |
| 789 |
+ |
|
| 790 |
+ |
|
| 791 |
|
void jzb_negative_generator_study() { |
| 792 |
|
write_warning(__FUNCTION__,"We are going to use a t5zz scan file, and \033[1;31m WON'T \033[1;35m cut on MassGlu/MassLSP in order to improve statistics. This is ok for small studies, for a real study you'll need to look at points individually ..."); |
| 793 |
|
|
| 940 |
|
void compare_LM4_and_SMS() { |
| 941 |
|
TFile *f1 = new TFile("/shome/lbaeni/jzb/LM4_SMS/SMS_LM4_JZB.root"); |
| 942 |
|
TTree *LM4events = (TTree*)f1->Get("events"); |
| 943 |
< |
TFile *f2 = new TFile("/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_TestingSMS_v5/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root "); |
| 943 |
> |
TFile *f2 = new TFile("/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_TestingSMS_v5/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root"); |
| 944 |
|
TTree *SMSevents = (TTree*)f2->Get("events"); |
| 945 |
|
|
| 946 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "mll",100,50,150,cutOSSF&&cutnJets&&basiccut,"mll",true); |
| 954 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "genMET",100,0.0,600.0,cutOSSF&&basiccut,"genMET",false); |
| 955 |
|
compare_lm4_sms_variable(LM4events, SMSevents, "genNjets",10,-0.5,9.5,cutOSSF&&basiccut,"genNjets",true); |
| 956 |
|
} |
| 957 |
+ |
|
| 958 |
+ |
|
| 959 |
+ |
void compare_onpeak_offpeak_signal_distributions(string mcjzb, string datajzb, int sampleindex) { |
| 960 |
+ |
cout << (signalsamples.collection)[sampleindex].filename << endl; |
| 961 |
+ |
} |
| 962 |
+ |
|
| 963 |
+ |
void compare_onpeak_offpeak_signal_distributions(string mcjzb, string datajzb) { |
| 964 |
+ |
for(int i=0;i<signalsamples.collection.size();i++) compare_onpeak_offpeak_signal_distributions(mcjzb,datajzb,i); |
| 965 |
+ |
} |
| 966 |
+ |
|
