Changeset 511 in svn for trunk/Resolutions_ATLAS.cpp

Timestamp:

Jul 23, 2009, 3:41:27 PM (15 years ago)

Author:

Xavier Rouby

Message:

print the number of tau-jets

File:

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

trunk/Resolutions_ATLAS.cpp

@@ -82 +82 @@
             if((d1 < array.GetEntries()) && (d1 > 0) && (d2 < array.GetEntries()) && (d2 > 0))
                {
-                 tauhad=true;
+                 //if(fabs(gen1->Eta)<2.5) tauhad=true;
+                 if(gen1->PT<60 && gen1->PT>30)tauhad=true;
                  for(int d=d1; d < d2+1; d++)
                     {
@@ -240 +241 @@
   // 1.3 Reads the trees in  input file 
   TChain chainGEN("GEN");
-  //chainGEN.Add(inputfilename.c_str());
-  chainGEN.Add("all_dijets_atlas_kt.root");
-  chainGEN.Add("all_dijets_atlas_kt_2.root");
+  chainGEN.Add(inputfilename.c_str());
+  //chainGEN.Add("all_dijets_atlas_kt.root");
+  //chainGEN.Add("all_dijets_atlas_kt_2.root");
   ExRootTreeReader *treeReaderGEN = new ExRootTreeReader(&chainGEN);
   TChain chain("Analysis");
-  //chain.Add(inputfilename.c_str());
-  chain.Add("all_dijets_atlas_kt.root");
-  chain.Add("all_dijets_atlas_kt_2.root");
+  chain.Add(inputfilename.c_str());
+  //chain.Add("all_dijets_atlas_kt.root");
+  //chain.Add("all_dijets_atlas_kt_2.root");
   ExRootTreeReader *treeReader = new ExRootTreeReader(&chain);
   const TClonesArray *branchJet = treeReader->UseBranch("Jet");
@@ -445 +446 @@
                     elementTaujet->NumTrack = NumTracks(branchTracks,DET->TAU_track_pt,JETT.Eta(),JETT.Phi(),DET->TAU_track_scone);
                     if( (EnergySmallCone(towers,JETT.Eta(),JETT.Phi(),DET->TAU_energy_scone,DET->JET_seed)/JETT.E()) > 0.95 
-                     &&  (NumTracks(branchTracks,DET->TAU_track_pt,JETT.Eta(),JETT.Phi(),DET->TAU_track_scone))==1)numTauRec++;
+                     &&  (NumTracks(branchTracks,DET->TAU_track_pt,JETT.Eta(),JETT.Phi(),DET->TAU_track_scone))==1) if(TausHadr[i]->PT>30 && TausHadr[i]->PT<60) numTauRec++;
                  }
               }
@@ -458 +459 @@
   globalwatch.Stop();
 
+  cout << "This sample contains " << numTau << " taus that decayed hadronically.\n";
+  cout << numTauRec << " have been reconstructed.\n";
 
   // Screen output