]>
Commit | Line | Data |
---|---|---|
3a9a3487 | 1 | //-------------------------------------------------------------------------- |
2 | // Test macro for reconstruction and analysis of D0->Kpi | |
3 | // | |
4 | // Andrea Dainese, andrea.dainese@pd.infn.it | |
5 | //-------------------------------------------------------------------------- | |
6 | ||
7 | void AliD0toKpiReco() { | |
8 | ||
9 | //============== R E C O N S T R U C T I O N ============================== | |
10 | ||
11 | // Look for field value in galice.root | |
12 | Double_t field = 0.4; | |
13 | if (gAlice) { | |
14 | delete gAlice->GetRunLoader(); | |
15 | delete gAlice; | |
16 | gAlice=0; | |
17 | } | |
18 | if(!gSystem->AccessPathName("galice.root",kFileExists)) { | |
19 | AliRunLoader *rl = AliRunLoader::Open("galice.root"); | |
20 | rl->LoadgAlice(); | |
c84a5e9e | 21 | field=gAlice->Field(); |
22 | Double_t bz=field->SolenoidField()/10.; | |
23 | printf("B = %3.1f T read from gAlice and set\n",bz); | |
3a9a3487 | 24 | delete gAlice->GetRunLoader(); |
25 | delete gAlice; | |
26 | gAlice=0; | |
27 | } else { | |
28 | printf(" File galice.root not found: default %3.1f T being used!\n",field); | |
29 | } | |
30 | ||
31 | AliD0toKpiAnalysis *analysis = new AliD0toKpiAnalysis(); | |
32 | //--- set magnetic field | |
33 | analysis->SetBz(field); | |
34 | // set simulation to take info on PDG codes from kine tree | |
35 | //analysis->SetSimulation(); | |
36 | //--- set this is you want only signal candidates in output file | |
37 | //analysis->SetOnlySignal(); | |
38 | //--- set this if you want to compute primary vertex D0 by D0 using | |
39 | // other tracks in the event (for pp, broad interaction region); | |
40 | // it is time-consuming procedure, so it can be done after a | |
41 | // preselection on invariant mass | |
42 | //analysis->SetVertexOnTheFly(); | |
43 | //analysis->SetMassCut(0.1); // GeV | |
44 | //--- set single-track preselections | |
45 | analysis->SetPtCut(0.5); // GeV | |
46 | analysis->Setd0Cut(50.); // micron | |
47 | //--- set cuts on D0 candidates to be written to file | |
48 | // (see AliD0toKpiAnalysis.h for a description and for the defaults) | |
49 | analysis->SetD0Cuts(0.1,1000.,1.1,0.,0.,10000.,10000.,0.,.5); | |
50 | //analysis->SetD0Cuts(); | |
51 | ||
52 | //--- check the settings | |
53 | analysis->PrintStatus(); | |
54 | ||
55 | Int_t evFirst = 0; | |
5e653232 | 56 | Int_t evLast = 1000000; |
3a9a3487 | 57 | //analysis->SetDebug(); |
58 | //analysis->FindCandidates(evFirst,evLast,"AliD0toKpi.root"); | |
59 | analysis->FindCandidatesESD(evFirst,evLast,"AliD0toKpi.root"); | |
60 | delete analysis; | |
61 | ||
62 | return; | |
63 | } | |
64 | //========================================================================== | |
65 | void AliD0toKpiSele() { | |
66 | ||
67 | //======================== S E L E C T I O N ============================ | |
68 | ||
69 | AliD0toKpiAnalysis *analysis = new AliD0toKpiAnalysis(); | |
70 | analysis->SetSimulation(); | |
71 | analysis->SetOnlySignal(); | |
72 | analysis->SetD0Cuts(0.1,1000.,1.1,0.,0.,10000.,10000.,0.,.5); | |
73 | analysis->ApplySelection("AliD0toKpi.root","AliD0toKpi_sele.root"); | |
74 | delete analysis; | |
75 | ||
76 | return; | |
77 | } | |
78 | //========================================================================== | |
79 | ||
80 | ||
81 | ||
82 | ||
83 | ||
84 |