Difference: OliverGutscheElsChargeMisidentification ( vs. 1)

Revision 127 Apr 2009 - Main.OliverGutsche

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="OliverGutsche"

Electron Charge Misidentification

Input, 10k single particle gun electrons (2 electrons per event, antiPartilce = true)

uaf-4:/data/tmp/yanjuntu/singleElectron_ntuple_v3.root

All electrons in els_p4

      for (unsigned int els = 0; 
           els < els_p4().size(); els++) {

     els_pt->Fill(els_mc_p4().at(els).Pt());
     els_eta->Fill(els_mc_p4().at(els).eta());

     // correct charge identified
     if ( (els_charge().at(els) == -1 && els_mc_id().at(els) == 11) || (els_charge().at(els) == 1 && els_mc_id().at(els) == -11) ) {
       els_pt_corCharge->Fill(els_mc_p4().at(els).Pt());
       els_eta_corCharge->Fill(els_mc_p4().at(els).eta());
       // incorrect charge identified
     } else {
       els_pt_incorCharge->Fill(els_mc_p4().at(els).Pt());
       els_eta_incorCharge->Fill(els_mc_p4().at(els).eta());
     }
      }

eta and pt


electron misidentification percentage in eta and pt

  els_pt_incorCharge->Divide(els_pt);
  els_eta_incorCharge->Divide(els_eta);

  • overall: 429/16946=2.5%


good electrons without isolation

      for (unsigned int els = 0; 
           els < els_p4().size(); els++) {

   if ( goodElectronWithoutIsolation(els) ) {

     els_pt->Fill(els_mc_p4().at(els).Pt());
     els_eta->Fill(els_mc_p4().at(els).eta());

     // correct charge identified
     if ( (els_charge().at(els) == -1 && els_mc_id().at(els) == 11) || (els_charge().at(els) == 1 && els_mc_id().at(els) == -11) ) {
       els_pt_corCharge->Fill(els_mc_p4().at(els).Pt());
       els_eta_corCharge->Fill(els_mc_p4().at(els).eta());
       // incorrect charge identified
     } else {
       els_pt_incorCharge->Fill(els_mc_p4().at(els).Pt());
       els_eta_incorCharge->Fill(els_mc_p4().at(els).eta());
     }
   }
      }

eta and pt


electron misidentification percentage in eta and pt

  els_pt_incorCharge->Divide(els_pt);
  els_eta_incorCharge->Divide(els_eta);

  • overall: 327/15972=2.0%


 
This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback