[u/mrichter/AliRoot.git] / CORRFW / AliCFAcceptanceCuts.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

///////////////////////////////////////////////////////////////////////////
//          ----   CORRECTION FRAMEWORK   ----
// AliCFAcceptanceCuts implementation
// Class to cut on the number of AliTrackReference's 
// for each detector
///////////////////////////////////////////////////////////////////////////
// author : R. Vernet (renaud.vernet@cern.ch)
///////////////////////////////////////////////////////////////////////////

#include "AliLog.h"
#include "AliMCParticle.h"
#include "AliCFAcceptanceCuts.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TList.h"
#include "TBits.h"

ClassImp(AliCFAcceptanceCuts)

//______________________________
AliCFAcceptanceCuts::AliCFAcceptanceCuts() : 
  AliCFCutBase(),
  fMCInfo(0x0),
  fMinNHitITS(0),
  fMinNHitTPC(0),
  fMinNHitTRD(0),
  fMinNHitTOF(0),
  fMinNHitMUON(0),
  fhCutStatistics(0x0),
  fhCutCorrelation(0x0),
  fBitmap(new TBits(0))
{
  //
  //ctor
  //
  for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=0x0;
}

//______________________________
AliCFAcceptanceCuts::AliCFAcceptanceCuts(const Char_t* name, const Char_t* title) : 
  AliCFCutBase(name,title),
  fMCInfo(0x0),
  fMinNHitITS(0),
  fMinNHitTPC(0),
  fMinNHitTRD(0),
  fMinNHitTOF(0),
  fMinNHitMUON(0),
  fhCutStatistics(0x0),
  fhCutCorrelation(0x0),
  fBitmap(new TBits(0))
{
  //
  //ctor
  //
  for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=0x0;
}

//______________________________
AliCFAcceptanceCuts::AliCFAcceptanceCuts(const AliCFAcceptanceCuts& c) : 
  AliCFCutBase(c),
  fMCInfo(c.fMCInfo),
  fMinNHitITS(c.fMinNHitITS),
  fMinNHitTPC(c.fMinNHitTPC),
  fMinNHitTRD(c.fMinNHitTRD),
  fMinNHitTOF(c.fMinNHitTOF),
  fMinNHitMUON(c.fMinNHitMUON),
  fhCutStatistics(c.fhCutStatistics),
  fhCutCorrelation(c.fhCutCorrelation),
  fBitmap(c.fBitmap)
{
  //
  //copy ctor
  //
  for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=c.fhQA[i][j];
}

//______________________________
AliCFAcceptanceCuts& AliCFAcceptanceCuts::operator=(const AliCFAcceptanceCuts& c)
{
  //
  // Assignment operator
  //
  if (this != &c) {
    AliCFCutBase::operator=(c) ;
    fMCInfo=c.fMCInfo;
    fMinNHitITS=c.fMinNHitITS;
    fMinNHitTPC=c.fMinNHitTPC;
    fMinNHitTRD=c.fMinNHitTRD;
    fMinNHitTOF=c.fMinNHitTOF;
    fMinNHitMUON=c.fMinNHitMUON;
    fhCutStatistics  = c.fhCutStatistics ;
    fhCutCorrelation = c.fhCutCorrelation ;
    fBitmap          = c.fBitmap ;
    for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=c.fhQA[i][j];
  }
  return *this ;
}

//______________________________________________________________
Bool_t AliCFAcceptanceCuts::IsSelected(TObject *obj) {
  //
  // check if selections on 'obj' are passed
  // 'obj' must be an AliMCParticle
  //
  
  SelectionBitMap(obj);

  if (fIsQAOn) FillHistograms(obj,kFALSE);

  for (UInt_t icut=0; icut<fBitmap->GetNbits(); icut++)
    if (!fBitmap->TestBitNumber(icut)) return kFALSE ;
  
  if (fIsQAOn) FillHistograms(obj,kTRUE);
  return kTRUE;
}

