Changes in modules/ClusterCounting.cc [6216fb2:d4fb268] in git

File:

• modules/ClusterCounting.cc (modified) (4 diffs)

modules/ClusterCounting.cc

@@ -1 @@
-   /*
+  /*
  *  Delphes: a framework for fast simulation of a generic collider experiment
  *  Copyright (C) 2020  Universite catholique de Louvain (UCLouvain), Belgium
@@ -28 +28 @@
  */
 
+//FIXME add reference to Bedeschi-code
+//FIXME make sure about units of P, X
+//FIXME fix pt > 200 GeV issue and angle > 6.41
+
 #include "modules/ClusterCounting.h"
+
 #include "classes/DelphesClasses.h"
 #include "TrackCovariance/TrkUtil.h"
@@ -103 +108 @@
 {
   Candidate *candidate, *mother, *particle;
-  Double_t mass, trackLength, Ncl;
+  Double_t mass, Ncl;
 
   fItInputArray->Reset();
@@ -119 +124 @@
     mass = candidateMomentum.M();
 
-    trackLength = fTrackUtil->TrkLen(Par);
-
     mother    = candidate;
     candidate = static_cast<Candidate*>(candidate->Clone());
 
-    Ncl = 0.;
-    if (fTrackUtil->IonClusters(Ncl, mass, Par))
+    Ncl = -999;
+    // computation of Nclusters is not supported for electrons
+    if (TMath::Abs(particle->PID) == 11)
     {
       candidate->Nclusters = Ncl;
-      candidate->dNdx = (trackLength > 0.) ? Ncl/trackLength : -1;
     }
+    else if (fTrackUtil->IonClusters(Ncl, mass, Par))
+    {
+      candidate->Nclusters = Ncl;
+    }
+    //cout<<candidate->PID<<", "<<mass<<", "<<candidate->Nclusters<<endl;
 
     candidate->AddCandidate(mother);