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

Comparing UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C (file contents):
Revision 1.83 by buchmann, Mon Feb 11 13:32:18 2013 UTC vs.
Revision 1.84 by buchmann, Mon Feb 18 18:55:58 2013 UTC

# Line 1606 | Line 1606 | void ProduceMCSystematicPlots() {
1606    cout << "Kicking cutWeight  " << (const char*) cutWeight << endl;
1607    cout << "Keeping OSSF cut " << (const char*)cutOSSF << endl;
1608    cutWeight="1.0";
1609 –  TCut bkpessential=essentialcut;
1610 –  essentialcut=TCut("genPt1>0");
1609    
1610 <  QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",TCut("abs(genId1)==abs(genId2) && abs(genEta1)<1.4 && abs(genEta2)<1.4 && genPt1>20 && genPt2>20 &&  genMll>20&&genMll<70"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE__NOPURW___GENLEVEL");
1613 <  QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",TCut("abs(genId1)!=abs(genId2) && abs(genEta1)<1.4 && abs(genEta2)<1.4 && genPt1>20 && genPt2>20 &&  genMll>20&&genMll<70"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE__NOPURW___GENLEVEL");
1614 <  QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",TCut("abs(genId1)==abs(genId2) && abs(genEta1)<1.4 && abs(genEta2)<1.4 && genPt1>20 && genPt2>20 &&  genMll>120"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE_HighMass__NOPURW___GENLEVEL");
1615 <  QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",TCut("abs(genId1)!=abs(genId2) && abs(genEta1)<1.4 && abs(genEta2)<1.4 && genPt1>20 && genPt2>20 &&  genMll>120"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE_HighMass__NOPURW___GENLEVEL");
1616 <  essentialcut=bkpessential;
1617 <  assert(0);
1618 <  
1619 < /*  QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE__NOPURW");
1610 >  QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE__NOPURW");
1611    QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE__NOPURW");  
1612    QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE_HighMass__NOPURW");
1613    QuickProduceMCSystematicPlot("genMET",20,0,200,"genMET [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE_HighMass__NOPURW");  
# Line 1625 | Line 1616 | void ProduceMCSystematicPlots() {
1616    QuickProduceMCSystematicPlot("met[4]",20,0,200,"PFMET [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_OppositeFlavor_INCLUSIVE__NOPURW");
1617    QuickProduceMCSystematicPlot("met[4]",20,0,200,"PFMET [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/MET/ETH_SameFlavor_INCLUSIVE_HighMass__NOPURW");
1618    QuickProduceMCSystematicPlot("met[4]",20,0,200,"PFMET [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/MET/ETH_OppositeFlavor_INCLUSIVE_HighMass__NOPURW");
1619 <  */
1619 >  
1620    
1621    // --------------------------------------    ***** AACHEN *****      --------------------------------------    
1622    
# Line 1765 | Line 1756 | void do_kinematic_plots(string mcjzb, st
1756    bool dolog=true;
1757    bool nolog=false;
1758  
1759 <  bool doOFSF = true;
1759 >  bool doOFSF = false;
1760    bool doKin  = true;
1761    bool doDataComp = false;
1762    bool MakeTwoThreeJetComparison = false;
# Line 1784 | Line 1775 | void do_kinematic_plots(string mcjzb, st
1775    //TauQuestion();
1776  
1777    //ProduceJanPlots();
1778 <  ProduceMCSystematicPlots();
1788 <  assert(0);
1778 >  //ProduceMCSystematicPlots();
1779   //  make_plain_kin_plot("pt",Cut2Str(cutOSSF&&TCut("mll>20&&pfJetGoodNumBtag30>=2")),40,0,200,nolog,"Z p_{T}","PlotForKostas",doPF,true);
1780    
1781    if ( doOFSF ) {
# Line 1911 | Line 1901 | void do_kinematic_plots(string mcjzb, st
1901  
1902    if ( doKin ) {
1903      string mllCut("");
1904 <    /*
1904 >    
1905      
1906    
1907      make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_0j",doPF,true);
# Line 1963 | Line 1953 | void do_kinematic_plots(string mcjzb, st
1953      make_kin_plot("MetFactor","",20,0,2,nolog,"MetFactor","MetFactor",doPF,true);
1954      
1955      make_ttbar_comparison("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_30to70_2j__ALLMCSAMPLES");
1956 <     */
1956 >    /*
1957      DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),dolog,"PFMET","Debunking/MET_OF_30to70_0j");
1958      DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),dolog,"PFMET","Debunking/MET_OF_30to70_1j");
1959      DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_30to70_2j");
# Line 2058 | Line 2048 | void do_kinematic_plots(string mcjzb, st
2048      
2049      //   make_kin_plot("mll","met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV] (MET>100GeV)","mll_met100_ll",doPF,true);
2050      //make_kin_plot("mll","met[4]>150&&id1==0",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ee} [GeV] (MET>150GeV)","mll_met150_ee",doPF,true);
2051 <    //make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true);
2051 >    //make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true);*/
2052    }
2053 <    
2053 >    /*
2054     make_special_obs_pred_mll_plot(datajzb,mcjzb,0);
2055     make_special_obs_pred_mll_plot(datajzb,mcjzb,50);
2056     make_special_obs_pred_mll_plot(datajzb,mcjzb,80);
# Line 2071 | Line 2061 | void do_kinematic_plots(string mcjzb, st
2061  
2062     make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf");
2063     make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof");
2064 <
2064 > */
2065    switch_overunderflow(false);
2066   }
2067  
# Line 3713 | Line 3703 | void find_correction_factors(string &jzb
3703    
3704   }
3705  
3706 < void pick_up_events(string cut) {
3707 <  dout << "Picking up events with cut " << cut << endl;
3718 <  allsamples.PickUpEvents(cut);
3706 > void pick_up_events(string cut, string filename) {
3707 >  allsamples.PickUpEvents(cut,filename);
3708   }
3709  
3710   void save_template(string mcjzb, string datajzb,vector<float> jzb_cuts,float MCPeakError,float DataPeakError, vector<float> jzb_shape_limit_bins) {

Diff Legend

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