Fix leaks in AliIsolationCut and AlianaOmegaToPi0Gamma

[u/mrichter/AliRoot.git] / FASTSIM / fastGenQuarkonia.C

diff --git a/FASTSIM/fastGenQuarkonia.C b/FASTSIM/fastGenQuarkonia.C
index a1c365cf94f6d886a7f1bbd3dc4e811c8fff2c44..ee6512c77978dd08914eb0a41739399fd7bf8146 100644 (file)
--- a/FASTSIM/fastGenQuarkonia.C
+++ b/FASTSIM/fastGenQuarkonia.C
@@ -110,7 +110,7 @@ void fastGenQuarkonia(Int_t nev = 1, char* filename = "galice.root")
 // Smear
            // the mu+
            Float_t thetas1, phis1, ps1, thetas2, phis2, ps2, pts1, pts2, etas1, etas2;
-           Float_t trigeffpL, trigeffpH, trigeffnL, trigeffnH;
+           Float_t trigeffpL, trigeffpH, trigeffnL, trigeffnH, trigeffnA;
            
            res->SetCharge(1);
            eff->SetCharge(1);
@@ -120,7 +120,7 @@ void fastGenQuarkonia(Int_t nev = 1, char* filename = "galice.root")
            
            Float_t effp     = eff->Evaluate(pt1, thetad1, phid1);
            Float_t accp     = acc->Evaluate(pt1, thetad1, phid1);
-           trigeff->Evaluate(1, pt1, thetad1, phid1, trigeffpL, trigeffpH);
+           trigeff->Evaluate(1, pt1, thetad1, phid1, trigeffpL, trigeffpH, trigeffnA);
            thetas1 *= TMath::Pi()/180.;
            phis1 *= TMath::Pi()/180.;
            Float_t etas;
@@ -138,7 +138,7 @@ void fastGenQuarkonia(Int_t nev = 1, char* filename = "galice.root")
            res->Evaluate(p2, thetad2, phid2, ps2, thetas2, phis2);
            Float_t effn     = eff->Evaluate(pt2, thetad2, phid2);
            Float_t accn     = acc->Evaluate(pt2, thetad2, phid2);
-           trigeff->Evaluate(-1, pt2, thetad2, phid2, trigeffnL, trigeffnH);
+           trigeff->Evaluate(-1, pt2, thetad2, phid2, trigeffnL, trigeffnH, trigeffnA);
            thetas2 *= TMath::Pi()/180.;
            phis2 *= TMath::Pi()/180.;
            if (TMath::Tan(thetas2/2.) > 0) {
@@ -169,6 +169,7 @@ void fastGenQuarkonia(Int_t nev = 1, char* filename = "galice.root")
        rl->TreeE()->Fill();
        header->SetNprimary(stack->GetNprimary());
        header->SetNtrack(stack->GetNtrack());  
+       rl->WriteKinematics("OVERWRITE");
     } // event loop
     TCanvas *c1 = new TCanvas("c1","Canvas 1",400,10,600,700);
     massHU->Draw();
@@ -197,7 +198,7 @@ AliGenerator*  CreateGenerator()
         
     gener->SetMomentumRange(0,999);
     gener->SetPtRange(0,100.);
-    gener->SetPhiRange(-180, 180);
+    gener->SetPhiRange(0., 360.);
     gener->SetYRange(2.5,4);
     gener->SetCutOnChild(1);
     gener->SetChildThetaRange(2.0,9.);