BDSIM: /home/cern/BDSIM_new/utils/display/powerELEMENT.C Source File

00001 powerELEMENT(char *gmadFile, char *rootFile)
00002 {
00003         gStyle->SetOptTitle(0);
00004         BDSCanvas *b = new BDSCanvas(gmadFile);
00005    TFile *f = TFile::Open(rootFile);
00006         TNtuple *ntuple = (TNtuple*)f->Get("ElossNtuple");
00007         Int_t i=1,n=1;
00008         Double_t coeff = 0.061;
00009         char temp[256];
00010         TH1F *myhisto= new TH1F("myhisto","Histo from the ntuple projection",
00011         b->GetNumberElement(),b->GetEndPositionElement());
00012         sprintf(temp,"E*%lf",coeff);
00013         ntuple->Project("myhisto","z",temp);
00014         for(Int_t i=1;i<=b->GetNumberElement();i++) 
00015 //      cout << "sample,range="<< b->GetNameElement()[i] << ";" << endl;
00016          cout << b->GetNameElement()[i] << " " << i << " " << b->GetEndPositionElement()[i] <<"  " << 
00017          myhisto->GetBinContent(i)<< endl;   
00018         
00019         for(Int_t i=1;i<b->GetNumberElement();i++)
00020         {
00021         myhisto->SetBinContent(i,myhisto->GetBinContent(i)/(b->GetEndPositionElement()[i]-b->GetEndPositionElement()[i-1]));
00022         cout << b->GetNameElement()[i]  << " " << b->GetTypeElement()[i] << endl;
00023         //++<< (b->GetEndPositionElement()[i]-b->GetEndPositionElement()[i-1]) << " " 
00024         //<< myhisto->GetBinContent(i)   << " "  
00025         //<<myhisto->GetBinContent(i)/(b->GetEndPositionElement()[i]-b->GetEndPositionElement()[i-1]) << endl;
00026         }
00027         
00028 //for(Int_t i=1;i<b->GetNumberElement();i++) 
00029 //      cout << b->GetNameElement()[i] << " " << i << " " << b->GetEndPositionElement()[i] <<"  " << 
00030 //      myhisto->GetBinContent(i)<< endl;       
00031         myhisto->GetXaxis()->SetLimits(BDSPad::xmin_axis,BDSPad::xmax_axis);
00032         myhisto->Draw();
00033         b->SetHisto(myhisto);
00034         b->cd();
00035 }
00036