Fork me on GitHub

Changeset c9e7466 in git for examples/Example6.C


Ignore:
Timestamp:
Mar 1, 2021, 4:21:41 PM (4 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
master
Children:
a6e8029
Parents:
825beb7
Message:

fixed trailing printouts

File:
1 edited

Legend:

Unmodified
Added
Removed
  • examples/Example6.C

    r825beb7 rc9e7466  
    140140                                Double_t pullCtg = (obsCtg - genCtg) / trk->ErrorCtgTheta;
    141141
    142         if (fabs(pullZ0) > 5)
    143                                 {
    144                                         // Get associated generated particle
    145                                         Int_t unique = gp->GetUniqueID();
    146                                         //cout<<unique<<endl;
    147 
    148                                         if (branchGenPart->GetEntries() > 0)
    149                                         {
    150 
    151                                         // Loop on tracks
    152                                                 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 
    174142                                //
    175143                                histDpull->Fill(pullD0);
Note: See TracChangeset for help on using the changeset viewer.