Changeset 71648c2 in git

Timestamp:

May 16, 2013, 4:25:05 PM (11 years ago)

Author:

pavel <pavel@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

8608ef8

Parents:

809e3a9

Message:

add Rho branch

Files:

• classes/ClassesLinkDef.h (modified) (1 diff)
• classes/DelphesClasses.h (modified) (1 diff)
• examples/delphes_card_CMS_PileUp.tcl (modified) (1 diff)
• modules/TreeWriter.cc (modified) (2 diffs)
• modules/TreeWriter.h (modified) (1 diff)

classes/ClassesLinkDef.h

@@ -36 +36 @@
 #pragma link C++ class MissingET+;
 #pragma link C++ class ScalarHT+;
+#pragma link C++ class Rho+;
 #pragma link C++ class Photon+;
 #pragma link C++ class Electron+;

classes/DelphesClasses.h

@@ -165 +165 @@
 //---------------------------------------------------------------------------
 
+class Rho: public TObject
+{
+public:
+  Float_t Rho; // rho energy density
+
+  ClassDef(Rho, 1)
+};
+
+//---------------------------------------------------------------------------
+
 class Photon: public SortableObject
 {

examples/delphes_card_CMS_PileUp.tcl

@@ -626 +626 @@
   add Branch MissingET/momentum MissingET MissingET
   add Branch ScalarHT/energy ScalarHT ScalarHT
-}
-
+  add Branch Rho/rho Rho Rho
+}
+

modules/TreeWriter.cc

@@ -64 +64 @@
   fClassMap[MissingET::Class()] = &TreeWriter::ProcessMissingET;
   fClassMap[ScalarHT::Class()] = &TreeWriter::ProcessScalarHT;
+  fClassMap[Rho::Class()] = &TreeWriter::ProcessRho;
 
   TBranchMap::iterator itBranchMap;
@@ -519 +520 @@
 //------------------------------------------------------------------------------
 
+void TreeWriter::ProcessRho(ExRootTreeBranch *branch, TObjArray *array)
+{
+  Candidate *candidate = 0;
+  Rho *entry = 0;
+
+  // get the first entry
+  if((candidate = static_cast<Candidate*>(array->At(0))))
+  {
+    const TLorentzVector &momentum = candidate->Momentum;
+
+    entry = static_cast<Rho*>(branch->NewEntry());
+
+    entry->Rho = momentum.Pt();
+  }
+}
+
+//------------------------------------------------------------------------------
+
 void TreeWriter::Process()
 {

modules/TreeWriter.h

@@ -50 +50 @@
   void ProcessMissingET(ExRootTreeBranch *branch, TObjArray *array);
   void ProcessScalarHT(ExRootTreeBranch *branch, TObjArray *array);
+  void ProcessRho(ExRootTreeBranch *branch, TObjArray *array);
 
 #ifndef __CINT__