[u/mrichter/AliRoot.git] / ITS / AliITSDigitizer.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
 
/*
$Log$
Revision 1.3  2002/03/25 10:48:55  nilsen
New ITS SDigit merging with region of interest cut. Update for changes in
AliDigitizer. Additional optimization should be done.

Revision 1.2  2002/03/15 17:26:40  nilsen
New SDigit version of ITS Digitizer.

Revision 1.1  2001/11/27 16:27:28  nilsen
Adding AliITSDigitizer class to do merging and digitization . Based on the
TTask method. AliITSDigitizer class added to the Makefile and ITSLinkDef.h
file. The following files required minor changes. AliITS, added functions
SetHitsAddressBranch, MakeBranchInTreeD and modified MakeBranchD.
AliITSsimulationSDD.cxx needed a Tree independent way of returning back to
the original Root Directory in function Compress1D. Now it uses gDirectory.

Revision 1.2  2002/03/01  E. Lopez
Digitization changed to start from SDigits instead of Hits.
The SDigits are reading as TClonesArray of AliITSpListItem
*/

#include <stdlib.h>
#include <iostream.h>
#include <TObjArray.h>
#include <TTree.h>
#include <TBranch.h>
#include <TFile.h>

#include <AliRun.h>
#include <AliRunDigitizer.h>

#include "AliITSDigitizer.h"
#include "AliITSpList.h"
#include "AliITSmodule.h"
#include "AliITSsimulation.h"
#include "AliITSDetType.h"
#include "AliITSgeom.h"

ClassImp(AliITSDigitizer)

