ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/GeneratorLevelStudyModule.C

Comparing UserCode/cbrown/Development/Plotting/Modules/GeneratorLevelStudyModule.C (file contents):
Revision 1.4 by buchmann, Tue Apr 10 13:45:00 2012 UTC vs.
Revision 1.5 by buchmann, Mon Apr 16 17:22:33 2012 UTC

#	Line 50 | Line 50 | namespace GenLevelStudy {
50		void DrawOnly100to150inJetBand();
51		void DrawOnly100to150inJetBand(float massglu,float masslsp);
52		void ImpactOfGluinoChi2MassDifference();
53	+	void color_histos(TH1F *a, TH1F *b, TH1F *c, string xlabel);
54
55		//      string LM4sample="/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_PrimeTime__SAMPLES__Gen_V11/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root";
56		string LM4sample="/scratch/buchmann/buchmann/SignalWithLeptonsStartingAt5GeV/LM4_SUSY_sftsht_7TeV-pythia.root";
57		//      string LM4sample="/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_PrimeTime__SAMPLES__Gen_V16_sumJet_with_LSPsum_LM8only_fixed/LM8_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2_novtx_fixed3.root";
58	<	string SMSsample="/shome/buchmann/JellyfishCBAF/various_assignments/SimPa/Tests/FullSMS.root";
58	>	string SMSsample="/shome/buchmann/JellyfishCBAF/various_assignments/SimPa_LOCAL/Tests/FullSMS.root";
59		//scratch/buchmann/buchmann/GeneratorInformationInJZB___JZBplusSamples_PrimeTime__SAMPLES__Gen_V12_with_zchi2angle_bf1/SMS-T5zz_x-05_Mgluino-150to1200_mLSP-50to1150_7TeV-Pythia6Z__Summer11-PU_START42_V11_FastSim-v2__AODSIM___inclindex_v2_complet.root";
60		string histoname(int i);
61		Color_t histocolor(int i);
#	Line 197 | Line 198 | void GenLevelStudy::MomentumFractionScen
198		stringstream cut;
199		cut << "genNjets>2&&abs(pureGeneratorJZB)>0";
200		if(massglu>0&&masslsp>0) cut << "&&abs(MassGlu-" << massglu << ")<5&&abs(MassLSP-" << masslsp << ")";
201	<	events->Draw("(LSP1pt/LSP1Mopt):pureGeneratorJZB",cut.str().c_str(),"PROF");
202	<	TH1F *h = new TH1F("h","h",100,-500,500);
203	<	h->SetLineColor(kBlack);
204	<	TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
205	<	p->GetXaxis()->SetTitle("Generator JZB");
206	<	p->GetXaxis()->CenterTitle();
201	>
202	>	TProfile *fLSP = new TProfile("fLSP","fLSP",500,-500,1000,0,1.5);
203	>	TProfile *fZ = new TProfile("fZ","fZ",500,-500,1000,0,1.5);
204	>	events->Draw("(LSP1pt/LSP1Mopt):pureGeneratorJZB>>fLSP",cut.str().c_str(),"PROF");
205	>	//  TH1F *h = new TH1F("h","h",100,-500,500);
206	>	//  h->SetLineColor(kBlack);
207	>
208	>	fLSP->GetXaxis()->SetTitle("Generator JZB");
209	>	fLSP->GetXaxis()->CenterTitle();
210	>	fLSP->SetMarkerSize(0);
211		//  p->GetYaxis()->SetTitle("( LSP p_{T} ) / ( LSP mother p_{T} )");
212	<	p->GetYaxis()->SetTitle("f");
213	<	p->GetYaxis()->CenterTitle();
214	<	p->SetLineColor(kBlue);
215	<	p->GetXaxis()->SetRangeUser(-500,1000);
216	<	p->GetYaxis()->SetRangeUser(0,2);
212	>	fLSP->GetYaxis()->SetTitle("f");
213	>	fLSP->GetYaxis()->CenterTitle();
214	>	fLSP->SetLineColor(kBlue);
215	>	fLSP->GetYaxis()->SetRangeUser(0,1.5);
216	>	//  p->GetYaxis()->SetRangeUser(0,2);
217		/*  TLegend* leg = make_legend("", 0.20, 1.03, false);
218		leg->SetNColumns(2);
219		leg->AddEntry(p,"f_{LSP}","l");
#	Line 216 | Line 221 | void GenLevelStudy::MomentumFractionScen
221		leg->Draw();*/
222		//  TText *title = write_title("JZB as a function of the first LSP's momentum transfer");
223		//  title->Draw();
224	<	events->Draw("(genZPt/LSP1Mopt):pureGeneratorJZB","genNjets>2","PROF,same");
224	>	events->Draw("(genZPt/LSP1Mopt):pureGeneratorJZB>>fZ","genNjets>2","PROF,same");
225	>	fZ->SetMarkerSize(0);
226	>	fLSP->Draw();
227	>	fZ->Draw("same");
228
229		stringstream saveas;
230		saveas << "GeneratorStudies/MomentumFraction";
#	Line 295 | Line 303 | void GenLevelStudy::AngleLSPLSPScenario(
303		stringstream cut;
304		cut << "genNjets>2&&abs(pureGeneratorJZB)>0";
305		if(massglu>0&&masslsp>0) cut << "&&abs(MassGlu-" << massglu << ")<5&&abs(MassLSP-" << masslsp << ")";
306	<	events->Draw("angleLSPLSP/TMath::Pi():pureGeneratorJZB",cut.str().c_str(),"PROF");
307	<	TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
308	<	p->GetXaxis()->SetTitle("Generator JZB");
309	<	p->GetXaxis()->CenterTitle();
310	<	p->GetYaxis()->SetTitle("#angle (LSP,LSP) / #pi");
311	<	p->GetYaxis()->CenterTitle();
312	<	p->SetLineColor(kBlue);
313	<	p->GetXaxis()->SetRangeUser(-500,1000);
314	<	p->GetYaxis()->SetRangeUser(0,1);
306	>	TProfile *angLSP = new TProfile("angLSP","angLSP",500,-500,1000,0,3.2);
307	>	events->Draw("angleLSPLSP:pureGeneratorJZB>>angLSP",cut.str().c_str(),"PROF");
308	>	//  TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
309	>	angLSP->GetXaxis()->SetTitle("Generator JZB");
310	>	angLSP->GetXaxis()->CenterTitle();
311	>	angLSP->GetYaxis()->SetTitle("#angle (LSP,LSP)");
312	>	angLSP->GetYaxis()->CenterTitle();
313	>	angLSP->SetLineColor(kBlue);
314	>	angLSP->SetMarkerSize(0);
315	>	//  p->GetXaxis()->SetRangeUser(-500,1000);
316	>	angLSP->GetYaxis()->SetRangeUser(0,3.1415);
317	>
318		stringstream saveas;
319		saveas << "GeneratorStudies/AngleLSPLSP";
320		if(massglu>0&&masslsp>0) saveas << "__MassGlu" << massglu << "_MassLSP_" << masslsp;
#	Line 543 | Line 554 | void GenLevelStudy::DeltaLSPmomentumScen
554		stringstream cut;
555		cut << "genNjets>2&&abs(pureGeneratorJZB)>0";
556		if(massglu>0&&masslsp>0) cut << "&&abs(MassGlu-" << massglu << ")<5&&abs(MassLSP-" << masslsp << ")";
557	<	events->Draw("(LSP1pt-LSP2pt):pureGeneratorJZB",cut.str().c_str(),"PROF");
558	<	TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
559	<	p->GetXaxis()->SetTitle("Generator JZB");
560	<	p->GetXaxis()->CenterTitle();
561	<	p->GetYaxis()->SetTitle("p_{T}(LSP1)-p_{T}(LSP2)");
562	<	p->GetYaxis()->CenterTitle();
563	<	p->GetXaxis()->SetRangeUser(-500,1000);
564	<	p->GetYaxis()->SetRangeUser(-800,800);
565	<	p->SetLineColor(kBlue);
557	>	TProfile *DeltaLSP = new TProfile("DeltaLSP","DeltaLSP",500,-500,1000,-800,800);
558	>	events->Draw("(LSP1pt-LSP2pt):pureGeneratorJZB>>DeltaLSP",cut.str().c_str(),"PROF");
559	>	//  TProfile *p = (TProfile*)jcan->GetPrimitive("htemp");
560	>	DeltaLSP->GetXaxis()->SetTitle("Generator JZB");
561	>	DeltaLSP->GetXaxis()->CenterTitle();
562	>	DeltaLSP->GetYaxis()->SetTitle("p_{T}(LSP1)-p_{T}(LSP2)");
563	>	DeltaLSP->GetYaxis()->CenterTitle();
564	>	DeltaLSP->GetXaxis()->SetRangeUser(-500,1000);
565	>	DeltaLSP->GetYaxis()->SetRangeUser(-800,800);
566	>	DeltaLSP->SetLineColor(kBlue);
567	>	DeltaLSP->SetMarkerSize(0);
568
569		stringstream saveas;
570		saveas << "GeneratorStudies/DeltaLSPPt";
#	Line 771 | Line 784 | void GenLevelStudy::DrawJetBand() {
784		GenLevelStudy::DrawJetBand(625,425);
785		}
786
787	+	void GenLevelStudy::color_histos(TH1F *a, TH1F *b, TH1F *c, string xlabel) {
788	+	a->GetXaxis()->SetTitle(xlabel.c_str());
789	+	a->GetYaxis()->SetTitle("events");
790	+	a->GetXaxis()->CenterTitle();
791	+	a->GetYaxis()->CenterTitle();
792	+	b->SetLineColor(diversehistocolor(1));
793	+
794	+	if(!Contains(((string)(a->GetName())),"LSP")) {
795	+	a->SetLineColor(diversehistocolor(0));
796	+	c->SetLineColor(diversehistocolor(2));
797	+	} else {
798	+	a->SetLineColor(diversehistocolor(2));
799	+	c->SetLineColor(diversehistocolor(0));
800	+	}
801	+	}
802	+
803	+
804		void GenLevelStudy::ImpactOfGluinoChi2MassDifference() {
805	+
806	+	// STEP 1 : Loading everything
807		TCanvas *can = new TCanvas("can","can");
808		can->SetLogy(1);
809		TChain smsp025events("events");
#	Line 783 | Line 815 | void GenLevelStudy::ImpactOfGluinoChi2Ma
815		stringstream filenamep25;
816		stringstream filenamep50;
817		stringstream filenamep75;
818	<	filenamep25 << "/shome/buchmann/ntuples/SMS_T5zzh_newslots/SMS_clean_splitup_" << i << "_" << j << ".root";
818	>	filenamep25 << "/shome/lbaeni/jzb/SMS_T5zzh_newslots/SMS_clean_splitup_" << i << "_" << j << ".root";
819		filenamep50 << "/shome/buchmann/ntuples/SMS_T5zz/SMS_clean_splitup_" << i << "_" << j << ".root";
820		filenamep75 << "/shome/buchmann/ntuples/SMS_T5zzl/SMS_clean_splitup_" << i << "_" << j << ".root";
821		smsp025events.Add(filenamep25.str().c_str());
#	Line 792 | Line 824 | void GenLevelStudy::ImpactOfGluinoChi2Ma
824		}
825		}
826
795	–	stringstream cutp25m2;
796	–	cutp25m2 << "abs(MassGlu-700)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
797	–	stringstream cutp50m2;
798	–	cutp50m2 << "abs(MassGlu-400)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
799	–	stringstream cutp75m2;
800	–	cutp75m2 << "abs(MassGlu-300)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
801	–
802	–	stringstream cutp25m1;
803	–	cutp25m1 << "abs(MassGlu-700)<5&&abs(MassLSP-500)&&pfJetGoodNum>0";
804	–	stringstream cutp50m1;
805	–	cutp50m1 << "abs(MassGlu-700)<5&&abs(MassLSP-400)&&pfJetGoodNum>0";
806	–	stringstream cutp75m1;
807	–	cutp75m1 << "abs(MassGlu-700)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
808	–
809	–	TH1F *nJets25m2 = new TH1F("nJets25m2","nJets25m2",20,-0.5,19.5);
810	–	TH1F *nJets50m2 = new TH1F("nJets50m2","nJets50m2",20,-0.5,19.5);
811	–	TH1F *nJets75m2 = new TH1F("nJets75m2","nJets75m2",20,-0.5,19.5);
812	–
813	–	TH1F *sumJetpt25m2 = new TH1F("sumJetpt25m2","sumJetpt25m2",100,0,300);
814	–	TH1F *sumJetpt50m2 = new TH1F("sumJetpt50m2","sumJetpt50m2",100,0,300);
815	–	TH1F *sumJetpt75m2 = new TH1F("sumJetpt75m2","sumJetpt75m2",100,0,300);
816	–
817	–	TH1F *Zpt25m2 = new TH1F("Zpt25m2","Zpt25m2",100,0,200);
818	–	TH1F *Zpt50m2 = new TH1F("Zpt50m2","Zpt50m2",100,0,200);
819	–	TH1F *Zpt75m2 = new TH1F("Zpt75m2","Zpt75m2",100,0,200);
820	–
821	–	TH1F *nJets25m1 = new TH1F("nJets25m1","nJets25m1",20,-0.5,19.5);
822	–	TH1F *nJets50m1 = new TH1F("nJets50m1","nJets50m1",20,-0.5,19.5);
823	–	TH1F *nJets75m1 = new TH1F("nJets75m1","nJets75m1",20,-0.5,19.5);
824	–
825	–	TH1F *sumJetpt25m1 = new TH1F("sumJetpt25m1","sumJetpt25m1",100,0,300);
826	–	TH1F *sumJetpt50m1 = new TH1F("sumJetpt50m1","sumJetpt50m1",100,0,300);
827	–	TH1F *sumJetpt75m1 = new TH1F("sumJetpt75m1","sumJetpt75m1",100,0,300);
828	–
829	–	TH1F *Zpt25m1 = new TH1F("Zpt25m1","Zpt25m1",100,0,200);
830	–	TH1F *Zpt50m1 = new TH1F("Zpt50m1","Zpt50m1",100,0,200);
831	–	TH1F *Zpt75m1 = new TH1F("Zpt75m1","Zpt75m1",100,0,200);
832	–
833	–	smsp025events.Draw("pfJetGoodNum>>nJets25m2",("abs(mll-91.2)&&id1==id2&&"+cutp25m2.str()).c_str());
834	–	smsp050events.Draw("pfJetGoodNum>>nJets50m2",("abs(mll-91.2)&&id1==id2&&"+cutp50m2.str()).c_str());
835	–	smsp075events.Draw("pfJetGoodNum>>nJets75m2",("abs(mll-91.2)&&id1==id2&&"+cutp75m2.str()).c_str());
836	–
837	–	smsp025events.Draw("sumJetPt[1]>>sumJetpt25m2",("abs(mll-91.2)&&id1==id2&&"+cutp25m2.str()).c_str());
838	–	smsp050events.Draw("sumJetPt[1]>>sumJetpt50m2",("abs(mll-91.2)&&id1==id2&&"+cutp50m2.str()).c_str());
839	–	smsp075events.Draw("sumJetPt[1]>>sumJetpt75m2",("abs(mll-91.2)&&id1==id2&&"+cutp75m2.str()).c_str());
840	–
841	–	smsp025events.Draw("pt>>Zpt25m2",("abs(mll-91.2)&&id1==id2&&"+cutp25m2.str()).c_str());
842	–	smsp050events.Draw("pt>>Zpt50m2",("abs(mll-91.2)&&id1==id2&&"+cutp50m2.str()).c_str());
843	–	smsp075events.Draw("pt>>Zpt75m2",("abs(mll-91.2)&&id1==id2&&"+cutp75m2.str()).c_str());
844	–
845	–	smsp025events.Draw("pfJetGoodNum>>nJets25m1",("abs(mll-91.2)&&id1==id2&&"+cutp25m1.str()).c_str());
846	–	smsp050events.Draw("pfJetGoodNum>>nJets50m1",("abs(mll-91.2)&&id1==id2&&"+cutp50m1.str()).c_str());
847	–	smsp075events.Draw("pfJetGoodNum>>nJets75m1",("abs(mll-91.2)&&id1==id2&&"+cutp75m1.str()).c_str());
848	–
849	–	smsp025events.Draw("sumJetPt[1]>>sumJetpt25m1",("abs(mll-91.2)&&id1==id2&&"+cutp25m1.str()).c_str());
850	–	smsp050events.Draw("sumJetPt[1]>>sumJetpt50m1",("abs(mll-91.2)&&id1==id2&&"+cutp50m1.str()).c_str());
851	–	smsp075events.Draw("sumJetPt[1]>>sumJetpt75m1",("abs(mll-91.2)&&id1==id2&&"+cutp75m1.str()).c_str());
852	–
853	–	smsp025events.Draw("pt>>Zpt25m1",("abs(mll-91.2)&&id1==id2&&"+cutp25m1.str()).c_str());
854	–	smsp050events.Draw("pt>>Zpt50m1",("abs(mll-91.2)&&id1==id2&&"+cutp50m1.str()).c_str());
855	–	smsp075events.Draw("pt>>Zpt75m1",("abs(mll-91.2)&&id1==id2&&"+cutp75m1.str()).c_str());
856	–
857	–	sumJetpt25m1->GetXaxis()->SetTitle("Sum of Jet Pt (GeV)");
858	–	sumJetpt25m1->GetYaxis()->SetTitle("events");
859	–	sumJetpt25m1->GetXaxis()->CenterTitle();
860	–	sumJetpt25m1->GetYaxis()->CenterTitle();
861	–
862	–	sumJetpt25m1->SetLineColor(diversehistocolor(0));
863	–	sumJetpt50m1->SetLineColor(diversehistocolor(1));
864	–	sumJetpt75m1->SetLineColor(diversehistocolor(2));
865	–
866	–	TLegend *leg = make_legend("",0.6,0.6,false);
867	–	leg->AddEntry(sumJetpt25m1,"#Delta m_{2} = 50 GeV","l");
868	–	leg->AddEntry(sumJetpt50m1,"#Delta m_{2} = 150 GeV","l");
869	–	leg->AddEntry(sumJetpt75m1,"#Delta m_{2} = 450 GeV","l");
870	–
871	–	sumJetpt25m1->Draw("histo");
872	–	sumJetpt50m1->Draw("histo,same");
873	–	sumJetpt75m1->Draw("histo,same");
874	–	leg->Draw("same");
875	–	CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM2/ComparingSumJetPt_Of_Three_Different_Gluino_Chi2_MassDifferences");
876	–
877	–	nJets25m1->SetLineColor(diversehistocolor(0));
878	–	nJets50m1->SetLineColor(diversehistocolor(1));
879	–	nJets75m1->SetLineColor(diversehistocolor(2));
880	–	nJets25m1->GetXaxis()->SetTitle("N(Jets)");
881	–	nJets25m1->GetXaxis()->CenterTitle();
882	–	nJets25m1->GetYaxis()->SetTitle("events");
883	–	nJets25m1->GetYaxis()->CenterTitle();
884	–	nJets25m1->Draw("histo");
885	–	nJets50m1->Draw("histo,same");
886	–	nJets75m1->Draw("histo,same");
887	–	leg->Draw("same");
888	–	CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM2/ComparingNJets_Of_Three_Different_Gluino_Chi2_MassDifferences");
889	–
890	–	Zpt25m1->GetXaxis()->SetTitle("Z p_{T} (GeV)");
891	–	Zpt25m1->GetXaxis()->CenterTitle();
892	–	Zpt25m1->GetYaxis()->SetTitle("events");
893	–	Zpt25m1->GetYaxis()->CenterTitle();
894	–	Zpt25m1->SetLineColor(diversehistocolor(0));
895	–	Zpt50m1->SetLineColor(diversehistocolor(1));
896	–	Zpt75m1->SetLineColor(diversehistocolor(2));
897	–	Zpt25m1->Draw("histo");
898	–	Zpt50m1->Draw("histo,same");
899	–	Zpt75m1->Draw("histo,same");
900	–	leg->Draw("same");
901	–	CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM2/ComparingZPt_Of_Three_Different_Gluino_Chi2_MassDifferences");
902	–
903	–	//------------
904	–
905	–	sumJetpt25m2->GetXaxis()->SetTitle("Sum of Jet Pt (GeV)");
906	–	sumJetpt25m2->GetYaxis()->SetTitle("events");
907	–	sumJetpt25m2->GetXaxis()->CenterTitle();
908	–	sumJetpt25m2->GetYaxis()->CenterTitle();
909	–
910	–	sumJetpt25m2->SetLineColor(diversehistocolor(0));
911	–	sumJetpt50m2->SetLineColor(diversehistocolor(1));
912	–	sumJetpt75m2->SetLineColor(diversehistocolor(2));
913	–
914	–	TLegend *leg2 = make_legend("",0.6,0.6,false);
915	–	leg2->AddEntry(sumJetpt75m2,"#Delta m_{1} = 50 GeV","l");
916	–	leg2->AddEntry(sumJetpt50m2,"#Delta m_{1} = 150 GeV","l");
917	–	leg2->AddEntry(sumJetpt25m2,"#Delta m_{1} = 450 GeV","l");
918	–
919	–	sumJetpt25m2->Draw("histo");
920	–	sumJetpt50m2->Draw("histo,same");
921	–	sumJetpt75m2->Draw("histo,same");
922	–	leg2->Draw("same");
923	–	CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM1/ComparingSumJetPt_Of_Three_Different_Chi1_Chi2_MassDifferences");
924	–
925	–	nJets25m2->SetLineColor(diversehistocolor(0));
926	–	nJets50m2->SetLineColor(diversehistocolor(1));
927	–	nJets75m2->SetLineColor(diversehistocolor(2));
928	–	nJets25m2->GetXaxis()->SetTitle("N(Jets)");
929	–	nJets25m2->GetXaxis()->CenterTitle();
930	–	nJets25m2->GetYaxis()->SetTitle("events");
931	–	nJets25m2->GetYaxis()->CenterTitle();
932	–
933	–	nJets25m2->Draw("histo");
934	–	nJets50m2->Draw("histo,same");
935	–	nJets75m2->Draw("histo,same");
936	–	leg2->Draw("same");
937	–	CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM1/ComparingNJets_Of_Three_Different_Chi1_Chi2_MassDifferences");
938	–
939	–	Zpt25m2->GetXaxis()->SetTitle("Z p_{T} (GeV)");
940	–	Zpt25m2->GetXaxis()->CenterTitle();
941	–	Zpt25m2->GetYaxis()->SetTitle("events");
942	–	Zpt25m2->GetYaxis()->CenterTitle();
943	–
944	–	Zpt25m2->SetLineColor(diversehistocolor(0));
945	–	Zpt50m2->SetLineColor(diversehistocolor(1));
946	–	Zpt75m2->SetLineColor(diversehistocolor(2));
947	–	Zpt25m2->Draw("histo");
948	–	Zpt50m2->Draw("histo,same");
949	–	Zpt75m2->Draw("histo,same");
950	–	leg2->Draw("same");
951	–	CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM1/ComparingZPt_Of_Three_Different_Chi1_Chi2_MassDifferences");
827
828	<	//--------------
954	<	delete can;
828	>	// STEP 2: Define cuts
829
830	<	delete nJets25m1;
831	<	delete nJets50m1;
832	<	delete nJets75m1;
833	<
834	<	delete nJets25m2;
835	<	delete nJets50m2;
836	<	delete nJets75m2;
837	<
838	<	delete sumJetpt25m1;
839	<	delete sumJetpt50m1;
840	<	delete sumJetpt75m1;
841	<
842	<	delete sumJetpt25m2;
843	<	delete sumJetpt50m2;
844	<	delete sumJetpt75m2;
830	>	stringstream cut_VaryingGluiChi2_MDiff_x25;
831	>	cut_VaryingGluiChi2_MDiff_x25 << "abs(MassGlu-700)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";   // config: (700,250,100)
832	>	stringstream cut_VaryingGluiChi2_MDiff_x50;
833	>	cut_VaryingGluiChi2_MDiff_x50  << "abs(MassGlu-400)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";  // config: (400,250,100)
834	>	stringstream cut_VaryingGluiChi2_MDiff_x75;
835	>	cut_VaryingGluiChi2_MDiff_x75  << "abs(MassGlu-300)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";  // config: (300,250,100)
836	>
837	>
838	>	stringstream cut_VaryingChi2LSP_MDiff_x25;
839	>	cut_VaryingChi2LSP_MDiff_x25 << "abs(MassGlu-700)<5&&abs(MassLSP-500)&&pfJetGoodNum>0";   // config: (700,550,500)
840	>	stringstream cut_VaryingChi2LSP_MDiff_x50;
841	>	cut_VaryingChi2LSP_MDiff_x50  << "abs(MassGlu-700)<5&&abs(MassLSP-550)&&pfJetGoodNum>0";  // config: (700,550,400)
842	>	stringstream cut_VaryingChi2LSP_MDiff_x75;
843	>	cut_VaryingChi2LSP_MDiff_x75  << "abs(MassGlu-700)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";  // config: (700,550,100)
844	>
845	>
846	>	// STEP 3: Define histos
847	>	TH1F *nJets25GluChi2 = new TH1F("nJets25GluChi2","nJets25GluChi2",20,-0.5,19.5);
848	>	TH1F *nJets50GluChi2 = new TH1F("nJets50GluChi2","nJets50GluChi2",20,-0.5,19.5);
849	>	TH1F *nJets75GluChi2 = new TH1F("nJets75GluChi2","nJets75GluChi2",20,-0.5,19.5);
850	>
851	>	TH1F *sumJetpt25GluChi2 = new TH1F("sumJetpt25GluChi2","sumJetpt25GluChi2",100,0,300);
852	>	TH1F *sumJetpt50GluChi2 = new TH1F("sumJetpt50GluChi2","sumJetpt50GluChi2",100,0,300);
853	>	TH1F *sumJetpt75GluChi2 = new TH1F("sumJetpt75GluChi2","sumJetpt75GluChi2",100,0,300);
854	>
855	>	TH1F *Zpt25GluChi2 = new TH1F("Zpt25GluChi2","Zpt25GluChi2",100,0,200);
856	>	TH1F *Zpt50GluChi2 = new TH1F("Zpt50GluChi2","Zpt50GluChi2",100,0,200);
857	>	TH1F *Zpt75GluChi2 = new TH1F("Zpt75GluChi2","Zpt75GluChi2",100,0,200);
858	>
859	>	TH1F *MET25GluChi2 = new TH1F("MET25GluChi2","MET25GluChi2",100,0,500);
860	>	TH1F *MET50GluChi2 = new TH1F("MET50GluChi2","MET50GluChi2",100,0,500);
861	>	TH1F *MET75GluChi2 = new TH1F("MET75GluChi2","MET75GluChi2",100,0,500);
862	>
863	>	TH1F *JZB25GluChi2 = new TH1F("JZB25GluChi2","JZB25GluChi2",100,-400,800);
864	>	TH1F *JZB50GluChi2 = new TH1F("JZB50GluChi2","JZB50GluChi2",100,-400,800);
865	>	TH1F *JZB75GluChi2 = new TH1F("JZB75GluChi2","JZB75GluChi2",100,-400,800);
866	>
867	>	TH1F *nJets25Chi2LSP = new TH1F("nJets25Chi2LSP","nJets25Chi2LSP",20,-0.5,19.5);
868	>	TH1F *nJets50Chi2LSP = new TH1F("nJets50Chi2LSP","nJets50Chi2LSP",20,-0.5,19.5);
869	>	TH1F *nJets75Chi2LSP = new TH1F("nJets75Chi2LSP","nJets75Chi2LSP",20,-0.5,19.5);
870	>
871	>	TH1F *sumJetpt25Chi2LSP = new TH1F("sumJetpt25Chi2LSP","sumJetpt25Chi2LSP",100,0,300);
872	>	TH1F *sumJetpt50Chi2LSP = new TH1F("sumJetpt50Chi2LSP","sumJetpt50Chi2LSP",100,0,300);
873	>	TH1F *sumJetpt75Chi2LSP = new TH1F("sumJetpt75Chi2LSP","sumJetpt75Chi2LSP",100,0,300);
874	>
875	>	TH1F *Zpt25Chi2LSP = new TH1F("Zpt25Chi2LSP","Zpt25Chi2LSP",100,0,200);
876	>	TH1F *Zpt50Chi2LSP = new TH1F("Zpt50Chi2LSP","Zpt50Chi2LSP",100,0,200);
877	>	TH1F *Zpt75Chi2LSP = new TH1F("Zpt75Chi2LSP","Zpt75Chi2LSP",100,0,200);
878	>
879	>	TH1F *MET25Chi2LSP = new TH1F("MET25Chi2LSP","MET25Chi2LSP",100,0,500);
880	>	TH1F *MET50Chi2LSP = new TH1F("MET50Chi2LSP","MET50Chi2LSP",100,0,500);
881	>	TH1F *MET75Chi2LSP = new TH1F("MET75Chi2LSP","MET75Chi2LSP",100,0,500);
882	>
883	>	TH1F *JZB25Chi2LSP = new TH1F("JZB25Chi2LSP","JZB25Chi2LSP",100,-400,800);
884	>	TH1F *JZB50Chi2LSP = new TH1F("JZB50Chi2LSP","JZB50Chi2LSP",100,-400,800);
885	>	TH1F *JZB75Chi2LSP = new TH1F("JZB75Chi2LSP","JZB75Chi2LSP",100,-400,800);
886	>
887	>	// STEP 4: Filling the histos
888	>	smsp025events.Draw("pfJetGoodNum>>nJets25GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x25.str()).c_str());
889	>	smsp050events.Draw("pfJetGoodNum>>nJets50GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x50.str()).c_str());
890	>	smsp075events.Draw("pfJetGoodNum>>nJets75GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x75.str()).c_str());
891	>
892	>	smsp025events.Draw("pfJetGoodNum>>nJets25Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x25.str()).c_str());
893	>	smsp050events.Draw("pfJetGoodNum>>nJets50Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x50.str()).c_str());
894	>	smsp075events.Draw("pfJetGoodNum>>nJets75Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x75.str()).c_str());
895	>
896	>	smsp025events.Draw("sumJetPt[1]>>sumJetpt25GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x25.str()).c_str());
897	>	smsp050events.Draw("sumJetPt[1]>>sumJetpt50GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x50.str()).c_str());
898	>	smsp075events.Draw("sumJetPt[1]>>sumJetpt75GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x75.str()).c_str());
899	>
900	>	smsp025events.Draw("sumJetPt[1]>>sumJetpt25Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x25.str()).c_str());
901	>	smsp050events.Draw("sumJetPt[1]>>sumJetpt50Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x50.str()).c_str());
902	>	smsp075events.Draw("sumJetPt[1]>>sumJetpt75Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x75.str()).c_str());
903	>
904	>	smsp025events.Draw("pt>>Zpt25GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x25.str()).c_str());
905	>	smsp050events.Draw("pt>>Zpt50GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x50.str()).c_str());
906	>	smsp075events.Draw("pt>>Zpt75GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x75.str()).c_str());
907	>
908	>	smsp025events.Draw("pt>>Zpt25Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x25.str()).c_str());
909	>	smsp050events.Draw("pt>>Zpt50Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x50.str()).c_str());
910	>	smsp075events.Draw("pt>>Zpt75Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x75.str()).c_str());
911	>
912	>	smsp025events.Draw("met[4]>>MET25GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x25.str()).c_str());
913	>	smsp050events.Draw("met[4]>>MET50GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x50.str()).c_str());
914	>	smsp075events.Draw("met[4]>>MET75GluChi2",("id1==id2&&"+cut_VaryingGluiChi2_MDiff_x75.str()).c_str());
915	>
916	>	smsp025events.Draw("met[4]>>MET25Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x25.str()).c_str());
917	>	smsp050events.Draw("met[4]>>MET50Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x50.str()).c_str());
918	>	smsp075events.Draw("met[4]>>MET75Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x75.str()).c_str());
919	>
920	>	smsp025events.Draw("(((jzb[1]+0.043*pt)-3.05933))>>JZB25GluChi2",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x25.str()).c_str());
921	>	smsp050events.Draw("(((jzb[1]+0.043*pt)-3.05933))>>JZB50GluChi2",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x50.str()).c_str());
922	>	smsp075events.Draw("(((jzb[1]+0.043*pt)-3.05933))>>JZB75GluChi2",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x75.str()).c_str());
923	>
924	>	smsp025events.Draw("(((jzb[1]+0.043*pt)-3.05933))>>JZB25Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x25.str()).c_str());
925	>	smsp050events.Draw("(((jzb[1]+0.043*pt)-3.05933))>>JZB50Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x50.str()).c_str());
926	>	smsp075events.Draw("(((jzb[1]+0.043*pt)-3.05933))>>JZB75Chi2LSP",("id1==id2&&"+cut_VaryingChi2LSP_MDiff_x75.str()).c_str());
927	>
928	>	// STEP 4: Cosmetics
929	>	color_histos(nJets25GluChi2,nJets50GluChi2,nJets75GluChi2,"N(jets)");
930	>	color_histos(nJets25Chi2LSP,nJets50Chi2LSP,nJets75Chi2LSP,"N(jets)");
931	>	color_histos(sumJetpt25GluChi2,sumJetpt50GluChi2,sumJetpt75GluChi2,"Sum of Jet Pt (GeV)");
932	>	color_histos(sumJetpt25Chi2LSP,sumJetpt50Chi2LSP,sumJetpt75Chi2LSP,"Sum of Jet Pt (GeV)");
933	>	color_histos(Zpt25GluChi2,Zpt50GluChi2,Zpt75GluChi2,"Z p_{T} (GeV)");
934	>	color_histos(Zpt25Chi2LSP,Zpt50Chi2LSP,Zpt75Chi2LSP,"Z p_{T} (GeV)");
935	>	color_histos(MET25GluChi2,MET50GluChi2,MET75GluChi2,"Missing E_{T} (GeV)");
936	>	color_histos(MET25Chi2LSP,MET50Chi2LSP,MET75Chi2LSP,"Missing E_{T} (GeV)");
937	>	color_histos(JZB25GluChi2,JZB50GluChi2,JZB75GluChi2,"JZB (GeV)");
938	>	color_histos(JZB25Chi2LSP,JZB50Chi2LSP,JZB75Chi2LSP,"JZB (GeV)");
939	>
940	>	TLegend *legGluChi2 = make_legend("",0.55,0.6,false);
941	>	legGluChi2->AddEntry(nJets75GluChi2,"m_{#tilde{g}}-m_{#chi_{2}} = 50 GeV","l");
942	>	legGluChi2->AddEntry(nJets50GluChi2,"m_{#tilde{g}}-m_{#chi_{2}} = 150 GeV","l");
943	>	legGluChi2->AddEntry(nJets25GluChi2,"m_{#tilde{g}}-m_{#chi_{2}} = 450 GeV","l");
944	>
945	>	TLegend *legChi2LSP = make_legend("",0.55,0.6,false);
946	>	legChi2LSP->AddEntry(nJets25Chi2LSP,"m_{#chi_{2}}-m_{LSP} = 50 GeV","l");
947	>	legChi2LSP->AddEntry(nJets50Chi2LSP,"m_{#chi_{2}}-m_{LSP} = 150 GeV","l");
948	>	legChi2LSP->AddEntry(nJets75Chi2LSP,"m_{#chi_{2}}-m_{LSP} = 450 GeV","l");
949	>
950	>	// STEP 4: Draw!
951	>
952	>	nJets25GluChi2->DrawNormalized("histo");
953	>	nJets50GluChi2->DrawNormalized("histo,same");
954	>	nJets75GluChi2->DrawNormalized("histo,same");
955	>	legGluChi2->Draw("same");
956	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingNJets_Of_Three_Different_Gluino_Chi2_MassDifferences");
957	>
958	>	sumJetpt25GluChi2->SetMaximum(sumJetpt25GluChi2->GetMaximum()*15);
959	>	sumJetpt25GluChi2->DrawNormalized("histo");
960	>	sumJetpt50GluChi2->DrawNormalized("histo,same");
961	>	sumJetpt75GluChi2->DrawNormalized("histo,same");
962	>	legGluChi2->Draw("same");
963	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingSumJetPt_Of_Three_Different_Gluino_Chi2_MassDifferences");
964	>
965	>	Zpt25GluChi2->SetMaximum(Zpt25GluChi2->GetMaximum()*15);
966	>	Zpt25GluChi2->DrawNormalized("histo");
967	>	Zpt50GluChi2->DrawNormalized("histo,same");
968	>	Zpt75GluChi2->DrawNormalized("histo,same");
969	>	legGluChi2->Draw("same");
970	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingZPt_Of_Three_Different_Gluino_Chi2_MassDifferences");
971	>
972	>	MET25GluChi2->SetMaximum(MET25GluChi2->GetMaximum()*15);
973	>	MET25GluChi2->DrawNormalized("histo");
974	>	MET50GluChi2->DrawNormalized("histo,same");
975	>	MET75GluChi2->DrawNormalized("histo,same");
976	>	legGluChi2->Draw("same");
977	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingMET_Of_Three_Different_Gluino_Chi2_MassDifferences");
978	>
979	>	JZB25GluChi2->SetMaximum(JZB25GluChi2->GetMaximum()*10);
980	>	JZB25GluChi2->DrawNormalized("histo");
981	>	JZB50GluChi2->DrawNormalized("histo,same");
982	>	JZB75GluChi2->DrawNormalized("histo,same");
983	>	//  legGluChi2->Draw("same");
984	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingJZB_Of_Three_Different_Gluino_Chi2_MassDifferences");
985	>
986	>	//----------------------
987	>
988	>
989	>	sumJetpt25Chi2LSP->SetMaximum(sumJetpt25Chi2LSP->GetMaximum()*15);
990	>	sumJetpt25Chi2LSP->DrawNormalized("histo");
991	>	sumJetpt50Chi2LSP->DrawNormalized("histo,same");
992	>	sumJetpt75Chi2LSP->DrawNormalized("histo,same");
993	>	legChi2LSP->Draw("same");
994	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingSumJetPt_Of_Three_Different_Chi1_Chi2_MassDifferences");
995	>
996	>	nJets25Chi2LSP->DrawNormalized("histo");
997	>	nJets50Chi2LSP->DrawNormalized("histo,same");
998	>	nJets75Chi2LSP->DrawNormalized("histo,same");
999	>	legChi2LSP->Draw("same");
1000	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingNJets_Of_Three_Different_Chi1_Chi2_MassDifferences");
1001	>
1002	>	Zpt25Chi2LSP->SetMaximum(Zpt25Chi2LSP->GetMaximum()*30);
1003	>	Zpt25Chi2LSP->DrawNormalized("histo");
1004	>	Zpt50Chi2LSP->DrawNormalized("histo,same");
1005	>	Zpt75Chi2LSP->DrawNormalized("histo,same");
1006	>	legChi2LSP->Draw("same");
1007	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingZPt_Of_Three_Different_Chi1_Chi2_MassDifferences");
1008	>
1009	>	MET25Chi2LSP->SetMaximum(MET25Chi2LSP->GetMaximum()*30);
1010	>	MET25Chi2LSP->DrawNormalized("histo");
1011	>	MET50Chi2LSP->DrawNormalized("histo,same");
1012	>	MET75Chi2LSP->DrawNormalized("histo,same");
1013	>	legChi2LSP->Draw("same");
1014	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingMET_Of_Three_Different_Chi1_Chi2_MassDifferences");
1015	>
1016	>	JZB25Chi2LSP->SetMaximum(JZB25Chi2LSP->GetMaximum()*30);
1017	>	JZB25Chi2LSP->DrawNormalized("histo");
1018	>	JZB50Chi2LSP->DrawNormalized("histo,same");
1019	>	JZB75Chi2LSP->DrawNormalized("histo,same");
1020	>	//  legChi2LSP->Draw("same");
1021	>	CompleteSave(can,"GeneratorStudies/MassDifferences/ComparingJZB_Of_Three_Different_Chi1_Chi2_MassDifferences");
1022		}
972	–
973	–

Diff Legend

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