ROOT logo
/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id:$ */

#include "TChain.h"
#include "TList.h"
#include "TFile.h"
#include "TTree.h"
#include "TH1D.h"
#include "TH2D.h"
#include "TH3D.h"
#include "THnSparse.h"
#include "TCanvas.h"

#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include <TList.h>
#include <TMath.h>
#include <TTree.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TH3F.h>
#include <TProfile.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TH3D.h>
#include <TRandom.h>
#include "AliAnalysisManager.h"

#include "AliAODHandler.h"
#include "AliAODInputHandler.h"
#include "AliInputEventHandler.h"
#include "AliLog.h"
#include "AliESDEvent.h"
#include "AliESDInputHandler.h"
#include "AliMultiplicity.h"
#include "AliCentrality.h"
#include "AliDptDptInMC.h"

#include "AliPID.h"
#include "AliPIDResponse.h"

#include "AliESDVertex.h"
#include "AliESDEvent.h"
#include "AliESDInputHandler.h"
#include "AliAODEvent.h"
#include "AliAODTrack.h"
#include "AliAODInputHandler.h"
#include "AliESD.h"
#include "AliESDEvent.h"
#include "AliAODEvent.h"
#include "AliStack.h"
#include "AliESDtrackCuts.h"
#include "AliAODMCHeader.h"
#include "AliAODMCParticle.h"
#include "AliAODMCParticle.h"
#include "AliAODMCHeader.h"
#include "AliAODHeader.h"

#include "AliGenHijingEventHeader.h"
#include "AliGenEventHeader.h"
#include "AliPID.h"
#include "AliAODPid.h"
#include "AliPIDResponse.h"
#include "AliAODpidUtil.h"
#include "AliPIDCombined.h"

using std::cout;
using std::endl;

ClassImp(AliDptDptInMC)

AliDptDptInMC::AliDptDptInMC()
: AliAnalysisTaskSE(),
fAODEvent(0), 
fESDEvent(0),             //! ESD Event
fInputHandler(0),
fPIDResponse(0x0),
_outputHistoList(0),
fArrayMC (0),
_twoPi         ( 2.0 * 3.1415927),
_eventCount    ( 0), 
_debugLevel    ( 0),
_singlesOnly   ( 0), 
_useWeights    ( 0), 
_sameFilter    ( false),
_rejectPileup  ( 1), 
_rejectPairConversion ( 0), 
_vertexZMin           ( -10), 
_vertexZMax           (  10), 
_vertexXYMin          ( -10),
_vertexXYMax          (  10),
_centralityMethod     (  4),
_centralityMin        (  0.),
_centralityMax        (  0.),
_requestedCharge_1    (   1),
_requestedCharge_2    (  -1),
_dcaZMin              ( -3), 
_dcaZMax              (  3.), 
_dcaXYMin             ( -2.4), 
_dcaXYMax             (  2.4),
_dedxMin              ( 0),
_dedxMax              ( 100000),
_nClusterMin          ( 80), 
_trackFilterBit       (0),
fNSigmaCut            (3.),
fAnalysisType ("MCAOD"),
fExcludeResonancesInMC (kFALSE),
fExcludeElectronsInMC (kFALSE), 
_field    ( 1.),
_nTracks  ( 0 ),
_mult0    ( 0 ),
_mult1    ( 0 ),
_mult2    ( 0 ),
_mult3    ( 0 ),
_mult4    ( 0 ),
_mult4a    ( 0 ),
_mult5    ( 0 ),
_mult6    ( 0 ),
arraySize ( 2000),
_id_1(0),       
_charge_1(0),    
_iEtaPhi_1(0),    
_iPt_1(0),     
_pt_1(0),       
_px_1(0),      
_py_1(0),      
_pz_1(0),      
_correction_1(0),   
_dedx_1(0),   
_id_2(0),      
_charge_2(0),    
_iEtaPhi_2(0),    
_iPt_2(0),     
_pt_2(0),      
_px_2(0),      
_py_2(0),      
_pz_2(0),      
_correction_2(0),   
_dedx_2(0),   
_correctionWeight_1(0),   
_correctionWeight_2(0),
_nBins_M0(500),       _min_M0(0),        _max_M0(10000),          _width_M0(20),
_nBins_M1(500),       _min_M1(0),        _max_M1(10000),          _width_M1(20),
_nBins_M2(500),       _min_M2(0),        _max_M2(10000),          _width_M2(20),
_nBins_M3(500),       _min_M3(0),        _max_M3(10000),          _width_M3(20),
_nBins_M4(100),       _min_M4(0),        _max_M4(1),              _width_M4(0.01),
_nBins_M5(100),       _min_M5(0),        _max_M5(1),              _width_M5(0.01),
_nBins_M6(100),       _min_M6(0),        _max_M6(1),              _width_M6(0.01),
_nBins_vertexZ(40),   _min_vertexZ(-10), _max_vertexZ(10),        _width_vertexZ(0.5),

_nBins_pt_1(18),      _min_pt_1(0.2),    _max_pt_1(2.0),          _width_pt_1(0.1),
_nBins_phi_1(72),     _min_phi_1(0),     _max_phi_1(2.*3.1415927),_width_phi_1(2.*3.1415927/72.),
//_nBins_eta_1(18),    
//_min_eta_1(-0.9),  _max_eta_1(0.9),           _width_eta_1(0.1),
  _nBins_eta_1(0),    
  _min_eta_1(0),  _max_eta_1(0),           _width_eta_1(0.1),

_nBins_etaPhi_1(0), 
_nBins_etaPhiPt_1(0),
_nBins_zEtaPhiPt_1(0),

_nBins_pt_2(18),     _min_pt_2(0.2),     _max_pt_2(2.0),          _width_pt_2(0.1),
_nBins_phi_2(72),    _min_phi_2(0),      _max_phi_2(2.*3.1415927),_width_phi_2(2.*3.1415927/72),
//_nBins_eta_2(18),    
//_min_eta_2(-0.9),     _max_eta_2(0.9),           _width_eta_2(0.1),
  _nBins_eta_2(0),    
  _min_eta_2(0),     _max_eta_2(0),           _width_eta_2(0.1),



_nBins_etaPhi_2(0), 
_nBins_etaPhiPt_2(0),
_nBins_zEtaPhiPt_2(0),
_nBins_etaPhi_12(0),
__n1_1(0),
__n1_2(0),
__n2_12(0),   
__s1pt_1(0),
__s1pt_2(0),
__s2ptpt_12(0),
__s2NPt_12(0),
__s2PtN_12(0),
__n1Nw_1(0),
__n1Nw_2(0),
__n2Nw_12(0),   
__s1ptNw_1(0),
__s1ptNw_2(0),
__s2ptptNw_12(0),
__s2NPtNw_12(0),
__s2PtNNw_12(0),
__n1_1_vsPt(0),
__n1_1_vsEtaPhi(0), 
__s1pt_1_vsEtaPhi(0),
__n1_1_vsZEtaPhiPt(0),
__n1_2_vsPt(0),
__n1_2_vsEtaPhi(0), 
__s1pt_2_vsEtaPhi(0),
__n1_2_vsZEtaPhiPt(0),
__n2_12_vsPtPt(0),
__n2_12_vsEtaPhi(0),
__s2ptpt_12_vsEtaPhi(0),
__s2PtN_12_vsEtaPhi(0),
__s2NPt_12_vsEtaPhi(0),
_weight_1      ( 0    ),
_weight_2      ( 0    ),
_eventAccounting ( 0),
_m0 ( 0),
_m1 ( 0),
_m2 ( 0),
_m3 ( 0),
_m4 ( 0),
_m5 ( 0),
_m6 ( 0),
_vertexZ ( 0),
  
_Ncluster1  ( 0),
_Ncluster2  ( 0),
_etadis ( 0),
_phidis ( 0),
_dcaz   ( 0),
_dcaxy  ( 0),
_etadis2 ( 0),
_phidis2 ( 0),
_dcaz2   ( 0),
_dcaxy2  ( 0),
_etadis3 ( 0),
_phidis3 ( 0),
_dcaz3   ( 0),
_dcaxy3  ( 0),
_n1_1_vsPt         ( 0),         
_n1_1_vsEtaVsPhi   ( 0),
_s1pt_1_vsEtaVsPhi ( 0), 
_n1_1_vsZVsEtaVsPhiVsPt ( 0),
_n1_1_vsM          ( 0),  // w/ weight
_s1pt_1_vsM        ( 0),
_n1Nw_1_vsM        ( 0), // w/o weight
_s1ptNw_1_vsM      ( 0),
_dedxVsP_1         ( 0),
_corrDedxVsP_1     ( 0),
_betaVsP_1         ( 0),
_n1_2_vsPt         ( 0),       
_n1_2_vsEtaVsPhi   ( 0),
_s1pt_2_vsEtaVsPhi ( 0),
_n1_2_vsZVsEtaVsPhiVsPt ( 0), 
_n1_2_vsM          ( 0),
_s1pt_2_vsM        ( 0),
_n1Nw_2_vsM        ( 0),
_s1ptNw_2_vsM      ( 0),
_dedxVsP_2         ( 0),
_corrDedxVsP_2     ( 0),
_betaVsP_2         ( 0),
_n2_12_vsEtaPhi    ( 0),  
_n2_12_vsPtVsPt    ( 0),
_s2PtPt_12_vsEtaPhi( 0),    
_s2PtN_12_vsEtaPhi ( 0),       
_s2NPt_12_vsEtaPhi ( 0),     
_n2_12_vsM         ( 0),        
_s2PtPt_12_vsM     ( 0),    
_s2PtN_12_vsM      ( 0),       
_s2NPt_12_vsM      ( 0), 
_n2Nw_12_vsM       ( 0),        
_s2PtPtNw_12_vsM   ( 0),    
_s2PtNNw_12_vsM    ( 0),       
_s2NPtNw_12_vsM    ( 0), 
_invMass           ( 0),
_invMassElec       ( 0),
n1Name("NA"),
n1NwName("NA"),
n2Name("NA"),
n2NwName("NA"),
n3Name("NA"),
n1n1Name("NA"),
n1n1n1Name("NA"),
n2n1Name("NA"),
r1Name("NA"),
r2Name("NA"),
r3Name("NA"),
r2r1Name("NA"),
c2Name("NA"),
c3Name("NA"),
d3Name("NA"),
p3Name("NA"),
cName("NA"),

intR2Name("NA"),
binCorrName("NA"),
intBinCorrName("NA"),

countsName("NA"),
part_1_Name("NA"),
part_2_Name("NA"),
part_3_Name("NA"),
pair_12_Name("NA"),
pair_13_Name("NA"),
pair_23_Name("NA"),
tripletName("NA"),

avg("NA"),
avgName("NA"),
sumName("NA"),
s1ptName("NA"),
s1ptNwName("NA"),
s1DptName("NA"),

s2PtPtName("NA"),
s2NPtName("NA"),
s2PtNName("NA"),
s2DptDptName("NA"),

s2PtPtNwName("NA"),
s2NPtNwName("NA"),
s2PtNNwName("NA"),

ptName("NA"),
ptptName("NA"),
pt1pt1Name("NA"),
DptName("NA"),
DptDptName("NA"),
RDptDptName("NA"),
nPtName("NA"),
ptNName("NA"),
seanName("NA"),

_title_counts("NA"),

_title_m0("NA"),
_title_m1("NA"),
_title_m2("NA"),
_title_m3("NA"),
_title_m4("NA"),
_title_m5("NA"),
_title_m6("NA"),

_title_eta_1("NA"),
_title_phi_1("NA"),
_title_pt_1("NA"),
_title_etaPhi_1("NA"),
_title_n_1("NA"),
_title_SumPt_1("NA"),
_title_AvgPt_1("NA"),
_title_AvgN_1("NA"),
_title_AvgSumPt_1("NA"),

_title_eta_2("NA"),
_title_phi_2("NA"),
_title_pt_2("NA"),
_title_etaPhi_2("NA"),
_title_n_2("NA"),
_title_SumPt_2("NA"),
_title_AvgPt_2("NA"),
_title_AvgN_2("NA"),
_title_AvgSumPt_2("NA"),

_title_etaPhi_12("NA"),

_title_AvgN2_12("NA"),
_title_AvgSumPtPt_12("NA"),
_title_AvgSumPtN_12("NA"),
_title_AvgNSumPt_12("NA"),

vsZ("NA"),
vsM("NA"),
vsPt("NA"),
vsPhi("NA"), 
vsEta("NA"), 
vsEtaPhi("NA"), 
vsPtVsPt("NA")
{
  printf("Default constructor called \n");
  
  printf("passed \n ");
  
}

AliDptDptInMC::AliDptDptInMC(const TString & name)
: AliAnalysisTaskSE(name),
fAODEvent(0), 
fESDEvent(0),  
fInputHandler(0),
fPIDResponse(0),
_outputHistoList(0),
fArrayMC (0),
_twoPi         ( 2.0 * 3.1415927),
_eventCount    ( 0), 
_debugLevel    ( 0),
_singlesOnly   ( 0), 
_useWeights    ( 0), 
_sameFilter    ( false),
_rejectPileup  ( 1), 
_rejectPairConversion ( 0), 
_vertexZMin           ( -10.), 
_vertexZMax           (  10.), 
_vertexXYMin          ( -10.),
_vertexXYMax          (  10.),
_centralityMethod     (  4),
_centralityMin        (  0.),
_centralityMax        (  1.),
_requestedCharge_1    (   1),
_requestedCharge_2    (  -1),
_dcaZMin              ( -3), 
_dcaZMax              (  3.), 
_dcaXYMin             ( -2.4), 
_dcaXYMax             (  2.4),
_dedxMin              ( 0),
_dedxMax              ( 100000),
_nClusterMin          ( 80), 
_trackFilterBit       ( 0),
fNSigmaCut            ( 3.),

fAnalysisType ("MCAOD"),
fExcludeResonancesInMC (kFALSE),
fExcludeElectronsInMC (kFALSE), 

_field    ( 1.),
_nTracks  ( 0 ),
_mult0    ( 0 ),
_mult1    ( 0 ),
_mult2    ( 0 ),
_mult3    ( 0 ),
_mult4    ( 0 ),
_mult4a    ( 0 ),
_mult5    ( 0 ),
_mult6    ( 0 ),
arraySize ( 2000),
_id_1(0),       
_charge_1(0),    
_iEtaPhi_1(0),    
_iPt_1(0),     
_pt_1(0),       
_px_1(0),      
_py_1(0),      
_pz_1(0),      
_correction_1(0),   
_dedx_1(0),   
_id_2(0),      
_charge_2(0),    
_iEtaPhi_2(0),    
_iPt_2(0),     
_pt_2(0),      
_px_2(0),      
_py_2(0),      
_pz_2(0),      
_correction_2(0),   
_dedx_2(0),   
_correctionWeight_1(0),   
_correctionWeight_2(0),
_nBins_M0(500),       _min_M0(0),        _max_M0(10000),          _width_M0(20),
_nBins_M1(500),       _min_M1(0),        _max_M1(10000),          _width_M1(20),
_nBins_M2(500),       _min_M2(0),        _max_M2(10000),          _width_M2(20),
_nBins_M3(500),       _min_M3(0),        _max_M3(10000),          _width_M3(20),
_nBins_M4(100),       _min_M4(0),        _max_M4(1),              _width_M4(0.01),
_nBins_M5(100),       _min_M5(0),        _max_M5(1),              _width_M5(0.01),
_nBins_M6(100),       _min_M6(0),        _max_M6(1),              _width_M6(0.01),
_nBins_vertexZ(40),   _min_vertexZ(-10), _max_vertexZ(10),        _width_vertexZ(0.5),

_nBins_pt_1(18),      _min_pt_1(0.2),    _max_pt_1(2.0),          _width_pt_1(0.1),
_nBins_phi_1(72),     _min_phi_1(0),     _max_phi_1(2.*3.1415927),_width_phi_1(2.*3.1415927/72.),
  //_nBins_eta_1(18),     
  //_min_eta_1(-0.9),    _max_eta_1(0.9),           _width_eta_1(0.1),
  _nBins_eta_1(0),     
  _min_eta_1(0),    _max_eta_1(0),           _width_eta_1(0.1),

_nBins_etaPhi_1(0), 
_nBins_etaPhiPt_1(0),
_nBins_zEtaPhiPt_1(0),

_nBins_pt_2(18),     _min_pt_2(0.2),     _max_pt_2(2.0),          _width_pt_2(0.1),
_nBins_phi_2(72),    _min_phi_2(0),      _max_phi_2(2.*3.1415927),_width_phi_2(2.*3.1415927/72),
  //_nBins_eta_2(18),   
  //_min_eta_2(-0.9),     _max_eta_2(0.9),           _width_eta_2(0.1),
  _nBins_eta_2(0),   
  _min_eta_2(0),     _max_eta_2(0),           _width_eta_2(0.1),


_nBins_etaPhi_2(0), 
_nBins_etaPhiPt_2(0),
_nBins_zEtaPhiPt_2(0),
_nBins_etaPhi_12(0),
__n1_1(0),
__n1_2(0),
__n2_12(0),   
__s1pt_1(0),
__s1pt_2(0),
__s2ptpt_12(0),
__s2NPt_12(0),
__s2PtN_12(0),
__n1Nw_1(0),
__n1Nw_2(0),
__n2Nw_12(0),   
__s1ptNw_1(0),
__s1ptNw_2(0),
__s2ptptNw_12(0),
__s2NPtNw_12(0),
__s2PtNNw_12(0),
__n1_1_vsPt(0),
__n1_1_vsEtaPhi(0), 
__s1pt_1_vsEtaPhi(0),
__n1_1_vsZEtaPhiPt(0),
__n1_2_vsPt(0),
__n1_2_vsEtaPhi(0), 
__s1pt_2_vsEtaPhi(0),
__n1_2_vsZEtaPhiPt(0),
__n2_12_vsPtPt(0),
__n2_12_vsEtaPhi(0),
__s2ptpt_12_vsEtaPhi(0),
__s2PtN_12_vsEtaPhi(0),
__s2NPt_12_vsEtaPhi(0),
_weight_1        ( 0    ),
_weight_2        ( 0    ),
_eventAccounting ( 0),
_m0 ( 0),
_m1 ( 0),
_m2 ( 0),
_m3 ( 0),
_m4 ( 0),
_m5 ( 0),
_m6 ( 0),
_vertexZ ( 0),
_Ncluster1  ( 0),
_Ncluster2  ( 0),
_etadis ( 0),
_phidis ( 0),

_dcaz ( 0),
_dcaxy ( 0),
_etadis2 ( 0),
_phidis2 ( 0),
_dcaz2   ( 0),
_dcaxy2  ( 0),
_etadis3 ( 0),
_phidis3 ( 0),
_dcaz3   ( 0),
_dcaxy3  ( 0),
_n1_1_vsPt         ( 0),         
_n1_1_vsEtaVsPhi   ( 0),
_s1pt_1_vsEtaVsPhi ( 0), 
_n1_1_vsZVsEtaVsPhiVsPt ( 0),
_n1_1_vsM          ( 0),  // w/ weight
_s1pt_1_vsM        ( 0),
_n1Nw_1_vsM        ( 0), // w/o weight
_s1ptNw_1_vsM      ( 0),
_dedxVsP_1         ( 0),
_corrDedxVsP_1     ( 0),
_betaVsP_1         ( 0),
_n1_2_vsPt         ( 0),       
_n1_2_vsEtaVsPhi   ( 0),
_s1pt_2_vsEtaVsPhi ( 0),
_n1_2_vsZVsEtaVsPhiVsPt ( 0), 
_n1_2_vsM          ( 0),
_s1pt_2_vsM        ( 0),
_n1Nw_2_vsM        ( 0),
_s1ptNw_2_vsM      ( 0),
_dedxVsP_2         ( 0),
_corrDedxVsP_2     ( 0),
_betaVsP_2         ( 0),
_n2_12_vsEtaPhi    ( 0),  
_n2_12_vsPtVsPt    ( 0),
_s2PtPt_12_vsEtaPhi( 0),    
_s2PtN_12_vsEtaPhi ( 0),       
_s2NPt_12_vsEtaPhi ( 0),     
_n2_12_vsM         ( 0),        
_s2PtPt_12_vsM     ( 0),    
_s2PtN_12_vsM      ( 0),       
_s2NPt_12_vsM      ( 0), 
_n2Nw_12_vsM       ( 0),        
_s2PtPtNw_12_vsM   ( 0),    
_s2PtNNw_12_vsM    ( 0),       
_s2NPtNw_12_vsM    ( 0), 
_invMass           ( 0),
_invMassElec       ( 0),
n1Name("NA"),
n1NwName("NA"),
n2Name("NA"),
n2NwName("NA"),
n3Name("NA"),
n1n1Name("NA"),
n1n1n1Name("NA"),
n2n1Name("NA"),
r1Name("NA"),
r2Name("NA"),
r3Name("NA"),
r2r1Name("NA"),
c2Name("NA"),
c3Name("NA"),
d3Name("NA"),
p3Name("NA"),
cName("NA"),

intR2Name("NA"),
binCorrName("NA"),
intBinCorrName("NA"),

countsName("NA"),
part_1_Name("NA"),
part_2_Name("NA"),
part_3_Name("NA"),
pair_12_Name("NA"),
pair_13_Name("NA"),
pair_23_Name("NA"),
tripletName("NA"),

avg("NA"),
avgName("NA"),
sumName("NA"),
s1ptName("NA"),
s1ptNwName("NA"),
s1DptName("NA"),

s2PtPtName("NA"),
s2NPtName("NA"),
s2PtNName("NA"),
s2DptDptName("NA"),

s2PtPtNwName("NA"),
s2NPtNwName("NA"),
s2PtNNwName("NA"),

ptName("NA"),
ptptName("NA"),
pt1pt1Name("NA"),
DptName("NA"),
DptDptName("NA"),
RDptDptName("NA"),
nPtName("NA"),
ptNName("NA"),
seanName("NA"),

_title_counts("NA"),

_title_m0("NA"),
_title_m1("NA"),
_title_m2("NA"),
_title_m3("NA"),
_title_m4("NA"),
_title_m5("NA"),
_title_m6("NA"),

_title_eta_1("NA"),
_title_phi_1("NA"),
_title_pt_1("NA"),
_title_etaPhi_1("NA"),
_title_n_1("NA"),
_title_SumPt_1("NA"),
_title_AvgPt_1("NA"),
_title_AvgN_1("NA"),
_title_AvgSumPt_1("NA"),

_title_eta_2("NA"),
_title_phi_2("NA"),
_title_pt_2("NA"),
_title_etaPhi_2("NA"),
_title_n_2("NA"),
_title_SumPt_2("NA"),
_title_AvgPt_2("NA"),
_title_AvgN_2("NA"),
_title_AvgSumPt_2("NA"),

_title_etaPhi_12("NA"),

_title_AvgN2_12("NA"),
_title_AvgSumPtPt_12("NA"),
_title_AvgSumPtN_12("NA"),
_title_AvgNSumPt_12("NA"),

vsZ("NA"),
vsM("NA"),
vsPt("NA"),
vsPhi("NA"), 
vsEta("NA"), 
vsEtaPhi("NA"), 
vsPtVsPt("NA")
{
  printf("2nd constructor called ");
  
  DefineOutput(0, TList::Class());
  
  printf("passed  ");
  
}

AliDptDptInMC::~AliDptDptInMC()
{
  /*
   delete _id_1;
   delete _charge_1; 
   delete _iEtaPhi_1; 
   delete _iPt_1;  
   delete _pt_1;    
   delete _px_1;   
   delete _py_1;   
   delete _pz_1;   
   delete _correction_1;
   delete _dedx_1;
   delete __n1_1_vsPt;
   delete __n1_1_vsEtaPhi;
   delete __s1pt_1_vsEtaPhi;
   delete __n1_1_vsZEtaPhiPt; 
   if (_correctionWeight_1) delete _correctionWeight_1;
   
   if (!_sameFilter)
   {
   delete _id_2;
   delete _charge_2; 
   delete _iEtaPhi_2; 
   delete _iPt_2;  
   delete _pt_2;    
   delete _px_2;   
   delete _py_2;   
   delete _pz_2;   
   delete _correction_2;
   delete _dedx_2;
   delete __n1_2_vsPt;
   delete __n1_2_vsEtaPhi;
   delete __s1pt_2_vsEtaPhi;
   delete __n1_2_vsZEtaPhiPt; 
   if (_correctionWeight_2) delete _correctionWeight_2;
   }
   
   if (!_singlesOnly)
   {
   delete __n2_12_vsPtPt;
   delete __n2_12_vsEtaPhi;
   delete __s2ptpt_12_vsEtaPhi;
   delete __s2PtN_12_vsEtaPhi;
   delete __s2NPt_12_vsEtaPhi;
   }
   */
}

