Changes in classes/DelphesSTDHEPReader.cc [341014c:1716933] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
classes/DelphesSTDHEPReader.cc
r341014c r1716933 17 17 */ 18 18 19 19 20 /** \class DelphesSTDHEPReader 20 21 * … … 27 28 #include "classes/DelphesSTDHEPReader.h" 28 29 30 #include <stdexcept> 29 31 #include <iostream> 30 32 #include <sstream> 31 #include <stdexcept> 32 33 34 #include <stdio.h> 33 35 #include <errno.h> 34 36 #include <stdint.h> 35 #include <stdio.h>36 37 #include <string.h> 37 38 39 #include "TObjArray.h" 40 #include "TStopwatch.h" 38 41 #include "TDatabasePDG.h" 42 #include "TParticlePDG.h" 39 43 #include "TLorentzVector.h" 40 #include "TObjArray.h"41 #include "TParticlePDG.h"42 #include "TStopwatch.h"43 44 44 45 #include "classes/DelphesClasses.h" … … 50 51 using namespace std; 51 52 52 static const int kBufferSize = 1000000;53 static const int kBufferSize = 1000000; 53 54 54 55 //--------------------------------------------------------------------------- … … 57 58 fInputFile(0), fBuffer(0), fPDG(0), fBlockType(-1) 58 59 { 59 fBuffer = new uint8_t[kBufferSize *96 + 24];60 fBuffer = new uint8_t[kBufferSize*96 + 24]; 60 61 61 62 fPDG = TDatabasePDG::Instance(); … … 116 117 ReadEventHeader(); 117 118 } 118 else if(fBlockType == MCFIO_STDHEPBEG || fBlockType == MCFIO_STDHEPEND) 119 else if(fBlockType == MCFIO_STDHEPBEG || 120 fBlockType == MCFIO_STDHEPEND) 119 121 { 120 122 ReadSTDCM1(); … … 168 170 uint32_t size; 169 171 fReader[0].ReadValue(&size, 4); 170 SkipBytes(size *elsize);172 SkipBytes(size*elsize); 171 173 } 172 174 … … 176 178 { 177 179 uint32_t i; 178 enum STDHEPVersion 179 { 180 UNKNOWN, 181 V1, 182 V2, 183 V21 184 } version; 180 enum STDHEPVersion {UNKNOWN, V1, V2, V21} version; 185 181 186 182 // version 187 183 fReader[0].ReadString(fBuffer, 100); 188 if(fBuffer[0] == '\0' || fBuffer[1] == '\0') 189 version = UNKNOWN; 190 else if(fBuffer[0] == '1') 191 version = V1; 192 else if(strncmp((char *)fBuffer, "2.01", 4) == 0) 193 version = V21; 194 else if(fBuffer[0] == '2') 195 version = V2; 196 else 197 version = UNKNOWN; 184 if(fBuffer[0] == '\0' || fBuffer[1] == '\0') version = UNKNOWN; 185 else if(fBuffer[0] == '1') version = V1; 186 else if(strncmp((char *)fBuffer, "2.01", 4) == 0) version = V21; 187 else if(fBuffer[0] == '2') version = V2; 188 else version = UNKNOWN; 198 189 199 190 if(version == UNKNOWN) … … 372 363 // 4 + 4 + 4 + 4 + 4 + 4 = 96*n + 24 373 364 374 fReader[0].ReadRaw(fBuffer, 96 *fEventSize + 24);365 fReader[0].ReadRaw(fBuffer, 96*fEventSize + 24); 375 366 376 367 fReader[1].SetBuffer(fBuffer); 377 fReader[2].SetBuffer(fBuffer + 4 * 1 + 4 * 1 *fEventSize);378 fReader[3].SetBuffer(fBuffer + 4 * 2 + 4 * 2 *fEventSize);379 fReader[4].SetBuffer(fBuffer + 4 * 3 + 4 * 4 *fEventSize);380 fReader[5].SetBuffer(fBuffer + 4 * 4 + 4 * 6 *fEventSize);381 fReader[6].SetBuffer(fBuffer + 4 * 5 + 4 * 16 *fEventSize);368 fReader[2].SetBuffer(fBuffer + 4*1 + 4*1*fEventSize); 369 fReader[3].SetBuffer(fBuffer + 4*2 + 4*2*fEventSize); 370 fReader[4].SetBuffer(fBuffer + 4*3 + 4*4*fEventSize); 371 fReader[5].SetBuffer(fBuffer + 4*4 + 4*6*fEventSize); 372 fReader[6].SetBuffer(fBuffer + 4*5 + 4*16*fEventSize); 382 373 383 374 fReader[1].ReadValue(&idhepSize, 4); … … 400 391 fAlphaQCD = 0.0; 401 392 fScaleSize = 0; 402 memset(fScale, 0, 10 *sizeof(double));393 memset(fScale, 0, 10*sizeof(double)); 403 394 } 404 395 … … 503 494 504 495 pdgParticle = fPDG->GetParticle(pid); 505 candidate->Charge = pdgParticle ? int(pdgParticle->Charge() /3.0) : -999;496 candidate->Charge = pdgParticle ? int(pdgParticle->Charge()/3.0) : -999; 506 497 candidate->Mass = mass; 507 498
Note:
See TracChangeset
for help on using the changeset viewer.