ROOT logo
//#include <TClonesArray.h> 

#include <TDatabasePDG.h>
#include <TFile.h>
#include "AliConst.h"
#include "AliGenMUONLMR.h" 
#include "AliMC.h" 
#include "AliRun.h" 
#include "AliLog.h" 
#include "AliGenEventHeader.h"

ClassImp(AliGenMUONLMR)

  AliGenMUONLMR::AliGenMUONLMR () : AliGenMC(), 
				    fNMuMin(2), 
				    fCMSEnergy(kNCMSEnergies),
				    fGenSingleProc(-1),
				    fYCM(0), 
				    fCosTheta(0x0), 
				    fRhoLineShape(0x0),  
				    fHMultMu(0x0), 
				    fHNProc(0x0) { 
    //
    // default constructor 
    //
    for (int i=0; i<fgkNpart; i++) {
      fPDG[i] = 0;
      fScaleMult[i] = 1.;
      fPt[i] = NULL;
      fY[i] = NULL;
      fMult[i] = NULL;
      fParticle[i] = NULL;
    }
    for (int i=0; i<2; i++) {
      fMu[i]     = NULL;
      fDecay[i] = NULL;
    }						 

    for (int i=0; i<3; i++) {
      fDalitz[i] = NULL;
    }						 

  }

//-----------------------------------------------------------

void AliGenMUONLMR::SetCMSEnergy(CMSEnergies energy){
  fCMSEnergy = energy;
  // initialize pt and y distributions according to a fit to 
  // Pythia simulation at sqrt(s) = 7 TeV
  for (Int_t ipart=0; ipart < fgkNpart; ipart++) fScaleMult[ipart] = 1; 
  fScaleMult[kPionLMR] = 0; // set pion multiplicity to zero 
  fScaleMult[kKaonLMR] = 0; // set kaon multiplicity to zero
  const char* fdname[2] = {"fDecPion","fDecKaon"};
  Double_t ctau[2] = {7.8045, 3.712};  
  Int_t pdg[7] = {211, 321, 221, 113, 223, 333, 331}; 
  const char* fptname[7] = {"fPtPion","fPtKaon","fPtEta","fPtRho","fPtOmega","fPtPhi","fPtEtaPrime"};
  const char* fyname[7] = {"fYPion","fYKaon","fYEta","fYRho","fYOmega","fYPhi","fYEtaPrime"}; 
  const char* fnname[7] = {"fMultPion","fMultKaon","fMultEta","fMultRho","fMultOmega","fMultPhi","fMultEtaPrime"};
  Double_t ptparam[7][9];
  Double_t yparam[7][9];
  Double_t nparam[7][9];

  // parameters for 8 TeV generation
  if (fCMSEnergy==kCMS8000GeV) {
    AliInfo ("Using pp parameterization at 8 TeV\n");

  // Parameters of transverse momentum spectra
  Double_t ptparam8000[7][9] = {{1,0.427,2.52,0,0,0,0,0,0},   // pions from 7 TeV section of code
				{1,0.58,2.57,0,0,0,0,0,0},    // kaons from 7 TeV section of code
				{1,0.657,2.685,0,0,0,0,0,0},  // eta from PYTHIA6.4 ATLAS-CSC at 8 TeV
				{1,1.44,3.16,0,0,0,0,0,0},    // rho+omega from 7 TeV section of code
				{1,1.44,3.16,0,0,0,0,0,0},    // rho+omega from 7 TeV section of code
  				{1,1.16,2.74,0,0,0,0,0,0},    // phi from 7 TeV section of code
				{1,0.755,2.578,0,0,0,0,0,0}}; // etaPrime from PYTHIA 6.4 ATLAS-CSC at 8 TeV

  // Parameters of rapidity spectra
  Double_t yparam8000[7][9] = {{1,0.8251,3.657,0,0,0,0,0,0},  // pions from 7 TeV section of code
			       {1,1.83,2.698,0,0,0,0,0,0},    // kaons from 7 TeV section of code
			       {1,0.0509,3.96,0,0,0,0,0,0},   // eta from PYTHIA6.4 ATLAS-CSC at 8 TeV
			       {1,0.0489,3.961,0,0,0,0,0,0},  // rho from PYTHIA6.4 ATLAS-CSC at 8 TeV
			       {1,0.0650,3.966,0,0,0,0,0,0},  // omega from PYTHIA 6.4 ATLAS-CSC at 8 TeV
			       {1,1.279,2.745,0,0,0,0,0,0},   // phi from from PYTHIA6.4 ATLAS-CSC at 8 TeV
			       {1,0.1627,3.883,0,0,0,0,0,0}}; // eta prime from PYTHIA6.4 ATLAS-CSC at 8 TeV

  // Parameters of multiplicity spectra
  Double_t nparam8000[7][9] = {{353.582, 6.76263, 1.66979, 998.445, 9.73281, 12.6704, 175.187, 29.08, 40.2531}, //pions from 7 TeV section of code
                               {1.e4,    0.2841, 0,0,0,0,0,0,0},     // kaons from 7 TeV section of code
                               {2.279e4, 0.2622, 0,0,0,0,0,0,0},  // eta from PYTHIA6.4 ATLAS-CSC at 8 TeV
                               {1.564e4, 0.1713, 0,0,0,0,0,0,0},  // rho from PYTHIA6.4 ATLAS-CSC at 8 TeV
                               {1.662e4, 0.183,  0,0,0,0,0,0,0},   // omega from PYTHIA6.4 ATLAS-CSC at 8 TeV 
                               {6.723e4, 1.121,  0,0,0,0,0,0,0},  // phi from PYTHIA6.4 ATLAS-CSC at 8 TeV
                               {5.005e4, 0.6971, 0,0,0,0,0,0,0}}; // eta prime from PYTHIA6.4 ATLAS-CSC at 8 TeV

    for (Int_t i=0; i<fgkNpart; i++) {
      for (Int_t j=0; j<9; j++) {
	ptparam[i][j] = ptparam8000[i][j];
	yparam[i][j] = yparam8000[i][j];
	nparam[i][j] = nparam8000[i][j];
      }
    }
  }
	
  // parameters for 7 TeV generation
  else if (fCMSEnergy==kCMS7000GeV) {
    AliInfo ("Using pp parameterization at 7 TeV\n");  
    Double_t ptparam7000[7][9] = {{1,0.427,2.52,0,0,0,0,0,0}, // pions from Pythia
				  {1,0.58,2.57,0,0,0,0,0,0},  // kaons from Pythia
				  {1,0.641,2.62,0,0,0,0,0,0}, // eta from Pythia
				  {1,1.44,3.16,0,0,0,0,0,0},  // rho+omega from ALICE muon  
				  {1,1.44,3.16,0,0,0,0,0,0},  // rho+omega from ALICE muon  
				  {1,1.16,2.74,0,0,0,0,0,0},  // phi from ALICE muon  
				  {1,0.72,2.5,0,0,0,0,0,0}};  // etaPrime from Pythia    
		
    Double_t yparam7000[7][9] = {{1,0.8251,3.657,0,0,0,0,0,0}, // pions from pythia
				 {1,1.83,2.698,0,0,0,0,0,0},   // kaons from pythia
				 {1,1.169,3.282,0,0,0,0,0,0},  // eta from pythia
				 {1,1.234,3.264,0,0,0,0,0,0},  // rho from pythia
				 {1,1.311,3.223,0,0,0,0,0,0},  // omega from pythia
				 {1,2.388,2.129,0,0,0,0,0,0},  // phi from pythia
				 {1,1.13,3.3,0,0,0,0,0,0}};    // eta prime from pythia
		
    // multiplicity parameters from pythia
    Double_t nparam7000[7][9] = {{353.582, 6.76263, 1.66979, 998.445, 9.73281, 12.6704, 175.187, 29.08, 40.2531},
				 {1.e4,  0.2841, 0,0,0,0,0,0,0},
				 {1.e4,  0.2647, 0,0,0,0,0,0,0},
				 {7055,  0.1786, 0,0,0,0,0,0,0},
				 {7500,  0.1896, 0,0,0,0,0,0,0},
				 {5.e4,  1.167,  0,0,0,0,0,0,0}, 
				 {2.9e4, 0.714,  0,0,0,0,0,0,0}};
		
    for (Int_t i=0; i<fgkNpart; i++) { 
      for (Int_t j=0; j<9; j++) {
	ptparam[i][j] = ptparam7000[i][j];
	yparam[i][j] = yparam7000[i][j];
	nparam[i][j] = nparam7000[i][j];	
      }
    }
  }  
  else if (fCMSEnergy==kCMS5020GeVpPb || fCMSEnergy==kCMS5020GeVPbp) {
    AliInfo ("Using pPb parameterization at 5.02 TeV\n");  
    Double_t ptparam5020[7][9] = {{1,0.427,2.52,0,0,0,0,0,0}, // pions from Pythia at 7 TeV
				  {1,0.58,2.57,0,0,0,0,0,0},  // kaons from Pythia at 7 TeV
				  {1,0.665,2.796,0,0,0,0,0,0}, // eta from Pythia at 5.02 TeV
				  {1,1.66,3.12,0,0,0,0,0,0},  // rho+omega from ALICE muon  
				  {1,1.66,3.12,0,0,0,0,0,0},  // rho+omega from ALICE muon  
				  {1,2.03,3.13,0,0,0,0,0,0},  // phi from ALICE muon  
				  {1,0.767,2.713,0,0,0,0,0,0}};  // etaPrime from Pythia at 5.02 TeV
		
    Double_t yparam5020[7][9] = {{1,0.8251,3.657,0,0,0,0,0,0}, // pions from pythia at 7 TeV
				 {1,1.83,2.698,0,0,0,0,0,0},   // kaons from pythia at 7 TeV
				 {1,1.169,3.282,0,0,0,0,0,0},  // eta from pythia at 7 TeV
				 {1,1.234,3.264,0,0,0,0,0,0},  // rho from pythia at 7 TeV
				 {1,1.311,3.223,0,0,0,0,0,0},  // omega from pythia at 7 TeV
				 {1,2.388,2.129,0,0,0,0,0,0},  // phi from pythia at 7 TeV
				 {1,1.13,3.3,0,0,0,0,0,0}};    // eta prime from pythia at 7 TeV
		
    // multiplicity parameters from pythia at 7 TeV
    Double_t nparam5020[7][9] = {{353.582, 6.76263, 1.66979, 998.445, 9.73281, 12.6704, 175.187, 29.08, 40.2531},
				 {1.e4,  0.2841, 0,0,0,0,0,0,0},
				 {1.e4,  0.2647, 0,0,0,0,0,0,0},
				 {7055,  0.1786, 0,0,0,0,0,0,0},
				 {7500,  0.1896, 0,0,0,0,0,0,0},
				 {5.e4,  1.167,  0,0,0,0,0,0,0}, 
				 {2.9e4, 0.714,  0,0,0,0,0,0,0}};
		
    for (Int_t i=0; i<fgkNpart; i++) { 
      for (Int_t j=0; j<9; j++) {
	ptparam[i][j] = ptparam5020[i][j];
	yparam[i][j] = yparam5020[i][j];
	nparam[i][j] = nparam5020[i][j];	
      }
    }
    if (fCMSEnergy==kCMS5020GeVpPb) fYCM = -0.4654;   
    else fYCM = 0.4654;
  }  
  else if (fCMSEnergy==kCMS2760GeV){
    // parameters for 2.76 generation
    // pt params has been determined as <pt>ALICE_2.76 = <pt>ALICE_7 * <pt>PYTHIA_2.76 / <pt>PYTHIA_7
    AliInfo ("Using pp parameterization at 2.76 TeV\n");  
    Double_t yparam2760[7][9] = {{1,0.8251,3.657,0,0,0,0,0,0},      // pions from pythia
				 {1,1.83,2.698,0,0,0,0,0,0},        // kaons from pythia
				 {1,0.011,3.474,0,0,0,0,0,0},       // eta from pythia
				 {1,-0.01,3.409,0,0,0,0,0,0},       // rho from pythia
				 {1,-0.037,3.294,0,0,0,0,0,0},      // omega from pythia
				 {1,-0.016,2.717,0,0,0,0,0,0},      // phi from pythia
				 {1,-0.010,3.312,0,0,0,0,0,0}};     // eta prime from pythia  
		
    Double_t ptparam2760[7][9] = {{1,0.1665,8.878,0,0,0,0,0,0},  // pions from Pythia
				  {1,0.1657,8.591,0,0,0,0,0,0},  // kaons from Pythia
				  {1,0.641,2.62,0,0,0,0,0,0},    // eta from ALICE 7 TeV
				  {1,1.3551,3.16,0,0,0,0,0,0},   // rho with <pt> scaled 
				  {1,1.3551,3.16,0,0,0,0,0,0},   // omega with <pt> scaled 
				  {1,1.0811,2.74,0,0,0,0,0,0},   // phi with <pt> scaled 
				  {1,0.72,2.5,0,0,0,0,0,0}};     // etaPrime from ALICE 7 TeV
		
    Double_t nparam2760[7][9] = {{9752,-2.693,3.023,9.5e9,-84.68,16.75,-14.06,635.3,-423.2}, // pions
				 {1.e5, 1.538,  0,0,0,0,0,0,0},                              // kaons
				 {1.e4, 0.351,  0,0,0,0,0,0,0},                              // eta
				 {1.e4, 0.2471, 0,0,0,0,0,0,0},                              // rho
				 {1.e4, 0.2583, 0,0,0,0,0,0,0},                              // omega
				 {1.e5, 1.393,  0,0,0,0,0,0,0},                              // phi
				 {1.e4, 0.9005, 0,0,0,0,0,0,0}};                             // etaPrime
		
    for (Int_t i=0; i<fgkNpart; i++) { 
      for (Int_t j=0; j<9; j++) {
	ptparam[i][j] = ptparam2760[i][j];
	yparam[i][j] = yparam2760[i][j];
	nparam[i][j] = nparam2760[i][j];	
      }
    }
  } 
  else AliFatal("Energy not correctly defined");
	
  for (Int_t i=0; i<fgkNpart; i++) { 
    fPDG[i] = pdg[i]; 
    if (i!=0) { 
      fMult[i] = new TF1(fnname[i],"[0]*exp(-[1]*x)",0,30);
      fMult[i]->SetParameters(nparam[i][0],nparam[i][1]);  
    }
    else { 
      fMult[i] = new TF1(fnname[i],"gaus(0)+gaus(3)+gaus(6)",0,150);
      for (Int_t j=0; j<9; j++) fMult[i]->SetParameter(j,nparam[i][j]);
    }
		
    fPt[i] = new TF1(fptname[i],AliGenMUONLMR::PtDistr,0,20,3);
    fPt[i]->SetParameters(ptparam[i][0], ptparam[i][1], ptparam[i][2]);  
    fY[i] = new TF1(fyname[i],AliGenMUONLMR::YDistr,-10,10,4);
    fY[i]->SetParameters(yparam[i][0], yparam[i][1], yparam[i][2],fYCM); 
  }
	
  for(Int_t i = 0; i<2; i++){
    fDecay[i] = new TF1(fdname[i],"exp(-x/[0])",0,150);
    fDecay[i]->SetParameter(0,ctau[i]);
  }
	
  for (Int_t ipart = 0; ipart < fgkNpart; ipart++) { 
    fParticle[ipart] = new TParticle(); 
    fParticle[ipart]->SetPdgCode(fPDG[ipart]); 
  }
	
  TDatabasePDG *pdgdb = TDatabasePDG::Instance(); 
  Double_t mumass = pdgdb->GetParticle(13)->Mass();
  fMu[0] = new TParticle(); 
  fMu[0]->SetPdgCode(-13); 
  fMu[0]->SetCalcMass(mumass); 
  fMu[1] = new TParticle(); 
  fMu[1]->SetPdgCode(13); 
  fMu[1]->SetCalcMass(mumass); 
	
  // function for polarized theta distributions
  fCosTheta = new TF1 ("fCosTheta","1+[0]*x*x",-1,1);
  fCosTheta->SetParameter(0,1);
	
  // Dalitz decays 
  Int_t nbins = 1000;
  Double_t xmin = 0, xmax = 2; 
  fDalitz[0] = new TH1F("hDalitzEta","",nbins,xmin,xmax);
  fDalitz[1] = new TH1F("hDalitzOmega","",nbins,xmin,xmax);
  fDalitz[2] = new TH1F("hDalitzEtaPrime","",nbins,xmin,xmax);
	
  Double_t meta   = pdgdb->GetParticle("eta")->Mass(); 
  Double_t momega = pdgdb->GetParticle("omega")->Mass(); 
  Double_t metaPrime = pdgdb->GetParticle("eta'")->Mass(); 
  Double_t mpi0   = pdgdb->GetParticle("pi0")->Mass(); 
  Double_t md3 = 0, mres = 0; 
	
  for (Int_t index = 0; index < 3; index++) { 
    if (index == 0) { 
      mres = meta; 
      md3 = 0; 
    }
    else if (index == 1) { 
      mres = momega; 
      md3 = mpi0; 
    }
    else if (index == 2) { 
      mres = metaPrime; 
      md3 = 0; 
    }
    Double_t delta   = md3 * md3 / (mres * mres);
    Double_t epsilon = mumass * mumass / (mres * mres);
    nbins = fDalitz[index]->GetNbinsX();
    xmin = fDalitz[index]->GetXaxis()->GetXmin(); 
    Double_t deltax =  fDalitz[index]->GetBinWidth(1);
    Double_t xd = xmin - deltax/2.; 
    for (Int_t ibin = 0; ibin< nbins; ibin++) { 
      Double_t dalval = 0; 
      xd += deltax; 
      if (xd > 4. *epsilon) { 
	Double_t bracket = TMath::Power(1. + xd/(1. - delta),2)      
	  - 4. * xd / ((1. - delta) * (1. - delta));
	if (bracket > 0) { 
	  dalval = TMath::Power(bracket,1.5) /xd *
	    TMath::Sqrt(1 - 4 * epsilon / xd) * (1 + 2 * epsilon / xd) * 
	    FormFactor(xd * mres * mres, index);
	  fDalitz[index]->Fill(xd,dalval); 
	}
      }
    }
  }
	
  fRhoLineShape = new TF1("fRhoLineShape",RhoLineShapeNew,0,2,2); 
  fHMultMu = new TH1D("fHMultMu","Muon multiplicity",20,-0.5,19.5); 
  fHNProc = new TH1D("fHNProc","Number of gen. evts. per process in 4 pi",9,-0.5,8.5); 
	
}
//-----------------------------------------------------------