void AliDptDptInMC::UserCreateOutputObjects()
{
  //cout<< "AliDptDptInMC::CreateOutputObjects() Starting " << endl;
  OpenFile(0);
  _outputHistoList = new TList();
  _outputHistoList->SetOwner();
  
  //if (_useWeights) DefineInput(2, TList::Class());
  
  //Setup the parameters of the histograms
  _nBins_M0 = 500; _min_M0   = 0.;    _max_M0    = 5000.;  _width_M0 = (_max_M0-_min_M0)/_nBins_M0;
  _nBins_M1 = 500; _min_M1   = 0.;    _max_M1    = 5000.;  _width_M1 = (_max_M1-_min_M1)/_nBins_M1;
  _nBins_M2 = 500; _min_M2   = 0.;    _max_M2    = 5000.;  _width_M2 = (_max_M2-_min_M2)/_nBins_M2;
  _nBins_M3 = 500; _min_M3   = 0.;    _max_M3    = 5000.;  _width_M3 = (_max_M3-_min_M3)/_nBins_M3;
  _nBins_M4 = 100; _min_M4   = 0.;    _max_M4    = 100.;   _width_M4 = (_max_M4-_min_M4)/_nBins_M4;
  _nBins_M5 = 100; _min_M5   = 0.;    _max_M5    = 100.;   _width_M5 = (_max_M5-_min_M5)/_nBins_M5;
  _nBins_M6 = 100; _min_M6   = 0.;    _max_M6    = 100.;   _width_M6 = (_max_M6-_min_M6)/_nBins_M6;
  
  _min_vertexZ       = _vertexZMin;  
  _max_vertexZ       = _vertexZMax;  
  _width_vertexZ     = 0.5;
  _nBins_vertexZ     = int(0.5+ (_max_vertexZ - _min_vertexZ)/_width_vertexZ); 
  _nBins_pt_1        = int(0.5+ (_max_pt_1 -_min_pt_1 )/_width_pt_1); 
  _nBins_eta_1       = int(0.5+ (_max_eta_1-_min_eta_1)/_width_eta_1);  
  _width_phi_1       = (_max_phi_1  - _min_phi_1)  /_nBins_phi_1;
  _nBins_etaPhi_1    = _nBins_phi_1    * _nBins_eta_1;
  _nBins_etaPhiPt_1  = _nBins_etaPhi_1 * _nBins_pt_1;
  _nBins_zEtaPhiPt_1 = _nBins_vertexZ  * _nBins_etaPhiPt_1;
  
  _nBins_pt_2        =  int(0.5+ (_max_pt_2 -_min_pt_2 )/_width_pt_2);  
  _nBins_eta_2       = int(0.5+ (_max_eta_2-_min_eta_2)/_width_eta_2); 
  _width_phi_2       = (_max_phi_2  - _min_phi_2)  /_nBins_phi_2;
  _nBins_etaPhi_2    = _nBins_phi_2    * _nBins_eta_2;
  _nBins_etaPhiPt_2  = _nBins_etaPhi_2 * _nBins_pt_2;
  _nBins_zEtaPhiPt_2 = _nBins_vertexZ  * _nBins_etaPhiPt_2;
  _nBins_etaPhi_12   = _nBins_etaPhi_1 * _nBins_etaPhi_2;
  
  //setup the work arrays
  
  _id_1       = new int[arraySize];   
  _charge_1   = new int[arraySize]; 
  //_iPhi_1 = new int[arraySize]; 
  //_iEta_1 = new int[arraySize]; 
  _iEtaPhi_1  = new int[arraySize]; 
  _iPt_1      = new int[arraySize];  
  _pt_1       = new float[arraySize];    
  _px_1       = new float[arraySize];   
  _py_1       = new float[arraySize];   
  _pz_1       = new float[arraySize];   
  //_phi_1 = new float[arraySize];  
  //_eta_1 = new float[arraySize];  
  _correction_1 = new float[arraySize];
  //_dedx_1     = new float[arraySize];
  
  __n1_1_vsPt              = getDoubleArray(_nBins_pt_1,        0.);

  __n1_1_vsEtaPhi          = getDoubleArray(_nBins_etaPhi_1,    0.);
  __s1pt_1_vsEtaPhi        = getDoubleArray(_nBins_etaPhi_1,    0.);
  __n1_1_vsZEtaPhiPt       = getFloatArray(_nBins_zEtaPhiPt_1,  0.);
  
  //cout << "==========================================================================================" << endl;
  //cout << "=============== Booking for particle 1 done." << endl;
  //cout << "_requestedCharge_1: " << _requestedCharge_1 << endl;
  //cout << "_requestedCharge_2: " << _requestedCharge_2 << endl;
  
  if (_requestedCharge_2!=_requestedCharge_1)
    {
      //cout << " creating arrays for particle 2 with size: " << arraySize << endl;
    _sameFilter = 0;
    //particle 2
    _id_2       = new int[arraySize];   
    _charge_2   = new int[arraySize]; 
    //_iPhi_2   = new int[arraySize]; 
    //_iEta_2   = new int[arraySize]; 
    _iEtaPhi_2  = new int[arraySize]; 
    _iPt_2      = new int[arraySize];  
    _pt_2       = new float[arraySize];   
    _px_2       = new float[arraySize];   
    _py_2       = new float[arraySize];   
    _pz_2       = new float[arraySize];   
    //_phi_2    = new float[arraySize];  
    //_eta_2    = new float[arraySize];  
    _correction_2 = new float[arraySize];
    //_dedx_2       = new float[arraySize];
    
    __n1_2_vsPt              = getDoubleArray(_nBins_pt_2,        0.);
    __n1_2_vsEtaPhi          = getDoubleArray(_nBins_etaPhi_2,    0.);
    __s1pt_2_vsEtaPhi        = getDoubleArray(_nBins_etaPhi_2,    0.);
    __n1_2_vsZEtaPhiPt       = getFloatArray(_nBins_zEtaPhiPt_2, 0.);
    
    }
  
  __n2_12_vsPtPt           = getDoubleArray(_nBins_pt_1*_nBins_pt_2,0.);
  __n2_12_vsEtaPhi         = getFloatArray(_nBins_etaPhi_12,       0.);
  __s2ptpt_12_vsEtaPhi     = getFloatArray(_nBins_etaPhi_12,       0.);
  __s2PtN_12_vsEtaPhi      = getFloatArray(_nBins_etaPhi_12,       0.);
  __s2NPt_12_vsEtaPhi      = getFloatArray(_nBins_etaPhi_12,       0.);
  
  // Setup all the labels needed.
  
  part_1_Name   = "_1";
  part_2_Name   = "_2";
  pair_12_Name  = "_12";
  
  n1Name     = "n1";
  n2Name     = "n2";
  n1NwName   = "n1Nw";
  n2NwName   = "n2Nw";
  r1Name     = "r1";
  r2Name     = "r2";
  r3Name     = "r3";
  r2r1Name   = "r2r1";
  c2Name     = "c2";
  c3Name     = "c3";
  d3Name     = "d3";
  p3Name     = "p3";
  cName      = "sean";
  
  intR2Name       = "intR2";
  binCorrName     = "binCorr";
  intBinCorrName  = "intBinCorr";
  
  avgName      = "avg";
  sumName      = "sum";
  s1ptName     = "sumPt";
  s1ptNwName   = "sumPtNw";
  s1DptName    = "sumDpt";
  s2PtPtName   = "sumPtPt";
  s2PtPtNwName = "sumPtPtNw";
  s2DptDptName = "sumDptDpt";
  s2NPtName    = "sumNPt";
  s2NPtNwName  = "sumNPtNw";
  s2PtNName    = "sumPtN";
  s2NPtNwName  = "sumNPtNw";
  s2PtNNwName  = "sumPtNNw";
  ptName       = "avgPt";
  ptptName     = "avgPtPt";
  pt1pt1Name   = "avgPtavgPt";
  DptName      = "avgDpt";
  DptDptName   = "avgDptDpt";
  RDptDptName  = "relDptDpt"; // ratio of avgDptDpt by avgPt*avgPt
  nPtName      = "avgNpt";
  ptNName      = "avgPtN";
  seanName     = "seanC";
  
  _title_counts = "yield";
  
  _title_m0     = "M_{0}";
  _title_m1     = "M_{1}";
  _title_m2     = "M_{2}";
  _title_m3     = "M_{3}";
  _title_m4     = "V0Centrality";
  _title_m5     = "TrkCentrality";
  _title_m6     = "SpdCentrality";
  
  _title_eta_1       = "#eta_{1}";
  _title_phi_1       = "#varphi_{1} (radian)";
  _title_etaPhi_1    = "#eta_{1}#times#varphi_{1}";
  _title_pt_1        = "p_{t,1} (GeV/c)";
  _title_n_1         = "n_{1}";
  _title_SumPt_1     = "#Sigma p_{t,1} (GeV/c)";
  _title_AvgPt_1     = "#LT p_{t,1} #GT (GeV/c)";
  _title_AvgN_1      = "#LT n_{1} #GT";
  _title_AvgSumPt_1  = "#LT #Sigma p_{t,1} #GT (GeV/c)";
  
  _title_eta_2       = "#eta_{2}";
  _title_phi_2       = "#varphi_{2} (radian)";
  _title_etaPhi_2    = "#eta_{2}#times#varphi_{2}";
  _title_pt_2        = "p_{t,2} (GeV/c)";
  _title_n_2         = "n_{2}";
  _title_SumPt_2     = "#Sigma p_{t,1} (GeV/c)";
  _title_AvgPt_2     = "#LT p_{t,2} #GT (GeV/c)";
  _title_AvgN_2      = "#LT n_{2} #GT";
  _title_AvgSumPt_2  = "#LT #Sigma p_{t,2} #GT (GeV/c)";
  
  _title_etaPhi_12   = "#eta_{1}#times#varphi_{1}#times#eta_{2}#times#varphi_{2}";
  
  _title_AvgN2_12       = "#LT n_{2} #GT";;
  _title_AvgSumPtPt_12  = "#LT #Sigma p_{t,1}p_{t,2} #GT";;
  _title_AvgSumPtN_12   = "#LT #Sigma p_{t,1}N #GT";;
  _title_AvgNSumPt_12   = "#LT N#Sigma p_{t,2} #GT";;
  
  
  vsZ         = "_vsZ";
  vsM         = "_vsM";
  vsPt         = "_vsPt";
  vsPhi        = "_vsPhi"; 
  vsEta        = "_vsEta"; 
  vsEtaPhi     = "_vsEtaPhi"; 
  vsPtVsPt     = "_vsPtVsPt";
  
  
  if (_useWeights)
    {
    int iZ, iEtaPhi, iPt;
    int iZ1,iEtaPhi1,iPt1;
    int a, b;
    if (_weight_1)
      {
      _correctionWeight_1 = new float[_nBins_vertexZ*_nBins_etaPhi_1*_nBins_pt_1];
      a = _nBins_pt_1;
      b = _nBins_etaPhi_1*_nBins_pt_1;
      for (iZ=0,iZ1=1; iZ<_nBins_vertexZ; iZ++, iZ1++)
        {
        for (iEtaPhi=0,iEtaPhi1=1; iEtaPhi<_nBins_etaPhi_1; iEtaPhi++, iEtaPhi1++)
          {
          for (iPt=0,iPt1=1; iPt<_nBins_pt_1; iPt++, iPt1++)
            {
            _correctionWeight_1[iZ*b+iEtaPhi*a+iPt] = _weight_1->GetBinContent(iZ1,iEtaPhi1,iPt1);
            }      
          }
        }
      } // _weight_1
    else
      {
      AliError("AliDptDptInMC:: _weight_1 is a null pointer.");
      return;
      }
    if (!_sameFilter) 
      {
      if (_weight_2)
        {
        _correctionWeight_2 = new float[_nBins_vertexZ*_nBins_etaPhi_2*_nBins_pt_2];
        a = _nBins_pt_2;
        b = _nBins_etaPhi_2*_nBins_pt_2;
        for (iZ=0,iZ1=1; iZ<_nBins_vertexZ; iZ++, iZ1++)
          {
          for (iEtaPhi=0,iEtaPhi1=1; iEtaPhi<_nBins_etaPhi_2; iEtaPhi++, iEtaPhi1++)
            {
            for (iPt=0,iPt1=1; iPt<_nBins_pt_2; iPt++, iPt1++)
              {
              _correctionWeight_2[iZ*b+iEtaPhi*a+iPt] = _weight_2->GetBinContent(iZ1,iEtaPhi1,iPt1);
              }      
            }
          }
        } // _weight_2
      else
        {
        AliError("AliDptDptInMC:: _weight_1 is a null pointer.");
        return;
        }
      }
    }
  
  createHistograms();
  PostData(0,_outputHistoList);
  
  //cout<< "AliDptDptInMC::CreateOutputObjects() DONE " << endl;
  
}

void  AliDptDptInMC::createHistograms()
{
  AliInfo(" AliDptDptInMC::createHistoHistograms() Creating Event Histos");
  TString name;
  
  name = "eventAccounting";
  
  // bin index : what it is...
  //         0 :  number of event submitted
  //         1 :  number accepted by centrality cut
  //         2 :  number accepted by centrality cut and z cut
  //         3 :  total number of particles that satisfied filter 1
  //         4 :  total number of particles that satisfied filter 2
  _eventAccounting      = createHisto1D(name,name,10, -0.5, 9.5, "event Code", _title_counts);
  
  name = "m0"; _m0      = createHisto1D(name,name,_nBins_M1, _min_M1, _max_M1, _title_m0, _title_counts);
  name = "m1"; _m1      = createHisto1D(name,name,_nBins_M1, _min_M1, _max_M1, _title_m1, _title_counts);
  name = "m2"; _m2      = createHisto1D(name,name,_nBins_M2, _min_M2, _max_M2, _title_m2, _title_counts);
  name = "m3"; _m3      = createHisto1D(name,name,_nBins_M3, _min_M3, _max_M3, _title_m3, _title_counts);
  name = "m4"; _m4      = createHisto1D(name,name,_nBins_M4, _min_M4, _max_M4, _title_m4, _title_counts);
  name = "m5"; _m5      = createHisto1D(name,name,_nBins_M5, _min_M5, _max_M5, _title_m5, _title_counts);
  name = "m6"; _m6      = createHisto1D(name,name,_nBins_M6, _min_M6, _max_M6, _title_m6, _title_counts);
  name = "zV"; _vertexZ = createHisto1D(name,name,_nBins_vertexZ, _min_vertexZ, _max_vertexZ, "z-Vertex (cm)", _title_counts);
  

  name = "Eta";     _etadis   = createHisto1F(name,name, 200, -1.0, 1.0, "#eta","counts");
  name = "Phi";     _phidis   = createHisto1F(name,name, 360, 0.0, 6.4, "#phi","counts");
  name = "DCAz";    _dcaz     = createHisto1F(name,name, 340, -3.3, 3.3, "dcaZ","counts");
  name = "DCAxy";   _dcaxy    = createHisto1F(name,name, 100, -0.1, 2.5, "dcaXY","counts");

  name = "Eta2";     _etadis2   = createHisto1F(name,name, 200, -1.0, 1.0, "#eta","counts");
  name = "Phi2";     _phidis2   = createHisto1F(name,name, 360, 0.0, 6.4, "#phi","counts");
  name = "DCAz2";    _dcaz2     = createHisto1F(name,name, 340, -3.3, 3.3, "dcaZ","counts");
  name = "DCAxy2";   _dcaxy2    = createHisto1F(name,name, 100, -0.1, 2.5, "dcaXY","counts");

  name = "Eta3";     _etadis3   = createHisto1F(name,name, 200, -1.0, 1.0, "#eta","counts");
  name = "Phi3";     _phidis3   = createHisto1F(name,name, 360, 0.0, 6.4, "#phi","counts");
  name = "DCAz3";    _dcaz3     = createHisto1F(name,name, 340, -3.3, 3.3, "dcaZ","counts");
  name = "DCAxy3";   _dcaxy3    = createHisto1F(name,name, 100, -0.1, 2.5, "dcaXY","counts");

  //name = "Eta";     _etadis   = createHisto1F(name,name, 250, 0.0, 2.5, "#eta","counts"); //temporaryly
  //name = "Phi";     _phidis   = createHisto1F(name,name, 250, 0.0, 2.5, "#phi","counts");
  //name = "DCAz";    _dcaz     = createHisto1F(name,name, 250, 0.0, 2.5, "dcaZ","counts");

  //name = "Nclus1";   _Ncluster1    = createHisto1F(name,name, 200, 0, 200, "Ncluster1","counts");
  //name = "Nclus2";   _Ncluster2    = createHisto1F(name,name, 200, 0, 200, "Ncluster2","counts");
  
  if (_singlesOnly)
    {
    name = n1Name+part_1_Name+vsPt;              _n1_1_vsPt              = createHisto1F(name,name, _nBins_pt_1,  _min_pt_1,  _max_pt_1,   _title_pt_1,  _title_AvgN_1);

    name = n1Name+part_1_Name+vsZ+vsEtaPhi+vsPt; _n1_1_vsZVsEtaVsPhiVsPt = createHisto3F(name,name, _nBins_vertexZ,_min_vertexZ,_max_vertexZ, _nBins_etaPhi_1, 0., double(_nBins_etaPhi_1), _nBins_pt_1, _min_pt_1, _max_pt_1, "zVertex", _title_etaPhi_1,  _title_pt_1);
    //name = "dedxVsP_1";                          _dedxVsP_1              = createHisto2F(name,name,400,-2.,2.,120,0.,120.,"p (GeV/c)", "dedx", "counts");
    //name = "corrDedxVsP_1";                      _corrDedxVsP_1          = createHisto2F(name,name,400,-2.,2.,120,0.,120.,"p (GeV/c)", "dedx", "counts");
    //name = "betaVsP_1";                          _betaVsP_1              = createHisto2F(name,name,400,-2.,2.,120,0.5,1.1,"p (GeV/c)", "beta", "counts");

    name = n1Name+part_2_Name+vsPt;              _n1_2_vsPt              = createHisto1F(name,name, _nBins_pt_2,  _min_pt_2,  _max_pt_2,   _title_pt_2,  _title_AvgN_2);
    name = n1Name+part_2_Name+vsZ+vsEtaPhi+vsPt; _n1_2_vsZVsEtaVsPhiVsPt = createHisto3F(name,name, _nBins_vertexZ,_min_vertexZ,_max_vertexZ, _nBins_etaPhi_2, 0., double(_nBins_etaPhi_2), _nBins_pt_2, _min_pt_2, _max_pt_2, "zVertex", _title_etaPhi_2,  _title_pt_2);
    //name = "dedxVsP_2";                          _dedxVsP_2              = createHisto2F(name,name,400,-2.,2.,120,0.,120.,"p (GeV/c)", "dedx", "counts");
    //name = "corrDedxVsP_2";                      _corrDedxVsP_2          = createHisto2F(name,name,400,-2.,2.,120,0.,120.,"p (GeV/c)", "dedx", "counts");
    //name = "betaVsP_2";                          _betaVsP_2              = createHisto2F(name,name,400,-2.,2.,120,0.5,1.1,"p (GeV/c)", "beta", "counts");

    }
  else
    {
    name = n1Name+part_1_Name+vsEtaPhi;       _n1_1_vsEtaVsPhi      = createHisto2F(name,name, _nBins_eta_1, _min_eta_1, _max_eta_1,  _nBins_phi_1, _min_phi_1, _max_phi_1,  _title_eta_1,  _title_phi_1,  _title_AvgN_1);
    name = s1ptName+part_1_Name+vsEtaPhi;     _s1pt_1_vsEtaVsPhi    = createHisto2F(name,name, _nBins_eta_1, _min_eta_1, _max_eta_1,  _nBins_phi_1, _min_phi_1, _max_phi_1,  _title_eta_1,  _title_phi_1,  _title_AvgSumPt_1);
    name = n1Name+part_1_Name+vsM;            _n1_1_vsM             = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgN_1);
    name = s1ptName+part_1_Name+vsM;          _s1pt_1_vsM           = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPt_1);
    name = n1NwName+part_1_Name+vsM;          _n1Nw_1_vsM           = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgN_1);
    name = s1ptNwName+part_1_Name+vsM;        _s1ptNw_1_vsM         = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPt_1);

    name = n1Name+part_2_Name+vsEtaPhi;       _n1_2_vsEtaVsPhi      = createHisto2F(name,name, _nBins_eta_2, _min_eta_2, _max_eta_2,  _nBins_phi_2, _min_phi_2, _max_phi_2,  _title_eta_2,  _title_phi_2,  _title_AvgN_2);
    name = s1ptName+part_2_Name+vsEtaPhi;     _s1pt_2_vsEtaVsPhi    = createHisto2F(name,name, _nBins_eta_2, _min_eta_2, _max_eta_2,  _nBins_phi_2, _min_phi_2, _max_phi_2,  _title_eta_2,  _title_phi_2,  _title_AvgSumPt_2);
    name = n1Name+part_2_Name + vsM;          _n1_2_vsM             = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgN_2);
    name = s1ptName+part_2_Name + vsM;        _s1pt_2_vsM           = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPt_2);
    name = n1NwName+part_2_Name+vsM;          _n1Nw_2_vsM           = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgN_1);
    name = s1ptNwName+part_2_Name+vsM;        _s1ptNw_2_vsM         = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPt_1);

    name = n2Name+pair_12_Name+vsEtaPhi;      _n2_12_vsEtaPhi       = createHisto1F(name,name, _nBins_etaPhi_12, 0.,        double(_nBins_etaPhi_12), _title_etaPhi_12, _title_AvgN2_12);        
    name = s2PtPtName+pair_12_Name + vsEtaPhi;_s2PtPt_12_vsEtaPhi   = createHisto1F(name,name, _nBins_etaPhi_12, 0.,        double(_nBins_etaPhi_12), _title_etaPhi_12,  _title_AvgSumPtPt_12);    
    name = s2PtNName+pair_12_Name + vsEtaPhi; _s2PtN_12_vsEtaPhi    = createHisto1F(name,name, _nBins_etaPhi_12, 0.,        double(_nBins_etaPhi_12), _title_etaPhi_12,  _title_AvgSumPtN_12);    
    name = s2NPtName+pair_12_Name + vsEtaPhi; _s2NPt_12_vsEtaPhi    = createHisto1F(name,name, _nBins_etaPhi_12, 0.,        double(_nBins_etaPhi_12), _title_etaPhi_12,  _title_AvgNSumPt_12);    
    name = n2Name+pair_12_Name+vsPtVsPt;      _n2_12_vsPtVsPt       = createHisto2F(name,name, _nBins_pt_1, _min_pt_1, _max_pt_1, _nBins_pt_2, _min_pt_2, _max_pt_2, _title_pt_1, _title_pt_2, _title_AvgN2_12);        
    
    name = n2Name+pair_12_Name + vsM;         _n2_12_vsM            = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgN2_12);        
    name = s2PtPtName+pair_12_Name + vsM;     _s2PtPt_12_vsM        = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPtPt_12);    
    name = s2PtNName+pair_12_Name + vsM;      _s2PtN_12_vsM         = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPtN_12);       
    name = s2NPtName+pair_12_Name + vsM;      _s2NPt_12_vsM         = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgNSumPt_12);     
    
    name = n2NwName+pair_12_Name + vsM;       _n2Nw_12_vsM          = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgN2_12);        
    name = s2PtPtNwName+pair_12_Name + vsM;   _s2PtPtNw_12_vsM      = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPtPt_12);    
    name = s2PtNNwName+pair_12_Name + vsM;    _s2PtNNw_12_vsM       = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgSumPtN_12);       
    name = s2NPtNwName+pair_12_Name + vsM;    _s2NPtNw_12_vsM       = createProfile(name,name, _nBins_M4, _min_M4, _max_M4, _title_m4, _title_AvgNSumPt_12);     
    
    name = "mInv";     _invMass     = createHisto1F(name,name, 500, 0., 1.000, "M_{inv}","counts"); 
    name = "mInvElec"; _invMassElec = createHisto1F(name,name, 500, 0., 1.000, "M_{inv}","counts"); 
    }
  
  AliInfo(" AliDptDptInMC::createHistoHistograms() All Done"); 
}
//-----------------------//

void  AliDptDptInMC::finalizeHistograms()
{
  
  AliInfo("AliDptDptInMC::finalizeHistograms() starting");
  AliInfo(Form("CorrelationAnalyzers::finalizeHistograms()   _eventCount : %d",int(_eventCount)));
  if (_singlesOnly)
    {
    if (_sameFilter)
      {
      fillHistoWithArray(_n1_1_vsPt,              __n1_1_vsPt,        _nBins_pt_1);
      fillHistoWithArray(_n1_1_vsZVsEtaVsPhiVsPt, __n1_1_vsZEtaPhiPt, _nBins_vertexZ, _nBins_etaPhi_1, _nBins_pt_1);
      fillHistoWithArray(_n1_2_vsPt,              __n1_1_vsPt,        _nBins_pt_1);
      fillHistoWithArray(_n1_2_vsZVsEtaVsPhiVsPt, __n1_1_vsZEtaPhiPt, _nBins_vertexZ, _nBins_etaPhi_1, _nBins_pt_1);
      }
    else
      {
      fillHistoWithArray(_n1_1_vsPt,              __n1_1_vsPt,        _nBins_pt_1);
      fillHistoWithArray(_n1_1_vsZVsEtaVsPhiVsPt, __n1_1_vsZEtaPhiPt, _nBins_vertexZ, _nBins_etaPhi_1, _nBins_pt_1);
      fillHistoWithArray(_n1_2_vsPt,              __n1_2_vsPt,        _nBins_pt_2);
      fillHistoWithArray(_n1_2_vsZVsEtaVsPhiVsPt, __n1_2_vsZEtaPhiPt, _nBins_vertexZ, _nBins_etaPhi_2, _nBins_pt_2);
      }
    }
  else
    {
    if (_sameFilter)
      {
      fillHistoWithArray(_n1_1_vsEtaVsPhi,        __n1_1_vsEtaPhi,    _nBins_eta_1,   _nBins_phi_1);
      fillHistoWithArray(_s1pt_1_vsEtaVsPhi,      __s1pt_1_vsEtaPhi,  _nBins_eta_1,   _nBins_phi_1);
      fillHistoWithArray(_n1_2_vsEtaVsPhi,        __n1_1_vsEtaPhi,    _nBins_eta_1,   _nBins_phi_1);
      fillHistoWithArray(_s1pt_2_vsEtaVsPhi,      __s1pt_1_vsEtaPhi,  _nBins_eta_1,   _nBins_phi_1);
      }
    else
      {
      fillHistoWithArray(_n1_1_vsEtaVsPhi,        __n1_1_vsEtaPhi,    _nBins_eta_1,   _nBins_phi_1);
      fillHistoWithArray(_s1pt_1_vsEtaVsPhi,      __s1pt_1_vsEtaPhi,  _nBins_eta_1,   _nBins_phi_1);
      fillHistoWithArray(_n1_2_vsEtaVsPhi,        __n1_2_vsEtaPhi,    _nBins_eta_2,   _nBins_phi_2);
      fillHistoWithArray(_s1pt_2_vsEtaVsPhi,      __s1pt_2_vsEtaPhi,  _nBins_eta_2,   _nBins_phi_2);
      }
    fillHistoWithArray(_n2_12_vsEtaPhi,     __n2_12_vsEtaPhi,     _nBins_etaPhi_12);
    fillHistoWithArray(_s2PtPt_12_vsEtaPhi, __s2ptpt_12_vsEtaPhi, _nBins_etaPhi_12);
    fillHistoWithArray(_s2PtN_12_vsEtaPhi,  __s2PtN_12_vsEtaPhi,  _nBins_etaPhi_12);
    fillHistoWithArray(_s2NPt_12_vsEtaPhi,  __s2NPt_12_vsEtaPhi,  _nBins_etaPhi_12);
    fillHistoWithArray(_n2_12_vsPtVsPt,     __n2_12_vsPtPt,       _nBins_pt_1,    _nBins_pt_2);

    }  
  AliInfo("AliDptDptInMC::finalizeHistograms()  Done ");
}
//--------------//


