]>
Commit | Line | Data |
---|---|---|
34ce8d31 | 1 | /* $Id$ */ |
2 | /* $Log$ */ | |
3 | ||
4 | //------------------------------------ | |
5 | // Configuration macro example: | |
6 | // | |
7 | // Do prompt photon analysis with kinematics | |
8 | // | |
9 | // Author : Gustavo Conesa Balbastre (INFN-LNF) | |
10 | //------------------------------------ | |
11 | AliAnaGamma* ConfigGammaAnalysis() | |
12 | { | |
13 | // | |
14 | // Configuration goes here | |
15 | // | |
16 | printf("ConfigGammaAnalysis() \n"); | |
17 | ||
18 | //----------------------------------------------------------- | |
19 | // Reader | |
20 | //----------------------------------------------------------- | |
21 | ||
22 | AliGammaMCReader *reader = new AliGammaMCReader(); | |
23 | // Switch "on" or "off" detectors used in study | |
24 | reader->SwitchOnEMCAL(kTRUE); | |
25 | reader->SwitchOnPHOS(kFALSE) ; | |
26 | reader->SwitchOnCTS(kTRUE) ; | |
27 | ||
28 | //Set detectors acceptance (in real and montecarlo data) | |
29 | reader->SetCTSEtaCut(1); reader->SetEMCALEtaCut(0.7); //reader->SetPHOSEtaCut(0.2); | |
30 | reader->SetPhiEMCALCut(40*TMath::DegToRad(), 200*TMath::DegToRad()); | |
31 | // reader->SetPhiPHOSCut(200*TMath::DegToRad(), 350*TMath::DegToRad()); | |
32 | ||
33 | //Set minimum pt for particles in analysis (in real and montecarlo data) | |
34 | reader->SetNeutralPtCut(0.5); reader->SetChargedPtCut(0.3); | |
35 | ||
36 | reader->SetDecayPi0Flag(AliGammaMCReader::kGeantDecay) ; //Options | |
37 | //kNoDecay :Do not decay pi0, keep them in the list as they are | |
38 | //kGeantDecay: Look for gamma decayed by GEANT | |
39 | //kDecay: Decay pi0 by hand (geant was not used) | |
40 | //kDecayGamma: Pi0 is decayed by PYTHIA, pi0 is not final check if photons overlapp | |
41 | ||
42 | //parameters to study if decay is overlapped: | |
43 | reader->SetEMCALIPDistance(450.); reader->SetPHOSIPDistance(460.); | |
44 | reader->SetEMCALMinAngle(2.5 * TMath::DegToRad() ); | |
45 | reader->SetPHOSMinAngle(0.45 * TMath::DegToRad() ); //Minimum overlapp distance | |
46 | reader->SetCheckOverlapping(kTRUE); | |
47 | ||
48 | //============================ | |
49 | //Prompt photon algoritm | |
50 | //============================== | |
51 | AliAnaGammaDirect *gd = new AliAnaGammaDirect(); | |
52 | gd->SetMinGammaPt(5.); | |
53 | gd->SetConeSize(0.5); gd->SetPtThreshold(1); gd->SetPtSumThreshold(.); | |
54 | gd->SetICMethod(AliAnaGammaDirect::kPtIC) ;//Options: | |
55 | //kNoIC: Accept all photons, no isolation used | |
56 | //kPtIC: IC with cut on pT | |
57 | //kSumPtIC: IC with cut on pT sum in cone | |
58 | ||
59 | //--------------------------------------------------------------------- | |
60 | // Finally: Set analysis algorithm and reader | |
61 | //--------------------------------------------------------------------- | |
62 | ana = new AliAnaGamma(); | |
63 | ana->SetReader(reader);//pointer to reader | |
64 | ana->SetAnalysisType(AliAnaGamma::kPrompt); //set kPrompt, kCorrelation | |
65 | ana->SetGammaDirect(gd);//pointer to direct photon algorithm | |
66 | ana->SetCalorimeter("EMCAL"); //Prompt photon calorimeter | |
67 | ||
68 | // | |
69 | return ana ; | |
70 | } |