[u/mrichter/AliRoot.git] / MUON / AliMUONReconstructor.cxx

/**************************************************************************
 * 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$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// class for MUON reconstruction                                             //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////
#include <TParticle.h>
#include <TArrayF.h>

#include "AliRunLoader.h"
#include "AliHeader.h"
#include "AliGenEventHeader.h"
#include "AliESD.h"
#include "AliMUONReconstructor.h"

#include "AliMUONData.h"
#include "AliMUONEventReconstructor.h"
#include "AliMUONClusterReconstructor.h"
#include "AliMUONTriggerDecision.h"
#include "AliMUONClusterFinderVS.h"
#include "AliMUONTrack.h"
#include "AliMUONTrackParam.h"
#include "AliMUONTriggerTrack.h"
#include "AliESDMuonTrack.h"
#include "AliRawReader.h"

ClassImp(AliMUONReconstructor)
//_____________________________________________________________________________
AliMUONReconstructor::AliMUONReconstructor()
  : AliReconstructor()
{
}
//_____________________________________________________________________________
AliMUONReconstructor::~AliMUONReconstructor()
{
}
//_____________________________________________________________________________
void AliMUONReconstructor::Reconstruct(AliRunLoader* runLoader) const
{
//  AliLoader
  AliLoader* loader = runLoader->GetLoader("MUONLoader");
  Int_t nEvents = runLoader->GetNumberOfEvents();

// used local container for each method
// passing fLoader as argument, could be avoided ???
  AliMUONEventReconstructor* recoEvent = new AliMUONEventReconstructor(loader);
  AliMUONData* dataEvent = recoEvent->GetMUONData();

  AliMUONClusterReconstructor* recoCluster = new AliMUONClusterReconstructor(loader);
  AliMUONData* dataCluster = recoCluster->GetMUONData();

  AliMUONTriggerDecision* trigDec = new AliMUONTriggerDecision(loader,0,dataCluster);
  //  AliMUONData* dataTrig = trigDec->GetMUONData();


  for (Int_t i = 0; i < 10; i++) {
    AliMUONClusterFinderVS *recModel = new AliMUONClusterFinderVS();
    recModel->SetGhostChi2Cut(10);
    recoCluster->SetReconstructionModel(i,recModel);
  } 

  loader->LoadDigits("READ");
  loader->LoadRecPoints("RECREATE");
  loader->LoadTracks("RECREATE");
  
  //   Loop over events              
  for(Int_t ievent = 0; ievent < nEvents; ievent++) {
    printf("Event %d\n",ievent);
    runLoader->GetEvent(ievent);

    //----------------------- digit2cluster & Trigger2Trigger -------------------
    if (!loader->TreeR()) loader->MakeRecPointsContainer();
     
    // tracking branch
    dataCluster->MakeBranch("RC");
    dataCluster->SetTreeAddress("D,RC");
    recoCluster->Digits2Clusters(); 
    dataCluster->Fill("RC"); 

    // trigger branch
    dataCluster->MakeBranch("TC");
    dataCluster->SetTreeAddress("TC");
    trigDec->Trigger2Trigger(); 
    dataCluster->Fill("TC");

    loader->WriteRecPoints("OVERWRITE");

    //---------------------------- Track & TriggerTrack ---------------------
    if (!loader->TreeT()) loader->MakeTracksContainer();

    // trigger branch
    dataEvent->MakeBranch("RL"); //trigger track
    dataEvent->SetTreeAddress("RL");
    recoEvent->EventReconstructTrigger();
    dataEvent->Fill("RL");

    // tracking branch
    dataEvent->MakeBranch("RT"); //track
    dataEvent->SetTreeAddress("RT");
    recoEvent->EventReconstruct();
    dataEvent->Fill("RT");

    loader->WriteTracks("OVERWRITE");  
  
    //--------------------------- Resetting branches -----------------------
    dataCluster->ResetDigits();
    dataCluster->ResetRawClusters();
    dataCluster->ResetTrigger();

    dataEvent->ResetRawClusters();
    dataEvent->ResetTrigger();
    dataEvent->ResetRecTracks();
    dataEvent->ResetRecTriggerTracks();
  
  }
  loader->UnloadDigits();
  loader->UnloadRecPoints();
  loader->UnloadTracks();

  delete recoCluster;
  delete recoEvent;
  delete trigDec;
}

//_____________________________________________________________________________
void AliMUONReconstructor::Reconstruct(AliRunLoader* runLoader, AliRawReader* rawReader) const
{
//  AliLoader
  AliLoader* loader = runLoader->GetLoader("MUONLoader");

// used local container for each method
// passing fLoader as argument, could be avoided ???
  AliMUONEventReconstructor* recoEvent = new AliMUONEventReconstructor(loader);
  AliMUONData* dataEvent = recoEvent->GetMUONData();

  AliMUONClusterReconstructor* recoCluster = new AliMUONClusterReconstructor(loader);
  AliMUONData* dataCluster = recoCluster->GetMUONData();

  AliMUONTriggerDecision* trigDec = new AliMUONTriggerDecision(loader,0,dataCluster);
 
  for (Int_t i = 0; i < 10; i++) {
    AliMUONClusterFinderVS *recModel = new AliMUONClusterFinderVS();
    recModel->SetGhostChi2Cut(10);
    recoCluster->SetReconstructionModel(i,recModel);
  } 

  loader->LoadRecPoints("RECREATE");
  loader->LoadTracks("RECREATE");
  
  //   Loop over events  
  Int_t iEvent = 0;
            
  while (rawReader->NextEvent()) {
    printf("Event %d\n",iEvent);
    runLoader->GetEvent(iEvent++);

    //----------------------- digit2cluster & Trigger2Trigger -------------------
    if (!loader->TreeR()) loader->MakeRecPointsContainer();
     
    // tracking branch
    dataCluster->MakeBranch("RC");
    dataCluster->SetTreeAddress("D,RC");
    recoCluster->Digits2Clusters(rawReader); 
    dataCluster->Fill("RC"); 

    // trigger branch
    dataCluster->MakeBranch("TC");
    dataCluster->SetTreeAddress("TC");
    trigDec->Trigger2Trigger(rawReader); 
    dataCluster->Fill("TC");

    loader->WriteRecPoints("OVERWRITE");

    //---------------------------- Track & TriggerTrack ---------------------
    if (!loader->TreeT()) loader->MakeTracksContainer();

    // trigger branch
    dataEvent->MakeBranch("RL"); //trigger track
    dataEvent->SetTreeAddress("RL");
    recoEvent->EventReconstructTrigger();
    dataEvent->Fill("RL");

    // tracking branch
    dataEvent->MakeBranch("RT"); //track
    dataEvent->SetTreeAddress("RT");
    recoEvent->EventReconstruct();
    dataEvent->Fill("RT");

    loader->WriteTracks("OVERWRITE");  
  
    //--------------------------- Resetting branches -----------------------
    dataCluster->ResetDigits();
    dataCluster->ResetRawClusters();
    dataCluster->ResetTrigger();

    dataEvent->ResetRawClusters();
    dataEvent->ResetTrigger();
    dataEvent->ResetRecTracks();
    dataEvent->ResetRecTriggerTracks();
  
  }
  loader->UnloadRecPoints();
  loader->UnloadTracks();

  delete recoCluster;
  delete recoEvent;
  delete trigDec;
}

//_____________________________________________________________________________
void AliMUONReconstructor::FillESD(AliRunLoader* runLoader, AliESD* esd) const
{
  TClonesArray* recTracksArray = 0;
  TClonesArray* recTrigTracksArray = 0;
  
  AliLoader* loader = runLoader->GetLoader("MUONLoader");
  loader->LoadTracks("READ");
  AliMUONData* muonData = new AliMUONData(loader,"MUON","MUON");

   // declaration  
  Int_t iEvent, nPart;
  Int_t nTrackHits, nPrimary;
  Double_t fitFmin;
  TArrayF vertex(3);

  Double_t bendingSlope, nonBendingSlope, inverseBendingMomentum;
  Double_t xRec, yRec, zRec, chi2MatchTrigger;
  Bool_t matchTrigger;

  // setting pointer for tracks, triggertracks & trackparam at vertex
  AliMUONTrack* recTrack = 0;
  AliMUONTrackParam* trackParam = 0;
  AliMUONTriggerTrack* recTriggerTrack = 0;
  TParticle* particle = new TParticle();
  AliGenEventHeader* header = 0;
  iEvent = runLoader->GetEventNumber(); 
  runLoader->GetEvent(iEvent);

  // vertex calculation (maybe it exists already somewhere else)
  vertex[0] = vertex[1] = vertex[2] = 0.;
  nPrimary = 0;
  if ( (header = runLoader->GetHeader()->GenEventHeader()) ) {
    header->PrimaryVertex(vertex);
  } else {
    runLoader->LoadKinematics("READ");
    runLoader->TreeK()->GetBranch("Particles")->SetAddress(&particle);
    nPart = (Int_t)runLoader->TreeK()->GetEntries();
    for(Int_t iPart = 0; iPart < nPart; iPart++) {
      runLoader->TreeK()->GetEvent(iPart);
      if (particle->GetFirstMother() == -1) {
	vertex[0] += particle->Vx();
	vertex[1] += particle->Vy();
	vertex[2] += particle->Vz();
	nPrimary++;
      }
      if (nPrimary) {
	vertex[0] /= (double)nPrimary;
	vertex[1] /= (double)nPrimary;
	vertex[2] /= (double)nPrimary;
      }
    }
  }
  // setting ESD MUON class
  AliESDMuonTrack* theESDTrack = new  AliESDMuonTrack() ;

  //-------------------- trigger tracks-------------
  Long_t trigPat = 0;
  muonData->SetTreeAddress("RL");
  muonData->GetRecTriggerTracks();
  recTrigTracksArray = muonData->RecTriggerTracks();

  // ready global trigger pattern from first track
  if (recTrigTracksArray) 
    recTriggerTrack = (AliMUONTriggerTrack*) recTrigTracksArray->First();
  if (recTriggerTrack) trigPat = recTriggerTrack->GetGTPattern();

  //printf(">>> Event %d Number of Recconstructed tracks %d \n",iEvent, nrectracks);
 
  // -------------------- tracks-------------
  muonData->SetTreeAddress("RT");
  muonData->GetRecTracks();
  recTracksArray = muonData->RecTracks();
        
  Int_t nRecTracks = 0;
  if (recTracksArray)
    nRecTracks = (Int_t) recTracksArray->GetEntriesFast(); //
  
  // loop over tracks
  for (Int_t iRecTracks = 0; iRecTracks <  nRecTracks;  iRecTracks++) {

    // reading info from tracks
    recTrack = (AliMUONTrack*) recTracksArray->At(iRecTracks);

    trackParam = (AliMUONTrackParam*) (recTrack->GetTrackParamAtHit())->First();
    trackParam->ExtrapToVertex(vertex[0],vertex[1],vertex[2]);

    bendingSlope            = trackParam->GetBendingSlope();
    nonBendingSlope         = trackParam->GetNonBendingSlope();
    inverseBendingMomentum = trackParam->GetInverseBendingMomentum();
    xRec  = trackParam->GetNonBendingCoor();
    yRec  = trackParam->GetBendingCoor();
    zRec  = trackParam->GetZ();

    nTrackHits       = recTrack->GetNTrackHits();
    fitFmin          = recTrack->GetFitFMin();
    matchTrigger     = recTrack->GetMatchTrigger();
    chi2MatchTrigger = recTrack->GetChi2MatchTrigger();

    // setting data member of ESD MUON
    theESDTrack->SetInverseBendingMomentum(inverseBendingMomentum);
    theESDTrack->SetThetaX(TMath::ATan(nonBendingSlope));
    theESDTrack->SetThetaY(TMath::ATan(bendingSlope));
    theESDTrack->SetZ(zRec);
    theESDTrack->SetBendingCoor(yRec); // calculate vertex at ESD or Tracking level ?
    theESDTrack->SetNonBendingCoor(xRec);
    theESDTrack->SetChi2(fitFmin);
    theESDTrack->SetNHit(nTrackHits);
    theESDTrack->SetMatchTrigger(matchTrigger);
    theESDTrack->SetChi2MatchTrigger(chi2MatchTrigger);

    // storing ESD MUON Track into ESD Event 
    if (nRecTracks != 0)  
      esd->AddMuonTrack(theESDTrack);
  } // end loop tracks

  // add global trigger pattern
  if (nRecTracks != 0)  
    esd->SetTrigger(trigPat);

  // reset muondata
  muonData->ResetRecTracks();
  muonData->ResetRecTriggerTracks();

  //} // end loop on event  
  loader->UnloadTracks(); 
  if (!header)
    runLoader->UnloadKinematics();
  delete theESDTrack;
  delete muonData;
  // delete particle;
}
Commit	Line	Data
cac2eb58	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
cac2eb58	15	/* $Id$ */
	16
	17	///////////////////////////////////////////////////////////////////////////////
	18	// //
