From 5a622da07efaf0663ba78e9c9c325b90daf7d0ee Mon Sep 17 00:00:00 2001 From: nemu Date: Tue, 15 Dec 2009 08:12:03 +0000 Subject: [PATCH] newly added --- src/tests/t0NotEqFirstGoodData/testPoisson.C | 45 ++++++++++++++++++++ 1 file changed, 45 insertions(+) create mode 100644 src/tests/t0NotEqFirstGoodData/testPoisson.C diff --git a/src/tests/t0NotEqFirstGoodData/testPoisson.C b/src/tests/t0NotEqFirstGoodData/testPoisson.C new file mode 100644 index 00000000..d0045bd8 --- /dev/null +++ b/src/tests/t0NotEqFirstGoodData/testPoisson.C @@ -0,0 +1,45 @@ +/*************************************************************************** + + testPoisson.C + + Author: Andreas Suter + e-mail: andreas.suter@psi.ch + + $Id$ + +***************************************************************************/ + +/*************************************************************************** + * Copyright (C) 2009 by Andreas Suter * + * andreas.suter@psi.ch * + * * + * This program is free software; you can redistribute it and/or modify * + * it under the terms of the GNU General Public License as published by * + * the Free Software Foundation; either version 2 of the License, or * + * (at your option) any later version. * + * * + * This program is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * + * GNU General Public License for more details. * + * * + * You should have received a copy of the GNU General Public License * + * along with this program; if not, write to the * + * Free Software Foundation, Inc., * + * 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. * + ***************************************************************************/ + +{ + gSystem->Load("$ROOTSYS/lib/libPAddPoissonNoise"); + + PAddPoissonNoise *apn = new PAddPoissonNoise(); + + const Double_t mean = 10.0; + TH1F *h = new TH1F("h", "h", 31, -0.5, 30.5); + + for (UInt_t i=0; i<1e6; i++) { + h->Fill(apn->AddNoise(mean)); + } + + h->Draw(); +}