AliGenMUONLMR::AliGenMUONLMR (AliGenMUONLMR &gen) : AliGenMC(), 
						    fNMuMin(gen.fNMuMin), 
						    fCMSEnergy(gen.fCMSEnergy), 
						    fGenSingleProc(gen.fGenSingleProc),
						    fYCM(gen.fYCM),
						    fCosTheta(gen.fCosTheta), 
						    fRhoLineShape(gen.fRhoLineShape),  
						    fHMultMu(gen.fHMultMu), 
						    fHNProc(gen.fHNProc) {  
  for (Int_t i=0; i < fgkNpart; i++) { 
    fPDG[i] = gen.fPDG[i]; 
    fScaleMult[i] = gen.fScaleMult[i]; 
    fPt[i] = (TF1*) gen.fPt[i]->Clone(); 
    fY[i] = (TF1*) gen.fY[i]->Clone();  
    fMult[i] = (TF1*) gen.fMult[i]->Clone(); 
    fParticle[i] = (TParticle*) gen.fParticle[i]->Clone(); 
  }
  
  for(Int_t i = 0; i<2; i++) fDecay[i] = (TF1*) gen.fDecay[i]->Clone(); 
  for(Int_t i = 0; i<3; i++) fDalitz[i] = (TH1F*) gen.fDalitz[i]->Clone(); 
  for(Int_t i = 0; i<2; i++) fMu[i] = (TParticle*) gen.fMu[i]->Clone(); 
}

//-----------------------------------------------------------

AliGenMUONLMR& AliGenMUONLMR::operator=(const AliGenMUONLMR &gen) {
  // Assignment operator
  if (this!=&gen) {
    fNMuMin = gen.fNMuMin; 
    fCMSEnergy = gen.fCMSEnergy; 
    fGenSingleProc = gen.fGenSingleProc; 
    fYCM = gen.fYCM; 
    fCosTheta = (TF1*) gen.fCosTheta->Clone();  
    fRhoLineShape = (TF1*) gen.fRhoLineShape->Clone();
    fHMultMu = (TH1D*) gen.fHMultMu->Clone();
    fHNProc = (TH1D*) gen.fHNProc->Clone();  
    
    for (Int_t i=0; i < fgkNpart; i++) { 
      fPDG[i] = gen.fPDG[i]; 
      fScaleMult[i] = gen.fScaleMult[i]; 
      fPt[i] = (TF1*) gen.fPt[i]->Clone(); 
      fY[i] = (TF1*) gen.fY[i]->Clone();  
      fMult[i] = (TF1*) gen.fMult[i]->Clone(); 
      fParticle[i] = (TParticle*) gen.fParticle[i]->Clone(); 
    }
    
    for(Int_t i = 0; i<2; i++) fDecay[i] = (TF1*) gen.fDecay[i]->Clone(); 
    for(Int_t i = 0; i<3; i++) fDalitz[i] = (TH1F*) gen.fDalitz[i]->Clone(); 
    for(Int_t i = 0; i<2; i++) fMu[i] = (TParticle*) gen.fMu[i]->Clone(); 
  }
  return *this; 
}

 
//-----------------------------------------------------------

