Changeset b4cd1e7 in git

Timestamp:

Oct 12, 2014, 11:12:55 PM (10 years ago)

Author:

Christophe Delaere <christophe.delaere@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

a844c07

Parents:

c4d879d

Message:

Added MET to the display

First version. The MET vector appears properly on the 3D view but not on
the projections... Why?

File:

• examples/geometry.C (modified) (7 diffs)

examples/geometry.C

@@ -74 +74 @@
 TClonesArray *gBranchTrack = 0;
 TClonesArray *gBranchJet = 0;
+TClonesArray *gBranchMet = 0;
 
 DelphesCaloData *gCaloData = 0;
 TEveElementList *gJetList = 0;
+TEveArrow *gMet = 0;
 TEveTrackList *gTrackList = 0;
 
@@ -529 +531 @@
    gBranchTrack = gTreeReader->UseBranch("Track");
    gBranchJet   = gTreeReader->UseBranch("Jet");
+   gBranchMet   = gTreeReader->UseBranch("MissingET");
 
 //TODO make it configurable, for more objects (or can we guess from the config?)
@@ -563 +566 @@
    gJetList = new TEveElementList("Jets");
    gEve->AddElement(gJetList);
+
+   gMet = new TEveArrow(1., 0., 0., 0., 0., 0.);
+   gMet->SetMainColor(kViolet);
+   gMet->SetTubeR(0.02);
+   gMet->SetPickable(kTRUE);
+   gMet->SetName("Missing Et");
+   gEve->GetCurrentEvent()->AddElement(gMet);
 
    gTrackList = new TEveTrackList("Tracks");
@@ -634 +644 @@
   TIter itTrack(gBranchTrack);
   TIter itJet(gBranchJet);
+  TIter itMet(gBranchMet);
 
   Tower *tower;
   Track *track;
   Jet *jet;
+  MissingET *MET;
 
   TEveJetCone *eveJetCone;
@@ -643 +655 @@
 
   Int_t counter;
+  Float_t maxPt = 0.;
 
   TEveTrackPropagator *trkProp = gTrackList->GetPropagator();
@@ -686 +699 @@
     gTrackList->AddElement(eveTrack);
     eveTrack->MakeTrack();
+    maxPt = maxPt > track->PT ? maxPt : track->PT;
   }
   // Loop over all jets
@@ -700 +714 @@
     eveJetCone->AddEllipticCone(jet->Eta, jet->Phi, jet->DeltaEta, jet->DeltaPhi);
     gJetList->AddElement(eveJetCone);
+    maxPt = maxPt > jet->PT ? maxPt : jet->PT;
+  }
+
+  // Missing Et
+  // recipe: gRadius * MET/maxpt(tracks, jets)
+  itMet.Reset();
+  while((MET = (MissingET*) itMet.Next())) {
+    delete gMet;
+    gMet = new TEveArrow((gRadius * MET->MET/maxPt)*cos(MET->Phi), (gRadius * MET->MET/maxPt)*sin(MET->Phi), 0., 0., 0., 0.);
+    gMet->SetMainColor(kViolet);
+//    gMet->SetTubeR(0.04);
+    gMet->SetPickable(kTRUE);
+    gMet->SetName("Missing Et");
+    gMet->SetTitle(Form("Missing Et (%.1f GeV)",MET->MET));
+    gMet->ProjectAllChildren();
+    gEve->GetCurrentEvent()->AddElement(gMet);
   }
 }