- Timestamp:
- May 6, 2021, 10:40:18 AM (3 years ago)
- Branches:
- master
- Children:
- d489660
- Parents:
- 65776c0
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/ClusterCounting.cc
r65776c0 r781af69 28 28 */ 29 29 30 //FIXME add reference to Bedeschi-code31 //FIXME make sure about units of P, X32 //FIXME fix pt > 200 GeV issue and angle > 6.4133 34 30 #include "modules/ClusterCounting.h" 35 36 31 #include "classes/DelphesClasses.h" 37 32 #include "TrackCovariance/TrkUtil.h" … … 108 103 { 109 104 Candidate *candidate, *mother, *particle; 110 Double_t mass, Ncl;105 Double_t mass, trackLength, Ncl; 111 106 112 107 fItInputArray->Reset(); … … 124 119 mass = candidateMomentum.M(); 125 120 121 trackLength = fTrackUtil->TrkLen(Par); 122 126 123 mother = candidate; 127 124 candidate = static_cast<Candidate*>(candidate->Clone()); 128 125 129 126 Ncl = -999; 127 128 130 129 // computation of Nclusters is not supported for electrons 130 /* 131 131 if (TMath::Abs(particle->PID) == 11) 132 132 { 133 133 candidate->Nclusters = Ncl; 134 candidate->dNdx = -999; 134 135 } 135 else if (fTrackUtil->IonClusters(Ncl, mass, Par)) 136 */ 137 if (fTrackUtil->IonClusters(Ncl, mass, Par)) 136 138 { 137 139 candidate->Nclusters = Ncl; 140 candidate->dNdx = (trackLength > 0.) ? Ncl/trackLength : -1; 138 141 } 139 //cout<<candidate->PID<<", "<<mass<<", "<<candidate->Nclusters<<endl;140 142 141 143 candidate->AddCandidate(mother); -
modules/TreeWriter.cc
r65776c0 r781af69 356 356 entry->DZ = candidate->DZ; 357 357 entry->Nclusters = candidate->Nclusters; 358 entry->dNdx = candidate->dNdx; 358 359 359 360 entry->ErrorP = candidate->ErrorP; … … 510 511 entry->DZ = candidate->DZ; 511 512 entry->Nclusters = candidate->Nclusters; 513 entry->dNdx = candidate->dNdx; 512 514 513 515 entry->ErrorP = candidate->ErrorP;
Note:
See TracChangeset
for help on using the changeset viewer.