| 25 |
|
// |
| 26 |
|
|
| 27 |
|
|
| 28 |
< |
pGlobalChi2_=new TH1F( "globalChi2","Global chi2", 60, 0, 20 ); |
| 28 |
> |
pGlobalChi2_=new TH1F( "globalChi2","Global chi2", 60, 0, 6 ); |
| 29 |
|
pGlobalChi2_->SetDirectory(pDirectory); |
| 30 |
|
|
| 31 |
+ |
pNumberOfValidPixel_=new TH1F( "numberOfValidPixel", "Number of valid pixel", 11, -0.5, 10.5 ); |
| 32 |
+ |
pNumberOfValidPixel_->SetDirectory(pDirectory); |
| 33 |
+ |
|
| 34 |
|
pNumberOfPixelHits_=new TH1F( "numberOfPixelHits", "Number of pixel hits", 11, -0.5, 10.5 ); |
| 35 |
|
pNumberOfPixelHits_->SetDirectory(pDirectory); |
| 36 |
|
|
| 52 |
|
pEta_=new TH1F( "eta", "Eta", 60, -3, 3 ); |
| 53 |
|
pEta_->SetDirectory(pDirectory); |
| 54 |
|
|
| 55 |
+ |
pPhi_=new TH1F( "phi", "phi", 60, -3.15, 3.15 ); |
| 56 |
+ |
pPhi_->SetDirectory(pDirectory); |
| 57 |
+ |
|
| 58 |
|
pPt_=new TH1F( "pT", "pT", 60, 0, 250 ); |
| 59 |
|
pPt_->SetDirectory(pDirectory); |
| 60 |
|
|
| 91 |
|
pRho25CorrectedIsolation_=new TH1F( "rho25CorrectedIsolation", "rho25CorrectedIsolation", 120, 0, 2 ); |
| 92 |
|
pRho25CorrectedIsolation_->SetDirectory(pDirectory); |
| 93 |
|
|
| 94 |
+ |
pThirdRho25CorrectedIsolation_=new TH1F( "thirdRho25CorrectedIsolation", "thirdRho25CorrectedIsolation", 120, 0, 2 ); |
| 95 |
+ |
pThirdRho25CorrectedIsolation_->SetDirectory(pDirectory); |
| 96 |
+ |
|
| 97 |
+ |
pDeltaBetaCorrectedIsolationNoZeroing_=new TH1F( "deltaBetaCorrectedIsolationNoZeroing", "deltaBetaCorrectedIsolationNoZeroing", 120, -1, 1 ); |
| 98 |
+ |
pDeltaBetaCorrectedIsolationNoZeroing_->SetDirectory(pDirectory); |
| 99 |
+ |
|
| 100 |
+ |
pRho25CorrectedIsolationNoZeroing_=new TH1F( "rho25CorrectedIsolationNoZeroing", "rho25CorrectedIsolationNoZeroing", 120, -1, 1 ); |
| 101 |
+ |
pRho25CorrectedIsolationNoZeroing_->SetDirectory(pDirectory); |
| 102 |
+ |
|
| 103 |
+ |
pThirdRho25CorrectedIsolationNoZeroing_=new TH1F( "thirdRho25CorrectedIsolationNoZeroing", "thirdRho25CorrectedIsolationNoZeroing", 120, -1, 1 ); |
| 104 |
+ |
pThirdRho25CorrectedIsolationNoZeroing_->SetDirectory(pDirectory); |
| 105 |
+ |
|
| 106 |
|
if( createIsolationTree_ ) |
| 107 |
|
{ |
| 108 |
|
pIsolationTree_=new TTree("isolationTree","Data about the muon isolation"); |
| 139 |
|
pNumberOfMatches_->Fill( muon.nMatches ); |
| 140 |
|
pIPDB_->Fill( muon.ipDb ); |
| 141 |
|
pEta_->Fill( muon.p4.Eta() ); |
| 142 |
+ |
pPhi_->Fill( muon.p4.Phi() ); |
| 143 |
|
pPt_->Fill( muon.p4.Pt() ); |
| 144 |
|
pChargedIsolation_->Fill( muon.pfChaIso ); |
| 145 |
|
pPhotonIsolation_->Fill( muon.pfPhoIso ); |
| 151 |
|
pRelativePileupIsolation_->Fill( muon.pfChaPUIso/muon.p4.Pt() ); |
| 152 |
|
pRelativeIsolation_->Fill( combinedRelativeIsolation(muon) ); |
| 153 |
|
pDeltaBetaCorrectedIsolation_->Fill( deltaBetaCorrectedIsolation(muon) ); |
| 154 |
< |
if( pAuxInfo ) pRho25CorrectedIsolation_->Fill( rho25CorrectedIsolation(muon,rho25) ); // Don't bother filling if I couldn't get the rho25 info |
| 154 |
> |
pDeltaBetaCorrectedIsolationNoZeroing_->Fill( deltaBetaCorrectedIsolationNoZeroing(muon) ); |
| 155 |
> |
if( pAuxInfo ) |
| 156 |
> |
{ |
| 157 |
> |
// Don't bother filling if I couldn't get the rho25 info |
| 158 |
> |
pRho25CorrectedIsolation_->Fill( rho25CorrectedIsolation(muon,rho25) ); |
| 159 |
> |
pRho25CorrectedIsolationNoZeroing_->Fill( rho25CorrectedIsolationNoZeroing(muon,rho25) ); |
| 160 |
> |
pThirdRho25CorrectedIsolation_->Fill( thirdRho25CorrectedIsolation(muon,rho25) ); |
| 161 |
> |
pThirdRho25CorrectedIsolationNoZeroing_->Fill( thirdRho25CorrectedIsolationNoZeroing(muon,rho25) ); |
| 162 |
> |
} |
| 163 |
|
|
| 164 |
|
// If the TTree is not null then isolation data has been requested |
| 165 |
|
if( pIsolationTree_ ) |
| 193 |
|
float trkupgradeanalysis::MuonInfoPlotSet::rho25CorrectedIsolation( const VHbbEvent::MuonInfo& muon, float rho25 ) |
| 194 |
|
{ |
| 195 |
|
float uncorrectedIsolation=combinedRelativeIsolation(muon); |
| 196 |
< |
float correctedIsolation=uncorrectedIsolation-(M_PI*0.4*0.4*rho25)/muon.p4.Pt(); |
| 196 |
> |
float correctedIsolation=uncorrectedIsolation-(0.4*0.4*M_PI*rho25)/muon.p4.Pt(); |
| 197 |
|
if( correctedIsolation<0 ) correctedIsolation=0; |
| 198 |
|
|
| 199 |
|
return correctedIsolation; |
| 200 |
|
} |
| 201 |
+ |
|
| 202 |
+ |
float trkupgradeanalysis::MuonInfoPlotSet::thirdRho25CorrectedIsolation( const VHbbEvent::MuonInfo& muon, float rho25 ) |
| 203 |
+ |
{ |
| 204 |
+ |
float uncorrectedIsolation=combinedRelativeIsolation(muon); |
| 205 |
+ |
float correctedIsolation=uncorrectedIsolation-(0.4*0.4*M_PI*rho25)/3.0/muon.p4.Pt(); |
| 206 |
+ |
if( correctedIsolation<0 ) correctedIsolation=0; |
| 207 |
+ |
|
| 208 |
+ |
return correctedIsolation; |
| 209 |
+ |
} |
| 210 |
+ |
|
| 211 |
+ |
float trkupgradeanalysis::MuonInfoPlotSet::deltaBetaCorrectedIsolationNoZeroing( const VHbbEvent::MuonInfo& muon, float deltaBetaFactor ) |
| 212 |
+ |
{ |
| 213 |
+ |
// deltaBetaFactor defaults to -0.5 |
| 214 |
+ |
float correctedNeutralIsolation=muon.pfPhoIso+muon.pfNeuIso+deltaBetaFactor*muon.pfChaPUIso; |
| 215 |
+ |
|
| 216 |
+ |
return (muon.pfChaIso+correctedNeutralIsolation)/muon.p4.Pt(); |
| 217 |
+ |
} |
| 218 |
+ |
|
| 219 |
+ |
float trkupgradeanalysis::MuonInfoPlotSet::rho25CorrectedIsolationNoZeroing( const VHbbEvent::MuonInfo& muon, float rho25 ) |
| 220 |
+ |
{ |
| 221 |
+ |
float uncorrectedIsolation=combinedRelativeIsolation(muon); |
| 222 |
+ |
float correctedIsolation=uncorrectedIsolation-(0.4*0.4*M_PI*rho25)/muon.p4.Pt(); |
| 223 |
+ |
|
| 224 |
+ |
|
| 225 |
+ |
return correctedIsolation; |
| 226 |
+ |
} |
| 227 |
+ |
|
| 228 |
+ |
float trkupgradeanalysis::MuonInfoPlotSet::thirdRho25CorrectedIsolationNoZeroing( const VHbbEvent::MuonInfo& muon, float rho25 ) |
| 229 |
+ |
{ |
| 230 |
+ |
float uncorrectedIsolation=combinedRelativeIsolation(muon); |
| 231 |
+ |
float correctedIsolation=uncorrectedIsolation-(0.4*0.4*M_PI*rho25)/3.0/muon.p4.Pt(); |
| 232 |
+ |
|
| 233 |
+ |
|
| 234 |
+ |
return correctedIsolation; |
| 235 |
+ |
} |
