doc/examples/example3.C

//
//      Example to loop over events, reading more than one file.
//
//      Input parameters are:
//  - dir:      name of the directory where Level2 files are stored
//  - list:     name of a text file with the list of file names to be processed.
//    If list="" all ROOT files inside the directory are processed (ROOT files that does not
//    match Pamela Level2 structure are automatically discarded)
//  - treelist: string containing the list of trees to be loaded (es: "-ORB-ND" or "-ALL+TRK+CAL+TOF").
//    Possible options are:
//    +ALL +TRK +CAL +TRG +TOF +ND +AC +S4 +ORB
//    -ALL -TRK -CAL -TRG -TOF -ND -AC -S4 -ORB
//      - nmax:     maximum number of events to be processed
//
//  All selected trees in all the files are chained and made friends.
//  The whole "blob" can be accessed as if it was a single tree.
//
//  --- Modified on January 2007 ---
//
//      see the comment at the beginning of example1.C      
//
void example3(TString dir,TString list="", TString treelist="+ALL", Int_t nmax=500000000){

    //   Create a timer object to benchmark this loop
    TStopwatch timer;
    timer.Start();
    TString wd=gSystem->WorkingDirectory();
    //
    PamLevel2*  event  = new PamLevel2();       // << create pamela event
        
    // =========================================================
    // << First a TList of root files is created
    // << (either from a specific list or from the whole input directory)
    TList *l  = event->GetListOfLevel2Files(dir,list);
    // << Hence a TChain is created from the list. 
    // << All the trees required by the user are made friend.
    TChain *T = event->GetPamTree(l,treelist);
    // << In the following the TChain can be used similarly to a TTree.

    //
    TChain *RUNT = event->GetRunTree(dir,list); // << get the run information and store it in RUNT

    //
    // NB! It is still possible to read a single file and get a TTree, by means of the methods:
    // TTree* LoadPamTrees(TFile *f);
    // TTree* LoadPamTrees(TFile *f,TString treelist);
    // =========================================================
        
    Int_t nevent = T->GetEntries();     
                
    Int_t ntrk=0;
    if(nevent < nmax)nmax=nevent;
    cout << endl<<" Start loop over events   "<<endl;
    for (Int_t i=0; i<nmax;i++){
        //
        event->Clear();
        T->GetEntry(i);

        if ( event->UpdateRunInfo(RUNT,i) ) printf(" New RUN at entry = %llu RUN ID is %u NEVENTS is %u \n",i,event->GetRunInfo()->ID,event->GetRunInfo()->NEVENTS);
                
        if(event->GetTrkLevel2()->GetNTracks() > 1)cout << " Ev "<<i << " multiple tracks: "<< event->GetTrkLevel2()->GetNTracks() << endl;
                
        if(event->GetTrkLevel2()->ntrk() >0){
            ntrk++;
        };
                
        for(Int_t it=0; it<event->GetTrkLevel2()->GetNTracks(); it++){
                        
            PamTrack *track = event->GetTrack(it);
            PamTrack *image = 0;
                        
            if(track->GetTrkTrack()->HasImage())image = event->GetTrackImage(it);
                        
        };

    };                                       // end loop over the events
    cout << endl << endl << " Done "<< endl<<endl;
    cout << ntrk <<" tracks over "<<nmax<<" events ("<< 100*ntrk/nmax<<"%)"<<endl;

//  Stop timer and print results
    timer.Stop();
    Double_t rtime = timer.RealTime();
    Double_t ctime = timer.CpuTime();
    printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);

    gSystem->ChangeDirectory(wd);
        
};
1	pam-fi	1.1	//
2			// Example to loop over events, reading more than one file.
3			//
4			// Input parameters are:
5			// - dir: name of the directory where Level2 files are stored
6			// - list: name of a text file with the list of file names to be processed.
7			// If list="" all ROOT files inside the directory are processed (ROOT files that does not
8			// match Pamela Level2 structure are automatically discarded)
9			// - treelist: string containing the list of trees to be loaded (es: "-ORB-ND" or "-ALL+TRK+CAL+TOF").
10			// Possible options are:
11			// +ALL +TRK +CAL +TRG +TOF +ND +AC +S4 +ORB
12			// -ALL -TRK -CAL -TRG -TOF -ND -AC -S4 -ORB
13			// - nmax: maximum number of events to be processed
14			//
15			// All selected trees in all the files are chained and made friends.
16			// The whole "blob" can be accessed as if it was a single tree.
17			//
18	pam-fi	1.3	// --- Modified on January 2007 ---
19			//
20			// see the comment at the beginning of example1.C
21			//
22	pam-fi	1.1	void example3(TString dir,TString list="", TString treelist="+ALL", Int_t nmax=500000000){
23
24	pam-fi	1.3	// Create a timer object to benchmark this loop
25			TStopwatch timer;
26			timer.Start();
27			TString wd=gSystem->WorkingDirectory();
28			//
29			PamLevel2* event = new PamLevel2(); // << create pamela event
30	pam-fi	1.1
31	pam-fi	1.3	// =========================================================
32			// << First a TList of root files is created
33			// << (either from a specific list or from the whole input directory)
34			TList *l = event->GetListOfLevel2Files(dir,list);
35			// << Hence a TChain is created from the list.
36			// << All the trees required by the user are made friend.
37			TChain *T = event->GetPamTree(l,treelist);
38			// << In the following the TChain can be used similarly to a TTree.
39	mocchiut	1.4
40			//
41			TChain *RUNT = event->GetRunTree(dir,list); // << get the run information and store it in RUNT
42
43	pam-fi	1.3	//
44			// NB! It is still possible to read a single file and get a TTree, by means of the methods:
45			// TTree* LoadPamTrees(TFile *f);
46			// TTree* LoadPamTrees(TFile *f,TString treelist);
47			// =========================================================
48	pam-fi	1.1
49	pam-fi	1.3	Int_t nevent = T->GetEntries();
50	pam-fi	1.1
51	pam-fi	1.3	Int_t ntrk=0;
52			if(nevent < nmax)nmax=nevent;
53			cout << endl<<" Start loop over events "<<endl;
54	pam-fi	1.1	for (Int_t i=0; i<nmax;i++){
55			//
56	pam-fi	1.3	event->Clear();
57			T->GetEntry(i);
58	mocchiut	1.4
59			if ( event->UpdateRunInfo(RUNT,i) ) printf(" New RUN at entry = %llu RUN ID is %u NEVENTS is %u \n",i,event->GetRunInfo()->ID,event->GetRunInfo()->NEVENTS);
60	pam-fi	1.1
61	pam-fi	1.3	if(event->GetTrkLevel2()->GetNTracks() > 1)cout << " Ev "<<i << " multiple tracks: "<< event->GetTrkLevel2()->GetNTracks() << endl;
62	pam-fi	1.1
63	pam-fi	1.3	if(event->GetTrkLevel2()->ntrk() >0){
64			ntrk++;
65			};
66	pam-fi	1.1
67	pam-fi	1.3	for(Int_t it=0; it<event->GetTrkLevel2()->GetNTracks(); it++){
68	pam-fi	1.1
69	pam-fi	1.3	PamTrack *track = event->GetTrack(it);
70			PamTrack *image = 0;
71	pam-fi	1.1
72	pam-fi	1.3	if(track->GetTrkTrack()->HasImage())image = event->GetTrackImage(it);
73	pam-fi	1.1
74	pam-fi	1.3	};
75	pam-fi	1.1
76			}; // end loop over the events
77			cout << endl << endl << " Done "<< endl<<endl;
78	pam-fi	1.3	cout << ntrk <<" tracks over "<<nmax<<" events ("<< 100*ntrk/nmax<<"%)"<<endl;
79	pam-fi	1.1
80			// Stop timer and print results
81	pam-fi	1.3	timer.Stop();
82			Double_t rtime = timer.RealTime();
83			Double_t ctime = timer.CpuTime();
84			printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
85	pam-fi	1.1
86	pam-fi	1.3	gSystem->ChangeDirectory(wd);
87	pam-fi	1.1
88			};