Changeset 481 in svn

Timestamp:

Jul 14, 2009, 12:54:32 PM (15 years ago)

Author:

Xavier Rouby

Message:

new leaves in RESOLJET needed for atlas resolution

Location:

trunk

Files:

• Resolutions_ATLAS.cpp (modified) (3 diffs)
• interface/TreeClasses.h (modified) (1 diff)

trunk/Resolutions_ATLAS.cpp

@@ -237 +237 @@
   TFile *outputFile = TFile::Open(outputfilename.c_str(), "RECREATE");// Creates the file, but should be closed just after
   outputFile->Close();
- 
+
   // 1.3 Reads the trees in  input file 
   TChain chainGEN("GEN");
-  chainGEN.Add(inputfilename.c_str());
+  //chainGEN.Add(inputfilename.c_str());
+  //chainGEN.Add("all_dijets_atlas_kt.root");
+  //chainGEN.Add("all_dijets_atlas_kt_2.root");
+  chainGEN.Add("all_dijets_atlas_jetclu_new.root");
+  chainGEN.Add("all_dijets_atlas_jetclu_new_2.root");
   ExRootTreeReader *treeReaderGEN = new ExRootTreeReader(&chainGEN);
   TChain chain("Analysis");
-  chain.Add(inputfilename.c_str());
+  //chain.Add(inputfilename.c_str());
+  //chain.Add("all_dijets_atlas_kt.root");
+  //chain.Add("all_dijets_atlas_kt_2.root");
+  chain.Add("all_dijets_atlas_jetclu_new.root");
+  chain.Add("all_dijets_atlas_jetclu_new_2.root");
   ExRootTreeReader *treeReader = new ExRootTreeReader(&chain);
   const TClonesArray *branchJet = treeReader->UseBranch("Jet");
@@ -299 +307 @@
   cout <<"*********************************************************************"<< endl;
   cout <<"**                                                                 **"<< endl;
-  cout <<"**                                                                 **"<< endl;
+  cout <<"**        ####### RESOLUTION STUDIES FOR ATLAS   ########          **"<< endl;
   cout <<"**                                                                 **"<< endl;
   cout <<"**        ####### Start resolution processing    ########          **"<< endl;
@@ -417 +425 @@
             elementJet->dE = (JETreco.E()-JETgen.E())/JETgen.E() ;
             elementJet->dE2 = pow( (JETreco.E()-JETgen.E())/JETgen.E() , 2.) ;
+            elementJet->Eta = JETgen.Eta();
+            elementJet->dE_reco = (JETreco.E()-JETgen.E())/JETreco.E() ;
+            elementJet->dE2_reco = pow( (JETreco.E()-JETgen.E())/JETreco.E() , 2.) ;
           }
       }

trunk/interface/TreeClasses.h

@@ -45 +45 @@
 {
 public:
-  float PT;         // generator  PT
-  float SmearedPT;  // recoPT / genPT
+  float SmearedPT;
+  float PT;
   float E;          // generator E
   float dE;         // (recoE-genE)/genE
   float dE2;        // ( (recoE-genE)/ genE )^2
+  float Eta;        // generator-level eta
+  float dE_reco;    // (recoE-genE)/recoE
+  float dE2_reco;   // ( (recoE-genE)/ recoE )^2
+
 
   static TCompare *fgCompare; //!