ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/GeneratorLevelStudyModule.C
(Generate patch)

Comparing UserCode/cbrown/Development/Plotting/Modules/GeneratorLevelStudyModule.C (file contents):
Revision 1.1 by buchmann, Mon Feb 20 17:06:18 2012 UTC vs.
Revision 1.2 by buchmann, Thu Feb 23 15:19:55 2012 UTC

# Line 48 | Line 48 | namespace GenLevelStudy {
48          void fill_pstar_histos(string scandir, float x, TH1F* pstarhistos[], int npstarhistos, TLegend *leg, bool addlegentry);
49          void DrawOnly100to150inJetBand();
50          void DrawOnly100to150inJetBand(float massglu,float masslsp);
51 +        void ImpactOfGluinoChi2MassDifference();
52  
53 <        string LM4sample="/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_PrimeTime__SAMPLES__Gen_V11/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root";
54 <        string SMSsample="/scratch/buchmann/ntuples/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.root";
53 > //      string LM4sample="/scratch/buchmann/ntuples/GeneratorInformationInJZB___JZBplusSamples_PrimeTime__SAMPLES__Gen_V11/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2__withIndex.root";
54 >        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";
55 >        string SMSsample="/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";
56          string histoname(int i);
57          Color_t histocolor(int i);
58          Color_t diversehistocolor(int i);
# Line 434 | Line 436 | void GenLevelStudy::ZDecayIllustration()
436  
437  
438   void GenLevelStudy::pStarIllustration(float x) {
439 <  TH2F *pstartmap = new TH2F("limitmap", "",(int)((mgluend-mglustart)/mglustep+1),mglustart-0.5*mglustep,mgluend+0.5*mglustep,(int)((mLSPend-mLSPstart)/mLSPstep+1),mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
439 >  TH2F *pstarlspmap = new TH2F("pstarlspmap", "",(int)((mgluend-mglustart)/mglustep+1),mglustart-0.5*mglustep,mgluend+0.5*mglustep,(int)((mLSPend-mLSPstart)/mLSPstep+1),mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
440 >  TH2F *pstarchimap = new TH2F("pstarchimap", "",(int)((mgluend-mglustart)/mglustep+1),mglustart-0.5*mglustep,mgluend+0.5*mglustep,(int)((mLSPend-mLSPstart)/mLSPstep+1),mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
441 >  TH2F *pstardiffmap = new TH2F("pstardiffmap", "",(int)((mgluend-mglustart)/mglustep+1),mglustart-0.5*mglustep,mgluend+0.5*mglustep,(int)((mLSPend-mLSPstart)/mLSPstep+1),mLSPstart-0.5*mLSPstep,mLSPend+0.5*mLSPstep);
442  
443    for(int mglu=(int)mglustart;mglu<=mgluend;mglu+=(int)mglustep) {
444      for (int mlsp=(int)mLSPstart;mlsp<=mLSPend&&mlsp<=mglu;mlsp+=(int)mLSPstep)
445      {
446 <        float pstar = pSTAR(mglu,mlsp,mlsp+x*(mglu-mlsp));
447 <        pstartmap->Fill(mglu,mlsp,pstar);
446 >        float pstarlsp = pSTAR(mglu,mlsp,mlsp+x*(mglu-mlsp));
447 >        pstarlspmap->Fill(mglu,mlsp,pstarlsp);
448 >        float pstarchi = generalizedpSTAR(mglu,mlsp+x*(mglu-mlsp),1.0);//the second daughter is a jet; we're using 1 GeV here.
449 >        pstarchimap->Fill(mglu,mlsp,pstarchi);
450 >        pstardiffmap->Fill(mglu,mlsp,pstarchi-pstarlsp);
451      }
452    }
453  
454    stringstream pstartitle;
455    pstartitle << "pstar_histo_x_equal_0p" << 100*x << "____NOT_efficiencymap";
456 <  pstartmap->SetTitle(pstartitle.str().c_str());
457 <  pstartmap->SetName(pstartitle.str().c_str());
458 <  pstartmap->GetXaxis()->SetTitle("m_{#tilde{g}} (GeV/c^{2})");
459 <  pstartmap->GetXaxis()->CenterTitle();
460 <  pstartmap->GetYaxis()->SetTitle("m_{LSP} (GeV/c^{2})");
461 <  pstartmap->GetYaxis()->CenterTitle();
462 <  pstartmap->GetZaxis()->SetTitle("p^{*} (GeV/c)");
463 <  pstartmap->GetZaxis()->CenterTitle();
456 >  pstarlspmap->SetTitle(pstartitle.str().c_str());
457 >  pstarlspmap->SetName(pstartitle.str().c_str());
458 >  pstarlspmap->GetXaxis()->SetTitle("m_{#tilde{g}} (GeV/c^{2})");
459 >  pstarlspmap->GetXaxis()->CenterTitle();
460 >  pstarlspmap->GetYaxis()->SetTitle("m_{LSP} (GeV/c^{2})");
461 >  pstarlspmap->GetYaxis()->CenterTitle();
462 >  pstarlspmap->GetZaxis()->SetTitle("p^{*}(LSP) (GeV/c)");
463 >  pstarlspmap->GetZaxis()->CenterTitle();
464  
465    float corr=91.2/x;
466    TLine *l = new TLine(corr,0,1200,1200-corr);
# Line 465 | Line 472 | void GenLevelStudy::pStarIllustration(fl
472  
473    TCanvas *can = new TCanvas("can","can");
474    can->SetRightMargin(0.16);
475 <  pstartmap->Draw("COLZ");
475 >  pstarlspmap->Draw("COLZ");
476    l->Draw("same");
477    l2->Draw("same");
478    stringstream saveas;
479 <  saveas << "GeneratorStudies/PStar_x_0p" << x*100;
479 >  saveas << "GeneratorStudies/PStarLSP_x_0p" << x*100;
480    CompleteSave(can,saveas.str());  
481 +  
482 +  stringstream pstartitle2;
483 +  pstartitle2 << "pstar_histo_x_equal_0p" << 100*x << "____NOT_efficiencymap";
484 +  pstarchimap->SetTitle(pstartitle2.str().c_str());
485 +  pstarchimap->SetName(pstartitle2.str().c_str());
486 +  pstarchimap->GetXaxis()->SetTitle("m_{#tilde{g}} (GeV/c^{2})");
487 +  pstarchimap->GetXaxis()->CenterTitle();
488 +  pstarchimap->GetYaxis()->SetTitle("m_{LSP} (GeV/c^{2})");
489 +  pstarchimap->GetYaxis()->CenterTitle();
490 +  pstarchimap->GetZaxis()->SetTitle("p^{*}(#chi_{2}^{0}) (GeV/c)");
491 +  pstarchimap->GetZaxis()->CenterTitle();
492 +  pstarchimap->Draw("COLZ");
493 +  l->Draw("same");
494 +  l2->Draw("same");
495 +  stringstream saveas2;
496 +  saveas2 << "GeneratorStudies/PStarCHI_x_0p" << x*100;
497 +  CompleteSave(can,saveas2.str());  
498 +  
499 +  stringstream pstartitle3;
500 +  pstartitle3 << "pstar_histo_x_equal_0p" << 100*x << "____NOT_efficiencymap";
501 +  pstardiffmap->SetTitle(pstartitle3.str().c_str());
502 +  pstardiffmap->SetName(pstartitle3.str().c_str());
503 +  pstardiffmap->GetXaxis()->SetTitle("m_{#tilde{g}} (GeV/c^{2})");
504 +  pstardiffmap->GetXaxis()->CenterTitle();
505 +  pstardiffmap->GetYaxis()->SetTitle("m_{LSP} (GeV/c^{2})");
506 +  pstardiffmap->GetYaxis()->CenterTitle();
507 +  pstardiffmap->GetZaxis()->SetTitle("p^{*}(#chi_{2}^{0}) -p^{*}(LSP) (GeV/c)");
508 +  pstardiffmap->GetZaxis()->CenterTitle();
509 + //  pstardiffmap->GetZaxis()->SetRangeUser(0,250);
510 +  pstardiffmap->Draw("COLZ");
511 +  l->Draw("same");
512 +  l2->Draw("same");
513 +  stringstream saveas3;
514 +  saveas3 << "GeneratorStudies/PStar_CHI_minus_LSP__x_0p" << x*100;
515 +  CompleteSave(can,saveas3.str());  
516 +  
517 +  
518    TFile *g = new TFile("pstarshapes.root","UPDATE");
519 <  pstartmap->Write();
519 >  pstarlspmap->Write();
520 >  pstarchimap->Write();
521 >  pstardiffmap->Write();
522    g->Close();
523    delete can;
524 <  delete pstartmap;
524 >  delete pstarlspmap;
525 >  delete pstarchimap;
526 >  delete pstardiffmap;
527   }
528  
529   void GenLevelStudy::DeltaLSPmomentumScenario(float massglu, float masslsp) {
# Line 703 | Line 751 | void GenLevelStudy::DrawJetBand() {
751    GenLevelStudy::DrawJetBand(625,425);
752   }
753  
754 + void GenLevelStudy::ImpactOfGluinoChi2MassDifference() {
755 +  TCanvas *can = new TCanvas("can","can");
756 +  can->SetLogy(1);
757 +  TChain smsp025events("events");
758 +  TChain smsp050events("events");
759 +  TChain smsp075events("events");
760 +  
761 +  for(int i=0;i<15;i++) {
762 +    for(int j=0;j<15;j++) {
763 +       stringstream filenamep25;
764 +       stringstream filenamep50;
765 +       stringstream filenamep75;
766 +       filenamep25 << "/shome/buchmann/ntuples/SMS_T5zzh_newslots/SMS_clean_splitup_" << i << "_" << j << ".root";
767 +       filenamep50 << "/shome/buchmann/ntuples/SMS_T5zz/SMS_clean_splitup_" << i << "_" << j << ".root";
768 +       filenamep75 << "/shome/buchmann/ntuples/SMS_T5zzl/SMS_clean_splitup_" << i << "_" << j << ".root";
769 +       smsp025events.Add(filenamep25.str().c_str());
770 +       smsp050events.Add(filenamep50.str().c_str());
771 +       smsp075events.Add(filenamep75.str().c_str());
772 +    }
773 +   }
774 +  
775 +  stringstream cutp25m2;
776 +  cutp25m2 << "abs(MassGlu-700)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
777 +  stringstream cutp50m2;
778 +  cutp50m2 << "abs(MassGlu-400)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
779 +  stringstream cutp75m2;
780 +  cutp75m2 << "abs(MassGlu-300)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
781 +  
782 +  stringstream cutp25m1;
783 +  cutp25m1 << "abs(MassGlu-700)<5&&abs(MassLSP-500)&&pfJetGoodNum>0";
784 +  stringstream cutp50m1;
785 +  cutp50m1 << "abs(MassGlu-700)<5&&abs(MassLSP-400)&&pfJetGoodNum>0";
786 +  stringstream cutp75m1;
787 +  cutp75m1 << "abs(MassGlu-700)<5&&abs(MassLSP-100)&&pfJetGoodNum>0";
788 +  
789 +  TH1F *nJets25m2 = new TH1F("nJets25m2","nJets25m2",20,-0.5,19.5);
790 +  TH1F *nJets50m2 = new TH1F("nJets50m2","nJets50m2",20,-0.5,19.5);
791 +  TH1F *nJets75m2 = new TH1F("nJets75m2","nJets75m2",20,-0.5,19.5);
792 +  
793 +  TH1F *sumJetpt25m2 = new TH1F("sumJetpt25m2","sumJetpt25m2",100,0,300);
794 +  TH1F *sumJetpt50m2 = new TH1F("sumJetpt50m2","sumJetpt50m2",100,0,300);
795 +  TH1F *sumJetpt75m2 = new TH1F("sumJetpt75m2","sumJetpt75m2",100,0,300);
796 +  
797 +  TH1F *Zpt25m2 = new TH1F("Zpt25m2","Zpt25m2",100,0,200);
798 +  TH1F *Zpt50m2 = new TH1F("Zpt50m2","Zpt50m2",100,0,200);
799 +  TH1F *Zpt75m2 = new TH1F("Zpt75m2","Zpt75m2",100,0,200);
800 +  
801 +  TH1F *nJets25m1 = new TH1F("nJets25m1","nJets25m1",20,-0.5,19.5);
802 +  TH1F *nJets50m1 = new TH1F("nJets50m1","nJets50m1",20,-0.5,19.5);
803 +  TH1F *nJets75m1 = new TH1F("nJets75m1","nJets75m1",20,-0.5,19.5);
804 +  
805 +  TH1F *sumJetpt25m1 = new TH1F("sumJetpt25m1","sumJetpt25m1",100,0,300);
806 +  TH1F *sumJetpt50m1 = new TH1F("sumJetpt50m1","sumJetpt50m1",100,0,300);
807 +  TH1F *sumJetpt75m1 = new TH1F("sumJetpt75m1","sumJetpt75m1",100,0,300);
808 +  
809 +  TH1F *Zpt25m1 = new TH1F("Zpt25m1","Zpt25m1",100,0,200);
810 +  TH1F *Zpt50m1 = new TH1F("Zpt50m1","Zpt50m1",100,0,200);
811 +  TH1F *Zpt75m1 = new TH1F("Zpt75m1","Zpt75m1",100,0,200);
812 +  
813 +  smsp025events.Draw("pfJetGoodNum>>nJets25m2",("abs(mll-91.2)&&id1==id2&&"+cutp25m2.str()).c_str());
814 +  smsp050events.Draw("pfJetGoodNum>>nJets50m2",("abs(mll-91.2)&&id1==id2&&"+cutp50m2.str()).c_str());
815 +  smsp075events.Draw("pfJetGoodNum>>nJets75m2",("abs(mll-91.2)&&id1==id2&&"+cutp75m2.str()).c_str());
816 +  
817 +  smsp025events.Draw("sumJetPt[1]>>sumJetpt25m2",("abs(mll-91.2)&&id1==id2&&"+cutp25m2.str()).c_str());
818 +  smsp050events.Draw("sumJetPt[1]>>sumJetpt50m2",("abs(mll-91.2)&&id1==id2&&"+cutp50m2.str()).c_str());
819 +  smsp075events.Draw("sumJetPt[1]>>sumJetpt75m2",("abs(mll-91.2)&&id1==id2&&"+cutp75m2.str()).c_str());
820 +  
821 +  smsp025events.Draw("pt>>Zpt25m2",("abs(mll-91.2)&&id1==id2&&"+cutp25m2.str()).c_str());
822 +  smsp050events.Draw("pt>>Zpt50m2",("abs(mll-91.2)&&id1==id2&&"+cutp50m2.str()).c_str());
823 +  smsp075events.Draw("pt>>Zpt75m2",("abs(mll-91.2)&&id1==id2&&"+cutp75m2.str()).c_str());
824 +  
825 +  smsp025events.Draw("pfJetGoodNum>>nJets25m1",("abs(mll-91.2)&&id1==id2&&"+cutp25m1.str()).c_str());
826 +  smsp050events.Draw("pfJetGoodNum>>nJets50m1",("abs(mll-91.2)&&id1==id2&&"+cutp50m1.str()).c_str());
827 +  smsp075events.Draw("pfJetGoodNum>>nJets75m1",("abs(mll-91.2)&&id1==id2&&"+cutp75m1.str()).c_str());
828 +  
829 +  smsp025events.Draw("sumJetPt[1]>>sumJetpt25m1",("abs(mll-91.2)&&id1==id2&&"+cutp25m1.str()).c_str());
830 +  smsp050events.Draw("sumJetPt[1]>>sumJetpt50m1",("abs(mll-91.2)&&id1==id2&&"+cutp50m1.str()).c_str());
831 +  smsp075events.Draw("sumJetPt[1]>>sumJetpt75m1",("abs(mll-91.2)&&id1==id2&&"+cutp75m1.str()).c_str());
832 +  
833 +  smsp025events.Draw("pt>>Zpt25m1",("abs(mll-91.2)&&id1==id2&&"+cutp25m1.str()).c_str());
834 +  smsp050events.Draw("pt>>Zpt50m1",("abs(mll-91.2)&&id1==id2&&"+cutp50m1.str()).c_str());
835 +  smsp075events.Draw("pt>>Zpt75m1",("abs(mll-91.2)&&id1==id2&&"+cutp75m1.str()).c_str());
836 +  
837 +  sumJetpt25m1->GetXaxis()->SetTitle("Sum of Jet Pt (GeV)");
838 +  sumJetpt25m1->GetYaxis()->SetTitle("events");
839 +  sumJetpt25m1->GetXaxis()->CenterTitle();
840 +  sumJetpt25m1->GetYaxis()->CenterTitle();
841 +  
842 +  sumJetpt25m1->SetLineColor(diversehistocolor(0));
843 +  sumJetpt50m1->SetLineColor(diversehistocolor(1));
844 +  sumJetpt75m1->SetLineColor(diversehistocolor(2));
845 +  
846 +  TLegend *leg = make_legend("",0.6,0.6,false);
847 +  leg->AddEntry(sumJetpt25m1,"#Delta m_{2} = 50 GeV","l");
848 +  leg->AddEntry(sumJetpt50m1,"#Delta m_{2} = 150 GeV","l");
849 +  leg->AddEntry(sumJetpt75m1,"#Delta m_{2} = 450 GeV","l");
850 +  
851 +  sumJetpt25m1->Draw("histo");
852 +  sumJetpt50m1->Draw("histo,same");
853 +  sumJetpt75m1->Draw("histo,same");
854 +  leg->Draw("same");
855 +  CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM2/ComparingSumJetPt_Of_Three_Different_Gluino_Chi2_MassDifferences");
856 +  
857 +  nJets25m1->SetLineColor(diversehistocolor(0));
858 +  nJets50m1->SetLineColor(diversehistocolor(1));
859 +  nJets75m1->SetLineColor(diversehistocolor(2));  
860 +  nJets25m1->GetXaxis()->SetTitle("N(Jets)");
861 +  nJets25m1->GetXaxis()->CenterTitle();
862 +  nJets25m1->GetYaxis()->SetTitle("events");
863 +  nJets25m1->GetYaxis()->CenterTitle();
864 +  nJets25m1->Draw("histo");
865 +  nJets50m1->Draw("histo,same");
866 +  nJets75m1->Draw("histo,same");
867 +  leg->Draw("same");
868 +  CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM2/ComparingNJets_Of_Three_Different_Gluino_Chi2_MassDifferences");
869 +  
870 +  Zpt25m1->GetXaxis()->SetTitle("Z p_{T} (GeV)");
871 +  Zpt25m1->GetXaxis()->CenterTitle();
872 +  Zpt25m1->GetYaxis()->SetTitle("events");
873 +  Zpt25m1->GetYaxis()->CenterTitle();
874 +  Zpt25m1->SetLineColor(diversehistocolor(0));
875 +  Zpt50m1->SetLineColor(diversehistocolor(1));
876 +  Zpt75m1->SetLineColor(diversehistocolor(2));
877 +  Zpt25m1->Draw("histo");
878 +  Zpt50m1->Draw("histo,same");
879 +  Zpt75m1->Draw("histo,same");
880 +  leg->Draw("same");
881 +  CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM2/ComparingZPt_Of_Three_Different_Gluino_Chi2_MassDifferences");
882 +  
883 + //------------
884 +  
885 +  sumJetpt25m2->GetXaxis()->SetTitle("Sum of Jet Pt (GeV)");
886 +  sumJetpt25m2->GetYaxis()->SetTitle("events");
887 +  sumJetpt25m2->GetXaxis()->CenterTitle();
888 +  sumJetpt25m2->GetYaxis()->CenterTitle();
889 +  
890 +  sumJetpt25m2->SetLineColor(diversehistocolor(0));
891 +  sumJetpt50m2->SetLineColor(diversehistocolor(1));
892 +  sumJetpt75m2->SetLineColor(diversehistocolor(2));
893 +  
894 +  TLegend *leg2 = make_legend("",0.6,0.6,false);
895 +  leg2->AddEntry(sumJetpt75m2,"#Delta m_{1} = 50 GeV","l");
896 +  leg2->AddEntry(sumJetpt50m2,"#Delta m_{1} = 150 GeV","l");
897 +  leg2->AddEntry(sumJetpt25m2,"#Delta m_{1} = 450 GeV","l");
898 +  
899 +  sumJetpt25m2->Draw("histo");
900 +  sumJetpt50m2->Draw("histo,same");
901 +  sumJetpt75m2->Draw("histo,same");
902 +  leg2->Draw("same");
903 +  CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM1/ComparingSumJetPt_Of_Three_Different_Chi1_Chi2_MassDifferences");
904 +  
905 +  nJets25m2->SetLineColor(diversehistocolor(0));
906 +  nJets50m2->SetLineColor(diversehistocolor(1));
907 +  nJets75m2->SetLineColor(diversehistocolor(2));  
908 +  nJets25m2->GetXaxis()->SetTitle("N(Jets)");
909 +  nJets25m2->GetXaxis()->CenterTitle();
910 +  nJets25m2->GetYaxis()->SetTitle("events");
911 +  nJets25m2->GetYaxis()->CenterTitle();
912 +
913 +  nJets25m2->Draw("histo");
914 +  nJets50m2->Draw("histo,same");
915 +  nJets75m2->Draw("histo,same");
916 +  leg2->Draw("same");
917 +  CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM1/ComparingNJets_Of_Three_Different_Chi1_Chi2_MassDifferences");
918 +  
919 +  Zpt25m2->GetXaxis()->SetTitle("Z p_{T} (GeV)");
920 +  Zpt25m2->GetXaxis()->CenterTitle();
921 +  Zpt25m2->GetYaxis()->SetTitle("events");
922 +  Zpt25m2->GetYaxis()->CenterTitle();
923 +  
924 +  Zpt25m2->SetLineColor(diversehistocolor(0));
925 +  Zpt50m2->SetLineColor(diversehistocolor(1));
926 +  Zpt75m2->SetLineColor(diversehistocolor(2));  
927 +  Zpt25m2->Draw("histo");
928 +  Zpt50m2->Draw("histo,same");
929 +  Zpt75m2->Draw("histo,same");
930 +  leg2->Draw("same");
931 +  CompleteSave(can,"GeneratorStudies/MassDifferences/VaryingDeltaM1/ComparingZPt_Of_Three_Different_Chi1_Chi2_MassDifferences");
932 +  
933 + //--------------
934 +  delete can;
935 +  
936 +  delete nJets25m1;
937 +  delete nJets50m1;
938 +  delete nJets75m1;
939 +  
940 +  delete nJets25m2;
941 +  delete nJets50m2;
942 +  delete nJets75m2;
943 +  
944 +  delete sumJetpt25m1;
945 +  delete sumJetpt50m1;
946 +  delete sumJetpt75m1;
947 +  
948 +  delete sumJetpt25m2;
949 +  delete sumJetpt50m2;
950 +  delete sumJetpt75m2;
951 + }
952 +
953  

Diff Legend

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