Fork me on GitHub

Changeset 781af69 in git


Ignore:
Timestamp:
May 6, 2021, 10:40:18 AM (4 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
master
Children:
d489660
Parents:
65776c0
Message:

added dNdx member

Files:
4 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesClasses.cc

    r65776c0 r781af69  
    231231  Xd(0), Yd(0), Zd(0),
    232232  Nclusters(0.0),
     233  dNdx(0.0),
    233234  TrackResolution(0),
    234235  NCharged(0),
     
    408409  object.Zd = Zd;
    409410  object.Nclusters = Nclusters;
     411  object.dNdx = dNdx;
    410412  object.TrackResolution = TrackResolution;
    411413  object.NCharged = NCharged;
     
    541543  Zd = 0.0;
    542544  Nclusters = 0.0;
     545  dNdx = 0.0;
    543546  TrackResolution = 0.0;
    544547  NCharged = 0;
  • classes/DelphesClasses.h

    r65776c0 r781af69  
    462462  Float_t DZ; // track longitudinal impact parameter
    463463  Float_t Nclusters; // Number of ionization clusters
     464  Float_t dNdx; // Number of ionization clusters
    464465
    465466  Float_t ErrorP; // track momentum error
     
    567568  Float_t DZ; // track longitudinal impact parameter
    568569  Float_t Nclusters; // Number of ionization clusters
     570  Float_t dNdx; // Number of ionization clusters
    569571
    570572  Float_t ErrorP; // track momentum error
     
    697699
    698700  Float_t Nclusters; // Number of ionization clusters
     701  Float_t dNdx; // Number of ionization clusters per unit length
    699702
    700703  Float_t Xd;
  • modules/ClusterCounting.cc

    r65776c0 r781af69  
    2828 */
    2929
    30 //FIXME add reference to Bedeschi-code
    31 //FIXME make sure about units of P, X
    32 //FIXME fix pt > 200 GeV issue and angle > 6.41
    33 
    3430#include "modules/ClusterCounting.h"
    35 
    3631#include "classes/DelphesClasses.h"
    3732#include "TrackCovariance/TrkUtil.h"
     
    108103{
    109104  Candidate *candidate, *mother, *particle;
    110   Double_t mass, Ncl;
     105  Double_t mass, trackLength, Ncl;
    111106
    112107  fItInputArray->Reset();
     
    124119    mass = candidateMomentum.M();
    125120
     121    trackLength = fTrackUtil->TrkLen(Par);
     122
    126123    mother    = candidate;
    127124    candidate = static_cast<Candidate*>(candidate->Clone());
    128125
    129126    Ncl = -999;
     127
     128
    130129    // computation of Nclusters is not supported for electrons
     130    /*
    131131    if (TMath::Abs(particle->PID) == 11)
    132132    {
    133133      candidate->Nclusters = Ncl;
     134      candidate->dNdx = -999;
    134135    }
    135     else if (fTrackUtil->IonClusters(Ncl, mass, Par))
     136    */
     137    if (fTrackUtil->IonClusters(Ncl, mass, Par))
    136138    {
    137139      candidate->Nclusters = Ncl;
     140      candidate->dNdx = (trackLength > 0.) ? Ncl/trackLength : -1;
    138141    }
    139     //cout<<candidate->PID<<", "<<mass<<", "<<candidate->Nclusters<<endl;
    140142
    141143    candidate->AddCandidate(mother);
  • modules/TreeWriter.cc

    r65776c0 r781af69  
    356356    entry->DZ = candidate->DZ;
    357357    entry->Nclusters = candidate->Nclusters;
     358    entry->dNdx = candidate->dNdx;
    358359
    359360    entry->ErrorP = candidate->ErrorP;
     
    510511    entry->DZ = candidate->DZ;
    511512    entry->Nclusters = candidate->Nclusters;
     513    entry->dNdx = candidate->dNdx;
    512514
    513515    entry->ErrorP = candidate->ErrorP;
Note: See TracChangeset for help on using the changeset viewer.