From ec7f2c5f212f7c134296f0fc09619f837a757d76 Mon Sep 17 00:00:00 2001 From: morsch Date: Fri, 24 Oct 2003 12:54:46 +0000 Subject: [PATCH] Redefinition of ptMin and ptMax corrrected. Can be run in compiled mode now (A. De Falco). --- FASTSIM/uncorrBg.C | 33 +++++++++++++++++++++++++-------- 1 file changed, 25 insertions(+), 8 deletions(-) diff --git a/FASTSIM/uncorrBg.C b/FASTSIM/uncorrBg.C index 5f87e9462d1..7d1f3df5be6 100644 --- a/FASTSIM/uncorrBg.C +++ b/FASTSIM/uncorrBg.C @@ -1,3 +1,20 @@ +#if !defined(__CINT__) || defined(__MAKECINT__) +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include "AliFastGlauber.h" +#include "AliFastMuonTrackingRes.h" +#include "AliFastMuonTrackingAcc.h" +#include "AliFastMuonTrackingEff.h" +#include "AliFastMuonTriggerEff.h" +#endif + void uncorrBg(Int_t nev = 1000000, Double_t bmin = 0., Double_t bmax = 5.) { // @@ -17,17 +34,17 @@ void uncorrBg(Int_t nev = 1000000, Double_t bmin = 0., Double_t bmax = 5.) printf("Impact parameter range: %10.3f - %10.3f fm \n", bmin, bmax); printf("Luminosity: %10.3e cm^-2 s^-1 \n", lumi); printf("Rate: %10.3f Hz\n", rate); - printf("Fraction of hard cross-section: %10.3f %\n", fhard * 100.); - printf("Fraction of geom. cross-section: %10.3f %\n", fgeo * 100.); - printf("Events in 10^6 s: %10.3e %\n", events); + printf("Fraction of hard cross-section: %10.3f \n", fhard * 100.); + printf("Fraction of geom. cross-section: %10.3f \n", fgeo * 100.); + printf("Events in 10^6 s: %10.3e \n", events); printf("-------------------------------------------------------------------\n"); // // - Float_t ptMinCut = 3.; + Float_t ptMinCut = 3.; // GeV Float_t etamin = 2.543; Float_t etar = 1.457; - Float_t ptUp = 20.; // GeV + Float_t ptUp = 20.; // GeV Float_t dpt = 0.01; // GeV // // For b = 0 @@ -106,7 +123,7 @@ void uncorrBg(Int_t nev = 1000000, Double_t bmin = 0., Double_t bmax = 5.) // // pi/K -> mu // - f = new TFile("$(ALICE_ROOT)/FASTSIM/data/pikmu.root"); + TFile* f = new TFile("$(ALICE_ROOT)/FASTSIM/data/pikmu.root"); TH2F* etaptPiK = (TH2F*) f->Get("etaptH"); TAxis* etaAxis = etaptPiK->GetXaxis(); TAxis* ptAxis = etaptPiK->GetYaxis(); @@ -239,8 +256,8 @@ void uncorrBg(Int_t nev = 1000000, Double_t bmin = 0., Double_t bmax = 5.) Float_t ptMax = pT1; Float_t ptMin = pT2; if (pT2 > pT1) { - Float_t ptMax = pT2; - Float_t ptMin = pT1; + ptMax = pT2; + ptMin = pT1; } if (ptMin > ptMinCut && p1 > 4. && p2 > 4.) { -- 2.43.0