- Timestamp:
- Mar 1, 2021, 4:21:41 PM (4 years ago)
- Branches:
- master
- Children:
- a6e8029
- Parents:
- 825beb7
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
examples/Example6.C
r825beb7 rc9e7466 140 140 Double_t pullCtg = (obsCtg - genCtg) / trk->ErrorCtgTheta; 141 141 142 if (fabs(pullZ0) > 5)143 {144 // Get associated generated particle145 Int_t unique = gp->GetUniqueID();146 //cout<<unique<<endl;147 148 if (branchGenPart->GetEntries() > 0)149 {150 151 // Loop on tracks152 for (Int_t itp = 0; itp < branchGenPart->GetEntries(); itp++)153 {154 GenParticle* part = (GenParticle*)branchGenPart->At(itp);155 if (part->GetUniqueID() == unique)156 {157 cout<<" ------------ found: ------------------------"<<unique<<endl;158 159 cout<<gp->X<<","<<part->X<<","<<TMath::Power((gp->X - part->X),2)*1e6<<endl;160 cout<<gp->Y<<","<<part->Y<<","<<TMath::Power((gp->Y - part->Y),2)*1e6<<endl;161 cout<<gp->Z<<","<<part->Z<<","<<TMath::Power((gp->Z - part->Z),2)*1e6<<endl;162 cout<<gp->Px<<","<<part->Px<<","<<TMath::Power((gp->Px - part->Px),2)*1e6<<endl;163 cout<<gp->Py<<","<<part->Py<<","<<TMath::Power((gp->Py - part->Py),2)*1e6<<endl;164 cout<<gp->Pz<<","<<part->Pz<<","<<TMath::Power((gp->Pz - part->Pz),2)*1e6<<endl;165 166 }167 }168 }169 170 171 }172 173 174 142 // 175 143 histDpull->Fill(pullD0);
Note:
See TracChangeset
for help on using the changeset viewer.