00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 inline void MyTruthParticle::init(const TLorentzVector & p,
00016 const double & charge,
00017 const int & pdgId,
00018 const int & barcode,
00019 const int & status,
00020 const int & nParents,
00021 const int & nDaughters,
00022 const int & mother0_barcode,
00023 const int & mother0_pdgId,
00024 const int & index
00025 ) {
00026
00027 m_tlv = p;
00028 m_charge = charge;
00029 m_pdgId = pdgId;
00030 m_barcode = barcode;
00031 m_status = status;
00032 m_nParents = nParents;
00033 m_nDaughters = nDaughters;
00034 m_mother0_barcode = mother0_barcode;
00035 m_mother0_pdgId = mother0_pdgId;
00036 m_index = index;
00037 m_type = "truth_particle";
00038 m_matchingflag = 0;
00039
00040 for (unsigned int k=0; k<5; k++) {
00041 m_Et_cone[k] = 0.0;
00042 }
00043
00044 return;
00045
00046 }
00047
00048
00049
00050
00051
00052
00053
00054 inline int MyTruthParticle::barcode(void) const {
00055
00056 return m_barcode;
00057
00058 }
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078 inline int MyTruthParticle::status(void) const {
00079
00080 return m_status;
00081
00082 }
00083
00084
00085
00086
00087
00088
00089
00090 inline int MyTruthParticle::nParents(void) const {
00091
00092 return m_nParents;
00093
00094 }
00095
00096
00097
00098
00099
00100
00101
00102 inline int MyTruthParticle::nDaughters(void) const {
00103
00104 return m_nDaughters;
00105
00106 }
00107
00108
00109
00110
00111
00112
00113
00114 inline int MyTruthParticle::mother0_barcode(void) const {
00115
00116 return m_mother0_barcode;
00117
00118 }
00119
00120
00121
00122
00123
00124
00125
00126 inline int MyTruthParticle::mother0_pdgId(void) const {
00127
00128 return m_mother0_pdgId;
00129
00130 }
00131
00132
00133
00134
00135
00136
00137
00138 inline double MyTruthParticle::Et_in_cone(const int & delta_R_index) const {
00139
00140
00141
00142 if (delta_R_index<0 || delta_R_index>4) {
00143 std::cerr << "\n"
00144 << "Class MyTruthParticle, method Et_in_cone: ERROR!\n"
00145 << "Wrong delta_R_index " << delta_R_index
00146 << ", must be >=0 and <=7.\n";
00147 exit(1);
00148 }
00149
00150
00151 return m_Et_cone[delta_R_index];
00152
00153 }
00154
00155
00156
00157
00158
00159
00160
00161 inline void MyTruthParticle::set_Et_in_cone(const double Et_cone[5]) {
00162
00163 for (unsigned int k=0; k<5; k++) {
00164 m_Et_cone[k] = Et_cone[k];
00165 }
00166
00167 return;
00168
00169 }
00170
00171
00172
00173
00174
00175
00176
00177 inline void MyTruthParticle::PrintParticle(std::string option) {
00178
00179 return;
00180 }