//______________________________
void AliCFAcceptanceCuts::SelectionBitMap(TObject* obj) {
  //
  // checks the number of track references associated to 'obj'
  // 'obj' must be an AliMCParticle
  //

  for (UInt_t i=0; i<kNCuts; i++) fBitmap->SetBitNumber(i,kFALSE);

  if (!obj) return;
  TString className(obj->ClassName());
  if (className.CompareTo("AliMCParticle") != 0) {
    AliError("obj must point to an AliMCParticle !");
    return ;
  }

  AliMCParticle * part = dynamic_cast<AliMCParticle*>(obj);
  if(!part) return ;

  Int_t nHitsITS=0, nHitsTPC=0, nHitsTRD=0, nHitsTOF=0, nHitsMUON=0 ;
  for (Int_t iTrackRef=0; iTrackRef<part->GetNumberOfTrackReferences(); iTrackRef++) {
    AliTrackReference * trackRef = part->GetTrackReference(iTrackRef);
    if(trackRef){
      Int_t detectorId = trackRef->DetectorId();
      switch(detectorId) {
      case AliTrackReference::kITS  : nHitsITS++  ; break ;
      case AliTrackReference::kTPC  : nHitsTPC++  ; break ;
      case AliTrackReference::kTRD  : nHitsTRD++  ; break ;
      case AliTrackReference::kTOF  : nHitsTOF++  ; break ;
      case AliTrackReference::kMUON : nHitsMUON++ ; break ;
      default : break ;
      }
    }
  }
  
  Int_t iCutBit = 0;

  if (nHitsITS  >= fMinNHitITS  ) fBitmap->SetBitNumber(iCutBit,kTRUE);
  iCutBit++;
  if (nHitsTPC  >= fMinNHitTPC  ) fBitmap->SetBitNumber(iCutBit,kTRUE);
  iCutBit++;
  if (nHitsTRD  >= fMinNHitTRD  ) fBitmap->SetBitNumber(iCutBit,kTRUE);
  iCutBit++;
  if (nHitsTOF  >= fMinNHitTOF  ) fBitmap->SetBitNumber(iCutBit,kTRUE);
  iCutBit++;
  if (nHitsMUON  >= fMinNHitMUON) fBitmap->SetBitNumber(iCutBit,kTRUE);
}


void AliCFAcceptanceCuts::SetMCEventInfo(const TObject* mcInfo) {
  //
  // Sets pointer to MC event information (AliMCEvent)
  //

  if (!mcInfo) {
    AliError("Pointer to AliMCEvent !");
    return;
  }
  
  TString className(mcInfo->ClassName());
  if (className.CompareTo("AliMCEvent") != 0) {
    AliError("argument must point to an AliMCEvent !");
    return ;
  }
  
  fMCInfo = (AliMCEvent*) mcInfo ;
}


//__________________________________________________________________________________
void AliCFAcceptanceCuts::AddQAHistograms(TList *qaList) {
  //
  // saves the histograms in a TList
  //

  DefineHistograms();

  qaList->Add(fhCutStatistics);
  qaList->Add(fhCutCorrelation);

  for (Int_t j=0; j<kNStepQA; j++) {
    for(Int_t i=0; i<kNCuts; i++)
      qaList->Add(fhQA[i][j]);
  }
}

//__________________________________________________________________________________
void AliCFAcceptanceCuts::DefineHistograms() {
  //
  // histograms for cut variables, cut statistics and cut correlations
  //
  Int_t color = 2;

  // book cut statistics and cut correlation histograms
  fhCutStatistics = new TH1F(Form("%s_cut_statistics",GetName()),"",kNCuts,0.5,kNCuts+0.5);
  fhCutStatistics->SetLineWidth(2);
  int k = 1;
  fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits ITS") ; k++;
  fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits TPC") ; k++;
  fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits TRD") ; k++;
  fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits TOF") ; k++;
  fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits MUON"); k++;


  fhCutCorrelation = new TH2F(Form("%s_cut_correlation",GetName()),"",kNCuts,0.5,kNCuts+0.5,kNCuts,0.5,kNCuts+0.5);
  fhCutCorrelation->SetLineWidth(2);
  for (k=1; k<=kNCuts; k++) {
    fhCutCorrelation->GetXaxis()->SetBinLabel(k,fhCutStatistics->GetXaxis()->GetBinLabel(k));
    fhCutCorrelation->GetYaxis()->SetBinLabel(k,fhCutStatistics->GetXaxis()->GetBinLabel(k));
  }

  Char_t str[5];
  for (int i=0; i<kNStepQA; i++) {
    if (i==0) snprintf(str,5," ");
    else snprintf(str,5,"_cut");
    fhQA[kCutHitsITS] [i] = new TH1F(Form("%s_HitsITS%s"  ,GetName(),str),"",10,0,10);
    fhQA[kCutHitsTPC] [i] = new TH1F(Form("%s_HitsTPC%s"  ,GetName(),str),"",5,0,5);
    fhQA[kCutHitsTRD] [i] = new TH1F(Form("%s_HitsTRD%s"  ,GetName(),str),"",20,0,20);
    fhQA[kCutHitsTOF] [i] = new TH1F(Form("%s_HitsTOF%s"  ,GetName(),str),"",5,0,5);
    fhQA[kCutHitsMUON][i] = new TH1F(Form("%s_HitsMUON%s" ,GetName(),str),"",5,0,5);
  }
  for(Int_t i=0; i<kNCuts; i++) fhQA[i][1]->SetLineColor(color);
}

