TList * AliAnaCaloTrackCorrMaker::FillAndGetAODBranchList()
{
- // Get any new output AOD branches from analysis and put them in a list
- // The list is filled in the maker, and new branch passed to the analysis frame
- // AliAnalysisTaskCaloTrackCorrelation
+ // Get any new output AOD branches from analysis and put them in a list
+ // The list is filled in the maker, and new branch passed to the analysis frame
+ // AliAnalysisTaskCaloTrackCorrelation
- TList *aodBranchList = fReader->GetAODBranchList() ;
+ TList *aodBranchList = fReader->GetAODBranchList() ;
+
+ for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
+ {
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ if(ana->NewOutputAOD()) aodBranchList->Add(ana->GetCreateOutputAODBranch());
+ }
+
+ return aodBranchList ;
- for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
- {
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
- if(ana->NewOutputAOD()) aodBranchList->Add(ana->GetCreateOutputAODBranch());
- }
-
- return aodBranchList ;
-
}
//_________________________________________________________
TList * AliAnaCaloTrackCorrMaker::GetListOfAnalysisCuts()
{
- // Get the list of the cuts used for the analysis
- // The list is filled in the maker, called by the task in LocalInit() and posted there
+ // Get the list of the cuts used for the analysis
+ // The list is filled in the maker, called by the task in LocalInit() and posted there
- for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
- {
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
- TObjString * objstring = ana->GetAnalysisCuts();
-
- if(objstring)fCuts->Add(objstring);
- }
+ for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
+ {
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ TObjString * objstring = ana->GetAnalysisCuts();
+
+ if(objstring)fCuts->Add(objstring);
+ }
- return fCuts ;
+ return fCuts ;
}
printf("Number of analysis tasks = %d\n", fAnalysisContainer->GetEntries()) ;
if(!strcmp("all",opt))
- {
- printf("Print analysis Tasks settings :\n") ;
- for(Int_t iana = 0; iana<fAnalysisContainer->GetEntries(); iana++)
{
- ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana))->Print("");
- }
-
- printf("Print analysis Reader settings :\n") ;
- fReader->Print("");
- printf("Print analysis Calorimeter Utils settings :\n") ;
- fCaloUtils->Print("");
-
+ printf("Print analysis Tasks settings :\n") ;
+ for(Int_t iana = 0; iana<fAnalysisContainer->GetEntries(); iana++)
+ {
+ ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana))->Print("");
+ }
+
+ printf("Print analysis Reader settings :\n") ;
+ fReader->Print("");
+ printf("Print analysis Calorimeter Utils settings :\n") ;
+ fCaloUtils->Print("");
+
}
}
//Process analysis for this event
if(fMakeHisto && !fOutputContainer)
- {
- printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Histograms not initialized\n");
- abort();
- }
-
+ {
+ printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Histograms not initialized\n");
+ abort();
+ }
+
if(fAnaDebug >= 0 )
- {
- printf("*** AliAnaCaloTrackCorrMaker::ProcessEvent() Event %d *** \n",iEntry);
- if(fAnaDebug > 1 )
{
- printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Current File Name : %s\n", currentFileName);
- //printf("fAODBranchList %p, entries %d\n",fAODBranchList,fAODBranchList->GetEntries());
- }
- }
+ printf("*** AliAnaCaloTrackCorrMaker::ProcessEvent() Event %d *** \n",iEntry);
+ if(fAnaDebug > 1 )
+ {
+ printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Current File Name : %s\n", currentFileName);
+ //printf("fAODBranchList %p, entries %d\n",fAODBranchList,fAODBranchList->GetEntries());
+ }
+ }
//Each event needs an empty branch
TList * aodList = fReader->GetAODBranchList();
//Set the AODB calibration, bad channels etc. parameters at least once
fCaloUtils->AccessOADB(fReader->GetInputEvent());
-
//Tell the reader to fill the data in the 3 detector lists
Bool_t ok = fReader->FillInputEvent(iEntry, currentFileName);
if(!ok)
{
- if(fAnaDebug >= 1 )printf("*** Skip event *** %d \n",iEntry);
+ if(fAnaDebug >= 1 )printf("*** Skip event *** %d \n",iEntry);
fReader->ResetLists();
- return ;
+ return ;
}
//Magic line to write events to file
if (!outputList)
{
- Error("Terminate", "No output list");
- return;
+ Error("Terminate", "No output list");
+ return;
}
-
+
for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
- {
-
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
- if(ana->MakePlotsOn())ana->Terminate(outputList);
-
- }//Loop on analysis defined
+ {
+
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ if(ana->MakePlotsOn())ana->Terminate(outputList);
+
+ }//Loop on analysis defined
}
fEventTriggerAtSE, fEventTriggerMask,fMixEventTriggerMask);
if(fComparePtHardAndClusterPt)
- printf("Compare jet pt and pt hard to accept event, factor = %2.2f",fPtHardAndJetPtFactor);
+ printf("Compare jet pt and pt hard to accept event, factor = %2.2f",fPtHardAndJetPtFactor);
if(fComparePtHardAndClusterPt)
- printf("Compare cluster pt and pt hard to accept event, factor = %2.2f",fPtHardAndClusterPtFactor);
+ printf("Compare cluster pt and pt hard to accept event, factor = %2.2f",fPtHardAndClusterPtFactor);
printf("Read Kine from, stack? %d, AOD ? %d \n", fReadStack, fReadAODMCParticles) ;
printf("Delta AOD File Name = %s\n", fDeltaAODFileName.Data()) ;
//fCurrentFileName = TString(currentFileName);
if(!fInputEvent)
{
- if(fDebug >= 0) printf("AliCaloTrackReader::FillInputEvent() - Input event not available, skip event analysis\n");
- return kFALSE;
+ if(fDebug >= 0) printf("AliCaloTrackReader::FillInputEvent() - Input event not available, skip event analysis\n");
+ return kFALSE;
}
//Select events only fired by a certain trigger configuration if it is provided
Int_t eventType = 0;
if(fInputEvent->GetHeader())
- eventType = ((AliVHeader*)fInputEvent->GetHeader())->GetEventType();
+ eventType = ((AliVHeader*)fInputEvent->GetHeader())->GetEventType();
if (GetFiredTriggerClasses().Contains("FAST") && !GetFiredTriggerClasses().Contains("ALL") && !fAcceptFastCluster)
{
if ( !IsEMCALCluster(clus)) continue ;
- //Skip clusters with too low energy to be triggering
+ //Skip clusters with too low energy to be triggering
if ( clus->E() < fTriggerEventThreshold / 2. ) continue ;
Float_t frac = -1;
Double_t tofPatchMaxUS = TMath::Abs(tofPatchMax);
- if (tofPatchMaxUS < 25 ) fTriggerClusterBC = 0 ;
+ if (tofPatchMaxUS < 28 ) fTriggerClusterBC = 0 ;
else if(tofPatchMaxUS < 75 ) fTriggerClusterBC = 1 ;
else if(tofPatchMaxUS < 125) fTriggerClusterBC = 2 ;
else if(tofPatchMaxUS < 175) fTriggerClusterBC = 3 ;