[u/mrichter/AliRoot.git] / PWG0 / multPbPb / AliAnalysisTaskTriggerStudy.cxx

// AliAnalysisTaskTriggerStudy

// Author: Michele Floris, CERN
// TODO:
// - Add chi2/cluster plot for primary, secondaries and fakes


#include "AliAnalysisTaskTriggerStudy.h"
#include "AliESDInputHandler.h"
#include "AliHistoListWrapper.h"
#include "AliAnalysisManager.h"
#include "AliMCEvent.h"
#include "AliStack.h"
#include "TH1I.h"
#include "TH3D.h"
#include "AliMCParticle.h"
#include "AliGenEventHeader.h"
#include "AliESDCentrality.h"

#include <iostream>
#include "AliTriggerAnalysis.h"
#include "AliMultiplicity.h"
#include "TFile.h"
#include "AliLog.h"

using namespace std;

ClassImp(AliAnalysisTaskTriggerStudy)

//const char * AliAnalysisTaskTriggerStudy::kVDNames[] = {"C0SM1","C0SM2","C0VBA","C0VBC","C0OM2"};       
const char * AliAnalysisTaskTriggerStudy::kVDNames[] = {"C0SM1","C0SM2","C0VBA","C0VBC"};//,"C0OM2"};       

AliAnalysisTaskTriggerStudy::AliAnalysisTaskTriggerStudy()
: AliAnalysisTaskSE("TaskTriggerStudy"),
  fESD(0),fHistoList(0),fIsMC(0),fTriggerAnalysis(0),fHistoSuffix(""),fNTrackletsCut(1000000)
{
  // constructor

  DefineOutput(1, AliHistoListWrapper::Class());

}
AliAnalysisTaskTriggerStudy::AliAnalysisTaskTriggerStudy(const char * name)
  : AliAnalysisTaskSE(name),
    fESD(0),fHistoList(0),fIsMC(0),fTriggerAnalysis(0),fHistoSuffix(""),fNTrackletsCut(1000000)
{
  //
  // Standard constructur which should be used
  //

  DefineOutput(1, AliHistoListWrapper::Class());

}

AliAnalysisTaskTriggerStudy::AliAnalysisTaskTriggerStudy(const AliAnalysisTaskTriggerStudy& obj) : 
  AliAnalysisTaskSE(obj) ,fESD (0), fIsMC(0), fTriggerAnalysis(0),fHistoSuffix("")
{
  //copy ctor
  fESD = obj.fESD ;
  fHistoList = obj.fHistoList;
  fTriggerAnalysis = obj.fTriggerAnalysis;
  fHistoSuffix = obj.fHistoSuffix;
}

AliAnalysisTaskTriggerStudy::~AliAnalysisTaskTriggerStudy(){
  // destructor

  if(!AliAnalysisManager::GetAnalysisManager()->IsProofMode()) {
    if(fHistoList) {
      delete fHistoList;
      fHistoList = 0;
    }
    if(fTriggerAnalysis) {
      delete fTriggerAnalysis;
      fHistoList = 0;
    }
  }
  // Histo list should not be destroyed: fListWrapper is owner!

}
void AliAnalysisTaskTriggerStudy::UserCreateOutputObjects()
{
  // Called once
  fHistoList = new AliHistoListWrapper("histoList","histogram list for trigger studies");
  fTriggerAnalysis = new AliTriggerAnalysis();
}


