Line: 1 to 1 | ||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Continue Electron Charge Misidentification studies 7/14/09: comparison to 31X
| ||||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < | Base: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 22X | |||||||||||||||||||||||||||||||||||
> > | Basic selection | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < |
| |||||||||||||||||||||||||||||||||||
> > | ||||||||||||||||||||||||||||||||||||
Added: | ||||||||||||||||||||||||||||||||||||
> > | Code | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < | Base: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31X | |||||||||||||||||||||||||||||||||||
> > | Show code Hide code <--/twistyPlugin twikiMakeVisibleInline-->for( unsigned int event = 0; event < nEvents; ++event) { cms2.GetEntry(event); ++nEventsTotal; | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < |
| |||||||||||||||||||||||||||||||||||
> > | if ( nEventsTotal%10000 == 0 ) { std::cout << "Event: " << nEventsTotal << endl; } | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < | Yanjun_2: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun (new conversion removal with ≥ 2) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31X | |||||||||||||||||||||||||||||||||||
> > | // loop over true electrons for ( unsigned int els = 0; els < genps_p4().size(); ++els ) { | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < |
| |||||||||||||||||||||||||||||||||||
> > | // check for true electron if ( TMath::Abs(genps_id()[els]) = 11 ) continue; | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < | Yanjun_1: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun (new conversion removal with ≥ 1) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31X | |||||||||||||||||||||||||||||||||||
> > | // check for pT >= 10 GeV if ( genps_p4()[els].pt() < 10. ) continue; | |||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < |
| |||||||||||||||||||||||||||||||||||
> > | // fill true histrograms
els_pt_sim->Fill(genps_p4()[els].pt());
els_eta_sim->Fill(genps_p4()[els].eta());
}
// loop over reco electrons
for (unsigned int els = 0;
els < els_p4().size();
++els) {
// cuts
// check for pT >= 10 GeV
if ( els_p4()[els].pt() < 10. ) continue;
//if ( goodElectronWithoutIsolation(els) ) continue;
if ( goodElectronIsolated(els,true) ) continue;
if ( conversionElectron(els) ) continue;
// check how many electrons don't have an associated track
els_trkId->Fill(els_trkidx().at(els));
// tmp charge variable
double charge = els_charge().at(els);
// if electron has associated track and track charge is not equal to electron charge, veto
int trk = els_trkidx().at(els);
if ( (trk >= 0) && (charge = trks_charge().at(trk)) ) {
continue;
}
// exclude reco electron which is not a true electron
if ( abs(els_mc_id()[els]) = 11 ) continue;
// further reduce misidentifcation rate using the ratio of super cluster energy
// to gsf track momentum
if ( els_eSC()[els]/els_trk_p4()[els].P() < 0.3 ) continue;
// further redyce misidentification rate using the ratio of gsf track momentum
// over ctf track momentum and only veto high tail
if ( trk >= 0 && els_trk_p4()[els].P()/trks_trk_p4()[trk].P() > 3.5 ) continue;
// further redyce misidentification rate using the ratio of gsf track momentum
// over ctf track momentum and also veto lower tail
if ( trk >= 0 && els_trk_p4()[els].P()/trks_trk_p4()[trk].P() < 0.5 ) continue;
// fill reco
els_pt_reco->Fill(els_p4().at(els).Pt());
els_eta_reco->Fill(els_p4().at(els).eta());
// fill reco_corCharge
if ( (charge = -1 && els_mc_id().at(els) = 11) || (charge = 1 && els_mc_id().at(els) = -11) ) {
els_pt_reco_corCharge->Fill(els_p4().at(els).Pt());
els_eta_reco_corCharge->Fill(els_p4().at(els).eta());
}
// fill recosim
els_pt_recosim->Fill(els_mc_p4().at(els).Pt());
els_eta_recosim->Fill(els_mc_p4().at(els).eta());
// correct charge identified
if ( (charge = -1 && els_mc_id().at(els) = 11) || (charge = 1 && els_mc_id().at(els) = -11) ) {
els_pt_recosim_corCharge->Fill(els_mc_p4().at(els).Pt());
els_eta_recosim_corCharge->Fill(els_mc_p4().at(els).eta());
els_ScOvGsfTrkP_corId->Fill(els_eSC()[els]/els_trk_p4()[els].P());
if ( trk >= 0 ) {
els_GfsTrkPOvCtfTrkP_corId->Fill(els_trk_p4()[els].P()/trks_trk_p4()[trk].P());
}
// incorrect charge identified
} else {
els_pt_recosim_incorCharge->Fill(els_mc_p4().at(els).Pt());
els_eta_recosim_incorCharge->Fill(els_mc_p4().at(els).eta());
els_ScOvGsfTrkP_incorId->Fill(els_eSC()[els]/els_trk_p4()[els].P());
if ( trk >= 0 ) {
els_GfsTrkPOvCtfTrkP_incorId->Fill(els_trk_p4()[els].P()/trks_trk_p4()[trk].P());
}
}
}
}
}
<--/twistyPlugin--> Base: GoodIsolatedElectrons veto'ing conversions (Puneeth) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 22XBase: GoodIsolatedElectrons veto'ing conversions (Puneeth) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31X | |||||||||||||||||||||||||||||||||||
GoodIsolatedElectrons veto'ing conversions that don't have a track associated in 22X | ||||||||||||||||||||||||||||||||||||
Line: 37 to 146 | ||||||||||||||||||||||||||||||||||||
Deleted: | ||||||||||||||||||||||||||||||||||||
< < | ||||||||||||||||||||||||||||||||||||
GoodIsolatedElectrons veto'ing conversions that don't have a track associated in 31X | ||||||||||||||||||||||||||||||||||||
Changed: | ||||||||||||||||||||||||||||||||||||
< < | ||||||||||||||||||||||||||||||||||||
> > |
Further reduce charge misidentification rate (starting from base selection)
Base selection with
|
Line: 1 to 1 | ||||||||
---|---|---|---|---|---|---|---|---|
Added: | ||||||||
> > |
Continue Electron Charge Misidentification studies 7/14/09: comparison to 31X
Base: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 22XBase: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31XYanjun_2: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun (new conversion removal with ≥ 2) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31XYanjun_1: GoodIsolatedElectrons veto'ing conversions (Puneeth and Yanjun (new conversion removal with ≥ 1) and veto'ing electrons with inconsistent charge rel. to their assoc. track in 31XGoodIsolatedElectrons veto'ing conversions that don't have a track associated in 22XGoodIsolatedElectrons veto'ing conversions that don't have a track associated in 31X |