AliGenMUONLMR::~AliGenMUONLMR()
{
  // Default destructor
  for (Int_t i=0; i<7; i++) { 
    delete fPt[i]; 
    delete fY[i]; 
    delete fMult[i]; 
    delete fParticle[i]; 
  }    
  
  for (Int_t i=0; i<2; i++) { 
    delete fDecay[i]; 
    delete fMu[i]; 
  }

  for (Int_t i=0; i<3; i++) delete fDalitz[i]; 

  delete fCosTheta; fCosTheta = 0;  
  delete fRhoLineShape; fRhoLineShape = 0;  
  delete fHMultMu; fHMultMu = 0;   
  delete fHNProc;  fHNProc = 0;   
}

//-----------------------------------------------------------

void AliGenMUONLMR::FinishRun(){ 
  // save some histograms to an output file 
  Int_t nbins = fHNProc->GetNbinsX(); 
  for (Int_t ibin=1; ibin <= nbins; ibin++) AliInfo (Form("ibin = %d nEvProc = %g",
							  ibin,fHNProc->GetBinContent(ibin)));
  TFile *fout = new TFile("AliGenMUONLMR_histos.root","recreate"); 
  fHMultMu->Write(); 
  fHNProc->Write(); 
  fout->Close(); 
}

//-----------------------------------------------------------

Double_t AliGenMUONLMR::YDistr(Double_t *px, Double_t *par){ 
  // function for rapidity distribution: plateau at par[0] +
  // gaussian tails centered at par[1] and with par[2]=sigma  
  Double_t ylab = px[0]; 
  Double_t y0 = par[3]; // center of mass rapidity  
  Double_t func = 0;
  if (ylab<y0+par[1] && ylab>y0-par[1]) func = par[0]; 
  else if (ylab>y0+par[1]) { 
    Double_t z = (ylab-(par[1]+y0) )/(par[2]); 
    func = par[0] * TMath::Exp(-0.5 * z * z); 
  }
  else {
    Double_t z = (ylab-(-par[1]+y0) )/(par[2]); 
    func = par[0] * TMath::Exp(-0.5 * z * z); 
  }
  return func; 
}

//-----------------------------------------------------------

Double_t AliGenMUONLMR::PtDistr(Double_t *px, Double_t *par){
  // pt distribution: power law 
  Double_t x = px[0];
  Double_t func = par[0] * x / TMath::Power((1+(x/par[1])*(x/par[1])),par[2]); 
  return func; 
}

//-----------------------------------------------------------

void AliGenMUONLMR::Generate() {
  //
  // generate the low mass resonances and their decays according to  
  // the multiplicity parameterized by pythia and BR from PDG  
  // rapidity distributions parametrized from pythia 
  // pt distributions from data (or pythia for etaprime) 
  //
  Double_t pxPushed[100], pyPushed[100], pzPushed[100], ePushed[100]; 
  Int_t nmuons = -1, npartPushed = 0, pdgPushed[100]; 
  Double_t polar[3]= {0,0,0};  // Polarisation of the parent particle (for GEANT tracking)
  Double_t origin0[3];         // Origin of the generated parent particle (for GEANT tracking)
  // Calculating vertex position per event
  for (Int_t j=0;j<3;j++) origin0[j]=fOrigin[j];
  if(fVertexSmear==kPerEvent) {
    Vertex();
    for (Int_t j=0;j<3;j++) origin0[j]=fVertex[j];
  }
  
  TParticle *mother; 
  TDatabasePDG* pdg = TDatabasePDG::Instance();

  Double_t pt, y, phi, mass, px, py, pz, ene, mt; 

  const Int_t nproc = 9; 
  Int_t idRes[nproc] = {kEtaLMR, kEtaLMR, kRhoLMR, kOmegaLMR, kOmegaLMR, kPhiLMR, kEtaPrimeLMR, kPionLMR, kKaonLMR}; 
  Double_t BR[nproc] = {5.8e-6, 3.1e-4, 4.55e-5, 7.28e-5, 1.3e-4, 2.86e-4, 1.04e-4, 1, 0.6344};
  //  Double_t BR[nproc] = {1, 1, 1, 1, 1, 1, 1, 1, 1};
  Int_t idDec[nproc] = {0, 1, 0, 0, 1, 0, 1, 2, 2}; // 0:2body, 1:Dalitz, 2:pi/K 
  Int_t mult[nproc] = {0,0,0,0,0,0,0,0,0}; 

  while (nmuons < fNMuMin) { 

    nmuons = 0; 
    npartPushed = 0; 
    for (Int_t iproc=0; iproc<nproc; iproc++) { 
      if (fGenSingleProc == -1) { 
	mult[iproc] = Int_t(fMult[idRes[iproc]]->GetRandom()*fScaleMult[idRes[iproc]]); 
      }
      else { 
	if (iproc==fGenSingleProc) { 
	  mult[iproc] = 1; 
	  BR[iproc] = 1;
	} 
	else { 
	  mult[iproc] = 0; 
	  BR[iproc] = 0;
	}
      }
    }
    
    if (fGenSingleProc == -1) { 
      mult[1] = mult[0]; 
      mult[4] = mult[3]; 
    }
    
    for (Int_t iproc = 0; iproc < nproc; iproc++) { 
      //       printf ("Multiplicity for process %d is %d\n",iproc,mult[iproc]); 
      for (Int_t imult=0; imult<mult[iproc]; imult++) { 
	if (gRandom->Rndm() < BR[iproc]) { 
	  fHNProc->Fill(iproc); 
	  Int_t ipart = idRes[iproc]; 
	  pt  = fPt[ipart]->GetRandom(); 
	  y   = fY[ipart]->GetRandom(); 
	  phi = gRandom->Rndm() * 2 * TMath::Pi(); 
	  mass = pdg->GetParticle(fPDG[ipart])->Mass(); 
	  px  = pt * TMath::Cos(phi); 
	  py  = pt * TMath::Sin(phi); 
	  mt  = TMath::Sqrt(pt * pt + mass * mass);
	  pz  = mt * TMath::SinH(y); 
	  ene = mt * TMath::CosH(y); 
	
	  mother = fParticle[ipart]; 
	  mother->SetMomentum(px,py,pz,ene); 
	  mother->SetCalcMass(mass);
	  if (!KinematicSelection(mother,0)) continue; 

	  Bool_t hasDecayed = kTRUE;
	  if (idDec[iproc] == 0) Decay2Body(mother);
	  else if (idDec[iproc] == 1) DalitzDecay(mother); 
	  else DecayPiK(mother,hasDecayed); 
	  if (!hasDecayed) continue; 
	  Bool_t isMu0Acc = KinematicSelection(fMu[0],1); 
	  Bool_t isMu1Acc = KinematicSelection(fMu[1],1); 
	  Bool_t isMuFromPiKAcc = kTRUE; 

	  if (idDec[iproc] == 2) isMuFromPiKAcc = (mother->GetPdgCode()>0) ? isMu0Acc : isMu1Acc;
	  // mother 
	  if ((idDec[iproc]  < 2 && (isMu0Acc || isMu1Acc)) || 
	      (idDec[iproc] == 2 && isMuFromPiKAcc)) { 
	    pdgPushed[npartPushed] = mother->GetPdgCode(); 
	    pxPushed[npartPushed] = mother->Px(); 
	    pyPushed[npartPushed] = mother->Py(); 
	    pzPushed[npartPushed] = mother->Pz();
	    ePushed[npartPushed] = mother->Energy(); 
	    npartPushed++; 
	    if (isMu0Acc && (idDec[iproc] < 2 || mother->GetPdgCode() > 0)) { 
	      pdgPushed[npartPushed] = fMu[0]->GetPdgCode(); 
	      pxPushed[npartPushed] = fMu[0]->Px(); 
	      pyPushed[npartPushed] = fMu[0]->Py(); 
	      pzPushed[npartPushed] = fMu[0]->Pz();
	      ePushed[npartPushed] = fMu[0]->Energy(); 
	      npartPushed++; 
	      nmuons++; 
	    }
	    
	    if (isMu1Acc && (idDec[iproc] < 2 || mother->GetPdgCode() < 0)) { 
	      pdgPushed[npartPushed] = fMu[1]->GetPdgCode(); 
	      pxPushed[npartPushed] = fMu[1]->Px(); 
	      pyPushed[npartPushed] = fMu[1]->Py(); 
	      pzPushed[npartPushed] = fMu[1]->Pz();
	      ePushed[npartPushed] = fMu[1]->Energy(); 
	      npartPushed++; 
	      nmuons++; 
	    }
	  } 
	} // end if BR
      } // end loop on multiplicity 
    }  // end loop on process 
    fHMultMu->Fill(nmuons); 
  } // keep on generating until at least a muon is created in the event
  
  Int_t ntmother = 0, ntchild =0; 
  for (Int_t ipart = 0; ipart < npartPushed; ipart++) { 
    if (TMath::Abs(pdgPushed[ipart]) != 13) { // particle is not a muon, hence it's a mother
      PushTrack(0,-1,pdgPushed[ipart],
		pxPushed[ipart],pyPushed[ipart],pzPushed[ipart],ePushed[ipart],
		origin0[0],origin0[1],origin0[2],0.,
		polar[0],polar[1],polar[2],
		kPPrimary,ntmother,1,11);
      KeepTrack(ntmother); 
    }
    else { 
      PushTrack(1,ntmother,pdgPushed[ipart],
		pxPushed[ipart],pyPushed[ipart],pzPushed[ipart],ePushed[ipart],
		origin0[0],origin0[1],origin0[2],0.,
		polar[0],polar[1],polar[2],
		kPDecay,ntchild,1,1);
      KeepTrack(ntchild); 
    }
  }
  SetHighWaterMark(ntchild); 
  AliGenEventHeader* header = new AliGenEventHeader("LMR");
  header->SetPrimaryVertex(fVertex);
  header->SetNProduced(fNprimaries);
  AddHeader(header); 
}

//------------------------------------------------------------------

