Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesSTDHEPReader.cc

    r341014c r1716933  
    1717 */
    1818
     19
    1920/** \class DelphesSTDHEPReader
    2021 *
     
    2728#include "classes/DelphesSTDHEPReader.h"
    2829
     30#include <stdexcept>
    2931#include <iostream>
    3032#include <sstream>
    31 #include <stdexcept>
    32 
     33
     34#include <stdio.h>
    3335#include <errno.h>
    3436#include <stdint.h>
    35 #include <stdio.h>
    3637#include <string.h>
    3738
     39#include "TObjArray.h"
     40#include "TStopwatch.h"
    3841#include "TDatabasePDG.h"
     42#include "TParticlePDG.h"
    3943#include "TLorentzVector.h"
    40 #include "TObjArray.h"
    41 #include "TParticlePDG.h"
    42 #include "TStopwatch.h"
    4344
    4445#include "classes/DelphesClasses.h"
     
    5051using namespace std;
    5152
    52 static const int kBufferSize = 1000000;
     53static const int kBufferSize  = 1000000;
    5354
    5455//---------------------------------------------------------------------------
     
    5758  fInputFile(0), fBuffer(0), fPDG(0), fBlockType(-1)
    5859{
    59   fBuffer = new uint8_t[kBufferSize * 96 + 24];
     60  fBuffer = new uint8_t[kBufferSize*96 + 24];
    6061
    6162  fPDG = TDatabasePDG::Instance();
     
    116117    ReadEventHeader();
    117118  }
    118   else if(fBlockType == MCFIO_STDHEPBEG || fBlockType == MCFIO_STDHEPEND)
     119  else if(fBlockType == MCFIO_STDHEPBEG ||
     120          fBlockType == MCFIO_STDHEPEND)
    119121  {
    120122    ReadSTDCM1();
     
    168170  uint32_t size;
    169171  fReader[0].ReadValue(&size, 4);
    170   SkipBytes(size * elsize);
     172  SkipBytes(size*elsize);
    171173}
    172174
     
    176178{
    177179  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;
    185181
    186182  // version
    187183  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;
    198189
    199190  if(version == UNKNOWN)
     
    372363  // 4 + 4 + 4 + 4 + 4 + 4 = 96*n + 24
    373364
    374   fReader[0].ReadRaw(fBuffer, 96 * fEventSize + 24);
     365  fReader[0].ReadRaw(fBuffer, 96*fEventSize + 24);
    375366
    376367  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);
    382373
    383374  fReader[1].ReadValue(&idhepSize, 4);
     
    400391  fAlphaQCD = 0.0;
    401392  fScaleSize = 0;
    402   memset(fScale, 0, 10 * sizeof(double));
     393  memset(fScale, 0, 10*sizeof(double));
    403394}
    404395
     
    503494
    504495    pdgParticle = fPDG->GetParticle(pid);
    505     candidate->Charge = pdgParticle ? int(pdgParticle->Charge() / 3.0) : -999;
     496    candidate->Charge = pdgParticle ? int(pdgParticle->Charge()/3.0) : -999;
    506497    candidate->Mass = mass;
    507498
Note: See TracChangeset for help on using the changeset viewer.