- for (Int_t i=0; i<12; i++) {\r
- if( time[i] != 0 ){\r
- printf(" %i time %f \n", i, time[i] );\r
- fCFD[i]->Fill( time[i]);\r
- if( time[0] != 0 ) \r
- fTimeDiff[i]->Fill( time[i]-time[0]);\r
- }\r
- }\r
- for (Int_t i=12; i<24; i++) {\r
- if( time[i] != 0) {\r
- fCFD[i]->Fill( time[i]);\r
- if( time[12] != 0 ) \r
- fTimeDiff[i]->Fill( time[i]-time[12]);\r
+ const Double32_t* amp = fESD->GetT0amplitude();\r
+ \r
+ Double32_t diff;\r
+ for (Int_t i=0; i<24; i++) {\r
+ if( time[i] > 0 && amp[i]>0.1 ){\r
+ if (eq) {\r
+ fCFD[i]->Fill( time[i] );\r
+ if( time[fRefPMTC] > 0 && i<12) {\r
+ diff = time[i]-time[fRefPMTC];\r
+ fTimeDiff[i]->Fill( diff);\r
+ }\r
+ if( time[fRefPMTA] >0 && i>11) {\r
+ diff = time[i]-time[fRefPMTA] ;\r
+ fTimeDiff[i]->Fill( diff);\r
+ }\r
+ } //eq=1\r
+ else {\r
+ fCFD[i]->Fill( time[i] + fCDBdelays[i] );\r
+ if( time[fRefPMTC] > 0 && i<12) {\r
+ diff = time[i]-time[fRefPMTC] + fCDBdelays[i];\r
+ fTimeDiff[i]->Fill( diff);\r
+ } //C\r
+ if( time[fRefPMTA] >0 && i>11) {\r
+ diff = time[i]-time[fRefPMTA] + fCDBdelays[i];\r
+ fTimeDiff[i]->Fill( diff);\r
+ } //A\r
+ } //eq=0\r