void AliGenMUONLMR::Decay2Body(TParticle *mother){ 
  // performs decay in two muons of the low mass resonances
  Double_t md1 = fMu[0]->GetMass(); 
  Int_t pdg = mother->GetPdgCode(); 
  Double_t mres =0; 
  // if mother is a rho, extract the mass from its line shape
  // otherwise consider the resonance mass 
  if (pdg == 113) mres = fRhoLineShape->GetRandom(); 
  else mres = mother->GetCalcMass(); 
  //  while (mres < md1 + md2) mres =  fDsigmaDm[res]->GetRandom();
  // energies and momenta in rest frame 
  Double_t e1 = mres / 2.;
  Double_t p1 = TMath::Sqrt((e1 + md1)*(e1 - md1)); 
  // orientation in decaying particle rest frame
  Double_t costheta = gRandom->Rndm() * 2 - 1;
  Double_t sintheta = TMath::Sqrt((1. + costheta)*(1. - costheta));
  Double_t phi      = 2. * TMath::Pi() * gRandom->Rndm(); 
  Double_t px1      = p1 * sintheta * TMath::Cos(phi); 
  Double_t py1      = p1 * sintheta * TMath::Sin(phi); 
  Double_t pz1      = p1 * costheta; 

  // boost muons into lab frame 

  TLorentzVector vmother, v1, v2;
  //  TLorentzVector boosted1, boosted2;   
  vmother.SetPxPyPzE(mother->Px(),mother->Py(),mother->Pz(),mother->Energy());
  v1.SetPxPyPzE(px1,py1,pz1,e1); 
  v2.SetPxPyPzE(-px1,-py1,-pz1,e1); 

  TVector3 betaParent = (1./vmother.E())*vmother.Vect(); // beta = p/E
  v1.Boost(betaParent);
  v2.Boost(betaParent);

  //   TLorentzVector vtot = v1 + v2; 
  //   printf ("mother: %g   %g    %g     %g\n",vmother.Px(), vmother.Py(), vmother.Pz(), vmother.E());
  //   printf ("vtot  : %g   %g    %g     %g\n",vtot.Px(), vtot.Py(), vtot.Pz(), vtot.E());

  fMu[0]->SetMomentum(v1.Px(),v1.Py(),v1.Pz(),v1.E());
  fMu[1]->SetMomentum(v2.Px(),v2.Py(),v2.Pz(),v2.E());
} 

//------------------------------------------------------------------

void AliGenMUONLMR::DecayPiK(TParticle *mother, Bool_t &hasDecayed){ 
  // performs decays of pions and kaons
  Double_t md1 = fMu[0]->GetMass(); 
  // extract the mass from the resonance's line shape
  Double_t mres = mother->GetMass(); 
  // choose the pi/k sign, assuming 50% probabilities for both signs
  Int_t sign = (gRandom->Rndm() > 0.5) ? 1 : -1;
  mother->SetPdgCode(sign * TMath::Abs(mother->GetPdgCode())); 

  // energies and momenta in rest frame 
  Double_t e1 = (mres*mres + md1*md1)/(2*mres);
  Double_t p1 = TMath::Sqrt((e1 + md1)*(e1 - md1)); 
  // orientation in decaying particle rest frame
  Double_t costheta = gRandom->Rndm() * 2 - 1;
  Double_t sintheta = TMath::Sqrt((1. + costheta)*(1. - costheta));
  Double_t phi      = 2. * TMath::Pi() * gRandom->Rndm(); 
  Double_t px1      = p1 * sintheta * TMath::Cos(phi); 
  Double_t py1      = p1 * sintheta * TMath::Sin(phi); 
  Double_t pz1      = p1 * costheta;  

  // boost muons into lab frame 
  TLorentzVector vmother, v1;
  vmother.SetPxPyPzE(mother->Px(),mother->Py(),mother->Pz(),mother->Energy());
  v1.SetPxPyPzE(px1,py1,pz1,e1); 

  TVector3 betaParent = (1./vmother.E())*vmother.Vect(); // beta = p/E
  v1.Boost(betaParent);  
  if (mother->GetPdgCode()>0) fMu[0]->SetMomentum(v1.Px(),v1.Py(),v1.Pz(),v1.E());
  else fMu[1]->SetMomentum(v1.Px(),v1.Py(),v1.Pz(),v1.E());

  Int_t idmother = 0; 
  if (TMath::Abs(mother->GetPdgCode())== 211) idmother = 0; 
  if (TMath::Abs(mother->GetPdgCode())== 321) idmother = 1; 
  Double_t gammaRes = mother->Energy()/mres;
  Double_t zResCM = fDecay[idmother]->GetRandom();
  Double_t zResLab = gammaRes*zResCM;  
  if(zResLab > 0.938) hasDecayed = 0; // 0.938: distance from IP to absorber + lambda_i
  else hasDecayed = 1;
} 

//-------------------------------------------------------------------

void AliGenMUONLMR::DalitzDecay(TParticle *mother){
  //
  // perform dalitz decays of eta, omega and etaprime 
  //
  //in the rest frame of the virtual photon:
  Double_t mres = mother->GetCalcMass(); 
  Double_t mumass  = fMu[0]->GetMass(); 
  Double_t md3  = 0;  // unless differently specified, third particle is a photon 
  if (mother->GetPdgCode() == 223) md3 = 0.134977; // if mother is an omega, third particle is a pi0
  Int_t index = 0; 
  if (mother->GetPdgCode() == 221) index = 0;  // eta
  else if (mother->GetPdgCode() == 223) index = 1; // omega  
  else if (mother->GetPdgCode() == 331) index = 2; // etaPrime  
  Int_t flag = 0; 
  Double_t xd=0, mvirt2=0; 
  Double_t countIt = 0;
  while (flag==0) {  
    xd       = fDalitz[index]->GetRandom(); 
    mvirt2   = xd * mres * mres;   // mass of virtual photon 
    // check kinematics 
    if (mres - md3 > TMath::Sqrt(mvirt2) && TMath::Sqrt(mvirt2)/2. > mumass) flag=1;
    if (++countIt>1E11) {
      mvirt2 =  mres * mres * 0.998; 
      break;
    }
  }  

  //
  //        Generate muons in virtual photon rest frame. 
  //        z axis is the virt. photon direction (before boost)  
  //

  Double_t e1 = TMath::Sqrt(mvirt2)/2.; // energy of mu1 in the virtual photon frame
  Double_t psquare = (e1 + mumass)*(e1 - mumass); 
  if (psquare<0) {
    AliError(Form("sqrt of psquare = %f put to 0\n",psquare)); 
    psquare = 0;
  }
  Double_t p1 = TMath::Sqrt(psquare);
  //theta angle between the pos. muon and the virtual photon 
  Double_t costheta = fCosTheta->GetRandom();
  if (costheta>1)  costheta = 1; 
  if (costheta<-1) costheta = -1; 
  Double_t sintheta = TMath::Sqrt((1. + costheta)*(1. - costheta));
  Double_t phi      = 2 * TMath::Pi() * gRandom->Rndm();
  Double_t sinphi   = TMath::Sin(phi);
  Double_t cosphi   = TMath::Cos(phi);

  // fill 4-vectors of leptons in the virtual photon frame

  Double_t px1 = p1*sintheta*cosphi; 
  Double_t py1 = p1*sintheta*sinphi; 
  Double_t pz1 = p1*costheta; 
  Double_t px2 = -p1*sintheta*cosphi; 
  Double_t py2 = -p1*sintheta*sinphi; 
  Double_t pz2 = -p1*costheta; 
  Double_t e2  = e1; 

  fMu[0]->SetMomentum(px1,py1,pz1,e1); 
  fMu[1]->SetMomentum(px2,py2,pz2,e2); 

  // calculate components of non-dilepton in CMS of parent resonance 

  Double_t e3 = (mres * mres + md3 * md3 - mvirt2) / (2.*mres);
  Double_t psquare3 = (e3 + md3)*(e3 - md3); 
  if (psquare3<0) {
    AliError(Form("Sqrt of psquare3 = %f put to 0\n",psquare3)); 
    psquare3 = 0;
  }
  Double_t p3 = TMath::Sqrt(psquare3);
  Double_t costheta2 = 2.* gRandom->Rndm() - 1.;   // angle between virtual photon and resonance
  if (costheta2>1)  costheta2 = 1; 
  if (costheta2<-1) costheta2 = -1; 
  Double_t sintheta2 = TMath::Sqrt((1. + costheta2)*(1. - costheta2));
  Double_t phi2      = 2 * TMath::Pi() * gRandom->Rndm();
  Double_t sinphi2   = TMath::Sin(phi2);
  Double_t cosphi2   = TMath::Cos(phi2);
  Double_t px3 = p3*sintheta2*cosphi2; 
  Double_t py3 = p3*sintheta2*sinphi2; 
  Double_t pz3 = p3*costheta2; 
  TLorentzVector v3(px3,py3,pz3,e3); 

  sintheta2 = -sintheta2;
  cosphi2   = -cosphi2;
  sinphi2   = -sinphi2;

  Double_t px1new = px1*costheta2*cosphi2 - py1*sinphi2 + pz1*sintheta2*cosphi2; 
  Double_t py1new = px1*costheta2*sinphi2 + py1*cosphi2 + pz1*sintheta2*sinphi2; 
  Double_t pz1new =-px1*sintheta2                       + pz1*costheta2; 
  Double_t px2new = px2*costheta2*cosphi2 - py2*sinphi2 + pz2*sintheta2*cosphi2; 
  Double_t py2new = px2*costheta2*sinphi2 + py2*cosphi2 + pz2*sintheta2*sinphi2; 
  Double_t pz2new =-px2*sintheta2                       + pz2*costheta2; 

  fMu[0]->SetMomentum(px1new,py1new,pz1new,e1); 
  fMu[1]->SetMomentum(px2new,py2new,pz2new,e2); 

  Double_t evirt = mres - e3; 
  Double_t pxvirt = -px3;
  Double_t pyvirt = -py3;
  Double_t pzvirt = -pz3;
  TLorentzVector vvirt(pxvirt,pyvirt,pzvirt,evirt); 
  TVector3 betaVirt = (1./evirt) * vvirt.Vect(); // virtual photon beta in res frame

  TLorentzVector v1(px1,py1,pz1,e1); 
  TLorentzVector v2(px2,py2,pz2,e2);

  // boost the muons in the frame where the resonance is at rest 

  v1.Boost(betaVirt); 
  v2.Boost(betaVirt); 

  // boost muons and third particle in lab frame

  TLorentzVector vmother(mother->Px(), mother->Py(), mother->Pz(), mother->Energy());  
  TVector3 resBetaLab = (1./vmother.E())*vmother.Vect(); // eta beta in lab frame
  v1.Boost(resBetaLab); 
  v2.Boost(resBetaLab); 
  v3.Boost(resBetaLab); 
  vvirt.Boost(resBetaLab); 

  fMu[0]->SetMomentum(v1.Px(),v1.Py(),v1.Pz(),v1.E());
  fMu[1]->SetMomentum(v2.Px(),v2.Py(),v2.Pz(),v2.E());
  //   part3->SetMomentum(v3.Px(),v3.Py(),v3.Pz(),v3.E());

  //   TLorentzVector vtot = v1 + v2 + v3; 
  //   TLorentzVector vdimu = v1 + v2; 
  //   printf ("mother: %g   %g    %g     %g\n",vmother.Px(), vmother.Py(), vmother.Pz(), vmother.E());
  //   printf ("vtot  : %g   %g    %g     %g\n",vtot.Px(), vtot.Py(), vtot.Pz(), vtot.E());
  //   printf ("vvirt : %g   %g    %g     %g\n",vvirt.Px(), vvirt.Py(), vvirt.Pz(), vvirt.E());
  //   printf ("vdimu : %g   %g    %g     %g\n",vdimu.Px(), vdimu.Py(), vdimu.Pz(), vdimu.E());

}