void AliAnalysisTaskTriggerStudy::UserExec(Option_t *)
{
  // User code

  /* PostData(0) is taken care of by AliAnalysisTaskSE */
  PostData(1,fHistoList);

  fESD = dynamic_cast<AliESDEvent*>(fInputEvent);
  if (strcmp(fESD->ClassName(),"AliESDEvent")) {
    AliFatal("Not processing ESDs");
  }

  
  // get the multiplicity object
  const AliMultiplicity* mult = fESD->GetMultiplicity();
  Int_t ntracklets = mult->GetNumberOfTracklets();

  if(ntracklets > fNTrackletsCut) return;

  // Reset histo suffix and fill reference histograms without any suffix
  fHistoSuffix = "";
  GetHistoTracklets("all","All events")->Fill(ntracklets);

  // Fast or in the outer layer  
  Int_t nFastOrOnline  = fTriggerAnalysis->SPDFiredChips(fESD, 1, 0, 2); // offline
  Int_t nFastOrOffline = fTriggerAnalysis->SPDFiredChips(fESD, 0, 0, 2); // online
  
  Bool_t c0sm1 = nFastOrOffline >= 1;
  Bool_t c0sm2 = nFastOrOffline >= 2;
  Bool_t c0sm3 = nFastOrOffline >= 3;
  Bool_t c0sm4 = nFastOrOffline >= 4;
  Bool_t c0sm5 = nFastOrOffline >= 5;
  
  // V0 triggers
  Bool_t c0v0A       = fTriggerAnalysis->IsOfflineTriggerFired(fESD, AliTriggerAnalysis::kV0A);
  Bool_t c0v0C       = fTriggerAnalysis->IsOfflineTriggerFired(fESD, AliTriggerAnalysis::kV0C);
  Bool_t v0AHW     = (fTriggerAnalysis->V0Trigger(fESD, AliTriggerAnalysis::kASide, kTRUE) == AliTriggerAnalysis::kV0BB);// should replay hw trigger
  Bool_t v0CHW     = (fTriggerAnalysis->V0Trigger(fESD, AliTriggerAnalysis::kCSide, kTRUE) == AliTriggerAnalysis::kV0BB);// should replay hw trigger

  // TOF triggers 
  // FIXME: move to triggeranalysis?
  AliESDHeader*h = fESD->GetHeader(); // taken the header from AliESDEvent 
  Bool_t c0OM2 = h->IsTriggerInputFired("0OM2"); // thr >= 2 (input 19)
  Bool_t c0OM3 = h->IsTriggerInputFired("0OM3"); // thr >= 3 (input 20)

  // ZDC triggers
  Bool_t zdcA   = fTriggerAnalysis->ZDCTrigger(fESD, AliTriggerAnalysis::kASide) ;
  Bool_t zdcC   = fTriggerAnalysis->ZDCTrigger(fESD, AliTriggerAnalysis::kCSide) ;
  Bool_t zdcBar = fTriggerAnalysis->ZDCTrigger(fESD, AliTriggerAnalysis::kCentralBarrel) ;

  // Some macros for the online triggers
  Bool_t cMBS2A = c0sm2 && c0v0A;
  Bool_t cMBS2C = c0sm2 && c0v0C;
  Bool_t cMBAC  = c0v0A && c0v0C;
  

  Bool_t vdArray[kNVDEntries];
  vdArray[kVDC0MBS1] = c0sm1;
  vdArray[kVDC0MBS2] = c0sm2;
  vdArray[kVDC0VBA]  = c0v0A;
  vdArray[kVDC0VBC]  = c0v0C;
  //vdArray[kVDC0OM2]  = c0OM2;

  FillTriggerOverlaps("All", "All Events",vdArray);
  
  // Fill some combination of trigger classes
  Bool_t cmbs1aOnline = fESD->IsTriggerClassFired("CMBS1A-B-NOPF-ALL");
  Bool_t cmbs1cOnline = fESD->IsTriggerClassFired("CMBS1C-B-NOPF-ALL");
  Bool_t cmbacOnline  = fESD->IsTriggerClassFired("CMBAC-B-NOPF-ALL");
  
  if (cmbs1aOnline || cmbs1cOnline ||cmbacOnline) GetHistoTracklets("TwoOutOfThree" ,"Events 2-out-of-3 online" )->Fill(ntracklets);


  // loop over trigger classes in the event
  TObjArray * tokens = 0;
  if(fIsMC) {
    // in case of montecarlo I override the trigger class
    tokens = new TObjArray;
    tokens->SetOwner();
    //    tokens->Add(new TObjString("CINT1B-ABCE-NOPF-ALL")); 
    tokens->Add(new TObjString("MC")); 
  }
  else {  
    TString trgClasses = fESD->GetFiredTriggerClasses();
    tokens = trgClasses.Tokenize("  ");
  }
  TIterator  * iter = (TIterator*) tokens->MakeIterator();
  
  TString classes = fESD->GetFiredTriggerClasses();

  // if (classes.Contains("SMH")) {
  //    tokens->Print();
  //    cout << classes.Data() << endl;
  // }
  //  iter->Reset();
  Int_t itoken = 0;
  TObjString * tok=0;
  while((tok = (TObjString*) iter->Next())){
    // clean up trigger name
    TString trg = tok->GetString();
    trg.Strip(TString::kTrailing, ' ');
    trg.Strip(TString::kLeading, ' ');

    fHistoSuffix = "_";
    fHistoSuffix += trg;

    //    cout << itoken++ << " " << trg.Data() << endl;
    // continue;
    //    if (trg.Contains("SMH")) cout << itoken++ << " " << trg.Data() << endl;
    
    // Fill tracklets 
    GetHistoTracklets("All" ,"Events no offline trigger" )->Fill(ntracklets);


    // Fill histograms mismatchs
    // TODO: check mismatch trigger class 
    if(nFastOrOffline != nFastOrOnline) {
      GetHistoTracklets("mismatchingFastOr", "Events where fast or offline differs from fast-or online")->Fill(ntracklets);
    }
    
    if (c0v0A != v0AHW){
      GetHistoTracklets("mismatchingV0A", "Events where V0A offline differs from V0A online")->Fill(ntracklets);
    }
    
    if (c0v0C != v0CHW){
      GetHistoTracklets("mismatchingV0C", "Events where V0C offline differs from V0C online")->Fill(ntracklets);
    }    
    
    // Fill a tracklet histo for each trigger type
    if(c0sm1)  GetHistoTracklets("c0sm1" ,"Events were trigger c0sm1 fired" )->Fill(ntracklets);
    if(c0sm2)  GetHistoTracklets("c0sm2" ,"Events were trigger c0sm2 fired" )->Fill(ntracklets);
    if(c0sm3)  GetHistoTracklets("c0sm3" ,"Events were trigger c0sm3 fired" )->Fill(ntracklets);
    if(c0sm4)  GetHistoTracklets("c0sm4" ,"Events were trigger c0sm4 fired" )->Fill(ntracklets);
    if(c0sm5)  GetHistoTracklets("c0sm5" ,"Events were trigger c0sm5 fired" )->Fill(ntracklets);
    if(c0OM2)  GetHistoTracklets("c0OM2" ,"Events were trigger c0OM2 fired" )->Fill(ntracklets);
    if(c0OM3)  GetHistoTracklets("c0OM3" ,"Events were trigger c0OM3 fired" )->Fill(ntracklets);
    if(c0v0A)  GetHistoTracklets("c0v0A" ,"Events were trigger c0v0A fired" )->Fill(ntracklets);
    if(c0v0C)  GetHistoTracklets("c0v0C" ,"Events were trigger c0v0C fired" )->Fill(ntracklets);
    if(cMBS2A) GetHistoTracklets("cMBS2A","Events were trigger cMBS2A fired")->Fill(ntracklets);
    if(cMBS2C) GetHistoTracklets("cMBS2C","Events were trigger cMBS2C fired")->Fill(ntracklets);
    if(cMBAC ) GetHistoTracklets("cMBAC", "Events were trigger cMBAC  fired")->Fill(ntracklets);
    if(zdcA )  GetHistoTracklets("cZDCA", "Events were trigger cZDCA  fired")->Fill(ntracklets);
    if(zdcC )  GetHistoTracklets("cZDCC", "Events were trigger cZDCC  fired")->Fill(ntracklets);
    if(zdcBar) GetHistoTracklets("cZDCBar","Events were trigger cZDCB  fired")->Fill(ntracklets);
    //  if() GetHistoTracklets("","Events were trigger  fired");
    
    // Fill trigger overlaps
    FillTriggerOverlaps("All", "All Events in trigger class",vdArray);

  }
  delete tokens;
    
    // if (fIsMC) {
    

  //   if (!fMCEvent) {
  //     AliError("No MC info found");
  //   } else {
      
  //     //loop on the MC event
  //     //      Int_t nMCTracks = fMCEvent->GetNumberOfTracks();
  //     Int_t offset    = fMCEvent->GetPrimaryOffset();
  //     Int_t nMCTracks = fMCEvent->GetNumberOfPrimaries()+offset;
  //     for (Int_t ipart=offset; ipart<nMCTracks; ipart++) { 
	
  // 	AliMCParticle *mcPart  = (AliMCParticle*)fMCEvent->GetTrack(ipart);
	
  // 	// We don't care about neutrals and non-physical primaries
  // 	if(mcPart->Charge() == 0) continue;

  // PHYSICAL PRIMARY
  // 	// Get MC vertex
    // 	TArrayF   vertex;
  // 	fMCEvent->GenEventHeader()->PrimaryVertex(vertex);
  // 	Float_t zv = vertex[2];
  // 	//	Float_t zv = vtxESD->GetZ();
  // 	// Fill generated histo
  // 	hTracks[AliAnalysisMultPbTrackHistoManager::kHistoGen]->Fill(mcPart->Pt(),mcPart->Eta(),zv);
	
  //     }
  //   }
  // }
  

}

