Changeset 341014c in git for classes/DelphesLHEFReader.cc
- Timestamp:
- Feb 12, 2019, 9:29:17 PM (6 years ago)
- Branches:
- ImprovedOutputFile, Timing, llp, master
- Children:
- 6455202
- Parents:
- 45e58be
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
classes/DelphesLHEFReader.cc
r45e58be r341014c 17 17 */ 18 18 19 20 19 /** \class DelphesLHEFReader 21 20 * … … 28 27 #include "classes/DelphesLHEFReader.h" 29 28 30 #include <stdexcept>31 29 #include <iostream> 32 30 #include <sstream> 31 #include <stdexcept> 33 32 34 33 #include <stdio.h> 35 34 35 #include "TDatabasePDG.h" 36 #include "TLorentzVector.h" 36 37 #include "TObjArray.h" 38 #include "TParticlePDG.h" 37 39 #include "TStopwatch.h" 38 #include "TDatabasePDG.h"39 #include "TParticlePDG.h"40 #include "TLorentzVector.h"41 40 42 41 #include "classes/DelphesClasses.h" … … 54 53 DelphesLHEFReader::DelphesLHEFReader() : 55 54 fInputFile(0), fBuffer(0), fPDG(0), 56 fEventReady(kFALSE), fEventCounter(-1), fParticleCounter(-1), 55 fEventReady(kFALSE), fEventCounter(-1), fParticleCounter(-1), fCrossSection(1) 57 56 58 57 { … … 124 123 if(!rc) 125 124 { 126 cerr << "** ERROR: " << "invalid event format" << endl; 125 cerr << "** ERROR: " 126 << "invalid event format" << endl; 127 127 return kFALSE; 128 128 } … … 148 148 if(!rc) 149 149 { 150 cerr << "** ERROR: " << "invalid particle format" << endl; 150 cerr << "** ERROR: " 151 << "invalid particle format" << endl; 151 152 return kFALSE; 152 153 } … … 162 163 if(!pch) 163 164 { 164 cerr << "** ERROR: " << "invalid weight format" << endl; 165 cerr << "** ERROR: " 166 << "invalid weight format" << endl; 165 167 return kFALSE; 166 168 } … … 172 174 if(!pch) 173 175 { 174 cerr << "** ERROR: " << "invalid weight format" << endl; 176 cerr << "** ERROR: " 177 << "invalid weight format" << endl; 175 178 return kFALSE; 176 179 } … … 181 184 if(!rc) 182 185 { 183 cerr << "** ERROR: " << "invalid weight format" << endl; 186 cerr << "** ERROR: " 187 << "invalid weight format" << endl; 184 188 return kFALSE; 185 189 } … … 192 196 if(!pch) 193 197 { 194 cerr << "** ERROR: " << "invalid cross section format" << endl; 198 cerr << "** ERROR: " 199 << "invalid cross section format" << endl; 195 200 return kFALSE; 196 201 } … … 199 204 if(!pch) 200 205 { 201 cerr << "** ERROR: " << "invalid cross section format" << endl; 206 cerr << "** ERROR: " 207 << "invalid cross section format" << endl; 202 208 return kFALSE; 203 209 } … … 208 214 if(!rc) 209 215 { 210 cerr << "** ERROR: " << "invalid cross section format" << endl; 216 cerr << "** ERROR: " 217 << "invalid cross section format" << endl; 211 218 return kFALSE; 212 219 } … … 249 256 { 250 257 LHEFWeight *element; 251 vector< pair< int, double >>::const_iterator itWeightList;258 vector<pair<int, double>>::const_iterator itWeightList; 252 259 253 260 for(itWeightList = fWeightList.begin(); itWeightList != fWeightList.end(); ++itWeightList) … … 279 286 280 287 pdgParticle = fPDG->GetParticle(fPID); 281 candidate->Charge = pdgParticle ? int(pdgParticle->Charge() /3.0) : -999;288 candidate->Charge = pdgParticle ? int(pdgParticle->Charge() / 3.0) : -999; 282 289 candidate->Mass = fMass; 283 290
Note:
See TracChangeset
for help on using the changeset viewer.