//__________________________________________________________________________________
void AliCFAcceptanceCuts::FillHistograms(TObject* obj, Bool_t afterCuts)
{
  //
  // fill the QA histograms
  //
  AliMCParticle* part = dynamic_cast<AliMCParticle *>(obj);
  if (!part) {
    AliError("casting failed");
    return;
  }

  Int_t nHitsITS=0, nHitsTPC=0, nHitsTRD=0, nHitsTOF=0, nHitsMUON=0 ;
  for (Int_t iTrackRef=0; iTrackRef<part->GetNumberOfTrackReferences(); iTrackRef++) {
    AliTrackReference * trackRef = part->GetTrackReference(iTrackRef);
    if(trackRef){
      Int_t detectorId = trackRef->DetectorId();
      switch(detectorId) {
      case AliTrackReference::kITS  : nHitsITS++  ; break ;
      case AliTrackReference::kTPC  : nHitsTPC++  ; break ;
      case AliTrackReference::kTRD  : nHitsTRD++  ; break ;
      case AliTrackReference::kTOF  : nHitsTOF++  ; break ;
      case AliTrackReference::kMUON : nHitsMUON++ ; break ;
      default : break ;
      }
    }
  }
  
  fhQA[kCutHitsITS ][afterCuts]->Fill(nHitsITS );
  fhQA[kCutHitsTPC ][afterCuts]->Fill(nHitsTPC );
  fhQA[kCutHitsTRD ][afterCuts]->Fill(nHitsTRD );
  fhQA[kCutHitsTOF ][afterCuts]->Fill(nHitsTOF );
  fhQA[kCutHitsMUON][afterCuts]->Fill(nHitsMUON);

  // fill cut statistics and cut correlation histograms with information from the bitmap
  if (afterCuts) return;

  // Number of single cuts in this class
  UInt_t ncuts = fBitmap->GetNbits();
  for(UInt_t bit=0; bit<ncuts;bit++) {
    if (!fBitmap->TestBitNumber(bit)) {
      fhCutStatistics->Fill(bit+1);
      for (UInt_t bit2=bit; bit2<ncuts;bit2++) {
        if (!fBitmap->TestBitNumber(bit2)) 
          fhCutCorrelation->Fill(bit+1,bit2+1);
      }
    }
  }
}
Commit	Line	Data
563113d0	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	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	**************************************************************************/
	15
	16	///////////////////////////////////////////////////////////////////////////
	17	// ---- CORRECTION FRAMEWORK ----
	18	// AliCFAcceptanceCuts implementation
	19	// Class to cut on the number of AliTrackReference's
	20	// for each detector
	21	///////////////////////////////////////////////////////////////////////////
	22	// author : R. Vernet (renaud.vernet@cern.ch)
	23	///////////////////////////////////////////////////////////////////////////
	24
	25	#include "AliLog.h"
	26	#include "AliMCParticle.h"
	27	#include "AliCFAcceptanceCuts.h"
107a3100	28	#include "TH1F.h"
	29	#include "TH2F.h"
	30	#include "TList.h"
	31	#include "TBits.h"
563113d0	32
	33	ClassImp(AliCFAcceptanceCuts)
	34
	35	//______________________________
	36	AliCFAcceptanceCuts::AliCFAcceptanceCuts() :
	37	AliCFCutBase(),
	38	fMCInfo(0x0),
	39	fMinNHitITS(0),
	40	fMinNHitTPC(0),
	41	fMinNHitTRD(0),
	42	fMinNHitTOF(0),