//------------------------------------------------------------------

Double_t AliGenMUONLMR::FormFactor(Double_t q2, Int_t decay){ 
  //  Calculates the form factor for Dalitz decays A->B+l+l
  //  Returns: |F(q^2)|^2
  //
  //  References: L.G. Landsberg, Physics Reports 128 No.6 (1985) 301-376. 
 
  Double_t ff2, mass2;
  Double_t n2, n4, m2; 
  // Lepton-G
  
  Double_t lambda2inv = 0;
  switch (decay) { 
  case 0:   // eta -> mu mu gamma  
    // eta   -> l+ l- gamma: pole approximation
    lambda2inv = 1.95; 
    mass2 = fParticle[kEtaLMR]->GetMass() * fParticle[kEtaLMR]->GetMass(); 
    if (q2 < mass2) ff2 = TMath::Power(1./(1.-lambda2inv*q2),2);
    else ff2 = 0; 
    break;
  case 1:   // omega -> mu mu pi0 
    // omega -> l+ l- pi0: pole approximation
    mass2 = fParticle[kOmegaLMR]->GetMass() * fParticle[kOmegaLMR]->GetMass(); 
    lambda2inv = 2.26; 
    if (q2 < mass2) ff2 = TMath::Power(1./(1.-lambda2inv*q2),2);
    else ff2 = 0; 
    break;
  case 2:   // etaPrime -> mu mu gamma 
    mass2 = fParticle[kEtaPrimeLMR]->GetMass() * fParticle[kEtaPrimeLMR]->GetMass(); 
    // eta'  -> l+ l- gamma: Breit-Wigner fitted to data
    n2 = 0.764 * 0.764; 
    n4 = n2 * n2; 
    m2 = 0.1020 * 0.1020;
    if (q2 < mass2) ff2 = n4 / (TMath::Power(n2-q2,2) + m2 * n2); 
    else ff2 = 0; 
    break;
  default:
    AliError ("FormFactor: Decay not found"); 
    return 0; 
    break; 
  }
  return ff2; 
}

//____________________________________________________________