b002ac70	19	// class for MUON reconstruction //
cac2eb58	20	// //
cac2eb58	21	///////////////////////////////////////////////////////////////////////////////
b002ac70	22	#include <TParticle.h>
30178c30	23	#include <TArrayF.h>
30178c30	24
cac2eb58	25	#include "AliRunLoader.h"
b002ac70	26	#include "AliHeader.h"
b002ac70	27	#include "AliGenEventHeader.h"
cac2eb58	28	#include "AliESD.h"
cf464691	29	#include "AliMUONReconstructor.h"
30178c30	30
cac2eb58	31	#include "AliMUONData.h"
	32	#include "AliMUONEventReconstructor.h"
	33	#include "AliMUONClusterReconstructor.h"
	34	#include "AliMUONTriggerDecision.h"
	35	#include "AliMUONClusterFinderVS.h"
	36	#include "AliMUONTrack.h"
	37	#include "AliMUONTrackParam.h"
	38	#include "AliMUONTriggerTrack.h"
	39	#include "AliESDMuonTrack.h"
cf464691	40	#include "AliRawReader.h"
cac2eb58	41
	42	ClassImp(AliMUONReconstructor)
	43	//_____________________________________________________________________________
	44	AliMUONReconstructor::AliMUONReconstructor()
30178c30	45	: AliReconstructor()
8789635b	46	{
	47	}
	48	//_____________________________________________________________________________
	49	AliMUONReconstructor::~AliMUONReconstructor()
