Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/TreeWriter.cc

    r2600216 r839deb7  
    215215    entry->Pz = momentum.Pz();
    216216
    217     entry->D0            = candidate->D0;
    218     entry->DZ            = candidate->DZ;
    219     entry->P             = candidate->P;
    220     entry->PT            = candidate->PT;
    221     entry->CtgTheta      = candidate->CtgTheta;
    222     entry->Phi           = candidate->Phi;
    223 
    224217    entry->Eta = eta;
    225218    entry->Phi = momentum.Phi();
     
    244237
    245238  const Double_t c_light = 2.99792458E8;
    246  
    247   Double_t x, y, z, t, xError, yError, zError, sigma, sumPT2, btvSumPT2, genDeltaZ, genSumPT2;
    248   UInt_t index, ndf;
    249  
     239
    250240  // loop over all vertices
    251241  iterator.Reset();
    252242  while((candidate = static_cast<Candidate*>(iterator.Next())))
    253243  {
    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 ();
     244    const TLorentzVector &position = candidate->Position;
    271245
    272246    entry = static_cast<Vertex*>(branch->NewEntry());
    273247
    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 
     248    entry->X = position.X();
     249    entry->Y = position.Y();
     250    entry->Z = position.Z();
     251    entry->T = position.T()*1.0E-3/c_light;
     252  }
     253}
    293254
    294255//------------------------------------------------------------------------------
     
    330291    entry->ZOuter = position.Z();
    331292    entry->TOuter = position.T()*1.0E-3/c_light;
    332  
    333     entry->L = candidate->L;
    334  
    335     entry->D0            = candidate->D0;
    336     entry->ErrorD0       = candidate->ErrorD0;
    337     entry->DZ            = candidate->DZ;
    338     entry->ErrorDZ       = candidate->ErrorDZ;
    339     entry->P             = candidate->P;
    340     entry->ErrorP        = candidate->ErrorP;
    341     entry->PT            = candidate->PT;
    342     entry->ErrorPT       = candidate->ErrorPT;
    343     entry->CtgTheta      = candidate->CtgTheta;
    344     entry->ErrorCtgTheta = candidate->ErrorCtgTheta;
    345     entry->Phi           = candidate->Phi;
    346     entry->ErrorPhi      = candidate->ErrorPhi;
    347    
     293
     294    entry->Dxy = candidate->Dxy;
     295    entry->SDxy = candidate->SDxy ;
    348296    entry->Xd = candidate->Xd;
    349297    entry->Yd = candidate->Yd;
     
    359307
    360308    entry->Eta = eta;
    361    
     309    entry->Phi = momentum.Phi();
     310    entry->PT = pt;
     311
    362312    particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));
    363313    const TLorentzVector &initialPosition = particle->Position;
     
    369319
    370320    entry->Particle = particle;
    371 
    372     entry->VertexIndex = candidate->ClusterIndex;
    373 
    374321  }
    375322}
Note: See TracChangeset for help on using the changeset viewer.