1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
19 // Macro for checking trigger integrated efficiency for dimuons.
20 // The efficiency is calculated with respect to the 3/4 coincidence.
21 // Author: Fabien Guerin, LPC Clermont-Ferrand, Jan. 2006
24 #if !defined(__CINT__) || defined(__MAKECINT__)
27 #include "TClonesArray.h"
33 #include "TParticle.h"
37 #include "AliRunLoader.h"
38 #include "AliHeader.h"
39 #include "AliLoader.h"
44 #include "AliMUONData.h"
45 #include "AliMUONHit.h"
46 #include "AliMUONConstants.h"
47 #include "AliMUONDigit.h"
48 #include "AliMUONRawCluster.h"
49 #include "AliMUONGlobalTrigger.h"
50 #include "AliMUONLocalTrigger.h"
51 #include "AliMUONTrack.h"
57 void MUONTriggerEfficiency (char filename[10]="galice.root"){
61 char currentfile[100];
62 sprintf(digitdat,"MUONTriggerEfficiency.out");
64 FILE *fdat=fopen(digitdat,"w");
66 Int_t nevents,coincmuon,muonapt,muonlpt,muonhpt;
67 Int_t CoincMuPlus,CoincMuMinus;
75 // Creating Run Loader and openning file containing Hits
76 AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","READ");
78 if (RunLoader ==0x0) {
79 printf(">>> Error : Error Opening %s file \n",currentfile);
83 nevents = RunLoader->GetNumberOfEvents();
86 AliLoader * MUONHits = RunLoader->GetLoader("MUONLoader");
87 MUONHits->LoadHits("READ");
88 AliMUONData data_hits(MUONHits,"MUON","MUON");
91 AliLoader * MUONDigits = RunLoader->GetLoader("MUONLoader");
92 MUONDigits->LoadDigits("READ");
93 AliMUONData data_digits(MUONDigits,"MUON","MUON");
95 TClonesArray * globalTrigger;
96 AliMUONGlobalTrigger * gloTrg;
99 for (Int_t ievent=0; ievent<nevents; ievent++) { // event loop
104 RunLoader->GetEvent(ievent);
106 if (ievent%1000==0) printf("\t Event = %d\n",ievent);
109 data_hits.SetTreeAddress("H");
111 Int_t itrack, ntracks, NbHits[2][4];
113 for (Int_t j=0; j<2; j++) {
114 for (Int_t jj=0; jj<4; jj++) {
118 ntracks = (Int_t) data_hits.GetNtracks();
119 for (itrack=0; itrack<ntracks; itrack++) { // Track loop
120 data_hits.GetTrack(itrack);
123 nhits = (Int_t) data_hits.Hits()->GetEntriesFast();
125 for(ihit=0; ihit<nhits; ihit++) {
126 mHit = static_cast<AliMUONHit*>(data_hits.Hits()->At(ihit));
127 Int_t Nch = mHit->Chamber();
128 Int_t hittrack = mHit->Track();
129 Float_t IdPart = mHit->Particle();
131 for (Int_t j=0;j<4;j++) {
133 if (hittrack==1 || hittrack==2) {
134 if (Nch==kch && IdPart==-13 && NbHits[0][j]==0) NbHits[0][j]++;
135 if (Nch==kch && IdPart==13 && NbHits[1][j]==0) NbHits[1][j]++;
140 data_hits.ResetHits();
145 SumNbHits=NbHits[0][0]+NbHits[0][1]+NbHits[0][2]+NbHits[0][3];
146 if (SumNbHits==3 || SumNbHits==4) CoincMuPlus=1;
148 SumNbHits=NbHits[1][0]+NbHits[1][1]+NbHits[1][2]+NbHits[1][3];
149 if (SumNbHits==3 || SumNbHits==4) CoincMuMinus=1;
151 if (CoincMuPlus==1 && CoincMuMinus==1) coincmuon++;
154 data_digits.SetTreeAddress("D,GLT");
155 data_digits.GetTriggerD();
157 globalTrigger = data_digits.GlobalTrigger();
159 Int_t nglobals = (Int_t) globalTrigger->GetEntriesFast(); // should be 1
161 for (Int_t iglobal=0; iglobal<nglobals; iglobal++) { // Global Trigger
162 gloTrg = static_cast<AliMUONGlobalTrigger*>(globalTrigger->At(iglobal));
164 if (gloTrg->PairUnlikeApt()>=1) muonapt++;
165 if (gloTrg->PairUnlikeLpt()>=1) muonlpt++;
166 if (gloTrg->PairUnlikeHpt()>=1) muonhpt++;
168 } // end of loop on Global Trigger
170 data_digits.ResetTrigger();
172 } // end loop on event
174 MUONHits->UnloadHits();
175 MUONDigits->UnloadDigits();
178 // calculate efficiency with as a ref. at least 3/4 planes fired
179 Float_t efficiencyapt,efficiencylpt,efficiencyhpt;
180 Double_t coincmu,aptmu,lptmu,hptmu;
182 coincmu=Double_t(coincmuon);
183 aptmu=Double_t(muonapt);
184 lptmu=Double_t(muonlpt);
185 hptmu=Double_t(muonhpt);
190 fprintf(fdat," Number of events = %d \n",nevents);
191 fprintf(fdat," Number of events with 3/4 coinc = %d \n",coincmuon);
192 fprintf(fdat," Number of dimuons with 3/4 coinc Apt cut = %d \n",muonapt);
193 fprintf(fdat," Nomber of dimuons with 3/4 coinc Lpt cut = %d \n",muonlpt);
194 fprintf(fdat," Number of dimuons with 3/4 coinc Hpt cut = %d \n",muonhpt);
197 efficiencyapt=aptmu/coincmu;
198 error=efficiencyapt*TMath::Sqrt((aptmu+coincmu)/(aptmu*coincmu));
199 fprintf(fdat," Efficiency Apt cut = %4.4f +/- %4.4f\n",efficiencyapt,error);
201 efficiencylpt=lptmu/coincmu;
202 error=efficiencylpt*TMath::Sqrt((lptmu+coincmu)/(lptmu*coincmu));
203 fprintf(fdat," Efficiency Lpt cut = %4.4f +/- %4.4f\n",efficiencylpt,error);
205 efficiencyhpt=hptmu/coincmu;
206 error=efficiencyhpt*TMath::Sqrt((hptmu+coincmu)/(hptmu*coincmu));
207 fprintf(fdat," Efficiency Hpt cut = %4.4f +/- %4.4f\n",efficiencyhpt,error);