//______________________________________________________________________
AliITSDigitizer::AliITSDigitizer() : AliDigitizer(){
    // Default constructor. Assign fITS since it is never written out from
    // here. 
    // Inputs:
    //      Option_t * opt   Not used
    // Outputs:
    //      none.
    // Return:
    //      A blank AliITSDigitizer class.

    fITS      = 0;
    fActive   = 0;
    fRoif     = -1;
    fRoiifile = 0;
    fInit     = kFALSE;
}
//______________________________________________________________________
AliITSDigitizer::AliITSDigitizer(AliRunDigitizer *mngr) : AliDigitizer(mngr){
    // Standard constructor. Assign fITS since it is never written out from
    // here. 
    // Inputs:
    //      Option_t * opt   Not used
    // Outputs:
    //      none.
    // Return:
    //      An AliItSDigitizer class.

    fITS      = 0;
    fActive   = 0;
    fRoif     = -1;
    fRoiifile = 0;
    fInit     = kFALSE;
}
//______________________________________________________________________
AliITSDigitizer::~AliITSDigitizer(){
    // Default destructor. 
    // Inputs:
    //      Option_t * opt   Not used
    // Outputs:
    //      none.
    // Return:
    //      none.

    fITS = 0; // don't delete fITS. Done else where.
    if(fActive) delete[] fActive;
}
//______________________________________________________________________
Bool_t AliITSDigitizer::Init(){
    // Initialization. Set up region of interest, if switched on, and
    // loads ITS and ITSgeom.
    // Inputs:
    //      none.
    // Outputs:
    //      none.
    // Return:
    //      none.

    fInit = kTRUE; // Assume for now init will work.
    if(!gAlice) {
	fITS      = 0;
	fActive   = 0;
	fRoif     = -1;
	fRoiifile = 0;
	fInit     = kFALSE;
	Warning("Init","gAlice not found");
	return fInit;
    } // end if
    fITS = (AliITS *)(gAlice->GetDetector("ITS"));
    if(!fITS){
	fActive   = 0;
	fRoif     = -1;
	fRoiifile = 0;
	fInit     = kFALSE;
	Warning("Init","ITS not found");
	return fInit;
    } else if(fITS->GetITSgeom()){
	//cout << "fRoif,fRoiifile="<<fRoif<<" "<<fRoiifile<<endl;
	fActive = new Bool_t[fITS->GetITSgeom()->GetIndexMax()];
    } else{
	fActive   = 0;
	fRoif     = -1;
	fRoiifile = 0;
	fInit     = kFALSE;
	Warning("Init","ITS geometry not found");
	return fInit;
    } // end if
    // fActive needs to be set to a default all kTRUE value
    for(Int_t i=0;i<fITS->GetITSgeom()->GetIndexMax();i++) fActive[i] = kTRUE;
/*  This will not work from Init. ts is aways returned as zero.
    TTree *ts;
    if(fRoif>=0 && fRoiifile>=0 && fRoiifile<GetManager()->GetNinputs()){
	ts = GetManager()->GetInputTreeS(fRoiifile);
	if(!ts){
	    if(!gAlice) ts = gAlice->TreeS();
	    if(!ts){
		cout <<"The TTree TreeS needed to set by region not found."
		    " No region of interest cut will be applied."<< endl;
		return fInit;
	    } // end if
	} // end if
	cout << "calling SetByReionOfInterest ts="<< ts <<endl;
	SetByRegionOfInterest(ts);
    } // end if
*/
    return fInit;
}
//______________________________________________________________________
void AliITSDigitizer::Exec(Option_t* opt){
    // Main digitization function. 
    // Inputs:
    //      Option_t * opt   list of sub detector to digitize. =0 all.
    // Outputs:
    //      none.
    // Return:
    //      none.

    char name[20] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
    char *all;
    const char *det[3] = {strstr(opt,"SPD"),strstr(opt,"SDD"),
                          strstr(opt,"SSD")};
    if( !det[0] && !det[1] && !det[2] ) all = "All";
    else all = 0;
    AliITSsimulation *sim      = 0;
    AliITSDetType    *iDetType = 0;
    static Bool_t    setDef    = kTRUE;

    if(!fInit){
	Error("Exec","Init not succesfull, aborting.");
	return;
    } // end if

    if( setDef ) fITS->SetDefaultSimulation();
    setDef = kFALSE;
    sprintf(name,"%s",fITS->GetName());

    Int_t nfiles = GetManager()->GetNinputs();
    Int_t event  = GetManager()->GetOutputEventNr();
    Int_t size   = fITS->GetITSgeom()->GetIndexMax();
    Int_t module,id,ifiles,mask;
    Bool_t lmod;
    Int_t *fl = new Int_t[nfiles];
    fl[0] = fRoiifile;
    mask = 1;
    for(id=0;id<nfiles;id++) if(id!=fRoiifile){
	// just in case fRoiifile!=0.
	fl[mask] = id;
	mask++;
    } // end for,if
    TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
    
    // Digitize
    fITS->MakeBranchInTreeD(GetManager()->GetTreeD());

    for(module=0; module<size; module++ ){
	if(fActive && fRoif!=0) if(!fActive[module]) continue;
        id = fITS->GetITSgeom()->GetModuleType(module);
        if(!all && !det[id]) continue;
        iDetType = fITS->DetType( id );
        sim      = (AliITSsimulation*)iDetType->GetSimulationModel();
        if(!sim) {
            Error( "Exec", "The simulation class was not instanciated!" );
            exit(1);
        } // end if !sim

        // Fill the module with the sum of SDigits
        sim->InitSimulationModule(module, event);
	//cout << "Module=" << module;
        for(ifiles=0; ifiles<nfiles; ifiles++ ){
	    if(fActive && fRoif!=0) if(!fActive[module]) continue;
	    //cout <<" fl[ifiles=" << ifiles << "]=" << fl[ifiles];
            TTree *treeS = GetManager()->GetInputTreeS(fl[ifiles]);
            if( !(treeS && fITS->GetSDigits()) ) continue;   
            TBranch *brchSDigits = treeS->GetBranch( name );
            if( brchSDigits ) {
                brchSDigits->SetAddress( &sdig ); 
            } else {
                Error( "Exec", "branch ITS not found in TreeS, input file %d ",
                       ifiles );
		return;
            } // end if brchSDigits
            sdig->Clear();
            mask = GetManager()->GetMask(ifiles);
            // add summable digits to module
            brchSDigits->GetEvent( module );
            lmod = sim->AddSDigitsToModule(sdig,mask);
	    if(ifiles==0){
		fActive[module] = lmod;
	    } // end if
	    //cout << " fActive["<<module<<"]=";
	    //if(fActive[module]) cout << "kTRUE";
	    //else cout << "kFALSE";
        } // end for ifiles
	//cout << " end ifiles loop" << endl;
        // Digitize current module sum(SDigits)->Digits
        sim->FinishSDigitiseModule();

        // fills all branches - wasted disk space
        GetManager()->GetTreeD()->Fill();
        fITS->ResetDigits();
    } // end for module
    //cout << "end modules loop"<<endl;

    GetManager()->GetTreeD()->AutoSave();

    delete[] fl;
    sdig->Clear();
    delete sdig;
    for(Int_t i=0;i<fITS->GetITSgeom()->GetIndexMax();i++) fActive[i] = kTRUE;
    return;
}
//______________________________________________________________________
void AliITSDigitizer::SetByRegionOfInterest(TTree *ts){
    // Scans through the ITS branch of the SDigits tree, ts, for modules
    // which have SDigits in them. For these modules, a flag is set to
    // digitize only these modules. The value of fRoif determines how many
    // neighboring modules will also be turned on. fRoif=0 will turn on only
    // those modules with SDigits in them. fRoif=1 will turn on, in addition,
    // those modules that are +-1 module from the one with the SDigits. And
    // So on. This last feature is not supported yet.
    // Inputs:
    //      TTree *ts  The tree in which the existing SDigits will define the
    //                 region of interest.
    // Outputs:
    //      none.
    // Return:
    //      none.
    Int_t m,nm,i;

    if(fRoif==0) return;
    if(ts==0) return;
    TBranch *brchSDigits = ts->GetBranch(fITS->GetName());
    TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
    //cout << "Region of Interest ts="<<ts<<" brchSDigits="<<brchSDigits<<" sdig="<<sdig<<endl;

    if( brchSDigits ) {
	brchSDigits->SetAddress( &sdig );
    } else {
	Error( "SetByRegionOfInterest","branch ITS not found in TreeS");
	return;
    } // end if brchSDigits

    nm = fITS->GetITSgeom()->GetIndexMax();
    for(m=0;m<nm;m++){
	//cout << " fActive["<<m<<"]=";
	fActive[m] = kFALSE; // Not active by default
	sdig->Clear();
	brchSDigits->GetEvent(m);
	if(sdig->GetLast()>=0) for(i=0;i<sdig->GetLast();i++){
	    // activate the necessary modules
	    if(((AliITSpList*)sdig->At(m))->GetpListItem(i)->GetSignal()>0.0){ // Must have non zero signal.
		fActive[m] = kTRUE;
		break;
	    } // end if
	} // end if. end for i.
	//cout << fActive[m];
	//cout << endl;
    } // end for m
    sdig->Clear();
    delete sdig;
    return;
}
Commit	Line	Data
9ad8b5dd	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/*
	17	$Log$
fd04285a	18	Revision 1.3 2002/03/25 10:48:55 nilsen
	19	New ITS SDigit merging with region of interest cut. Update for changes in
	20	AliDigitizer. Additional optimization should be done.
	21
fca85276	22	Revision 1.2 2002/03/15 17:26:40 nilsen
	23	New SDigit version of ITS Digitizer.
	24
f1888ca5	25	Revision 1.1 2001/11/27 16:27:28 nilsen
	26	Adding AliITSDigitizer class to do merging and digitization . Based on the
	27	TTask method. AliITSDigitizer class added to the Makefile and ITSLinkDef.h
	28	file. The following files required minor changes. AliITS, added functions
	29	SetHitsAddressBranch, MakeBranchInTreeD and modified MakeBranchD.
fca85276	30	AliITSsimulationSDD.cxx needed a Tree independent way of returning back to
f1888ca5	31	the original Root Directory in function Compress1D. Now it uses gDirectory.
	32
	33	Revision 1.2 2002/03/01 E. Lopez
fca85276	34	Digitization changed to start from SDigits instead of Hits.
f1888ca5	35	The SDigits are reading as TClonesArray of AliITSpListItem
9ad8b5dd	36	*/
	37
	38	#include <stdlib.h>
	39	#include <iostream.h>
	40	#include <TObjArray.h>
	41	#include <TTree.h>
	42	#include <TBranch.h>
