Changes in modules/RunPUPPI.cc [d6ce231:659c7b6] in git

File:

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

modules/RunPUPPI.cc

@@ -166 +166 @@
   Candidate *pv = static_cast<Candidate*>(fPVItInputArray->Next());
   if (pv) PVZ = pv->Position.Z();
+
   // Fill input particles for puppi
   std::vector<RecoObj> puppiInputVector;
@@ -178 +179 @@
       curRecoObj.phi = momentum.Phi();
       curRecoObj.m   = momentum.M();  
-      particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));//if(fApplyNoLep && TMath::Abs(candidate->PID) == 11) continue; //Dumb cut to minimize the nolepton on electron
+      particle = static_cast<Candidate*>(candidate->GetCandidates()->Last());
+      //if(fApplyNoLep && TMath::Abs(candidate->PID) == 11) continue; //Dumb cut to minimize the nolepton on electron
       //if(fApplyNoLep && TMath::Abs(candidate->PID) == 13) continue;
       if (candidate->IsRecoPU and candidate->Charge !=0) { // if it comes fromPU vertexes after the resolution smearing and the dZ matching within resolution
@@ -190 +192 @@
         curRecoObj.dZ = particle->Position.Z()-PVZ;
       } 
-      else if(!candidate->IsRecoPU && candidate->Charge !=0) {
+      else if(!candidate->IsRecoPU and candidate->Charge !=0) {
         curRecoObj.id    = 1;  // charge from LV
         curRecoObj.vtxId = 1; // from PV
@@ -217 +219 @@
       curRecoObj.m   = momentum.M();
       curRecoObj.charge = 0;
-      particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));
+      particle = static_cast<Candidate*>(candidate->GetCandidates()->Last());
+
       if(candidate->Charge == 0){
         curRecoObj.id    = 0; // neutrals have id==0