| 1 |
// |
| 2 |
// Example to integrate all detector info by means of the PamLevel2 class |
| 3 |
// |
| 4 |
// An object of this class has access to all the public members of all the detector classes. |
| 5 |
// Some local methods allows to sort tracks and solve the tracker y-view ambiguity (reject track images) |
| 6 |
// |
| 7 |
void example1(TString file){ |
| 8 |
// |
| 9 |
// histos |
| 10 |
// |
| 11 |
TH1F* qtot = new TH1F("qtot","Total energy in the calorimeter",1000,0.,2000.); |
| 12 |
TH1F* npaddle = new TH1F("npaddle","Total number of hit ToF paddles",16,-0.5,15.5); |
| 13 |
TH1F* ntrack = new TH1F("ntr","Number of fitted tracks",6,-0.5,5.5); |
| 14 |
// |
| 15 |
TH1F* rig = new TH1F("rig","Particle rigidity",100,0.,100.); |
| 16 |
TH1F* resxs = new TH1F("resxs","Spatial residual (X) on the 1^ calo plane (sorted tracks)",100,-10.,10.); |
| 17 |
TH1F* resys = new TH1F("resys","Spatial residual (Y) on the 1^ calo plane (sorted tracks)",100,-10.,10.); |
| 18 |
TH1F* resxi = new TH1F("resxi","Spatial residual (X) on the 1^ calo plane (image tracks)",100,-10.,10.); |
| 19 |
TH1F* resyi = new TH1F("resyi","Spatial residual (Y) on the 1^ calo plane (image tracks)",100,-10.,10.); |
| 20 |
// |
| 21 |
PamLevel2* pam_event = new PamLevel2(); // << create pamela event |
| 22 |
// |
| 23 |
TFile f(file); |
| 24 |
// |
| 25 |
TTree *T = pam_event->LoadPamTrees(&f); // << load Pamela trees from file f |
| 26 |
Int_t nevent = T->GetEntries(); |
| 27 |
// |
| 28 |
cout << endl<<" Start loop over events "; |
| 29 |
for (Int_t i=0; i<nevent;i++){ |
| 30 |
// |
| 31 |
T->GetEntry(i); |
| 32 |
//================================================================================ |
| 33 |
// some general quantities |
| 34 |
//================================================================================ |
| 35 |
// tracker |
| 36 |
ntrack->Fill( pam_event->GetNTracks()); |
| 37 |
// calorimeter |
| 38 |
qtot->Fill(pam_event->qtot); |
| 39 |
// ToF |
| 40 |
Int_t npa=0; |
| 41 |
for(Int_t ipa=0; ipa<6; ipa++)npa = npa + pam_event->GetNHitPaddles(ipa); |
| 42 |
npaddle->Fill(npa); |
| 43 |
//================================================================================ |
| 44 |
// track related variables |
| 45 |
//================================================================================ |
| 46 |
for(Int_t it=0; it<pam_event->GetNTracks(); it++){ // << loop over the "physical" tracks (no track images) |
| 47 |
// |
| 48 |
// << get the it-th physical pamela track |
| 49 |
// << PamTrack combines the tracker, calorimeter and ToF track-related variables |
| 50 |
// << (in case of image, choose the best track by means of calorimeter and ToF data) |
| 51 |
// |
| 52 |
PamTrack *track = pam_event->GetTrack(it); |
| 53 |
// |
| 54 |
// << if the track fit is good, fill some histos |
| 55 |
// |
| 56 |
if( track->chi2 > 0 && track->chi2 < 100){ |
| 57 |
// |
| 58 |
// << spatial residuals on the first calorimeter plane for the sorted track |
| 59 |
// |
| 60 |
Float_t rxs = track->tbar[0][0] - pam_event->cbar[0][0]; |
| 61 |
Float_t rys = track->tbar[0][1] - pam_event->cbar[0][1]; |
| 62 |
resxs->Fill(rxs); |
| 63 |
resys->Fill(rys); |
| 64 |
// |
| 65 |
rig->Fill( track->GetRigidity() ); |
| 66 |
// |
| 67 |
// << spatial residuals on the first calorimeter plane for the image track |
| 68 |
// |
| 69 |
if(track->HasImage()){ // << if the sorted track has an image... |
| 70 |
// |
| 71 |
PamTrack *image = pam_event->GetTrackImage(it); // << ...get the image track |
| 72 |
// |
| 73 |
Float_t rxi = image->tbar[0][0] - pam_event->cbar[0][0]; |
| 74 |
Float_t ryi = image->tbar[0][1] - pam_event->cbar[0][1]; |
| 75 |
resxi->Fill(rxi); |
| 76 |
resyi->Fill(ryi); |
| 77 |
}; |
| 78 |
}; |
| 79 |
}; // end loop over tracks |
| 80 |
//================================================================================ |
| 81 |
|
| 82 |
}; // end loop over the events |
| 83 |
cout << endl << endl << " Done "<< endl<<endl; |
| 84 |
// |
| 85 |
// close file |
| 86 |
// |
| 87 |
f.Close(); |
| 88 |
|
| 89 |
}; |
Parent Directory
| 
Revision Log