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.5 by buchmann, Mon Apr 16 17:22:33 2012 UTC vs.
Revision 1.6 by buchmann, Wed Sep 19 09:16:24 2012 UTC

# Line 483 | Line 483 | void GenLevelStudy::pStarIllustration(fl
483    pstarlspmap->GetZaxis()->SetTitle("p^{*}(LSP) (GeV/c)");
484    pstarlspmap->GetZaxis()->CenterTitle();
485  
486 <  float corr=91.2/x;
486 >  float corr=91/x;
487    TLine *l = new TLine(corr,0,1200,1200-corr);
488    l->SetLineWidth(3);
489    l->SetLineColor(kWhite);
# Line 620 | Line 620 | void GenLevelStudy::fill_pstar_histos(st
620     float stepsize=(pstarhigh-pstarlow)/(float)npstarhistos;
621  
622     stringstream accdrawcommand;
623 <   accdrawcommand << "1+((" << x << "*(MassGlu+MassLSP)-MassLSP)>91.2)";
623 >   accdrawcommand << "1+((" << x << "*(MassGlu+MassLSP)-MassLSP)>91)";
624    
625     TCanvas *ca = new TCanvas("ca","ca");
626     events->Draw(accdrawcommand.str().c_str(),"","COLZ");
# Line 631 | Line 631 | void GenLevelStudy::fill_pstar_histos(st
631     for(int i=0;i<npstarhistos;i++) {
632          float curr=(i+1)*stepsize;
633          stringstream cut;
634 <        cut << "pfJetGoodNum>=3&&TMath::Abs(mll-91.2)<20&&id1==id2&&";
635 <        cut << "(abs(TMath::Sqrt((("<< x*x << ")*(MassGlu+MassLSP)*(MassGlu+MassLSP)-(MassLSP+91.2)*(MassLSP+91.2))*((" << x*x << ")*(MassGlu+MassLSP)*(MassGlu+MassLSP)-(MassLSP-91.2)*(MassLSP-91.2)))/(2*(" << x << ")*(MassGlu+MassLSP))-" << curr-stepsize/2.0 << ")<" << stepsize/2.0 << ")";
636 <        cut << "&&((" << x << "*(MassGlu+MassLSP)-MassLSP)>91.2)"; // cleaning all events out that dont' have the imposed x
634 >        cut << "pfJetGoodNum>=3&&TMath::Abs(mll-91)<10&&id1==id2&&";
635 >        cut << "(abs(TMath::Sqrt((("<< x*x << ")*(MassGlu+MassLSP)*(MassGlu+MassLSP)-(MassLSP+91)*(MassLSP+91))*((" << x*x << ")*(MassGlu+MassLSP)*(MassGlu+MassLSP)-(MassLSP-91)*(MassLSP-91)))/(2*(" << x << ")*(MassGlu+MassLSP))-" << curr-stepsize/2.0 << ")<" << stepsize/2.0 << ")";
636 >        cut << "&&((" << x << "*(MassGlu+MassLSP)-MassLSP)>91)"; // cleaning all events out that dont' have the imposed x
637          stringstream drawcommand;
638          drawcommand << "jzb[1]>>tempie";
639          TH1F *temp = (TH1F*)pstarhistos[i]->Clone("tempie");

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines