[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.4.1  2002/06/10 17:51:14  hristov
Merged with v3-08-02

Revision 1.4  2002/04/24 22:08:12  nilsen
New ITS Digitizer/merger with two macros. One to make SDigits (old way) and
one to run the  merger (modified for Jiri).

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