Changeset fec809d in git

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

Files:

• cards/FCC/FCChh_PileUpVtx.tcl (modified) (1 diff)
• modules/EnergyLoss.cc (modified) (7 diffs)
• modules/MomentumSmearing.cc (modified) (1 diff)

cards/FCC/FCChh_PileUpVtx.tcl

@@ -258 +258 @@
 
   # absolute resolution per measurement (normalized in MeV/cm)
-  set Resolution 0.4
+  # CMS pixel detector performance is reproduceable with r = 0.4
+  # dedicated dEdX detector can achieve r = 0.0 or below (i.e better than Landau)
+  
+  #set Resolution 0.4
+  set Resolution 0.2
 
   # fraction of measurements to ignore when computing truncated mean
+  # suggested range [0.4-0.6]
+  
   set TruncatedMeanFraction 0.5
 
   # detector properties (active fraction = nhits*thickness/L)
-  set Thickness 200E-6
-  set ActiveFraction 0.00133
+  set Thickness 100E-6
+  set ActiveFraction 0.0006666
 
   # Silicon properties, for other materials:

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

modules/MomentumSmearing.cc

@@ -120 +120 @@
     candidate->ErrorPT = res*pt;
     candidate->PT = pt;
-    
+        
     candidate->TrackResolution = res;
     candidate->AddCandidate(mother);