Option kPyCharmPbMNR added. Produce charm pairs in agreement with MNR
[u/mrichter/AliRoot.git] / EVGEN / AliGenPythia.cxx
CommitLineData
4c039060 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 **************************************************************************/
15
16/*
17$Log$
fbd1348b 18Revision 1.49 2002/02/08 16:50:50 morsch
19Add name and title in constructor.
20
8b31bfac 21Revision 1.48 2001/12/20 11:44:28 morsch
22Add kinematic bias for direct gamma production.
23
b2bb450d 24Revision 1.47 2001/12/19 14:45:00 morsch
25Store number of trials in header.
26
fd8fb861 27Revision 1.46 2001/12/19 10:36:19 morsch
28Add possibility if jet kinematic biasing.
29
988b118b 30Revision 1.45 2001/11/28 08:06:52 morsch
31Use fMaxLifeTime parameter.
32
8cdfd303 33Revision 1.44 2001/11/27 13:13:07 morsch
34Maximum lifetime for long-lived particles to be put on the stack is parameter.
35It can be set via SetMaximumLifetime(..).
36
47fc6bd5 37Revision 1.43 2001/10/21 18:35:56 hristov
38Several pointers were set to zero in the default constructors to avoid memory management problems
39
2685bf00 40Revision 1.42 2001/10/15 08:21:55 morsch
41Vertex truncation settings moved to AliGenMC.
42
4c07486d 43Revision 1.41 2001/10/08 08:45:42 morsch
44Possibility of vertex cut added.
45
a0b2b296 46Revision 1.40 2001/09/25 11:30:23 morsch
47Pass event vertex to header.
48
c5d36e84 49Revision 1.39 2001/07/27 17:09:36 morsch
50Use local SetTrack, KeepTrack and SetHighWaterMark methods
51to delegate either to local stack or to stack owned by AliRun.
52(Piotr Skowronski, A.M.)
53
a99cf51f 54Revision 1.38 2001/07/13 10:58:54 morsch
55- Some coded moved to AliGenMC
56- Improved handling of secondary vertices.
57
fff02fee 58Revision 1.37 2001/06/28 11:17:28 morsch
59SetEventListRange setter added. Events in specified range are listed for
60debugging. (Yuri Kharlov)
61
11dfaf8d 62Revision 1.36 2001/03/30 07:05:49 morsch
63Final print-out in finish run.
64Write parton system for jet-production (preliminary solution).
65
9be6226f 66Revision 1.35 2001/03/09 13:03:40 morsch
67Process_t and Struc_Func_t moved to AliPythia.h
68
f1a48a38 69Revision 1.34 2001/02/14 15:50:40 hristov
70The last particle in event marked using SetHighWaterMark
71
28adac24 72Revision 1.33 2001/01/30 09:23:12 hristov
73Streamers removed (R.Brun)
74
a8a6107b 75Revision 1.32 2001/01/26 19:55:51 hristov
76Major upgrade of AliRoot code
77
2ab0c725 78Revision 1.31 2001/01/17 10:54:31 hristov
79Better protection against FPE
80
0f9e52a3 81Revision 1.30 2000/12/18 08:55:35 morsch
82Make AliPythia dependent generartors work with new scheme of random number generation
83
3356c022 84Revision 1.29 2000/12/04 11:22:03 morsch
85Init of sRandom as in 1.15
86
097ec7c1 87Revision 1.28 2000/12/02 11:41:39 morsch
88Use SetRandom() to initialize random number generator in constructor.
89
1ed25e11 90Revision 1.27 2000/11/30 20:29:02 morsch
91Initialise static variable sRandom in constructor: sRandom = fRandom;
92
ad1df918 93Revision 1.26 2000/11/30 07:12:50 alibrary
94Introducing new Rndm and QA classes
95
65fb704d 96Revision 1.25 2000/10/18 19:11:27 hristov
97Division by zero fixed
98
919c2096 99Revision 1.24 2000/09/18 10:41:35 morsch
100Add possibility to use nuclear structure functions from PDF library V8.
101
811826d8 102Revision 1.23 2000/09/14 14:05:40 morsch
103dito
104
819f84ad 105Revision 1.22 2000/09/14 14:02:22 morsch
106- Correct conversion from mm to cm when passing particle vertex to MC.
107- Correct handling of fForceDecay == all.
108
4451ef92 109Revision 1.21 2000/09/12 14:14:55 morsch
110Call fDecayer->ForceDecay() at the beginning of Generate().
111
00d6ce7d 112Revision 1.20 2000/09/06 14:29:33 morsch
113Use AliPythia for event generation an AliDecayPythia for decays.
114Correct handling of "nodecay" option
115
18edb254 116Revision 1.19 2000/07/11 18:24:56 fca
117Coding convention corrections + few minor bug fixes
118
aee8290b 119Revision 1.18 2000/06/30 12:40:34 morsch
120Pythia takes care of vertex smearing. Correct conversion from Pythia units (mm) to
121Geant units (cm).
122
1512e357 123Revision 1.17 2000/06/09 20:34:07 morsch
124All coding rule violations except RS3 corrected
125
f87cfe57 126Revision 1.16 2000/05/15 15:04:20 morsch
127The full event is written for fNtrack = -1
128Coding rule violations corrected.
129
5ddeb374 130Revision 1.15 2000/04/26 10:14:24 morsch
131Particles array has one entry more than pythia particle list. Upper bound of
132particle loop changed to np-1 (R. Guernane, AM)
133
2a0e6f5b 134Revision 1.14 2000/04/05 08:36:13 morsch
135Check status code of particles in Pythia event
136to avoid double counting as partonic state and final state particle.
137
23211d3e 138Revision 1.13 1999/11/09 07:38:48 fca
139Changes for compatibility with version 2.23 of ROOT
140
084c1b4a 141Revision 1.12 1999/11/03 17:43:20 fca
142New version from G.Martinez & A.Morsch
143
886b6f73 144Revision 1.11 1999/09/29 09:24:14 fca
145Introduction of the Copyright and cvs Log
4c039060 146*/
147
fe4da5cc 148#include "AliGenPythia.h"
fff02fee 149#include "AliGenPythiaEventHeader.h"
18edb254 150#include "AliDecayerPythia.h"
fe4da5cc 151#include "AliRun.h"
152#include "AliPythia.h"
18edb254 153#include "AliPDG.h"
1578254f 154#include <TParticle.h>
18edb254 155#include <TSystem.h>
f87cfe57 156
fe4da5cc 157 ClassImp(AliGenPythia)
158
159AliGenPythia::AliGenPythia()
fff02fee 160 :AliGenMC()
fe4da5cc 161{
18edb254 162// Default Constructor
2685bf00 163 fParticles = 0;
164 fPythia = 0;
988b118b 165 fDecayer = new AliDecayerPythia();
11dfaf8d 166 SetEventListRange();
988b118b 167 SetJetPhiRange();
168 SetJetEtaRange();
fe4da5cc 169}
170
171AliGenPythia::AliGenPythia(Int_t npart)
fff02fee 172 :AliGenMC(npart)
fe4da5cc 173{
174// default charm production at 5. 5 TeV
175// semimuonic decay
176// structure function GRVHO
177//
8b31bfac 178 fName = "Pythia";
179 fTitle= "Particle Generator using PYTHIA";
fe4da5cc 180 fXsection = 0.;
811826d8 181 fNucA1=0;
182 fNucA2=0;
fe4da5cc 183 SetProcess();
184 SetStrucFunc();
886b6f73 185 SetForceDecay();
fe4da5cc 186 SetPtHard();
187 SetEnergyCMS();
18edb254 188 fDecayer = new AliDecayerPythia();
1ed25e11 189 // Set random number generator
097ec7c1 190 sRandom=fRandom;
fff02fee 191 fFlavorSelect = 0;
192 // Produced particles
193 fParticles = new TClonesArray("TParticle",1000);
c5d36e84 194 fEventVertex.Set(3);
988b118b 195 SetEventListRange();
196 SetJetPhiRange();
197 SetJetEtaRange();
fe4da5cc 198}
199
f87cfe57 200AliGenPythia::AliGenPythia(const AliGenPythia & Pythia)
201{
202// copy constructor
203}
204
fe4da5cc 205AliGenPythia::~AliGenPythia()
206{
5ddeb374 207// Destructor
fe4da5cc 208}
209
11dfaf8d 210void AliGenPythia::SetEventListRange(Int_t eventFirst, Int_t eventLast)
211{
212 // Set a range of event numbers, for which a table
213 // of generated particle will be printed
214 fDebugEventFirst = eventFirst;
215 fDebugEventLast = eventLast;
216 if (fDebugEventLast==-1) fDebugEventLast=fDebugEventFirst;
217}
218
fe4da5cc 219void AliGenPythia::Init()
220{
5ddeb374 221// Initialisation
18edb254 222 SetMC(AliPythia::Instance());
fe4da5cc 223 fPythia=(AliPythia*) fgMCEvGen;
224//
225 fParentWeight=1./Float_t(fNpart);
226//
227// Forward Paramters to the AliPythia object
00d6ce7d 228 fDecayer->SetForceDecay(fForceDecay);
18edb254 229 fDecayer->Init();
00d6ce7d 230
18edb254 231
fe4da5cc 232 fPythia->SetCKIN(3,fPtHardMin);
233 fPythia->SetCKIN(4,fPtHardMax);
811826d8 234 if (fNucA1 > 0 && fNucA2 > 0) fPythia->SetNuclei(fNucA1, fNucA2);
fe4da5cc 235 fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc);
18edb254 236
237 // fPythia->Pylist(0);
238 // fPythia->Pystat(2);
fe4da5cc 239// Parent and Children Selection
240 switch (fProcess)
241 {
f1a48a38 242 case kPyCharm:
f1a48a38 243 case kPyCharmUnforced:
fbd1348b 244 case kPyCharmPbMNR:
b2bb450d 245 fParentSelect[0] = 411;
246 fParentSelect[1] = 421;
247 fParentSelect[2] = 431;
248 fParentSelect[3] = 4122;
8b31bfac 249 fFlavorSelect = 4;
fe4da5cc 250 break;
f1a48a38 251 case kPyBeauty:
fff02fee 252 fParentSelect[0]= 511;
253 fParentSelect[1]= 521;
254 fParentSelect[2]= 531;
255 fParentSelect[3]= 5122;
256 fParentSelect[4]= 5132;
257 fParentSelect[5]= 5232;
258 fParentSelect[6]= 5332;
259 fFlavorSelect = 5;
fe4da5cc 260 break;
f1a48a38 261 case kPyBeautyUnforced:
fff02fee 262 fParentSelect[0] = 511;
263 fParentSelect[1] = 521;
264 fParentSelect[2] = 531;
265 fParentSelect[3] = 5122;
266 fParentSelect[4] = 5132;
267 fParentSelect[5] = 5232;
268 fParentSelect[6] = 5332;
269 fFlavorSelect = 5;
fe4da5cc 270 break;
f1a48a38 271 case kPyJpsiChi:
272 case kPyJpsi:
fff02fee 273 fParentSelect[0] = 443;
fe4da5cc 274 break;
f1a48a38 275 case kPyMb:
276 case kPyJets:
277 case kPyDirectGamma:
5be1fe76 278 break;
fe4da5cc 279 }
fff02fee 280 AliGenMC::Init();
fe4da5cc 281}
282
283void AliGenPythia::Generate()
284{
5ddeb374 285// Generate one event
00d6ce7d 286 fDecayer->ForceDecay();
a0b2b296 287
9be6226f 288 Float_t polar[3] = {0,0,0};
289 Float_t origin[3] = {0,0,0};
fff02fee 290 Float_t p[3];
5be1fe76 291// converts from mm/c to s
292 const Float_t kconv=0.001/2.999792458e8;
fe4da5cc 293//
fe4da5cc 294 Int_t nt=0;
295 Int_t jev=0;
7921755b 296 Int_t j, kf;
fe4da5cc 297 fTrials=0;
fff02fee 298
fd8fb861 299 // Set collision vertex position
aee8290b 300 if(fVertexSmear==kPerEvent) {
1512e357 301 fPythia->SetMSTP(151,1);
fe4da5cc 302 for (j=0;j<3;j++) {
fff02fee 303 fPythia->SetPARP(151+j, fOsigma[j]*10.);
fe4da5cc 304 }
aee8290b 305 } else if (fVertexSmear==kPerTrack) {
fe4da5cc 306 fPythia->SetMSTP(151,0);
fe4da5cc 307 }
fff02fee 308// event loop
fe4da5cc 309 while(1)
310 {
084c1b4a 311 fPythia->Pyevnt();
11dfaf8d 312 if (gAlice->GetEvNumber()>=fDebugEventFirst &&
313 gAlice->GetEvNumber()<=fDebugEventLast) fPythia->Pylist(1);
fe4da5cc 314 fTrials++;
47fc6bd5 315
fff02fee 316 fPythia->ImportParticles(fParticles,"All");
317
318//
319//
320//
321 Int_t i;
322
323 Int_t np = fParticles->GetEntriesFast();
a0b2b296 324 if (np == 0 ) continue;
325// Get event vertex and discard the event if the Z coord. is too big
4c07486d 326 TParticle *iparticle = (TParticle *) fParticles->At(0);
327 Float_t distz = iparticle->Vz()/10.;
328 if(TMath::Abs(distz)>fCutVertexZ*fOsigma[2]) continue;
a0b2b296 329//
4c07486d 330 fEventVertex[0] = iparticle->Vx()/10.+fOrigin.At(0);
331 fEventVertex[1] = iparticle->Vy()/10.+fOrigin.At(1);
332 fEventVertex[2] = iparticle->Vz()/10.+fOrigin.At(2);
a0b2b296 333//
fff02fee 334 Int_t* pParent = new Int_t[np];
335 Int_t* pSelected = new Int_t[np];
336 Int_t* trackIt = new Int_t[np];
337 for (i=0; i< np-1; i++) {
338 pParent[i] = -1;
339 pSelected[i] = 0;
340 }
fe4da5cc 341 printf("\n **************************************************%d\n",np);
fff02fee 342 Int_t nc = 0;
f1a48a38 343 if (fProcess != kPyMb && fProcess != kPyJets && fProcess != kPyDirectGamma) {
fd8fb861 344
fff02fee 345 for (i = 0; i<np-1; i++) {
c5d36e84 346 iparticle = (TParticle *) fParticles->At(i);
23211d3e 347 Int_t ks = iparticle->GetStatusCode();
a8228d85 348 kf = CheckPDGCode(iparticle->GetPdgCode());
fff02fee 349// No initial state partons
2a0e6f5b 350 if (ks==21) continue;
fe4da5cc 351//
fff02fee 352// Heavy Flavor Selection
5be1fe76 353//
fff02fee 354 // quark ?
355 kf = TMath::Abs(kf);
356 Int_t kfl = kf;
357 // meson ?
358 if (kfl > 10) kfl/=100;
359 // baryon
360 if (kfl > 10) kfl/=10;
361 if (kfl > 10) kfl/=10;
362
363 Int_t ipa = iparticle->GetFirstMother()-1;
364 Int_t kfMo = 0;
365
366 if (ipa > -1) {
367 TParticle * mother = (TParticle *) fParticles->At(ipa);
368 kfMo = TMath::Abs(mother->GetPdgCode());
369 }
370// printf("\n particle (all) %d %d %d", i, pSelected[i], kf);
371 if (kfl >= fFlavorSelect) {
5be1fe76 372//
fff02fee 373// Heavy flavor hadron or quark
fe4da5cc 374//
fff02fee 375// Kinematic seletion on final state heavy flavor mesons
376 if (ParentSelected(kf) && !KinematicSelection(iparticle, 0))
377 {
378 nc = -1;
379 break;
380 }
381 pSelected[i] = 1;
382// printf("\n particle (HF) %d %d %d", i, pSelected[i], kf);
383 } else {
384// Kinematic seletion on decay products
385 if (fCutOnChild && ParentSelected(kfMo) && ChildSelected(kf)
386 && !KinematicSelection(iparticle, 1))
387 {
388 nc = -1;
389 break;
390 }
fe4da5cc 391//
fff02fee 392// Decay products
393// Select if mother was selected and is not tracked
394
395 if (pSelected[ipa] &&
396 !trackIt[ipa] && // mother will be tracked ?
397 kfMo != 5 && // mother is b-quark, don't store fragments
398 kfMo != 4 && // mother is c-quark, don't store fragments
399 kf != 92) // don't store string
400 {
fe4da5cc 401//
fff02fee 402// Semi-stable or de-selected: diselect decay products:
fe4da5cc 403//
5be1fe76 404//
47fc6bd5 405 if (pSelected[i] == -1 || fDecayer->GetLifetime(kf) > fMaxLifeTime)
fff02fee 406 {
407 Int_t ipF = iparticle->GetFirstDaughter();
408 Int_t ipL = iparticle->GetLastDaughter();
409 if (ipF > 0) for (j = ipF-1; j < ipL; j++) pSelected[j] = -1;
410 }
411// printf("\n particle (decay) %d %d %d", i, pSelected[i], kf);
412 pSelected[i] = (pSelected[i] == -1) ? 0 : 1;
413 }
414 }
415 if (pSelected[i] == -1) pSelected[i] = 0;
416 if (!pSelected[i]) continue;
417 nc++;
418// Decision on tracking
419 trackIt[i] = 0;
420//
421// Track final state particle
422 if (ks == 1) trackIt[i] = 1;
423// Track semi-stable particles
8cdfd303 424 if ((ks ==1) || (fDecayer->GetLifetime(kf) > fMaxLifeTime)) trackIt[i] = 1;
fff02fee 425// Track particles selected by process if undecayed.
426 if (fForceDecay == kNoDecay) {
427 if (ParentSelected(kf)) trackIt[i] = 1;
428 } else {
429 if (ParentSelected(kf)) trackIt[i] = 0;
430 }
431//
432//
433
434 } // particle selection loop
435 if (nc > -1) {
436 for (i = 0; i<np-1; i++) {
437 if (!pSelected[i]) continue;
438 TParticle * iparticle = (TParticle *) fParticles->At(i);
439 kf = CheckPDGCode(iparticle->GetPdgCode());
a0b2b296 440 p[0] = iparticle->Px();
441 p[1] = iparticle->Py();
442 p[2] = iparticle->Pz();
443 origin[0] = fOrigin[0]+iparticle->Vx()/10.;
444 origin[1] = fOrigin[1]+iparticle->Vy()/10.;
445 origin[2] = fOrigin[2]+iparticle->Vz()/10.;
446 Float_t tof = kconv*iparticle->T();
447 Int_t ipa = iparticle->GetFirstMother()-1;
fff02fee 448 Int_t iparent = (ipa > -1) ? pParent[ipa] : -1;
a99cf51f 449 SetTrack(fTrackIt*trackIt[i] ,
fff02fee 450 iparent, kf, p, origin, polar, tof, kPPrimary, nt, 1.);
451 pParent[i] = nt;
a99cf51f 452 KeepTrack(nt);
fff02fee 453 } // SetTrack loop
454 }
455 } else {
456 nc = GenerateMB();
5be1fe76 457 } // mb ?
fff02fee 458
fe4da5cc 459 if (nc > 0) {
5be1fe76 460 jev+=nc;
5ddeb374 461 if (jev >= fNpart || fNpart == -1) {
fe4da5cc 462 fKineBias=Float_t(fNpart)/Float_t(fTrials);
23211d3e 463 printf("\n Trials: %i %i %i\n",fTrials, fNpart, jev);
fd8fb861 464 MakeHeader();
fe4da5cc 465 break;
466 }
467 }
468 } // event loop
a99cf51f 469 SetHighWaterMark(nt);
fe4da5cc 470// adjust weight due to kinematic selection
471 AdjustWeights();
472// get cross-section
473 fXsection=fPythia->GetPARI(1);
474}
475
fff02fee 476Int_t AliGenPythia::GenerateMB()
fe4da5cc 477{
fff02fee 478 Int_t i, kf, nt, iparent;
479 Int_t nc = 0;
480 Float_t p[3];
481 Float_t polar[3] = {0,0,0};
482 Float_t origin[3] = {0,0,0};
fff02fee 483// converts from mm/c to s
484 const Float_t kconv=0.001/2.999792458e8;
4451ef92 485
fff02fee 486 Int_t np = fParticles->GetEntriesFast();
487 Int_t* pParent = new Int_t[np];
488 for (i=0; i< np-1; i++) pParent[i] = -1;
b2bb450d 489 if (fProcess == kPyJets || fProcess == kPyDirectGamma) {
988b118b 490 TParticle* jet1 = (TParticle *) fParticles->At(6);
491 TParticle* jet2 = (TParticle *) fParticles->At(7);
492 if (!CheckTrigger(jet1, jet2)) return 0;
493 }
494
fff02fee 495 for (i = 0; i<np-1; i++) {
496 Int_t trackIt = 0;
497 TParticle * iparticle = (TParticle *) fParticles->At(i);
498 kf = CheckPDGCode(iparticle->GetPdgCode());
499 Int_t ks = iparticle->GetStatusCode();
500 Int_t km = iparticle->GetFirstMother();
501// printf("\n Particle: %d %d %d", i, kf, ks);
502
503 if ((ks == 1 && kf!=0 && KinematicSelection(iparticle, 0)) ||
504 (ks != 1) ||
505 (fProcess == kPyJets && ks == 21 && km == 0 && i>1)) {
506 nc++;
507 if (ks == 1) trackIt = 1;
508 Int_t ipa = iparticle->GetFirstMother()-1;
fe4da5cc 509
fff02fee 510 iparent = (ipa > -1) ? pParent[ipa] : -1;
fe4da5cc 511
512//
fff02fee 513// store track information
a0b2b296 514 p[0] = iparticle->Px();
515 p[1] = iparticle->Py();
516 p[2] = iparticle->Pz();
517 origin[0] = fOrigin[0]+iparticle->Vx()/10.;
518 origin[1] = fOrigin[1]+iparticle->Vy()/10.;
519 origin[2] = fOrigin[2]+iparticle->Vz()/10.;
fff02fee 520 Float_t tof=kconv*iparticle->T();
a99cf51f 521 SetTrack(fTrackIt*trackIt, iparent, kf, p, origin, polar,
a0b2b296 522 tof, kPPrimary, nt);
a99cf51f 523 KeepTrack(nt);
fff02fee 524 pParent[i] = nt;
525 } // select particle
526 } // particle loop
527
528 if (pParent) delete[] pParent;
fe4da5cc 529
fff02fee 530 printf("\n I've put %i particles on the stack \n",nc);
fff02fee 531 return nc;
532}
fe4da5cc 533
fe4da5cc 534
fff02fee 535void AliGenPythia::FinishRun()
536{
537// Print x-section summary
538 fPythia->Pystat(1);
fe4da5cc 539}
fff02fee 540
fe4da5cc 541void AliGenPythia::AdjustWeights()
542{
5ddeb374 543// Adjust the weights after generation of all events
544//
5ddeb374 545 TParticle *part;
fe4da5cc 546 Int_t ntrack=gAlice->GetNtrack();
547 for (Int_t i=0; i<ntrack; i++) {
2ab0c725 548 part= gAlice->Particle(i);
5ddeb374 549 part->SetWeight(part->GetWeight()*fKineBias);
fe4da5cc 550 }
551}
811826d8 552
553void AliGenPythia::SetNuclei(Int_t a1, Int_t a2)
554{
555// Treat protons as inside nuclei with mass numbers a1 and a2
556 fNucA1 = a1;
557 fNucA2 = a2;
558}
fff02fee 559
560
561void AliGenPythia::MakeHeader()
562{
563// Builds the event header, to be called after each event
564 AliGenEventHeader* header = new AliGenPythiaEventHeader("Pythia");
c5d36e84 565 ((AliGenPythiaEventHeader*) header)->SetProcessType(fPythia->GetMSTI(1));
fd8fb861 566 ((AliGenPythiaEventHeader*) header)->SetTrials(fTrials);
c5d36e84 567 header->SetPrimaryVertex(fEventVertex);
568 gAlice->SetGenEventHeader(header);
fff02fee 569}
811826d8 570
988b118b 571
572Bool_t AliGenPythia::CheckTrigger(TParticle* jet1, TParticle* jet2)
573{
574// Check the kinematic trigger condition
575//
b2bb450d 576 Double_t eta[2];
577 eta[0] = jet1->Eta();
578 eta[1] = jet2->Eta();
579 Double_t phi[2];
580 phi[0] = jet1->Phi();
581 phi[1] = jet2->Phi();
582 Int_t pdg[2];
583 pdg[0] = jet1->GetPdgCode();
584 pdg[1] = jet2->GetPdgCode();
988b118b 585 Bool_t triggered = kFALSE;
b2bb450d 586
587 if (fProcess == kPyJets) {
588 //Check eta range first...
589 if ((eta[0] < fEtaMaxJet && eta[0] > fEtaMinJet) ||
590 (eta[1] < fEtaMaxJet && eta[1] > fEtaMinJet))
591 {
592 //Eta is okay, now check phi range
593 if ((phi[0] < fPhiMaxJet && phi[0] > fPhiMinJet) ||
594 (phi[1] < fPhiMaxJet && phi[1] > fPhiMinJet))
595 {
596 triggered = kTRUE;
597 }
598 }
599 } else {
600 Int_t ij = 0;
601 Int_t ig = 1;
602 if (pdg[0] == kGamma) {
603 ij = 1;
604 ig = 0;
605 }
606 //Check eta range first...
607 if ((eta[ij] < fEtaMaxJet && eta[ij] > fEtaMinJet) &&
608 (eta[ig] < fEtaMaxGamma && eta[ig] > fEtaMinGamma))
609 {
610 //Eta is okay, now check phi range
611 if ((phi[ij] < fPhiMaxJet && phi[ij] > fPhiMinJet) &&
612 (phi[ig] < fPhiMaxGamma && phi[ig] > fPhiMinGamma))
613 {
614 triggered = kTRUE;
615 }
616 }
988b118b 617 }
b2bb450d 618
619
988b118b 620 return triggered;
621}
811826d8 622
f87cfe57 623AliGenPythia& AliGenPythia::operator=(const AliGenPythia& rhs)
624{
625// Assignment operator
626 return *this;
627}
fe4da5cc 628
629
811826d8 630
a8a6107b 631#ifdef never
18edb254 632void AliGenPythia::Streamer(TBuffer &R__b)
633{
634 // Stream an object of class AliGenPythia.
635
636 if (R__b.IsReading()) {
637 Version_t R__v = R__b.ReadVersion(); if (R__v) { }
638 AliGenerator::Streamer(R__b);
639 R__b >> (Int_t&)fProcess;
640 R__b >> (Int_t&)fStrucFunc;
641 R__b >> (Int_t&)fForceDecay;
642 R__b >> fEnergyCMS;
643 R__b >> fKineBias;
644 R__b >> fTrials;
645 fParentSelect.Streamer(R__b);
646 fChildSelect.Streamer(R__b);
647 R__b >> fXsection;
648// (AliPythia::Instance())->Streamer(R__b);
649 R__b >> fPtHardMin;
650 R__b >> fPtHardMax;
651// if (fDecayer) fDecayer->Streamer(R__b);
652 } else {
653 R__b.WriteVersion(AliGenPythia::IsA());
654 AliGenerator::Streamer(R__b);
655 R__b << (Int_t)fProcess;
656 R__b << (Int_t)fStrucFunc;
657 R__b << (Int_t)fForceDecay;
658 R__b << fEnergyCMS;
659 R__b << fKineBias;
660 R__b << fTrials;
661 fParentSelect.Streamer(R__b);
662 fChildSelect.Streamer(R__b);
663 R__b << fXsection;
664// R__b << fPythia;
665 R__b << fPtHardMin;
666 R__b << fPtHardMax;
667 // fDecayer->Streamer(R__b);
668 }
669}
a8a6107b 670#endif
18edb254 671