Changeset fec809d in git for modules/EnergyLoss.cc

Timestamp:

Jan 27, 2020, 12:26:15 PM (5 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

Timing

Children:

32b2ff1

Parents:

6777565

Message:

fixed bug in energy loss, removed printouts, and added comments in the card

File:

• modules/EnergyLoss.cc (modified) (7 diffs)

modules/EnergyLoss.cc

@@ -81 +81 @@
   fX1   =  GetDouble("x1", 2.8716);
   fI    =  GetDouble("I", 173.0); // mean excitation potential in (eV)
-  fX0   =  GetDouble("c0", 4.4355);
+  fC0   =  GetDouble("c0", 4.4355);
 
   // import arrays with output from other modules
@@ -120 +120 @@
 void EnergyLoss::Process()
 {
-  Candidate *candidate, *particle;
+  Candidate *candidate;
   vector<TIterator *>::iterator itInputList;
   TIterator *iterator;
@@ -131 +131 @@
   vector<Double_t> elosses;
 
-  cout<<"---------------- new event -------------------"<<endl;
+  //cout<<"---------------- new event -------------------"<<endl;
 
 
@@ -143 +143 @@
     while((candidate = static_cast<Candidate *>(iterator->Next())))
     {
-      cout<<"    ---------------- new candidate -------------------"<<endl;
+      //cout<<"    ---------------- new candidate -------------------"<<endl;
       const TLorentzVector &candidateMomentum = candidate->Momentum;
 
@@ -149 +149 @@
       gamma     = candidateMomentum.Gamma();
       charge    = TMath::Abs(candidate->Charge);
+
 
       // length of the track normalized by the fraction of active material and the charge collection efficiency in the tracker (in cm)
@@ -161 +162 @@
       nhits = Int_t(L*fActiveFraction/dx);
 
+
+      //beta = 0.999945;
+      //gamma = 95.6446;
+      //charge = 1.;
+      //nhits = 100;      
+
       //cout<<L<<","<<fActiveFraction<<","<<dx<<","<<nhits<<endl;
 
@@ -180 +187 @@
       dP = chi*( TMath::Log(Wmax/I) + TMath::Log(chi/I) + 0.2 - beta*beta - delta);
 
-      if (candidateMomentum.Pt() > 5) {
-        cout<<"Nhits: "<<nhits<<", dx: "<<dx<<", Charge: "<<charge<<", Beta: "<< beta<<", Gamma: "<<gamma<<", PT: "<<candidateMomentum.Pt()<<endl;
+      //cout<<"L: "<<L<<", PT: "<<candidateMomentum.Pt()<<", Eta: "<<candidateMomentum.Eta()<<", Phi: "<< candidateMomentum.Phi()<<endl;
+      //cout<<"Nhits: "<<nhits<<", dx: "<<dx<<", Charge: "<<charge<<", Beta: "<< beta<<", Gamma: "<<gamma<<", PT: "<<candidateMomentum.Pt()<<endl;
       //cout<<x<<","<<kappa<<endl;
-        cout<<"    Wmax: "<<Wmax<<", Chi: "<<chi<<", delta: "<<delta<<", DeDx: "<<avdE<<", DeltaP: "<<dP<<endl;
-      }
+      //cout<<"    Wmax: "<<Wmax<<", Chi: "<<chi<<", delta: "<<delta<<", DeDx: "<<avdE<<", DeltaP: "<<dP<<endl;
 
       // simulate Nhits energy loss measurements