cac2eb58	50	{
	51	}
	52	//_____________________________________________________________________________
	53	void AliMUONReconstructor::Reconstruct(AliRunLoader* runLoader) const
	54	{
	55	// AliLoader
	56	AliLoader* loader = runLoader->GetLoader("MUONLoader");
	57	Int_t nEvents = runLoader->GetNumberOfEvents();
	58
	59	// used local container for each method
	60	// passing fLoader as argument, could be avoided ???
	61	AliMUONEventReconstructor* recoEvent = new AliMUONEventReconstructor(loader);
	62	AliMUONData* dataEvent = recoEvent->GetMUONData();
	63
	64	AliMUONClusterReconstructor* recoCluster = new AliMUONClusterReconstructor(loader);
	65	AliMUONData* dataCluster = recoCluster->GetMUONData();
	66
ce3e25a8	67	AliMUONTriggerDecision* trigDec = new AliMUONTriggerDecision(loader,0,dataCluster);
ce3e25a8	68	// AliMUONData* dataTrig = trigDec->GetMUONData();
cac2eb58	69
	70
	71	for (Int_t i = 0; i < 10; i++) {
30178c30	72	AliMUONClusterFinderVS *recModel = new AliMUONClusterFinderVS();
	73	recModel->SetGhostChi2Cut(10);
	74	recoCluster->SetReconstructionModel(i,recModel);
cac2eb58	75	}
	76
	77	loader->LoadDigits("READ");
	78	loader->LoadRecPoints("RECREATE");
	79	loader->LoadTracks("RECREATE");
	80
	81	// Loop over events
	82	for(Int_t ievent = 0; ievent < nEvents; ievent++) {
	83	printf("Event %d\n",ievent);
	84	runLoader->GetEvent(ievent);
	85
8210cb40	86	//----------------------- digit2cluster & Trigger2Trigger -------------------
cac2eb58	87	if (!loader->TreeR()) loader->MakeRecPointsContainer();
	88
	89	// tracking branch
	90	dataCluster->MakeBranch("RC");
	91	dataCluster->SetTreeAddress("D,RC");
	92	recoCluster->Digits2Clusters();
	93	dataCluster->Fill("RC");
	94
	95	// trigger branch
ce3e25a8	96	dataCluster->MakeBranch("TC");
	97	dataCluster->SetTreeAddress("TC");
	98	trigDec->Trigger2Trigger();
	99	dataCluster->Fill("TC");
cac2eb58	100
	101	loader->WriteRecPoints("OVERWRITE");
	102
	103	//---------------------------- Track & TriggerTrack ---------------------
	104	if (!loader->TreeT()) loader->MakeTracksContainer();
	105
	106	// trigger branch
	107	dataEvent->MakeBranch("RL"); //trigger track
	108	dataEvent->SetTreeAddress("RL");
	109	recoEvent->EventReconstructTrigger();
	110	dataEvent->Fill("RL");
	111
	112	// tracking branch
	113	dataEvent->MakeBranch("RT"); //track
	114	dataEvent->SetTreeAddress("RT");
	115	recoEvent->EventReconstruct();
	116	dataEvent->Fill("RT");
	117
	118	loader->WriteTracks("OVERWRITE");
	119
	120	//--------------------------- Resetting branches -----------------------
	121	dataCluster->ResetDigits();
	122	dataCluster->ResetRawClusters();
ce3e25a8	123	dataCluster->ResetTrigger();
cac2eb58	124
	125	dataEvent->ResetRawClusters();
	126	dataEvent->ResetTrigger();
	127	dataEvent->ResetRecTracks();
	128	dataEvent->ResetRecTriggerTracks();
	129
	130	}
	131	loader->UnloadDigits();
	132	loader->UnloadRecPoints();
	133	loader->UnloadTracks();
	134
	135	delete recoCluster;
	136	delete recoEvent;
	137	delete trigDec;
	138	}
