remove analysis classes new implemention
[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigKineGammaDirect.C
CommitLineData
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//------------------------------------
11AliAnaGamma* 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}