628982481fd44db94c496a4f0066eb1559ffe1bc
[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigKineGammaDirect.C
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 }