| 1 |
|
#include "TheLimits.h" |
| 2 |
|
#include "SusyScan.h" |
| 3 |
|
#include "GeneratorMasses.h" |
| 4 |
– |
#include "IsoMassLine.h" |
| 5 |
– |
#include "TGraph.h" |
| 6 |
– |
|
| 7 |
– |
#include "TH2.h" |
| 8 |
– |
#include "TH2F.h" |
| 9 |
– |
|
| 10 |
– |
void TheLimits::plot(TH2*h, double(*x)(const SusyScan*), double(*y)(const SusyScan*), double(*f)(const SusyScan*) ) |
| 11 |
– |
{ |
| 12 |
– |
for (std::vector<SusyScan*>::const_iterator it=_scan.begin();it!=_scan.end();++it ) |
| 13 |
– |
//h->Fill( x(*it), y(*it), f(*it) ); |
| 14 |
– |
h->SetBinContent( h->GetXaxis()->FindBin(x(*it)), |
| 15 |
– |
h->GetYaxis()->FindBin(y(*it)), f(*it) ); |
| 16 |
– |
} |
| 4 |
|
|
| 5 |
|
void TheLimits::match(const std::vector<GeneratorMasses> genm) |
| 6 |
|
{ |
| 7 |
|
for (std::vector<SusyScan*>::iterator it=_scan.begin();it!=_scan.end();++it){ |
| 8 |
|
bool match = false; |
| 9 |
|
for (std::vector<GeneratorMasses>::const_iterator gt=genm.begin();gt!=genm.end();++gt){ |
| 10 |
< |
if ((*it)->Mzero==gt->MZERO && |
| 11 |
< |
(*it)->Mhalf==gt->MHALF && |
| 12 |
< |
(*it)->Azero==gt->AZERO && |
| 13 |
< |
(*it)->TanBeta==gt->TANB && |
| 14 |
< |
(*it)->Mu==gt->SGNMU ) { |
| 10 |
> |
if ((*it)->Mzero==gt->Mzero && |
| 11 |
> |
(*it)->Mhalf==gt->Mhalf && |
| 12 |
> |
(*it)->Azero==gt->Azero && |
| 13 |
> |
(*it)->TanBeta==gt->TanBeta && |
| 14 |
> |
(*it)->Mu==gt->Mu ) { |
| 15 |
|
(*it)->M1 = gt->M1; |
| 16 |
|
(*it)->M2 = gt->M2; |
| 17 |
|
(*it)->M3 = gt->M3; |