107a3100	43	fMinNHitMUON(0),
	44	fhCutStatistics(0x0),
	45	fhCutCorrelation(0x0),
	46	fBitmap(new TBits(0))
563113d0	47	{
	48	//
	49	//ctor
	50	//
107a3100	51	for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=0x0;
563113d0	52	}
	53
	54	//______________________________
	55	AliCFAcceptanceCuts::AliCFAcceptanceCuts(const Char_t* name, const Char_t* title) :
	56	AliCFCutBase(name,title),
	57	fMCInfo(0x0),
	58	fMinNHitITS(0),
	59	fMinNHitTPC(0),
	60	fMinNHitTRD(0),
	61	fMinNHitTOF(0),
107a3100	62	fMinNHitMUON(0),
	63	fhCutStatistics(0x0),
	64	fhCutCorrelation(0x0),
	65	fBitmap(new TBits(0))
563113d0	66	{
	67	//
	68	//ctor
	69	//
107a3100	70	for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=0x0;
563113d0	71	}
	72
	73	//______________________________
	74	AliCFAcceptanceCuts::AliCFAcceptanceCuts(const AliCFAcceptanceCuts& c) :
	75	AliCFCutBase(c),
	76	fMCInfo(c.fMCInfo),
	77	fMinNHitITS(c.fMinNHitITS),
	78	fMinNHitTPC(c.fMinNHitTPC),
	79	fMinNHitTRD(c.fMinNHitTRD),
	80	fMinNHitTOF(c.fMinNHitTOF),
107a3100	81	fMinNHitMUON(c.fMinNHitMUON),
	82	fhCutStatistics(c.fhCutStatistics),
	83	fhCutCorrelation(c.fhCutCorrelation),
	84	fBitmap(c.fBitmap)
563113d0	85	{
	86	//
	87	//copy ctor
	88	//
107a3100	89	for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=c.fhQA[i][j];
563113d0	90	}
	91
	92	//______________________________
	93	AliCFAcceptanceCuts& AliCFAcceptanceCuts::operator=(const AliCFAcceptanceCuts& c)
	94	{
	95	//
	96	// Assignment operator
	97	//
	98	if (this != &c) {
	99	AliCFCutBase::operator=(c) ;
	100	fMCInfo=c.fMCInfo;
	101	fMinNHitITS=c.fMinNHitITS;
	102	fMinNHitTPC=c.fMinNHitTPC;
	103	fMinNHitTRD=c.fMinNHitTRD;
	104	fMinNHitTOF=c.fMinNHitTOF;
	105	fMinNHitMUON=c.fMinNHitMUON;
107a3100	106	fhCutStatistics = c.fhCutStatistics ;
	107	fhCutCorrelation = c.fhCutCorrelation ;
	108	fBitmap = c.fBitmap ;
	109	for (int i=0; i<kNCuts; i++) for (int j=0; j<kNStepQA; j++) fhQA[i][j]=c.fhQA[i][j];
563113d0	110	}
	111	return *this ;
	112	}
	113
107a3100	114	//______________________________________________________________
9eeae5d5	115	Bool_t AliCFAcceptanceCuts::IsSelected(TObject *obj) {
107a3100	116	//
	117	// check if selections on 'obj' are passed
	118	// 'obj' must be an AliMCParticle
	119	//
	120
	121	SelectionBitMap(obj);
	122
	123	if (fIsQAOn) FillHistograms(obj,kFALSE);
	124
	125	for (UInt_t icut=0; icut<fBitmap->GetNbits(); icut++)
	126	if (!fBitmap->TestBitNumber(icut)) return kFALSE ;
	127
	128	if (fIsQAOn) FillHistograms(obj,kTRUE);
	129	return kTRUE;
	130	}
	131