void  AliDptDptInMC::UserExec(Option_t */*option*/)
{
  

  int    k1,k2;
  int    iPhi, iEta, iEtaPhi, iPt, charge;
  float  q, phi, pt, eta, corr, corrPt, px, py, pz;
  int    ij;
  int    id_1, q_1, iEtaPhi_1, iPt_1;
  float  pt_1, px_1, py_1, pz_1, corr_1;
  int    id_2, q_2, iEtaPhi_2, iPt_2;
  float  pt_2, px_2, py_2, pz_2, corr_2;
  float  ptpt;
  int    iVertex, iVertexP1, iVertexP2;
  int    iZEtaPhiPt;
  float  massElecSq = 2.5e-7;
  const  AliAODVertex*  vertex;
  //int    nClus;                                                                                                                           
  bool   bitOK;

  AliAnalysisManager* manager = AliAnalysisManager::GetAnalysisManager();
  if (!manager) {
    //cout<<"ERROR: Analysis manager not found."<<endl;                                                                                     
    return;
  }
  //coneect to the inputHandler------------                                                                                                 
  AliAODInputHandler* inputHandler = dynamic_cast<AliAODInputHandler*> (manager->GetInputEventHandler());
  if (!inputHandler) {
    return;
  }

  fAODEvent = dynamic_cast<AliAODEvent*>(InputEvent());
  if (!fAODEvent)
    {
      return;
    }
  fPIDResponse =inputHandler->GetPIDResponse();
  if (!fPIDResponse){
    AliFatal("This Task needs the PID response attached to the inputHandler");
    return;
  }
  // count all events looked at here                                                                                                        
  _eventCount++;

  if (_eventAccounting)
    {
      _eventAccounting->Fill(0);// count all calls to this function                                                                         
    }
  else
    {
      return;
    }

  _eventAccounting->Fill(1);// count all calls to this function with a valid pointer                                                        
  //reset single particle counters                                                                                                          
  k1 = k2 = 0;
  __n1_1 = __n1_2 = __s1pt_1 = __s1pt_2 = __n1Nw_1 = __n1Nw_2 = __s1ptNw_1 = __s1ptNw_2 = 0;

  float v0Centr  = -999.;
  float v0ACentr  = -999.;
  float trkCentr = -999.;
  float spdCentr = -999.;
  //float vertexX  = -999;
  //float vertexY  = -999;
  float vertexZ  = -999;
  //float vertexXY = -999;

  float centrality = -999;
  //Double_t nSigma =-999;  
  if(fAODEvent)
    {

      //Centrality                                                                                                                          
      //AliAODHeader* centralityObject = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP();                                                          
      AliCentrality* centralityObject =  ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP();
      if (centralityObject)
        {
          v0Centr  = centralityObject->GetCentralityPercentile("V0M");
          v0ACentr  = centralityObject->GetCentralityPercentile("V0A");
          trkCentr = centralityObject->GetCentralityPercentile("TRK");
          spdCentr = centralityObject->GetCentralityPercentile("CL1");

        }

      _mult4    = v0Centr;
      _mult4a   = v0ACentr;
      _mult5    = trkCentr;
      _mult6    = spdCentr;
      _field    = fAODEvent->GetMagneticField();
      //_centralityMethod                                                                                                                   
      switch (_centralityMethod)
        {
        case 0: centrality = _mult0; break;
        case 1: centrality = _mult1; break;
        case 2: centrality = _mult2; break;
        case 3: centrality = _mult3; break;
	case 4: centrality = _mult4; break;
        case 5: centrality = _mult5; break;
        case 6: centrality = _mult6; break;
        case 7: centrality = _mult4a; break;
        }

      if ( centrality < _centralityMin ||  centrality > _centralityMax ) return;

      _eventAccounting->Fill(2);// count all events with right centrality                                                                   

      vertex = (AliAODVertex*) fAODEvent->GetPrimaryVertex();

      if(vertex->GetNContributors() > 0)
        {
          //vertexX = vertex->GetX();
          //vertexY = vertex->GetY();
          vertexZ = vertex->GetZ();
          //vertexXY = sqrt(vertexX*vertexX+vertexY*vertexY);
        }
      if (vertexZ  < _vertexZMin  || vertexZ  > _vertexZMax ) return;

      iVertex = int((vertexZ-_min_vertexZ)/_width_vertexZ);
      iVertexP1 = iVertex*_nBins_etaPhiPt_1;
      iVertexP2 = iVertex*_nBins_etaPhiPt_2;
      if (iVertex<0 || iVertex>=_nBins_vertexZ)
        {
          AliError("AliTaskDptCorrMC::Exec(Option_t *option) iVertex<0 || iVertex>=_nBins_vertexZ ");
          return;
        }
      _eventAccounting->Fill(3);// count all calls to this function with a valid pointer                                                    
      //======================                 
      //***********************************            
      //MC AOD Truth 
      if(fAnalysisType == "MCAOD")
        { 
	  AliMCEvent* mcEvent = MCEvent();
	  _nTracks = mcEvent->GetNumberOfTracks();
          _mult3    = _nTracks;
          //loop over tracks starts here                                                                                                       
          for (int iTrack=0; iTrack< _nTracks; iTrack++)
            {
              AliAODMCParticle *aodTrack = (AliAODMCParticle*) mcEvent->GetTrack(iTrack);

              if (!aodTrack)
                {
                  AliError(Form("AliTaskDptCorrMC::Exec(Option_t *option) No track ofr iTrack=%d", iTrack));
                  continue;
                }
	      
	      if(!aodTrack->IsPhysicalPrimary()) continue;

	      q      = aodTrack->Charge();
              charge = int(q);
              phi    = aodTrack->Phi();
              pt     = aodTrack->Pt();
              px     = aodTrack->Px();
              py     = aodTrack->Py();
              pz     = aodTrack->Pz();
              eta    = aodTrack->Eta();
	      // Kinematics cuts from ESD track cuts                                                           
	      if( pt < 0.2 || pt > 2.0)      continue;
	      if( eta < -0.8 || eta > 0.8)  continue;

	      if(q == 0) continue;  
	      _dcaz->Fill(pt); //AllCh 

	      if(TMath::Abs(aodTrack->GetPdgCode()) == 211)
		{
		  _etadis->Fill(pt); //pion
		} 
	      if(TMath::Abs(aodTrack->GetPdgCode()) == 321)
		{
		  _phidis->Fill(pt); //kaon
		} 

	      //_etadis->Fill(eta);
              //_phidis->Fill(phi); 

	      if(fExcludeResonancesInMC)
                {
                  //cout<<"***************Prabhat on Weak Decay Particles ************"<<endl;                               
		  Int_t gMotherIndex = aodTrack->GetMother();
                  if(gMotherIndex != -1) {
                    AliAODMCParticle* motherTrack = dynamic_cast<AliAODMCParticle *>(mcEvent->GetTrack(gMotherIndex));
                    if(motherTrack) {
                      Int_t pdgCodeOfMother = motherTrack->GetPdgCode();
		      
                      if(pdgCodeOfMother == 311  ||
                         pdgCodeOfMother == -311 ||
                         pdgCodeOfMother == 310  ||
                         pdgCodeOfMother == 3122 ||
                         pdgCodeOfMother == -3122 ||
			 pdgCodeOfMother == 111 ||
			 pdgCodeOfMother == 22 ) continue;
                    }
                  }
                }
	      //Exclude electrons with PDG                                                                                   
                                                                                                                               
              if(fExcludeElectronsInMC) {
                if(TMath::Abs(aodTrack->GetPdgCode()) == 11) continue;
              }

	      //Particle loop 1st particle
              if (q > 0)
                {
		  
                  iPhi   = int( phi/_width_phi_1);
		  
                  if (iPhi<0 || iPhi>=_nBins_phi_1 )
                    {
                      AliWarning("AliTaskDptCorrMC::analyze() iPhi<0 || iPhi>=_nBins_phi_1");
                      return;
                    }

                  iEta    = int((eta-_min_eta_1)/_width_eta_1);
                  if (iEta<0 || iEta>=_nBins_eta_1)
                    {
                      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iEta: %d", iEta));
                      continue;
                    }

		  iPt     = int((pt -_min_pt_1 )/_width_pt_1 );
                  if (iPt<0  || iPt >=_nBins_pt_1)
                    {
                      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iPt: %d",iPt));
                      continue;
                    }
                  iEtaPhi = iEta*_nBins_phi_1+iPhi;
                  iZEtaPhiPt = iVertexP1 + iEtaPhi*_nBins_pt_1 + iPt;

                  if (_correctionWeight_1)
                    corr = _correctionWeight_1[iZEtaPhiPt];
                  else
                    corr = 1;
                  if (iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_1)
                    {
                      AliWarning("AliTaskDptCorrMC::analyze(AliceEvent * event) iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_1");
                      continue;
                    }

		  if (_singlesOnly)
                    {

                      __n1_1_vsPt[iPt]               += corr;          //cout << "step 15" << endl;                                            
                      
                      __n1_1_vsZEtaPhiPt[iZEtaPhiPt] += corr;       //cout << "step 12" << endl;                                               

                    }

                  else
                    {
                      corrPt                      = corr*pt;
                      _id_1[k1]                   = iTrack;
                      _charge_1[k1]               = charge;
                      _iEtaPhi_1[k1]              = iEtaPhi;
                      _iPt_1[k1]                  = iPt;
                      _pt_1[k1]                   = pt;
                      _px_1[k1]                   = px;
                      _py_1[k1]                   = py;
                      _pz_1[k1]                   = pz;
                      _correction_1[k1]           = corr;
                      __n1_1                      += corr;
                      __n1_1_vsEtaPhi[iEtaPhi]    += corr;
                      __s1pt_1                    += corrPt;
                      __s1pt_1_vsEtaPhi[iEtaPhi]  += corrPt;
                      __n1Nw_1                    += 1;
                      __s1ptNw_1                  += pt;
                      ++k1;
                      if (k1>=arraySize)
                        {
                          AliError(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) k1 >=arraySize; arraySize: %d",arraySize));
                          return;
                        }
                    }
                }

	      if (q < 0)
                {
                  iPhi   = int( phi/_width_phi_2);

                  if (iPhi<0 || iPhi>=_nBins_phi_2 )
                    {
                      AliWarning("AliTaskDptCorrMC::analyze() iPhi<0 || iPhi>=_nBins_phi_1");
                      return;
                    }

                  iEta    = int((eta-_min_eta_2)/_width_eta_2);
                  if (iEta<0 || iEta>=_nBins_eta_2)
                    {
                      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iEta: %d", iEta));
                      continue;
                    }

		  iPt     = int((pt -_min_pt_2 )/_width_pt_2 );
                  if (iPt<0  || iPt >=_nBins_pt_2)
                    {
                      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iPt: %d",iPt));
                      continue;
                    }

                  iEtaPhi = iEta*_nBins_phi_2+iPhi;
                  iZEtaPhiPt = iVertexP2 + iEtaPhi*_nBins_pt_2 + iPt;
                  if (iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_2)
                    {
                      AliWarning("AliTaskDptCorrMC::analyze(AliceEvent * event) iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_2");
                      continue;
                    }

		  if (_correctionWeight_2)
                    corr = _correctionWeight_2[iZEtaPhiPt];
                  else
                    corr = 1;
                  //dpt = pt - (charge>0) ? _avgPt_vsEtaPhi_2p[iEtaPhi] : _avgPt_vsEtaPhi_2m[iEtaPhi];                                         

                  if (_singlesOnly)
                    {
                      __n1_2_vsPt[iPt]               += corr;          //cout << "step 15" << endl;                                            
                      __n1_2_vsZEtaPhiPt[iZEtaPhiPt] += corr;       //cout << "step 12" << endl;                                               
                    }
		  else
                    {
                      corrPt                      = corr*pt;
                      _id_2[k2]                   = iTrack;         //cout << "step 1" << endl;                                                
                      _charge_2[k2]               = charge;         //cout << "step 2" << endl;                                                
                      _iEtaPhi_2[k2]              = iEtaPhi;        //cout << "step 3" << endl;                                                
                      _iPt_2[k2]                  = iPt;            //cout << "step 4" << endl;                                                
                      _pt_2[k2]                   = pt;             //cout << "step 5" << endl;                                                
                      _px_2[k2]                   = px;             //cout << "step 6" << endl;                                                
                      _py_2[k2]                   = py;             //cout << "step 7" << endl;                                                
                      _pz_2[k2]                   = pz;             //cout << "step 8" << endl;                                                
                      _correction_2[k2]           = corr;           //cout << "step 9" << endl;                                                
                      __n1_2                      += corr;          //cout << "step 10" << endl;                                               
                      __s1pt_2                    += corrPt;        //cout << "step 13" << endl;                                               
                      __n1Nw_2                    += 1;
                      __n1_2_vsEtaPhi[iEtaPhi]    += corr;          //cout << "step 11" << endl;                                               
                      __s1pt_2_vsEtaPhi[iEtaPhi]  += corrPt;        //cout << "step 14" << endl;                                               
                      __s1ptNw_2                  += pt;
                      ++k2;
                      if (k2>=arraySize)
                        {
                          AliWarning(Form("-W- k2 >=arraySize; arraySize: %d",arraySize));
                          return;
                        }
                    }

                  //cout << "done with track" << endl;                                                                                         
                } //iTrack                                                                                                                     
            } //data aod loop                                                                                                                  
	  } //MC AOD loop ends                                                                                                                   

      //***********************************************
      //MC AOD Reconstructed tracks
      if(fAnalysisType == "MCAODreco")
	{
	  
	  fArrayMC = dynamic_cast<TClonesArray*>(fAODEvent->FindListObject(AliAODMCParticle::StdBranchName()));
	  if (!fArrayMC) {
	    AliError("No array of MC particles found !!!");
	  }
	  
	  AliMCEvent* mcEvent = MCEvent();
	  if (!mcEvent) {
	    AliError("ERROR: Could not retrieve MC event");
	  }

	  TExMap *trackMap = new TExMap();//Mapping matrix---- 
	  //1st loop track for Global tracks                                                                                                    
	  for(Int_t i = 0; i < fAODEvent->GetNumberOfTracks(); i++)
	    {
	      AliAODTrack* aodTrack = dynamic_cast<AliAODTrack *>(fAODEvent->GetTrack(i));
	      if(!aodTrack) {
		AliError(Form("ERROR: Could not retrieve AODtrack %d",i));
		continue;
	      }
	      Int_t gID = aodTrack->GetID();
	      if (aodTrack->TestFilterBit(1)) trackMap->Add(gID, i);//Global tracks                                                             
	    }
	  
	  AliAODTrack* newAodTrack;
	  
	  for (int iTrack=0; iTrack< fAODEvent->GetNumberOfTracks(); iTrack++)
	    {
	      
	      AliAODTrack *t = dynamic_cast<AliAODTrack *>(fAODEvent->GetTrack(iTrack));
	      
	      if(!t) {
		AliError(Form("ERROR: Could not retrieve AODtrack %d",iTrack));
		continue;
	      }
	      
	      bitOK  = t->TestFilterBit(_trackFilterBit);
	      if (!bitOK) continue;
	      Int_t gID = t->GetID();
	      newAodTrack = gID >= 0 ?t : dynamic_cast<AliAODTrack*>(fAODEvent->GetTrack(trackMap->GetValue(-1-gID)));
              if(!newAodTrack) {
                AliFatal("Not a standard AOD?");
              }
	      
	      q      = t->Charge();
	      charge = int(q);
	      phi    = t->Phi();
	      pt     = t->Pt();
	      px     = t->Px();
	      py     = t->Py();
	      pz     = t->Pz();
	      eta    = t->Eta();
	      //Float_t dcaXY = t->DCA();  
	      //Float_t dcaZ  = t->ZAtDCA(); 
	      
	      // get the electron nsigma                                                                                                
	      Double_t nSigma = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
	      Double_t nSigmaPions   = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));
              Double_t nSigmaKaons   = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kKaon));
              Double_t nSigmaProtons = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kProton));
	      
	      //Make the decision based on the n-sigma of electrons exclusively 
	      if(nSigma < fNSigmaCut
                 && nSigmaPions   > fNSigmaCut
                 && nSigmaKaons   > fNSigmaCut
                 && nSigmaProtons > fNSigmaCut ) continue;
	      

	      if(q == 0) continue;	      
	      // Kinematics cuts                                                                                 
	      if( pt < 0.2 || pt > 2.0)      continue;
	      if( eta < _min_eta_1 || eta > _max_eta_1)  continue;
	      
	      
	      /*
		_dcaz->Fill(pt);
		if (nSigmaPions < fNSigmaCut)
		{
		_etadis->Fill(pt);
		}
		if (nSigmaKaons < fNSigmaCut)
		{
		_phidis->Fill(pt);
		} 
	      */ //particle ratio calculation
	      
	      Double_t pos[3];
	      newAodTrack->GetXYZ(pos);
	      
	      //Double_t DCAX = pos[0] - vertexX;
	      //Double_t DCAY = pos[1] - vertexY;
	      Double_t DCAZ = pos[2] - vertexZ;
	      
	      //Double_t DCAXY = TMath::Sqrt((DCAX*DCAX) + (DCAY*DCAY));
	      
	      /*if (DCAZ     <  _dcaZMin ||
		DCAZ     >  _dcaZMax ||
		DCAXY    >  _dcaXYMax ) continue;
	      */
	      
	      Int_t label = TMath::Abs(t->GetLabel());                                                                                              
	      AliAODMCParticle *AODmcTrack = (AliAODMCParticle*) fArrayMC->At(label); 
	      //------------Test of tracks QAs------------	      
	      if (AODmcTrack->IsPhysicalPrimary())
		{
		  _dcaz->Fill(DCAZ);                                                           
		  //_dcaxy->Fill(DCAXY);                                                         
		  //_etadis->Fill(eta);                                                          
		  _phidis->Fill(phi); 
		}
	      if (AODmcTrack->IsSecondaryFromWeakDecay())
		{
		  _dcaz2->Fill(DCAZ);                                                           
		  //_dcaxy2->Fill(DCAXY);                                                         
		  //_etadis2->Fill(eta);                                                          
		  _phidis2->Fill(phi); 
		}
	      
	      if (AODmcTrack->IsSecondaryFromMaterial())
		{
		  _dcaz3->Fill(DCAZ);                                                           
		  //_dcaxy3->Fill(DCAXY);                                                         
		  //_etadis3->Fill(eta);                                                          
		  _phidis3->Fill(phi); 
		}
	      //---------------------------------
	      
	      /* //W/Wo Secondaries
	      //if (!AODmcTrack->IsPhysicalPrimary()) continue;
	      //cout<<"***************Prabhat on Weak Decay Particles ************"<<endl;
	      if(fExcludeResonancesInMC)
	      {
	      Int_t label = TMath::Abs(t->GetLabel());
	      AliAODMCParticle *AODmcTrack = (AliAODMCParticle*) fArrayMC->At(label);
	      
	      Int_t gMotherIndex = AODmcTrack->GetMother();
	      if(gMotherIndex != -1) {
	      AliAODMCParticle* motherTrack = dynamic_cast<AliAODMCParticle *>(mcEvent->GetTrack(gMotherIndex));
	      if(motherTrack) {
	      Int_t pdgCodeOfMother = motherTrack->GetPdgCode();
	      
	      if(pdgCodeOfMother == 311  ||
	      pdgCodeOfMother == -311 ||
	      pdgCodeOfMother == 310  ||
	      pdgCodeOfMother == 3122 ||
	      pdgCodeOfMother == -3122 ||
	      pdgCodeOfMother == 111 ||
	      pdgCodeOfMother == 22 ) continue;
	      }
	      }
	      } */
	      
	      
	      //Int_t label = TMath::Abs(t->GetLabel());
	      //AliAODMCParticle *AODmcTrack = (AliAODMCParticle*) fArrayMC->At(label);
	      if (AODmcTrack)
		{
		  if(TMath::Abs(AODmcTrack->GetPdgCode()) == 11) continue;
		}
	      
	      //Particle 1                                                                                                                  
	      if (t->Charge() > 0)
		{
		  
		  iPhi   = int( phi/_width_phi_1);

		  if (iPhi<0 || iPhi>=_nBins_phi_1 )
		    {
		      AliWarning("AliTaskDptCorrMC::analyze() iPhi<0 || iPhi>=_nBins_phi_1");
		      return;
		    }

		  iEta    = int((eta-_min_eta_1)/_width_eta_1);
		  if (iEta<0 || iEta>=_nBins_eta_1)
		    {
		      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iEta: %d", iEta));
		      continue;
		    }
		  iPt     = int((pt -_min_pt_1 )/_width_pt_1 );
		  if (iPt<0  || iPt >=_nBins_pt_1)
		    {
		      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iPt: %d",iPt));
		      continue;
		    }
		  iEtaPhi = iEta*_nBins_phi_1+iPhi;
		  iZEtaPhiPt = iVertexP1 + iEtaPhi*_nBins_pt_1 + iPt;

		  if (_correctionWeight_1)
		    corr = _correctionWeight_1[iZEtaPhiPt];
		  else
		    corr = 1;
		  if (iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_1)
		    {
		      AliWarning("AliTaskDptCorrMC::analyze(AliceEvent * event) iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_1");
		      continue;
		    }


		  if (_singlesOnly)
		    {

		      __n1_1_vsPt[iPt]               += corr;          //cout << "step 15" << endl;                                           

		      __n1_1_vsZEtaPhiPt[iZEtaPhiPt] += corr;       //cout << "step 12" << endl;                                              

		    }

		  else
		    {
		      corrPt                      = corr*pt;
		      _id_1[k1]                   = iTrack;
		      _charge_1[k1]               = charge;
		      _iEtaPhi_1[k1]              = iEtaPhi;
		      _iPt_1[k1]                  = iPt;
		      _pt_1[k1]                   = pt;
		      _px_1[k1]                   = px;
		      _py_1[k1]                   = py;
		      _pz_1[k1]                   = pz;
		      _correction_1[k1]           = corr;
		      __n1_1                      += corr;
		      __n1_1_vsEtaPhi[iEtaPhi]    += corr;
		      __s1pt_1                    += corrPt;
		      __s1pt_1_vsEtaPhi[iEtaPhi]  += corrPt;
		      __n1Nw_1                    += 1;
		      __s1ptNw_1                  += pt;
		      ++k1;
		      if (k1>=arraySize)
			{
			  AliError(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) k1 >=arraySize; arraySize: %d",arraySize));
			  return;
			}
		    }
		}

	      if (t->Charge() < 0 )
		{
		  //===================================                                                                                       

		  iPhi   = int( phi/_width_phi_2);

		  if (iPhi<0 || iPhi>=_nBins_phi_2 )
		    {
		      AliWarning("AliTaskDptCorrMC::analyze() iPhi<0 || iPhi>=_nBins_phi_1");
		      return;
		    }

		  iEta    = int((eta-_min_eta_2)/_width_eta_2);
		  if (iEta<0 || iEta>=_nBins_eta_2)
		    {
		      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iEta: %d", iEta));
		      continue;
		    }
		  iPt     = int((pt -_min_pt_2 )/_width_pt_2 );
		  if (iPt<0  || iPt >=_nBins_pt_2)
		    {
		      AliWarning(Form("AliTaskDptCorrMC::analyze(AliceEvent * event) Mismatched iPt: %d",iPt));
		      continue;
		    }

		  iEtaPhi = iEta*_nBins_phi_2+iPhi;
		  iZEtaPhiPt = iVertexP2 + iEtaPhi*_nBins_pt_2 + iPt;
		  if (iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_2)
		    {
		      AliWarning("AliTaskDptCorrMC::analyze(AliceEvent * event) iZEtaPhiPt<0 || iZEtaPhiPt>=_nBins_zEtaPhiPt_2");
		      continue;
		    }


		  if (_correctionWeight_2)
		    corr = _correctionWeight_2[iZEtaPhiPt];
		  else
		    corr = 1;
		  //dpt = pt - (charge>0) ? _avgPt_vsEtaPhi_2p[iEtaPhi] : _avgPt_vsEtaPhi_2m[iEtaPhi];                                        

		  if (_singlesOnly)
		    {
		      __n1_2_vsPt[iPt]               += corr;          //cout << "step 15" << endl;                                           
		      __n1_2_vsZEtaPhiPt[iZEtaPhiPt] += corr;       //cout << "step 12" << endl;                                              
		    }

		  else
		    {
		      corrPt                      = corr*pt;
		      _id_2[k2]                   = iTrack;         //cout << "step 1" << endl;                                               
		      _charge_2[k2]               = charge;         //cout << "step 2" << endl;                                               
		      _iEtaPhi_2[k2]              = iEtaPhi;        //cout << "step 3" << endl;                                               
		      _iPt_2[k2]                  = iPt;            //cout << "step 4" << endl;                                               
		      _pt_2[k2]                   = pt;             //cout << "step 5" << endl;                                               
		      _px_2[k2]                   = px;             //cout << "step 6" << endl;                                               
		      _py_2[k2]                   = py;             //cout << "step 7" << endl;                                               
		      _pz_2[k2]                   = pz;             //cout << "step 8" << endl;                                               
		      _correction_2[k2]           = corr;           //cout << "step 9" << endl;                                               
		      __n1_2                      += corr;          //cout << "step 10" << endl;                                              
		      __s1pt_2                    += corrPt;        //cout << "step 13" << endl;                                              
		      __n1Nw_2                    += 1;
		      __n1_2_vsEtaPhi[iEtaPhi]    += corr;          //cout << "step 11" << endl;                                              
		      __s1pt_2_vsEtaPhi[iEtaPhi]  += corrPt;        //cout << "step 14" << endl;                                              
		      __s1ptNw_2                  += pt;
		      ++k2;

		      if (k2>=arraySize)
			{
			  AliWarning(Form("-W- k2 >=arraySize; arraySize: %d",arraySize));
			  return;
			}
		    }

		  //cout << "done with track" << endl;                                                                                        
		} //iTrack                                                                                                                    
	    } //data aod loop                                                                                                                 
	} //MC AOD loop ends                                                                                                                  

      //************************************************
    } //AOD events             


  //  _m0->Fill(_mult0);
  //_m1->Fill(_mult1);
  //_m2->Fill(_mult2);
  //_m3->Fill(_mult3);
  //_m4->Fill(_mult4);
  //_m5->Fill(_mult5);
  //_m6->Fill(_mult6);
  _vertexZ->Fill(vertexZ);

  if (_singlesOnly)
    {
      // nothing to do here.                                                                                                                
    }
  else
    {
      if (_sameFilter)
        {
          _n1_1_vsM->Fill(centrality,      __n1_1);
          _s1pt_1_vsM->Fill(centrality,    __s1pt_1);
          _n1Nw_1_vsM->Fill(centrality,    __n1Nw_1);
          _s1ptNw_1_vsM->Fill(centrality,  __s1ptNw_1);
          _n1_2_vsM->Fill(centrality,      __n1_1);
          _s1pt_2_vsM->Fill(centrality,    __s1pt_1);
          _n1Nw_2_vsM->Fill(centrality,    __n1Nw_1);
          _s1ptNw_2_vsM->Fill(centrality,  __s1ptNw_1);
          // reset pair counters                                                                                                            
          __n2_12   = __s2ptpt_12   = __s2NPt_12    = __s2PtN_12    = 0;
          __n2Nw_12 = __s2ptptNw_12 = __s2NPtNw_12  = __s2PtNNw_12  = 0;
	  if (_field>0)
            {
              for (int i1=0; i1<k1; i1++)
                {
                  ////cout << "         i1:" << i1 << endl;                                                                                 
                  id_1      = _id_1[i1];           ////cout << "       id_1:" << id_1 << endl;                                              
                  q_1       = _charge_1[i1];       ////cout << "        q_1:" << q_1 << endl;                                               
                  iEtaPhi_1 = _iEtaPhi_1[i1];      ////cout << "  iEtaPhi_1:" << iEtaPhi_1 << endl;                                         
                  iPt_1     = _iPt_1[i1];          ////cout << "      iPt_1:" << iPt_1 << endl;                                             
                  corr_1    = _correction_1[i1];   ////cout << "     corr_1:" << corr_1 << endl;                                            
                  pt_1      = _pt_1[i1];           ////cout << "       pt_1:" << pt_1 << endl;                                              


		  //1 and 2                                                                                                                 
                  for (int i2=i1+1; i2<k1; i2++)
                    {
                      ////cout << "         i2:" << i2 << endl;                                                                             
                      id_2      = _id_1[i2];              ////cout << "       id_2:" << id_2 << endl;                                       
                      if (id_1!=id_2)
                        {
                          q_2       = _charge_1[i2];     ////cout << "        q_1:" << q_1 << endl;                                         
                          iEtaPhi_2 = _iEtaPhi_1[i2];    ////cout << "  iEtaPhi_1:" << iEtaPhi_1 << endl;                                   
                          iPt_2     = _iPt_1[i2];        ////cout << "      iPt_1:" << iPt_1 << endl;                                       
                          corr_2    = _correction_1[i2]; ////cout << "     corr_1:" << corr_1 << endl;                                      
                          pt_2      = _pt_1[i2];         ////cout << "       pt_1:" << pt_1 << endl;                                        

                          corr      = corr_1*corr_2;
                          if (q_2>q_1 || (q_1>0 && q_2>0 && pt_2<=pt_1) || (q_1<0 && q_2<0 && pt_2>=pt_1))
                            {
                              ij = iEtaPhi_1*_nBins_etaPhi_1 + iEtaPhi_2;   ////cout << " ij:" << ij<< endl;                                
                            }
                          else // swap particles                                                                                            
                            {
                              ij = iEtaPhi_2*_nBins_etaPhi_1 + iEtaPhi_1;   ////cout << " ij:" << ij<< endl;                                
                            }

			  __n2_12                  += corr;
                          __n2_12_vsEtaPhi[ij]     += corr;
                          ptpt                     = pt_1*pt_2;
                          __s2ptpt_12              += corr*ptpt;
                          __s2PtN_12               += corr*pt_1;
                          __s2NPt_12               += corr*pt_2;
                          __s2ptpt_12_vsEtaPhi[ij] += corr*ptpt;
                          __s2PtN_12_vsEtaPhi[ij]  += corr*pt_1;
                          __s2NPt_12_vsEtaPhi[ij]  += corr*pt_2;
                          __n2_12_vsPtPt[iPt_1*_nBins_pt_2 + iPt_2] += corr;

                          __n2Nw_12                  += 1;
                          __s2ptptNw_12              += ptpt;
                          __s2PtNNw_12               += pt_1;
                          __s2NPtNw_12               += pt_2;

                        }
                    } //i2                                                                                                                  
                } //i1                                                                                                                      
            }

	  else // field<0                                                                                                                   
            {
              for (int i1=0; i1<k1; i1++)
                {
                  ////cout << "         i1:" << i1 << endl;                                                                                 
                  id_1      = _id_1[i1];           ////cout << "       id_1:" << id_1 << endl;                                              
                  q_1       = _charge_1[i1];       ////cout << "        q_1:" << q_1 << endl;                                               
                  iEtaPhi_1 = _iEtaPhi_1[i1];      ////cout << "  iEtaPhi_1:" << iEtaPhi_1 << endl;                                         
                  iPt_1     = _iPt_1[i1];          ////cout << "      iPt_1:" << iPt_1 << endl;                                             
                  corr_1    = _correction_1[i1];   ////cout << "     corr_1:" << corr_1 << endl;                                            
                  pt_1      = _pt_1[i1];           ////cout << "       pt_1:" << pt_1 << endl;                                              

		  //1 and 2                                                                                                                 
                  for (int i2=i1+1; i2<k1; i2++)
                    {
                      ////cout << "         i2:" << i2 << endl;                                                                             
                      id_2      = _id_1[i2];              ////cout << "       id_2:" << id_2 << endl;                                       
                      if (id_1!=id_2)
                        {
                          q_2       = _charge_1[i2];     ////cout << "        q_2:" << q_2 << endl;                                         
                          iEtaPhi_2 = _iEtaPhi_1[i2];    ////cout << "  iEtaPhi_2:" << iEtaPhi_2 << endl;                                   
                          iPt_2     = _iPt_1[i2];        ////cout << "      iPt_2:" << iPt_2 << endl;                                       
                          corr_2    = _correction_1[i2]; ////cout << "     corr_2:" << corr_2 << endl;                                      
                          pt_2      = _pt_1[i2];         ////cout << "       pt_2:" << pt_2 << endl;                                        
                          corr      = corr_1*corr_2;
                          if ( q_2<q_1 || (q_1>0 && q_2>0 && pt_2>=pt_1) || (q_1<0 && q_2<0 && pt_2<=pt_1))
                            {
                              ij = iEtaPhi_1*_nBins_etaPhi_1 + iEtaPhi_2;   ////cout << " ij:" << ij<< endl;                                
                            }
                          else // swap particles                                                                                            
                            {
                              ij = iEtaPhi_2*_nBins_etaPhi_1 + iEtaPhi_1;   ////cout << " ij:" << ij<< endl;                                
                            }

			  __n2_12                  += corr;
                          __n2_12_vsEtaPhi[ij]     += corr;
                          ptpt                     = pt_1*pt_2;
                          __s2ptpt_12              += corr*ptpt;
                          __s2PtN_12               += corr*pt_1;
                          __s2NPt_12               += corr*pt_2;
                          __s2ptpt_12_vsEtaPhi[ij] += corr*ptpt;
                          __s2PtN_12_vsEtaPhi[ij]  += corr*pt_1;
                          __s2NPt_12_vsEtaPhi[ij]  += corr*pt_2;
                          __n2_12_vsPtPt[iPt_1*_nBins_pt_2 + iPt_2] += corr;

                          __n2Nw_12                  += 1;
                          __s2ptptNw_12              += ptpt;
                          __s2PtNNw_12               += pt_1;
                          __s2NPtNw_12               += pt_2;

			}
                    } //i2                                                                                                                  
                } //i1                                                                                                                      
            }
        }
      else  // filter 1 and 2 are different -- must do all particle pairs...                                                                
        {
          _n1_1_vsM->Fill(centrality,      __n1_1);
          _s1pt_1_vsM->Fill(centrality,    __s1pt_1);
          _n1Nw_1_vsM->Fill(centrality,    __n1Nw_1);
          _s1ptNw_1_vsM->Fill(centrality,  __s1ptNw_1);
          _n1_2_vsM->Fill(centrality,      __n1_2);
          _s1pt_2_vsM->Fill(centrality,    __s1pt_2);
          _n1Nw_2_vsM->Fill(centrality,    __n1Nw_2);
          _s1ptNw_2_vsM->Fill(centrality,  __s1ptNw_2);
          // reset pair counters                                  
	  __n2_12   = __s2ptpt_12   = __s2NPt_12    = __s2PtN_12    = 0;
          __n2Nw_12 = __s2ptptNw_12 = __s2NPtNw_12  = __s2PtNNw_12  = 0;
          for (int i1=0; i1<k1; i1++)
            {
              ////cout << "         i1:" << i1 << endl;                                                                                     
              id_1      = _id_1[i1];           ////cout << "       id_1:" << id_1 << endl;                                                  
              q_1       = _charge_1[i1];       ////cout << "        q_1:" << q_1 << endl;                                                   
              iEtaPhi_1 = _iEtaPhi_1[i1];      ////cout << "  iEtaPhi_1:" << iEtaPhi_1 << endl;                                             
              iPt_1     = _iPt_1[i1];          ////cout << "      iPt_1:" << iPt_1 << endl;                                                 
              corr_1    = _correction_1[i1];   ////cout << "     corr_1:" << corr_1 << endl;                                                
              pt_1      = _pt_1[i1];           ////cout << "       pt_1:" << pt_1 << endl;                                                  
              px_1      = _px_1[i1];          ////cout << "      px_1:" << px_1 << endl;                                                    
              py_1      = _py_1[i1];          ////cout << "      py_1:" << py_1 << endl;                                                    
              pz_1      = _pz_1[i1];          ////cout << "      pz_1:" << pz_1 << endl;                        
	      //1 and 2                                                                                                                     
              for (int i2=0; i2<k2; i2++)
                {
                  ////cout << "         i2:" << i2 << endl;                                                                                 
                  id_2   = _id_2[i2];              ////cout << "       id_2:" << id_2 << endl;                                              
                  if (id_1!=id_2)  // exclude auto correlation                                                                              
                    {
                      q_2       = _charge_2[i2];     ////cout << "        q_2:" << q_2 << endl;                                             
                      iEtaPhi_2 = _iEtaPhi_2[i2];    ////cout << "  iEtaPhi_2:" << iEtaPhi_2 << endl;                                       
                      iPt_2     = _iPt_2[i2];        ////cout << "      iPt_2:" << iPt_2 << endl;                                           
                      corr_2    = _correction_2[i2]; ////cout << "     corr_2:" << corr_2 << endl;                                          
                      pt_2      = _pt_2[i2];         ////cout << "       pt_2:" << pt_2 << endl;                                            
                      px_2      = _px_2[i2];          ////cout << "      px_2:" << px_2 << endl;                                            
                      py_2      = _py_2[i2];          ////cout << "      py_2:" << py_2 << endl;                                            
                      pz_2      = _pz_2[i2];          ////cout << "      pz_2:" << pz_2 << endl;                                            


		      if (_rejectPairConversion)
                        {
                          float e1Sq = massElecSq + pt_1*pt_1 + pz_1*pz_1;
                          float e2Sq = massElecSq + pt_2*pt_2 + pz_2*pz_2;
                          float mInvSq = 2*(massElecSq + sqrt(e1Sq*e2Sq) - px_1*px_2 - py_1*py_2 - pz_1*pz_2 );
                          float mInv = sqrt(mInvSq);
                          _invMass->Fill(mInv);
                        }

		      corr      = corr_1*corr_2;
                      ij        = iEtaPhi_1*_nBins_etaPhi_1 + iEtaPhi_2;   ////cout << " ij:" << ij<< endl;                                 
                      __n2_12                  += corr;
                      __n2_12_vsEtaPhi[ij]     += corr;
                      ptpt                     = pt_1*pt_2;
                      __s2ptpt_12              += corr*ptpt;
                      __s2PtN_12               += corr*pt_1;
                      __s2NPt_12               += corr*pt_2;
                      __s2ptpt_12_vsEtaPhi[ij] += corr*ptpt;
                      __s2PtN_12_vsEtaPhi[ij]  += corr*pt_1;
                      __s2NPt_12_vsEtaPhi[ij]  += corr*pt_2;
                      __n2_12_vsPtPt[iPt_1*_nBins_pt_2 + iPt_2] += corr;
                      __n2Nw_12                  += 1;
                      __s2ptptNw_12              += ptpt;
                      __s2PtNNw_12               += pt_1;
                      __s2NPtNw_12               += pt_2;

                    }
                } //i2                                                                                                                      
            } //i1                                                                                                                          
        }

      _n2_12_vsM->Fill(centrality,     __n2_12);
      _s2PtPt_12_vsM->Fill(centrality, __s2ptpt_12);
      _s2PtN_12_vsM->Fill(centrality,  __s2NPt_12);
      _s2NPt_12_vsM->Fill(centrality,  __s2PtN_12);

      _n2Nw_12_vsM->Fill(centrality,     __n2Nw_12);
      _s2PtPtNw_12_vsM->Fill(centrality, __s2ptptNw_12);
      _s2PtNNw_12_vsM->Fill(centrality,  __s2NPtNw_12);
      _s2NPtNw_12_vsM->Fill(centrality,  __s2PtNNw_12);

    }


  AliInfo("AliTaskDptCorrMC::UserExec()   -----------------Event Done ");
  PostData(0,_outputHistoList);

} //End of UserExec     




