Fork me on GitHub

Changeset 511 in svn


Ignore:
Timestamp:
Jul 23, 2009, 3:41:27 PM (15 years ago)
Author:
Xavier Rouby
Message:

print the number of tau-jets

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/Resolutions.cpp

    r506 r511  
    8282            if((d1 < array.GetEntries()) && (d1 > 0) && (d2 < array.GetEntries()) && (d2 > 0))
    8383               {
     84                 //cout << "eta du tau = " << gen1->Eta << endl;
     85                 //if(fabs(gen1->Eta)<2.5)tauhad=true;
    8486                 tauhad=true;
    8587                 for(int d=d1; d < d2+1; d++)
     
    445447  globalwatch.Stop();
    446448
     449  cout << "This sample contains " << numTau << " taus that decayed hadronically.\n";
     450  cout << numTauRec << " have been reconstructed.\n";
    447451
    448452    // Screen output
  • trunk/Resolutions_ATLAS.cpp

    r485 r511  
    8282            if((d1 < array.GetEntries()) && (d1 > 0) && (d2 < array.GetEntries()) && (d2 > 0))
    8383               {
    84                  tauhad=true;
     84                 //if(fabs(gen1->Eta)<2.5) tauhad=true;
     85                 if(gen1->PT<60 && gen1->PT>30)tauhad=true;
    8586                 for(int d=d1; d < d2+1; d++)
    8687                    {
     
    240241  // 1.3 Reads the trees in  input file
    241242  TChain chainGEN("GEN");
    242   //chainGEN.Add(inputfilename.c_str());
    243   chainGEN.Add("all_dijets_atlas_kt.root");
    244   chainGEN.Add("all_dijets_atlas_kt_2.root");
     243  chainGEN.Add(inputfilename.c_str());
     244  //chainGEN.Add("all_dijets_atlas_kt.root");
     245  //chainGEN.Add("all_dijets_atlas_kt_2.root");
    245246  ExRootTreeReader *treeReaderGEN = new ExRootTreeReader(&chainGEN);
    246247  TChain chain("Analysis");
    247   //chain.Add(inputfilename.c_str());
    248   chain.Add("all_dijets_atlas_kt.root");
    249   chain.Add("all_dijets_atlas_kt_2.root");
     248  chain.Add(inputfilename.c_str());
     249  //chain.Add("all_dijets_atlas_kt.root");
     250  //chain.Add("all_dijets_atlas_kt_2.root");
    250251  ExRootTreeReader *treeReader = new ExRootTreeReader(&chain);
    251252  const TClonesArray *branchJet = treeReader->UseBranch("Jet");
     
    445446                    elementTaujet->NumTrack = NumTracks(branchTracks,DET->TAU_track_pt,JETT.Eta(),JETT.Phi(),DET->TAU_track_scone);
    446447                    if( (EnergySmallCone(towers,JETT.Eta(),JETT.Phi(),DET->TAU_energy_scone,DET->JET_seed)/JETT.E()) > 0.95
    447                      &&  (NumTracks(branchTracks,DET->TAU_track_pt,JETT.Eta(),JETT.Phi(),DET->TAU_track_scone))==1)numTauRec++;
     448                     &&  (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++;
    448449                 }
    449450              }
     
    458459  globalwatch.Stop();
    459460
     461  cout << "This sample contains " << numTau << " taus that decayed hadronically.\n";
     462  cout << numTauRec << " have been reconstructed.\n";
    460463
    461464  // Screen output
Note: See TracChangeset for help on using the changeset viewer.