cf464691	139
	140	//_____________________________________________________________________________
	141	void AliMUONReconstructor::Reconstruct(AliRunLoader* runLoader, AliRawReader* rawReader) const
	142	{
	143	// AliLoader
	144	AliLoader* loader = runLoader->GetLoader("MUONLoader");
	145
	146	// used local container for each method
	147	// passing fLoader as argument, could be avoided ???
	148	AliMUONEventReconstructor* recoEvent = new AliMUONEventReconstructor(loader);
	149	AliMUONData* dataEvent = recoEvent->GetMUONData();
	150
	151	AliMUONClusterReconstructor* recoCluster = new AliMUONClusterReconstructor(loader);
	152	AliMUONData* dataCluster = recoCluster->GetMUONData();
	153
	154	AliMUONTriggerDecision* trigDec = new AliMUONTriggerDecision(loader,0,dataCluster);
8210cb40	155
cf464691	156	for (Int_t i = 0; i < 10; i++) {
	157	AliMUONClusterFinderVS *recModel = new AliMUONClusterFinderVS();
	158	recModel->SetGhostChi2Cut(10);
	159	recoCluster->SetReconstructionModel(i,recModel);
	160	}
	161
cf464691	162	loader->LoadRecPoints("RECREATE");
	163	loader->LoadTracks("RECREATE");
	164
	165	// Loop over events
8210cb40	166	Int_t iEvent = 0;
cf464691	167
8210cb40	168	while (rawReader->NextEvent()) {
cf464691	169	printf("Event %d\n",iEvent);
	170	runLoader->GetEvent(iEvent++);
	171
8210cb40	172	//----------------------- digit2cluster & Trigger2Trigger -------------------
cf464691	173	if (!loader->TreeR()) loader->MakeRecPointsContainer();
	174
	175	// tracking branch
	176	dataCluster->MakeBranch("RC");
	177	dataCluster->SetTreeAddress("D,RC");
	178	recoCluster->Digits2Clusters(rawReader);
	179	dataCluster->Fill("RC");
	180
	181	// trigger branch
	182	dataCluster->MakeBranch("TC");
	183	dataCluster->SetTreeAddress("TC");
8210cb40	184	trigDec->Trigger2Trigger(rawReader);
cf464691	185	dataCluster->Fill("TC");
	186
	187	loader->WriteRecPoints("OVERWRITE");
	188
	189	//---------------------------- Track & TriggerTrack ---------------------
	190	if (!loader->TreeT()) loader->MakeTracksContainer();
	191
	192	// trigger branch
	193	dataEvent->MakeBranch("RL"); //trigger track
	194	dataEvent->SetTreeAddress("RL");
	195	recoEvent->EventReconstructTrigger();
	196	dataEvent->Fill("RL");
	197
	198	// tracking branch
	199	dataEvent->MakeBranch("RT"); //track
	200	dataEvent->SetTreeAddress("RT");
	201	recoEvent->EventReconstruct();
	202	dataEvent->Fill("RT");
	203
	204	loader->WriteTracks("OVERWRITE");
	205
	206	//--------------------------- Resetting branches -----------------------
	207	dataCluster->ResetDigits();
	208	dataCluster->ResetRawClusters();
	209	dataCluster->ResetTrigger();
	210
	211	dataEvent->ResetRawClusters();
	212	dataEvent->ResetTrigger();
	213	dataEvent->ResetRecTracks();
	214	dataEvent->ResetRecTriggerTracks();
	215
	216	}
cf464691	217	loader->UnloadRecPoints();
	218	loader->UnloadTracks();
	219
	220	delete recoCluster;
	221	delete recoEvent;
	222	delete trigDec;
	223	}
	224
