double acceptance2(double accrecoil, double n,
double fbeta=0.27, double frecoil=0.71)
{
TF1* f = (TF1*)(gROOT->GetFunction("acc"));
if(!f)
f = new TF1("acc","acceptance([0],[1],x)",0,1);
f->SetParameters(n,frecoil);
double fcut = f->GetX(accrecoil);
return acceptance(n,fbeta,fcut);
}
void acceptanceR(){
//getFraction();
//acceptance("TTbarMcatnlo");
//acceptance("TTbarTuneZ2");
acceptance("TTbarFullLepMGDecays");
//acceptance("TTbarPowheg");
//acceptance("TTbar_matchingdown");
//acceptance("TTbar_matchingup");
//acceptance("TTbar_scaledown");
//acceptance("TTbar_scaleup");
}
int main(int argc, const char** argv){
acceptance(argv[1],argv[2]);
}
