1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
17 // MC reader for jet analysis
18 // Author: Mercedes Lopez Noriega (mercedes.lopez.noriega@cern.ch)
21 #include <TClonesArray.h>
23 #include <TParticle.h>
24 #include <TParticlePDG.h>
26 #include <TLorentzVector.h>
30 #include "AliJetMCReader.h"
31 #include "AliJetMCReaderHeader.h"
32 #include "AliESDEvent.h"
33 #include "AliESDtrack.h"
35 ClassImp(AliJetMCReader)
38 AliJetMCReader::AliJetMCReader():
45 //____________________________________________________________________________
47 AliJetMCReader::~AliJetMCReader()
53 //____________________________________________________________________________
56 Bool_t AliJetMCReader::FillMomentumArray(Int_t event)
58 // Fill momentum array
59 TClonesArray &arrayMC = *fArrayMC;
67 // get event from chains
68 fChain->GetEntry(event);
69 fChainMC->GetEntry(event);
70 // get number of tracks in event (for the loop)
71 nt = fESD->GetNumberOfTracks();
73 // get cuts set by user
74 Double_t ptMin = ((AliJetMCReaderHeader*) fReaderHeader)->GetPtCut();
77 for (Int_t it = 0; it < nt; it++) {
78 AliESDtrack *track = fESD->GetTrack(it); //track
79 UInt_t status = track->GetStatus();
80 if ((status & AliESDtrack::kITSrefit) == 0) continue; // quality check
81 // track->GetImpactParameters(dca,z);
82 // if (dca > dcaMax) continue; // check track is reasonable
83 Int_t label = TMath::Abs(track->GetLabel());
84 TParticle *part = (TParticle*)arrayMC[label]; //particle
85 pt = part->Pt(); // pt of the particle
86 if (pt < ptMin) continue; //check cuts
89 // fill momentum array
90 new ((*fMomentumArray)[goodTrack]) TLorentzVector(p.X(), p.Y(), p.Z(), e);
93 printf("\nNumber of good tracks %d \n", goodTrack);