+ aux=profnameResXvsPhi;
+ aux+=nhist;
+ fProfResXvsPhi[nhist]=new TProfile("histname","histname",40,-40.,40.0); // binning: range: -40°- 40°
+ fProfResXvsPhi[nhist]->SetName(aux.Data()); // bin width: 2°
+ fProfResXvsPhi[nhist]->SetTitle(aux.Data());
+
+ aux=profnameResZvsDip;
+ aux+=nhist;
+ fProfResZvsDip[nhist]=new TProfile("histname","histname",48,-72.,72.0); // binning: range: -70°-4°
+ fProfResZvsDip[nhist]->SetName(aux.Data()); // bin width: 3°
+ fProfResZvsDip[nhist]->SetTitle(aux.Data());
+
+ fProfResXvsPhiclu[nhist]=new TProfile*[kNclu];
+ fProfResZvsDipclu[nhist]=new TProfile*[kNclu];
+ for(Int_t clu=0; clu<kNclu; clu++) { // clu=0 --> cluster size 1
+ aux=profnameResXvsPhiclu;
+ aux+=nhist;
+ aux+="_clu_";
+ aux+=clu+1; // clu=0 --> cluster size 1
+ fProfResXvsPhiclu[nhist][clu]=new TProfile("histname","histname",40,-40.,40.0); // binning: range: -40°- 40
+ fProfResXvsPhiclu[nhist][clu]->SetName(aux.Data()); // bin width: 2°
+ fProfResXvsPhiclu[nhist][clu]->SetTitle(aux.Data());
+
+ aux=profnameResZvsDipclu;
+ aux+=nhist;
+ aux+="_clu_";
+ aux+=clu+1; // clu=0 --> cluster size 1
+ fProfResZvsDipclu[nhist][clu]= new TProfile("histname","histname",48,-72.,72.0); // binning: range: -70°-7°
+ fProfResZvsDipclu[nhist][clu]->SetName(aux.Data()); // bin width: 3°
+ fProfResZvsDipclu[nhist][clu]->SetTitle(aux.Data());
+ }
+
+ fHisTrackXFOtrue[nhist]=new TH1F*[kNClockPhase];
+ fHisTrackZFOtrue[nhist]=new TH1F*[kNClockPhase];
+ fHisTrackXFOfalse[nhist]=new TH1F*[kNClockPhase];
+ fHisTrackZFOfalse[nhist]=new TH1F*[kNClockPhase];
+ fHisTrackXZFOtrue[nhist]=new TH2F*[kNClockPhase];
+ fHisTrackXZFOfalse[nhist]=new TH2F*[kNClockPhase];
+ for(Int_t phas=0; phas<kNClockPhase;phas++){
+ aux=histnameTrackXFOtrue;
+ aux+=nhist;
+ aux+="_BCmod4_";
+ aux+=phas;
+ fHisTrackXFOtrue[nhist][phas]=new TH1F("histname","histname",128,-0.64,0.64); // +- 6.4 mm; 1 bin=0.1 mm
+ fHisTrackXFOtrue[nhist][phas]->SetName(aux.Data());
+ fHisTrackXFOtrue[nhist][phas]->SetTitle(aux.Data());
+
+ aux=histnameTrackZFOtrue;
+ aux+=nhist;
+ aux+="_BCmod4_";
+ aux+=phas;
+ fHisTrackZFOtrue[nhist][phas]=new TH1F("histname","histname",350,-3.5,3.5); // +- 35. mm; 1 bin=0.2 mm
+ fHisTrackZFOtrue[nhist][phas]->SetName(aux.Data());
+ fHisTrackZFOtrue[nhist][phas]->SetTitle(aux.Data());
+
+ aux=histnameTrackXFOfalse;
+ aux+=nhist;
+ aux+="_BCmod4_";
+ aux+=phas;
+ fHisTrackXFOfalse[nhist][phas]=new TH1F("histname","histname",128,-0.64,0.64); // +- 6.4 mm; 1 bin=0.1 mm
+ fHisTrackXFOfalse[nhist][phas]->SetName(aux.Data());
+ fHisTrackXFOfalse[nhist][phas]->SetTitle(aux.Data());
+
+ aux=histnameTrackZFOfalse;
+ aux+=nhist;
+ aux+="_BCmod4_";
+ aux+=phas;
+ fHisTrackZFOfalse[nhist][phas]=new TH1F("histname","histname",350,-3.5,3.5); // +- 35. mm; 1 bin=0.2 mm
+ fHisTrackZFOfalse[nhist][phas]->SetName(aux.Data());
+ fHisTrackZFOfalse[nhist][phas]->SetTitle(aux.Data());
+
+ aux=histnameTrackXZFOtrue;
+ aux+=nhist;
+ aux+="_BCmod4_";
+ aux+=phas;
+ fHisTrackXZFOtrue[nhist][phas]=new TH2F("histname","histname",22,-3.5,3.5,32,-0.64,0.64); // localZ +- 35. mm; 1 bin=3.2 mm
+ fHisTrackXZFOtrue[nhist][phas]->SetName(aux.Data()); // localX +- 6.4 mm; 1 bin=0.4 mm
+ fHisTrackXZFOtrue[nhist][phas]->SetTitle(aux.Data());
+
+ aux=histnameTrackXZFOfalse;
+ aux+=nhist;
+ aux+="_BCmod4_";
+ aux+=phas;
+ fHisTrackXZFOfalse[nhist][phas]=new TH2F("histname","histname",22,-3.5,3.5,32,-0.64,0.64); // localZ +- 35. mm; 1 bin=3.2 mm
+ fHisTrackXZFOfalse[nhist][phas]->SetName(aux.Data()); // localX +- 6.4 mm; 1 bin=0.4 mm
+ fHisTrackXZFOfalse[nhist][phas]->SetTitle(aux.Data());
+ }
+ } // end loop on module
+
+ TH1::AddDirectory(kTRUE);
+