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); |
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"); |
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"); |