void   AliDptDptInMC::FinishTaskOutput()
{
  AliInfo("AliDptDptInMC::FinishTaskOutput() Starting.");
  Printf("= 0 ====================================================================");
  finalizeHistograms();
  AliInfo("= 1 ====================================================================");
  PostData(0,_outputHistoList);
  AliInfo("= 2 ====================================================================");
  AliInfo("AliDptDptInMC::FinishTaskOutput() Done.");
}

void   AliDptDptInMC::Terminate(Option_t* /*option*/)
{
  AliInfo("AliDptDptInMC::Terminate() Starting/Done.");
}


//Tools
//===================================================================================================
void  AliDptDptInMC::fillHistoWithArray(TH1 * h, float * array, int size)
{
  int i, i1;
  float v1, ev1, v2, ev2, sum, esum;
  for (i=0, i1=1; i<size; ++i,++i1)
    {
    v1  = array[i]; ev1 = sqrt(v1);
    v2  = h->GetBinContent(i1);
    ev2 = h->GetBinError(i1);
    sum = v1 + v2;
    esum = sqrt(ev1*ev1+ev2*ev2);
    h->SetBinContent(i1,sum);
    h->SetBinError(i1,esum);
    }
}

void  AliDptDptInMC::fillHistoWithArray(TH2 * h, float * array, int size1, int size2)
{
  int i, i1;
  int j, j1;
  float v1, ev1, v2, ev2, sum, esum;
  for (i=0, i1=1; i<size1; ++i,++i1)
    {
    for (j=0, j1=1; j<size2; ++j,++j1)
      {
      v1  = array[i*size2+j]; ev1 = sqrt(v1);
      v2  = h->GetBinContent(i1,j1);
      ev2 = h->GetBinError(i1,j1);
      sum = v1 + v2;
      esum = sqrt(ev1*ev1+ev2*ev2);
      h->SetBinContent(i1,j1,sum);
      h->SetBinError(i1,j1,esum);
      }
    }
}

void  AliDptDptInMC::fillHistoWithArray(TH3 * h, float * array, int size1, int size2, int size3)
{
  int i, i1;
  int j, j1;
  int k, k1;
  float v1, ev1, v2, ev2, sum, esum;
  int size23 = size2*size3;
  for (i=0, i1=1; i<size1; ++i,++i1)
    {
    for (j=0, j1=1; j<size2; ++j,++j1)
      {
      for (k=0, k1=1; k<size3; ++k,++k1)
        {
        v1  = array[i*size23+j*size3+k]; ev1 = sqrt(v1);
        v2  = h->GetBinContent(i1,j1,k1);
        ev2 = h->GetBinError(i1,j1,k1);
        sum = v1 + v2;
        esum = sqrt(ev1*ev1+ev2*ev2);
        h->SetBinContent(i1,j1,k1,sum);
        h->SetBinError(i1,j1,k1,esum);
        }
      }
    }
}

void  AliDptDptInMC::fillHistoWithArray(TH1 * h, double * array, int size)
{
  int i, i1;
  double v1, ev1, v2, ev2, sum, esum;
  for (i=0, i1=1; i<size; ++i,++i1)
    {
    v1  = array[i]; ev1 = sqrt(v1);
    v2  = h->GetBinContent(i1);
    ev2 = h->GetBinError(i1);
    sum = v1 + v2;
    esum = sqrt(ev1*ev1+ev2*ev2);
    h->SetBinContent(i1,sum);
    h->SetBinError(i1,esum);
    }
}

void  AliDptDptInMC::fillHistoWithArray(TH2 * h, double * array, int size1, int size2)
{
  int i, i1;
  int j, j1;
  double v1, ev1, v2, ev2, sum, esum;
  for (i=0, i1=1; i<size1; ++i,++i1)
    {
    for (j=0, j1=1; j<size2; ++j,++j1)
      {
      v1  = array[i*size2+j]; ev1 = sqrt(v1);
      v2  = h->GetBinContent(i1,j1);
      ev2 = h->GetBinError(i1,j1);
      sum = v1 + v2;
      esum = sqrt(ev1*ev1+ev2*ev2);
      h->SetBinContent(i1,j1,sum);
      h->SetBinError(i1,j1,esum);
      }
    }
}

void  AliDptDptInMC::fillHistoWithArray(TH3 * h, double * array, int size1, int size2, int size3)
{
  int i, i1;
  int j, j1;
  int k, k1;
  double v1, ev1, v2, ev2, sum, esum;
  int size23 = size2*size3;
  for (i=0, i1=1; i<size1; ++i,++i1)
    {
    for (j=0, j1=1; j<size2; ++j,++j1)
      {
      for (k=0, k1=1; k<size3; ++k,++k1)
        {
        v1  = array[i*size23+j*size3+k]; ev1 = sqrt(v1);
        v2  = h->GetBinContent(i1,j1,k1);
        ev2 = h->GetBinError(i1,j1,k1);
        sum = v1 + v2;
        esum = sqrt(ev1*ev1+ev2*ev2);
        h->SetBinContent(i1,j1,k1,sum);
        h->SetBinError(i1,j1,k1,esum);
        }
      }
    }
}

//________________________________________________________________________
double *  AliDptDptInMC::getDoubleArray(int size, double v)
{
  /// Allocate an array of type double with n values
  /// Initialize the array to the given value
  double * array = new double [size];
  for (int i=0;i<size;++i) array[i]=v;
  return array;
}

//________________________________________________________________________
float *  AliDptDptInMC::getFloatArray(int size, float v)
{
  /// Allocate an array of type float with n values
  /// Initialize the array to the given value
  float * array = new float [size];
  for (int i=0;i<size;++i) array[i]=v;
  return array;
}