void   AliAnalysisTaskTriggerStudy::Terminate(Option_t *){
  // terminate
  // Save output in a more friendly format
  fHistoList = dynamic_cast<AliHistoListWrapper*> (GetOutputData(1));
  if (!fHistoList){
    Printf("ERROR: fHistoList not available");
    return;
  }
  TFile * f = new TFile("trigger_study.root", "recreate");
  fHistoList->GetList()->Write();
  f->Close();

}

TH1 *   AliAnalysisTaskTriggerStudy::GetHistoTracklets(const char * name, const char * title){
  // Book histo of events vs ntracklets, if needed

  TString hname = "hTracklets_";
  hname+=name;  
  hname+=fHistoSuffix;
  TH1 * h = (TH1*) fHistoList->GetList()->FindObject(hname.Data());
  
  if(!h) {
    AliInfo(Form("Booking histo %s",hname.Data()));
    Bool_t oldStatus = TH1::AddDirectoryStatus();
    TH1::AddDirectory(kFALSE);
    // h = new TH1F (hname.Data(), title, 1000, -0.5, 10000-0.5);
    h = new TH1F (hname.Data(), title, 4000, -0.5, 4000-0.5);
    h->Sumw2();
    h->SetXTitle("ntracklets");
    fHistoList->GetList()->Add(h);
    TH1::AddDirectory(oldStatus);
  }
  return h;
}

