Changeset cd75093 in git for modules/IdentificationMap.cc

Timestamp:

Dec 12, 2014, 2:58:04 PM (10 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

0c871e5

Parents:

17826f2 (diff), 3d10d1f (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge pull request #8 from selvaggi/master

Added external fast jet examples.

File:

• modules/IdentificationMap.cc (modified) (1 diff)

modules/IdentificationMap.cc

@@ -181 +181 @@
       Pi = formula->Eval(pt, eta);
     
-     // cout<<"PdgIn: "<<pdgIn<<", PdgOut: "<<pdgOut<<", Pi = "<<Pi<<endl;
-    
       // check that sum of probabilities does not exceed 1.
-      
       P = (P - Pi)/P;
         
-      if( P < 0.0 ) 
-      {
-       // cout<<"Sum of probabilities exceeds 1 for particle ID "<<pdgIn<<", check formulas in IdentificationMap module"<<endl;
-        continue;
-      }
+      if( P < 0.0 ) continue;
       else
       {
       
-        //randomly assign a PID to particle according to map
-      
-        Double_t rndm = gRandom->Uniform();
-       // cout<<"Random number vs P //  "<<rndm<<"  "<<P<<endl;
-     
-        if(rndm > P)
-        //if(gRandom->Uniform() > P)
-        {
-         // cout<<"particle "<<candidate->PID<<" --> ";
-        
-          //change PID of particle
-          if(pdgOut != 0) candidate->PID = charge*pdgOut;
-         // cout<<"Final PID: "<<candidate->PID<<endl;
-        
-          fOutputArray->Add(candidate);
-          break;
-      }
-     }  
+       //randomly assign a PID to particle according to map
+       Double_t rndm = gRandom->Uniform();
+       
+       if(rndm > P)
+       {
+         //change PID of particle
+         if(pdgOut != 0) candidate->PID = charge*pdgOut;
+         fOutputArray->Add(candidate);
+         break;
+       }
+      } 
         
     }