Changeset 481 in svn for trunk/Resolutions_ATLAS.cpp

Timestamp:

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

Author:

Xavier Rouby

Message:

new leaves in RESOLJET needed for atlas resolution

File:

• trunk/Resolutions_ATLAS.cpp (modified) (3 diffs)

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.) ;
           }
       }