void AliAnalysisTaskTriggerStudy::FillTriggerOverlaps (const char * name, const char * title, Bool_t * vdArray){
  //Fills a histo with the different trigger statistics in a venn like diagramm. Books it if needed.

  // Get or book histo
  TString hname = "hTrigStat_";
  hname+=name;  
  hname+=fHistoSuffix;
  TH1 * h = (TH1*) fHistoList->GetList()->FindObject(hname.Data());
  
  if(!h) {
    AliInfo(Form("Booking histo %s",hname.Data()));
    Bool_t oldStatus = TH1::AddDirectoryStatus();
    TH1::AddDirectory(kFALSE);
    Int_t nbins = 0;
    for(Int_t ientry = 0; ientry < kNVDEntries; ientry++){
      nbins = nbins | (1<<ientry);
    }
    //    cout << "NBINS " << nbins << endl;
    nbins = nbins+1;
    h = new TH1I (hname, title, nbins, -0.5, nbins-0.5);
    fHistoList->GetList()->Add(h);
    TH1::AddDirectory(oldStatus);
  
    // we look at the combinations of n triggers
    // We set a bit for each trigger to fill the diagram
    // This is much simpler and faster than any recursive function
    h->GetXaxis()->SetBinLabel(1,"NONE"); 
    for(Int_t ibin = 1; ibin < nbins; ibin++){
      TString binname = "";
      Bool_t first = kTRUE;
      for(Int_t ivdentry = 0; ivdentry < kNVDEntries; ivdentry++){
	if (ibin & (1<<ivdentry)) {
	  if(!first) binname += " & ";
	  binname += kVDNames[ivdentry];
	  first=kFALSE;
	}
      }
      h->GetXaxis()->SetBinLabel(ibin+1,binname.Data());
    }
    
  }

  UInt_t mask = 0;
  for(Int_t ivdentry = 0; ivdentry < kNVDEntries; ivdentry++){
    if(vdArray[ivdentry]) {
      mask  = mask | (1<<ivdentry);
      //      cout << " 1 "   ;
    } //else cout << " 0 ";
  }
  //  cout << hex << " = " << mask << endl;
  
  h->Fill(mask);

}
Commit	Line	Data
9d173aad	1	// AliAnalysisTaskTriggerStudy
	2
	3	// Author: Michele Floris, CERN
	4	// TODO:
	5	// - Add chi2/cluster plot for primary, secondaries and fakes
	6
	7
	8	#include "AliAnalysisTaskTriggerStudy.h"
	9	#include "AliESDInputHandler.h"
	10	#include "AliHistoListWrapper.h"
	11	#include "AliAnalysisManager.h"
	12	#include "AliMCEvent.h"
	13	#include "AliStack.h"
	14	#include "TH1I.h"
	15	#include "TH3D.h"
	16	#include "AliMCParticle.h"
	17	#include "AliGenEventHeader.h"
	18	#include "AliESDCentrality.h"
	19
	20	#include <iostream>
	21	#include "AliTriggerAnalysis.h"
	22	#include "AliMultiplicity.h"
	23	#include "TFile.h"
	24	#include "AliLog.h"
	25
	26	using namespace std;
	27
	28	ClassImp(AliAnalysisTaskTriggerStudy)
	29
