doc/examples/example1.C

//
//  Example to integrate all detector info by means of the PamLevel2 class
//
//  An object of this class has access to all the public members of all the detector classes.
//  Some local methods allows to sort tracks and solve the tracker y-view ambiguity (reject track images)
//
void example1(TString file){
    //
    // histos
    //
    TH1F* qtot    =  new TH1F("qtot","Total energy in the calorimeter",1000,0.,2000.);
    TH1F* npaddle =  new TH1F("npaddle","Total number of hit ToF paddles",16,-0.5,15.5);
    TH1F* ntrack  =  new TH1F("ntr","Number of fitted tracks",6,-0.5,5.5);
    //
    TH1F* rig   = new TH1F("rig","Particle rigidity",100,0.,100.);
    TH1F* resxs = new TH1F("resxs","Spatial residual (X) on the 1^ calo plane (sorted tracks)",100,-10.,10.);
    TH1F* resys = new TH1F("resys","Spatial residual (Y) on the 1^ calo plane (sorted tracks)",100,-10.,10.);
    TH1F* resxi = new TH1F("resxi","Spatial residual (X) on the 1^ calo plane (image tracks)",100,-10.,10.);
    TH1F* resyi = new TH1F("resyi","Spatial residual (Y) on the 1^ calo plane (image tracks)",100,-10.,10.); 
    //
    PamLevel2*  pam_event  = new PamLevel2();       // << create pamela event
    //
    TFile f(file);
    //
    TTree *T = pam_event->LoadPamTrees(&f);         // << load Pamela trees from file f
    Int_t nevent = T->GetEntries();
    //
    cout << endl<<" Start loop over events   ";
    for (Int_t i=0; i<nevent;i++){
        //
        T->GetEntry(i);
        //================================================================================
        // some general quantities
        //================================================================================
        // tracker
        ntrack->Fill( pam_event->GetNTracks()); 
        // calorimeter
        qtot->Fill(pam_event->qtot);
        // ToF
        Int_t npa=0;
        for(Int_t ipa=0; ipa<6; ipa++)npa = npa + pam_event->GetNHitPaddles(ipa);
        npaddle->Fill(npa);
        //================================================================================
        // track related variables
        //================================================================================
        for(Int_t it=0; it<pam_event->GetNTracks(); it++){     // << loop over the "physical" tracks (no track images)
            //
            // << get the it-th physical pamela track 
            // << PamTrack combines the tracker, calorimeter and ToF track-related variables
            // << (in case of image, choose the best track by means of calorimeter and ToF data)
            //
            PamTrack *track = pam_event->GetTrack(it);
            //
            // << if the track fit is good, fill some histos
            //
            if( track->chi2 > 0 && track->chi2 < 100){
                // 
                // << spatial residuals on the first calorimeter plane for the sorted track 
                //
                Float_t rxs = track->tbar[0][0] - pam_event->cbar[0][0];
                Float_t rys = track->tbar[0][1] - pam_event->cbar[0][1];
                resxs->Fill(rxs);
                resys->Fill(rys);
                //
                rig->Fill( track->GetRigidity() );
                //
                // << spatial residuals on the first calorimeter plane for the image track 
                //
                if(track->HasImage()){                                // << if the sorted track has an image...
                    //
                    PamTrack *image = pam_event->GetTrackImage(it);   // << ...get the image track
                    //
                    Float_t rxi = image->tbar[0][0] - pam_event->cbar[0][0];
                    Float_t ryi = image->tbar[0][1] - pam_event->cbar[0][1];
                    resxi->Fill(rxi);
                    resyi->Fill(ryi);
                };
            };
        };  // end loop over tracks
        //================================================================================

    };                                       // end loop over the events
    cout << endl << endl << " Done "<< endl<<endl;
    //
    // close file
    //
    f.Close();  
    
};
1	pam-fi	1.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			};