cac2eb58	225	//_____________________________________________________________________________
	226	void AliMUONReconstructor::FillESD(AliRunLoader* runLoader, AliESD* esd) const
	227	{
	228	TClonesArray* recTracksArray = 0;
	229	TClonesArray* recTrigTracksArray = 0;
	230
	231	AliLoader* loader = runLoader->GetLoader("MUONLoader");
	232	loader->LoadTracks("READ");
cac2eb58	233	AliMUONData* muonData = new AliMUONData(loader,"MUON","MUON");
	234
	235	// declaration
b002ac70	236	Int_t iEvent, nPart;
b002ac70	237	Int_t nTrackHits, nPrimary;
cac2eb58	238	Double_t fitFmin;
b002ac70	239	TArrayF vertex(3);
cac2eb58	240
	241	Double_t bendingSlope, nonBendingSlope, inverseBendingMomentum;
	242	Double_t xRec, yRec, zRec, chi2MatchTrigger;
	243	Bool_t matchTrigger;
	244
b002ac70	245	// setting pointer for tracks, triggertracks & trackparam at vertex
cac2eb58	246	AliMUONTrack* recTrack = 0;
	247	AliMUONTrackParam* trackParam = 0;
	248	AliMUONTriggerTrack* recTriggerTrack = 0;
b002ac70	249	TParticle* particle = new TParticle();
b002ac70	250	AliGenEventHeader* header = 0;
cac2eb58	251	iEvent = runLoader->GetEventNumber();
	252	runLoader->GetEvent(iEvent);
	253
b002ac70	254	// vertex calculation (maybe it exists already somewhere else)
	255	vertex[0] = vertex[1] = vertex[2] = 0.;
	256	nPrimary = 0;
	257	if ( (header = runLoader->GetHeader()->GenEventHeader()) ) {
	258	header->PrimaryVertex(vertex);
	259	} else {
	260	runLoader->LoadKinematics("READ");
	261	runLoader->TreeK()->GetBranch("Particles")->SetAddress(&particle);
	262	nPart = (Int_t)runLoader->TreeK()->GetEntries();
	263	for(Int_t iPart = 0; iPart < nPart; iPart++) {
	264	runLoader->TreeK()->GetEvent(iPart);
	265	if (particle->GetFirstMother() == -1) {
	266	vertex[0] += particle->Vx();
	267	vertex[1] += particle->Vy();
	268	vertex[2] += particle->Vz();
	269	nPrimary++;
	270	}
	271	if (nPrimary) {
	272	vertex[0] /= (double)nPrimary;
	273	vertex[1] /= (double)nPrimary;
	274	vertex[2] /= (double)nPrimary;
	275	}
	276	}
	277	}
cac2eb58	278	// setting ESD MUON class
30178c30	279	AliESDMuonTrack* theESDTrack = new AliESDMuonTrack() ;
cac2eb58	280
	281	//-------------------- trigger tracks-------------
	282	Long_t trigPat = 0;
	283	muonData->SetTreeAddress("RL");
	284	muonData->GetRecTriggerTracks();
	285	recTrigTracksArray = muonData->RecTriggerTracks();
	286
	287	// ready global trigger pattern from first track
	288	if (recTrigTracksArray)
	289	recTriggerTrack = (AliMUONTriggerTrack*) recTrigTracksArray->First();
	290	if (recTriggerTrack) trigPat = recTriggerTrack->GetGTPattern();
	291
	292	//printf(">>> Event %d Number of Recconstructed tracks %d \n",iEvent, nrectracks);
	293
	294	// -------------------- tracks-------------
	295	muonData->SetTreeAddress("RT");
	296	muonData->GetRecTracks();
	297	recTracksArray = muonData->RecTracks();
	298
	299	Int_t nRecTracks = 0;
	300	if (recTracksArray)
	301	nRecTracks = (Int_t) recTracksArray->GetEntriesFast(); //
	302
	303	// loop over tracks
	304	for (Int_t iRecTracks = 0; iRecTracks < nRecTracks; iRecTracks++) {
	305
	306	// reading info from tracks
	307	recTrack = (AliMUONTrack*) recTracksArray->At(iRecTracks);
	308
889a0215	309	trackParam = (AliMUONTrackParam*) (recTrack->GetTrackParamAtHit())->First();
889a0215	310	trackParam->ExtrapToVertex(vertex[0],vertex[1],vertex[2]);
cac2eb58	311
	312	bendingSlope = trackParam->GetBendingSlope();
	313	nonBendingSlope = trackParam->GetNonBendingSlope();
	314	inverseBendingMomentum = trackParam->GetInverseBendingMomentum();
	315	xRec = trackParam->GetNonBendingCoor();
	316	yRec = trackParam->GetBendingCoor();
	317	zRec = trackParam->GetZ();
	318
	319	nTrackHits = recTrack->GetNTrackHits();
	320	fitFmin = recTrack->GetFitFMin();
	321	matchTrigger = recTrack->GetMatchTrigger();
	322	chi2MatchTrigger = recTrack->GetChi2MatchTrigger();
	323
	324	// setting data member of ESD MUON
30178c30	325	theESDTrack->SetInverseBendingMomentum(inverseBendingMomentum);
	326	theESDTrack->SetThetaX(TMath::ATan(nonBendingSlope));
	327	theESDTrack->SetThetaY(TMath::ATan(bendingSlope));
889a0215	328	theESDTrack->SetZ(zRec);
	329	theESDTrack->SetBendingCoor(yRec); // calculate vertex at ESD or Tracking level ?
	330	theESDTrack->SetNonBendingCoor(xRec);
30178c30	331	theESDTrack->SetChi2(fitFmin);
	332	theESDTrack->SetNHit(nTrackHits);
	333	theESDTrack->SetMatchTrigger(matchTrigger);
	334	theESDTrack->SetChi2MatchTrigger(chi2MatchTrigger);
cac2eb58	335
	336	// storing ESD MUON Track into ESD Event
	337	if (nRecTracks != 0)
30178c30	338	esd->AddMuonTrack(theESDTrack);
cac2eb58	339	} // end loop tracks
	340
	341	// add global trigger pattern
	342	if (nRecTracks != 0)
	343	esd->SetTrigger(trigPat);
	344
	345	// reset muondata
	346	muonData->ResetRecTracks();
	347	muonData->ResetRecTriggerTracks();
	348
	349	//} // end loop on event
	350	loader->UnloadTracks();
b002ac70	351	if (!header)
b002ac70	352	runLoader->UnloadKinematics();
30178c30	353	delete theESDTrack;
cac2eb58	354	delete muonData;
b002ac70	355	// delete particle;
cac2eb58	356	}