7f5f2e0c	30	//const char * AliAnalysisTaskTriggerStudy::kVDNames[] = {"C0SM1","C0SM2","C0VBA","C0VBC","C0OM2"};
7f5f2e0c	31	const char * AliAnalysisTaskTriggerStudy::kVDNames[] = {"C0SM1","C0SM2","C0VBA","C0VBC"};//,"C0OM2"};
52d405b5	32
9d173aad	33	AliAnalysisTaskTriggerStudy::AliAnalysisTaskTriggerStudy()
9d173aad	34	: AliAnalysisTaskSE("TaskTriggerStudy"),
52d405b5	35	fESD(0),fHistoList(0),fIsMC(0),fTriggerAnalysis(0),fHistoSuffix(""),fNTrackletsCut(1000000)
9d173aad	36	{
	37	// constructor
	38
	39	DefineOutput(1, AliHistoListWrapper::Class());
	40
	41	}
	42	AliAnalysisTaskTriggerStudy::AliAnalysisTaskTriggerStudy(const char * name)
	43	: AliAnalysisTaskSE(name),
52d405b5	44	fESD(0),fHistoList(0),fIsMC(0),fTriggerAnalysis(0),fHistoSuffix(""),fNTrackletsCut(1000000)
9d173aad	45	{
	46	//
	47	// Standard constructur which should be used
	48	//
	49
	50	DefineOutput(1, AliHistoListWrapper::Class());
	51
	52	}
	53
	54	AliAnalysisTaskTriggerStudy::AliAnalysisTaskTriggerStudy(const AliAnalysisTaskTriggerStudy& obj) :
