Changeset 290 in svn for trunk/Delphes.cpp

Timestamp:

Mar 3, 2009, 12:48:15 AM (15 years ago)

Author:

Xavier Rouby

Message:

tracks and Taujets have charge

File:

• trunk/Delphes.cpp (modified) (5 diffs)

trunk/Delphes.cpp

@@ -318 +318 @@
   vector<fastjet::PseudoJet> input_particles;//for FastJet algorithm
   vector<fastjet::PseudoJet> sorted_jets;
-  vector<TLorentzVector> TrackCentral;  
+  //vector<TLorentzVector> TrackCentral;  
+  vector<TRootTracks> TrackCentral;  
   vector<PhysicsTower> towers;
   //vector<ParticleUtil> electron;
@@ -326 +327 @@
   //vector<ParticleUtil> gamma;
   vector<D_Particle> gamma;
+
+bool FLAG_lhco = true;
 
   TSimpleArray<TRootGenParticle> NFCentralQ;
@@ -489 +492 @@
                    if( (rand()%100) < DET->TRACK_eff && sET!=0) {
                      elementTrack = (TRootTracks*) branchTrack->NewEntry();
-                     elementTrack->Set(particle->Eta, particle->Phi, particle->EtaCalo, particle->PhiCalo, sET);
-                     TrackCentral.push_back(elementTrack->GetFourVector()); // tracks at vertex!
-                     // TODO!!! associates the tracks to the calo where it points to
+                     elementTrack->Set(particle->Eta, particle->Phi, particle->EtaCalo, particle->PhiCalo, sET, particle->Charge);
+                     //TrackCentral.push_back(elementTrack->GetFourVector()); // tracks at vertex!
+                     TrackCentral.push_back(*elementTrack); // tracks at vertex!
                      // TODO!!! apply a smearing on the position of the origin of the track
                      // TODO!!! elementTracks->SetPositionOut(Xout,Yout,Zout);
@@ -577 +580 @@
         elementElec->Charge = sign(electron[i].PID());
         elementElec->IsolFlag = DET->Isolation(electron[i].Phi(),electron[i].Eta(),TrackCentral,2.0);//isolation based on tracks
-      }
+      }                                                                                 /////////////// HARDCODING
       DET->SortedVector(muon);
       for(unsigned int i=0; i < muon.size(); i++) {
@@ -585 +588 @@
         elementMu->EtaCalo = muon[i].EtaCalo(); 
         elementMu->PhiCalo = muon[i].PhiCalo();
-        elementMu->IsolFlag = DET->Isolation(muon[i].Phi(),muon[i].Eta(),TrackCentral,2.0);
+        elementMu->IsolFlag = DET->Isolation(muon[i].Phi(),muon[i].Eta(),TrackCentral,2.0); /////////////// HARDCODING
       }
       DET->SortedVector(gamma);