563113d0	132	//______________________________
107a3100	133	void AliCFAcceptanceCuts::SelectionBitMap(TObject* obj) {
563113d0	134	//
	135	// checks the number of track references associated to 'obj'
	136	// 'obj' must be an AliMCParticle
	137	//
	138
107a3100	139	for (UInt_t i=0; i<kNCuts; i++) fBitmap->SetBitNumber(i,kFALSE);
563113d0	140
107a3100	141	if (!obj) return;
563113d0	142	TString className(obj->ClassName());
563113d0	143	if (className.CompareTo("AliMCParticle") != 0) {
107a3100	144	AliError("obj must point to an AliMCParticle !");
107a3100	145	return ;
563113d0	146	}
563113d0	147
107a3100	148	AliMCParticle * part = dynamic_cast<AliMCParticle*>(obj);
107a3100	149	if(!part) return ;
563113d0	150
	151	Int_t nHitsITS=0, nHitsTPC=0, nHitsTRD=0, nHitsTOF=0, nHitsMUON=0 ;
	152	for (Int_t iTrackRef=0; iTrackRef<part->GetNumberOfTrackReferences(); iTrackRef++) {
	153	AliTrackReference * trackRef = part->GetTrackReference(iTrackRef);
	154	if(trackRef){
	155	Int_t detectorId = trackRef->DetectorId();
	156	switch(detectorId) {
	157	case AliTrackReference::kITS : nHitsITS++ ; break ;
	158	case AliTrackReference::kTPC : nHitsTPC++ ; break ;
	159	case AliTrackReference::kTRD : nHitsTRD++ ; break ;
	160	case AliTrackReference::kTOF : nHitsTOF++ ; break ;
	161	case AliTrackReference::kMUON : nHitsMUON++ ; break ;
	162	default : break ;
	163	}
	164	}
	165	}
	166
107a3100	167	Int_t iCutBit = 0;
563113d0	168
107a3100	169	if (nHitsITS >= fMinNHitITS ) fBitmap->SetBitNumber(iCutBit,kTRUE);
	170	iCutBit++;
	171	if (nHitsTPC >= fMinNHitTPC ) fBitmap->SetBitNumber(iCutBit,kTRUE);
	172	iCutBit++;
	173	if (nHitsTRD >= fMinNHitTRD ) fBitmap->SetBitNumber(iCutBit,kTRUE);
	174	iCutBit++;
	175	if (nHitsTOF >= fMinNHitTOF ) fBitmap->SetBitNumber(iCutBit,kTRUE);
	176	iCutBit++;
	177	if (nHitsMUON >= fMinNHitMUON) fBitmap->SetBitNumber(iCutBit,kTRUE);
563113d0	178	}
	179
	180