a82dc581	55	AliAnalysisTaskSE(obj) ,fESD (0), fIsMC(0), fTriggerAnalysis(0),fHistoSuffix("")
9d173aad	56	{
	57	//copy ctor
	58	fESD = obj.fESD ;
	59	fHistoList = obj.fHistoList;
	60	fTriggerAnalysis = obj.fTriggerAnalysis;
a82dc581	61	fHistoSuffix = obj.fHistoSuffix;
9d173aad	62	}
	63
	64	AliAnalysisTaskTriggerStudy::~AliAnalysisTaskTriggerStudy(){
	65	// destructor
	66
	67	if(!AliAnalysisManager::GetAnalysisManager()->IsProofMode()) {
	68	if(fHistoList) {
	69	delete fHistoList;
	70	fHistoList = 0;
	71	}
	72	if(fTriggerAnalysis) {
	73	delete fTriggerAnalysis;
	74	fHistoList = 0;
	75	}
	76	}
	77	// Histo list should not be destroyed: fListWrapper is owner!
	78
	79	}
	80	void AliAnalysisTaskTriggerStudy::UserCreateOutputObjects()
	81	{
	82	// Called once
	83	fHistoList = new AliHistoListWrapper("histoList","histogram list for trigger studies");
	84	fTriggerAnalysis = new AliTriggerAnalysis();
	85	}
	86
	87
	88	void AliAnalysisTaskTriggerStudy::UserExec(Option_t *)
	89	{
	90	// User code
	91
	92	/* PostData(0) is taken care of by AliAnalysisTaskSE */
	93	PostData(1,fHistoList);
	94
	95	fESD = dynamic_cast<AliESDEvent*>(fInputEvent);
	96	if (strcmp(fESD->ClassName(),"AliESDEvent")) {
	97	AliFatal("Not processing ESDs");
	98	}
	99
9d173aad	100
	101	// get the multiplicity object
	102	const AliMultiplicity* mult = fESD->GetMultiplicity();
	103	Int_t ntracklets = mult->GetNumberOfTracklets();
	104
52d405b5	105	if(ntracklets > fNTrackletsCut) return;
52d405b5	106
a82dc581	107	// Reset histo suffix and fill reference histograms without any suffix
a82dc581	108	fHistoSuffix = "";
9d173aad	109	GetHistoTracklets("all","All events")->Fill(ntracklets);
	110
	111	// Fast or in the outer layer
	112	Int_t nFastOrOnline = fTriggerAnalysis->SPDFiredChips(fESD, 1, 0, 2); // offline
	113	Int_t nFastOrOffline = fTriggerAnalysis->SPDFiredChips(fESD, 0, 0, 2); // online
9d173aad	114
	115	Bool_t c0sm1 = nFastOrOffline >= 1;
	116	Bool_t c0sm2 = nFastOrOffline >= 2;
	117	Bool_t c0sm3 = nFastOrOffline >= 3;
	118	Bool_t c0sm4 = nFastOrOffline >= 4;
	119	Bool_t c0sm5 = nFastOrOffline >= 5;
	120
	121	// V0 triggers
	122	Bool_t c0v0A = fTriggerAnalysis->IsOfflineTriggerFired(fESD, AliTriggerAnalysis::kV0A);
	123	Bool_t c0v0C = fTriggerAnalysis->IsOfflineTriggerFired(fESD, AliTriggerAnalysis::kV0C);
	124	Bool_t v0AHW = (fTriggerAnalysis->V0Trigger(fESD, AliTriggerAnalysis::kASide, kTRUE) == AliTriggerAnalysis::kV0BB);// should replay hw trigger
a82dc581	125	Bool_t v0CHW = (fTriggerAnalysis->V0Trigger(fESD, AliTriggerAnalysis::kCSide, kTRUE) == AliTriggerAnalysis::kV0BB);// should replay hw trigger
9d173aad	126
9d173aad	127	// TOF triggers
52d405b5	128	// FIXME: move to triggeranalysis?
	129	AliESDHeader*h = fESD->GetHeader(); // taken the header from AliESDEvent
	130	Bool_t c0OM2 = h->IsTriggerInputFired("0OM2"); // thr >= 2 (input 19)
	131	Bool_t c0OM3 = h->IsTriggerInputFired("0OM3"); // thr >= 3 (input 20)
9d173aad	132
5ec1c2f5	133	// ZDC triggers
	134	Bool_t zdcA = fTriggerAnalysis->ZDCTrigger(fESD, AliTriggerAnalysis::kASide) ;
	135	Bool_t zdcC = fTriggerAnalysis->ZDCTrigger(fESD, AliTriggerAnalysis::kCSide) ;
	136	Bool_t zdcBar = fTriggerAnalysis->ZDCTrigger(fESD, AliTriggerAnalysis::kCentralBarrel) ;
	137
9d173aad	138	// Some macros for the online triggers
	139	Bool_t cMBS2A = c0sm2 && c0v0A;
	140	Bool_t cMBS2C = c0sm2 && c0v0C;
	141	Bool_t cMBAC = c0v0A && c0v0C;
52d405b5	142
9d173aad	143
52d405b5	144	Bool_t vdArray[kNVDEntries];
5ec1c2f5	145	vdArray[kVDC0MBS1] = c0sm1;
52d405b5	146	vdArray[kVDC0MBS2] = c0sm2;
	147	vdArray[kVDC0VBA] = c0v0A;
	148	vdArray[kVDC0VBC] = c0v0C;
7f5f2e0c	149	//vdArray[kVDC0OM2] = c0OM2;
52d405b5	150
52d405b5	151	FillTriggerOverlaps("All", "All Events",vdArray);
a82dc581	152
7f5f2e0c	153	// Fill some combination of trigger classes
	154	Bool_t cmbs1aOnline = fESD->IsTriggerClassFired("CMBS1A-B-NOPF-ALL");
	155	Bool_t cmbs1cOnline = fESD->IsTriggerClassFired("CMBS1C-B-NOPF-ALL");
	156	Bool_t cmbacOnline = fESD->IsTriggerClassFired("CMBAC-B-NOPF-ALL");
	157
	158	if (cmbs1aOnline \|\| cmbs1cOnline \|\|cmbacOnline) GetHistoTracklets("TwoOutOfThree" ,"Events 2-out-of-3 online" )->Fill(ntracklets);
	159
a82dc581	160
	161	// loop over trigger classes in the event
	162	TObjArray * tokens = 0;
	163	if(fIsMC) {
	164	// in case of montecarlo I override the trigger class
	165	tokens = new TObjArray;
	166	tokens->SetOwner();
	167	// tokens->Add(new TObjString("CINT1B-ABCE-NOPF-ALL"));
	168	tokens->Add(new TObjString("MC"));
	169	}
	170	else {
	171	TString trgClasses = fESD->GetFiredTriggerClasses();
7f5f2e0c	172	tokens = trgClasses.Tokenize(" ");
a82dc581	173	}
7f5f2e0c	174	TIterator * iter = (TIterator*) tokens->MakeIterator();
	175
	176	TString classes = fESD->GetFiredTriggerClasses();
	177
	178	// if (classes.Contains("SMH")) {
	179	// tokens->Print();
	180	// cout << classes.Data() << endl;
	181	// }
	182	// iter->Reset();
	183	Int_t itoken = 0;
	184	TObjString * tok=0;
	185	while((tok = (TObjString*) iter->Next())){
a82dc581	186	// clean up trigger name
	187	TString trg = tok->GetString();
	188	trg.Strip(TString::kTrailing, ' ');
	189	trg.Strip(TString::kLeading, ' ');
	190
	191	fHistoSuffix = "_";
	192	fHistoSuffix += trg;
	193
7f5f2e0c	194	// cout << itoken++ << " " << trg.Data() << endl;
	195	// continue;
	196	// if (trg.Contains("SMH")) cout << itoken++ << " " << trg.Data() << endl;
	197
	198	// Fill tracklets
	199	GetHistoTracklets("All" ,"Events no offline trigger" )->Fill(ntracklets);
	200
	201
a82dc581	202	// Fill histograms mismatchs
	203	// TODO: check mismatch trigger class
	204	if(nFastOrOffline != nFastOrOnline) {
	205	GetHistoTracklets("mismatchingFastOr", "Events where fast or offline differs from fast-or online")->Fill(ntracklets);
	206	}
	207
	208	if (c0v0A != v0AHW){
	209	GetHistoTracklets("mismatchingV0A", "Events where V0A offline differs from V0A online")->Fill(ntracklets);
	210	}
	211
	212	if (c0v0C != v0CHW){
	213	GetHistoTracklets("mismatchingV0C", "Events where V0C offline differs from V0C online")->Fill(ntracklets);
	214	}
	215
	216	// Fill a tracklet histo for each trigger type
	217	if(c0sm1) GetHistoTracklets("c0sm1" ,"Events were trigger c0sm1 fired" )->Fill(ntracklets);
	218	if(c0sm2) GetHistoTracklets("c0sm2" ,"Events were trigger c0sm2 fired" )->Fill(ntracklets);
	219	if(c0sm3) GetHistoTracklets("c0sm3" ,"Events were trigger c0sm3 fired" )->Fill(ntracklets);
	220	if(c0sm4) GetHistoTracklets("c0sm4" ,"Events were trigger c0sm4 fired" )->Fill(ntracklets);
	221	if(c0sm5) GetHistoTracklets("c0sm5" ,"Events were trigger c0sm5 fired" )->Fill(ntracklets);
	222	if(c0OM2) GetHistoTracklets("c0OM2" ,"Events were trigger c0OM2 fired" )->Fill(ntracklets);
	223	if(c0OM3) GetHistoTracklets("c0OM3" ,"Events were trigger c0OM3 fired" )->Fill(ntracklets);
	224	if(c0v0A) GetHistoTracklets("c0v0A" ,"Events were trigger c0v0A fired" )->Fill(ntracklets);
	225	if(c0v0C) GetHistoTracklets("c0v0C" ,"Events were trigger c0v0C fired" )->Fill(ntracklets);
	226	if(cMBS2A) GetHistoTracklets("cMBS2A","Events were trigger cMBS2A fired")->Fill(ntracklets);
	227	if(cMBS2C) GetHistoTracklets("cMBS2C","Events were trigger cMBS2C fired")->Fill(ntracklets);
5ec1c2f5	228	if(cMBAC ) GetHistoTracklets("cMBAC", "Events were trigger cMBAC fired")->Fill(ntracklets);
	229	if(zdcA ) GetHistoTracklets("cZDCA", "Events were trigger cZDCA fired")->Fill(ntracklets);
	230	if(zdcC ) GetHistoTracklets("cZDCC", "Events were trigger cZDCC fired")->Fill(ntracklets);
	231	if(zdcBar) GetHistoTracklets("cZDCBar","Events were trigger cZDCB fired")->Fill(ntracklets);
a82dc581	232	// if() GetHistoTracklets("","Events were trigger fired");
52d405b5	233
	234	// Fill trigger overlaps
	235	FillTriggerOverlaps("All", "All Events in trigger class",vdArray);
	236
a82dc581	237	}
7f5f2e0c	238	delete tokens;
a82dc581	239
a82dc581	240	// if (fIsMC) {
9d173aad	241
	242
	243	// if (!fMCEvent) {
	244	// AliError("No MC info found");
	245	// } else {
	246
	247	// //loop on the MC event
	248	// // Int_t nMCTracks = fMCEvent->GetNumberOfTracks();
	249	// Int_t offset = fMCEvent->GetPrimaryOffset();
	250	// Int_t nMCTracks = fMCEvent->GetNumberOfPrimaries()+offset;
	251	// for (Int_t ipart=offset; ipart<nMCTracks; ipart++) {
	252
	253	// AliMCParticle mcPart = (AliMCParticle)fMCEvent->GetTrack(ipart);
	254
	255	// // We don't care about neutrals and non-physical primaries
	256	// if(mcPart->Charge() == 0) continue;
	257
e0376287	258	// PHYSICAL PRIMARY
9d173aad	259	// // Get MC vertex
e0376287	260	// TArrayF vertex;
9d173aad	261	// fMCEvent->GenEventHeader()->PrimaryVertex(vertex);
	262	// Float_t zv = vertex[2];
	263	// // Float_t zv = vtxESD->GetZ();
	264	// // Fill generated histo
	265	// hTracks[AliAnalysisMultPbTrackHistoManager::kHistoGen]->Fill(mcPart->Pt(),mcPart->Eta(),zv);
	266
	267	// }
	268	// }
	269	// }
	270
	271
	272
	273
	274	}
	275
	276	void AliAnalysisTaskTriggerStudy::Terminate(Option_t *){
	277	// terminate
	278	// Save output in a more friendly format
	279	fHistoList = dynamic_cast<AliHistoListWrapper*> (GetOutputData(1));
	280	if (!fHistoList){
	281	Printf("ERROR: fHistoList not available");
	282	return;
	283	}
	284	TFile * f = new TFile("trigger_study.root", "recreate");
	285	fHistoList->GetList()->Write();
	286	f->Close();
	287
	288	}
	289
	290	TH1 * AliAnalysisTaskTriggerStudy::GetHistoTracklets(const char * name, const char * title){
a82dc581	291	// Book histo of events vs ntracklets, if needed
9d173aad	292
a82dc581	293	TString hname = "hTracklets_";
	294	hname+=name;
	295	hname+=fHistoSuffix;
9d173aad	296	TH1 * h = (TH1*) fHistoList->GetList()->FindObject(hname.Data());
	297
	298	if(!h) {
	299	AliInfo(Form("Booking histo %s",hname.Data()));
	300	Bool_t oldStatus = TH1::AddDirectoryStatus();
	301	TH1::AddDirectory(kFALSE);
7f5f2e0c	302	// h = new TH1F (hname.Data(), title, 1000, -0.5, 10000-0.5);
7f5f2e0c	303	h = new TH1F (hname.Data(), title, 4000, -0.5, 4000-0.5);
9d173aad	304	h->Sumw2();
	305	h->SetXTitle("ntracklets");
	306	fHistoList->GetList()->Add(h);
	307	TH1::AddDirectory(oldStatus);
	308	}
	309	return h;
	310	}
	311
