+Bool_t AliJetKineReader::GetGenJets(AliJet* genJets)
+{
+ // Get generated jets from mc header
+ AliHeader* alih = GetAliHeader();
+ if (alih == 0) return kFALSE;
+ AliGenEventHeader * genh = alih->GenEventHeader();
+ if (genh == 0) return kFALSE;
+ Int_t nj =((AliGenPythiaEventHeader*)genh)->NTriggerJets();
+ Int_t* m = new Int_t[nj];
+ Int_t* k = new Int_t[nj];
+ for (Int_t i=0; i< nj; i++) {
+ Float_t p[4];
+ ((AliGenPythiaEventHeader*)genh)->TriggerJet(i,p);
+ genJets->AddJet(p[0],p[1],p[2],p[3]);
+ m[i]=1;
+ k[i]=i;
+ }
+ genJets->SetNinput(nj);
+ genJets->SetMultiplicities(m);
+ genJets->SetInJet(k);
+ return kTRUE;
+}