]>
Commit | Line | Data |
---|---|---|
ec908cae | 1 | void AliRsnSelect |
2 | ( | |
3 | Int_t event = 1, | |
4 | Option_t *inRootPath = "/home/pulvir/resonances/aliroot-v4-04-Rev-08/pythia/events", | |
5 | Option_t *outRootPath = "/home/pulvir/resonances/aliroot-v4-04-Rev-08/pythia/selections", | |
6 | Double_t ptLimit4PID = 10000.0 | |
7 | ) | |
8 | { | |
9 | // load ANALYSIS library for output objects | |
10 | gSystem->Load("libPWG2.so"); | |
11 | ||
12 | // instantiate reading manager | |
13 | AliRsnReader *reader = new AliRsnReader; | |
14 | ||
15 | // Define prior probabilities (only for ESD PID) | |
16 | reader->SetPriorProbability(AliPID::kElectron, 0.0339947); //0.0339947 | |
17 | reader->SetPriorProbability(AliPID::kMuon, 0.0192307); //0.0192307 | |
18 | reader->SetPriorProbability(AliPID::kPion, 0.822957); //0.822957 | |
19 | reader->SetPriorProbability(AliPID::kKaon, 0.0751355); //0.0751355 | |
20 | reader->SetPriorProbability(AliPID::kProton, 0.0486821); //0.0486821 | |
21 | reader->SetProbabilityThreshold(0.5); | |
22 | ||
23 | // Define PID method | |
24 | reader->SetPIDMethod(AliRsnReader::kPerfectPID); | |
25 | ||
26 | // create input/output names | |
27 | Text_t inPath[200], outFileName[200]; | |
28 | sprintf(inPath, "%s/%d", inRootPath, event); | |
29 | sprintf(outFileName, "%s/%d.root", outRootPath, event); | |
30 | ||
31 | // do event reading | |
32 | cout << "Reading data in " << inPath << endl; | |
33 | TTree *events = reader->ReadTracksAndParticles(inPath, "R"); | |
34 | if (!events) return; | |
35 | cout << endl; | |
36 | ||
37 | // open output file | |
38 | TFile *fileOut = TFile::Open(outFileName, "RECREATE"); | |
39 | cout << "Saving data in " << outFileName << endl; | |
40 | events->Write(events->GetName(), TObject::kOverwrite); | |
41 | fileOut->Close(); | |
42 | } |