/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/root/HistoFileBuilder.C

Go to the documentation of this file.

//void HistoFileBuilder(TString path ; Tstring fileName )
{
gROOT->Reset();
gSystem->Load("./libMicro.so");

TString path = "./";
TString fileName= "analyse.root"; 

TFile *f = new TFile(path+fileName);
cout<<"reading File : "<<path<<fileName<<endl;
TTree *tree;
tree  = (TTree*) f.Get("t1");cout<<"tree built\n";
int nEntries = tree.GetEntries();
MTEvent *evt = new MTEvent();
TBranch *branch= tree->GetBranch("MTEvent");
branch->SetAddress(&evt);
//nEntries=500;
cout <<"nEntries="<<nEntries<<endl;


  //construire NC Histogrammes qui vont contenir le signal de chaque canal
  //THStack hs("HStack","Histo for each pad");
  THStack *hs = new THStack("HStack","Histo for each channel");
  TH1F *hh = new TH1F("hh","single channel histo",1024,0.,1023.);
  TString trackname;
  int NC = 672;
  for(int j=0;j<NC;j++)
    {
      trackname = "channel";
      trackname+=j;
      hs->Add((TH1F*)hh->Clone(trackname));
    }
  cout<<"histo stack ready\n";
  //boucle sur les evt pour remplir les histo
  //nEntries=20000;

  for( int i=0;i<nEntries;i++)
    {    if(i%1000==0){cout<<i<<endl;}
      tree.GetEntry(i);
        for(int j=0;j<evt->GetNchannel();j++)
          {
            MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
            TH1F *h=(TH1F*)hs.GetStack().UncheckedAt(channel->GetHardId());
            //if(channel.GetOrder()==0||channel.GetOrder()==96||channel.GetOrder()==192||channel.GetOrder()==288)
            {
              //cout<<i<<"\t"<<j<<"\t"<<channel.GetValue()<<endl;
              h.Fill(channel.GetValue());
            }
          }
    }
  cout<<"->histo filled\n";
//////////////////////////////////////////////////////////////////////////////////////////////////////
fileName="Histo_"+fileName;
cout<<"root file to store histo"<<path<<fileName<<endl;
TFile *fh = new TFile(path+fileName,"recreate");
cout<<"Writing in rootfile"<<endl;
for(int i=0;i<672;i++)
{
  trackname = "channel";
  trackname+=i;
  
  TH1F *h=(TH1F*)hs.GetStack().UncheckedAt(i);
  
  h->Write(trackname);
}
fh->Close();
cout<<"Histos written\n";
////////////////////////////////////////////////////////////////////////////////////////////////////


}

---