1 #define AliAnalysisTaskMuonAODfromGeneral_cxx
6 // Class implementation for the specific muon AOD generation
7 // Extracts only muon tracks from a general AOD and builds dimuons
8 // Livio Bianchi, Universita' di Torino
18 #include "AliAODEvent.h"
19 #include "AliAnalysisTask.h"
20 #include "AliAnalysisTaskMuonAODfromGeneral.h"
21 #include "AliAODHandler.h"
22 #include "AliAnalysisManager.h"
23 #include "AliAODEventInfo.h"
24 #include "AliAODDimuon.h"
26 ClassImp(AliAnalysisTaskMuonAODfromGeneral)
28 //________________________________________________________________________
29 AliAnalysisTaskMuonAODfromGeneral::AliAnalysisTaskMuonAODfromGeneral(const char *name, Double_t BeamEnergy):
30 AliAnalysisTask(name, "AnalysisTaskMuonAODfromGeneral"),
39 // Input slot #0 works with an Ntuple
40 SetBeamEnergy(BeamEnergy);
41 DefineInput(0, TChain::Class());
42 // Output slot #0 writes into a TTree container
43 DefineOutput(0, TTree::Class());
46 //________________________________________________________________________
47 void AliAnalysisTaskMuonAODfromGeneral::ConnectInputData(Option_t *) {
48 printf(" ConnectInputData %s\n", GetName());
49 // New connection to the AOD
50 fChain = (TChain*)GetInputData(0);
51 fOrgAOD = new AliAODEvent();
52 fOrgAOD->ReadFromTree(fChain); // This checks also if the branch address is already set
55 //________________________________________________________________________
56 void AliAnalysisTaskMuonAODfromGeneral::CreateOutputObjects() {
57 printf("Creating output objects\n");
58 if(!fNewAOD)fNewAOD=new AliAODEvent();
59 fNewAOD->CreateStdContent();
62 ft = new TTree("AOD","Muon AOD tree",0);
64 ft->Branch(fNewAOD->GetList());
68 //Add new Objects: MuonInfos & Dimuons
69 fInfos = new AliAODEventInfo();
70 fNewAOD->AddObject(fInfos);
71 const char *nameInfos = "MuonInfos";
72 fInfos->SetName(nameInfos);
73 ft->Branch(nameInfos, &fInfos);
74 fDimuons = new TClonesArray("AliAODDimuon",0);
75 fNewAOD->AddObject(fDimuons);
76 const char *nameDimuons = "Dimuons";
77 fDimuons->SetName(nameDimuons);
78 ft->Branch(nameDimuons, &fDimuons);
82 //________________________________________________________________________
83 void AliAnalysisTaskMuonAODfromGeneral::Exec(Option_t *) {
85 // Task making a Muon AOD
87 TChain *chain = (TChain*)GetInputData(0);
88 Long64_t ientry = chain->GetReadEntry();
93 Int_t nTracks=fOrgAOD->GetNumberOfTracks();
96 Int_t mutrNumb[10]; for (Int_t i=0; i<10; mutrNumb[i++]=0) {}
97 for (Int_t iTrack=0; iTrack<nTracks; iTrack++){
98 const Double_t *trackpid=fOrgAOD->GetTrack(iTrack)->PID();
99 if (trackpid[AliAODTrack::kMuon]==1.) {
100 mutrNumb[nMuTracks]=iTrack;
103 if (fOrgAOD->GetTrack(iTrack)->Charge()> 0) nPosTracks++;
106 Bool_t ExistMuon=kFALSE;
107 if (nMuTracks>0) ExistMuon=kTRUE;
109 //--------------------------------------------------------------------
110 // set arrays and pointers
114 // Access to the header
115 AliAODHeader *header = fNewAOD->GetHeader();
118 header->SetRunNumber (fOrgAOD->GetRunNumber() );
119 header->SetBunchCrossNumber(fOrgAOD->GetBunchCrossNumber());
120 header->SetOrbitNumber (fOrgAOD->GetOrbitNumber() );
121 header->SetPeriodNumber (fOrgAOD->GetPeriodNumber() );
122 header->SetTriggerMask (fOrgAOD->GetTriggerMask() );
123 header->SetTriggerCluster (fOrgAOD->GetTriggerCluster() );
124 header->SetEventType (fOrgAOD->GetEventType() );
125 header->SetMagneticField (fOrgAOD->GetMagneticField() );
126 header->SetZDCN1Energy (fOrgAOD->GetZDCN1Energy() );
127 header->SetZDCP1Energy (fOrgAOD->GetZDCP1Energy() );
128 header->SetZDCN2Energy (fOrgAOD->GetZDCN2Energy() );
129 header->SetZDCP2Energy (fOrgAOD->GetZDCP2Energy() );
130 header->SetZDCEMEnergy (0,fOrgAOD->GetZDCEMEnergy(0) );
131 header->SetZDCEMEnergy (1,fOrgAOD->GetZDCEMEnergy(1) );
132 header->SetRefMultiplicity (nTracks);
133 header->SetRefMultiplicityPos(nPosTracks);
134 header->SetRefMultiplicityNeg(nTracks - nPosTracks);
135 header->SetMuonMagFieldScale(-999.); // FIXME
136 header->SetCentrality(0); // FIXME
138 const AliAODVertex *vtx = fOrgAOD->GetPrimaryVertex();
141 AliError("Primary vertex is not defined");
145 Int_t nDimuons=nMuTracks*(nMuTracks-1)/2;
147 // Access to the AOD container of vertices
148 TClonesArray &vertices = *(fNewAOD->GetVertices());
151 // Access to the AOD container of tracks
152 TClonesArray &tracks = *(fNewAOD->GetTracks());
157 vtx->GetCovMatrix(covVtx); //covariance matrix
159 AliAODVertex * primary = new(vertices[jVertices++])
160 AliAODVertex(pos, covVtx, vtx->GetChi2perNDF(), NULL, AliAODVertex::kPrimary);
163 static int numtracks=0;
164 for (Int_t iTrack=0; iTrack<nMuTracks; iTrack++) {
165 primary->AddDaughter(new(tracks[jTracks++])
166 AliAODTrack((*(fOrgAOD->GetTrack(mutrNumb[iTrack])))));
172 fDimuons->Expand(nDimuons);
173 TClonesArray &rDimuons = *fDimuons;
176 for(Int_t i=0; i<nMuTracks; i++){
177 for(Int_t j=i+1; j<nMuTracks; j++){
178 new(rDimuons[jDimuons++]) AliAODDimuon(tracks[i],tracks[j]);
182 fInfos->SetBeamEnergy(fBeamEnergy);
183 fInfos->SetEv(fNewAOD);
184 fInfos->SetHe(header);
185 fInfos->SetTr(fNewAOD->GetTracks());
186 fInfos->SetDi(fDimuons);
187 fInfos->SelectTriggerBits(0,1,2,3,4,5);
188 if(ExistMuon) ft->Fill();
191 // Post final data. It will be written to a file with option "RECREATE"
195 //________________________________________________________________________
196 void AliAnalysisTaskMuonAODfromGeneral::Terminate(Option_t *) {
197 // Draw some histogram at the end.
198 ft->GetUserInfo()->Add(fNewAOD);
199 ft->Write("",TObject::kOverwrite);
200 ft->ResetBranchAddresses();
201 if (!gROOT->IsBatch()) {
202 TCanvas *c1 = new TCanvas("c1","Pt",10,10,310,310);
203 c1->SetFillColor(10);
204 c1->SetHighLightColor(10);
206 c1->cd(1)->SetLeftMargin(0.15);
207 c1->cd(1)->SetBottomMargin(0.15);
208 c1->cd(1)->SetLogy();
209 ft->Draw("tracks.Pt()");
210 //ft->Draw("dimuons.M()");