Fork me on GitHub

Changeset 1367 in svn for trunk/modules


Ignore:
Timestamp:
Apr 16, 2014, 4:08:33 PM (11 years ago)
Author:
Michele Selvaggi
Message:

track counting btagging

Location:
trunk/modules
Files:
4 added
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/modules/ModulesLinkDef.h

    r1360 r1367  
    1919#include "modules/EnergySmearing.h"
    2020#include "modules/MomentumSmearing.h"
     21#include "modules/ImpactParameterSmearing.h"
    2122#include "modules/TimeSmearing.h"
    2223#include "modules/Calorimeter.h"
     
    2425#include "modules/EnergyScale.h"
    2526#include "modules/UniqueObjectFinder.h"
     27#include "modules/TrackCountingBTagging.h"
    2628#include "modules/BTagging.h"
    2729#include "modules/TauTagging.h"
     
    5355#pragma link C++ class EnergySmearing+;
    5456#pragma link C++ class MomentumSmearing+;
     57#pragma link C++ class ImpactParameterSmearing+;
    5558#pragma link C++ class TimeSmearing+;
    5659#pragma link C++ class Calorimeter+;
     
    5861#pragma link C++ class EnergyScale+;
    5962#pragma link C++ class UniqueObjectFinder+;
     63#pragma link C++ class TrackCountingBTagging+;
    6064#pragma link C++ class BTagging+;
    6165#pragma link C++ class TauTagging+;
  • trunk/modules/ParticlePropagator.cc

    r1342 r1367  
    106106  Double_t tmp, discr, discr2;
    107107  Double_t delta, gammam, omega, asinrho;
    108 
     108  Double_t ang_mom, rcu, rc2, dxy, xd, yd, zd;
     109 
    109110  const Double_t c_light = 2.99792458E8;
    110111
     
    214215      if(x_c < 0.0) phi += TMath::Pi();
    215216
     217      rcu = TMath::Abs(r);
     218      rc2 = r_c*r_c;
     219     
     220      // calculate coordinates of closest approach to track circle in transverse plane xd, yd, zd
     221      xd = x_c*x_c*x_c - x_c*rcu*r_c + x_c*y_c*y_c;
     222      xd  = ( rc2 > 0.0 ) ? xd / rc2 : -999;
     223      yd  = y_c*(-rcu*r_c + rc2);
     224      yd  = ( rc2 > 0.0 ) ? yd / rc2 : -999;
     225      zd  = z + (TMath::Sqrt(xd*xd+yd*yd) - TMath::Sqrt(x*x+y*y))*pz/pt;
     226
     227      // calculate impact paramater
     228      ang_mom = (xd*py - yd*px);
     229      dxy = ang_mom/pt;
     230   
     231         
    216232      // 3. time evaluation t = TMath::Min(t_r, t_z)
    217233      //    t_r : time to exit from the sides
     
    267283
    268284        candidate->Momentum = candidateMomentum;
    269         candidate->AddCandidate(mother);
     285        candidate->Xd = xd*1.0E3;
     286        candidate->Yd = yd*1.0E3;
     287        candidate->Zd = zd*1.0E3;
     288       
     289        candidate->AddCandidate(mother);
    270290
    271291        fOutputArray->Add(candidate);
     
    287307
    288308//------------------------------------------------------------------------------
     309
  • trunk/modules/TreeWriter.cc

    r1366 r1367  
    277277    entry->ZOuter = position.Z();
    278278    entry->TOuter = position.T()*1.0E-3/c_light;
    279 
     279   
     280    entry->Dxy = candidate->Dxy;
     281    entry->SDxy = candidate->SDxy ;
     282    entry->Xd = candidate -> Xd;
     283    entry->Yd = candidate -> Yd;
     284    entry->Zd = candidate -> Zd;
     285   
    280286    const TLorentzVector &momentum = candidate->Momentum;
    281287
Note: See TracChangeset for help on using the changeset viewer.