| 26 |
#include <stdlib.h> |
#include <stdlib.h> |
| 27 |
#include <math.h> |
#include <math.h> |
| 28 |
// |
// |
| 29 |
|
// RunInfo header |
| 30 |
|
// |
| 31 |
|
#include <RunInfo.h> |
| 32 |
|
#include <GLTables.h> |
| 33 |
|
// |
| 34 |
// YODA headers |
// YODA headers |
| 35 |
// |
// |
| 36 |
#include <PamelaRun.h> |
#include <PamelaRun.h> |
| 39 |
#include <EventHeader.h> |
#include <EventHeader.h> |
| 40 |
#include <physics/neutronDetector/NeutronEvent.h> |
#include <physics/neutronDetector/NeutronEvent.h> |
| 41 |
#include <physics/neutronDetector/NeutronRecord.h> |
#include <physics/neutronDetector/NeutronRecord.h> |
|
//#include <yodaUtility.h> |
|
| 42 |
// |
// |
| 43 |
// RunInfo header |
// RunInfo header |
| 44 |
// |
// |
| 45 |
#include <RunInfo.h> |
//#include <RunInfo.h> |
| 46 |
#include <GLTables.h> |
//#include <GLTables.h> |
| 47 |
// |
// |
| 48 |
// This program headers |
// This program headers |
| 49 |
// |
// |
| 58 |
// |
// |
| 59 |
// |
// |
| 60 |
int NDCore(UInt_t run, TFile *file, TSQLServer *dbc, Int_t NDargc, char *NDargv[]){ |
int NDCore(UInt_t run, TFile *file, TSQLServer *dbc, Int_t NDargc, char *NDargv[]){ |
|
// |
|
|
Int_t i = 0; |
|
|
// |
|
|
TString processFolder = "NDFolder"; |
|
| 61 |
// |
// |
| 62 |
// Set these to true to have a very verbose output. |
// Set these to true to have a verbose output. |
| 63 |
// |
// |
| 64 |
Bool_t debug = false; |
Bool_t debug = false; |
|
// |
|
| 65 |
Bool_t verbose = false; |
Bool_t verbose = false; |
| 66 |
// |
// |
| 67 |
|
// Output directory is the working directoy. |
| 68 |
|
// |
| 69 |
|
const char* outDir = gSystem->DirName(gSystem->DirName(file->GetPath())); |
| 70 |
|
// |
| 71 |
|
Int_t i = 0; |
| 72 |
|
TString processFolder = "NDFolder_"+run; |
| 73 |
if ( NDargc > 0 ){ |
if ( NDargc > 0 ){ |
| 74 |
i = 0; |
i = 0; |
| 75 |
while ( i < NDargc ){ |
while ( i < NDargc ){ |
| 90 |
}; |
}; |
| 91 |
}; |
}; |
| 92 |
// |
// |
| 93 |
const char* outDir = gSystem->DirName(gSystem->DirName(file->GetPath())); |
|
| 94 |
// Variables for level2 |
// Variables for level2 |
| 95 |
// |
// |
| 96 |
TTree *NDtr = 0; |
TTree *NDtr = 0; |
| 97 |
UInt_t nevents = 0; |
UInt_t nevents = 0; |
| 98 |
// |
// |
| 99 |
// variables needed to reprocess data |
// Variables needed to reprocess data |
| 100 |
// |
// |
| 101 |
TString NDversion; |
TString NDversion; |
| 102 |
ItoRunInfo *runinfo = 0; |
ItoRunInfo *runinfo = 0; |
| 121 |
Int_t code = 0; |
Int_t code = 0; |
| 122 |
Int_t sgnl; |
Int_t sgnl; |
| 123 |
// |
// |
| 124 |
|
Long64_t maxsize = 10000000000LL; |
| 125 |
|
TTree::SetMaxTreeSize(maxsize); |
| 126 |
|
// |
| 127 |
|
// |
| 128 |
// ND level2 classes |
// ND level2 classes |
| 129 |
// |
// |
| 130 |
NDLevel2 *nd = new NDLevel2(); |
NDLevel2 *nd = new NDLevel2(); |
| 247 |
} |
} |
| 248 |
} else { |
} else { |
| 249 |
// |
// |
| 250 |
|
NDtrclone->SetAutoSave(900000000000000LL); |
| 251 |
|
// |
| 252 |
// tree exists, we are reprocessing data. Are we reprocessing a single run or all the file? |
// tree exists, we are reprocessing data. Are we reprocessing a single run or all the file? |
| 253 |
// |
// |
| 254 |
reproc = true; |
reproc = true; |
| 295 |
// |
// |
| 296 |
file->cd(); |
file->cd(); |
| 297 |
NDtr = new TTree("NeutronD-new","PAMELA Level2 NeutronD data"); |
NDtr = new TTree("NeutronD-new","PAMELA Level2 NeutronD data"); |
| 298 |
|
NDtr->SetAutoSave(900000000000000LL); |
| 299 |
NDtr->Branch("NDLevel2","NDLevel2",&nd); |
NDtr->Branch("NDLevel2","NDLevel2",&nd); |
| 300 |
// |
// |
| 301 |
if ( reproc && !reprocall ){ |
if ( reproc && !reprocall ){ |
| 304 |
// |
// |
| 305 |
tempfile = new TFile(tempname.str().c_str(),"READ"); |
tempfile = new TFile(tempname.str().c_str(),"READ"); |
| 306 |
NDtrclone = (TTree*)tempfile->Get("NeutronD-old"); |
NDtrclone = (TTree*)tempfile->Get("NeutronD-old"); |
| 307 |
|
NDtrclone->SetAutoSave(900000000000000LL); |
| 308 |
NDtrclone->SetBranchAddress("NDLevel2",&ndclone); |
NDtrclone->SetBranchAddress("NDLevel2",&ndclone); |
| 309 |
// |
// |
| 310 |
if ( nobefrun > 0 ){ |
if ( nobefrun > 0 ){ |