- Timestamp:
- May 17, 2016, 3:13:43 PM (9 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 3f7ef54
- Parents:
- 4569b59
- Location:
- modules
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TreeWriter.cc
r4569b59 r2600216 244 244 245 245 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 247 250 // loop over all vertices 248 251 iterator.Reset(); 249 252 while((candidate = static_cast<Candidate*>(iterator.Next()))) 250 253 { 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 (); 252 271 253 272 entry = static_cast<Vertex*>(branch->NewEntry()); 254 273 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 261 293 262 294 //------------------------------------------------------------------------------ … … 337 369 338 370 entry->Particle = particle; 371 372 entry->VertexIndex = candidate->ClusterIndex; 373 339 374 } 340 375 } -
modules/VertexFinder.cc
r4569b59 r2600216 1 1 /** \class VertexFinder 2 2 * 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 9 4 * 10 5 * \author M. Selvaggi - UCL, Louvain-la-Neuve 11 6 * 12 7 */ 8 13 9 14 10 #include "modules/VertexFinder.h" … … 70 66 71 67 fOutputArray = ExportArray(GetString("OutputArray", "tracks")); 72 f ClusterOutputArray = ExportArray(GetString("ClusterOutputArray", "clusters"));68 fVertexOutputArray = ExportArray(GetString("VertexOutputArray", "vertices")); 73 69 } 74 70 … … 183 179 candidate->PositionError.SetXYZT(0.0, 0.0, clusterIDToDouble.at (cluster->first).at ("ez"), 0.0); 184 180 185 f ClusterOutputArray->Add(candidate);181 fVertexOutputArray->Add(candidate); 186 182 } 187 183 //////////////////////////////////////////////////////////////////////////////// -
modules/VertexFinder.h
r4569b59 r2600216 4 4 /** \class VertexFinder 5 5 * 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 12 7 * 13 8 * \author M. Selvaggi - UCL, Louvain-la-Neuve … … 63 58 64 59 TObjArray *fOutputArray; 65 TObjArray *f ClusterOutputArray;60 TObjArray *fVertexOutputArray; 66 61 67 62 map<unsigned, map<string, double> > trackIDToDouble;
Note:
See TracChangeset
for help on using the changeset viewer.