0c01ae65	181	void AliCFAcceptanceCuts::SetMCEventInfo(const TObject* mcInfo) {
563113d0	182	//
107a3100	183	// Sets pointer to MC event information (AliMCEvent)
563113d0	184	//
	185
	186	if (!mcInfo) {
107a3100	187	AliError("Pointer to AliMCEvent !");
563113d0	188	return;
	189	}
	190
	191	TString className(mcInfo->ClassName());
107a3100	192	if (className.CompareTo("AliMCEvent") != 0) {
107a3100	193	AliError("argument must point to an AliMCEvent !");
563113d0	194	return ;
	195	}
	196
107a3100	197	fMCInfo = (AliMCEvent*) mcInfo ;
	198	}
	199
	200
	201	//__________________________________________________________________________________
	202	void AliCFAcceptanceCuts::AddQAHistograms(TList *qaList) {
	203	//
	204	// saves the histograms in a TList
	205	//
	206
	207	DefineHistograms();
	208
	209	qaList->Add(fhCutStatistics);
	210	qaList->Add(fhCutCorrelation);
	211
	212	for (Int_t j=0; j<kNStepQA; j++) {
	213	for(Int_t i=0; i<kNCuts; i++)
	214	qaList->Add(fhQA[i][j]);
	215	}
	216	}
	217
	218	//__________________________________________________________________________________
	219	void AliCFAcceptanceCuts::DefineHistograms() {
	220	//
	221	// histograms for cut variables, cut statistics and cut correlations
	222	//
	223	Int_t color = 2;
	224
	225	// book cut statistics and cut correlation histograms
	226	fhCutStatistics = new TH1F(Form("%s_cut_statistics",GetName()),"",kNCuts,0.5,kNCuts+0.5);
	227	fhCutStatistics->SetLineWidth(2);
	228	int k = 1;
	229	fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits ITS") ; k++;
	230	fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits TPC") ; k++;
	231	fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits TRD") ; k++;
	232	fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits TOF") ; k++;
	233	fhCutStatistics->GetXaxis()->SetBinLabel(k,"hits MUON"); k++;
	234
	235
	236	fhCutCorrelation = new TH2F(Form("%s_cut_correlation",GetName()),"",kNCuts,0.5,kNCuts+0.5,kNCuts,0.5,kNCuts+0.5);
	237	fhCutCorrelation->SetLineWidth(2);
	238	for (k=1; k<=kNCuts; k++) {
	239	fhCutCorrelation->GetXaxis()->SetBinLabel(k,fhCutStatistics->GetXaxis()->GetBinLabel(k));
	240	fhCutCorrelation->GetYaxis()->SetBinLabel(k,fhCutStatistics->GetXaxis()->GetBinLabel(k));
	241	}
	242
4bc0c5f9	243	Char_t str[5];
107a3100	244	for (int i=0; i<kNStepQA; i++) {
4bc0c5f9	245	if (i==0) snprintf(str,5," ");
4bc0c5f9	246	else snprintf(str,5,"_cut");
107a3100	247	fhQA[kCutHitsITS] [i] = new TH1F(Form("%s_HitsITS%s" ,GetName(),str),"",10,0,10);
107a3100	248	fhQA[kCutHitsTPC] [i] = new TH1F(Form("%s_HitsTPC%s" ,GetName(),str),"",5,0,5);
177a1e83	249	fhQA[kCutHitsTRD] [i] = new TH1F(Form("%s_HitsTRD%s" ,GetName(),str),"",20,0,20);
107a3100	250	fhQA[kCutHitsTOF] [i] = new TH1F(Form("%s_HitsTOF%s" ,GetName(),str),"",5,0,5);
	251	fhQA[kCutHitsMUON][i] = new TH1F(Form("%s_HitsMUON%s" ,GetName(),str),"",5,0,5);
	252	}
	253	for(Int_t i=0; i<kNCuts; i++) fhQA[i][1]->SetLineColor(color);
	254	}
	255
	256	//__________________________________________________________________________________
	257	void AliCFAcceptanceCuts::FillHistograms(TObject* obj, Bool_t afterCuts)
	258	{
	259	//
	260	// fill the QA histograms
	261	//
	262	AliMCParticle* part = dynamic_cast<AliMCParticle *>(obj);
4bc0c5f9	263	if (!part) {
	264	AliError("casting failed");
	265	return;
	266	}
107a3100	267
	268	Int_t nHitsITS=0, nHitsTPC=0, nHitsTRD=0, nHitsTOF=0, nHitsMUON=0 ;
	269	for (Int_t iTrackRef=0; iTrackRef<part->GetNumberOfTrackReferences(); iTrackRef++) {
	270	AliTrackReference * trackRef = part->GetTrackReference(iTrackRef);
	271	if(trackRef){
	272	Int_t detectorId = trackRef->DetectorId();
	273	switch(detectorId) {
	274	case AliTrackReference::kITS : nHitsITS++ ; break ;
	275	case AliTrackReference::kTPC : nHitsTPC++ ; break ;
	276	case AliTrackReference::kTRD : nHitsTRD++ ; break ;
	277	case AliTrackReference::kTOF : nHitsTOF++ ; break ;
	278	case AliTrackReference::kMUON : nHitsMUON++ ; break ;
	279	default : break ;
	280	}
	281	}
	282	}
	283
	284	fhQA[kCutHitsITS ][afterCuts]->Fill(nHitsITS );
	285	fhQA[kCutHitsTPC ][afterCuts]->Fill(nHitsTPC );
	286	fhQA[kCutHitsTRD ][afterCuts]->Fill(nHitsTRD );
	287	fhQA[kCutHitsTOF ][afterCuts]->Fill(nHitsTOF );
	288	fhQA[kCutHitsMUON][afterCuts]->Fill(nHitsMUON);
	289
	290	// fill cut statistics and cut correlation histograms with information from the bitmap
	291	if (afterCuts) return;
	292
	293	// Number of single cuts in this class
	294	UInt_t ncuts = fBitmap->GetNbits();
	295	for(UInt_t bit=0; bit<ncuts;bit++) {
	296	if (!fBitmap->TestBitNumber(bit)) {
	297	fhCutStatistics->Fill(bit+1);
	298	for (UInt_t bit2=bit; bit2<ncuts;bit2++) {
	299	if (!fBitmap->TestBitNumber(bit2))
	300	fhCutCorrelation->Fill(bit+1,bit2+1);
	301	}
	302	}
	303	}
563113d0	304	}