Fork me on GitHub

Changeset 2600216 in git for modules


Ignore:
Timestamp:
May 17, 2016, 3:13:43 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
3f7ef54
Parents:
4569b59
Message:

added variables for vertexing

Location:
modules
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • modules/TreeWriter.cc

    r4569b59 r2600216  
    244244
    245245  const Double_t c_light = 2.99792458E8;
    246 
     246 
     247  Double_t x, y, z, t, xError, yError, zError, sigma, sumPT2, btvSumPT2, genDeltaZ, genSumPT2;
     248  UInt_t index, ndf;
     249 
    247250  // loop over all vertices
    248251  iterator.Reset();
    249252  while((candidate = static_cast<Candidate*>(iterator.Next())))
    250253  {
    251     const TLorentzVector &position = candidate->Position;
     254   
     255    index = candidate->ClusterIndex;
     256    ndf = candidate->ClusterNDF;
     257    sigma = candidate->ClusterSigma;
     258    sumPT2 = candidate->SumPT2;
     259    btvSumPT2 = candidate->BTVSumPT2;
     260    genDeltaZ = candidate->GenDeltaZ;
     261    genSumPT2 = candidate->GenSumPT2;
     262
     263    x = candidate->Position.X ();
     264    y = candidate->Position.Y ();
     265    z = candidate->Position.Z ();
     266    t = candidate->Position.T()*1.0E-3/c_light;
     267 
     268    xError = candidate->PositionError.X ();
     269    yError = candidate->PositionError.Y ();
     270    zError = candidate->PositionError.Z ();
    252271
    253272    entry = static_cast<Vertex*>(branch->NewEntry());
    254273
    255     entry->X = position.X();
    256     entry->Y = position.Y();
    257     entry->Z = position.Z();
    258     entry->T = position.T()*1.0E-3/c_light;
    259   }
    260 }
     274    entry->Index = index;
     275    entry->NDF = ndf;
     276    entry->Sigma = sigma;
     277    entry->SumPT2 = sumPT2;
     278    entry->BTVSumPT2 = btvSumPT2;
     279    entry->GenDeltaZ = genDeltaZ;
     280    entry->GenSumPT2 = genSumPT2;
     281 
     282    entry->X = x;
     283    entry->Y = y;
     284    entry->Z = z;
     285    entry->T = t;
     286   
     287    entry->ErrorX = xError;
     288    entry->ErrorY = yError;
     289    entry->ErrorZ = zError;
     290  }
     291}
     292
    261293
    262294//------------------------------------------------------------------------------
     
    337369
    338370    entry->Particle = particle;
     371
     372    entry->VertexIndex = candidate->ClusterIndex;
     373
    339374  }
    340375}
  • modules/VertexFinder.cc

    r4569b59 r2600216  
    11/** \class VertexFinder
    22 *
    3  *  Merges particles from pile-up sample into event
    4  *
    5  *
    6  *  $Date: 2013-02-12 15:13:59 +0100 (Tue, 12 Feb 2013) $
    7  *  $Revision: 907 $
    8  *
     3 *  Cluster vertices from tracks
    94 *
    105 *  \author M. Selvaggi - UCL, Louvain-la-Neuve
    116 *
    127 */
     8
    139
    1410#include "modules/VertexFinder.h"
     
    7066
    7167  fOutputArray = ExportArray(GetString("OutputArray", "tracks"));
    72   fClusterOutputArray = ExportArray(GetString("ClusterOutputArray", "clusters"));
     68  fVertexOutputArray = ExportArray(GetString("VertexOutputArray", "vertices"));
    7369}
    7470
     
    183179    candidate->PositionError.SetXYZT(0.0, 0.0, clusterIDToDouble.at (cluster->first).at ("ez"), 0.0);
    184180
    185     fClusterOutputArray->Add(candidate);
     181    fVertexOutputArray->Add(candidate);
    186182  }
    187183////////////////////////////////////////////////////////////////////////////////
  • modules/VertexFinder.h

    r4569b59 r2600216  
    44/** \class VertexFinder
    55 *
    6  *  Merges particles from pile-up sample into event
    7  *
    8  *
    9  *  $Date: 2013-02-12 15:13:59 +0100 (Tue, 12 Feb 2013) $
    10  *  $Revision: 907 $
    11  *
     6 *  Cluster vertices from tracks
    127 *
    138 *  \author M. Selvaggi - UCL, Louvain-la-Neuve
     
    6358
    6459  TObjArray *fOutputArray;
    65   TObjArray *fClusterOutputArray;
     60  TObjArray *fVertexOutputArray;
    6661
    6762  map<unsigned, map<string, double> > trackIDToDouble;
Note: See TracChangeset for help on using the changeset viewer.