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]); }