fca85276	43	#include <TFile.h>
9ad8b5dd	44
	45	#include <AliRun.h>
	46	#include <AliRunDigitizer.h>
	47
	48	#include "AliITSDigitizer.h"
fd04285a	49	#include "AliITSpList.h"
9ad8b5dd	50	#include "AliITSmodule.h"
	51	#include "AliITSsimulation.h"
	52	#include "AliITSDetType.h"
	53	#include "AliITSgeom.h"
	54
	55	ClassImp(AliITSDigitizer)
	56
	57	//______________________________________________________________________
	58	AliITSDigitizer::AliITSDigitizer() : AliDigitizer(){
	59	// Default constructor. Assign fITS since it is never written out from
	60	// here.
	61	// Inputs:
	62	// Option_t * opt Not used
	63	// Outputs:
	64	// none.
	65	// Return:
	66	// A blank AliITSDigitizer class.
	67
fd04285a	68	fITS = 0;
	69	fActive = 0;
	70	fRoif = -1;
fca85276	71	fRoiifile = 0;
fd04285a	72	fInit = kFALSE;
9ad8b5dd	73	}
	74	//______________________________________________________________________
	75	AliITSDigitizer::AliITSDigitizer(AliRunDigitizer *mngr) : AliDigitizer(mngr){
	76	// Standard constructor. Assign fITS since it is never written out from
	77	// here.
	78	// Inputs:
	79	// Option_t * opt Not used
	80	// Outputs:
	81	// none.
	82	// Return:
	83	// An AliItSDigitizer class.
	84
fd04285a	85	fITS = 0;
	86	fActive = 0;
	87	fRoif = -1;
	88	fRoiifile = 0;
	89	fInit = kFALSE;
9ad8b5dd	90	}
	91	//______________________________________________________________________
	92	AliITSDigitizer::~AliITSDigitizer(){
	93	// Default destructor.
	94	// Inputs:
	95	// Option_t * opt Not used
	96	// Outputs:
	97	// none.
	98	// Return:
	99	// none.
	100
	101	fITS = 0; // don't delete fITS. Done else where.
f1888ca5	102	if(fActive) delete[] fActive;
9ad8b5dd	103	}
9ad8b5dd	104	//______________________________________________________________________
9ad8b5dd	105	Bool_t AliITSDigitizer::Init(){
fd04285a	106	// Initialization. Set up region of interest, if switched on, and
fd04285a	107	// loads ITS and ITSgeom.
9ad8b5dd	108	// Inputs:
	109	// none.
	110	// Outputs:
	111	// none.
	112	// Return:
	113	// none.
fca85276	114
fd04285a	115	fInit = kTRUE; // Assume for now init will work.
fca85276	116	if(!gAlice) {
	117	fITS = 0;
	118	fActive = 0;
	119	fRoif = -1;
	120	fRoiifile = 0;
fd04285a	121	fInit = kFALSE;
	122	Warning("Init","gAlice not found");
	123	return fInit;
fca85276	124	} // end if
	125	fITS = (AliITS *)(gAlice->GetDetector("ITS"));
	126	if(!fITS){
	127	fActive = 0;
	128	fRoif = -1;
	129	fRoiifile = 0;
fd04285a	130	fInit = kFALSE;
	131	Warning("Init","ITS not found");
	132	return fInit;
fca85276	133	} else if(fITS->GetITSgeom()){
fd04285a	134	//cout << "fRoif,fRoiifile="<<fRoif<<" "<<fRoiifile<<endl;
fca85276	135	fActive = new Bool_t[fITS->GetITSgeom()->GetIndexMax()];
fca85276	136	} else{
fd04285a	137	fActive = 0;
fd04285a	138	fRoif = -1;
fca85276	139	fRoiifile = 0;
fd04285a	140	fInit = kFALSE;
	141	Warning("Init","ITS geometry not found");
	142	return fInit;
fca85276	143	} // end if
	144	// fActive needs to be set to a default all kTRUE value
	145	for(Int_t i=0;i<fITS->GetITSgeom()->GetIndexMax();i++) fActive[i] = kTRUE;
fd04285a	146	/* This will not work from Init. ts is aways returned as zero.
fd04285a	147	TTree *ts;
fca85276	148	if(fRoif>=0 && fRoiifile>=0 && fRoiifile<GetManager()->GetNinputs()){
fca85276	149	ts = GetManager()->GetInputTreeS(fRoiifile);
fd04285a	150	if(!ts){
	151	if(!gAlice) ts = gAlice->TreeS();
	152	if(!ts){
	153	cout <<"The TTree TreeS needed to set by region not found."
	154	" No region of interest cut will be applied."<< endl;
	155	return fInit;
	156	} // end if
	157	} // end if
	158	cout << "calling SetByReionOfInterest ts="<< ts <<endl;
fca85276	159	SetByRegionOfInterest(ts);
	160	} // end if
	161	*/
fd04285a	162	return fInit;
9ad8b5dd	163	}
	164	//______________________________________________________________________
	165	void AliITSDigitizer::Exec(Option_t* opt){
fca85276	166	// Main digitization function.
9ad8b5dd	167	// Inputs:
fca85276	168	// Option_t * opt list of sub detector to digitize. =0 all.
9ad8b5dd	169	// Outputs:
	170	// none.
	171	// Return:
	172	// none.
f1888ca5	173
9ad8b5dd	174	char name[20] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
	175	char *all;
	176	const char *det[3] = {strstr(opt,"SPD"),strstr(opt,"SDD"),
	177	strstr(opt,"SSD")};
f1888ca5	178	if( !det[0] && !det[1] && !det[2] ) all = "All";
9ad8b5dd	179	else all = 0;
	180	AliITSsimulation *sim = 0;
	181	AliITSDetType *iDetType = 0;
f1888ca5	182	static Bool_t setDef = kTRUE;
9ad8b5dd	183
fd04285a	184	if(!fInit){
fd04285a	185	Error("Exec","Init not succesfull, aborting.");
f1888ca5	186	return;
9ad8b5dd	187	} // end if
9ad8b5dd	188
f1888ca5	189	if( setDef ) fITS->SetDefaultSimulation();
f1888ca5	190	setDef = kFALSE;
fd04285a	191	sprintf(name,"%s",fITS->GetName());
9ad8b5dd	192
f1888ca5	193	Int_t nfiles = GetManager()->GetNinputs();
	194	Int_t event = GetManager()->GetOutputEventNr();
	195	Int_t size = fITS->GetITSgeom()->GetIndexMax();
fd04285a	196	Int_t module,id,ifiles,mask;
	197	Bool_t lmod;
	198	Int_t *fl = new Int_t[nfiles];
	199	fl[0] = fRoiifile;
	200	mask = 1;
	201	for(id=0;id<nfiles;id++) if(id!=fRoiifile){
	202	// just in case fRoiifile!=0.
	203	fl[mask] = id;
	204	mask++;
	205	} // end for,if
f1888ca5	206	TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
f1888ca5	207
fd04285a	208	// Digitize
fd04285a	209	fITS->MakeBranchInTreeD(GetManager()->GetTreeD());
f1888ca5	210
fd04285a	211	for(module=0; module<size; module++ ){
	212	if(fActive && fRoif!=0) if(!fActive[module]) continue;
	213	id = fITS->GetITSgeom()->GetModuleType(module);
	214	if(!all && !det[id]) continue;
f1888ca5	215	iDetType = fITS->DetType( id );
9ad8b5dd	216	sim = (AliITSsimulation*)iDetType->GetSimulationModel();
fca85276	217	if(!sim) {
f1888ca5	218	Error( "Exec", "The simulation class was not instanciated!" );
9ad8b5dd	219	exit(1);
9ad8b5dd	220	} // end if !sim
fca85276	221
f1888ca5	222	// Fill the module with the sum of SDigits
fd04285a	223	sim->InitSimulationModule(module, event);
	224	//cout << "Module=" << module;
	225	for(ifiles=0; ifiles<nfiles; ifiles++ ){
	226	if(fActive && fRoif!=0) if(!fActive[module]) continue;
	227	//cout <<" fl[ifiles=" << ifiles << "]=" << fl[ifiles];
	228	TTree *treeS = GetManager()->GetInputTreeS(fl[ifiles]);
f1888ca5	229	if( !(treeS && fITS->GetSDigits()) ) continue;
	230	TBranch *brchSDigits = treeS->GetBranch( name );
	231	if( brchSDigits ) {
	232	brchSDigits->SetAddress( &sdig );
	233	} else {
	234	Error( "Exec", "branch ITS not found in TreeS, input file %d ",
	235	ifiles );
fd04285a	236	return;
f1888ca5	237	} // end if brchSDigits
f1888ca5	238	sdig->Clear();
fd04285a	239	mask = GetManager()->GetMask(ifiles);
f1888ca5	240	// add summable digits to module
f1888ca5	241	brchSDigits->GetEvent( module );
fd04285a	242	lmod = sim->AddSDigitsToModule(sdig,mask);
	243	if(ifiles==0){
	244	fActive[module] = lmod;
	245	} // end if
	246	//cout << " fActive["<<module<<"]=";
	247	//if(fActive[module]) cout << "kTRUE";
	248	//else cout << "kFALSE";
f1888ca5	249	} // end for ifiles
fd04285a	250	//cout << " end ifiles loop" << endl;
fca85276	251	// Digitize current module sum(SDigits)->Digits
f1888ca5	252	sim->FinishSDigitiseModule();
f1888ca5	253
9ad8b5dd	254	// fills all branches - wasted disk space
	255	GetManager()->GetTreeD()->Fill();
	256	fITS->ResetDigits();
	257	} // end for module
fd04285a	258	//cout << "end modules loop"<<endl;
	259
	260	GetManager()->GetTreeD()->AutoSave();
	261
	262	delete[] fl;
f1888ca5	263	sdig->Clear();
f1888ca5	264	delete sdig;
fd04285a	265	for(Int_t i=0;i<fITS->GetITSgeom()->GetIndexMax();i++) fActive[i] = kTRUE;
fd04285a	266	return;
9ad8b5dd	267	}
fca85276	268	//______________________________________________________________________
	269	void AliITSDigitizer::SetByRegionOfInterest(TTree *ts){
	270	// Scans through the ITS branch of the SDigits tree, ts, for modules
	271	// which have SDigits in them. For these modules, a flag is set to
	272	// digitize only these modules. The value of fRoif determines how many
	273	// neighboring modules will also be turned on. fRoif=0 will turn on only
	274	// those modules with SDigits in them. fRoif=1 will turn on, in addition,
	275	// those modules that are +-1 module from the one with the SDigits. And
	276	// So on. This last feature is not supported yet.
	277	// Inputs:
	278	// TTree *ts The tree in which the existing SDigits will define the
	279	// region of interest.
	280	// Outputs:
	281	// none.
	282	// Return:
	283	// none.
fd04285a	284	Int_t m,nm,i;
	285
	286	if(fRoif==0) return;
	287	if(ts==0) return;
fca85276	288	TBranch *brchSDigits = ts->GetBranch(fITS->GetName());
fca85276	289	TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
fd04285a	290	//cout << "Region of Interest ts="<<ts<<" brchSDigits="<<brchSDigits<<" sdig="<<sdig<<endl;
fca85276	291
fca85276	292	if( brchSDigits ) {
fd04285a	293	brchSDigits->SetAddress( &sdig );
fca85276	294	} else {
	295	Error( "SetByRegionOfInterest","branch ITS not found in TreeS");
	296	return;
	297	} // end if brchSDigits
	298
	299	nm = fITS->GetITSgeom()->GetIndexMax();
	300	for(m=0;m<nm;m++){
fd04285a	301	//cout << " fActive["<<m<<"]=";
fca85276	302	fActive[m] = kFALSE; // Not active by default
	303	sdig->Clear();
	304	brchSDigits->GetEvent(m);
fd04285a	305	if(sdig->GetLast()>=0) for(i=0;i<sdig->GetLast();i++){
	306	// activate the necessary modules
	307	if(((AliITSpList*)sdig->At(m))->GetpListItem(i)->GetSignal()>0.0){ // Must have non zero signal.
fca85276	308	fActive[m] = kTRUE;
fd04285a	309	break;
fca85276	310	} // end if
fd04285a	311	} // end if. end for i.
	312	//cout << fActive[m];
	313	//cout << endl;
fca85276	314	} // end for m
	315	sdig->Clear();
	316	delete sdig;
fd04285a	317	return;
fca85276	318	}