//________________________________________________________________________
TH1D * AliDptDptInMC::createHisto1D(const TString &  name, const TString &  title, 
                                                      int n, double xMin, double xMax, 
                                                      const TString &  xTitle, const TString &  yTitle)
{
  //CreateHisto new 1D historgram
  AliInfo(Form("createHisto 1D histo %s   nBins: %d  xMin: %f   xMax: %f",name.Data(),n,xMin,xMax));
  TH1D * h = new TH1D(name,title,n,xMin,xMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  addToList(h);
  return h;
}


//________________________________________________________________________
TH1D * AliDptDptInMC::createHisto1D(const TString &  name, const TString &  title, 
                                                      int n, double * bins, 
                                                      const TString &  xTitle, const TString &  yTitle)
{
  AliInfo(Form("createHisto 1D histo %s   with %d non uniform nBins",name.Data(),n));
  TH1D * h = new TH1D(name,title,n,bins);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  addToList(h);
  return h;
}


//________________________________________________________________________
TH2D * AliDptDptInMC::createHisto2D(const TString &  name, const TString &  title, 
                                                      int nx, double xMin, double xMax, int ny, double yMin, double yMax, 
                                                      const TString &  xTitle, const TString &  yTitle, const TString &  zTitle)
{
  AliInfo(Form("createHisto 2D histo %s  nx: %d  xMin: %f10.4 xMax: %f10.4  ny: %d   yMin: %f10.4 yMax: %f10.4",name.Data(),nx,xMin,xMax,ny,yMin,yMax));
  TH2D * h = new TH2D(name,title,nx,xMin,xMax,ny,yMin,yMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  h->GetZaxis()->SetTitle(zTitle);
  addToList(h);
  return h;
}

//________________________________________________________________________
TH2D * AliDptDptInMC::createHisto2D(const TString &  name, const TString &  title, 
                                                      int nx, double* xbins, int ny, double yMin, double yMax, 
                                                      const TString &  xTitle, const TString &  yTitle, const TString &  zTitle)
{
  AliInfo(Form("createHisto 2D histo %s   with %d non uniform nBins",name.Data(),nx));
  TH2D * h;
  h = new TH2D(name,title,nx,xbins,ny,yMin,yMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  h->GetZaxis()->SetTitle(zTitle);
  addToList(h);
  return h;
}

//// F /////
//________________________________________________________________________
TH1F * AliDptDptInMC::createHisto1F(const TString &  name, const TString &  title, 
                                                        int n, double xMin, double xMax, 
                                                        const TString &  xTitle, const TString &  yTitle)
{
  //CreateHisto new 1D historgram
  AliInfo(Form("createHisto 1D histo %s   nBins: %d  xMin: %f   xMax: %f",name.Data(),n,xMin,xMax));
  TH1F * h = new TH1F(name,title,n,xMin,xMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  addToList(h);
  return h;
}


//________________________________________________________________________
TH1F * AliDptDptInMC::createHisto1F(const TString &  name, const TString &  title, 
                                                        int n, double * bins, 
                                                        const TString &  xTitle, const TString &  yTitle)
{
  AliInfo(Form("createHisto 1D histo %s   with %d non uniform nBins",name.Data(),n));
  TH1F * h = new TH1F(name,title,n,bins);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  addToList(h);
  return h;
}


//________________________________________________________________________
TH2F * AliDptDptInMC::createHisto2F(const TString &  name, const TString &  title, 
                                                        int nx, double xMin, double xMax, int ny, double yMin, double yMax, 
                                                        const TString &  xTitle, const TString &  yTitle, const TString &  zTitle)
{
  AliInfo(Form("createHisto 2D histo %s  nx: %d  xMin: %f10.4 xMax: %f10.4  ny: %d   yMin: %f10.4 yMax: %f10.4",name.Data(),nx,xMin,xMax,ny,yMin,yMax));
  TH2F * h = new TH2F(name,title,nx,xMin,xMax,ny,yMin,yMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  h->GetZaxis()->SetTitle(zTitle);
  addToList(h);
  return h;
}

//________________________________________________________________________
TH2F * AliDptDptInMC::createHisto2F(const TString &  name, const TString &  title, 
                                                        int nx, double* xbins, int ny, double yMin, double yMax, 
                                                        const TString &  xTitle, const TString &  yTitle, const TString &  zTitle)
{
  AliInfo(Form("createHisto 2D histo %s   with %d non uniform nBins",name.Data(),nx));
  TH2F * h;
  h = new TH2F(name,title,nx,xbins,ny,yMin,yMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  h->GetZaxis()->SetTitle(zTitle);
  addToList(h);
  return h;
}


//________________________________________________________________________
TH3F * AliDptDptInMC::createHisto3F(const TString &  name, const TString &  title, 
                                                      int nx, double xMin, double xMax, 
                                                      int ny, double yMin, double yMax, 
                                                      int nz, double zMin, double zMax, 
                                                      const TString &  xTitle, const TString &  yTitle, const TString &  zTitle)
{
  AliInfo(Form("createHisto 2D histo %s  nx: %d  xMin: %f10.4 xMax: %f10.4  ny: %d   yMin: %f10.4 yMax: %f10.4 nz: %d   zMin: %f10.4 zMax: %f10.4",name.Data(),nx,xMin,xMax,ny,yMin,yMax,nz,zMin,zMax));
  TH3F * h = new TH3F(name,title,nx,xMin,xMax,ny,yMin,yMax,nz,zMin,zMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  h->GetZaxis()->SetTitle(zTitle);
  addToList(h);
  return h;
}


//________________________________________________________________________
TProfile * AliDptDptInMC::createProfile(const TString & name, const TString & description, 
                                                            int nx,double xMin,double xMax,
                                                            const TString &  xTitle, const TString &  yTitle)
{
  AliInfo(Form("createHisto 1D profile %s   nBins: %d  xMin: %f10.4 xMax: %f10.4",name.Data(),nx,xMin,xMax));
  TProfile * h = new TProfile(name,description,nx,xMin,xMax);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  addToList(h);
  return h;  
}

//________________________________________________________________________
TProfile * AliDptDptInMC::createProfile(const TString &  name,const TString &  description, 
                                                            int nx,  double* bins,
                                                            const TString &  xTitle, const TString &  yTitle)
{
  AliInfo(Form("createHisto 1D profile %s  with %d non uniform bins",name.Data(),nx));
  TProfile * h = new TProfile(name,description,nx,bins);
  h->GetXaxis()->SetTitle(xTitle);
  h->GetYaxis()->SetTitle(yTitle);
  addToList(h);
  return h;
}


void   AliDptDptInMC::addToList(TH1 *h)
{
  if (_outputHistoList)
    {
    _outputHistoList->Add(h);
    }
  else
    AliInfo("addToList(TH1 *h) _outputHistoList is null!!!!! Should abort ship");

}



 AliDptDptInMC.cxx:1
 AliDptDptInMC.cxx:2
 AliDptDptInMC.cxx:3
 AliDptDptInMC.cxx:4
 AliDptDptInMC.cxx:5
 AliDptDptInMC.cxx:6
 AliDptDptInMC.cxx:7
 AliDptDptInMC.cxx:8
 AliDptDptInMC.cxx:9
 AliDptDptInMC.cxx:10
 AliDptDptInMC.cxx:11
 AliDptDptInMC.cxx:12
 AliDptDptInMC.cxx:13
 AliDptDptInMC.cxx:14
 AliDptDptInMC.cxx:15
 AliDptDptInMC.cxx:16
 AliDptDptInMC.cxx:17
 AliDptDptInMC.cxx:18
 AliDptDptInMC.cxx:19
 AliDptDptInMC.cxx:20
 AliDptDptInMC.cxx:21
 AliDptDptInMC.cxx:22
 AliDptDptInMC.cxx:23
 AliDptDptInMC.cxx:24
 AliDptDptInMC.cxx:25
 AliDptDptInMC.cxx:26
 AliDptDptInMC.cxx:27
 AliDptDptInMC.cxx:28
 AliDptDptInMC.cxx:29
 AliDptDptInMC.cxx:30
 AliDptDptInMC.cxx:31
 AliDptDptInMC.cxx:32
 AliDptDptInMC.cxx:33
 AliDptDptInMC.cxx:34
 AliDptDptInMC.cxx:35
 AliDptDptInMC.cxx:36
 AliDptDptInMC.cxx:37
 AliDptDptInMC.cxx:38
 AliDptDptInMC.cxx:39
 AliDptDptInMC.cxx:40
 AliDptDptInMC.cxx:41
 AliDptDptInMC.cxx:42
 AliDptDptInMC.cxx:43
 AliDptDptInMC.cxx:44
 AliDptDptInMC.cxx:45
 AliDptDptInMC.cxx:46
 AliDptDptInMC.cxx:47
 AliDptDptInMC.cxx:48
 AliDptDptInMC.cxx:49
 AliDptDptInMC.cxx:50
 AliDptDptInMC.cxx:51
 AliDptDptInMC.cxx:52
 AliDptDptInMC.cxx:53
 AliDptDptInMC.cxx:54
 AliDptDptInMC.cxx:55
 AliDptDptInMC.cxx:56
 AliDptDptInMC.cxx:57
 AliDptDptInMC.cxx:58
 AliDptDptInMC.cxx:59
 AliDptDptInMC.cxx:60
 AliDptDptInMC.cxx:61
 AliDptDptInMC.cxx:62
 AliDptDptInMC.cxx:63
 AliDptDptInMC.cxx:64
 AliDptDptInMC.cxx:65
 AliDptDptInMC.cxx:66
 AliDptDptInMC.cxx:67
 AliDptDptInMC.cxx:68
 AliDptDptInMC.cxx:69
 AliDptDptInMC.cxx:70
 AliDptDptInMC.cxx:71
 AliDptDptInMC.cxx:72
 AliDptDptInMC.cxx:73
 AliDptDptInMC.cxx:74
 AliDptDptInMC.cxx:75
 AliDptDptInMC.cxx:76
 AliDptDptInMC.cxx:77
 AliDptDptInMC.cxx:78
 AliDptDptInMC.cxx:79
 AliDptDptInMC.cxx:80
 AliDptDptInMC.cxx:81
 AliDptDptInMC.cxx:82
 AliDptDptInMC.cxx:83
 AliDptDptInMC.cxx:84
 AliDptDptInMC.cxx:85
 AliDptDptInMC.cxx:86
 AliDptDptInMC.cxx:87
 AliDptDptInMC.cxx:88
 AliDptDptInMC.cxx:89
 AliDptDptInMC.cxx:90
 AliDptDptInMC.cxx:91
 AliDptDptInMC.cxx:92
 AliDptDptInMC.cxx:93
 AliDptDptInMC.cxx:94
 AliDptDptInMC.cxx:95
 AliDptDptInMC.cxx:96
 AliDptDptInMC.cxx:97
 AliDptDptInMC.cxx:98
 AliDptDptInMC.cxx:99
 AliDptDptInMC.cxx:100
 AliDptDptInMC.cxx:101
 AliDptDptInMC.cxx:102
 AliDptDptInMC.cxx:103
 AliDptDptInMC.cxx:104
 AliDptDptInMC.cxx:105
 AliDptDptInMC.cxx:106
 AliDptDptInMC.cxx:107
 AliDptDptInMC.cxx:108
 AliDptDptInMC.cxx:109
 AliDptDptInMC.cxx:110
 AliDptDptInMC.cxx:111
 AliDptDptInMC.cxx:112
 AliDptDptInMC.cxx:113
 AliDptDptInMC.cxx:114
 AliDptDptInMC.cxx:115
 AliDptDptInMC.cxx:116
 AliDptDptInMC.cxx:117
 AliDptDptInMC.cxx:118
 AliDptDptInMC.cxx:119
 AliDptDptInMC.cxx:120
 AliDptDptInMC.cxx:121
 AliDptDptInMC.cxx:122
 AliDptDptInMC.cxx:123
 AliDptDptInMC.cxx:124
 AliDptDptInMC.cxx:125
 AliDptDptInMC.cxx:126
 AliDptDptInMC.cxx:127
 AliDptDptInMC.cxx:128
 AliDptDptInMC.cxx:129
 AliDptDptInMC.cxx:130
 AliDptDptInMC.cxx:131
 AliDptDptInMC.cxx:132
 AliDptDptInMC.cxx:133
 AliDptDptInMC.cxx:134
 AliDptDptInMC.cxx:135
 AliDptDptInMC.cxx:136
 AliDptDptInMC.cxx:137
 AliDptDptInMC.cxx:138
 AliDptDptInMC.cxx:139
 AliDptDptInMC.cxx:140
 AliDptDptInMC.cxx:141
 AliDptDptInMC.cxx:142
 AliDptDptInMC.cxx:143
 AliDptDptInMC.cxx:144
 AliDptDptInMC.cxx:145
 AliDptDptInMC.cxx:146
 AliDptDptInMC.cxx:147
 AliDptDptInMC.cxx:148
 AliDptDptInMC.cxx:149
 AliDptDptInMC.cxx:150
 AliDptDptInMC.cxx:151
 AliDptDptInMC.cxx:152
 AliDptDptInMC.cxx:153
 AliDptDptInMC.cxx:154
 AliDptDptInMC.cxx:155
 AliDptDptInMC.cxx:156
 AliDptDptInMC.cxx:157
 AliDptDptInMC.cxx:158
 AliDptDptInMC.cxx:159
 AliDptDptInMC.cxx:160
 AliDptDptInMC.cxx:161
 AliDptDptInMC.cxx:162
 AliDptDptInMC.cxx:163
 AliDptDptInMC.cxx:164
 AliDptDptInMC.cxx:165
 AliDptDptInMC.cxx:166
 AliDptDptInMC.cxx:167
 AliDptDptInMC.cxx:168
 AliDptDptInMC.cxx:169
 AliDptDptInMC.cxx:170
 AliDptDptInMC.cxx:171
 AliDptDptInMC.cxx:172
 AliDptDptInMC.cxx:173
 AliDptDptInMC.cxx:174
 AliDptDptInMC.cxx:175
 AliDptDptInMC.cxx:176
 AliDptDptInMC.cxx:177
 AliDptDptInMC.cxx:178
 AliDptDptInMC.cxx:179
 AliDptDptInMC.cxx:180
 AliDptDptInMC.cxx:181
 AliDptDptInMC.cxx:182
 AliDptDptInMC.cxx:183
 AliDptDptInMC.cxx:184
 AliDptDptInMC.cxx:185
 AliDptDptInMC.cxx:186
 AliDptDptInMC.cxx:187
 AliDptDptInMC.cxx:188
 AliDptDptInMC.cxx:189
 AliDptDptInMC.cxx:190
 AliDptDptInMC.cxx:191
 AliDptDptInMC.cxx:192
 AliDptDptInMC.cxx:193
 AliDptDptInMC.cxx:194
 AliDptDptInMC.cxx:195
 AliDptDptInMC.cxx:196
 AliDptDptInMC.cxx:197
 AliDptDptInMC.cxx:198
 AliDptDptInMC.cxx:199
 AliDptDptInMC.cxx:200
 AliDptDptInMC.cxx:201
 AliDptDptInMC.cxx:202
 AliDptDptInMC.cxx:203
 AliDptDptInMC.cxx:204
 AliDptDptInMC.cxx:205
 AliDptDptInMC.cxx:206
 AliDptDptInMC.cxx:207
 AliDptDptInMC.cxx:208
 AliDptDptInMC.cxx:209
 AliDptDptInMC.cxx:210
 AliDptDptInMC.cxx:211
 AliDptDptInMC.cxx:212
 AliDptDptInMC.cxx:213
 AliDptDptInMC.cxx:214
 AliDptDptInMC.cxx:215
 AliDptDptInMC.cxx:216
 AliDptDptInMC.cxx:217
 AliDptDptInMC.cxx:218
 AliDptDptInMC.cxx:219
 AliDptDptInMC.cxx:220
 AliDptDptInMC.cxx:221
 AliDptDptInMC.cxx:222
 AliDptDptInMC.cxx:223
 AliDptDptInMC.cxx:224
 AliDptDptInMC.cxx:225
 AliDptDptInMC.cxx:226
 AliDptDptInMC.cxx:227
 AliDptDptInMC.cxx:228
 AliDptDptInMC.cxx:229
 AliDptDptInMC.cxx:230
 AliDptDptInMC.cxx:231
 AliDptDptInMC.cxx:232
 AliDptDptInMC.cxx:233
 AliDptDptInMC.cxx:234
 AliDptDptInMC.cxx:235
 AliDptDptInMC.cxx:236
 AliDptDptInMC.cxx:237
 AliDptDptInMC.cxx:238
 AliDptDptInMC.cxx:239
 AliDptDptInMC.cxx:240
 AliDptDptInMC.cxx:241
 AliDptDptInMC.cxx:242
 AliDptDptInMC.cxx:243
 AliDptDptInMC.cxx:244
 AliDptDptInMC.cxx:245
 AliDptDptInMC.cxx:246
 AliDptDptInMC.cxx:247
 AliDptDptInMC.cxx:248
 AliDptDptInMC.cxx:249
 AliDptDptInMC.cxx:250
 AliDptDptInMC.cxx:251
 AliDptDptInMC.cxx:252
 AliDptDptInMC.cxx:253
 AliDptDptInMC.cxx:254
 AliDptDptInMC.cxx:255
 AliDptDptInMC.cxx:256
 AliDptDptInMC.cxx:257
 AliDptDptInMC.cxx:258
 AliDptDptInMC.cxx:259
 AliDptDptInMC.cxx:260
 AliDptDptInMC.cxx:261
 AliDptDptInMC.cxx:262
 AliDptDptInMC.cxx:263
 AliDptDptInMC.cxx:264
 AliDptDptInMC.cxx:265
 AliDptDptInMC.cxx:266
 AliDptDptInMC.cxx:267
 AliDptDptInMC.cxx:268
 AliDptDptInMC.cxx:269
 AliDptDptInMC.cxx:270
 AliDptDptInMC.cxx:271
 AliDptDptInMC.cxx:272
 AliDptDptInMC.cxx:273
 AliDptDptInMC.cxx:274
 AliDptDptInMC.cxx:275
 AliDptDptInMC.cxx:276
 AliDptDptInMC.cxx:277
 AliDptDptInMC.cxx:278
 AliDptDptInMC.cxx:279
 AliDptDptInMC.cxx:280
 AliDptDptInMC.cxx:281
 AliDptDptInMC.cxx:282
 AliDptDptInMC.cxx:283
 AliDptDptInMC.cxx:284
 AliDptDptInMC.cxx:285
 AliDptDptInMC.cxx:286
 AliDptDptInMC.cxx:287
 AliDptDptInMC.cxx:288
 AliDptDptInMC.cxx:289
 AliDptDptInMC.cxx:290
 AliDptDptInMC.cxx:291
 AliDptDptInMC.cxx:292
 AliDptDptInMC.cxx:293
 AliDptDptInMC.cxx:294
 AliDptDptInMC.cxx:295
 AliDptDptInMC.cxx:296
 AliDptDptInMC.cxx:297
 AliDptDptInMC.cxx:298
 AliDptDptInMC.cxx:299
 AliDptDptInMC.cxx:300
 AliDptDptInMC.cxx:301
 AliDptDptInMC.cxx:302
 AliDptDptInMC.cxx:303
 AliDptDptInMC.cxx:304
 AliDptDptInMC.cxx:305
 AliDptDptInMC.cxx:306
 AliDptDptInMC.cxx:307
 AliDptDptInMC.cxx:308
 AliDptDptInMC.cxx:309
 AliDptDptInMC.cxx:310
 AliDptDptInMC.cxx:311
 AliDptDptInMC.cxx:312
 AliDptDptInMC.cxx:313
 AliDptDptInMC.cxx:314
 AliDptDptInMC.cxx:315
 AliDptDptInMC.cxx:316
 AliDptDptInMC.cxx:317
 AliDptDptInMC.cxx:318
 AliDptDptInMC.cxx:319
 AliDptDptInMC.cxx:320
 AliDptDptInMC.cxx:321
 AliDptDptInMC.cxx:322
 AliDptDptInMC.cxx:323
 AliDptDptInMC.cxx:324
 AliDptDptInMC.cxx:325
 AliDptDptInMC.cxx:326
 AliDptDptInMC.cxx:327
 AliDptDptInMC.cxx:328
 AliDptDptInMC.cxx:329
 AliDptDptInMC.cxx:330
 AliDptDptInMC.cxx:331
 AliDptDptInMC.cxx:332
 AliDptDptInMC.cxx:333
 AliDptDptInMC.cxx:334
 AliDptDptInMC.cxx:335
 AliDptDptInMC.cxx:336
 AliDptDptInMC.cxx:337
 AliDptDptInMC.cxx:338
 AliDptDptInMC.cxx:339
 AliDptDptInMC.cxx:340
 AliDptDptInMC.cxx:341
 AliDptDptInMC.cxx:342
 AliDptDptInMC.cxx:343
 AliDptDptInMC.cxx:344
 AliDptDptInMC.cxx:345
 AliDptDptInMC.cxx:346
 AliDptDptInMC.cxx:347
 AliDptDptInMC.cxx:348
 AliDptDptInMC.cxx:349
 AliDptDptInMC.cxx:350
 AliDptDptInMC.cxx:351
 AliDptDptInMC.cxx:352
 AliDptDptInMC.cxx:353
 AliDptDptInMC.cxx:354
 AliDptDptInMC.cxx:355
 AliDptDptInMC.cxx:356
 AliDptDptInMC.cxx:357
 AliDptDptInMC.cxx:358
 AliDptDptInMC.cxx:359
 AliDptDptInMC.cxx:360
 AliDptDptInMC.cxx:361
 AliDptDptInMC.cxx:362
 AliDptDptInMC.cxx:363
 AliDptDptInMC.cxx:364
 AliDptDptInMC.cxx:365
 AliDptDptInMC.cxx:366
 AliDptDptInMC.cxx:367
 AliDptDptInMC.cxx:368
 AliDptDptInMC.cxx:369
 AliDptDptInMC.cxx:370
 AliDptDptInMC.cxx:371
 AliDptDptInMC.cxx:372
 AliDptDptInMC.cxx:373
 AliDptDptInMC.cxx:374
 AliDptDptInMC.cxx:375
 AliDptDptInMC.cxx:376
 AliDptDptInMC.cxx:377
 AliDptDptInMC.cxx:378
 AliDptDptInMC.cxx:379
 AliDptDptInMC.cxx:380
 AliDptDptInMC.cxx:381
 AliDptDptInMC.cxx:382
 AliDptDptInMC.cxx:383
 AliDptDptInMC.cxx:384
 AliDptDptInMC.cxx:385
 AliDptDptInMC.cxx:386
 AliDptDptInMC.cxx:387
 AliDptDptInMC.cxx:388
 AliDptDptInMC.cxx:389
 AliDptDptInMC.cxx:390
 AliDptDptInMC.cxx:391
 AliDptDptInMC.cxx:392
 AliDptDptInMC.cxx:393
 AliDptDptInMC.cxx:394
 AliDptDptInMC.cxx:395
 AliDptDptInMC.cxx:396
 AliDptDptInMC.cxx:397
 AliDptDptInMC.cxx:398
 AliDptDptInMC.cxx:399
 AliDptDptInMC.cxx:400
 AliDptDptInMC.cxx:401
 AliDptDptInMC.cxx:402
 AliDptDptInMC.cxx:403
 AliDptDptInMC.cxx:404
 AliDptDptInMC.cxx:405
 AliDptDptInMC.cxx:406
 AliDptDptInMC.cxx:407
 AliDptDptInMC.cxx:408
 AliDptDptInMC.cxx:409
 AliDptDptInMC.cxx:410
 AliDptDptInMC.cxx:411
 AliDptDptInMC.cxx:412
 AliDptDptInMC.cxx:413
 AliDptDptInMC.cxx:414
 AliDptDptInMC.cxx:415
 AliDptDptInMC.cxx:416
 AliDptDptInMC.cxx:417
 AliDptDptInMC.cxx:418
 AliDptDptInMC.cxx:419
 AliDptDptInMC.cxx:420
 AliDptDptInMC.cxx:421
 AliDptDptInMC.cxx:422
 AliDptDptInMC.cxx:423
 AliDptDptInMC.cxx:424
 AliDptDptInMC.cxx:425
 AliDptDptInMC.cxx:426
 AliDptDptInMC.cxx:427
 AliDptDptInMC.cxx:428
 AliDptDptInMC.cxx:429
 AliDptDptInMC.cxx:430
 AliDptDptInMC.cxx:431
 AliDptDptInMC.cxx:432
 AliDptDptInMC.cxx:433
 AliDptDptInMC.cxx:434
 AliDptDptInMC.cxx:435
 AliDptDptInMC.cxx:436
 AliDptDptInMC.cxx:437
 AliDptDptInMC.cxx:438
 AliDptDptInMC.cxx:439
 AliDptDptInMC.cxx:440
 AliDptDptInMC.cxx:441
 AliDptDptInMC.cxx:442
 AliDptDptInMC.cxx:443
 AliDptDptInMC.cxx:444
 AliDptDptInMC.cxx:445
 AliDptDptInMC.cxx:446
 AliDptDptInMC.cxx:447
 AliDptDptInMC.cxx:448
 AliDptDptInMC.cxx:449
 AliDptDptInMC.cxx:450
 AliDptDptInMC.cxx:451
 AliDptDptInMC.cxx:452
 AliDptDptInMC.cxx:453
 AliDptDptInMC.cxx:454
 AliDptDptInMC.cxx:455
 AliDptDptInMC.cxx:456
 AliDptDptInMC.cxx:457
 AliDptDptInMC.cxx:458
 AliDptDptInMC.cxx:459
 AliDptDptInMC.cxx:460
 AliDptDptInMC.cxx:461
 AliDptDptInMC.cxx:462
 AliDptDptInMC.cxx:463
 AliDptDptInMC.cxx:464
 AliDptDptInMC.cxx:465
 AliDptDptInMC.cxx:466
 AliDptDptInMC.cxx:467
 AliDptDptInMC.cxx:468
 AliDptDptInMC.cxx:469
 AliDptDptInMC.cxx:470
 AliDptDptInMC.cxx:471
 AliDptDptInMC.cxx:472
 AliDptDptInMC.cxx:473
 AliDptDptInMC.cxx:474
 AliDptDptInMC.cxx:475
 AliDptDptInMC.cxx:476
 AliDptDptInMC.cxx:477
 AliDptDptInMC.cxx:478
 AliDptDptInMC.cxx:479
 AliDptDptInMC.cxx:480
 AliDptDptInMC.cxx:481
 AliDptDptInMC.cxx:482
 AliDptDptInMC.cxx:483
 AliDptDptInMC.cxx:484
 AliDptDptInMC.cxx:485
 AliDptDptInMC.cxx:486
 AliDptDptInMC.cxx:487
 AliDptDptInMC.cxx:488
 AliDptDptInMC.cxx:489
 AliDptDptInMC.cxx:490
 AliDptDptInMC.cxx:491
 AliDptDptInMC.cxx:492
 AliDptDptInMC.cxx:493
 AliDptDptInMC.cxx:494
 AliDptDptInMC.cxx:495
 AliDptDptInMC.cxx:496
 AliDptDptInMC.cxx:497
 AliDptDptInMC.cxx:498
 AliDptDptInMC.cxx:499
 AliDptDptInMC.cxx:500
 AliDptDptInMC.cxx:501
 AliDptDptInMC.cxx:502
 AliDptDptInMC.cxx:503
 AliDptDptInMC.cxx:504
 AliDptDptInMC.cxx:505
 AliDptDptInMC.cxx:506
 AliDptDptInMC.cxx:507
 AliDptDptInMC.cxx:508
 AliDptDptInMC.cxx:509
 AliDptDptInMC.cxx:510
 AliDptDptInMC.cxx:511
 AliDptDptInMC.cxx:512
 AliDptDptInMC.cxx:513
 AliDptDptInMC.cxx:514
 AliDptDptInMC.cxx:515
 AliDptDptInMC.cxx:516
 AliDptDptInMC.cxx:517
 AliDptDptInMC.cxx:518
 AliDptDptInMC.cxx:519
 AliDptDptInMC.cxx:520
 AliDptDptInMC.cxx:521
 AliDptDptInMC.cxx:522
 AliDptDptInMC.cxx:523
 AliDptDptInMC.cxx:524
 AliDptDptInMC.cxx:525
 AliDptDptInMC.cxx:526
 AliDptDptInMC.cxx:527
 AliDptDptInMC.cxx:528
 AliDptDptInMC.cxx:529
 AliDptDptInMC.cxx:530
 AliDptDptInMC.cxx:531
 AliDptDptInMC.cxx:532
 AliDptDptInMC.cxx:533
 AliDptDptInMC.cxx:534
 AliDptDptInMC.cxx:535
 AliDptDptInMC.cxx:536
 AliDptDptInMC.cxx:537
 AliDptDptInMC.cxx:538
 AliDptDptInMC.cxx:539
 AliDptDptInMC.cxx:540
 AliDptDptInMC.cxx:541
 AliDptDptInMC.cxx:542
 AliDptDptInMC.cxx:543
 AliDptDptInMC.cxx:544
 AliDptDptInMC.cxx:545
 AliDptDptInMC.cxx:546
 AliDptDptInMC.cxx:547
 AliDptDptInMC.cxx:548
 AliDptDptInMC.cxx:549
 AliDptDptInMC.cxx:550
 AliDptDptInMC.cxx:551
 AliDptDptInMC.cxx:552
 AliDptDptInMC.cxx:553
 AliDptDptInMC.cxx:554
 AliDptDptInMC.cxx:555
 AliDptDptInMC.cxx:556
 AliDptDptInMC.cxx:557
 AliDptDptInMC.cxx:558
 AliDptDptInMC.cxx:559
 AliDptDptInMC.cxx:560
 AliDptDptInMC.cxx:561
 AliDptDptInMC.cxx:562
 AliDptDptInMC.cxx:563
 AliDptDptInMC.cxx:564
 AliDptDptInMC.cxx:565
 AliDptDptInMC.cxx:566
 AliDptDptInMC.cxx:567
 AliDptDptInMC.cxx:568
 AliDptDptInMC.cxx:569
 AliDptDptInMC.cxx:570
 AliDptDptInMC.cxx:571
 AliDptDptInMC.cxx:572
 AliDptDptInMC.cxx:573
 AliDptDptInMC.cxx:574
 AliDptDptInMC.cxx:575
 AliDptDptInMC.cxx:576
 AliDptDptInMC.cxx:577
 AliDptDptInMC.cxx:578
 AliDptDptInMC.cxx:579
 AliDptDptInMC.cxx:580
 AliDptDptInMC.cxx:581
 AliDptDptInMC.cxx:582
 AliDptDptInMC.cxx:583
 AliDptDptInMC.cxx:584
 AliDptDptInMC.cxx:585
 AliDptDptInMC.cxx:586
 AliDptDptInMC.cxx:587
 AliDptDptInMC.cxx:588
 AliDptDptInMC.cxx:589
 AliDptDptInMC.cxx:590
 AliDptDptInMC.cxx:591
 AliDptDptInMC.cxx:592
 AliDptDptInMC.cxx:593
 AliDptDptInMC.cxx:594
 AliDptDptInMC.cxx:595
 AliDptDptInMC.cxx:596
 AliDptDptInMC.cxx:597
 AliDptDptInMC.cxx:598
 AliDptDptInMC.cxx:599
 AliDptDptInMC.cxx:600
 AliDptDptInMC.cxx:601
 AliDptDptInMC.cxx:602
 AliDptDptInMC.cxx:603
 AliDptDptInMC.cxx:604
 AliDptDptInMC.cxx:605
 AliDptDptInMC.cxx:606
 AliDptDptInMC.cxx:607
 AliDptDptInMC.cxx:608
 AliDptDptInMC.cxx:609
 AliDptDptInMC.cxx:610
 AliDptDptInMC.cxx:611
 AliDptDptInMC.cxx:612
 AliDptDptInMC.cxx:613
 AliDptDptInMC.cxx:614
 AliDptDptInMC.cxx:615
 AliDptDptInMC.cxx:616
 AliDptDptInMC.cxx:617
 AliDptDptInMC.cxx:618
 AliDptDptInMC.cxx:619
 AliDptDptInMC.cxx:620
 AliDptDptInMC.cxx:621
 AliDptDptInMC.cxx:622
 AliDptDptInMC.cxx:623
 AliDptDptInMC.cxx:624
 AliDptDptInMC.cxx:625
 AliDptDptInMC.cxx:626
 AliDptDptInMC.cxx:627
 AliDptDptInMC.cxx:628
 AliDptDptInMC.cxx:629
 AliDptDptInMC.cxx:630
 AliDptDptInMC.cxx:631
 AliDptDptInMC.cxx:632
 AliDptDptInMC.cxx:633
 AliDptDptInMC.cxx:634
 AliDptDptInMC.cxx:635
 AliDptDptInMC.cxx:636
 AliDptDptInMC.cxx:637
 AliDptDptInMC.cxx:638
 AliDptDptInMC.cxx:639
 AliDptDptInMC.cxx:640
 AliDptDptInMC.cxx:641
 AliDptDptInMC.cxx:642
 AliDptDptInMC.cxx:643
 AliDptDptInMC.cxx:644
 AliDptDptInMC.cxx:645
 AliDptDptInMC.cxx:646
 AliDptDptInMC.cxx:647
 AliDptDptInMC.cxx:648
 AliDptDptInMC.cxx:649
 AliDptDptInMC.cxx:650
 AliDptDptInMC.cxx:651
 AliDptDptInMC.cxx:652
 AliDptDptInMC.cxx:653
 AliDptDptInMC.cxx:654
 AliDptDptInMC.cxx:655
 AliDptDptInMC.cxx:656
 AliDptDptInMC.cxx:657
 AliDptDptInMC.cxx:658
 AliDptDptInMC.cxx:659
 AliDptDptInMC.cxx:660
 AliDptDptInMC.cxx:661
 AliDptDptInMC.cxx:662
 AliDptDptInMC.cxx:663
 AliDptDptInMC.cxx:664
 AliDptDptInMC.cxx:665
 AliDptDptInMC.cxx:666
 AliDptDptInMC.cxx:667
 AliDptDptInMC.cxx:668
 AliDptDptInMC.cxx:669
 AliDptDptInMC.cxx:670
 AliDptDptInMC.cxx:671
 AliDptDptInMC.cxx:672
 AliDptDptInMC.cxx:673
 AliDptDptInMC.cxx:674
 AliDptDptInMC.cxx:675
 AliDptDptInMC.cxx:676
 AliDptDptInMC.cxx:677
 AliDptDptInMC.cxx:678
 AliDptDptInMC.cxx:679
 AliDptDptInMC.cxx:680
 AliDptDptInMC.cxx:681
 AliDptDptInMC.cxx:682
 AliDptDptInMC.cxx:683
 AliDptDptInMC.cxx:684
 AliDptDptInMC.cxx:685
 AliDptDptInMC.cxx:686
 AliDptDptInMC.cxx:687
 AliDptDptInMC.cxx:688
 AliDptDptInMC.cxx:689
 AliDptDptInMC.cxx:690
 AliDptDptInMC.cxx:691
 AliDptDptInMC.cxx:692
 AliDptDptInMC.cxx:693
 AliDptDptInMC.cxx:694
 AliDptDptInMC.cxx:695
 AliDptDptInMC.cxx:696
 AliDptDptInMC.cxx:697
 AliDptDptInMC.cxx:698
 AliDptDptInMC.cxx:699
 AliDptDptInMC.cxx:700
 AliDptDptInMC.cxx:701
 AliDptDptInMC.cxx:702
 AliDptDptInMC.cxx:703
 AliDptDptInMC.cxx:704
 AliDptDptInMC.cxx:705
 AliDptDptInMC.cxx:706
 AliDptDptInMC.cxx:707
 AliDptDptInMC.cxx:708
 AliDptDptInMC.cxx:709
 AliDptDptInMC.cxx:710
 AliDptDptInMC.cxx:711
 AliDptDptInMC.cxx:712
 AliDptDptInMC.cxx:713
 AliDptDptInMC.cxx:714
 AliDptDptInMC.cxx:715
 AliDptDptInMC.cxx:716
 AliDptDptInMC.cxx:717
 AliDptDptInMC.cxx:718
 AliDptDptInMC.cxx:719
 AliDptDptInMC.cxx:720
 AliDptDptInMC.cxx:721
 AliDptDptInMC.cxx:722
 AliDptDptInMC.cxx:723
 AliDptDptInMC.cxx:724
 AliDptDptInMC.cxx:725
 AliDptDptInMC.cxx:726
 AliDptDptInMC.cxx:727
 AliDptDptInMC.cxx:728
 AliDptDptInMC.cxx:729
 AliDptDptInMC.cxx:730
 AliDptDptInMC.cxx:731
 AliDptDptInMC.cxx:732
 AliDptDptInMC.cxx:733
 AliDptDptInMC.cxx:734
 AliDptDptInMC.cxx:735
 AliDptDptInMC.cxx:736
 AliDptDptInMC.cxx:737
 AliDptDptInMC.cxx:738
 AliDptDptInMC.cxx:739
 AliDptDptInMC.cxx:740
 AliDptDptInMC.cxx:741
 AliDptDptInMC.cxx:742
 AliDptDptInMC.cxx:743
 AliDptDptInMC.cxx:744
 AliDptDptInMC.cxx:745
 AliDptDptInMC.cxx:746
 AliDptDptInMC.cxx:747
 AliDptDptInMC.cxx:748
 AliDptDptInMC.cxx:749
 AliDptDptInMC.cxx:750
 AliDptDptInMC.cxx:751
 AliDptDptInMC.cxx:752
 AliDptDptInMC.cxx:753
 AliDptDptInMC.cxx:754
 AliDptDptInMC.cxx:755
 AliDptDptInMC.cxx:756
 AliDptDptInMC.cxx:757
 AliDptDptInMC.cxx:758
 AliDptDptInMC.cxx:759
 AliDptDptInMC.cxx:760
 AliDptDptInMC.cxx:761
 AliDptDptInMC.cxx:762
 AliDptDptInMC.cxx:763
 AliDptDptInMC.cxx:764
 AliDptDptInMC.cxx:765
 AliDptDptInMC.cxx:766
 AliDptDptInMC.cxx:767
 AliDptDptInMC.cxx:768
 AliDptDptInMC.cxx:769
 AliDptDptInMC.cxx:770
 AliDptDptInMC.cxx:771
 AliDptDptInMC.cxx:772
 AliDptDptInMC.cxx:773
 AliDptDptInMC.cxx:774
 AliDptDptInMC.cxx:775
 AliDptDptInMC.cxx:776
 AliDptDptInMC.cxx:777
 AliDptDptInMC.cxx:778
 AliDptDptInMC.cxx:779
 AliDptDptInMC.cxx:780
 AliDptDptInMC.cxx:781
 AliDptDptInMC.cxx:782
 AliDptDptInMC.cxx:783
 AliDptDptInMC.cxx:784
 AliDptDptInMC.cxx:785
 AliDptDptInMC.cxx:786
 AliDptDptInMC.cxx:787
 AliDptDptInMC.cxx:788
 AliDptDptInMC.cxx:789
 AliDptDptInMC.cxx:790
 AliDptDptInMC.cxx:791
 AliDptDptInMC.cxx:792
 AliDptDptInMC.cxx:793
 AliDptDptInMC.cxx:794
 AliDptDptInMC.cxx:795
 AliDptDptInMC.cxx:796
 AliDptDptInMC.cxx:797
 AliDptDptInMC.cxx:798
 AliDptDptInMC.cxx:799
 AliDptDptInMC.cxx:800
 AliDptDptInMC.cxx:801
 AliDptDptInMC.cxx:802
 AliDptDptInMC.cxx:803
 AliDptDptInMC.cxx:804
 AliDptDptInMC.cxx:805
 AliDptDptInMC.cxx:806
 AliDptDptInMC.cxx:807
 AliDptDptInMC.cxx:808
 AliDptDptInMC.cxx:809
 AliDptDptInMC.cxx:810
 AliDptDptInMC.cxx:811
 AliDptDptInMC.cxx:812
 AliDptDptInMC.cxx:813
 AliDptDptInMC.cxx:814
 AliDptDptInMC.cxx:815
 AliDptDptInMC.cxx:816
 AliDptDptInMC.cxx:817
 AliDptDptInMC.cxx:818
 AliDptDptInMC.cxx:819
 AliDptDptInMC.cxx:820
 AliDptDptInMC.cxx:821
 AliDptDptInMC.cxx:822
 AliDptDptInMC.cxx:823
 AliDptDptInMC.cxx:824
 AliDptDptInMC.cxx:825
 AliDptDptInMC.cxx:826
 AliDptDptInMC.cxx:827
 AliDptDptInMC.cxx:828
 AliDptDptInMC.cxx:829
 AliDptDptInMC.cxx:830
 AliDptDptInMC.cxx:831
 AliDptDptInMC.cxx:832
 AliDptDptInMC.cxx:833
 AliDptDptInMC.cxx:834
 AliDptDptInMC.cxx:835
 AliDptDptInMC.cxx:836
 AliDptDptInMC.cxx:837
 AliDptDptInMC.cxx:838
 AliDptDptInMC.cxx:839
 AliDptDptInMC.cxx:840
 AliDptDptInMC.cxx:841
 AliDptDptInMC.cxx:842
 AliDptDptInMC.cxx:843
 AliDptDptInMC.cxx:844
 AliDptDptInMC.cxx:845
 AliDptDptInMC.cxx:846
 AliDptDptInMC.cxx:847
 AliDptDptInMC.cxx:848
 AliDptDptInMC.cxx:849
 AliDptDptInMC.cxx:850
 AliDptDptInMC.cxx:851
 AliDptDptInMC.cxx:852
 AliDptDptInMC.cxx:853
 AliDptDptInMC.cxx:854
 AliDptDptInMC.cxx:855
 AliDptDptInMC.cxx:856
 AliDptDptInMC.cxx:857
 AliDptDptInMC.cxx:858
 AliDptDptInMC.cxx:859
 AliDptDptInMC.cxx:860
 AliDptDptInMC.cxx:861
 AliDptDptInMC.cxx:862
 AliDptDptInMC.cxx:863
 AliDptDptInMC.cxx:864
 AliDptDptInMC.cxx:865
 AliDptDptInMC.cxx:866
 AliDptDptInMC.cxx:867
 AliDptDptInMC.cxx:868
 AliDptDptInMC.cxx:869
 AliDptDptInMC.cxx:870
 AliDptDptInMC.cxx:871
 AliDptDptInMC.cxx:872
 AliDptDptInMC.cxx:873
 AliDptDptInMC.cxx:874
 AliDptDptInMC.cxx:875
 AliDptDptInMC.cxx:876
 AliDptDptInMC.cxx:877
 AliDptDptInMC.cxx:878
 AliDptDptInMC.cxx:879
 AliDptDptInMC.cxx:880
 AliDptDptInMC.cxx:881
 AliDptDptInMC.cxx:882
 AliDptDptInMC.cxx:883
 AliDptDptInMC.cxx:884
 AliDptDptInMC.cxx:885
 AliDptDptInMC.cxx:886
 AliDptDptInMC.cxx:887
 AliDptDptInMC.cxx:888
 AliDptDptInMC.cxx:889
 AliDptDptInMC.cxx:890
 AliDptDptInMC.cxx:891
 AliDptDptInMC.cxx:892
 AliDptDptInMC.cxx:893
 AliDptDptInMC.cxx:894
 AliDptDptInMC.cxx:895
 AliDptDptInMC.cxx:896
 AliDptDptInMC.cxx:897
 AliDptDptInMC.cxx:898
 AliDptDptInMC.cxx:899
 AliDptDptInMC.cxx:900
 AliDptDptInMC.cxx:901
 AliDptDptInMC.cxx:902
 AliDptDptInMC.cxx:903
 AliDptDptInMC.cxx:904
 AliDptDptInMC.cxx:905
 AliDptDptInMC.cxx:906
 AliDptDptInMC.cxx:907
 AliDptDptInMC.cxx:908
 AliDptDptInMC.cxx:909
 AliDptDptInMC.cxx:910
 AliDptDptInMC.cxx:911
 AliDptDptInMC.cxx:912
 AliDptDptInMC.cxx:913
 AliDptDptInMC.cxx:914
 AliDptDptInMC.cxx:915
 AliDptDptInMC.cxx:916
 AliDptDptInMC.cxx:917
 AliDptDptInMC.cxx:918
 AliDptDptInMC.cxx:919
 AliDptDptInMC.cxx:920
 AliDptDptInMC.cxx:921
 AliDptDptInMC.cxx:922
 AliDptDptInMC.cxx:923
 AliDptDptInMC.cxx:924
 AliDptDptInMC.cxx:925
 AliDptDptInMC.cxx:926
 AliDptDptInMC.cxx:927
 AliDptDptInMC.cxx:928
 AliDptDptInMC.cxx:929
 AliDptDptInMC.cxx:930
 AliDptDptInMC.cxx:931
 AliDptDptInMC.cxx:932
 AliDptDptInMC.cxx:933
 AliDptDptInMC.cxx:934
 AliDptDptInMC.cxx:935
 AliDptDptInMC.cxx:936
 AliDptDptInMC.cxx:937
 AliDptDptInMC.cxx:938
 AliDptDptInMC.cxx:939
 AliDptDptInMC.cxx:940
 AliDptDptInMC.cxx:941
 AliDptDptInMC.cxx:942
 AliDptDptInMC.cxx:943
 AliDptDptInMC.cxx:944
 AliDptDptInMC.cxx:945
 AliDptDptInMC.cxx:946
 AliDptDptInMC.cxx:947
 AliDptDptInMC.cxx:948
 AliDptDptInMC.cxx:949
 AliDptDptInMC.cxx:950
 AliDptDptInMC.cxx:951
 AliDptDptInMC.cxx:952
 AliDptDptInMC.cxx:953
 AliDptDptInMC.cxx:954
 AliDptDptInMC.cxx:955
 AliDptDptInMC.cxx:956
 AliDptDptInMC.cxx:957
 AliDptDptInMC.cxx:958
 AliDptDptInMC.cxx:959
 AliDptDptInMC.cxx:960
 AliDptDptInMC.cxx:961
 AliDptDptInMC.cxx:962
 AliDptDptInMC.cxx:963
 AliDptDptInMC.cxx:964
 AliDptDptInMC.cxx:965
 AliDptDptInMC.cxx:966
 AliDptDptInMC.cxx:967
 AliDptDptInMC.cxx:968
 AliDptDptInMC.cxx:969
 AliDptDptInMC.cxx:970
 AliDptDptInMC.cxx:971
 AliDptDptInMC.cxx:972
 AliDptDptInMC.cxx:973
 AliDptDptInMC.cxx:974
 AliDptDptInMC.cxx:975
 AliDptDptInMC.cxx:976
 AliDptDptInMC.cxx:977
 AliDptDptInMC.cxx:978
 AliDptDptInMC.cxx:979
 AliDptDptInMC.cxx:980
 AliDptDptInMC.cxx:981
 AliDptDptInMC.cxx:982
 AliDptDptInMC.cxx:983
 AliDptDptInMC.cxx:984
 AliDptDptInMC.cxx:985
 AliDptDptInMC.cxx:986
 AliDptDptInMC.cxx:987
 AliDptDptInMC.cxx:988
 AliDptDptInMC.cxx:989
 AliDptDptInMC.cxx:990
 AliDptDptInMC.cxx:991
 AliDptDptInMC.cxx:992
 AliDptDptInMC.cxx:993
 AliDptDptInMC.cxx:994
 AliDptDptInMC.cxx:995
 AliDptDptInMC.cxx:996
 AliDptDptInMC.cxx:997
 AliDptDptInMC.cxx:998
 AliDptDptInMC.cxx:999
 AliDptDptInMC.cxx:1000
 AliDptDptInMC.cxx:1001
 AliDptDptInMC.cxx:1002
 AliDptDptInMC.cxx:1003
 AliDptDptInMC.cxx:1004
 AliDptDptInMC.cxx:1005
 AliDptDptInMC.cxx:1006
 AliDptDptInMC.cxx:1007
 AliDptDptInMC.cxx:1008
 AliDptDptInMC.cxx:1009
 AliDptDptInMC.cxx:1010
 AliDptDptInMC.cxx:1011
 AliDptDptInMC.cxx:1012
 AliDptDptInMC.cxx:1013
 AliDptDptInMC.cxx:1014
 AliDptDptInMC.cxx:1015
 AliDptDptInMC.cxx:1016
 AliDptDptInMC.cxx:1017
 AliDptDptInMC.cxx:1018
 AliDptDptInMC.cxx:1019
 AliDptDptInMC.cxx:1020
 AliDptDptInMC.cxx:1021
 AliDptDptInMC.cxx:1022
 AliDptDptInMC.cxx:1023
 AliDptDptInMC.cxx:1024
 AliDptDptInMC.cxx:1025
 AliDptDptInMC.cxx:1026
 AliDptDptInMC.cxx:1027
 AliDptDptInMC.cxx:1028
 AliDptDptInMC.cxx:1029
 AliDptDptInMC.cxx:1030
 AliDptDptInMC.cxx:1031
 AliDptDptInMC.cxx:1032
 AliDptDptInMC.cxx:1033
 AliDptDptInMC.cxx:1034
 AliDptDptInMC.cxx:1035
 AliDptDptInMC.cxx:1036
 AliDptDptInMC.cxx:1037
 AliDptDptInMC.cxx:1038
 AliDptDptInMC.cxx:1039
 AliDptDptInMC.cxx:1040
 AliDptDptInMC.cxx:1041
 AliDptDptInMC.cxx:1042
 AliDptDptInMC.cxx:1043
 AliDptDptInMC.cxx:1044
 AliDptDptInMC.cxx:1045
 AliDptDptInMC.cxx:1046
 AliDptDptInMC.cxx:1047
 AliDptDptInMC.cxx:1048
 AliDptDptInMC.cxx:1049
 AliDptDptInMC.cxx:1050
 AliDptDptInMC.cxx:1051
 AliDptDptInMC.cxx:1052
 AliDptDptInMC.cxx:1053
 AliDptDptInMC.cxx:1054
 AliDptDptInMC.cxx:1055
 AliDptDptInMC.cxx:1056
 AliDptDptInMC.cxx:1057
 AliDptDptInMC.cxx:1058
 AliDptDptInMC.cxx:1059
 AliDptDptInMC.cxx:1060
 AliDptDptInMC.cxx:1061
 AliDptDptInMC.cxx:1062
 AliDptDptInMC.cxx:1063
 AliDptDptInMC.cxx:1064
 AliDptDptInMC.cxx:1065
 AliDptDptInMC.cxx:1066
 AliDptDptInMC.cxx:1067
 AliDptDptInMC.cxx:1068
 AliDptDptInMC.cxx:1069
 AliDptDptInMC.cxx:1070
 AliDptDptInMC.cxx:1071
 AliDptDptInMC.cxx:1072
 AliDptDptInMC.cxx:1073
 AliDptDptInMC.cxx:1074
 AliDptDptInMC.cxx:1075
 AliDptDptInMC.cxx:1076
 AliDptDptInMC.cxx:1077
 AliDptDptInMC.cxx:1078
 AliDptDptInMC.cxx:1079
 AliDptDptInMC.cxx:1080
 AliDptDptInMC.cxx:1081
 AliDptDptInMC.cxx:1082
 AliDptDptInMC.cxx:1083
 AliDptDptInMC.cxx:1084
 AliDptDptInMC.cxx:1085
 AliDptDptInMC.cxx:1086
 AliDptDptInMC.cxx:1087
 AliDptDptInMC.cxx:1088
 AliDptDptInMC.cxx:1089
 AliDptDptInMC.cxx:1090
 AliDptDptInMC.cxx:1091
 AliDptDptInMC.cxx:1092
 AliDptDptInMC.cxx:1093
 AliDptDptInMC.cxx:1094
 AliDptDptInMC.cxx:1095
 AliDptDptInMC.cxx:1096
 AliDptDptInMC.cxx:1097
 AliDptDptInMC.cxx:1098
 AliDptDptInMC.cxx:1099
 AliDptDptInMC.cxx:1100
 AliDptDptInMC.cxx:1101
 AliDptDptInMC.cxx:1102
 AliDptDptInMC.cxx:1103
 AliDptDptInMC.cxx:1104
 AliDptDptInMC.cxx:1105
 AliDptDptInMC.cxx:1106
 AliDptDptInMC.cxx:1107
 AliDptDptInMC.cxx:1108
 AliDptDptInMC.cxx:1109
 AliDptDptInMC.cxx:1110
 AliDptDptInMC.cxx:1111
 AliDptDptInMC.cxx:1112
 AliDptDptInMC.cxx:1113
 AliDptDptInMC.cxx:1114
 AliDptDptInMC.cxx:1115
 AliDptDptInMC.cxx:1116
 AliDptDptInMC.cxx:1117
 AliDptDptInMC.cxx:1118
 AliDptDptInMC.cxx:1119
 AliDptDptInMC.cxx:1120
 AliDptDptInMC.cxx:1121
 AliDptDptInMC.cxx:1122
 AliDptDptInMC.cxx:1123
 AliDptDptInMC.cxx:1124
 AliDptDptInMC.cxx:1125
 AliDptDptInMC.cxx:1126
 AliDptDptInMC.cxx:1127
 AliDptDptInMC.cxx:1128
 AliDptDptInMC.cxx:1129
 AliDptDptInMC.cxx:1130
 AliDptDptInMC.cxx:1131
 AliDptDptInMC.cxx:1132
 AliDptDptInMC.cxx:1133
 AliDptDptInMC.cxx:1134
 AliDptDptInMC.cxx:1135
 AliDptDptInMC.cxx:1136
 AliDptDptInMC.cxx:1137
 AliDptDptInMC.cxx:1138
 AliDptDptInMC.cxx:1139
 AliDptDptInMC.cxx:1140
 AliDptDptInMC.cxx:1141
 AliDptDptInMC.cxx:1142
 AliDptDptInMC.cxx:1143
 AliDptDptInMC.cxx:1144
 AliDptDptInMC.cxx:1145
 AliDptDptInMC.cxx:1146
 AliDptDptInMC.cxx:1147
 AliDptDptInMC.cxx:1148
 AliDptDptInMC.cxx:1149
 AliDptDptInMC.cxx:1150
 AliDptDptInMC.cxx:1151
 AliDptDptInMC.cxx:1152
 AliDptDptInMC.cxx:1153
 AliDptDptInMC.cxx:1154
 AliDptDptInMC.cxx:1155
 AliDptDptInMC.cxx:1156
 AliDptDptInMC.cxx:1157
 AliDptDptInMC.cxx:1158
 AliDptDptInMC.cxx:1159
 AliDptDptInMC.cxx:1160
 AliDptDptInMC.cxx:1161
 AliDptDptInMC.cxx:1162
 AliDptDptInMC.cxx:1163
 AliDptDptInMC.cxx:1164
 AliDptDptInMC.cxx:1165
 AliDptDptInMC.cxx:1166
 AliDptDptInMC.cxx:1167
 AliDptDptInMC.cxx:1168
 AliDptDptInMC.cxx:1169
 AliDptDptInMC.cxx:1170
 AliDptDptInMC.cxx:1171
 AliDptDptInMC.cxx:1172
 AliDptDptInMC.cxx:1173
 AliDptDptInMC.cxx:1174
 AliDptDptInMC.cxx:1175
 AliDptDptInMC.cxx:1176
 AliDptDptInMC.cxx:1177
 AliDptDptInMC.cxx:1178
 AliDptDptInMC.cxx:1179
 AliDptDptInMC.cxx:1180
 AliDptDptInMC.cxx:1181
 AliDptDptInMC.cxx:1182
 AliDptDptInMC.cxx:1183
 AliDptDptInMC.cxx:1184
 AliDptDptInMC.cxx:1185
 AliDptDptInMC.cxx:1186
 AliDptDptInMC.cxx:1187
 AliDptDptInMC.cxx:1188
 AliDptDptInMC.cxx:1189
 AliDptDptInMC.cxx:1190
 AliDptDptInMC.cxx:1191
 AliDptDptInMC.cxx:1192
 AliDptDptInMC.cxx:1193
 AliDptDptInMC.cxx:1194
 AliDptDptInMC.cxx:1195
 AliDptDptInMC.cxx:1196
 AliDptDptInMC.cxx:1197
 AliDptDptInMC.cxx:1198
 AliDptDptInMC.cxx:1199
 AliDptDptInMC.cxx:1200
 AliDptDptInMC.cxx:1201
 AliDptDptInMC.cxx:1202
 AliDptDptInMC.cxx:1203
 AliDptDptInMC.cxx:1204
 AliDptDptInMC.cxx:1205
 AliDptDptInMC.cxx:1206
 AliDptDptInMC.cxx:1207
 AliDptDptInMC.cxx:1208
 AliDptDptInMC.cxx:1209
 AliDptDptInMC.cxx:1210
 AliDptDptInMC.cxx:1211
 AliDptDptInMC.cxx:1212
 AliDptDptInMC.cxx:1213
 AliDptDptInMC.cxx:1214
 AliDptDptInMC.cxx:1215
 AliDptDptInMC.cxx:1216
 AliDptDptInMC.cxx:1217
 AliDptDptInMC.cxx:1218
 AliDptDptInMC.cxx:1219
 AliDptDptInMC.cxx:1220
 AliDptDptInMC.cxx:1221
 AliDptDptInMC.cxx:1222
 AliDptDptInMC.cxx:1223
 AliDptDptInMC.cxx:1224
 AliDptDptInMC.cxx:1225
 AliDptDptInMC.cxx:1226
 AliDptDptInMC.cxx:1227
 AliDptDptInMC.cxx:1228
 AliDptDptInMC.cxx:1229
 AliDptDptInMC.cxx:1230
 AliDptDptInMC.cxx:1231
 AliDptDptInMC.cxx:1232
 AliDptDptInMC.cxx:1233
 AliDptDptInMC.cxx:1234
 AliDptDptInMC.cxx:1235
 AliDptDptInMC.cxx:1236
 AliDptDptInMC.cxx:1237
 AliDptDptInMC.cxx:1238
 AliDptDptInMC.cxx:1239
 AliDptDptInMC.cxx:1240
 AliDptDptInMC.cxx:1241
 AliDptDptInMC.cxx:1242
 AliDptDptInMC.cxx:1243
 AliDptDptInMC.cxx:1244
 AliDptDptInMC.cxx:1245
 AliDptDptInMC.cxx:1246
 AliDptDptInMC.cxx:1247
 AliDptDptInMC.cxx:1248
 AliDptDptInMC.cxx:1249
 AliDptDptInMC.cxx:1250
 AliDptDptInMC.cxx:1251
 AliDptDptInMC.cxx:1252
 AliDptDptInMC.cxx:1253
 AliDptDptInMC.cxx:1254
 AliDptDptInMC.cxx:1255
 AliDptDptInMC.cxx:1256
 AliDptDptInMC.cxx:1257
 AliDptDptInMC.cxx:1258
 AliDptDptInMC.cxx:1259
 AliDptDptInMC.cxx:1260
 AliDptDptInMC.cxx:1261
 AliDptDptInMC.cxx:1262
 AliDptDptInMC.cxx:1263
 AliDptDptInMC.cxx:1264
 AliDptDptInMC.cxx:1265
 AliDptDptInMC.cxx:1266
 AliDptDptInMC.cxx:1267
 AliDptDptInMC.cxx:1268
 AliDptDptInMC.cxx:1269
 AliDptDptInMC.cxx:1270
 AliDptDptInMC.cxx:1271
 AliDptDptInMC.cxx:1272
 AliDptDptInMC.cxx:1273
 AliDptDptInMC.cxx:1274
 AliDptDptInMC.cxx:1275
 AliDptDptInMC.cxx:1276
 AliDptDptInMC.cxx:1277
 AliDptDptInMC.cxx:1278
 AliDptDptInMC.cxx:1279
 AliDptDptInMC.cxx:1280
 AliDptDptInMC.cxx:1281
 AliDptDptInMC.cxx:1282
 AliDptDptInMC.cxx:1283
 AliDptDptInMC.cxx:1284
 AliDptDptInMC.cxx:1285
 AliDptDptInMC.cxx:1286
 AliDptDptInMC.cxx:1287
 AliDptDptInMC.cxx:1288
 AliDptDptInMC.cxx:1289
 AliDptDptInMC.cxx:1290
 AliDptDptInMC.cxx:1291
 AliDptDptInMC.cxx:1292
 AliDptDptInMC.cxx:1293
 AliDptDptInMC.cxx:1294
 AliDptDptInMC.cxx:1295
 AliDptDptInMC.cxx:1296
 AliDptDptInMC.cxx:1297
 AliDptDptInMC.cxx:1298
 AliDptDptInMC.cxx:1299
 AliDptDptInMC.cxx:1300
 AliDptDptInMC.cxx:1301
 AliDptDptInMC.cxx:1302
 AliDptDptInMC.cxx:1303
 AliDptDptInMC.cxx:1304
 AliDptDptInMC.cxx:1305
 AliDptDptInMC.cxx:1306
 AliDptDptInMC.cxx:1307
 AliDptDptInMC.cxx:1308
 AliDptDptInMC.cxx:1309
 AliDptDptInMC.cxx:1310
 AliDptDptInMC.cxx:1311
 AliDptDptInMC.cxx:1312
 AliDptDptInMC.cxx:1313
 AliDptDptInMC.cxx:1314
 AliDptDptInMC.cxx:1315
 AliDptDptInMC.cxx:1316
 AliDptDptInMC.cxx:1317
 AliDptDptInMC.cxx:1318
 AliDptDptInMC.cxx:1319
 AliDptDptInMC.cxx:1320
 AliDptDptInMC.cxx:1321
 AliDptDptInMC.cxx:1322
 AliDptDptInMC.cxx:1323
 AliDptDptInMC.cxx:1324
 AliDptDptInMC.cxx:1325
 AliDptDptInMC.cxx:1326
 AliDptDptInMC.cxx:1327
 AliDptDptInMC.cxx:1328
 AliDptDptInMC.cxx:1329
 AliDptDptInMC.cxx:1330
 AliDptDptInMC.cxx:1331
 AliDptDptInMC.cxx:1332
 AliDptDptInMC.cxx:1333
 AliDptDptInMC.cxx:1334
 AliDptDptInMC.cxx:1335
 AliDptDptInMC.cxx:1336
 AliDptDptInMC.cxx:1337
 AliDptDptInMC.cxx:1338
 AliDptDptInMC.cxx:1339
 AliDptDptInMC.cxx:1340
 AliDptDptInMC.cxx:1341
 AliDptDptInMC.cxx:1342
 AliDptDptInMC.cxx:1343
 AliDptDptInMC.cxx:1344
 AliDptDptInMC.cxx:1345
 AliDptDptInMC.cxx:1346
 AliDptDptInMC.cxx:1347
 AliDptDptInMC.cxx:1348
 AliDptDptInMC.cxx:1349
 AliDptDptInMC.cxx:1350
 AliDptDptInMC.cxx:1351
 AliDptDptInMC.cxx:1352
 AliDptDptInMC.cxx:1353
 AliDptDptInMC.cxx:1354
 AliDptDptInMC.cxx:1355
 AliDptDptInMC.cxx:1356
 AliDptDptInMC.cxx:1357
 AliDptDptInMC.cxx:1358
 AliDptDptInMC.cxx:1359
 AliDptDptInMC.cxx:1360
 AliDptDptInMC.cxx:1361
 AliDptDptInMC.cxx:1362
 AliDptDptInMC.cxx:1363
 AliDptDptInMC.cxx:1364
 AliDptDptInMC.cxx:1365
 AliDptDptInMC.cxx:1366
 AliDptDptInMC.cxx:1367
 AliDptDptInMC.cxx:1368
 AliDptDptInMC.cxx:1369
 AliDptDptInMC.cxx:1370
 AliDptDptInMC.cxx:1371
 AliDptDptInMC.cxx:1372
 AliDptDptInMC.cxx:1373
 AliDptDptInMC.cxx:1374
 AliDptDptInMC.cxx:1375
 AliDptDptInMC.cxx:1376
 AliDptDptInMC.cxx:1377
 AliDptDptInMC.cxx:1378
 AliDptDptInMC.cxx:1379
 AliDptDptInMC.cxx:1380
 AliDptDptInMC.cxx:1381
 AliDptDptInMC.cxx:1382
 AliDptDptInMC.cxx:1383
 AliDptDptInMC.cxx:1384
 AliDptDptInMC.cxx:1385
 AliDptDptInMC.cxx:1386
 AliDptDptInMC.cxx:1387
 AliDptDptInMC.cxx:1388
 AliDptDptInMC.cxx:1389
 AliDptDptInMC.cxx:1390
 AliDptDptInMC.cxx:1391
 AliDptDptInMC.cxx:1392
 AliDptDptInMC.cxx:1393
 AliDptDptInMC.cxx:1394
 AliDptDptInMC.cxx:1395
 AliDptDptInMC.cxx:1396
 AliDptDptInMC.cxx:1397
 AliDptDptInMC.cxx:1398
 AliDptDptInMC.cxx:1399
 AliDptDptInMC.cxx:1400
 AliDptDptInMC.cxx:1401
 AliDptDptInMC.cxx:1402
 AliDptDptInMC.cxx:1403
 AliDptDptInMC.cxx:1404
 AliDptDptInMC.cxx:1405
 AliDptDptInMC.cxx:1406
 AliDptDptInMC.cxx:1407
 AliDptDptInMC.cxx:1408
 AliDptDptInMC.cxx:1409
 AliDptDptInMC.cxx:1410
 AliDptDptInMC.cxx:1411
 AliDptDptInMC.cxx:1412
 AliDptDptInMC.cxx:1413
 AliDptDptInMC.cxx:1414
 AliDptDptInMC.cxx:1415
 AliDptDptInMC.cxx:1416
 AliDptDptInMC.cxx:1417
 AliDptDptInMC.cxx:1418
 AliDptDptInMC.cxx:1419
 AliDptDptInMC.cxx:1420
 AliDptDptInMC.cxx:1421
 AliDptDptInMC.cxx:1422
 AliDptDptInMC.cxx:1423
 AliDptDptInMC.cxx:1424
 AliDptDptInMC.cxx:1425
 AliDptDptInMC.cxx:1426
 AliDptDptInMC.cxx:1427
 AliDptDptInMC.cxx:1428
 AliDptDptInMC.cxx:1429
 AliDptDptInMC.cxx:1430
 AliDptDptInMC.cxx:1431
 AliDptDptInMC.cxx:1432
 AliDptDptInMC.cxx:1433
 AliDptDptInMC.cxx:1434
 AliDptDptInMC.cxx:1435
 AliDptDptInMC.cxx:1436
 AliDptDptInMC.cxx:1437
 AliDptDptInMC.cxx:1438
 AliDptDptInMC.cxx:1439
 AliDptDptInMC.cxx:1440
 AliDptDptInMC.cxx:1441
 AliDptDptInMC.cxx:1442
 AliDptDptInMC.cxx:1443
 AliDptDptInMC.cxx:1444
 AliDptDptInMC.cxx:1445
 AliDptDptInMC.cxx:1446
 AliDptDptInMC.cxx:1447
 AliDptDptInMC.cxx:1448
 AliDptDptInMC.cxx:1449
 AliDptDptInMC.cxx:1450
 AliDptDptInMC.cxx:1451
 AliDptDptInMC.cxx:1452
 AliDptDptInMC.cxx:1453
 AliDptDptInMC.cxx:1454
 AliDptDptInMC.cxx:1455
 AliDptDptInMC.cxx:1456
 AliDptDptInMC.cxx:1457
 AliDptDptInMC.cxx:1458
 AliDptDptInMC.cxx:1459
 AliDptDptInMC.cxx:1460
 AliDptDptInMC.cxx:1461
 AliDptDptInMC.cxx:1462
 AliDptDptInMC.cxx:1463
 AliDptDptInMC.cxx:1464
 AliDptDptInMC.cxx:1465
 AliDptDptInMC.cxx:1466
 AliDptDptInMC.cxx:1467
 AliDptDptInMC.cxx:1468
 AliDptDptInMC.cxx:1469
 AliDptDptInMC.cxx:1470
 AliDptDptInMC.cxx:1471
 AliDptDptInMC.cxx:1472
 AliDptDptInMC.cxx:1473
 AliDptDptInMC.cxx:1474
 AliDptDptInMC.cxx:1475
 AliDptDptInMC.cxx:1476
 AliDptDptInMC.cxx:1477
 AliDptDptInMC.cxx:1478
 AliDptDptInMC.cxx:1479
 AliDptDptInMC.cxx:1480
 AliDptDptInMC.cxx:1481
 AliDptDptInMC.cxx:1482
 AliDptDptInMC.cxx:1483
 AliDptDptInMC.cxx:1484
 AliDptDptInMC.cxx:1485
 AliDptDptInMC.cxx:1486
 AliDptDptInMC.cxx:1487
 AliDptDptInMC.cxx:1488
 AliDptDptInMC.cxx:1489
 AliDptDptInMC.cxx:1490
 AliDptDptInMC.cxx:1491
 AliDptDptInMC.cxx:1492
 AliDptDptInMC.cxx:1493
 AliDptDptInMC.cxx:1494
 AliDptDptInMC.cxx:1495
 AliDptDptInMC.cxx:1496
 AliDptDptInMC.cxx:1497
 AliDptDptInMC.cxx:1498
 AliDptDptInMC.cxx:1499
 AliDptDptInMC.cxx:1500
 AliDptDptInMC.cxx:1501
 AliDptDptInMC.cxx:1502
 AliDptDptInMC.cxx:1503
 AliDptDptInMC.cxx:1504
 AliDptDptInMC.cxx:1505
 AliDptDptInMC.cxx:1506
 AliDptDptInMC.cxx:1507
 AliDptDptInMC.cxx:1508
 AliDptDptInMC.cxx:1509
 AliDptDptInMC.cxx:1510
 AliDptDptInMC.cxx:1511
 AliDptDptInMC.cxx:1512
 AliDptDptInMC.cxx:1513
 AliDptDptInMC.cxx:1514
 AliDptDptInMC.cxx:1515
 AliDptDptInMC.cxx:1516
 AliDptDptInMC.cxx:1517
 AliDptDptInMC.cxx:1518
 AliDptDptInMC.cxx:1519
 AliDptDptInMC.cxx:1520
 AliDptDptInMC.cxx:1521
 AliDptDptInMC.cxx:1522
 AliDptDptInMC.cxx:1523
 AliDptDptInMC.cxx:1524
 AliDptDptInMC.cxx:1525
 AliDptDptInMC.cxx:1526
 AliDptDptInMC.cxx:1527
 AliDptDptInMC.cxx:1528
 AliDptDptInMC.cxx:1529
 AliDptDptInMC.cxx:1530
 AliDptDptInMC.cxx:1531
 AliDptDptInMC.cxx:1532
 AliDptDptInMC.cxx:1533
 AliDptDptInMC.cxx:1534
 AliDptDptInMC.cxx:1535
 AliDptDptInMC.cxx:1536
 AliDptDptInMC.cxx:1537
 AliDptDptInMC.cxx:1538
 AliDptDptInMC.cxx:1539
 AliDptDptInMC.cxx:1540
 AliDptDptInMC.cxx:1541
 AliDptDptInMC.cxx:1542
 AliDptDptInMC.cxx:1543
 AliDptDptInMC.cxx:1544
 AliDptDptInMC.cxx:1545
 AliDptDptInMC.cxx:1546
 AliDptDptInMC.cxx:1547
 AliDptDptInMC.cxx:1548
 AliDptDptInMC.cxx:1549
 AliDptDptInMC.cxx:1550
 AliDptDptInMC.cxx:1551
 AliDptDptInMC.cxx:1552
 AliDptDptInMC.cxx:1553
 AliDptDptInMC.cxx:1554
 AliDptDptInMC.cxx:1555
 AliDptDptInMC.cxx:1556
 AliDptDptInMC.cxx:1557
 AliDptDptInMC.cxx:1558
 AliDptDptInMC.cxx:1559
 AliDptDptInMC.cxx:1560
 AliDptDptInMC.cxx:1561
 AliDptDptInMC.cxx:1562
 AliDptDptInMC.cxx:1563
 AliDptDptInMC.cxx:1564
 AliDptDptInMC.cxx:1565
 AliDptDptInMC.cxx:1566
 AliDptDptInMC.cxx:1567
 AliDptDptInMC.cxx:1568
 AliDptDptInMC.cxx:1569
 AliDptDptInMC.cxx:1570
 AliDptDptInMC.cxx:1571
 AliDptDptInMC.cxx:1572
 AliDptDptInMC.cxx:1573
 AliDptDptInMC.cxx:1574
 AliDptDptInMC.cxx:1575
 AliDptDptInMC.cxx:1576
 AliDptDptInMC.cxx:1577
 AliDptDptInMC.cxx:1578
 AliDptDptInMC.cxx:1579
 AliDptDptInMC.cxx:1580
 AliDptDptInMC.cxx:1581
 AliDptDptInMC.cxx:1582
 AliDptDptInMC.cxx:1583
 AliDptDptInMC.cxx:1584
 AliDptDptInMC.cxx:1585
 AliDptDptInMC.cxx:1586
 AliDptDptInMC.cxx:1587
 AliDptDptInMC.cxx:1588
 AliDptDptInMC.cxx:1589
 AliDptDptInMC.cxx:1590
 AliDptDptInMC.cxx:1591
 AliDptDptInMC.cxx:1592
 AliDptDptInMC.cxx:1593
 AliDptDptInMC.cxx:1594
 AliDptDptInMC.cxx:1595
 AliDptDptInMC.cxx:1596
 AliDptDptInMC.cxx:1597
 AliDptDptInMC.cxx:1598
 AliDptDptInMC.cxx:1599
 AliDptDptInMC.cxx:1600
 AliDptDptInMC.cxx:1601
 AliDptDptInMC.cxx:1602
 AliDptDptInMC.cxx:1603
 AliDptDptInMC.cxx:1604
 AliDptDptInMC.cxx:1605
 AliDptDptInMC.cxx:1606
 AliDptDptInMC.cxx:1607
 AliDptDptInMC.cxx:1608
 AliDptDptInMC.cxx:1609
 AliDptDptInMC.cxx:1610
 AliDptDptInMC.cxx:1611
 AliDptDptInMC.cxx:1612
 AliDptDptInMC.cxx:1613
 AliDptDptInMC.cxx:1614
 AliDptDptInMC.cxx:1615
 AliDptDptInMC.cxx:1616
 AliDptDptInMC.cxx:1617
 AliDptDptInMC.cxx:1618
 AliDptDptInMC.cxx:1619
 AliDptDptInMC.cxx:1620
 AliDptDptInMC.cxx:1621
 AliDptDptInMC.cxx:1622
 AliDptDptInMC.cxx:1623
 AliDptDptInMC.cxx:1624
 AliDptDptInMC.cxx:1625
 AliDptDptInMC.cxx:1626
 AliDptDptInMC.cxx:1627
 AliDptDptInMC.cxx:1628
 AliDptDptInMC.cxx:1629
 AliDptDptInMC.cxx:1630
 AliDptDptInMC.cxx:1631
 AliDptDptInMC.cxx:1632
 AliDptDptInMC.cxx:1633
 AliDptDptInMC.cxx:1634
 AliDptDptInMC.cxx:1635
 AliDptDptInMC.cxx:1636
 AliDptDptInMC.cxx:1637
 AliDptDptInMC.cxx:1638
 AliDptDptInMC.cxx:1639
 AliDptDptInMC.cxx:1640
 AliDptDptInMC.cxx:1641
 AliDptDptInMC.cxx:1642
 AliDptDptInMC.cxx:1643
 AliDptDptInMC.cxx:1644
 AliDptDptInMC.cxx:1645
 AliDptDptInMC.cxx:1646
 AliDptDptInMC.cxx:1647
 AliDptDptInMC.cxx:1648
 AliDptDptInMC.cxx:1649
 AliDptDptInMC.cxx:1650
 AliDptDptInMC.cxx:1651
 AliDptDptInMC.cxx:1652
 AliDptDptInMC.cxx:1653
 AliDptDptInMC.cxx:1654
 AliDptDptInMC.cxx:1655
 AliDptDptInMC.cxx:1656
 AliDptDptInMC.cxx:1657
 AliDptDptInMC.cxx:1658
 AliDptDptInMC.cxx:1659
 AliDptDptInMC.cxx:1660
 AliDptDptInMC.cxx:1661
 AliDptDptInMC.cxx:1662
 AliDptDptInMC.cxx:1663
 AliDptDptInMC.cxx:1664
 AliDptDptInMC.cxx:1665
 AliDptDptInMC.cxx:1666
 AliDptDptInMC.cxx:1667
 AliDptDptInMC.cxx:1668
 AliDptDptInMC.cxx:1669
 AliDptDptInMC.cxx:1670
 AliDptDptInMC.cxx:1671
 AliDptDptInMC.cxx:1672
 AliDptDptInMC.cxx:1673
 AliDptDptInMC.cxx:1674
 AliDptDptInMC.cxx:1675
 AliDptDptInMC.cxx:1676
 AliDptDptInMC.cxx:1677
 AliDptDptInMC.cxx:1678
 AliDptDptInMC.cxx:1679
 AliDptDptInMC.cxx:1680
 AliDptDptInMC.cxx:1681
 AliDptDptInMC.cxx:1682
 AliDptDptInMC.cxx:1683
 AliDptDptInMC.cxx:1684
 AliDptDptInMC.cxx:1685
 AliDptDptInMC.cxx:1686
 AliDptDptInMC.cxx:1687
 AliDptDptInMC.cxx:1688
 AliDptDptInMC.cxx:1689
 AliDptDptInMC.cxx:1690
 AliDptDptInMC.cxx:1691
 AliDptDptInMC.cxx:1692
 AliDptDptInMC.cxx:1693
 AliDptDptInMC.cxx:1694
 AliDptDptInMC.cxx:1695
 AliDptDptInMC.cxx:1696
 AliDptDptInMC.cxx:1697
 AliDptDptInMC.cxx:1698
 AliDptDptInMC.cxx:1699
 AliDptDptInMC.cxx:1700
 AliDptDptInMC.cxx:1701
 AliDptDptInMC.cxx:1702
 AliDptDptInMC.cxx:1703
 AliDptDptInMC.cxx:1704
 AliDptDptInMC.cxx:1705
 AliDptDptInMC.cxx:1706
 AliDptDptInMC.cxx:1707
 AliDptDptInMC.cxx:1708
 AliDptDptInMC.cxx:1709
 AliDptDptInMC.cxx:1710
 AliDptDptInMC.cxx:1711
 AliDptDptInMC.cxx:1712
 AliDptDptInMC.cxx:1713
 AliDptDptInMC.cxx:1714
 AliDptDptInMC.cxx:1715
 AliDptDptInMC.cxx:1716
 AliDptDptInMC.cxx:1717
 AliDptDptInMC.cxx:1718
 AliDptDptInMC.cxx:1719
 AliDptDptInMC.cxx:1720
 AliDptDptInMC.cxx:1721
 AliDptDptInMC.cxx:1722
 AliDptDptInMC.cxx:1723
 AliDptDptInMC.cxx:1724
 AliDptDptInMC.cxx:1725
 AliDptDptInMC.cxx:1726
 AliDptDptInMC.cxx:1727
 AliDptDptInMC.cxx:1728
 AliDptDptInMC.cxx:1729
 AliDptDptInMC.cxx:1730
 AliDptDptInMC.cxx:1731
 AliDptDptInMC.cxx:1732
 AliDptDptInMC.cxx:1733
 AliDptDptInMC.cxx:1734
 AliDptDptInMC.cxx:1735
 AliDptDptInMC.cxx:1736
 AliDptDptInMC.cxx:1737
 AliDptDptInMC.cxx:1738
 AliDptDptInMC.cxx:1739
 AliDptDptInMC.cxx:1740
 AliDptDptInMC.cxx:1741
 AliDptDptInMC.cxx:1742
 AliDptDptInMC.cxx:1743
 AliDptDptInMC.cxx:1744
 AliDptDptInMC.cxx:1745
 AliDptDptInMC.cxx:1746
 AliDptDptInMC.cxx:1747
 AliDptDptInMC.cxx:1748
 AliDptDptInMC.cxx:1749
 AliDptDptInMC.cxx:1750
 AliDptDptInMC.cxx:1751
 AliDptDptInMC.cxx:1752
 AliDptDptInMC.cxx:1753
 AliDptDptInMC.cxx:1754
 AliDptDptInMC.cxx:1755
 AliDptDptInMC.cxx:1756
 AliDptDptInMC.cxx:1757
 AliDptDptInMC.cxx:1758
 AliDptDptInMC.cxx:1759
 AliDptDptInMC.cxx:1760
 AliDptDptInMC.cxx:1761
 AliDptDptInMC.cxx:1762
 AliDptDptInMC.cxx:1763
 AliDptDptInMC.cxx:1764
 AliDptDptInMC.cxx:1765
 AliDptDptInMC.cxx:1766
 AliDptDptInMC.cxx:1767
 AliDptDptInMC.cxx:1768
 AliDptDptInMC.cxx:1769
 AliDptDptInMC.cxx:1770
 AliDptDptInMC.cxx:1771
 AliDptDptInMC.cxx:1772
 AliDptDptInMC.cxx:1773
 AliDptDptInMC.cxx:1774
 AliDptDptInMC.cxx:1775
 AliDptDptInMC.cxx:1776
 AliDptDptInMC.cxx:1777
 AliDptDptInMC.cxx:1778
 AliDptDptInMC.cxx:1779
 AliDptDptInMC.cxx:1780
 AliDptDptInMC.cxx:1781
 AliDptDptInMC.cxx:1782
 AliDptDptInMC.cxx:1783
 AliDptDptInMC.cxx:1784
 AliDptDptInMC.cxx:1785
 AliDptDptInMC.cxx:1786
 AliDptDptInMC.cxx:1787
 AliDptDptInMC.cxx:1788
 AliDptDptInMC.cxx:1789
 AliDptDptInMC.cxx:1790
 AliDptDptInMC.cxx:1791
 AliDptDptInMC.cxx:1792
 AliDptDptInMC.cxx:1793
 AliDptDptInMC.cxx:1794
 AliDptDptInMC.cxx:1795
 AliDptDptInMC.cxx:1796
 AliDptDptInMC.cxx:1797
 AliDptDptInMC.cxx:1798
 AliDptDptInMC.cxx:1799
 AliDptDptInMC.cxx:1800
 AliDptDptInMC.cxx:1801
 AliDptDptInMC.cxx:1802
 AliDptDptInMC.cxx:1803
 AliDptDptInMC.cxx:1804
 AliDptDptInMC.cxx:1805
 AliDptDptInMC.cxx:1806
 AliDptDptInMC.cxx:1807
 AliDptDptInMC.cxx:1808
 AliDptDptInMC.cxx:1809
 AliDptDptInMC.cxx:1810
 AliDptDptInMC.cxx:1811
 AliDptDptInMC.cxx:1812
 AliDptDptInMC.cxx:1813
 AliDptDptInMC.cxx:1814
 AliDptDptInMC.cxx:1815
 AliDptDptInMC.cxx:1816
 AliDptDptInMC.cxx:1817
 AliDptDptInMC.cxx:1818
 AliDptDptInMC.cxx:1819
 AliDptDptInMC.cxx:1820
 AliDptDptInMC.cxx:1821
 AliDptDptInMC.cxx:1822
 AliDptDptInMC.cxx:1823
 AliDptDptInMC.cxx:1824
 AliDptDptInMC.cxx:1825
 AliDptDptInMC.cxx:1826
 AliDptDptInMC.cxx:1827
 AliDptDptInMC.cxx:1828
 AliDptDptInMC.cxx:1829
 AliDptDptInMC.cxx:1830
 AliDptDptInMC.cxx:1831
 AliDptDptInMC.cxx:1832
 AliDptDptInMC.cxx:1833
 AliDptDptInMC.cxx:1834
 AliDptDptInMC.cxx:1835
 AliDptDptInMC.cxx:1836
 AliDptDptInMC.cxx:1837
 AliDptDptInMC.cxx:1838
 AliDptDptInMC.cxx:1839
 AliDptDptInMC.cxx:1840
 AliDptDptInMC.cxx:1841
 AliDptDptInMC.cxx:1842
 AliDptDptInMC.cxx:1843
 AliDptDptInMC.cxx:1844
 AliDptDptInMC.cxx:1845
 AliDptDptInMC.cxx:1846
 AliDptDptInMC.cxx:1847
 AliDptDptInMC.cxx:1848
 AliDptDptInMC.cxx:1849
 AliDptDptInMC.cxx:1850
 AliDptDptInMC.cxx:1851
 AliDptDptInMC.cxx:1852
 AliDptDptInMC.cxx:1853
 AliDptDptInMC.cxx:1854
 AliDptDptInMC.cxx:1855
 AliDptDptInMC.cxx:1856
 AliDptDptInMC.cxx:1857
 AliDptDptInMC.cxx:1858
 AliDptDptInMC.cxx:1859
 AliDptDptInMC.cxx:1860
 AliDptDptInMC.cxx:1861
 AliDptDptInMC.cxx:1862
 AliDptDptInMC.cxx:1863
 AliDptDptInMC.cxx:1864
 AliDptDptInMC.cxx:1865
 AliDptDptInMC.cxx:1866
 AliDptDptInMC.cxx:1867
 AliDptDptInMC.cxx:1868
 AliDptDptInMC.cxx:1869
 AliDptDptInMC.cxx:1870
 AliDptDptInMC.cxx:1871
 AliDptDptInMC.cxx:1872
 AliDptDptInMC.cxx:1873
 AliDptDptInMC.cxx:1874
 AliDptDptInMC.cxx:1875
 AliDptDptInMC.cxx:1876
 AliDptDptInMC.cxx:1877
 AliDptDptInMC.cxx:1878
 AliDptDptInMC.cxx:1879
 AliDptDptInMC.cxx:1880
 AliDptDptInMC.cxx:1881
 AliDptDptInMC.cxx:1882
 AliDptDptInMC.cxx:1883
 AliDptDptInMC.cxx:1884
 AliDptDptInMC.cxx:1885
 AliDptDptInMC.cxx:1886
 AliDptDptInMC.cxx:1887
 AliDptDptInMC.cxx:1888
 AliDptDptInMC.cxx:1889
 AliDptDptInMC.cxx:1890
 AliDptDptInMC.cxx:1891
 AliDptDptInMC.cxx:1892
 AliDptDptInMC.cxx:1893
 AliDptDptInMC.cxx:1894
 AliDptDptInMC.cxx:1895
 AliDptDptInMC.cxx:1896
 AliDptDptInMC.cxx:1897
 AliDptDptInMC.cxx:1898
 AliDptDptInMC.cxx:1899
 AliDptDptInMC.cxx:1900
 AliDptDptInMC.cxx:1901
 AliDptDptInMC.cxx:1902
 AliDptDptInMC.cxx:1903
 AliDptDptInMC.cxx:1904
 AliDptDptInMC.cxx:1905
 AliDptDptInMC.cxx:1906
 AliDptDptInMC.cxx:1907
 AliDptDptInMC.cxx:1908
 AliDptDptInMC.cxx:1909
 AliDptDptInMC.cxx:1910
 AliDptDptInMC.cxx:1911
 AliDptDptInMC.cxx:1912
 AliDptDptInMC.cxx:1913
 AliDptDptInMC.cxx:1914
 AliDptDptInMC.cxx:1915
 AliDptDptInMC.cxx:1916
 AliDptDptInMC.cxx:1917
 AliDptDptInMC.cxx:1918
 AliDptDptInMC.cxx:1919
 AliDptDptInMC.cxx:1920
 AliDptDptInMC.cxx:1921
 AliDptDptInMC.cxx:1922
 AliDptDptInMC.cxx:1923
 AliDptDptInMC.cxx:1924
 AliDptDptInMC.cxx:1925
 AliDptDptInMC.cxx:1926
 AliDptDptInMC.cxx:1927
 AliDptDptInMC.cxx:1928
 AliDptDptInMC.cxx:1929
 AliDptDptInMC.cxx:1930
 AliDptDptInMC.cxx:1931
 AliDptDptInMC.cxx:1932
 AliDptDptInMC.cxx:1933
 AliDptDptInMC.cxx:1934
 AliDptDptInMC.cxx:1935
 AliDptDptInMC.cxx:1936
 AliDptDptInMC.cxx:1937
 AliDptDptInMC.cxx:1938
 AliDptDptInMC.cxx:1939
 AliDptDptInMC.cxx:1940
 AliDptDptInMC.cxx:1941
 AliDptDptInMC.cxx:1942
 AliDptDptInMC.cxx:1943
 AliDptDptInMC.cxx:1944
 AliDptDptInMC.cxx:1945
 AliDptDptInMC.cxx:1946
 AliDptDptInMC.cxx:1947
 AliDptDptInMC.cxx:1948
 AliDptDptInMC.cxx:1949
 AliDptDptInMC.cxx:1950
 AliDptDptInMC.cxx:1951
 AliDptDptInMC.cxx:1952
 AliDptDptInMC.cxx:1953
 AliDptDptInMC.cxx:1954
 AliDptDptInMC.cxx:1955
 AliDptDptInMC.cxx:1956
 AliDptDptInMC.cxx:1957
 AliDptDptInMC.cxx:1958
 AliDptDptInMC.cxx:1959
 AliDptDptInMC.cxx:1960
 AliDptDptInMC.cxx:1961
 AliDptDptInMC.cxx:1962
 AliDptDptInMC.cxx:1963
 AliDptDptInMC.cxx:1964
 AliDptDptInMC.cxx:1965
 AliDptDptInMC.cxx:1966
 AliDptDptInMC.cxx:1967
 AliDptDptInMC.cxx:1968
 AliDptDptInMC.cxx:1969
 AliDptDptInMC.cxx:1970
 AliDptDptInMC.cxx:1971
 AliDptDptInMC.cxx:1972
 AliDptDptInMC.cxx:1973
 AliDptDptInMC.cxx:1974
 AliDptDptInMC.cxx:1975
 AliDptDptInMC.cxx:1976
 AliDptDptInMC.cxx:1977
 AliDptDptInMC.cxx:1978
 AliDptDptInMC.cxx:1979
 AliDptDptInMC.cxx:1980
 AliDptDptInMC.cxx:1981
 AliDptDptInMC.cxx:1982
 AliDptDptInMC.cxx:1983
 AliDptDptInMC.cxx:1984
 AliDptDptInMC.cxx:1985
 AliDptDptInMC.cxx:1986
 AliDptDptInMC.cxx:1987
 AliDptDptInMC.cxx:1988
 AliDptDptInMC.cxx:1989
 AliDptDptInMC.cxx:1990
 AliDptDptInMC.cxx:1991
 AliDptDptInMC.cxx:1992
 AliDptDptInMC.cxx:1993
 AliDptDptInMC.cxx:1994
 AliDptDptInMC.cxx:1995
 AliDptDptInMC.cxx:1996
 AliDptDptInMC.cxx:1997
 AliDptDptInMC.cxx:1998
 AliDptDptInMC.cxx:1999
 AliDptDptInMC.cxx:2000
 AliDptDptInMC.cxx:2001
 AliDptDptInMC.cxx:2002
 AliDptDptInMC.cxx:2003
 AliDptDptInMC.cxx:2004
 AliDptDptInMC.cxx:2005
 AliDptDptInMC.cxx:2006
 AliDptDptInMC.cxx:2007
 AliDptDptInMC.cxx:2008
 AliDptDptInMC.cxx:2009
 AliDptDptInMC.cxx:2010
 AliDptDptInMC.cxx:2011
 AliDptDptInMC.cxx:2012
 AliDptDptInMC.cxx:2013
 AliDptDptInMC.cxx:2014
 AliDptDptInMC.cxx:2015
 AliDptDptInMC.cxx:2016
 AliDptDptInMC.cxx:2017
 AliDptDptInMC.cxx:2018
 AliDptDptInMC.cxx:2019
 AliDptDptInMC.cxx:2020
 AliDptDptInMC.cxx:2021
 AliDptDptInMC.cxx:2022
 AliDptDptInMC.cxx:2023
 AliDptDptInMC.cxx:2024
 AliDptDptInMC.cxx:2025
 AliDptDptInMC.cxx:2026
 AliDptDptInMC.cxx:2027
 AliDptDptInMC.cxx:2028
 AliDptDptInMC.cxx:2029
 AliDptDptInMC.cxx:2030
 AliDptDptInMC.cxx:2031
 AliDptDptInMC.cxx:2032
 AliDptDptInMC.cxx:2033
 AliDptDptInMC.cxx:2034
 AliDptDptInMC.cxx:2035
 AliDptDptInMC.cxx:2036
 AliDptDptInMC.cxx:2037
 AliDptDptInMC.cxx:2038
 AliDptDptInMC.cxx:2039
 AliDptDptInMC.cxx:2040
 AliDptDptInMC.cxx:2041
 AliDptDptInMC.cxx:2042
 AliDptDptInMC.cxx:2043
 AliDptDptInMC.cxx:2044
 AliDptDptInMC.cxx:2045
 AliDptDptInMC.cxx:2046
 AliDptDptInMC.cxx:2047
 AliDptDptInMC.cxx:2048
 AliDptDptInMC.cxx:2049
 AliDptDptInMC.cxx:2050
 AliDptDptInMC.cxx:2051
 AliDptDptInMC.cxx:2052
 AliDptDptInMC.cxx:2053
 AliDptDptInMC.cxx:2054
 AliDptDptInMC.cxx:2055
 AliDptDptInMC.cxx:2056
 AliDptDptInMC.cxx:2057
 AliDptDptInMC.cxx:2058
 AliDptDptInMC.cxx:2059
 AliDptDptInMC.cxx:2060
 AliDptDptInMC.cxx:2061
 AliDptDptInMC.cxx:2062
 AliDptDptInMC.cxx:2063
 AliDptDptInMC.cxx:2064
 AliDptDptInMC.cxx:2065
 AliDptDptInMC.cxx:2066
 AliDptDptInMC.cxx:2067
 AliDptDptInMC.cxx:2068
 AliDptDptInMC.cxx:2069
 AliDptDptInMC.cxx:2070
 AliDptDptInMC.cxx:2071
 AliDptDptInMC.cxx:2072
 AliDptDptInMC.cxx:2073
 AliDptDptInMC.cxx:2074
 AliDptDptInMC.cxx:2075
 AliDptDptInMC.cxx:2076
 AliDptDptInMC.cxx:2077
 AliDptDptInMC.cxx:2078
 AliDptDptInMC.cxx:2079
 AliDptDptInMC.cxx:2080
 AliDptDptInMC.cxx:2081
 AliDptDptInMC.cxx:2082
 AliDptDptInMC.cxx:2083
 AliDptDptInMC.cxx:2084
 AliDptDptInMC.cxx:2085
 AliDptDptInMC.cxx:2086
 AliDptDptInMC.cxx:2087
 AliDptDptInMC.cxx:2088
 AliDptDptInMC.cxx:2089
 AliDptDptInMC.cxx:2090
 AliDptDptInMC.cxx:2091
 AliDptDptInMC.cxx:2092
 AliDptDptInMC.cxx:2093
 AliDptDptInMC.cxx:2094
 AliDptDptInMC.cxx:2095
 AliDptDptInMC.cxx:2096
 AliDptDptInMC.cxx:2097
 AliDptDptInMC.cxx:2098
 AliDptDptInMC.cxx:2099
 AliDptDptInMC.cxx:2100
 AliDptDptInMC.cxx:2101
 AliDptDptInMC.cxx:2102
 AliDptDptInMC.cxx:2103
 AliDptDptInMC.cxx:2104
 AliDptDptInMC.cxx:2105
 AliDptDptInMC.cxx:2106
 AliDptDptInMC.cxx:2107
 AliDptDptInMC.cxx:2108
 AliDptDptInMC.cxx:2109
 AliDptDptInMC.cxx:2110
 AliDptDptInMC.cxx:2111
 AliDptDptInMC.cxx:2112
 AliDptDptInMC.cxx:2113
 AliDptDptInMC.cxx:2114
 AliDptDptInMC.cxx:2115
 AliDptDptInMC.cxx:2116
 AliDptDptInMC.cxx:2117
 AliDptDptInMC.cxx:2118
 AliDptDptInMC.cxx:2119
 AliDptDptInMC.cxx:2120
 AliDptDptInMC.cxx:2121
 AliDptDptInMC.cxx:2122
 AliDptDptInMC.cxx:2123
 AliDptDptInMC.cxx:2124
 AliDptDptInMC.cxx:2125
 AliDptDptInMC.cxx:2126
 AliDptDptInMC.cxx:2127
 AliDptDptInMC.cxx:2128
 AliDptDptInMC.cxx:2129
 AliDptDptInMC.cxx:2130
 AliDptDptInMC.cxx:2131
 AliDptDptInMC.cxx:2132
 AliDptDptInMC.cxx:2133
 AliDptDptInMC.cxx:2134
 AliDptDptInMC.cxx:2135
 AliDptDptInMC.cxx:2136
 AliDptDptInMC.cxx:2137
 AliDptDptInMC.cxx:2138
 AliDptDptInMC.cxx:2139
 AliDptDptInMC.cxx:2140
 AliDptDptInMC.cxx:2141
 AliDptDptInMC.cxx:2142
 AliDptDptInMC.cxx:2143
 AliDptDptInMC.cxx:2144
 AliDptDptInMC.cxx:2145
 AliDptDptInMC.cxx:2146
 AliDptDptInMC.cxx:2147
 AliDptDptInMC.cxx:2148
 AliDptDptInMC.cxx:2149
 AliDptDptInMC.cxx:2150
 AliDptDptInMC.cxx:2151
 AliDptDptInMC.cxx:2152
 AliDptDptInMC.cxx:2153
 AliDptDptInMC.cxx:2154
 AliDptDptInMC.cxx:2155
 AliDptDptInMC.cxx:2156
 AliDptDptInMC.cxx:2157
 AliDptDptInMC.cxx:2158
 AliDptDptInMC.cxx:2159
 AliDptDptInMC.cxx:2160
 AliDptDptInMC.cxx:2161
 AliDptDptInMC.cxx:2162
 AliDptDptInMC.cxx:2163
 AliDptDptInMC.cxx:2164
 AliDptDptInMC.cxx:2165
 AliDptDptInMC.cxx:2166
 AliDptDptInMC.cxx:2167
 AliDptDptInMC.cxx:2168
 AliDptDptInMC.cxx:2169
 AliDptDptInMC.cxx:2170
 AliDptDptInMC.cxx:2171
 AliDptDptInMC.cxx:2172
 AliDptDptInMC.cxx:2173
 AliDptDptInMC.cxx:2174
 AliDptDptInMC.cxx:2175
 AliDptDptInMC.cxx:2176
 AliDptDptInMC.cxx:2177
 AliDptDptInMC.cxx:2178
 AliDptDptInMC.cxx:2179
 AliDptDptInMC.cxx:2180
 AliDptDptInMC.cxx:2181
 AliDptDptInMC.cxx:2182
 AliDptDptInMC.cxx:2183
 AliDptDptInMC.cxx:2184
 AliDptDptInMC.cxx:2185
 AliDptDptInMC.cxx:2186
 AliDptDptInMC.cxx:2187
 AliDptDptInMC.cxx:2188
 AliDptDptInMC.cxx:2189
 AliDptDptInMC.cxx:2190
 AliDptDptInMC.cxx:2191
 AliDptDptInMC.cxx:2192
 AliDptDptInMC.cxx:2193
 AliDptDptInMC.cxx:2194
 AliDptDptInMC.cxx:2195
 AliDptDptInMC.cxx:2196
 AliDptDptInMC.cxx:2197
 AliDptDptInMC.cxx:2198
 AliDptDptInMC.cxx:2199
 AliDptDptInMC.cxx:2200
 AliDptDptInMC.cxx:2201
 AliDptDptInMC.cxx:2202
 AliDptDptInMC.cxx:2203
 AliDptDptInMC.cxx:2204
 AliDptDptInMC.cxx:2205
 AliDptDptInMC.cxx:2206
 AliDptDptInMC.cxx:2207
 AliDptDptInMC.cxx:2208
 AliDptDptInMC.cxx:2209
 AliDptDptInMC.cxx:2210
 AliDptDptInMC.cxx:2211
 AliDptDptInMC.cxx:2212
 AliDptDptInMC.cxx:2213
 AliDptDptInMC.cxx:2214
 AliDptDptInMC.cxx:2215
 AliDptDptInMC.cxx:2216
 AliDptDptInMC.cxx:2217
 AliDptDptInMC.cxx:2218
 AliDptDptInMC.cxx:2219
 AliDptDptInMC.cxx:2220
 AliDptDptInMC.cxx:2221
 AliDptDptInMC.cxx:2222
 AliDptDptInMC.cxx:2223
 AliDptDptInMC.cxx:2224
 AliDptDptInMC.cxx:2225
 AliDptDptInMC.cxx:2226
 AliDptDptInMC.cxx:2227
 AliDptDptInMC.cxx:2228
 AliDptDptInMC.cxx:2229
 AliDptDptInMC.cxx:2230
 AliDptDptInMC.cxx:2231
 AliDptDptInMC.cxx:2232
 AliDptDptInMC.cxx:2233
 AliDptDptInMC.cxx:2234
 AliDptDptInMC.cxx:2235
 AliDptDptInMC.cxx:2236
 AliDptDptInMC.cxx:2237
 AliDptDptInMC.cxx:2238
 AliDptDptInMC.cxx:2239
 AliDptDptInMC.cxx:2240
 AliDptDptInMC.cxx:2241
 AliDptDptInMC.cxx:2242
 AliDptDptInMC.cxx:2243
 AliDptDptInMC.cxx:2244
 AliDptDptInMC.cxx:2245
 AliDptDptInMC.cxx:2246
 AliDptDptInMC.cxx:2247
 AliDptDptInMC.cxx:2248
 AliDptDptInMC.cxx:2249
 AliDptDptInMC.cxx:2250
 AliDptDptInMC.cxx:2251
 AliDptDptInMC.cxx:2252
 AliDptDptInMC.cxx:2253
 AliDptDptInMC.cxx:2254
 AliDptDptInMC.cxx:2255
 AliDptDptInMC.cxx:2256
 AliDptDptInMC.cxx:2257
 AliDptDptInMC.cxx:2258
 AliDptDptInMC.cxx:2259
 AliDptDptInMC.cxx:2260
 AliDptDptInMC.cxx:2261
 AliDptDptInMC.cxx:2262
 AliDptDptInMC.cxx:2263
 AliDptDptInMC.cxx:2264
 AliDptDptInMC.cxx:2265
 AliDptDptInMC.cxx:2266
 AliDptDptInMC.cxx:2267
 AliDptDptInMC.cxx:2268
 AliDptDptInMC.cxx:2269
 AliDptDptInMC.cxx:2270
 AliDptDptInMC.cxx:2271
 AliDptDptInMC.cxx:2272
 AliDptDptInMC.cxx:2273
 AliDptDptInMC.cxx:2274
 AliDptDptInMC.cxx:2275
 AliDptDptInMC.cxx:2276
 AliDptDptInMC.cxx:2277
 AliDptDptInMC.cxx:2278
 AliDptDptInMC.cxx:2279
 AliDptDptInMC.cxx:2280
 AliDptDptInMC.cxx:2281
 AliDptDptInMC.cxx:2282
 AliDptDptInMC.cxx:2283
 AliDptDptInMC.cxx:2284
 AliDptDptInMC.cxx:2285
 AliDptDptInMC.cxx:2286
 AliDptDptInMC.cxx:2287
 AliDptDptInMC.cxx:2288
 AliDptDptInMC.cxx:2289
 AliDptDptInMC.cxx:2290
 AliDptDptInMC.cxx:2291
 AliDptDptInMC.cxx:2292
 AliDptDptInMC.cxx:2293
 AliDptDptInMC.cxx:2294
 AliDptDptInMC.cxx:2295
 AliDptDptInMC.cxx:2296
 AliDptDptInMC.cxx:2297
 AliDptDptInMC.cxx:2298
 AliDptDptInMC.cxx:2299
 AliDptDptInMC.cxx:2300
 AliDptDptInMC.cxx:2301
 AliDptDptInMC.cxx:2302
 AliDptDptInMC.cxx:2303
 AliDptDptInMC.cxx:2304
 AliDptDptInMC.cxx:2305
 AliDptDptInMC.cxx:2306
 AliDptDptInMC.cxx:2307
 AliDptDptInMC.cxx:2308
 AliDptDptInMC.cxx:2309
 AliDptDptInMC.cxx:2310
 AliDptDptInMC.cxx:2311
 AliDptDptInMC.cxx:2312
 AliDptDptInMC.cxx:2313
 AliDptDptInMC.cxx:2314
 AliDptDptInMC.cxx:2315
 AliDptDptInMC.cxx:2316
 AliDptDptInMC.cxx:2317
 AliDptDptInMC.cxx:2318
 AliDptDptInMC.cxx:2319
 AliDptDptInMC.cxx:2320
 AliDptDptInMC.cxx:2321
 AliDptDptInMC.cxx:2322
 AliDptDptInMC.cxx:2323
 AliDptDptInMC.cxx:2324
 AliDptDptInMC.cxx:2325
 AliDptDptInMC.cxx:2326
 AliDptDptInMC.cxx:2327
 AliDptDptInMC.cxx:2328
 AliDptDptInMC.cxx:2329
 AliDptDptInMC.cxx:2330
 AliDptDptInMC.cxx:2331
 AliDptDptInMC.cxx:2332
 AliDptDptInMC.cxx:2333
 AliDptDptInMC.cxx:2334
 AliDptDptInMC.cxx:2335
 AliDptDptInMC.cxx:2336
 AliDptDptInMC.cxx:2337
 AliDptDptInMC.cxx:2338
 AliDptDptInMC.cxx:2339
 AliDptDptInMC.cxx:2340
 AliDptDptInMC.cxx:2341
 AliDptDptInMC.cxx:2342
 AliDptDptInMC.cxx:2343
 AliDptDptInMC.cxx:2344
 AliDptDptInMC.cxx:2345
 AliDptDptInMC.cxx:2346
 AliDptDptInMC.cxx:2347
 AliDptDptInMC.cxx:2348
 AliDptDptInMC.cxx:2349
 AliDptDptInMC.cxx:2350
 AliDptDptInMC.cxx:2351
 AliDptDptInMC.cxx:2352
 AliDptDptInMC.cxx:2353
 AliDptDptInMC.cxx:2354
 AliDptDptInMC.cxx:2355
 AliDptDptInMC.cxx:2356
 AliDptDptInMC.cxx:2357
 AliDptDptInMC.cxx:2358
 AliDptDptInMC.cxx:2359
 AliDptDptInMC.cxx:2360
 AliDptDptInMC.cxx:2361
 AliDptDptInMC.cxx:2362
 AliDptDptInMC.cxx:2363
 AliDptDptInMC.cxx:2364
 AliDptDptInMC.cxx:2365
 AliDptDptInMC.cxx:2366
 AliDptDptInMC.cxx:2367
 AliDptDptInMC.cxx:2368
 AliDptDptInMC.cxx:2369
 AliDptDptInMC.cxx:2370
 AliDptDptInMC.cxx:2371
 AliDptDptInMC.cxx:2372
 AliDptDptInMC.cxx:2373
 AliDptDptInMC.cxx:2374
 AliDptDptInMC.cxx:2375
 AliDptDptInMC.cxx:2376
 AliDptDptInMC.cxx:2377
 AliDptDptInMC.cxx:2378
 AliDptDptInMC.cxx:2379
 AliDptDptInMC.cxx:2380
 AliDptDptInMC.cxx:2381
 AliDptDptInMC.cxx:2382
 AliDptDptInMC.cxx:2383
 AliDptDptInMC.cxx:2384
 AliDptDptInMC.cxx:2385
 AliDptDptInMC.cxx:2386