Re: Script failing at unexpected point

You are creating a 3-D histogram with 4096*4096*8192*4 = 550 GBytes !!! in TH3F *hTD1GE1PE1 = new
TH3F("hTD1_GE1_PE1","",4096,10.5,4090,4096,10.5,4090,8192,10.5,8191.5);

Rene Brun

Russell Leslie wrote:
>
> Greetings all,
>
> I am trying to do a basic sort of data that was collected from an
> experiment last week (I haven't had problems with related sorts in the
> last few months). The data collection system used four germanium
> detectors (GE 1 to 4), five silicon particle detectors (PE 1 to 5) and
> four time discriminator circuits (TD 1 to 4).
>
> My first task is to generate a series of cuts of coincidences between
> each GE, PE and TD (80 in all - 4x5x4). The first step is to just try
> the single case with GE1, PE1 and TD1. I can do a related manual sort
> using the treeviewer, but I need to be able to automate the process
> (for some 20 run files each which would need 80 cuts/sorts).
>
> The script example below gives the following error message ""Error:
> Can't call
> TH3F::Draw("X.TD1:X.GE1:X.PE1>>hTD1_GE1_PE1",gt1g1p1||pt1g1p1) in
> current scope sortRIV-T1G1P1.C:90:".
>
> I can't figure out why the script falls over at this point and unless
> I can get past this point I can't do the next important steps (doing
> corrections for random coincidences and compton background).
>
> Any guidance you can give on what I am doing wrong would be greatly
> appreciated.
>
> Regards
>
> Russell Leslie
> Nuclear Physics
> R.S.Phys.S.E.
> ANU
>
> {
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%% Original data file from the experiment in ROOT format
> %%%%%%%%%%%%%%%%%//
>
> TString fname="PD1060309_R05.root";
>
> TString tree_name="TREE";
>
> //
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%%%%GERMANIUM
> DETECTORS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> Double_t GE1_min = 750; Double_t GE1_max = 1250; // Sensible energy
> range for GE1
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%%%%PARTICLE
> DETECTORS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> Double_t PE1_min = 300; Double_t PE1_max = 3300; // Sensible energy
> range for PE1
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%% TDC limits for Random subtraction %%%%%%%%%%%%%%%%%%%//
>
> Double_t TD1_GE1_PE1_min = 3900; Double_t TD1_GE1_PE1_max = 4100; //
> Peak limits
>
> Double_t TD1_GE1_PE1_L1 = 3400; Double_t TD1_GE1_PE1_L2 = 3800; //
> Left background limits
>
> Double_t TD1_GE1_PE1_R1 = 4150; Double_t TD1_GE1_PE1_R2 = 4450; //
> Right background limits
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%% Total Gamma energy Cut %%%%%%%%%%%%%%%%%%%//
>
> TCut ra_GE1 = Form("X.GE1>%g && X.GE1<%g",GE1_min,GE1_max);
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%% Total Particle Cut %%%%%%%%%%%%%%%%%%%//
>
> TCut ra_PE1 = Form("X.PE1>%g && X.PE1<%g",PE1_min,PE1_max);
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%% Gamma/Particle TDC Cuts for Random Subtraction %%%%%%%//
>
> TCut t1g1p1 = Form("X.TD1_GE1_PE1>%g &&
> X.TD1_GE1_PE1<%g",TD1_GE1_PE1_min,TD1_GE1_PE1_max); //Block begin
>
> TCut t1g1p1L = Form("X.TD1_GE1_PE1>%g &&
> X.TD1_GE1_PE1<%g",TD1_GE1_PE1_L1,TD1_GE1_PE1_L2);
>
> TCut t1g1p1R = Form("X.TD1_GE1_PE1>%g &&
> X.TD1_GE1_PE1<%g",TD1_GE1_PE1_R1,TD1_GE1_PE1_R2);
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%% Individual cuts for Gamma Detectors
> %%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> TCut gt1g1p1 = ra_GE1 && t1g1p1;
>
> TCut gt1g1p1L = ra_GE1 && t1g1p1L;
>
> TCut gt1g1p1R = ra_GE1 && t1g1p1R;
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%% Individual cuts for Gamma Detectors
> %%%%%%%%%%%%%%%%%%%%%%//
>
> TCut pt1g1p1 = ra_PE1 && t1g1p1;
>
> TCut pt1g1p1L = ra_PE1 && t1g1p1L;
>
> TCut pt1g1p1R = ra_PE1 && t1g1p1R;
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%% {Normalization Factor = (peak width in
> TDC)/(background width in TDC)}%%%%%%%%%%%%%
>
> //%%%% Top & Bottom %%%%%%%//
>
> Double_t nrt1g1p1 = ((TD1_GE1_PE1_max - TD1_GE1_PE1_min +1) /
> ((TD1_GE1_PE1_L2 - TD1_GE1_PE1_L1 +1) + (TD1_GE1_PE1_R2 -
> TD1_GE1_PE1_R1 +1))); //Begin block
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%% Sorted Matrices for TD1 GE1 PD1 %%%%%%%%%%%%%%%%//
>
> TString fTD1_GE1_PE1 = "Mat_TD1_GE1_PE1_RIV_PD106_R05.root"; TString
> hTD1_GE1_PE1 = "hTD1_GE1_PE1"; //BEGIN BLOCK
>
> TString fTD1_GE1_PE1_TC = "Mat_TD1_GE1_PE1_RIV_PD106_TC_R05.root";
> TString hTD1_GE1_PE1_TC = "hTD1_GE1_PE1_TC";
>
> TString fTD1_GE1_PE1_C = "Mat_TD1_GE1_PE1_RIV_PD106_C_R05.root";
> TString hTD1_GE1_PE1_C = "hTD1_GE1_PE1_C";
>
>
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //&&&&&&&&&&&&&&&&&&&& FINAL DATA SORTING &&&&&&&&&&&&&&&&&&&&&&&&&&&&&//
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%%%%%%% TD1 GE1 PD1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //
>
> //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //%%%%%%%%%%%%%%%%%%%%% True + Chance %%%%%%%%%%%%%%%%%%%%%%%%%%%%//
>
> //
>
> TFile* f=TFile::Open( fname.Data(), "READ" );
>
> TTree *TREE = (TTree*)f.Get("tree_name");
>
> //
>
> TCanvas *C1 = new TCanvas("C1","Fitting canvas",10,10,400,200);
>
> C1->Divide(2,2);
>
> //
>
> C1->cd(1);
>
> TH3F *hTD1GE1PE1 = new
> TH3F("hTD1_GE1_PE1","",4096,10.5,4090,4096,10.5,4090,8192,10.5,8191.5);
>
> //
>
> hTD1GE1PE1->Draw("X.TD1:X.GE1:X.PE1>>hTD1_GE1_PE1",gt1g1p1||pt1g1p1);
> //Problem occurs here
>
> //
>
> hTD1GE1PE1->GetXaxis()->SetTitle("PE1");
>
> hTD1GE1PE1->GetYaxis()->SetTitle("GE1");
>
> hTD1GE1PE1->GetZaxis()->SetTitle("TD1");
>
> //
>
> TFile* f_out=TFile::Open( fhTD1_GE1_PE1_TC.Data(), "RECREATE" );
>
> hTD1GE1PE1->SetName(hTD1_GE1_PE1_TC.Data());
>
> f_out->cd();
>
> hTD1GE1PE1->Write();
>
> f_out->Close();
>
> //
>
> }
>
Received on Mon Mar 16 2009 - 08:15:14 CET