Changeset 84 in svn for trunk/Examples/src/Analysis_Ex.cc

Timestamp:

Dec 4, 2008, 1:11:27 PM (16 years ago)

Author:

uid677

Message:

add output tree to example

File:

• trunk/Examples/src/Analysis_Ex.cc (modified) (4 diffs)

trunk/Examples/src/Analysis_Ex.cc

@@ -44 +44 @@
 }
 
-void Analysis_Ex::Run(ExRootTreeReader *treeReaderGen, ExRootTreeReader *treeReaderRec, ExRootTreeReader *treeReaderTrig)
+void Analysis_Ex::Run(ExRootTreeReader *treeReaderGen, ExRootTreeReader *treeReaderRec, ExRootTreeReader *treeReaderTrig, ExRootTreeWriter *treeWriter)
 { 
   total=0;//initialisation of total number of events
@@ -70 +70 @@
   //to get the trigger information
   const TClonesArray *TRIGGER  = treeReaderTrig->UseBranch("TrigResult");
-  
+
+  //Define the branches that will be filled during the analysis
+  ExRootTreeBranch *INVMASS = treeWriter->NewBranch("INVMass", TRootInvm::Class());
+  TRootInvm *inv_mass; 
   //*******************************************
   
@@ -296 +299 @@
      for(Int_t k = numElec; k < (numElec+mu.GetEntries()); k++)Lept[k].SetPxPyPzE(mu[k-numElec]->Px,mu[k-numElec]->Py,mu[k-numElec]->Pz,mu[k-numElec]->E);
     
+     //Example how to white a branch in the output file
+     inv_mass=(TRootInvm*) INVMASS->NewEntry();
+     inv_mass->M=(Lept[0]+Lept[1]).M();
+ 
      if((Lept[0]+Lept[1]).M() > INV_MASS_LL )continue;// the invariant mass should be < INV_MASS_LL
      cut_2++;//event accepted
+
+     treeWriter->Fill(); 
    }
+ treeWriter->Write();
 
 }
@@ -350 +360 @@
 }
 
-