4 //------------------------------------
5 // Configuration macro example:
7 // Do example analysis with MC
9 // Author : Gustavo Conesa Balbastre (INFN-LNF)
10 //------------------------------------
12 AliAnaPartCorrMaker* ConfigAnalysis()
15 // Configuration goes here
17 printf("======================== \n");
18 printf("ConfigAnalysis() \n");
19 printf("======================== \n");
22 //Detector Fiducial Cuts
23 AliFiducialCut * fidCut = new AliFiducialCut();
24 fidCut->DoCTSFiducialCut(kFALSE) ;
25 fidCut->DoEMCALFiducialCut(kTRUE) ;
26 fidCut->DoPHOSFiducialCut(kTRUE) ;
28 //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
29 fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
30 fidCut->SetSimplePHOSFiducialCut(0.12,220.,320.);
32 // //Fiducial cut EMCAL, 5 regions
33 // Float_t etamax[]={0.67,0.51,0.16,-0.21,-0.61};
34 // TArrayF etamaxarr(5,etamax);
35 // fidCut->AddEMCALFidCutMaxEtaArray(etamaxarr);
36 // Float_t etamin[]={0.61,0.21,-0.16,-0.51,-0.67};
37 // TArrayF etaminarr(5,etamin);
38 // fidCut->AddEMCALFidCutMinEtaArray(etaminarr);
39 // Float_t phimax[]={99.*TMath::DegToRad(), 119.*TMath::DegToRad(), 139.*TMath::DegToRad(),
40 // 159.*TMath::DegToRad(), 179.*TMath::DegToRad(), 189.*TMath::DegToRad()};
41 // TArrayF phimaxarr(6,phimax);
42 // fidCut->AddEMCALFidCutMaxPhiArray(phimaxarr);
43 // Float_t phimin[]={81.*TMath::DegToRad(), 101.*TMath::DegToRad(), 121.*TMath::DegToRad(),
44 // 141.*TMath::DegToRad(), 161.*TMath::DegToRad(), 181.*TMath::DegToRad()};
45 // TArrayF phiminarr(6,phimin);
46 // fidCut->AddEMCALFidCutMinPhiArray(phiminarr);
50 //-----------------------------------------------------------
52 //-----------------------------------------------------------
53 AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();
56 //Switch on or off the detectors information that you want
57 reader->SwitchOffEMCAL();
58 reader->SwitchOffCTS();
59 reader->SwitchOnPHOS();
60 reader->SwitchOffEMCALCells();
61 reader->SwitchOnPHOSCells();
64 reader->SetEMCALPtMin(0.);
65 reader->SetPHOSPtMin(0.);
66 reader->SetCTSPtMin(.2);
68 //Remove the temporal AODs we create.
69 reader->SwitchOnCleanStdAOD();
71 reader->SetFiducialCut(fidCut);
73 // //We want tracks fitted in the detectors:
74 // ULong_t status=AliAODTrack::kTPCrefit;
75 // status|=AliAODTrack::kITSrefit; //(default settings)
77 // We want tracks whose PID bit is set:
78 // ULong_t status =AliAODTrack::kITSpid;
79 // status|=AliAODTrack::kTPCpid;
81 // reader->SetTrackStatus(status);
86 //---------------------------------------------------------------------
88 //---------------------------------------------------------------------
90 //Detector Fiducial Cuts for analysis part
91 AliFiducialCut * fidCut2 = new AliFiducialCut();
92 fidCut2->DoCTSFiducialCut(kFALSE) ;
93 fidCut2->DoEMCALFiducialCut(kTRUE) ;
94 fidCut2->DoPHOSFiducialCut(kTRUE) ;
96 //fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
97 fidCut2->SetSimpleEMCALFiducialCut(0.5,100.,150.);
98 fidCut2->SetSimplePHOSFiducialCut(0.1,240.,280.);
100 AliCaloPID * pid = new AliCaloPID();
101 // use selection with simple weights
102 pid->SetPHOSPhotonWeight(0.7); pid->SetPHOSPi0Weight(0.7);
103 pid->SetEMCALPhotonWeight(0.7); pid->SetEMCALPi0Weight(0.7);
104 // use more complicated selection, particle weight depending on cluster energy
105 // pid->UsePHOSPIDWeightFormula(kTRUE);
106 // TFormula * photonF = new TFormula("photonWeight","0.98*(x<40)+ 0.68*(x>=100)+(x>=40 && x<100)*(0.98+x*(6e-3)-x*x*(2e-04)+x*x*x*(1.1e-06))");
107 // TFormula * pi0F = new TFormula("pi0Weight","0.98*(x<65)+ 0.915*(x>=100)+(x>=65 && x-x*(1.95e-3)-x*x*(4.31e-05)+x*x*x*(3.61e-07))");
108 // pid->SetPHOSPhotonWeightFormula(photonF);
109 // pid->SetPHOSPi0WeightFormula(pi0F);
112 AliAnaExample *ana = new AliAnaExample();
114 ana->SetCaloPID(pid);
115 ana->SetFiducialCut(fidCut2);
116 ana->SetDetector("PHOS");
117 //ana->SwitchOnDataMC();
119 ana->SetOutputAODName("Example");
120 ana->SetOutputAODClassName("AliAODPWG4Particle"); //Or AliAODPWG4ParticleCorrelation
121 //Set Histrograms bins and ranges
122 // ana->SetHistoPtRangeAndNBins(0, 50, 100) ;
123 // ana->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
124 // ana->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
126 //AliAnaExample *ana2 = new AliAnaExample();
128 //---------------------------------------------------------------------
129 // Set analysis algorithm and reader
130 //---------------------------------------------------------------------
131 maker = new AliAnaPartCorrMaker();
132 maker->SetReader(reader);//pointer to reader
133 maker->AddAnalysis(ana,0);
134 //maker->AddAnalysis(ana2,1);
135 maker->SetAnaDebug(-1) ;
136 maker->SwitchOnHistogramsMaker() ;
137 //maker->SwitchOffHistogramsMaker() ;
138 maker->SwitchOnAODsMaker() ;
139 //maker->SwitchOffAODsMaker() ;
143 printf("======================== \n");
144 printf("END ConfigAnalysis() \n");
145 printf("======================== \n");