52d405b5	312	void AliAnalysisTaskTriggerStudy::FillTriggerOverlaps (const char * name, const char * title, Bool_t * vdArray){
a82dc581	313	//Fills a histo with the different trigger statistics in a venn like diagramm. Books it if needed.
	314
	315	// Get or book histo
	316	TString hname = "hTrigStat_";
	317	hname+=name;
	318	hname+=fHistoSuffix;
	319	TH1 * h = (TH1*) fHistoList->GetList()->FindObject(hname.Data());
	320
	321	if(!h) {
	322	AliInfo(Form("Booking histo %s",hname.Data()));
	323	Bool_t oldStatus = TH1::AddDirectoryStatus();
	324	TH1::AddDirectory(kFALSE);
52d405b5	325	Int_t nbins = 0;
	326	for(Int_t ientry = 0; ientry < kNVDEntries; ientry++){
	327	nbins = nbins \| (1<<ientry);
	328	}
7f5f2e0c	329	// cout << "NBINS " << nbins << endl;
7f5f2e0c	330	nbins = nbins+1;
52d405b5	331	h = new TH1I (hname, title, nbins, -0.5, nbins-0.5);
a82dc581	332	fHistoList->GetList()->Add(h);
a82dc581	333	TH1::AddDirectory(oldStatus);
52d405b5	334
	335	// we look at the combinations of n triggers
	336	// We set a bit for each trigger to fill the diagram
	337	// This is much simpler and faster than any recursive function
	338	h->GetXaxis()->SetBinLabel(1,"NONE");
	339	for(Int_t ibin = 1; ibin < nbins; ibin++){
	340	TString binname = "";
	341	Bool_t first = kTRUE;
	342	for(Int_t ivdentry = 0; ivdentry < kNVDEntries; ivdentry++){
	343	if (ibin & (1<<ivdentry)) {
	344	if(!first) binname += " & ";
	345	binname += kVDNames[ivdentry];
	346	first=kFALSE;
	347	}
	348	}
	349	h->GetXaxis()->SetBinLabel(ibin+1,binname.Data());
	350	}
	351
a82dc581	352	}
a82dc581	353
52d405b5	354	UInt_t mask = 0;
	355	for(Int_t ivdentry = 0; ivdentry < kNVDEntries; ivdentry++){
	356	if(vdArray[ivdentry]) {
	357	mask = mask \| (1<<ivdentry);
	358	// cout << " 1 " ;
	359	} //else cout << " 0 ";
	360	}
	361	// cout << hex << " = " << mask << endl;
a82dc581	362
52d405b5	363	h->Fill(mask);
9d173aad	364
a82dc581	365	}