Double_t AliGenMUONLMR::RhoLineShapeNew(Double_t *x, Double_t */*para*/){
  //new parameterization implemented by Hiroyuki Sako (GSI)
  Double_t mass = *x;
  double r, GammaTot;
  Double_t mRho    = TDatabasePDG::Instance()->GetParticle("rho0")->Mass();
  Double_t mPi     = TDatabasePDG::Instance()->GetParticle("pi0")->Mass();
  Double_t mMu     = TDatabasePDG::Instance()->GetParticle("mu-")->Mass();
  Double_t Gamma0  = TDatabasePDG::Instance()->GetParticle("rho0")->Width();

  const double Norm = 0.0744416*1.01;  

  // 0.0744416 at m = 0.72297
  // is the max number with Norm=1 (for rho)
  
  double mThreshold = 2.*mPi;

  const double T = 0.170; // Assumption of pi+ temperature [GeV/c^2]
  //const double T = 0.11; // Taken from fit to pi+ temperature [GeV/c^2]
  // with Reference: LEBC-EHS collab., Z. Phys. C 50 (1991) 405

  if (mass < mThreshold) {
    r = 0.;
    return r;
  }

  double k = sqrt(0.25*mass*mass-(mThreshold/2)*(mThreshold/2));
  double k0 = sqrt(0.25*mRho*mRho-(mThreshold/2)*(mThreshold/2));

  GammaTot = (k/k0)*(k/k0)*(k/k0)*(mRho/mass)*(mRho/mass)*Gamma0;

  double FormFactor2 = 1/((mass*mass-mRho*mRho)*(mass*mass-mRho*mRho)+
			  mass*mass*GammaTot*GammaTot);

  r = pow(mass,1.5)*pow((1-mThreshold*mThreshold/(mass*mass)),1.5)*
    ((mass*mass+2*mMu*mMu)/(mass*mass))*(pow((mass*mass-4*mMu*mMu),0.5)/mass)*FormFactor2
    *exp(-mass/T)/Norm;

  return r;
}

 AliGenMUONLMR.cxx:1
 AliGenMUONLMR.cxx:2
 AliGenMUONLMR.cxx:3
 AliGenMUONLMR.cxx:4
 AliGenMUONLMR.cxx:5
 AliGenMUONLMR.cxx:6
 AliGenMUONLMR.cxx:7
 AliGenMUONLMR.cxx:8
 AliGenMUONLMR.cxx:9
 AliGenMUONLMR.cxx:10
 AliGenMUONLMR.cxx:11
 AliGenMUONLMR.cxx:12
 AliGenMUONLMR.cxx:13
 AliGenMUONLMR.cxx:14
 AliGenMUONLMR.cxx:15
 AliGenMUONLMR.cxx:16
 AliGenMUONLMR.cxx:17
 AliGenMUONLMR.cxx:18
 AliGenMUONLMR.cxx:19
 AliGenMUONLMR.cxx:20
 AliGenMUONLMR.cxx:21
 AliGenMUONLMR.cxx:22
 AliGenMUONLMR.cxx:23
 AliGenMUONLMR.cxx:24
 AliGenMUONLMR.cxx:25
 AliGenMUONLMR.cxx:26
 AliGenMUONLMR.cxx:27
 AliGenMUONLMR.cxx:28
 AliGenMUONLMR.cxx:29
 AliGenMUONLMR.cxx:30
 AliGenMUONLMR.cxx:31
 AliGenMUONLMR.cxx:32
 AliGenMUONLMR.cxx:33
 AliGenMUONLMR.cxx:34
 AliGenMUONLMR.cxx:35
 AliGenMUONLMR.cxx:36
 AliGenMUONLMR.cxx:37
 AliGenMUONLMR.cxx:38
 AliGenMUONLMR.cxx:39
 AliGenMUONLMR.cxx:40
 AliGenMUONLMR.cxx:41
 AliGenMUONLMR.cxx:42
 AliGenMUONLMR.cxx:43
 AliGenMUONLMR.cxx:44
 AliGenMUONLMR.cxx:45
 AliGenMUONLMR.cxx:46
 AliGenMUONLMR.cxx:47
 AliGenMUONLMR.cxx:48
 AliGenMUONLMR.cxx:49
 AliGenMUONLMR.cxx:50
 AliGenMUONLMR.cxx:51
 AliGenMUONLMR.cxx:52
 AliGenMUONLMR.cxx:53
 AliGenMUONLMR.cxx:54
 AliGenMUONLMR.cxx:55
 AliGenMUONLMR.cxx:56
 AliGenMUONLMR.cxx:57
 AliGenMUONLMR.cxx:58
 AliGenMUONLMR.cxx:59
 AliGenMUONLMR.cxx:60
 AliGenMUONLMR.cxx:61
 AliGenMUONLMR.cxx:62
 AliGenMUONLMR.cxx:63
 AliGenMUONLMR.cxx:64
 AliGenMUONLMR.cxx:65
 AliGenMUONLMR.cxx:66
 AliGenMUONLMR.cxx:67
 AliGenMUONLMR.cxx:68
 AliGenMUONLMR.cxx:69
 AliGenMUONLMR.cxx:70
 AliGenMUONLMR.cxx:71
 AliGenMUONLMR.cxx:72
 AliGenMUONLMR.cxx:73
 AliGenMUONLMR.cxx:74
 AliGenMUONLMR.cxx:75
 AliGenMUONLMR.cxx:76
 AliGenMUONLMR.cxx:77
 AliGenMUONLMR.cxx:78
 AliGenMUONLMR.cxx:79
 AliGenMUONLMR.cxx:80
 AliGenMUONLMR.cxx:81
 AliGenMUONLMR.cxx:82
 AliGenMUONLMR.cxx:83
 AliGenMUONLMR.cxx:84
 AliGenMUONLMR.cxx:85
 AliGenMUONLMR.cxx:86
 AliGenMUONLMR.cxx:87
 AliGenMUONLMR.cxx:88
 AliGenMUONLMR.cxx:89
 AliGenMUONLMR.cxx:90
 AliGenMUONLMR.cxx:91
 AliGenMUONLMR.cxx:92
 AliGenMUONLMR.cxx:93
 AliGenMUONLMR.cxx:94
 AliGenMUONLMR.cxx:95
 AliGenMUONLMR.cxx:96
 AliGenMUONLMR.cxx:97
 AliGenMUONLMR.cxx:98
 AliGenMUONLMR.cxx:99
 AliGenMUONLMR.cxx:100
 AliGenMUONLMR.cxx:101
 AliGenMUONLMR.cxx:102
 AliGenMUONLMR.cxx:103
 AliGenMUONLMR.cxx:104
 AliGenMUONLMR.cxx:105
 AliGenMUONLMR.cxx:106
 AliGenMUONLMR.cxx:107
 AliGenMUONLMR.cxx:108
 AliGenMUONLMR.cxx:109
 AliGenMUONLMR.cxx:110
 AliGenMUONLMR.cxx:111
 AliGenMUONLMR.cxx:112
 AliGenMUONLMR.cxx:113
 AliGenMUONLMR.cxx:114
 AliGenMUONLMR.cxx:115
 AliGenMUONLMR.cxx:116
 AliGenMUONLMR.cxx:117
 AliGenMUONLMR.cxx:118
 AliGenMUONLMR.cxx:119
 AliGenMUONLMR.cxx:120
 AliGenMUONLMR.cxx:121
 AliGenMUONLMR.cxx:122
 AliGenMUONLMR.cxx:123
 AliGenMUONLMR.cxx:124
 AliGenMUONLMR.cxx:125
 AliGenMUONLMR.cxx:126
 AliGenMUONLMR.cxx:127
 AliGenMUONLMR.cxx:128
 AliGenMUONLMR.cxx:129
 AliGenMUONLMR.cxx:130
 AliGenMUONLMR.cxx:131
 AliGenMUONLMR.cxx:132
 AliGenMUONLMR.cxx:133
 AliGenMUONLMR.cxx:134
 AliGenMUONLMR.cxx:135
 AliGenMUONLMR.cxx:136
 AliGenMUONLMR.cxx:137
 AliGenMUONLMR.cxx:138
 AliGenMUONLMR.cxx:139
 AliGenMUONLMR.cxx:140
 AliGenMUONLMR.cxx:141
 AliGenMUONLMR.cxx:142
 AliGenMUONLMR.cxx:143
 AliGenMUONLMR.cxx:144
 AliGenMUONLMR.cxx:145
 AliGenMUONLMR.cxx:146
 AliGenMUONLMR.cxx:147
 AliGenMUONLMR.cxx:148
 AliGenMUONLMR.cxx:149
 AliGenMUONLMR.cxx:150
 AliGenMUONLMR.cxx:151
 AliGenMUONLMR.cxx:152
 AliGenMUONLMR.cxx:153
 AliGenMUONLMR.cxx:154
 AliGenMUONLMR.cxx:155
 AliGenMUONLMR.cxx:156
 AliGenMUONLMR.cxx:157
 AliGenMUONLMR.cxx:158
 AliGenMUONLMR.cxx:159
 AliGenMUONLMR.cxx:160
 AliGenMUONLMR.cxx:161
 AliGenMUONLMR.cxx:162
 AliGenMUONLMR.cxx:163
 AliGenMUONLMR.cxx:164
 AliGenMUONLMR.cxx:165
 AliGenMUONLMR.cxx:166
 AliGenMUONLMR.cxx:167
 AliGenMUONLMR.cxx:168
 AliGenMUONLMR.cxx:169
 AliGenMUONLMR.cxx:170
 AliGenMUONLMR.cxx:171
 AliGenMUONLMR.cxx:172
 AliGenMUONLMR.cxx:173
 AliGenMUONLMR.cxx:174
 AliGenMUONLMR.cxx:175
 AliGenMUONLMR.cxx:176
 AliGenMUONLMR.cxx:177
 AliGenMUONLMR.cxx:178
 AliGenMUONLMR.cxx:179
 AliGenMUONLMR.cxx:180
 AliGenMUONLMR.cxx:181
 AliGenMUONLMR.cxx:182
 AliGenMUONLMR.cxx:183
 AliGenMUONLMR.cxx:184
 AliGenMUONLMR.cxx:185
 AliGenMUONLMR.cxx:186
 AliGenMUONLMR.cxx:187
 AliGenMUONLMR.cxx:188
 AliGenMUONLMR.cxx:189
 AliGenMUONLMR.cxx:190
 AliGenMUONLMR.cxx:191
 AliGenMUONLMR.cxx:192
 AliGenMUONLMR.cxx:193
 AliGenMUONLMR.cxx:194
 AliGenMUONLMR.cxx:195
 AliGenMUONLMR.cxx:196
 AliGenMUONLMR.cxx:197
 AliGenMUONLMR.cxx:198
 AliGenMUONLMR.cxx:199
 AliGenMUONLMR.cxx:200
 AliGenMUONLMR.cxx:201
 AliGenMUONLMR.cxx:202
 AliGenMUONLMR.cxx:203
 AliGenMUONLMR.cxx:204
 AliGenMUONLMR.cxx:205
 AliGenMUONLMR.cxx:206
 AliGenMUONLMR.cxx:207
 AliGenMUONLMR.cxx:208
 AliGenMUONLMR.cxx:209
 AliGenMUONLMR.cxx:210
 AliGenMUONLMR.cxx:211
 AliGenMUONLMR.cxx:212
 AliGenMUONLMR.cxx:213
 AliGenMUONLMR.cxx:214
 AliGenMUONLMR.cxx:215
 AliGenMUONLMR.cxx:216
 AliGenMUONLMR.cxx:217
 AliGenMUONLMR.cxx:218
 AliGenMUONLMR.cxx:219
 AliGenMUONLMR.cxx:220
 AliGenMUONLMR.cxx:221
 AliGenMUONLMR.cxx:222
 AliGenMUONLMR.cxx:223
 AliGenMUONLMR.cxx:224
 AliGenMUONLMR.cxx:225
 AliGenMUONLMR.cxx:226
 AliGenMUONLMR.cxx:227
 AliGenMUONLMR.cxx:228
 AliGenMUONLMR.cxx:229
 AliGenMUONLMR.cxx:230
 AliGenMUONLMR.cxx:231
 AliGenMUONLMR.cxx:232
 AliGenMUONLMR.cxx:233
 AliGenMUONLMR.cxx:234
 AliGenMUONLMR.cxx:235
 AliGenMUONLMR.cxx:236
 AliGenMUONLMR.cxx:237
 AliGenMUONLMR.cxx:238
 AliGenMUONLMR.cxx:239
 AliGenMUONLMR.cxx:240
 AliGenMUONLMR.cxx:241
 AliGenMUONLMR.cxx:242
 AliGenMUONLMR.cxx:243
 AliGenMUONLMR.cxx:244
 AliGenMUONLMR.cxx:245
 AliGenMUONLMR.cxx:246
 AliGenMUONLMR.cxx:247
 AliGenMUONLMR.cxx:248
 AliGenMUONLMR.cxx:249
 AliGenMUONLMR.cxx:250
 AliGenMUONLMR.cxx:251
 AliGenMUONLMR.cxx:252
 AliGenMUONLMR.cxx:253
 AliGenMUONLMR.cxx:254
 AliGenMUONLMR.cxx:255
 AliGenMUONLMR.cxx:256
 AliGenMUONLMR.cxx:257
 AliGenMUONLMR.cxx:258
 AliGenMUONLMR.cxx:259
 AliGenMUONLMR.cxx:260
 AliGenMUONLMR.cxx:261
 AliGenMUONLMR.cxx:262
 AliGenMUONLMR.cxx:263
 AliGenMUONLMR.cxx:264
 AliGenMUONLMR.cxx:265
 AliGenMUONLMR.cxx:266
 AliGenMUONLMR.cxx:267
 AliGenMUONLMR.cxx:268
 AliGenMUONLMR.cxx:269
 AliGenMUONLMR.cxx:270
 AliGenMUONLMR.cxx:271
 AliGenMUONLMR.cxx:272
 AliGenMUONLMR.cxx:273
 AliGenMUONLMR.cxx:274
 AliGenMUONLMR.cxx:275
 AliGenMUONLMR.cxx:276
 AliGenMUONLMR.cxx:277
 AliGenMUONLMR.cxx:278
 AliGenMUONLMR.cxx:279
 AliGenMUONLMR.cxx:280
 AliGenMUONLMR.cxx:281
 AliGenMUONLMR.cxx:282
 AliGenMUONLMR.cxx:283
 AliGenMUONLMR.cxx:284
 AliGenMUONLMR.cxx:285
 AliGenMUONLMR.cxx:286
 AliGenMUONLMR.cxx:287
 AliGenMUONLMR.cxx:288
 AliGenMUONLMR.cxx:289
 AliGenMUONLMR.cxx:290
 AliGenMUONLMR.cxx:291
 AliGenMUONLMR.cxx:292
 AliGenMUONLMR.cxx:293
 AliGenMUONLMR.cxx:294
 AliGenMUONLMR.cxx:295
 AliGenMUONLMR.cxx:296
 AliGenMUONLMR.cxx:297
 AliGenMUONLMR.cxx:298
 AliGenMUONLMR.cxx:299
 AliGenMUONLMR.cxx:300
 AliGenMUONLMR.cxx:301
 AliGenMUONLMR.cxx:302
 AliGenMUONLMR.cxx:303
 AliGenMUONLMR.cxx:304
 AliGenMUONLMR.cxx:305
 AliGenMUONLMR.cxx:306
 AliGenMUONLMR.cxx:307
 AliGenMUONLMR.cxx:308
 AliGenMUONLMR.cxx:309
 AliGenMUONLMR.cxx:310
 AliGenMUONLMR.cxx:311
 AliGenMUONLMR.cxx:312
 AliGenMUONLMR.cxx:313
 AliGenMUONLMR.cxx:314
 AliGenMUONLMR.cxx:315
 AliGenMUONLMR.cxx:316
 AliGenMUONLMR.cxx:317
 AliGenMUONLMR.cxx:318
 AliGenMUONLMR.cxx:319
 AliGenMUONLMR.cxx:320
 AliGenMUONLMR.cxx:321
 AliGenMUONLMR.cxx:322
 AliGenMUONLMR.cxx:323
 AliGenMUONLMR.cxx:324
 AliGenMUONLMR.cxx:325
 AliGenMUONLMR.cxx:326
 AliGenMUONLMR.cxx:327
 AliGenMUONLMR.cxx:328
 AliGenMUONLMR.cxx:329
 AliGenMUONLMR.cxx:330
 AliGenMUONLMR.cxx:331
 AliGenMUONLMR.cxx:332
 AliGenMUONLMR.cxx:333
 AliGenMUONLMR.cxx:334
 AliGenMUONLMR.cxx:335
 AliGenMUONLMR.cxx:336
 AliGenMUONLMR.cxx:337
 AliGenMUONLMR.cxx:338
 AliGenMUONLMR.cxx:339
 AliGenMUONLMR.cxx:340
 AliGenMUONLMR.cxx:341
 AliGenMUONLMR.cxx:342
 AliGenMUONLMR.cxx:343
 AliGenMUONLMR.cxx:344
 AliGenMUONLMR.cxx:345
 AliGenMUONLMR.cxx:346
 AliGenMUONLMR.cxx:347
 AliGenMUONLMR.cxx:348
 AliGenMUONLMR.cxx:349
 AliGenMUONLMR.cxx:350
 AliGenMUONLMR.cxx:351
 AliGenMUONLMR.cxx:352
 AliGenMUONLMR.cxx:353
 AliGenMUONLMR.cxx:354
 AliGenMUONLMR.cxx:355
 AliGenMUONLMR.cxx:356
 AliGenMUONLMR.cxx:357
 AliGenMUONLMR.cxx:358
 AliGenMUONLMR.cxx:359
 AliGenMUONLMR.cxx:360
 AliGenMUONLMR.cxx:361
 AliGenMUONLMR.cxx:362
 AliGenMUONLMR.cxx:363
 AliGenMUONLMR.cxx:364
 AliGenMUONLMR.cxx:365
 AliGenMUONLMR.cxx:366
 AliGenMUONLMR.cxx:367
 AliGenMUONLMR.cxx:368
 AliGenMUONLMR.cxx:369
 AliGenMUONLMR.cxx:370
 AliGenMUONLMR.cxx:371
 AliGenMUONLMR.cxx:372
 AliGenMUONLMR.cxx:373
 AliGenMUONLMR.cxx:374
 AliGenMUONLMR.cxx:375
 AliGenMUONLMR.cxx:376
 AliGenMUONLMR.cxx:377
 AliGenMUONLMR.cxx:378
 AliGenMUONLMR.cxx:379
 AliGenMUONLMR.cxx:380
 AliGenMUONLMR.cxx:381
 AliGenMUONLMR.cxx:382
 AliGenMUONLMR.cxx:383
 AliGenMUONLMR.cxx:384
 AliGenMUONLMR.cxx:385
 AliGenMUONLMR.cxx:386
 AliGenMUONLMR.cxx:387
 AliGenMUONLMR.cxx:388
 AliGenMUONLMR.cxx:389
 AliGenMUONLMR.cxx:390
 AliGenMUONLMR.cxx:391
 AliGenMUONLMR.cxx:392
 AliGenMUONLMR.cxx:393
 AliGenMUONLMR.cxx:394
 AliGenMUONLMR.cxx:395
 AliGenMUONLMR.cxx:396
 AliGenMUONLMR.cxx:397
 AliGenMUONLMR.cxx:398
 AliGenMUONLMR.cxx:399
 AliGenMUONLMR.cxx:400
 AliGenMUONLMR.cxx:401
 AliGenMUONLMR.cxx:402
 AliGenMUONLMR.cxx:403
 AliGenMUONLMR.cxx:404
 AliGenMUONLMR.cxx:405
 AliGenMUONLMR.cxx:406
 AliGenMUONLMR.cxx:407
 AliGenMUONLMR.cxx:408
 AliGenMUONLMR.cxx:409
 AliGenMUONLMR.cxx:410
 AliGenMUONLMR.cxx:411
 AliGenMUONLMR.cxx:412
 AliGenMUONLMR.cxx:413
 AliGenMUONLMR.cxx:414
 AliGenMUONLMR.cxx:415
 AliGenMUONLMR.cxx:416
 AliGenMUONLMR.cxx:417
 AliGenMUONLMR.cxx:418
 AliGenMUONLMR.cxx:419
 AliGenMUONLMR.cxx:420
 AliGenMUONLMR.cxx:421
 AliGenMUONLMR.cxx:422
 AliGenMUONLMR.cxx:423
 AliGenMUONLMR.cxx:424
 AliGenMUONLMR.cxx:425
 AliGenMUONLMR.cxx:426
 AliGenMUONLMR.cxx:427
 AliGenMUONLMR.cxx:428
 AliGenMUONLMR.cxx:429
 AliGenMUONLMR.cxx:430
 AliGenMUONLMR.cxx:431
 AliGenMUONLMR.cxx:432
 AliGenMUONLMR.cxx:433
 AliGenMUONLMR.cxx:434
 AliGenMUONLMR.cxx:435
 AliGenMUONLMR.cxx:436
 AliGenMUONLMR.cxx:437
 AliGenMUONLMR.cxx:438
 AliGenMUONLMR.cxx:439
 AliGenMUONLMR.cxx:440
 AliGenMUONLMR.cxx:441
 AliGenMUONLMR.cxx:442
 AliGenMUONLMR.cxx:443
 AliGenMUONLMR.cxx:444
 AliGenMUONLMR.cxx:445
 AliGenMUONLMR.cxx:446
 AliGenMUONLMR.cxx:447
 AliGenMUONLMR.cxx:448
 AliGenMUONLMR.cxx:449
 AliGenMUONLMR.cxx:450
 AliGenMUONLMR.cxx:451
 AliGenMUONLMR.cxx:452
 AliGenMUONLMR.cxx:453
 AliGenMUONLMR.cxx:454
 AliGenMUONLMR.cxx:455
 AliGenMUONLMR.cxx:456
 AliGenMUONLMR.cxx:457
 AliGenMUONLMR.cxx:458
 AliGenMUONLMR.cxx:459
 AliGenMUONLMR.cxx:460
 AliGenMUONLMR.cxx:461
 AliGenMUONLMR.cxx:462
 AliGenMUONLMR.cxx:463
 AliGenMUONLMR.cxx:464
 AliGenMUONLMR.cxx:465
 AliGenMUONLMR.cxx:466
 AliGenMUONLMR.cxx:467
 AliGenMUONLMR.cxx:468
 AliGenMUONLMR.cxx:469
 AliGenMUONLMR.cxx:470
 AliGenMUONLMR.cxx:471
 AliGenMUONLMR.cxx:472
 AliGenMUONLMR.cxx:473
 AliGenMUONLMR.cxx:474
 AliGenMUONLMR.cxx:475
 AliGenMUONLMR.cxx:476
 AliGenMUONLMR.cxx:477
 AliGenMUONLMR.cxx:478
 AliGenMUONLMR.cxx:479
 AliGenMUONLMR.cxx:480
 AliGenMUONLMR.cxx:481
 AliGenMUONLMR.cxx:482
 AliGenMUONLMR.cxx:483
 AliGenMUONLMR.cxx:484
 AliGenMUONLMR.cxx:485
 AliGenMUONLMR.cxx:486
 AliGenMUONLMR.cxx:487
 AliGenMUONLMR.cxx:488
 AliGenMUONLMR.cxx:489
 AliGenMUONLMR.cxx:490
 AliGenMUONLMR.cxx:491
 AliGenMUONLMR.cxx:492
 AliGenMUONLMR.cxx:493
 AliGenMUONLMR.cxx:494
 AliGenMUONLMR.cxx:495
 AliGenMUONLMR.cxx:496
 AliGenMUONLMR.cxx:497
 AliGenMUONLMR.cxx:498
 AliGenMUONLMR.cxx:499
 AliGenMUONLMR.cxx:500
 AliGenMUONLMR.cxx:501
 AliGenMUONLMR.cxx:502
 AliGenMUONLMR.cxx:503
 AliGenMUONLMR.cxx:504
 AliGenMUONLMR.cxx:505
 AliGenMUONLMR.cxx:506
 AliGenMUONLMR.cxx:507
 AliGenMUONLMR.cxx:508
 AliGenMUONLMR.cxx:509
 AliGenMUONLMR.cxx:510
 AliGenMUONLMR.cxx:511
 AliGenMUONLMR.cxx:512
 AliGenMUONLMR.cxx:513
 AliGenMUONLMR.cxx:514
 AliGenMUONLMR.cxx:515
 AliGenMUONLMR.cxx:516
 AliGenMUONLMR.cxx:517
 AliGenMUONLMR.cxx:518
 AliGenMUONLMR.cxx:519
 AliGenMUONLMR.cxx:520
 AliGenMUONLMR.cxx:521
 AliGenMUONLMR.cxx:522
 AliGenMUONLMR.cxx:523
 AliGenMUONLMR.cxx:524
 AliGenMUONLMR.cxx:525
 AliGenMUONLMR.cxx:526
 AliGenMUONLMR.cxx:527
 AliGenMUONLMR.cxx:528
 AliGenMUONLMR.cxx:529
 AliGenMUONLMR.cxx:530
 AliGenMUONLMR.cxx:531
 AliGenMUONLMR.cxx:532
 AliGenMUONLMR.cxx:533
 AliGenMUONLMR.cxx:534
 AliGenMUONLMR.cxx:535
 AliGenMUONLMR.cxx:536
 AliGenMUONLMR.cxx:537
 AliGenMUONLMR.cxx:538
 AliGenMUONLMR.cxx:539
 AliGenMUONLMR.cxx:540
 AliGenMUONLMR.cxx:541
 AliGenMUONLMR.cxx:542
 AliGenMUONLMR.cxx:543
 AliGenMUONLMR.cxx:544
 AliGenMUONLMR.cxx:545
 AliGenMUONLMR.cxx:546
 AliGenMUONLMR.cxx:547
 AliGenMUONLMR.cxx:548
 AliGenMUONLMR.cxx:549
 AliGenMUONLMR.cxx:550
 AliGenMUONLMR.cxx:551
 AliGenMUONLMR.cxx:552
 AliGenMUONLMR.cxx:553
 AliGenMUONLMR.cxx:554
 AliGenMUONLMR.cxx:555
 AliGenMUONLMR.cxx:556
 AliGenMUONLMR.cxx:557
 AliGenMUONLMR.cxx:558
 AliGenMUONLMR.cxx:559
 AliGenMUONLMR.cxx:560
 AliGenMUONLMR.cxx:561
 AliGenMUONLMR.cxx:562
 AliGenMUONLMR.cxx:563
 AliGenMUONLMR.cxx:564
 AliGenMUONLMR.cxx:565
 AliGenMUONLMR.cxx:566
 AliGenMUONLMR.cxx:567
 AliGenMUONLMR.cxx:568
 AliGenMUONLMR.cxx:569
 AliGenMUONLMR.cxx:570
 AliGenMUONLMR.cxx:571
 AliGenMUONLMR.cxx:572
 AliGenMUONLMR.cxx:573
 AliGenMUONLMR.cxx:574
 AliGenMUONLMR.cxx:575
 AliGenMUONLMR.cxx:576
 AliGenMUONLMR.cxx:577
 AliGenMUONLMR.cxx:578
 AliGenMUONLMR.cxx:579
 AliGenMUONLMR.cxx:580
 AliGenMUONLMR.cxx:581
 AliGenMUONLMR.cxx:582
 AliGenMUONLMR.cxx:583
 AliGenMUONLMR.cxx:584
 AliGenMUONLMR.cxx:585
 AliGenMUONLMR.cxx:586
 AliGenMUONLMR.cxx:587
 AliGenMUONLMR.cxx:588
 AliGenMUONLMR.cxx:589
 AliGenMUONLMR.cxx:590
 AliGenMUONLMR.cxx:591
 AliGenMUONLMR.cxx:592
 AliGenMUONLMR.cxx:593
 AliGenMUONLMR.cxx:594
 AliGenMUONLMR.cxx:595
 AliGenMUONLMR.cxx:596
 AliGenMUONLMR.cxx:597
 AliGenMUONLMR.cxx:598
 AliGenMUONLMR.cxx:599
 AliGenMUONLMR.cxx:600
 AliGenMUONLMR.cxx:601
 AliGenMUONLMR.cxx:602
 AliGenMUONLMR.cxx:603
 AliGenMUONLMR.cxx:604
 AliGenMUONLMR.cxx:605
 AliGenMUONLMR.cxx:606
 AliGenMUONLMR.cxx:607
 AliGenMUONLMR.cxx:608
 AliGenMUONLMR.cxx:609
 AliGenMUONLMR.cxx:610
 AliGenMUONLMR.cxx:611
 AliGenMUONLMR.cxx:612
 AliGenMUONLMR.cxx:613
 AliGenMUONLMR.cxx:614
 AliGenMUONLMR.cxx:615
 AliGenMUONLMR.cxx:616
 AliGenMUONLMR.cxx:617
 AliGenMUONLMR.cxx:618
 AliGenMUONLMR.cxx:619
 AliGenMUONLMR.cxx:620
 AliGenMUONLMR.cxx:621
 AliGenMUONLMR.cxx:622
 AliGenMUONLMR.cxx:623
 AliGenMUONLMR.cxx:624
 AliGenMUONLMR.cxx:625
 AliGenMUONLMR.cxx:626
 AliGenMUONLMR.cxx:627
 AliGenMUONLMR.cxx:628
 AliGenMUONLMR.cxx:629
 AliGenMUONLMR.cxx:630
 AliGenMUONLMR.cxx:631
 AliGenMUONLMR.cxx:632
 AliGenMUONLMR.cxx:633
 AliGenMUONLMR.cxx:634
 AliGenMUONLMR.cxx:635
 AliGenMUONLMR.cxx:636
 AliGenMUONLMR.cxx:637
 AliGenMUONLMR.cxx:638
 AliGenMUONLMR.cxx:639
 AliGenMUONLMR.cxx:640
 AliGenMUONLMR.cxx:641
 AliGenMUONLMR.cxx:642
 AliGenMUONLMR.cxx:643
 AliGenMUONLMR.cxx:644
 AliGenMUONLMR.cxx:645
 AliGenMUONLMR.cxx:646
 AliGenMUONLMR.cxx:647
 AliGenMUONLMR.cxx:648
 AliGenMUONLMR.cxx:649
 AliGenMUONLMR.cxx:650
 AliGenMUONLMR.cxx:651
 AliGenMUONLMR.cxx:652
 AliGenMUONLMR.cxx:653
 AliGenMUONLMR.cxx:654
 AliGenMUONLMR.cxx:655
 AliGenMUONLMR.cxx:656
 AliGenMUONLMR.cxx:657
 AliGenMUONLMR.cxx:658
 AliGenMUONLMR.cxx:659
 AliGenMUONLMR.cxx:660
 AliGenMUONLMR.cxx:661
 AliGenMUONLMR.cxx:662
 AliGenMUONLMR.cxx:663
 AliGenMUONLMR.cxx:664
 AliGenMUONLMR.cxx:665
 AliGenMUONLMR.cxx:666
 AliGenMUONLMR.cxx:667
 AliGenMUONLMR.cxx:668
 AliGenMUONLMR.cxx:669
 AliGenMUONLMR.cxx:670
 AliGenMUONLMR.cxx:671
 AliGenMUONLMR.cxx:672
 AliGenMUONLMR.cxx:673
 AliGenMUONLMR.cxx:674
 AliGenMUONLMR.cxx:675
 AliGenMUONLMR.cxx:676
 AliGenMUONLMR.cxx:677
 AliGenMUONLMR.cxx:678
 AliGenMUONLMR.cxx:679
 AliGenMUONLMR.cxx:680
 AliGenMUONLMR.cxx:681
 AliGenMUONLMR.cxx:682
 AliGenMUONLMR.cxx:683
 AliGenMUONLMR.cxx:684
 AliGenMUONLMR.cxx:685
 AliGenMUONLMR.cxx:686
 AliGenMUONLMR.cxx:687
 AliGenMUONLMR.cxx:688
 AliGenMUONLMR.cxx:689
 AliGenMUONLMR.cxx:690
 AliGenMUONLMR.cxx:691
 AliGenMUONLMR.cxx:692
 AliGenMUONLMR.cxx:693
 AliGenMUONLMR.cxx:694
 AliGenMUONLMR.cxx:695
 AliGenMUONLMR.cxx:696
 AliGenMUONLMR.cxx:697
 AliGenMUONLMR.cxx:698
 AliGenMUONLMR.cxx:699
 AliGenMUONLMR.cxx:700
 AliGenMUONLMR.cxx:701
 AliGenMUONLMR.cxx:702
 AliGenMUONLMR.cxx:703
 AliGenMUONLMR.cxx:704
 AliGenMUONLMR.cxx:705
 AliGenMUONLMR.cxx:706
 AliGenMUONLMR.cxx:707
 AliGenMUONLMR.cxx:708
 AliGenMUONLMR.cxx:709
 AliGenMUONLMR.cxx:710
 AliGenMUONLMR.cxx:711
 AliGenMUONLMR.cxx:712
 AliGenMUONLMR.cxx:713
 AliGenMUONLMR.cxx:714
 AliGenMUONLMR.cxx:715
 AliGenMUONLMR.cxx:716
 AliGenMUONLMR.cxx:717
 AliGenMUONLMR.cxx:718
 AliGenMUONLMR.cxx:719
 AliGenMUONLMR.cxx:720
 AliGenMUONLMR.cxx:721
 AliGenMUONLMR.cxx:722
 AliGenMUONLMR.cxx:723
 AliGenMUONLMR.cxx:724
 AliGenMUONLMR.cxx:725
 AliGenMUONLMR.cxx:726
 AliGenMUONLMR.cxx:727
 AliGenMUONLMR.cxx:728
 AliGenMUONLMR.cxx:729
 AliGenMUONLMR.cxx:730
 AliGenMUONLMR.cxx:731
 AliGenMUONLMR.cxx:732
 AliGenMUONLMR.cxx:733
 AliGenMUONLMR.cxx:734
 AliGenMUONLMR.cxx:735
 AliGenMUONLMR.cxx:736
 AliGenMUONLMR.cxx:737
 AliGenMUONLMR.cxx:738
 AliGenMUONLMR.cxx:739
 AliGenMUONLMR.cxx:740
 AliGenMUONLMR.cxx:741
 AliGenMUONLMR.cxx:742
 AliGenMUONLMR.cxx:743
 AliGenMUONLMR.cxx:744
 AliGenMUONLMR.cxx:745
 AliGenMUONLMR.cxx:746
 AliGenMUONLMR.cxx:747
 AliGenMUONLMR.cxx:748
 AliGenMUONLMR.cxx:749
 AliGenMUONLMR.cxx:750
 AliGenMUONLMR.cxx:751
 AliGenMUONLMR.cxx:752
 AliGenMUONLMR.cxx:753
 AliGenMUONLMR.cxx:754
 AliGenMUONLMR.cxx:755
 AliGenMUONLMR.cxx:756
 AliGenMUONLMR.cxx:757
 AliGenMUONLMR.cxx:758
 AliGenMUONLMR.cxx:759
 AliGenMUONLMR.cxx:760
 AliGenMUONLMR.cxx:761
 AliGenMUONLMR.cxx:762
 AliGenMUONLMR.cxx:763
 AliGenMUONLMR.cxx:764
 AliGenMUONLMR.cxx:765
 AliGenMUONLMR.cxx:766
 AliGenMUONLMR.cxx:767
 AliGenMUONLMR.cxx:768
 AliGenMUONLMR.cxx:769
 AliGenMUONLMR.cxx:770
 AliGenMUONLMR.cxx:771
 AliGenMUONLMR.cxx:772
 AliGenMUONLMR.cxx:773
 AliGenMUONLMR.cxx:774
 AliGenMUONLMR.cxx:775
 AliGenMUONLMR.cxx:776
 AliGenMUONLMR.cxx:777
 AliGenMUONLMR.cxx:778
 AliGenMUONLMR.cxx:779
 AliGenMUONLMR.cxx:780
 AliGenMUONLMR.cxx:781
 AliGenMUONLMR.cxx:782
 AliGenMUONLMR.cxx:783
 AliGenMUONLMR.cxx:784
 AliGenMUONLMR.cxx:785
 AliGenMUONLMR.cxx:786
 AliGenMUONLMR.cxx:787
 AliGenMUONLMR.cxx:788
 AliGenMUONLMR.cxx:789
 AliGenMUONLMR.cxx:790
 AliGenMUONLMR.cxx:791
 AliGenMUONLMR.cxx:792
 AliGenMUONLMR.cxx:793
 AliGenMUONLMR.cxx:794
 AliGenMUONLMR.cxx:795
 AliGenMUONLMR.cxx:796
 AliGenMUONLMR.cxx:797
 AliGenMUONLMR.cxx:798
 AliGenMUONLMR.cxx:799
 AliGenMUONLMR.cxx:800
 AliGenMUONLMR.cxx:801
 AliGenMUONLMR.cxx:802
 AliGenMUONLMR.cxx:803
 AliGenMUONLMR.cxx:804
 AliGenMUONLMR.cxx:805
 AliGenMUONLMR.cxx:806
 AliGenMUONLMR.cxx:807
 AliGenMUONLMR.cxx:808
 AliGenMUONLMR.cxx:809
 AliGenMUONLMR.cxx:810
 AliGenMUONLMR.cxx:811
 AliGenMUONLMR.cxx:812
 AliGenMUONLMR.cxx:813
 AliGenMUONLMR.cxx:814
 AliGenMUONLMR.cxx:815
 AliGenMUONLMR.cxx:816
 AliGenMUONLMR.cxx:817
 AliGenMUONLMR.cxx:818
 AliGenMUONLMR.cxx:819
 AliGenMUONLMR.cxx:820
 AliGenMUONLMR.cxx:821
 AliGenMUONLMR.cxx:822
 AliGenMUONLMR.cxx:823
 AliGenMUONLMR.cxx:824
 AliGenMUONLMR.cxx:825
 AliGenMUONLMR.cxx:826
 AliGenMUONLMR.cxx:827
 AliGenMUONLMR.cxx:828
 AliGenMUONLMR.cxx:829
 AliGenMUONLMR.cxx:830
 AliGenMUONLMR.cxx:831
 AliGenMUONLMR.cxx:832
 AliGenMUONLMR.cxx:833
 AliGenMUONLMR.cxx:834
 AliGenMUONLMR.cxx:835
 AliGenMUONLMR.cxx:836
 AliGenMUONLMR.cxx:837
 AliGenMUONLMR.cxx:838
 AliGenMUONLMR.cxx:839
 AliGenMUONLMR.cxx:840
 AliGenMUONLMR.cxx:841
 AliGenMUONLMR.cxx:842
 AliGenMUONLMR.cxx:843
 AliGenMUONLMR.cxx:844
 AliGenMUONLMR.cxx:845
 AliGenMUONLMR.cxx:846
 AliGenMUONLMR.cxx:847
 AliGenMUONLMR.cxx:848
 AliGenMUONLMR.cxx:849
 AliGenMUONLMR.cxx:850
 AliGenMUONLMR.cxx:851
 AliGenMUONLMR.cxx:852
 AliGenMUONLMR.cxx:853
 AliGenMUONLMR.cxx:854
 AliGenMUONLMR.cxx:855
 AliGenMUONLMR.cxx:856
 AliGenMUONLMR.cxx:857
 AliGenMUONLMR.cxx:858
 AliGenMUONLMR.cxx:859
 AliGenMUONLMR.cxx:860
 AliGenMUONLMR.cxx:861
 AliGenMUONLMR.cxx:862
 AliGenMUONLMR.cxx:863
 AliGenMUONLMR.cxx:864
 AliGenMUONLMR.cxx:865
 AliGenMUONLMR.cxx:866
 AliGenMUONLMR.cxx:867
 AliGenMUONLMR.cxx:868
 AliGenMUONLMR.cxx:869
 AliGenMUONLMR.cxx:870
 AliGenMUONLMR.cxx:871
 AliGenMUONLMR.cxx:872
 AliGenMUONLMR.cxx:873
 AliGenMUONLMR.cxx:874
 AliGenMUONLMR.cxx:875
 AliGenMUONLMR.cxx:876
 AliGenMUONLMR.cxx:877
 AliGenMUONLMR.cxx:878
 AliGenMUONLMR.cxx:879
 AliGenMUONLMR.cxx:880
 AliGenMUONLMR.cxx:881
 AliGenMUONLMR.cxx:882
 AliGenMUONLMR.cxx:883
 AliGenMUONLMR.cxx:884
 AliGenMUONLMR.cxx:885
 AliGenMUONLMR.cxx:886