New quenching algorithm.
[u/mrichter/AliRoot.git] / PYTHIA6 / AliGenPythiaEventHeader.cxx
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 /* $Id$ */
17
18 #include "AliGenPythiaEventHeader.h"
19 ClassImp(AliGenPythiaEventHeader);
20
21
22
23 AliGenPythiaEventHeader::AliGenPythiaEventHeader()
24 {
25 // Default Constructor
26     fNJets = 0;
27 }
28
29 AliGenPythiaEventHeader::AliGenPythiaEventHeader(const char* name):AliGenEventHeader(name)
30 {
31 // Constructor
32     fNJets = 0;
33 }
34
35 void AliGenPythiaEventHeader::AddJet(Float_t px, Float_t py, Float_t pz, Float_t e)
36 {
37 //
38 //  Add a jet 
39 //
40     if (fNJets < 10) {
41         fJets[0][fNJets] = px;
42         fJets[1][fNJets] = py;
43         fJets[2][fNJets] = pz;
44         fJets[3][fNJets] = e;
45         fNJets++;
46     } else {
47         printf("\nWarning: More than 10 jets triggered !!\n");
48     }
49 }
50
51 void AliGenPythiaEventHeader::TriggerJet(Int_t i, Float_t p[4])
52 {
53 //
54 // Give back jet #i
55 //
56     if (i >= fNJets) {
57         printf("\nWarning: TriggerJet, index out of Range!!\n");
58     } else {
59         p[0] = fJets[0][i];
60         p[1] = fJets[1][i];
61         p[2] = fJets[2][i];
62         p[3] = fJets[3][i];
63     }
64 }